Model Order Reduction for Quantum Molecular Dynamics
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
How is the orientation of molecular liquids ordered on cooling? What are the basic structures of molecular glasses, e.g., close to the crystalline structure or some special structures such as icosahedral cluster? These are long-standing questions in liquid and glass physics. Here, we have constructed a novel cryostat to prepare simple molecular glasses by vapor deposition and performed in situ synchrotron radiation x-ray diffraction experiments. The glassy state of a simple molecule CS 2 , which cannot be vitrified by normal liquid quenching, was successfully prepared with this instrument, and its diffraction data were collected in a wide Q-range of 0.16–25.7 Å -1 with a high-energy diffractometer at BL04B2, SPring-8. The diffraction data of liquid CS 2 were also recorded in a wide temperature range of 160–300 K. These diffraction data were analyzed with molecular dynamics simulations and reverse Monte Carlo modelings to investigate orientational correlation. From the obtained 3D structure models, the orientational correlation between neighboring CS 2 molecules was investigated quantitatively as a function of temperature. At room temperature, the parallel and T-shaped arrangements are preferred for the nearest neighbor correlation. On cooling, these arrangements are developed gradually, and its rate became prominent below the melting temperature (162 K). In the glassy state, the slipped-parallel arrangement is dominant as well as the T-shaped arrangement. Both arrangements appear in the CS 2 crystal, indicating that the structure of glassy CS 2 is close to that of crystalline CS 2 .
The structure–property relationships of small molecular acceptors (SMAs) are a key issue in the molecular design of new-generation acceptor materials for further improving the device efficiencies of polymer solar cells (PSCs). Herein, three couples of SMA isomers were synthesized, based on three central fused ring units and two 1,1-dicyanomethylene-3-indanone (IC) isomer electron-withdrawing terminal units with chlorine substitutions in different positions of its benzene ring: Cl-1 with chlorine on the same side as a C=O group of IC and Cl-2 with chlorine on the same side as the CN groups of IC. Through systematical investigation, we found that the chlorine substitution position of the terminal groups has a regular and significant influence on the molecular packing and photovoltaic performance of the SMAs. The molecular packing behavior of the SMAs is closely related to and determined by the configuration of their terminal groups, no matter which central fused ring of the SMAs is used. In particular, the Cl-1-based SMAs possess a stronger crystallinity with long range ordering packing in their molecular plane direction, while the more abundant and stereoscopic π–π intermolecular interaction in the Cl-2-based SMAs promotes the molecules to form three-dimensional charge transporting channels and leads to their red-shifted absorption and higher electron mobilities. Therefore, the Cl-2-based PSCs exhibit a higher power conversion efficiency (PCE) compared to that of the Cl-1-based devices, and the best PCE of a Cl-2 SMA-based PSC reached 16.42%. These results highlight the importance of the investigation of intermolecular interactions, packing and the arrangement of the SMAs in the solid-state, which may provide direct insights for exploring the relationship between the molecular structure and property of the photovoltaic materials. Moreover, we envision that if fragments such as end groups or side chains with more diverse molecular interactions are added into the design and the subsequent synthesis of the SMAs, this may be beneficial to promoting molecular π–π accumulation and further improving the molecular order, forming suitable molecular packing and morphology in the resulting blend films, and finally affecting the efficiency of the PSCs.
We present the framework and methodology for the new Los Alamos National Laboratory (LANL) G as chemical E quation O f S tate at H igher temperatures code (GEOSH) which aims to accurately model the behavior of chemically complex gaseous mixtures in equilibrium. Assuming the ideal gas approximation, GEOSH leverages the recursive nature of the Saha ionization and molecular equations in order to eliminate the molecular and ionic degrees of freedom, thereby reducing the problem size to the number of atomic species plus one for the free electrons if ions are included. This approach allows the chemical species, both molecular and ionic, of the mixture to be expressed in terms of the abundances of the elemental species. As a result, the GEOSH framework achieves a reduction in computational expense, increased processing speed, and the capability to efficiently model large-scale chemical networks. This report provides the necessary physical background and theoretical foundations for the GEOSH code, accompanied by benchmarking studies.
Elastic properties of amorphous solids remain a topic of intense interest due to their important roles in a wide range of applications. Prior works have focused on short-range, local structural features, such as number density and coordination number, to elucidate the underlying mechanism of elastic moduli in glasses. Here, we report strong correlations among the spatial extension of the interatomic potential, elastic moduli, and the coherence lengths of the medium-range structural order through molecular dynamics simulations for single element glasses. Our findings demonstrate an unconventional design principle to engineer elastic properties by extending the medium-range order and provide insights into the microscopic origin of elastic properties in amorphous solids.
Polyethylene is a promising low-cost alternative precursor material for carbon fiber production, but it has yet to show mechanical properties near or surpassing polyacrylonitrile-derived carbon fibers. The high molecular weight and order of ultra-high molecular weight polyethylene (UHMWPE) may offer a pathway to realizing this promise by enabling long-range graphitic structure formation and superior mechanical properties. The tension applied during the precursor stabilization process is crucial to maintaining the shape of the fibers during the conversion process, but no published study has yet probed the relationship between sulfonation tension and carbon fiber microstructure and mechanical properties. In this work, a logarithmic sweep of tensile stress was applied to UHMWPE fibers during the stabilization process followed by carbonization. Increasing tension significantly reduced fiber shrinkage, resulted in straighter fibers with less severe kink bands, and greatly improved the mechanical properties of the fibers. Raman spectroscopy and X-ray diffraction revealed that in all cases the carbon fibers were largely amorphous, but increasing tension resulted in increased size and alignment of the turbostratic crystallites with the fiber axis. Large voids were present in the sample fibers, so the Griffith-Irwin relation was employed to predict the potential ultimate tensile strength of the fibers with voids reduced to sizes comparable to commercially produced fibers. Furthermore, this work demonstrates the importance of tension applied during the stabilization of polyethylene fibers for carbon fiber production and establishes a framework for achieving high mechanical properties from these precursors.
ABSTRACT The performance of organic solar cells (OSCs) is governed by how molecular packing evolves into interconnected networks that facilitate exciton dissociation and charge transport. Using an all‐small‐molecule blend DR3TSBDT:Y6 as a model system, we study how local molecular stacking evolves into performance‐relevant morphology during solvent vapor annealing (SVA) and subsequent thermal annealing (TA). SVA promotes end‐to‐end stacking of amorphous acceptors to form interconnected fibrils, while TA compacts inter‐fibril spacing without disrupting favorable local order. Such molecular‐to‐morphological refinements broaden light absorption, enhance charge transport, and markedly improve device efficiency. Extending this approach to additional blend systems (D18:Y6, D18:L8‐BO, and DR3TSBDT:L8‐BO) yields similar structural evolution and performance gains, with the D18:L8‐BO system achieving up to 20.10% PCE. Our study establishes control over local stacking in amorphous acceptors into fibrillar networks as a general and effective route to realize high‐performance OSCs.
We review the role of nonperturbative phonon scattering in strongly anharmonic materials having ultralow lattice thermal conductivity with unusual temperature dependence. We take Tl 3 VSe 4 as an example and investigate its lattice dynamics using perturbation theory (PT) up to the fourth order and molecular dynamics (MD) with a machine-learning potential. We find distinct differences of phonon linewidth between PT and MD in the whole Brillouin zone. The comparison between the theoretical phonon linewidths and experiments suggests that PT severely underestimates the phonon scatterings, even when the fourth-order anharmonicity is included. Moreover, we extend our calculations to higher temperatures and evaluate the two-channel thermal conductivity based on the unified theory developed by Simoncelli et al. [Nat. Phys. 15, 809 (2019)]. We find a crucial coherence contribution to the total thermal conductivity at high temperatures. Our results pave the path for future studies of phonon properties and lattice thermal conductivities of strongly anharmonic crystals beyond the conventional PT realm.
We report in order to apply polymer semiconductors to stretchable electronics, they need to be easily deformed under strain without being damaged. A small number of conjugated polymers, typically with semicrystalline packing structures, have been reported to exhibit mechanical stretchability. Herein, a method is reported to modify polymer semiconductor packing-structure using a molecular additive, dioctyl phthalate (DOP), which is found to act as a molecular spacer, to be inserted between the amorphous chain networks and disrupt the crystalline packing. As a result, large-crystal growth is suppressed while short-range aggregations of conjugated polymers are promoted, which leads to an improved mechanical stretchability without affecting charge-carrier transport. Due to the reduced conjugated polymer intermolecular interactions, strain-induced chain alignment and crystallization are observed. By adding DOP to a well-known conjugated polymer, poly[2,5-bis(4-decyltetradecyl)pyrrolo[3,4-c]pyrrole-1,4-(2H,5H)-dione-(E)-1,2-di(2,2 '-bithiophen-5-yl)ethene] (DPPTVT), stretchable transistors are obtained with anisotropic charge-carrier mobilities under strain, and stable current output under strain up to 100%.
Covalent organic frameworks (COFs) are a new, growing class of crystalline polymers with exceptional characteristics such as low density, mechanical strength, and flexibility. Taking three-dimensional (3D) COF-300 derivatives as an example, we apply molecular dynamics simulations to show that ultrahigh thermal conductivity (>15 W/(mK)) is achieved at a pore size of d = 0.63 nm. This is the first observation of such an outstanding intrinsic thermal conductivity of a 3D polymer in all directions. Further temperature-dependent thermal conductivity calculations and X-ray diffraction (XRD) patterns highlight the importance of maintaining microscopic order for a large thermal conductivity. We also demonstrate that the tunability and anisotropy in COFs provide sufficient room to engineer and obtain desired thermal conductivity. These findings provide a fundamental understanding of structure-thermal conductivity relation in COFs and shine a light on high-thermal-conductivity COFs for thermal management applications.
Silicon (Si) atomic layer etching (ALE) by alternating exposure to chlorine gas (Cl 2 ) and argon ions (Ar + ) is studied by using molecular dynamics (MD) simulations and a reduced order model (ROM). Here, the purpose of this study is to elucidate the properties of the ALE window, a range of ion energies where the amount of Si etched over a series of cycles is nonzero and nearly independent of ion energy. Experimental studies of the Si–Cl 2 –Ar + ALE system report contradictory results related to the ALE window’s ion energy range. Both MD simulations and the ROM show that there is an ALE window present from approximately 15 to 20 eV for normal incidence argon ions. The Si–Cl 2 –Ar + system, therefore, exhibits a narrow ALE window. The amount of Si etched per cycle is less than one atomic layer because of the higher etch yield of Cl atoms relative to atomic Si and silicon chlorides. A modified version of the ROM with an artificially increased Si physical sputtering threshold energy expands the ALE window, illustrating the importance of the difference in chemical and physical sputtering threshold energies in the ALE window energy range. The ROM is also used to examine the dependence of the EPC on the Ar + ion fluence.
The competing and non-equilibrium phase transitions, involving dynamic tunability of cooperative electronic and magnetic states in strongly correlated materials, show great promise in quantum sensing and information technology. To date, the stabilization of transient states is still in the preliminary stage, particularly with respect to molecular electronic solids. Here, a dynamic and cooperative phase in potassium-7,7,8,8-tetracyanoquinodimethane (K-TCNQ) with the control of pulsed electromagnetic excitation is demonstrated. Simultaneous dynamic and coherent lattice perturbation with 8 ns pulsed laser (532 nm, 15 MW cm -2 , 10 Hz) in such a molecular electronic crystal initiates a stable long-lived (over 400 days) conducting paramagnetic state (≈42 Ωcm), showing the charge–spin bistability over a broad temperature range from 2 to 360 K. Comprehensive noise spectroscopy, in situ high-pressure measurements, electron spin resonance (ESR), theoretical model, and scanning tunneling microscopy/spectroscopy (STM/STS) studies provide further evidence that such a transition is cooperative, requiring a dedicated charge–spin–lattice decoupling to activate and subsequently stabilize nonequilibrium phase. Finally, the cooperativity triggered by ultrahigh-strain-rate (above 10 6 s -1 ) pulsed excitation offers a collective control toward the generation and stabilization of strongly correlated electronic and magnetic orders in molecular electronic solids and offers unique electro-magnetic phases with technological promises.
Complete tensors of the second- and third-order elastic constants of the organic molecular crystals acetaminophen, pentaerythritol tetranitrate (PETN), cyclotrimethylene trinitramine (RDX), cyclotetramethylene tetranitramine (HMX), 1,1-diamino-2,2-dinitroethylene (FOX-7), hexanitrohexaazaisowurtzitane (CL-20), and erythritol tetranitrate (ETN) have been calculated using dispersion-corrected density functional theory. The sets of second- and third-order elastic constants are expected to provide a more accurate and reliable description of the behavior of these materials under nonhydrostatic loads than pressure- and volume-dependent second-order elastic constants. The tensors of second-order constants have been compared with experimental data and/or other calculations when possible, and with the exception of results for CL-20 from Brillouin scattering experiments, we find good agreement. The calculated third-order elastic constants of PETN are in very good agreement with the subset of third-order constants derived from experimental wave speed measurements. The elastic anisotropies of the crystals have been estimated using the universal elastic anisotropy index, which shows that the crystals fall into three groups with low anisotropy (PETN, RDX, and CL-20), moderate anisotropy (acetaminophen, HMX, and ETN), and high elastic anisotropy (FOX-7).
Organophosphorus-based ligands represent a versatile set of solvent extraction reagents whose chemical makeup plays an important role in extraction mechanism. We hypothesize that the branching of the extractant hydrophobic tail and its oil-phase solvation affect the liquid/liquid interfacial structure. Understanding the structure mediated adsorption and interfacial ordering becomes key in designing ligands with enhanced selectivity and efficiency for targeted extractions. We employed vibrational sum frequency generation spectroscopy and interfacial tension measurements to extract thermodynamic adsorption energies, map interfacial ordering, and rationalize disparate behaviors of model di-(2-ethylhexyl) phosphoric acid and dioctyl phosphoric acid ligands at the hexadecane water interface. With increased surface loading, ligands with branched hydrophobic tails formed stable interfaces at much lower concentrations than those observed for ligands with linear alkyl tails. The lack of an oil phase and associated solvation results in markedly different interfacial properties, and thus measurements made at air/liquid surfaces cannot be assumed to correlate with the processes occurring at buried liquid/liquid interfaces. Here, we attribute these differences in the surface mediated self-assembly to key variations in hydrophobic interactions and tail solvation taking place in the oil phase demonstrating that interactions in both the polar and nonpolar phases are essential to understand self-assembly and function.
Abstract Here, a pair of A 1 –D–A 2 –D–A 1 unfused ring core‐based nonfullerene small molecule acceptors (NF‐SMAs), BO2FIDT‐4Cl and BT2FIDT‐4Cl is synthesized, which possess the same terminals (A 1 ) and indacenodithiophene unit (D), coupling with different fluorinated electron‐deficient central unit (difluorobenzoxadiazole or difluorobenzothiadiazole) (A 2 ). BT2FIDT‐4Cl exhibits a slightly smaller optical bandgap of 1.56 eV, upshifted highest occupied molecular orbital energy levels, much higher electron mobility, and slightly enhanced molecular packing order in neat thin films than that of BO2FIDT‐4Cl . The polymer solar cells (PSCs) based on BT2FIDT‐4Cl:PM7 yield the best power conversion efficiency (PCE) of 12.5% with a V oc of 0.97 V, which is higher than that of BO2FIDT‐4Cl ‐based devices (PCE of 10.4%). The results demonstrate that the subtle modification of A 2 unit would result in lower trap‐assisted recombination, more favorable morphology features, and more balanced electron and hole mobility in the PM7:BT2FIDT‐4Cl blend films. It is worth mentioning that the PCE of 12.5% is the highest value in nonfused ring NF‐SMA‐based binary PSCs with high V oc over 0.90 V. These results suggest that appropriate modulation of the quinoid electron‐deficient central unit is an effective approach to construct highly efficient unfused ring NF‐SMAs to boost PCE and V oc simultaneously.
The use of X-ray structures to determine and interpret the ferryl iron-oxygen bond order in molecular oxygen-activating heme enzymes has, in the past, been controversial. This has mainly stemmed from the susceptibility of ferryl species to X-ray-induced electronic state changes. In this work we establishe using time-resolved serial femtosecond X-ray crystallography (tr-SFX) on a dye-decolourising peroxidase that the ferryl intermediate species (Compounds I and II) captured following in situ mixing of microcrystals with H 2 O 2 have single, rather than the double bond character expected. X-ray emission validated tr-SFX data with quantum refinement, time-dependent-DFT calculations and QM/MM geometry optimizations together support the concept that the single iron-oxygen bond character is not an indication of ferryl reduction or a protonated form (Fe IV -OH) but is instead attributed to the existence of accessible excited states possessing ferric-oxyl (Fe III –O •– ) character. Such states offer insight into the nature of ferryl heme.
Employing dielectric spectroscopy, oscillatory shear rheology, and calorimetry, the present work explores the molecular dynamics of the widely used insecticide imidacloprid above and below its glass transition temperature. In its supercooled liquid regime, the applied techniques yield good agreement regarding the characteristic structural (alpha) relaxation times of this material. In addition, the generalized Gemant–DiMarzio–Bishop model provides a good conversion between the frequency-dependent dielectric and shear mechanical responses in its viscous state, allowing for an assessment of imidacloprid’s molecular hydrodynamic radius. In order to characterize the molecular dynamics in its glassy regime, we employ several approaches. These include the application of frequency–temperature superposition (FTS) to its isostructural dielectric and rheological responses as well as use of dielectric and calorimetric physical aging and the Adam–Gibbs–Vogel model. While the latter approach and dielectric FTS provide relaxation times that are close to each other, the other methods predict notably longer times that are closer to those reflecting a complete recovery of ergodicity. In conclusion, this seemingly conflicting dissimilarity demonstrates that the molecular dynamics of glassy imidacloprid strongly depends on its thermal history, with high relevance for the use of this insecticide as an active ingredient in technological applications.
Hund metals are multiorbital systems with 3d or 4d electrons exhibiting both an itinerant character and local moments, and they feature Kondo-like screenings of local orbital and spin moments, with suppressed coherence temperatures driven by Hund's coupling J H . They often exhibit magnetic order at low temperature, but how the interaction between the Kondo-like screening and long-range magnetic order is manifested in the quasiparticle spectrum remains an open question. Here, we present the spectroscopic signature of such an interaction in a Hund metal candidate MnSi exhibiting weak ferromagnetism. Our photoemission measurements reveal renormalized quasiparticle bands near the Fermi level with strong momentum dependence: The ferromagnetism manifests through possibly exchange-split bands (Q1) below T C , while the spin/orbital screenings lead to the gradual development of quasiparticles (Q2) upon cooling. Our results demonstrate how the characteristic spin/orbital coherence in a Hund metal could coexist and compete with the magnetic order to form a weak itinerant ferromagnet, via quasiparticle bands that are well separated in momentum space and exhibit distinct temperature dependence. Here, our results imply that the competition between the spin/orbital screening and the magnetic order in a Hund metal bears interesting similarities to the Kondo lattice systems