Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ab initio phasing”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Structure and Dynamics of CO 2 at the Air–Water Interface from Classical and Neural Network Potentials

The accurate description of the structure and dynamics of CO 2 at the instantaneous air–water interface, along with the effects of surface fluctuations on the CO 2 -transport processes, is essential for the development of negative emission technologies aimed at minimizing climate change. In this study, we performed molecular dynamics simulations of CO 2 at the air–water interface using neural network potentials (NNPs) trained on ab initio data generated through density-functional-theory-based molecular dynamics simulations. We compared these results with classical force fields to assess their performance in modeling interfacial CO 2 behavior. Our findings revealed that the asymmetric interactions, coupled with thermal surface fluctuations at the air–water interface, significantly influence CO 2 transport into the aqueous phase. The simulations demonstrate that classical force fields underestimate both the free energy of CO 2 transport and the strength of its interactions at the interface compared with the neural network potentials. In conclusion, the free energy and the interfacial dynamics of CO 2 are primarily influenced by the distribution of water within the instantaneous interfacial water layer, responsible for creating an asymmetric intermolecular interaction environment within the interfacial region.

Ab initio molecular dynamics↗
DOE OSTI2022

How arsenic makes amorphous GeSe a robust chalcogenide glass for advanced memory integration

The 3D integration technology in semiconductor fabrication requires a key component, the ovonic threshold switching (OTS) selector, to suppress the current leakage. The As doped amorphous (a-) GeSe glass is a commercialized OTS material in 3D phase-change memory, but the understanding of such a doping mechanism is still inadequate. Here we systematically explore the effect of As doping on the structural, bonding, and dynamics properties of a-GeAsSe using ab initio molecular dynamics simulations. The results reveal that As atoms form strong bonds with both Ge and Se atoms. The distorted octahedral structures and the 5-fold rings linked by atoms are increased. All of these structural features lead to a more disordered configuration. Moreover, as atoms have notably slowed down the atomic mobility, rendering a-GeAsSe a high stability. Overall, our studies offer insightful understanding of As-doping in OTS materials, paving the way for the design and application of advanced selector devices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Direct Extraction of Sodium Hydroxide by Calix[4]pyrrole-Based Ion-Pair Receptors

Here, described in this work, are calix[4]pyrrole-based ion-pair receptors, cis/trans-1 and cis/trans-2, designed for the extraction of sodium hydroxide. An X-ray diffraction analysis of a single crystal of the cis-1·NaOH isomer isolated from a mixture of cis/trans-1 revealed a unique dimeric supramolecular structure. An average dimer in toluene-d 8 solution was inferred on the basis of diffusion-ordered spectroscopy (DOSY). Support for the proposed stoichiometry came from density functional theory (DFT) calculations. The structural stability of the dimeric cis-1·NaOH complex in toluene solution was further confirmed by ab initio molecular dynamics (AIMD) simulation with explicit representation of solvent. Under conditions of liquid–liquid extraction (LLE), purified receptors cis- and trans-2 were both found to remove NaOH from a pH 11.01 aqueous source phase into toluene with extraction efficiencies (E%) of 50–60% when used equimolar to NaOH. However, in all cases, precipitation was observed. Complexities associated with precipitation could be avoided by immobilization of the receptors onto a chemically inert poly(styrene) resin by means of solvent impregnation. The use of solvent-impregnated resins (SIRs) eliminated precipitation in solution while retaining the extraction efficiency toward NaOH. This allowed both the pH and salinity of the alkaline source phase to be lowered.

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

Elucidating Norrish type I reactive pathways by ultrafast X-ray absorption spectroscopy

Norrish type I reactions selectively cleave carbon–carbon bonds directly adjacent to carbonyl groups. Despite their broad use in combination with aromatic carbonyls for additive manufacturing and dental UV curing applications, the nature of the photochemically active state and its population mechanism remain insufficiently understood. Detailed mechanistic insight requires mapping of the photoexcited population flow involving internal conversion and intersystem crossing. We present a time-domain study of gas phase acetophenone as a prototypical aromatic carbonyl combining soft X-ray time-resolved near-edge X-ray absorption fine structure (TR-NEXAFS) spectroscopy at the oxygen K-edge with ab initio multiple spawning (AIMS) simulations. Exploiting the specific sensitivity of TR-NEXAFS spectroscopy to states with nπ* character, we observe population transfer from the initially excited 1 ππ* state to the 1 nπ* state with a time constant of (0.13 ± 0.02) ps after an initial induction period of (0.12 ± 0.02) ps without population transfer, in quantitative agreement with the AIMS simulations. The population in the 1 nπ* state subsequently decays via intersystem crossing, likely mediated by a 3 ππ* state, within (3.17 ± 0.66) ps to a long-lived 3 nπ* state, which is presumed to be active towards Norrish type I chemistry.

Graßl, Martin [SLAC National Accelerator Laborator↗
DOE OSTI2024

Prediction challenge: First principles simulation of the ultrafast electron diffraction spectrum of cyclobutanone

Computer simulation has long been an essential partner of ultrafast experiments, allowing the assignment of microscopic mechanistic detail to low-dimensional spectroscopic data. However, the ability of theory to make a priori predictions of ultrafast experimental results is relatively untested. Herein, as a part of a community challenge, we attempt to predict the signal of an upcoming ultrafast photochemical experiment using state-of-the-art theory in the context of preexisting experimental data. Specifically, we employ ab initio Ehrenfest with collapse to a block mixed quantum–classical simulations to describe the real-time evolution of the electrons and nuclei of cyclobutanone following excitation to the 3s Rydberg state. The gas-phase ultrafast electron diffraction (GUED) signal is simulated for direct comparison to an upcoming experiment at the Stanford Linear Accelerator Laboratory. Following initial ring-opening, dissociation via two distinct channels is observed: the C3 dissociation channel, producing cyclopropane and CO, and the C2 channel, producing CH2CO and C2H4. Direct calculations of the GUED signal indicate how the ring-opened intermediate, the C2 products, and the C3 products can be discriminated in the GUED signal. We also report an a priori analysis of anticipated errors in our predictions: without knowledge of the experimental result, which features of the spectrum do we feel confident we have predicted correctly, and which might we have wrong?

Chemistry↗
DOE OSTI2020

Advances in ultrafast gas-phase x-ray scattering

Recent developments of x-ray free electron lasers and pulsed electron sources have enabled ultrafast scattering to become an increasingly powerful tool for exploring molecular dynamics. Here, this article describes our recent experimental and methodological advances in ultrafast gas-phase x-ray scattering experiments at the LCLS. A re-designed short-pathlength windowless diffractometer is coupled with careful optimization of sample density and independent normalization of x-ray intensity fluctuations to provide gas-phase scattering patterns with exceptionally high signal-to-noise ratios. These advances, coupled with careful geometry optimization and data treatment, provide both ground- and excited-state signals in excellent agreement with high level ab initio total scattering patterns.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Electron and Electron-Nuclear Dynamics of Open Quantum Subsystems

Leveraging an open-subsystem formulation of Density Functional Theory (DFT) we aim at describing periodic and molecular systems alike, including their electronic and nuclear dynamics. Subsystem DFT enables first principles simulations to approach realistic time- and length-scales, and most importantly sheds light on the dynamical behavior of complex systems. Taking subsystem DFT to the time domain allows us to inspect the electron dynamics of condensed-phase systems in real time. In liquids and interfaces, we observe all the relevant regimes proper of non-Markovian open quantum system dynamics, such as electronic energy transfer, and screening. In addition, the ab-initio modeling of system-bath interactions brought us to observe and justify the holographic time-dependent electron density theorem. Contrary to interactions between molecular (finite) systems, when molecules interact with metal or semiconductor surfaces the electron dynamics is strongly non-Markovian with dramatic repercussions to the molecule’s response to external perturbations. Metals and semiconductors typically have large polarizabilities, and even in a regime of low coupling their effect on impinging molecular species is significant–line broadening, peak shift, and intensity borrowing are observed, characterized, and explained in terms of inter-subsystem dynamical interactions and a many-body decomposition of the system’s density-density response function in a way that transcends the canons of Fermi Golden Rule.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Mobility assessment of the BCC and carbide phases in the C-Nb, C-U and Nb-U systems

Uranium carbides with refractory metal additions are considered for Gen IV nuclear reactors and nuclear thermal propulsion as fuels for their high-temperature and corrosion resistant properties. Understanding kinetic effects that dictate microstructural evolution during fabrication and operating conditions is essential to advance technological development of these fuels. This work presents the development of an atomic mobility database for C-Nb-U systems based off available experimental data supported with ab-initio methods. The mobility assessments and uncertainty quantification (using Markov chain Monte Carlo) were conducted in the Kawin software. Carbon diffusion is considered dominant, as metal diffusion is much slower, with niobium diffusion being even slower and rate limiting than uranium metal. We provide a comprehensive and self-consistent thermo-kinetic database that is validated by diffusion couple simulations through Kawin. In conclusion, this enables prediction of microstructural and phase evolution critical for the development and lifetime assessment of next generation nuclear fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Unveiling the structure and electronic characteristics of amorphous GeS for high performance threshold switching

Ovonic threshold switching selectors are widely studied owing to the essential application in high density phase-change memory. Amorphous GeS is proposed as a potential candidate for the excellent performance. However, the knowledge of amorphous GeS is still insufficient up to date. Here, we have studied the structure and electronic characteristics of GeS in the amorphization process, by using ab initio molecular dynamics simulations. The results indicate that the amorphous GeS is mainly made up of Ge–S bonds. The Ge- and S-centered clusters are dominantly in the form of octahedral structures in liquid GeS. Furthermore, during the amorphization process, most of Ge-centered clusters become highly coordinated octahedrons while a small number of Ge-centered clusters change to tetrahedrons, and the S-centered clusters deviate from the octahedral structure gradually. In addition, the large bandgap and the relatively small mid-gap states in amorphous GeS lead to a high switching voltage.

42 ENGINEERING↗
DOE OSTI2020

The Novel Hybrid Ab Initio Model of High-Performance Structural Alloys Design for Fossil Energy Power Plants

The current project developed a novel HT-CALPHAD/DFT approach, which can quickly design new high-performance structural alloys for the application of FE power plants. The PI will mainly take charge of high-throughput DFT simulations and computational thermodynamics of the selected multicomponent alloy systems for the FE power plant applications. At the end of the project, a novel hybrid model based on high throughput CALPHAD/DFT simulations and computational thermodynamics will be developed to provide guidance on how to identify multi-component new high-performance structural alloys with much less computational effort needed. It will address the extensive computation time needed for DFT on the new alloys design. In addition, it will also address the well-known headache of DFT, i.e. how to make the accurate prediction of the high-temperature equilibria. The novel hybrid model the PI proposed will not only be applied to the design of high-performance structural alloys in FE power plants but in many different applications, such as nuclear reactors. This hybrid modeling approach includes four sections: 1) HT-CALPHAD modeling of Al-Co-Cr-Ni-Fe system with FCC and BCC phase. In this section, 3561 non-equiatomic compositions were randomly generated in order to investigate the phase stability of single FCC and BCC phases. Meanwhile, we proposed a data screening procedure to screen out the good candidates within these compositions, considering the temperature range, average density, and melting temperature, etc. 2) Investigation of FCC-Cr lattice stability in Fe-Cr system. We systematically assessed the reliability of FCC-Cr lattice stability derived by DFT and CALPHAD approaches. Meanwhile, the Fe-Cr binary system was chosen as a case study to verify the Cr lattice stability obtained by both approaches. 3) High-throughput DFT modeling on elastic properties of Al-Co-Cr-Ni-Fe systems. We predicted and established the FCC quinary elastic constant database of the Al-Co-Cr-Fe-Ni systems at 0K by using special quasi-random structure (SQS) approach. The predictions will start with pure elements of Al-Co-Cr-Fe-Ni system and will be continued with binaries, ternaries, quaternaries, and finally the quinary compositions. 3) Modeling of temperature-dependent elastic properties in Al-Co-Cr-Ni-Fe systems. In this part, we predict the thermal expansion coefficient and elastic stiffness coefficient as a function of temperature by applying quasiharmonic approximation. With this approach, the elastic properties of HEAs at elevate temperature can be estimated.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

Prediction of chlorine and fluorine crystal structures at high pressure using symmetry driven structure search with geometric constraints

The high-pressure properties of fluorine and chlorine are not yet well understood because both are highly reactive and volatile elements, which have made conducting diamond anvil cell and x-ray diffraction experiments a challenge. Here, we use ab initio methods to search for stable crystal structures of both elements at megabar pressures. We demonstrate how symmetry and geometric constraints can be combined to efficiently generate crystal structures that are composed of diatomic molecules. Furthermore, our algorithm extends the symmetry driven structure search method by adding constraints for the bond length and the number of atoms in a molecule while still maintaining generality. As a method of validation, we have tested our approach for dense hydrogen and reproduced the known molecular structures of Cmca-12 and Cmca-4. We apply our algorithm to study chlorine and fluorine in the pressure range of 10 GPa–4000 GPa while considering crystal structures with up to 40 atoms per unit cell. We predict chlorine to follow the same series of phase transformations as elemental iodine from Cmca to Immm to Fm3¯m, but at substantially higher pressures. We predict fluorine to transition from a C2/c to Cmca structure at 70 GPa, to a novel orthorhombic and metallic structure with P42/mmc symmetry at 2500 GPa, and finally to its cubic analog form with Pm3¯n symmetry at 3000 GPa.

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

Multidimensional molecular high-harmonic spectroscopy: A road map for charge migration studies

This work discusses recent experimental and theoretical analyses of high-harmonic spectroscopy in small molecules, with the aim of characterizing charge migration. We discuss the formulation of molecular high-harmonic spectra followed by methods to extract molecular-target-specific information, both in experiments and ab initio simulations. We present measured and simulated high-harmonic spectra from carbon dioxide and carbonyl sulfide to illustrate the necessity for multidimensional analyses for high-harmonic spectroscopies that include both the spectral amplitude and phase. Leveraging these results, we examine how such multidimensional analyses pave the way for the study of charge migration with high-harmonic spectroscopy and illustrate the beneficial role a static molecular feature can play when probing dynamics. Finally, we briefly expand our scope with an outlook on the critical role of integrating theoretical and experimental approaches, beyond just high-harmonic spectroscopy, for the development of versatile harmonic spectroscopic probes of charge migration.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Insights into the Oxidation Mechanism of Vivianite to Metavivianite from First-Principles Calculations

Vivianite, a hydrous ferrous iron-phosphate mineral (Fe 3 (PO 4 ) 2 ·8H 2 O), readily oxidizes in contact with air yielding the less hydrous mixed-valent iron-phosphate mineral metavivianite. This topotactic transformation nominally occurs by oxidative dehydrogenation in which outgoing electrons from the iron sublattice are charge compensated by hydrolysis of structural water. However, the details of this internal charge balancing mechanism that allows the structure to remain electrostatically stable remain unknown. Here, in this study, we use density functional theory (DFT) calculations and ab initio thermodynamics (AIT) to evaluate the energetics of this process in terms of hydrogen release as a function of environmental variables such as partial pressure and temperature. The results show a thermodynamic driving force for vivianite phase transformation as oxidation progresses that is triggered by its rigid structure that has a limited accommodation for hydrogen vacancies. In contrast, metavivianite has a more flexible lattice and hydrogen bond network that stabilizes these hydrogen defects. Metavivianite is shown to be a stable intermediate for 66% residual Fe 2+ down to 33% Fe 2+ , below which other phases/structures such as santabarbarite should be more thermodynamically favorable. Our study provides a basis for experimental tests of our mechanistic findings and helps fill a basic knowledge gap about the solid-state process that defines how vivianite interacts with its surrounding environment.

Sassi, Michel [Pacific Northwest National Laborato↗
DOE OSTI2025

Designing a quantum-accurate machine-learning potential to enable large-scale simulations of deuterium under shock

Large-scale molecular dynamics of deuterium under shock can elucidate kinetic processes vital to the target design in inertial confinement fusion and high-energy-density experiments. However, modeling the complex evolution of this material from an insulating molecular state at ambient pressure to an ionized, atomic fluid under strong shock is beyond the capability of simple pair and even bond order potentials. We thus train a quantum-accurate and broadly transferable machine-learning interatomic potential for deuterium using the Chebyshev Interaction Model for Efficient Simulations framework. We show that due to an improved description of the molecular-to-atomic transition, our model is able to better reproduce the ab initio equation of state, radial distribution functions, and principal Hugoniot than bond order potentials. This represents an important step toward large-scale quantum-accurate and nonequilibrium simulations of complicated systems under dynamic changes including phase transitions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Imaginary phonon modes and phonon-mediated superconductivity in Y 2 C 3

For Y 2 C 3 with a superconducting critical temperature ($T_c$) ~18 K, zone-center imaginary optical phonon modes have been found for the high-symmetry I-43d structure due to C dimer wobbling motion and electronic instability from a flat band near Fermi energy. After lattice distortion to the more stable lowest-symmetry P1 structure, these stabilized low-energy phonon modes with a mixed C and Y character carry a strong electronphonon coupling to give rise to the observed sizable $T_c$. Our work shows that compounds with the calculated dynamical instability should not be simply excluded in the high-throughput search for new phonon-mediated superconductors. Moreover, we have studied the phase stability of the I-43d structure by calculating the enthalpy of different structural motifs of binary compounds containing group-IV elements at the 2:3 composition and also exploring the energy landscapes via ab initio molecular dynamics near and out of the I-43d structure. Here, our results show that the I-43d-type structures with C dimers are preferred in the low to medium pressure range. Because of the wobbling motion of the C dimers, there are many local energy minimums with degenerated energies. Thus, the ensemble average of many I-43d-distorted structures with C dimer wobbling motion at finite temperature still gives an overall I-43d structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Exploring non-adiabaticity to CO reduction reaction through ab initio molecular dynamics simulation

Non-adiabatic chemical reaction refers to the electronic excitation during reactions. This effect cannot be modeled by the ground-state Born-Oppenheimer molecular dynamics (BO-MD), where the electronic structure is at the ground state for every step of ions' movement. Although the non-adiabatic effect has been explored extensively in gas phase reactions, its role in electrochemical reactions, such as water splitting and CO 2 reduction, in electrolyte has been rarely explored. On the other hand, electrochemical reactions usually involve electron transport; thus, a non-adiabatic process can naturally play a significant role. In this work, using one-step CO 2 reduction as an example, we investigated the role of the non-adiabatic effect in the reaction. The reaction barriers were computed by adiabatic BO-MD and non-adiabatic real-time time dependent density functional theory (rt-TDDFT). We found that by including the non-adiabatic effect, rt-TDDFT could increase the reaction barrier up to 6% compared to the BO-MD calculated barrier when the solvent model is used to represent water. Simulations were carried out using explicit water molecules around the reaction site under different overpotentials, and similar non-adiabatic effects were found.

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

Unveiling the Finite Temperature Physics of Hydrogen Chains via Auxiliary Field Quantum Monte Carlo

The ability to accurately predict the finite temperature properties and phase diagrams of realistic quantum solids is central to uncovering new phases and engineering materials with novel properties ripe for device applications. Nonetheless, there remain comparatively few many-body techniques capable of elucidating the finite temperature physics of solids from first principles. Here, we take a significant step toward developing such a technique by generalizing our previous, exact fully ab initio finite temperature Auxiliary Field Quantum Monte Carlo (FT-AFQMC) method to model periodic solids and employing it to uncover the finite temperature physics of periodic hydrogen chains. Our chains’ unit cells consist of 10 hydrogen atoms modeled in a minimal basis, and we sample 5 k-points from the first Brillouin zone to arrive at a supercell consisting of 50 orbitals and 50 electrons. Based upon our calculations of these chains’ many-body energies, free energies, entropies, heat capacities, double and natural occupancies, and charge and spin correlation functions, we outline their metal–insulator and magnetic ordering as a function of both H–H bond distance and temperature. At low temperatures approaching the ground state, we observe both metal–insulator and ferromagnetic–antiferromagnetic crossovers at bond lengths between 0.5 and 0.75 Å. We then demonstrate how this low-temperature ordering evolves into a metallic phase with decreasing magnetic order at higher temperatures. In order to contextualize our results, we compare the features we observe to those previously seen in one-dimensional, half-filled Hubbard models at finite temperature and in ground state hydrogen chains. Interestingly, we identify signatures of the Pomeranchuk effect in hydrogen chains for the first time and show that spin and charge excitations that typically arise at distinct temperatures in the Hubbard model are indistinguishably coupled in these systems. Beyond qualitatively revealing the many-body phase behavior of hydrogen chains in a numerically exact manner without invoking the phaseless approximation, our efforts shed light on the further theoretical developments that will be required to construct the phase diagrams of the more complex transition metal, lanthanide, and actinide solids of longstanding interest to physicists.

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

Real-space visualization of short-range antiferromagnetic correlations in a magnetically enhanced thermoelectric

Short-range magnetic correlations can significantly increase the thermopower of magnetic semiconductors, representing a noteworthy development in the decades-long effort to develop high-performance thermoelectric materials. Here, we reveal the nature of the thermopower-enhancing magnetic correlations in the antiferromagnetic semiconductor MnTe. Using magnetic pair distribution function analysis of neutron scattering data, we obtain a detailed, real-space view of robust, nanometer-scale, antiferromagnetic correlations that persist into the paramagnetic phase above the Neel temperature $T_N$ = 307 K. In this work, the magnetic correlation length in the paramagnetic state is significantly longer along the crystallographic c axis than within the ab plane, pointing to anisotropic magnetic interactions. Ab initio calculations of the spin-spin correlations using density functional theory in the disordered local moment approach reproduce this result with quantitative accuracy. These findings constitute the first real-space picture of short-range spin correlations in a magnetically enhanced thermoelectric and inform future efforts to optimize thermoelectric performance by magnetic means.

36 MATERIALS SCIENCE↗