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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

High pressure chemical reactivity and structural study of the Na–P and Li–P systems

In this work, the Na–P and Li–P chemical systems were studied under pressure using synchrotron X-ray diffraction in a diamond anvil cell up to 20 GPa, combined with the AIRSS ab initio random structure searching technique. The results reveal an enhanced reactivity of both alkali metals with phosphorous at slightly elevated pressures. This enables the synthesis of Li 3 P and Na 3 P at room temperature (RT) starting from element precursors, bypassing the established chemical synthesis methods. Both compounds undergo a pressure-induced phase transition from the hexagonal Na 3 As-type structure (stable at ambient conditions) towards a Fm m (FCC) structure that remains stable up to 20 GPa. Attempts to synthesize compounds with higher alkali metal content (such as Li 5 P) using high-temperature and -pressure conditions (up to 2000+ K and 30 GPa), inspired by recent theoretical predictions, were not successful.

25 ENERGY STORAGE↗
DOE OSTI2023

Design of CO 2 Selective Type 3 Porous Liquids Through Porous Host Morphology

Direct air capture (DAC) of CO 2 is a negative emission technology under development to limit the impacts of climate change. The dilute concentration of CO 2 in the atmosphere (~400 ppm) requires new materials for carbon capture with increased CO 2 selectivity that is not met with current carbon capture materials. Porous liquids (PLs) are an emerging candidate for carbon capture and consists of a combination of solvents and porous hosts that creates a liquid with permanent porosity. The fundamental mechanisms of carbon capture in a PL are relatively unknown. To uncover these mechanisms, PLs were synthesized consisting of three different zeolitic-imidazolate framework (ZIF-8, ZIF-67, or ZIF-69) porous host in a water/glycol/2-methylimidazole solvent. The most stable composition was based on ZIF-8 and exhibited carbon capture following exposure to CO 2 . Density functional theory identified a three-step carbon capture mechanism based on (i) reaction of OH- with ethylene glycol in the solution followed by (ii) formation of 2-hydroxyethyl carbonate, which (iii) further react with OH- to form a carbonate species. This mechanism was validated with experimental nuclear magnetic resonance spectroscopy (NMR) to identify the dissolved carbonate phases and the decrease in the pH during CO 2 exposure. Deuterated samples of the ZIF-8 PLs were synthesized and analyzed via neutron diffraction at the Spallation Neutron Sources at Oak Ridge National Laboratory. Results identified differences in diffraction for PLs pre- and post-CO 2 exposure that will be combined with ab initio molecular dynamics data of the same PL composition to identify how the presence of a solvent-porous host interfaces results in carbon capture.

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

Quantifying uncertainties in α -nucleus reaction dynamics informed from first principles

The ab initio symmetry-adapted no-core shell model is a microscopic many-body method which naturally describes challenging collective and clustering features of atomic nuclei. Wave functions and observables computed with realistic nucleon-nucleon forces in this framework are tied to first principles, and are hence well-suited for rigorous uncertainty quantification. We discuss α-deuteron and α- 12 C cluster potentials informed by symmetry-adapted calculations, and propagate uncertainties in the effective binary cluster method as well as those in the nuclear interaction to reaction observables, namely scattering phase shifts, cross sections, partial widths and resonance energies. Here, we find that the overall uncertainties are dominated by those originating in the underlying nuclear force, speaking to the need for tighter constraints on realistic nucleon-nucleon interactions.

Ab initio↗
DOE OSTI2021

Nitromethane Decomposition via Automated Reaction Discovery and an Ab Initio Corrected Kinetic Model

In the explore the systematic construction of kinetic models from in silico reaction data for the decomposition of nitromethane. Our models are constructed in a computationally affordable manner by using reactions discovered through accelerated molecular dynamics simulations using the ReaxFF reactive force field. The reaction paths are then optimized to determine reaction rate parameters. We introduce a reaction barrier correction scheme that combines accurate thermochemical data from density functional theory with ReaxFF minimal energy paths. We validate our models across different thermodynamic regimes, showing predictions of gas phase CO and NO concentrations and high-pressure induction times that are similar to experimental data. The kinetic models are analyzed to find fundamental decomposition reactions in different thermodynamic regimes.

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

Photoinduced Electron–Nuclear Dynamics of Fullerene and Its Monolayer Networks in Solvated Environments

The recently synthesized monolayer fullerene network in a quasi-hexagonal phase (qHP-C 60 ) exhibits superior electron mobility and optoelectronic properties compared to molecular fullerene (C 60 ), making it highly promising for a variety of applications. However, the microscopic carrier dynamics of qHP-C 60 remain unclear, particularly in realistic environments, which are of significant importance for applications in optoelectronic devices. Unfortunately, traditional ab initio methods are prohibitive for capturing the real-time carrier dynamics of such large systems due to their high computational cost. In this work, we present the first real-time electron–nuclear dynamics study of qHP-C 60 using velocity-gauge density functional tight binding, which enables us to perform several picoseconds of excited-state electron–nuclear dynamics simulations for nanoscale systems with periodic boundary conditions. When applied to C 60 , qHP-C 60 , and their solvated counterparts, we demonstrate that water/moisture significantly increases the electron–hole recombination time in C 60 but has little impact on qHP-C 60 . Our excited-state electron–nuclear dynamics calculations show that qHP-C 60 is extremely unique and enable exploration of time-resolved dynamics for understanding excited-state processes of large systems in complex, solvated environments.

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

Probing the dynamics and bottleneck of the key atmospheric SO 2 oxidation reaction by the hydroxyl radical

SO 2 (Sulfur dioxide) is the major precursor to the production of sulfuric acid (H 2 SO 4 ), contributing to acid rain and atmospheric aerosols. Sulfuric acid formed from SO 2 generates light-reflecting sulfate aerosol particles in the atmosphere. This property has prompted recent geoengineering proposals to inject sulfuric acid or its precursors into the Earth’s atmosphere to increase the planetary albedo to counteract global warming. SO 2 oxidation in the atmosphere by the hydroxyl radical HO to form HOSO 2 is a key rate-limiting step in the mechanism for forming acid rain. However, the dynamics of the HO + SO 2 → HOSO 2 reaction and its slow rate in the atmosphere are poorly understood to date. Herein, we use photoelectron spectroscopy of cryogenically cooled HOSO 2 – anion to access the neutral HOSO 2 radical near the transition state of the HO + SO 2 reaction. Spectroscopic and dynamic calculations are conducted on the first ab initio-based full-dimensional potential energy surface to interpret the photoelectron spectra of HOSO 2 – and to probe the dynamics of the HO + SO 2 reaction. In addition to the finding of a unique pre-reaction complex (HO⋯SO 2 ) directly connected to the transition state, dynamic calculations reveal that the accessible phase space for the HO + SO 2 → HOSO 2 reaction is extremely narrow, forming a key reaction bottleneck and slowing the reaction rate in the atmosphere, despite the low reaction barrier. This study underlines the importance of understanding the full multidimensional potential energy surface to elucidate the dynamics of complex bimolecular reactions involving polyatomic reactants.

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

Photoelectron–residual-ion interaction in angle-resolved streaked shake-up ionization of helium

Streaked photoelectron emission spectra access the correlated dynamics of photoelectrons and residual target electrons with attosecond temporal resolution. Here, we calculated ab initio single-ionization spectra for photoemission from helium atoms by colinearly polarized ultrashort XUV and assisting few-femtosecond IR pulses. Distinguishing direct and shake-up ionization resulting in ground-state and excited (𝑛 = 2, 3) He + residual ions, respectively, we examined the effects of the correlated photoemission dynamics on the photoelectron phase accumulation as a function of the observable photoelectron detection direction and kinetic energy and XUV–IR pulse delay. We tracked the dynamical evolution of the residual ion in relative streaked photo- emission delays and found dominant contributions for shake-up emission from the residual-ion–photoelectron interaction. These are in very good and fair agreement, respectively, for 𝑛 = 2 and 𝑛 = 3 shake-up photoemission along the pulse-polarization directions, with previous experimental and theoretical investigations [Ossiander et al ., Nat. Phys. 13, 280 (2017)] and reveal a strong photoemission-direction dependence for shake-up ionization due to the coupling between the photoelectron and evolving residual-ion charge distribution in the IR-laser field.

Stark effect↗
DOE OSTI2020

TeraChem: A graphical processing unit-accelerated electronic structure package for large-scale ab initio molecular dynamics

TeraChem was born in 2008 with the goal of providing fast on-the-fly electronic structure calculations to facilitate ab initio molecular dynamics studies of large biochemical systems such as photoswitchable proteins and multichromophoric antenna complexes. Originally developed for videogaming applications, graphics processing units (GPUs) offered a low-cost parallel computer architecture that became more accessible for general-purpose GPU computing with the release of CUDA in 2007. The evaluation of the electron repulsion integrals (ERIs) is a major bottleneck in electronic structure codes and provides an attractive target for acceleration on GPUs. Thus, highly efficient routines for evaluation of and contractions between the ERIs and density matrices were implemented in TeraChem. Here, electronic structure methods were developed and implemented to leverage these integral contraction routines, resulting in the first quantum chemistry package designed from the ground up for GPUs. This GPU acceleration makes TeraChem capable of performing large-scale ground and excited state calculations in the gas and condensed phase. Today, TeraChem's speed forms the basis for a suite of quantum chemistry applications, including optimization and dynamics of proteins, automated and interactive chemical discovery tools, and large-scale nonadiabatic dynamics simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Modeling–Experiment–Theory Analysis of Reactions Initiated from Cl + Methyl Formate

Methyl formate (MF; CH 3 OCHO) is the smallest representative of esters, which are common components of biodiesel. The present study characterizes the thermal dissociation kinetics of the radicals formed by H atom abstraction from MF—CH 3 OCO and CH 2 OCHO—through a combination of modeling, experiment, and theory. For the experimental effort, excimer laser photolysis of Cl 2 was used as a source of Cl atoms to initiate reactions with MF in the gas phase. Time-resolved species profiles of MF, Cl 2 , HCl, CO 2 , CH 3 , CH 3 Cl, CH 2 O, and CH 2 ClOCHO were measured and quantified using photoionization mass spectrometry at temperatures of 400–750 K and 10 Torr. The experimental data were simulated using a kinetic model, which was informed by ab initio-based theoretical kinetics calculations and included chlorine chemistry and secondary reactions of radical decomposition products. Here, we calculated the rate coefficients for the H-abstraction reactions Cl + MF → HCl + CH 3 OCO (R1a) and Cl + MF → HCl + CH 2 OCHO (R1b): k 1a,theory = 6.71 × 10 –15 ·T 1.14 ·exp(—606/T) cm 3 /molecule·s; k 1b,theory = 4.67 × 10 –18 ·T 2.21 ·exp(—245/T) cm 3 /molecule·s over T = 200–2000 K. Electronic structure calculations indicate that the barriers to CH 3 OCO and CH 2 OCHO dissociation are 13.7 and 31.6 kcal/mol and lead to CH 3 + CO 2 (R3) and CH 2 O + HCO (R5), respectively. The master equation-based theoretical rate coefficients are k 3,theory (P = ∞) = 2.94 × 10 9 ·T 1.21 ·exp(—6209/T) s –1 and k 5,theory (P = ∞) = 8.45 × 10 8 ·T 1.39 ·exp(—15132/T) s –1 over T = 300–1500 K. The calculated branching fractions into R1a and R1b and the rate coefficient for R5 were validated by modeling of the experimental species time profiles and found to be in excellent agreement with theory. Additionally, we found that the bimolecular reactions CH 2 OCHO + Cl, CH 2 OCHO + Cl 2 , and CH 3 + Cl 2 were critical to accurately model the experimental data and constrain the kinetics of MF-radicals. Inclusion of the kinetic parameters determined in this study showed a significant impact on combustion simulations of larger methyl esters, which are considered as biodiesel surrogates.

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

Hole doping in a negative charge transfer insulator

RE NiO 3 is a negative charge transfer energy system and exhibits a temperature-driven metal-insulator transition (MIT), which is also accompanied by a bond disproportionation (BD) transition. In order to explore how hole doping affects the BD transition, we have investigated the electronic structure of single-crystalline thin films of Nd 1-x Ca x NiO 3 by synchrotron based experiments and ab-initio calculations. Here we show that for a small value of x, the doped holes are localized on one or more Ni sites around the dopant Ca 2+ ions, while the BD state for the rest of the lattice remains intact. The effective charge transfer energy (Δ) increases with Ca concentration and the formation of BD phase is not favored above a critical x, suppressing the insulating phase. Our present study firmly demonstrates that the appearance of BD mode is essential for the MIT of the RE NiO 3 series.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Second dome of superconductivity in YBa 2 Cu 3 O 7 at high pressure

Evidence is growing that a second dome of high-T c superconductivity can be accessed in the cuprates by increasing the doping beyond the first dome. Here, in this study, we use ab initio methods without invoking any free parameters, such as the Hubbard U, to reveal that pressure could turn YBa 2 ⁢Cu 3 ⁢O 7 into an ideal candidate for second dome superconductivity, displaying the predicted signature of strongly hybridized d x 2 -y 2 and d z 2 orbitals. Notably, pressure is found to induce a phase transition replacing the antiferromagnetic phases with an orbitally degenerate d-d phase. Our study suggests that the origin of the second dome is correlated with the oxygen-hole fraction in the CuO 2 planes and the collapse of the pseudogap phase.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Carrier Localization and Spontaneous Formation of Two-Dimensional Polarization Domain in Halide Perovskites

Halide perovskites are known for their rich phase diagram and superior performance in diverse optoelectronics applications. The latter property is often attributed to the long electron-hole recombination time, whose underlying physical mechanism has been a long-standing controversy. In this Letter, we investigate the transport and localization properties of electron and hole carriers in a prototypical halide perovskite (CsPbBr 3 ), through ab initio tight-binding nonadiabatic dynamics approach for large-scale (tens of nm size) supercell calculations, to simulate electron and ion dynamics on the same footing. Here, we found distinct structural, lattice polarization, and electron-phonon coupling properties at low (below 100 K) and high temperatures, consistent with experimental observations. In particular, at low temperature we find spontaneous formation of polar grain boundaries in the nonpolar bulk systems, which result in two-dimensional polarization patterns that serve to localize and separate electrons and holes. We reveal phonon-assisted variable-range hopping mostly responsible for low-temperature transport, and their characteristic frequency correlates with temperature-dependent phonon power spectrum and energy oscillation frequency in nonadiabatic dynamics. We answer the critical questions of long electron-hole recombination lifetime at low temperature and offer the correlation among polarization domains, electron-phonon couplings, and photocarrier dynamics.

Carrier generation & recombination↗
DOE OSTI2021

Orbital-resolved calculations of two-center interferences in linear triatomic molecules

In this report we perform ab initio calculations of high-harmonic spectroscopy (HHS) of two-center interference phenomena in oriented carbon-dichalcogen molecules, using time-dependent density functional theory (TDDFT). We show that by resolving the harmonic response into contributions from individual Kohn-Sham orbitals, we can extract target-specific characteristics for both the spectral amplitude and phase. We also discuss that this extraction is predicated on a careful analysis and normalization of the harmonic spectrum. Finally, we present field-free scattering calculations that mimic the recollision step in high-order-harmonic generation and allow the calculation of recombination dipole matrix elements without explicitly calculating the scattering states of a molecule. We show that the orbital-resolved TDDFT HHS results and results based on our field-free scattering calculations are in very good agreement with each other.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

The stability of FeH x and hydrogen transport at Earth’s core mantle boundary

Iron hydride in Earth’s interior can be formed by the reaction between hydrous minerals (water) and iron. Studying iron hydride improves our understanding of hydrogen transportation in Earth’s interior. Our high-pressure experiments found that face-centered cubic (fcc) FeH x (x≤1) is stable up to 165 GPa, and our ab initio molecular dynamics simulations predicted that fcc FeH x transforms to a superionic state under lower mantle conditions. In the superionic state, H-ions in fcc FeH become highly diffusive-like fluids with a high diffusion coefficient of ~3.7 × 10 -4 cm 2 s -1 , which is comparable to that in the liquid Fe-H phase. The densities and melting temperatures of fcc FeH x were systematically calculated. Similar to superionic ice, the extra entropy of diffusive H-ions increases the melting temperature of fcc FeH. Further, the wide stability field of fcc FeH enables hydrogen transport into the outer core to create a potential hydrogen reservoir in Earth’s interior, leaving oxygen-rich patches (ORP) above the core mantle boundary (CMB).

58 GEOSCIENCES↗
DOE OSTI2025

Predicting interface structure using the minima hopping method

Here, we adapt the minima hopping method (MHM) to the problem of interfacial structure prediction and apply it to study a canonical problem, the tilt grain boundaries in SrTiO 3 . Our method employs a hybrid approach by first exploring the potential energy surface (PES) of different grain boundary samplings with an empirical force field, among which the fifteen candidates with lower energies are then refined using ab initio density functional theory (DFT) calculations. During the exploratory stage, we bias the search using a local order parameter to primarily sample various reconstructions in the vicinity of the interface, while preserving the crystallinity of the bulk regions. We further enhance the search by incorporating initial structures with rigid body displacements to account for translational variations between bulk phases, enabling the MHM to effectively generate both stoichiometric and nonstoichiometric SrTiO 3 Σ⁢3(111)[110] and Σ⁢3(112)[110] grain boundaries. From an algorithmic standpoint, MHM outperforms earlier studies based on genetic algorithms (GA) by identifying more stable interfacial structures of several SrTiO 3 grain boundaries. The performance of the present implementation of the MHM approach is primarily limited by exploring an approximate description of the PES with a rather simple Buckingham potential. This limitation leads to variations in performance when compared to approaches utilizing more advanced surrogate PES models, such as direct DFT-PES sampling or GA with the embedded atom method (EAM). Despite the present limitations, the MHM approach is able to yield interfacial structures with comparable or lower interfacial energies in specific cases, such as Σ⁢3(111)[110] Γ=1, ±0.5 and Σ⁢3(112)[110] Γ= ±1, −2, underscoring the robustness of the MHM approach even with a simple approximation of the DFT PES. The MHM interfacial structure prediction method thus offers an efficient approach to understanding the grain boundaries and heterointerfaces at the atomic scale, providing an important prerequisite for effective materials design.

density functional theory↗
DOE OSTI2024

Breaking the Million-Electron and 1 EFLOP/s Barriers: Biomolecular-Scale Ab Initio Molecular Dynamics Using MP2 Potentials

The accurate simulation of complex biochemical phenomena has historically been hampered by the computational requirements of high-fidelity molecular-modeling techniques. Quantum mechanical methods, such as ab initio wave-function (WF) theory, deliver the desired accuracy, but have impractical scaling for modeling biosystems with thousands of atoms. Combining molecular fragmentation with MP2 perturbation theory, this study presents an innovative approach that enables biomolecular-scale ab initio molecular dynamics (AIMD) simulations at WF theory level. Leveraging the resolution-of-the-identity approximation for Hartree-Fock and MP2 gradients, our approach eliminates computationally intensive four-center integrals and their gradients, while achieving near-peak performance on modern GPU architectures. The introduction of asynchronous time steps minimizes time step latency, overlapping computational phases and effectively mitigating load imbalances. Utilizing up to 9,400 nodes of Frontier and achieving 59% (1006.7 PFLOP/s) of its double-precision floating-point peak, our method enables us to break the million-electron and 1EFLOP/s barriers for AIMD simulations with quantum accuracy.

Kurzak, Jakub↗
DOE OSTI2022

Mass transport in a highly immiscible alloy on extended shear deformation

Forced mixing to a single phase or supersaturated solid solution (SSS) and its prerequisite microstructure evolution in immiscible systems has been a focus of research for both fundamental science and a variety of applications. Controlling the formation of SSS by shear deformation assisted processing could enable a material design beyond conventional equilibrium microstructure in immiscible systems. Here, a highly immiscible (mixing enthalpy of ~ 20 kJ·mol -1 ) Cu-50 at. % Cr binary alloy was employed to investigate the microstructure evolution and localized tendencies of SSS during severe shear deformation. The present results demonstrate distinctive microstructure refinement process in each phase and how they lead to localized SSS as a function of shear strain. Preferential dynamic recrystallization occurs in the softer Cu phase due to strain localization, leading to substantial grain refinement. The refinement of Cr phase in the top-most layer, however, is enabled by the progressive evolution of grain lamination, splitting, and fragmentation as a function of shear strain. The eventual SSS is found to be strongly dependent on the local environments that affect the dislocation activity, including the level of microstructure refinement, the interfacial orientation relationship, hardness difference, and supposed stability of oxidation. Ab initio simulations comfirm that it is more favorable to oxidize Cr than Cu at incoherent Cu/Cr interfaces which then limits the mass transport on an incoherent boundary. Our results shed light on the underpinning mechanism for non-equilibrium mass transport in immiscible systems upon severe deformation that can be applicable to a variety of processing techniques aimed at producing immiscible alloys with superior mechanical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Origin of the magnetic and orbital ordering in α-Sr 2 CrO 4

Motivated by recent experimental progress in transition metal oxides with the K 2 NiF 4 structure, in this study we investigate the magnetic and orbital ordering in α-Sr 2 CrO 4 . Using first principles calculations, first we derive a three-orbital Hubbard model, which reproduces the ab initio band structure near the Fermi level. The unique reverse splitting of t 2g orbitals in α-Sr 2 CrO 4 , with the 3d 2 electronic configuration for the Cr 4+ oxidation state, opens up the possibility of orbital ordering in this material. Using real-space Hartree-Fock for multi-orbital systems, we constructed the ground state phase diagram for the two dimensional compound α-Sr 2 CrO 4 . We found stable ferromagnetic, antiferromagnetic, antiferro-orbital, and staggered orbital stripe ordering in robust regions of the phase diagram. Furthermore, using the density matrix renormalization group method for two-leg ladders with the realistic hopping parameters of α-Sr 2 CrO 4 , we explore magnetic and orbital ordering for experimentally relevant interaction parameters. Again, we find a clear signature of antiferromagnetic spin ordering along with antiferro-orbital ordering at moderate to large Hubbard interaction strength. We also explore the orbital-resolved density of states with Lanczos, predicting insulating behavior for the compound α-Sr 2 CrO 4 , in agreement with experiments. Finally, an intuitive understanding of the results is provided based on a hierarchy between orbitals, with d xy driving the spin order, while electronic repulsion and the effective one dimensionality of the movement within the d xz and d yz orbitals driving the orbital order.

36 MATERIALS SCIENCE↗