Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interface dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Exploring the Role of Inorganic and Organic Interfaces on CO 2 and CH 4 Partitioning: Case Study of Silica, Illite, Calcite, and Kerogen Nanopores on Gas Adsorption and Nanoscale Transport Behaviors

The adsorption, partitioning, and diffusion of CO 2 and CH 4 in organic–inorganic pores composed of kerogen with silica, illite, or calcite are studied using density functional theory and classical molecular dynamics simulations. The adsorption and partitioning behavior of CO 2 and CH 4 molecules is found to be a function of chemistry of the solid interface, pore size, and surface area. CO 2 molecules are preferentially adsorbed compared to CH 4 molecules on silica, illite, calcite, and kerogen surfaces due to significant contributions of electrostatic interactions with the atoms of the associated surfaces. CO 2 and CH 4 molecules display higher adsorption energy on calcite and silica compared to those on illite and kerogen because of enhanced interactions with the positively charged calcium ions on the calcite surface and the hydroxyl group (-OH) on the silica surface. Enhanced internal nanopores and the availability of adsorption active sites in illite and kerogen matrices aid in higher partitioning of the gas molecules into these pores. CO 2 has a significant influence on the clustering and swelling of kerogen fragments compared to CH 4 . Higher CO 2 adsorption on inorganic and kerogen surfaces results in lower overall self-diffusivities compared to CH 4 . This study illustrates surface and pore size controls on the adsorption, partitioning, and self-diffusivities of CO 2 and CH 4 at conditions relevant for subsurface energy applications. The chemistry of the surfaces along with the availability of internal pores and active adsorption sites influences the adsorption behavior of gases. Our results suggest that CO 2 stored in depleted gas reservoirs may preferentially inhabit nanoconfined pores.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Structural and dynamic mechanisms of GABA A receptor modulators with opposing activities

γ-Aminobutyric acid type A (GABA A ) receptors are pentameric ligand-gated ion channels abundant in the central nervous system and are prolific drug targets for treating anxiety, sleep disorders and epilepsy. Diverse small molecules exert a spectrum of effects on γ-aminobutyric acid type A (GABA A ) receptors by acting at the classical benzodiazepine site. They can potentiate the response to GABA, attenuate channel activity, or counteract modulation by other ligands. Structural mechanisms underlying the actions of these drugs are not fully understood. Here we present two high-resolution structures of GABA A receptors in complex with zolpidem, a positive allosteric modulator and heavily prescribed hypnotic, and DMCM, a negative allosteric modulator with convulsant and anxiogenic properties. These two drugs share the extracellular benzodiazepine site at the α/γ subunit interface and two transmembrane sites at β/α interfaces. Structural analyses reveal a basis for the subtype selectivity of zolpidem that underlies its clinical success. Molecular dynamics simulations provide insight into how DMCM switches from a negative to a positive modulator as a function of binding site occupancy. Together, these findings expand our understanding of how GABA A receptor allosteric modulators acting through a common site can have diverging activities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Spontaneous crack healing in calcite reveals the influence of dynamic strain evolution and surface chemistry

The mechanics of fracture healing in calcite remain poorly constrained yet are fundamental to managing fluid transport in geothermal reservoirs and hydrocarbon systems. Here, we apply microfocused synchrotron Laue X-ray diffraction and infrared spectroscopy to investigate subcritical crack healing in a 1 mm-thick calcite crystal subjected to controlled loading in a double-torsion device. Over a 44-hour period following load removal, we map the evolution of residual strain fields surrounding the crack tip and observe a progressive increase in compressive strain perpendicular to the crack plane accompanied by infrared spectroscopic signatures that reveal enhanced accumulation of water at the healed interface. The correlation between strain evolution and surface chemistry suggests that spontaneous crack healing in calcite is driven by dynamic anelastic relaxation coupled with irreversible fluid-mineral interactions. These findings offer insight into time-dependent crack closure processes in carbonates and highlight the role of chemically-mediated plasticity in subsurface fracture evolution.

materials science↗
DOE OSTI2021

Defect-Driven Oxide Transformations and the Electrochemical Interphase

The redox reaction pathway is crucial to the sustainable production of the fuels and chemicals required for a carbon-neutral society. Our society is becoming increasingly dependent on devices using batteries and electrolyzers, all of which rely on a series of redox reactions. The overall properties of oxide materials make them very well suited for such electrochemical and catalytic applications due to their associated cationic redox properties and the static site-adsorbate interactions. As these technologies have matured, it has become apparent that defect-driven redox reactions, defect-coupled diffusion, and structural transformations that are both time- and rate-dependent are also critical materials processes. This change in focus, considering not only redox properties but also more complex, dynamic behaviors, represents a new research frontier in the molecular sciences as they are strongly linked to device operation and degradation and lie at the heart of various phenomena that take place at electrochemical interfaces. Fundamental studies of the structural, electronic, and chemical transformation mechanisms are key to the advancement of materials and technological innovations that could be implemented in various electrochemical systems. In this Account, we focus on recent studies and advances in characterizing and understanding the dynamic redox evolution and structural transformations that take place in model perovskites and layered oxides under reactive conditions and correlate them with degradation mechanisms and operations in electrolyzers and batteries. We show that the dynamic evolution of oxygen vacancies and cationic migration in the surface or bulk occurs at the solid-liquid interface, using a combination of different synchrotron-based Xray spectroscopies and scattering probes. Detailed redox-structure-reactivity correlation studies show how defects and diffusion processes can be tailored to drive various physical and chemical transformations in electrolyzers and batteries. We also highlight a strong correlation between oxygen redox reactivity and structural reorganization in both model thin films and particles, helping to bridge the gap between fundamental studies of the reaction mechanism and device applications. Finally, on the basis of these findings, we discuss strategies to probe and tune the redox reactivity and structural stability of the redox-active oxide interphase toward devising efficient pathways for energy and chemical harvesting.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Studying Interfacial Dark Reactions of Glyoxal and Hydrogen Peroxide Using Vacuum Ultraviolet Single Photon Ionization Mass Spectrometry

Aqueous secondary organic aerosol (aqSOA) formation from volatile and semivolatile organic compounds at the air–liquid interface is considered as an important source of fine particles in the atmosphere. However, due to the lack of in situ detecting techniques, the detailed interfacial reaction mechanism and dynamics still remain uncertain. In this study, synchrotron-based vacuum ultraviolet single-photon ionization mass spectrometry (VUV SPI-MS) was coupled with the System for Analysis at the Liquid Vacuum Interface (SALVI) to investigate glyoxal dark oxidation products at the aqueous surface. Mass spectral analysis and determination of appearance energies (AEs) suggest that the main products of glyoxal dark interfacial aging are carboxylic acid related oligomers. Furthermore, the VUV SPI-MS results were compared and validated against those of in situ liquid time-of-flight secondary ion mass spectrometry (ToF-SIMS). The reaction mechanisms of the dark glyoxal interfacial oxidation, obtained using two different approaches, indicate that differences in ionization and instrument operation principles could contribute to their abilities to detect different oligomers. Therefore, the mechanistic differences revealed between the VUV SPI-MS and ToF-SIMS indicate that more in situ and real-time techniques are needed to investigate the contribution of the air–liquid interfacial reactions leading to aqSOA formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Quantifying motion blur by imaging shock front propagation with broadband and narrowband X-ray sources

Time-integrated radiography using MeV Bremsstrahlung X-ray sources is the norm for imaging during system-level testing of components and structures under dynamic condition. One source of error in the analysis of the time-integrated radiography data sets stems from motion blur which smears out sharp interfaces to a greater degree with longer exposure times, which become necessary to provide sufficient signal-to-noise with low X-ray penetration of objects of interest. To quantify motion blur, a 1D shock wave through PMMA was investigated experimentally at The Dynamic Compression Sector at The Advanced Photon Source (DCS@APS) with tapered broadband and 25.46 ± 1.06 keV narrowband X-rays. Four cameras with different exposure times were used for each experiment to compare the effect that exposure time has on motion blur. In addition, our methodology to accurately simulate motion blur in terms of transmission and shape is presented and compared to our experimental results and quantified. There is a high level of agreement between the experimental and simulation results across the range of data sets investigated in this study with a percent difference range of 0.29–1.31% for the four shots. The methodology of this work serves as a steppingstone towards a physically validated model that could be used in conjunction with experimental results to deconvolve physical parameters, densities, and interfaces of interest in a way that would not be possible with experimental results alone.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Experimental Setup and Learning-Based AI Model for Developing Accurate PV Inverter Models

The integration of power electronics-based interfaces presents challenges due to the absence of detailed models and the high computational complexity. Generic models used in system studies lack accuracy in capturing converter dynamics. This paper proposes a data-driven approach developed from experimental setup data. This approach enhances accuracy in photovoltaic inverter modeling. We used two types of PV inverters in the experiment. The recorded experimental data undergo processing through a machine learning model. Results from the model trained through machine learning is also presented.

artificial intelligence↗
DOE OSTI2020

Atomistic insights into the synergistic effects of tensile strain and hydroxyl group on increasing the thermal conductivity of monohydric alcohols as latent heat storage materials

The low thermal conductivity of linear monohydric alcohols significantly limits their performance when used as latent heat storage materials. It was found that imposing certain mechanical strains on nanomaterials with highly ordered structures can lead to appreciable thermal conductivity enhancement. To explore the potential of enhancing the intrinsic thermal conductivity of monohydric alcohols, a simplified ideal crystal model was employed for molecular arrangement in this work. The variation of atomistic heat transfer as a function of tensile strains was exploited by non-equilibrium molecular dynamics simulations. Furthermore, the influence of stretching on the vibration density of states and the corresponding molecular morphology were revealed. The results suggested that the thermal conductivity of monohydric alcohols with such simplified ideal crystalline structures can increase by about 170 % at the strain of 0.1. The distance between neighboring interfaces was found to determine the energy of inter-molecular interactions at the crystal grain boundaries. The hydrogen bonds energy has not reached the critical value when increasing tensile strains from 0 to 0.1. In contrast, the strongest van der Walls interactions for alkanes were obtained at the strain of 0.06 where the interfacial distances of alkanes have achieved their equilibrium distances. In addition, the interfacial heat transfer coefficient of ideal crystal monohydric alcohols at hydroxyl interfaces was found to be threefold than that between methyl interfaces, in consistent with the strength of inter-molecular interactions between functional groups at the interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction

This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressible solid structures for cases in which at least part of the solid boundary does not contact the incompressible fluid. Selected grid convergence studies demonstrate second-order convergence in volume conservation and in the pointwise discrepancies between corresponding positions of the two interface representations as well as between first and second-order convergence in the structural displacements. The time stepping scheme is also demonstrated to yield second-order convergence. To assess and validate the robustness and accuracy of the new algorithm, comparisons are made with computational and experimental FSI benchmarks. Test cases include both smooth and sharp geometries in various flow conditions. Furthermore, we also demonstrate the capabilities of this methodology by applying it to model the transport and capture of a geometrically realistic, deformable blood clot in an inferior vena cava filter.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Ultrafast Light-Induced Magnetoelectric Effect in van der Waals Magnetic Semiconductor Heterostructures

Atomic-scale heterostructures of van der Waals (vdW) magnets and semiconductors provide a unique environment for exploring magnetic dynamics. In contrast to typical photothermal excitation of precessional magnetization dynamics by a pump laser pulse, we find that ultrafast optical excitation of a WS2/CrGeTe3 (CGT) bilayer produces an opposite sign of magnetic torque compared to an isolated CGT film. Experimental observations by time-resolved magneto-optic Kerr effect (TR-MOKE) and theoretical analysis by density functional theory (DFT) and Landau-Lifshitz-Gilbert (LLG) simulations support a mechanism in which charge transfer of photoexcited carriers across the interface alters the perpendicular magnetic anisotropy, which in turn generates a torque on the magnetic layer to trigger precessional magnetization dynamics. These results provide new avenues for ultrafast manipulation of magnetization in vdW heterostructures with type-II band alignments. Lastly, we show that optically-generated spin currents from WS2 into CGT can also trigger precessional dynamics via angular momentum transfer.

Zhou, Wenyi [The Ohio State University]↗
DOE OSTI2022

An Engineered Multifunctional Composite for Passive Sensing, Power Harvesting, and In Situ Damage Identification with Enhanced Mechanical Performance

A multifunctional fiber-reinforced composite with passive self-sensing, energy-harvesting, and damage detection capabilities is presented. In this work, barium titanate piezoelectric microparticles are deposited on basalt fibers by a scalable, low-cost, environmentally friendly continuous feed-through process. The resulting composite derives a superior interlaminar shear strength from the microparticle-modified fiber–matrix interfaces. The composite also demonstrates passive self-sensing capabilities that produce electrical signals proportional to various dynamic loading events. Vibration and strain-controlled experiments are performed on composite beams to quantify the sensitivity and power output as a function of input acceleration and strain. Furthermore, these composite-generated electrical signals are used to identify in situ damage initiation for structural health monitoring to inform composite damage prior to structural failure. In brief, this truly multifunctional composite simultaneously displays a sensitivity of 0.5–2.6 mV g -1 at a resolution of 0.045–0.20g (g = gravitational acceleration), energy harvesting in the range of nW cc -1 , and prediction of early damage by exhibiting 0.017–1.17 mV peaks in voltage–time history profiles while assuring ≈20% improved interlaminar shear strength.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Temperature and density dependent pair potential for deuterium under shock

Large-scale classical molecular dynamics (CMD) simulations naturally include the microscopic physics necessary for atomistic modeling of shock release at the ablator-fuel interface in an inertial confinement fusion (ICF) capsule. Here, the multi-megabar shocks utilized in ICF experiments can drive the deuterium fuel from ambient to electron volt temperatures (T) and multi-fold compression. Modeling interatomic interactions over such an extreme range of conditions is challenging for empirical bond order potentials. We generate a pair potential for deuterium with explicit temperature and mass density dependence from ab initio density functional theory molecular dynamics using the iterative Boltzmann inversion method. This potential accurately reproduces the radial distribution functions and pressures from DFT in CMD equilibrium simulations across a wide range of thermodynamic conditions, yet fails to return the expected Hugoniot relations when used in direct CMD shock simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Probing Conformational Evolution and Associated Dynamics of Mg(N(SO 2 CF 3 ) 2 ) 2 ·Dimethoxyethane Adduct Using Solid-State 19 F and 1 H NMR

Bis(trifluoromethanesulfonimide) or TFSI is widely used as a counter anion in electrolyte design due to its structural flexibility and chemical stability. Here we studied the conformational variations and associated dynamics of TFSI in adduct of Mg(TFSI) 2 with dimethoxyethane (DME), a solvate crystalline material using solid-state 1 H and 19 F NMR. TFSI molecular motion in this solvate structure falls within the timescale of the 19 F NMR experiment, yielding spectroscopic signatures for unique TFSI conformers under the coordination environment of Mg 2+ cation. Within the temperature range of -5 to 82 °C, we observe nine distinct TFSI sites in both crystalline and disordered regions using 19 F NMR, reflecting complexity of structural and dynamics of TFS1 anions within solvate structure. The four distinguishable sites in the disordered region for the two CF 3 groups of the same TFSI molecule are identified using chemical shift analysis. The exchange rate constants from site to site are calculated through variable temperature 19 F NMR and two-dimensional (2D) exchange spectroscopy (EXSY) experiments, along with respective activation enthalpies using Eyring's formulation. The flip rate of CF 3 around the S-C bond is estimated as ~ 15 s -1 at 8 °C with ΔH ≠ ~ 22 kJ/mol, but the rotation of the entire TFSI is 4.8 s -1 at 8 °C with a significantly greater ΔH ≠ = 98 ± 10 kJ/mol. Furthermore, the slow conversion of trans to cis conformers at a lower temperature (T ≤ 1 °C) in the crystalline region is monitored, with a conversion rate of ~ 2 x 10 -5 s -1 at -5 °C. Density functional theory (DFT)-based calculations were performed to support further the assignment of experimental chemical shifts, and the activation energy E a = 21.1 kJ/mol obtained for the cis to trans transition is consistent with experimental values. The combined set of 19 F and 1 H under both one-dimensional (1D) and 2D NMR methods demonstrated here can be further used for examining electrode-electrolyte interfaces to probe the motions of various constituents that can enable detailed studies of interfacial processes and dynamics. Ultimately, such studies will aid in the design and discovery of interfacial constructs in which directed defect chemistry, chemical moiety distribution, and nanostructure are employed to drive efficient charge transport.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Quantifying transport and electrocatalytic reaction processes in a gastight rotating cylinder electrode reactor via integration of Computational Fluid Dynamics modeling and experiments

Understanding the complexity of the multiple processes of mass, momentum, charge, and heat transport, and how these affect reaction kinetics at the electrode/electrolyte interface is one of the major challenges in the field of energy and catalysis. The rapid and rational scale-up of electrocatalytic systems to industrial scales require a detailed understanding of nonlinear transport-reaction processes, accessible only through the building of multi-physics models that capture with high fidelity the complexity of real-world devices. The gastight rotating cylinder electrode (RCE) reactor is a promising lab-scale tool that can decouple transport from intrinsic kinetics to generate data for first-principle models useful in the design of industrial, electrochemical reactors. Computational Fluid Dynamics (CFD) studies have previously been used to investigate the bulk flow in RCE reactors for simple corrosion and electroplating processes. However, the quantification of changes in local concentration within the viscous layer where catalysis takes place requires capturing the correct flow conditions inside the hydrodynamic boundary layer near the surface of the electrode. Further, this requires simulations with spatial resolution in the nm and μm scale and temporal resolutions between ms and s scales that are similar to the timescales for reactions on the electrode surface. In this study, experimental electrocatalysis is combined with CFD modeling to elucidate and parameterize the hydrodynamics in a gastight RCE reactor. CFD simulations of the electrochemical ferricyanide reduction reaction under mass transport limited conditions are used to evaluate the validity of the CFD model parameters by comparing calculated dimensionless mass transport descriptors to dimensionless correlations obtained experimentally. Justifications for assumptions and details of the simulation methods used in this study are presented to provide a detailed understanding of the effect that each model parameter has on the ability to accurately simulate electrocatalysis in RCE systems. The simulation methodology reported here is a first step towards the development of multi-scale models for the study of transport dependent electrocatalytic processes, such as the electrochemical transformation of CO 2 to fuels and chemicals.

42 ENGINEERING↗
DOE OSTI2026

pH-Driven Restructuring of Hydration Layers and Cation Adsorption at the Alumina–Water Interface

Oxide−water interfaces underpin ion separations, catalysis, and electrochemical energy technologies, where the electrical double layer (EDL) controls adsorption, transport, and reactivity. However, the molecular-scale links between pHdependent surface protonation, hydration-layer structure, and counterion adsorption remain poorly defined. Here, we combine in situ crystal truncation rod (CTR) and resonant anomalous X-ray reflectivity (RAXR) with streaming potential measurements and ab initio molecular dynamics (AIMD) simulations to resolve the chemical and structural evolution of the EDL at the single-crystal alumina (012)−water interface in 10 mM Rb + over pH 3−12. CTR measurements reveal two distinct adsorbed water layers at ∼2.2 and ∼3.5 Å above the surface. Each water layer shifts toward the substrate at transition pHs near 6.5 and 10.6, respectively, reflecting changes in primary hydration layer structure in response to the deprotonation of bridging and terminal aluminol groups. RAXR shows a 10-fold increase in Rb + coverage and a decrease in mean adsorption height from ∼3.5 to ∼2.7 Å with increasing pH, indicating enhanced counterion binding accompanied by Stern layer contraction. Streaming potential measurements demonstrate that the zeta potential, i.e., the potential at the hydrodynamic shear plane, is positive at pH 3 and becomes negative at pH ≥ 3.5. The negative charge magnitude increases with increasing pH, consistent with progressive surface deprotonation at higher pH. AIMD identifies inner- and outer-sphere Rb+ complexes whose adsorption heights and coordination geometries depend sensitively on the protonation state of surface oxygens, providing atomistic support for the experimentally inferred trends. These measurements establish two discrete, site-specific pH transitions in hydration-layer structure that track aluminol (de)protonation and quantitatively link them to a pH-driven contraction of the Stern layer (increasing Rb + coverage and decreasing adsorption height). This provides a direct structural basis for connecting surface acid−base chemistry to ion binding distances at an oxide−water interface.

Adsorption↗
DOE OSTI2021

Imaging the spin chirality of ferrimagnetic Néel skyrmions stabilized on topological antiferromagnetic Mn 3 Sn

Néel skyrmions are generally realized in asymmetric multilayers made of heavy metals (HMs) and ultrathin ferromagnets possessing strong interfacial Dzyaloshinskii-Moriya interactions (iDMIs). Depending on the relative strengths of iDMIs at the interfaces, various types of Néel skyrmions have been suggested, which are typified with characteristically different topological properties and current-driven dynamics. This suggests the importance of a precise quantification of their spin chiralities. In this paper, we explore the possibility of realizing Néel skyrmions in magnetic multilayers without the direct usage of standard HMs. Specifically, through depositing a thin layer of ferrimagnetic (FIM) CoTb layer on top of an antiferromagnetic (AFM) quantum material of composition Mn 3 Sn, the AFM exchange interaction at the asymmetric interface provides an equivalent iDMI for stabilizing FIM Néel skyrmions. Secondly, through using advanced four-dimensional Lorentz scanning transmission electron microscopy (4D LSTEM), in combination with x-ray magnetic circular dichroism photoemission electron microscopy (XMCD-PEEM), we can directly determine the spin chirality of FIM Néel skyrmions. Finally, the present findings not only broaden the phase space for chiral interfacial magnetism but also provide a possibility for future applications of heavy-metal-free skyrmionic devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Diabatization with Electrostatic Embedding for Studying Photophysics in Organic Molecular Crystals

Highly emissive organic molecular crystals find applications in several areas, such as organic electronics, solar cells, and sensors. Understanding the excited-state mechanisms underlying these applications is essential for optimizing and controlling them effectively. Exciton models coupled with nonadiabatic dynamics, particularly quantum dynamics, provide crucial insights into photochemical and photophysical processes in molecular crystals. Nevertheless, there remains a lack of general tools and automated workflows to facilitate such simulations. In this paper, we present a computational strategy to investigate the photoactivated dynamics of organic molecular crystals, bridging methodologies traditionally used for molecular systems and materials science, with a particular focus on the interplay between local excitations and charge transfer (CT) processes. We have implemented an interface between the fromage and Overdia programs, enabling the construction of vibronic Hamiltonians for molecular crystals within an excited-state ONIOM(QM:QM′) framework, incorporating long-range electrostatics through a RESP-based Ewald summation. Fragment-based diabatization provides a route to quantum dynamics simulations in weak-to-intermediate coupling regimes. The method was applied to the photophysics of dibenzo[g,p]chrysene (DBC) crystals using time-dependent DFT. The fromage/ Overdia interface was employed to compute the couplings of local excitations and CT states for 18 unique DBC dimers in the crystal and to quantify the influence of electrostatic embedding, which was found to be modest (10−20%). Simulations on π-stacked dimers reproduced the small red shift observed experimentally from solution to crystal, attributed to electronic interactions among fixed monomers rather than crystal electrostatics. Quantum dynamics simulations revealed ultrafast population transfer from bright local excitations to CT states. This approach establishes a robust framework linking molecular and solid-state excited-state dynamics, with potential applications for studying excitations, defects, and impurities in molecular crystals.

Crystals↗
DOE OSTI2022

Easy-plane spin Hall nano-oscillators as spiking neurons for neuromorphic computing

Here we show analytically using a macrospin approximation that easy-plane spin Hall nano-oscillators excited by a spin current polarized perpendicularly to the easy plane have phase dynamics analogous to that of Josephson junctions. Similarly to Josephson junctions, they can reproduce the spiking behavior of biological neurons that is appropriate for neuromorphic computing. To take advantage of typical spin-orbit torques, we use a nanoconstriction geometry, in which the magnetostatic interaction and magnetocrystalline anisotropy are tuned to create an easy plane that includes the interface normal direction. We perform micromagnetic simulations of such oscillators realized in this geometry and show that the easy-plane spiking dynamics is preserved in this experimentally feasible architecture. Finally we simulate two elementary neural network blocks that implement operations essential for neuromorphic computing. First, we show that output spikes energies from two neurons can be summed and injected into a following layer neuron and second, we demonstrate that outputs can be multiplied by synaptic weights implemented by locally modifying the anisotropy.

36 MATERIALS SCIENCE↗