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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Theoretical analysis of the role of complex transition dipole phase in XUV transient-absorption probing of charge migration

We theoretically investigate the role of complex dipole phase in the attosecond probing of charge migration. The iodobromoacetylene ion (ICCBr + ) is considered as an example, in which one can probe charge migration by accessing both the iodine and bromine ends of the molecule with different spectral windows of an extreme-ultraviolet (XUV) pulse. The analytical expression for transient absorption shows that the site-specific information of charge migration is encoded in the complex phase of cross dipole products for XUV transitions between the I-4 d and Br-3 d spectral windows. Ab-initio quantum chemistry calculations on ICCBr + reveal that there is a constant π phase difference between the I-4 d and Br-3 d transient-absorption spectral windows, irrespective of the fine-structure energy splittings. Transient absorption spectra are simulated with a multistate model including the complex dipole phase, and the results correctly reconstruct the charge-migration dynamics via the quantum beats in the two element spectral windows, exhibiting out-of-phase oscillations.

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

Ordering in liquid and its heredity impact on phase transformation of Mg-Al-Ca alloys

It is a long-sought goal to achieve desired mechanical properties through tailoring phase formation in alloys, especially for complicated multi-phase alloys. In fact, unveiling nucleation of competitive crystalline phases during solidification hinges on the nature of liquid. Here we employ ab initio molecular dynamics simulations (AIMD) to reveal liquid configuration of the Mg-Al-Ca alloys and explore its effect on the transformation of Ca-containing Laves phase from Al 2 Ca to Mg 2 Ca with increasing Ca/Al ratio (r Ca/Al ). There is structural similarity between liquid and crystalline phase in terms of the local arrangement environment, and the connection schemes of polyhedras. The forming signature of Mg 2 Ca, as hinted by the topological and chemical short-range order originating from liquid, ascends monotonically with increasing r Ca/Al . However, Al 2 Ca crystal-like order increase at first and then decrease at the crossover of r Ca/Al = 0.74, corresponding to experimental composition of phase transition from Al 2 Ca to Mg 2 Ca. The origin of phase transformation across different compositions lies in the dense packing of atomic configurations and preferential bonding of chemical species in both liquid and solid. The present finding provides a feasible scenario for manipulating phase formation to achieve high performance alloys by tailoring the crystal-like order in liquid.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The Post‐Perovskite Transition in Fe‐ and Al‐Bearing Bridgmanite: Effects on Seismic Observables

Abstract The primary phase of the Earth’s lower mantle, (Al, Fe)‐bearing bridgmanite, transitions to the post‐perovskite (PPv) phase at Earth’s deep mantle conditions. Despite extensive experimental and ab initio investigations, there are still important aspects of this transformation that need clarification. Here, we address this transition in (Al 3+ , Fe 3+ )‐, (Al 3+ )‐, (Fe 2+ )‐, and (Fe 3+ )‐bearing bridgmanite using ab initio calculations and validate our results against experiments on similar compositions. Consistent with experiments, our results show that the onset transition pressure and the width of the two‐phase region depend distinctly on the chemical composition: (a) Fe 3+ ‐, Al 3+ ‐, or (Al 3+ , Fe 3+ )‐alloying increases the transition pressure, while Fe 2+ ‐alloying has the opposite effect; (b) in the absence of coexisting phases, the pressure‐depth range of the Pv‐PPv transition is likely too broad to cause a sharp D” discontinuity (<30 km); (c) the average Clapeyron slope of the two‐phase regions are consistent with previous measurements, calculations in MgSiO 3 , and inferences from seismic data. In addition, (d) we observe a softening of the bulk modulus in the two‐phase region. The consistency between our results and experiments gives us the confidence to proceed and examine this transition in aggregates with different compositions computationally, which will be fundamental for resolving the most likely chemical composition of the D" region by analyses of tomographic images.

Geochemistry & Geophysics↗
DOE OSTI2022

Revisit the VEC criterion in high entropy alloys (HEAs) with high-throughput ab initio calculations: A case study with Al-Co-Cr-Fe-Ni system

Valence electron concentration (VEC) was treated as a useful parameter to predict the stability of solid solution phases. However, the available experimental data to support this criterion is far from enough. In the current study, the high-throughput ab initio modeling is applied to investigate the relative stability of FCC and BCC single crystals of the Al-Co-Cr-Fe-Ni high entropy alloys (HEAs) by using the special quasi-random structure (SQS) approach. Furthermore, the predictions start with pure elements of the Al-Co-Cr-Fe-Ni system and are continued with binaries, ternaries, and quaternary compositions, which come up with 180 compositions (360 structures). After that, the reliability of the VEC criterion is testified. The results show that the VEC criterion not only works for the stable structure but also works effectively for metastable structure when both FCC and BCC are not thermodynamic stable. However, it is found that the old VEC criterion proposed by Guo et al. fails to work effectively for compositions containing high concentrations of light-weight metals such as Al at VEC<5. To solve this problem, the present work proposes a new VEC rule to define the stability of FCC and BCC structures at the ground state. With the implementation of the new VEC rule, the effectiveness of the VEC rule (EVEC) of both FCC and BCC structures is enhanced, especially for pure elements and binary compositions, indicating that this rule does not only work effectively for multicomponent systems but also works for low-order systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ab Initio Modeling of fcc Fe-Co-Cr-Ni High Entropy Alloys with Full Composition Range

The Fe-Co-Cr-Ni quaternary system has been studied extensively in the past decade, not only because of the superior properties achieved like high tensile ductility and fracture toughness, but also it is the foundation for the development of FeCrCoNi-based HEAs. However, most of the investigations are mainly focusing on the equiatomic and semi-equiatomic compositions. The physical properties of non-equiatomic compositions are barely explored. In the current work, the ab initio approach is adopted to predict the lattice parameter, structure stabilities, elastic properties, and enthalpy of formations of fcc Fe-Cr-Co-Ni HEA single crystals with the special quasirandom structure (SQS) method. To expand the design of Fe-Co-Cr-Ni HEAs into a full composition range, the current simulations cover all the binaries, ternaries, and the whole quaternary system. The 2-D composition-property contour diagrams are created to show the physical properties in the ternary and quaternary systems. Furthermore, our investigation shows that the fcc alloys' lattice parameter increases with the addition of Fe and Cr and decreases with Ni and Co. In addition to the enthalpy of formation at o K, the shape distortion rate is proposed to be another essential parameter to evaluate crystal stability. The results show the addition of Cr would destabilize the fcc lattice and cause symmetry breaking. Finally, the result of Pugh’s ratio shows that most of the fcc alloys show ductile behavior (P > 1.75), especially for the alloys around Fe 3 Ni, FeCo, and Fe 2 NiCo (P > 2.5). The brittle behavior (P < 1.75) is located around Fe 3 Cr and Fe 2 CoCr. Meanwhile, Cr 2 NiCo and Cr 2 FeNiCo are considered promising compositions to be tried and verified experimentally because of the higher bulk modulus and moderate shear modulus.

36 MATERIALS SCIENCE↗
DOE OSTI2024

HTESP (High-throughput electronic structure package): A package for high-throughput ab initio calculations

High-throughput ab initio calculations are the indispensable parts of data-driven discovery of new materials with desirable properties, as reflected in the establishment of several online material databases. The accumulation of extensive theoretical data through computations enables data-driven discovery by constructing machine learning and artificial intelligence models to predict novel compounds and forecast their properties. Efficient usage and extraction of data from these existing online material databases can accelerate the next stage materials discovery that targets different and more advanced properties, such as electron–phonon coupling for phonon-mediated superconductivity. However, extracting data from these databases, generating tailored input files for different ab initio calculations, performing such calculations, and analyzing new results can be demanding tasks. Here, in this work, we introduce a software package named “HTESP” (High-Throughput Electronic Structure Package) written in Python and Bash languages, which automates the entire workflow including data extraction, input file generation, calculation submission, result collection and plotting. Our HTESP will help speed up future computational materials discovery processes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Zinc-based cyclens containing pyridine and cross-bridges: X-ray and DFT structures, Lewis acidity, gas-phase acidity, and p K a values

Herein it is well known that catalytic centers containing the zinc(II) ion can act as both Lewis and Bronsted-Lowry acids. In addition to coordination number, the Zn coordination geometry can also strongly impact the acidity of the active site, no matter what measure of acidity is considered. Herein, we report the first pentacoordinate, zinc-ammonia complex containing a pyridine based tetraazamacrocycle, [(pyclen)Zn(NH 3 )](PF 6 ) 2 , that has applications in Lewis acid catalysis. From this structure, we obtain binding energies and acidities for a series of related pyclen and cross-bridged cyclen type tetraazamacrocycles comprising the pentacoordinate N 4 Zn(II)–OH 2 entity, collectively referred to as [(R-cyclen)Zn–OH 2 ] 2+ . Results from gas- and aqueous-phase density functional theory (M05-2X) and ab initio (MP2) calculations reported herein demonstrate that molecular geometry has a substantial influence on both the Zn–OH x binding strengths and the deprotonation energy of coordinated H 2 O, but not necessarily in a predicable way. While generally more constrained N–Zn–N coordination leads to greater Zn–OH 2 binding energies (Lewis acidities), the corresponding Lewis acidities of the complexes don’t always correlate with the Zn–OH 2 bond lengths nor the (Bronsted) acidity of the coordinated H 2 O. Additionally, the order of the Lewis acid strength of the [(R-cyclen)Zn(II)] complexes changes as the basic –OH 2 ligand is replaced with its –OH counterpart.

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

Experimental and Theoretical Study of Oxolan-3-one Thermal Decomposition

The thermal decomposition of oxolan-3-one, a common component of the bio-oil formed during biomass pyrolysis, has been studied using ab initio calculations and experiments employing pulsed gas-phase pyrolysis with matrix isolation FTIR product detection. Four pathways for unimolecular decomposition were predicted using computational methods. The dominant reaction channel led to carbon monoxide, formaldehyde, and ethylene, all of which were observed experimentally. Here, the other channels led to an assortment of products including ketene, water, propyne, and acetylene, which were all confirmed in the matrix-isolation FTIR spectra. There is also evidence for the production of substituted ketenes in pyrolysis, most likely hydroxyketene and methylketene.

09 BIOMASS FUELS↗
DOE OSTI2022

Revisiting the K-edge X-ray absorption fine structure of Si, Ge–Si alloys, and the isoelectronic series: CuBr, ZnSe, GaAs, and Ge

Extended X-ray absorption fine structure (EXAFS) has evolved into an unprecedented local-structure technique that is routinely used to study materials’ problems in the biological, chemical, and physical sciences. Like many other experimental techniques, EXAFS also requires that several key atomic parameters must be known a priori before structural information can be quantitatively determined. Utilizing current analytical methods, we revisit the isoelectronic series CuBr, ZnSe, GaAs, and Ge originally studied by Stern et al. during the early development of EXAFS. We demonstrate that the ab initio EXAFS code FEFF accurately predicts the atomic phase shifts and backscattering amplitudes that are primarily functions of the sum of atomic numbers Z along an EXAFS scattering path. We also investigate quantitative fitting and first- and second-shell phase transferability together with problems that arise if a backscattering atom is identified incorrectly in an EXAFS fitting model. Features in the near-edge region, on the other hand, are shown to require a comprehensive treatment of the band structure and density-of-states, including effects of the screened Coulomb interaction between the photoelectron and core hole. Here, we demonstrate that the Bethe–Salpeter equation (BSE) accurately captures the NEXAFS (or XANES) portion of the spectrum for the isoelectronic series in addition to Si and Ge–Si alloys, including within a few eV of the absorption edge, where band structure and excitonic effects are most important.

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

Ab initio equation of state for vanadium

Accurate materials’ equations of state (EOS) are essential for understanding materials properties as well as for use in multiphysics simulations. In particular, hydrodynamics simulations are based on three fundamental conservation laws (mass, momentum and energy) that form an under-determined system of equations. The equation of state serves as an additional closure relation between thermodynamic variables for a given material that enables numerical hydrodynamics simulation. In this report, we focus on the development of an ab initio EOS for the body centered cubic (BCC) phase of Vanadium (V) for eventual integration into a multiphase EOS in the OpenSesame EOS database.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Structure Factors for Hot Neutron Matter from Ab Initio Lattice Simulations with High-Fidelity Chiral Interactions

We present the first ab initio lattice calculations of spin and density correlations in hot neutron matter using high-fidelity interactions at next-to-next-to-next-to-leading order in chiral effective field theory. These correlations have a large impact on neutrino heating and shock revival in core-collapse supernovae and are encapsulated in functions called structure factors. Unfortunately, calculations of structure factors using high-fidelity chiral interactions were well out of reach using existing computational methods. In this Letter, we solve the problem using a computational approach called the rank-one operator (RO) method. The RO method is a general technique with broad applications to simulations of fermionic many-body systems. It solves the problem of exponential scaling of computational effort when using perturbation theory for higher-body operators and higher-order corrections. Using the RO method, we compute the vector and axial static structure factors for hot neutron matter as a function of temperature and density. Here, the ab initio lattice results are in good agreement with virial expansion calculations at low densities but are more reliable at higher densities. Random phase approximation codes used to estimate neutrino opacity in core-collapse supernovae simulations can now be calibrated with ab initio lattice calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

A Deep Potential model for liquid–vapor equilibrium and cavitation rates of water

Computational studies of liquid water and its phase transition into vapor have traditionally been performed using classical water models. Here, we utilize the Deep Potential methodology—a machine learning approach—to study this ubiquitous phase transition, starting from the phase diagram in the liquid–vapor coexistence regime. The machine learning model is trained on ab initio energies and forces based on the SCAN density functional, which has been previously shown to reproduce solid phases and other properties of water. Here, we compute the surface tension, saturation pressure, and enthalpy of vaporization for a range of temperatures spanning from 300 to 600 K and evaluate the Deep Potential model performance against experimental results and the semiempirical TIP4P/2005 classical model. Moreover, by employing the seeding technique, we evaluate the free energy barrier and nucleation rate at negative pressures for the isotherm of 296.4 K. Further, we find that the nucleation rates obtained from the Deep Potential model deviate from those computed for the TIP4P/2005 water model due to an underestimation in the surface tension from the Deep Potential model. From analysis of the seeding simulations, we also evaluate the Tolman length for the Deep Potential water model, which is (0.091 ± 0.008) nm at 296.4 K. Finally, we identify that water molecules display a preferential orientation in the liquid–vapor interface, in which H atoms tend to point toward the vapor phase to maximize the enthalpic gain of interfacial molecules. We find that this behavior is more pronounced for planar interfaces than for the curved interfaces in bubbles. This work represents the first application of Deep Potential models to the study of liquid–vapor coexistence and water cavitation.

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

Bi-directional tuning of thermal transport in SrCoO x with electrochemically induced phase transitions

Unlike the wide-ranging dynamic control of electrical conductivity, there does not exist an analogous ability to tune thermal conductivity by means of electric potential. The traditional picture assumes that atoms inserted into a material’s lattice act purely as a source of scattering for thermal carriers, which can only reduce thermal conductivity. In contrast, here we show that the electrochemical control of oxygen and proton concentration in an oxide provides a new ability to bi-directionally control thermal conductivity. On electrochemically oxygenating the brownmillerite SrCoO 2.5 to the perovskite SrCoO 3–δ , the thermal conductivity increases by a factor of 2.5, whereas protonating it to form hydrogenated SrCoO 2.5 effectively reduces the thermal conductivity by a factor of four. This bi-directional tuning of thermal conductivity across a nearly 10 ± 4-fold range at room temperature is achieved by using ionic liquid gating to trigger the ‘tri-state’ phase transitions in a single device. We elucidated the effects of these anionic and cationic species, and the resultant changes in lattice constants and lattice symmetry on thermal conductivity by combining chemical and structural information from X-ray absorption spectroscopy with thermoreflectance thermal conductivity measurements and ab initio calculations. This ability to control multiple ion types, multiple phase transitions and electronic conductivity that spans metallic through to insulating behaviour in oxides by electrical means provides a new framework for tuning thermal transport over a wide range.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modeling the thermodynamics of the FeTi hydrogenation under para-equilibrium: An ab-initio and experimental study

We report FeTi-based hydrides have recently re-attracted attention as stationary hydrogen storage materials due to favorable reversibility, good sorption kinetics and relatively low costs compared to alternative intermetallic hydrides. Employing the OpenCalphad software, the thermodynamics of the (FeTi) 1-x H x (0 ≤ x ≤ 1) system were assessed as a key basis for modeling hydrogenation of FeTi-based alloys. New thermodynamic data were acquired from our experimental pressure-composition-isotherm (PCI) curves, as well as first-principles calculations utilizing density functional theory (DFT). The thermodynamic phase models were carefully selected based on critical analysis of literature information and ab-initio investigations. Key thermodynamic properties such as dissociation pressure, formation enthalpies and phase diagrams were calculated in good agreement to our performed experiments and literature-reported data. This work provides an initial perspective, which can be extended to account for higher-order thermodynamic assessments and subsequently enables the design of novel FeTi-based hydrides. In addition, the assessed thermodynamic data can serve as key inputs for kinetic models and hydride microstructure simulations.

08 HYDROGEN↗
DOE OSTI2025

Ab Initio Molecular Dynamics Insights into Stress Corrosion Cracking and Dissolution of Metal Oxides

Oxide phases such as α-Fe2O3 (hematite) and α-Al2O3 (corundum) are highly insoluble in water; however, subcritical crack growth has been observed in humidity nonetheless. Chemically induced bond breaking at the crack tip appears unlikely due to sterically hindered molecular transport. The molecular mechanics of a crack in corundum with a reactive force field reveal minimal lattice trapping, leading to bond breaking before sufficient space opens for water transport. To address this, we model a pre-built blunt crack with space for H2O molecule adsorption at the tip and show that it reduces fracture toughness by lowering the critical J-integral. Then, we explore stress-enhanced dissolution to understand the mechanism of crack tip blunting in the oxide/water system. Density functional theory combined with metadynamics was employed to describe atomic dissolution from flat hematite and corundum surfaces in pure water. Strain accelerates dissolution, stabilizing intermediate states with broken bonds before full atom detachment, while the free energy profile of unstrained surfaces is almost monotonic. The atomistic calculations provided input for a kinetic model, predicting the shape evolution of a blunt crack tip, which displays three distinct regimes: (i) dissolution primarily away from the tip, (ii) enhanced blunting near but not at the apex, and (iii) sharpening near the apex. The transition between regimes occurs at a low strain, highlighting the critical role of water in the subcritical crack growth of oxide scales, with dissolution as the fundamental microscopic mechanism behind this process.

Chemistry↗
DOE OSTI2021

Gas-phase synthesis of corannulene – a molecular building block of fullerenes

Fullerenes (C 60 , C 70 ) detected in planetary nebulae and carbonaceous chondrites have been implicated to play a key role in the astrochemical evolution of the interstellar medium. However, the formation mechanism of even their simplest molecular building block—the corannulene molecule (C 20 H 10 )—has remained elusive. Here we demonstrate via a combined molecular beams and ab initio investigation that corannulene can be synthesized in the gas phase through the reactions of 7-fluoranthenyl (C 16 H 9 ˙) and benzo[ghi]fluoranthen-5-yl (C 18 H 9 ˙) radicals with acetylene (C 2 H 2 ) mimicking conditions in carbon-rich circumstellar envelopes. This reaction sequence reveals a reaction class in which a polycyclic aromatic hydrocarbon (PAH) radical undergoes ring expansion while simultaneously forming an out-of-plane carbon backbone central to 3D nanostructures such as buckybowls and buckyballs. In conclusion, these fundamental reaction mechanisms are critical in facilitating an intimate understanding of the origin and evolution of the molecular universe and, in particular, of carbon in our galaxy.

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

Generation of the TSL for Zirconium Hydrides from Ab Initio Methods

Zirconium hydride (ZrHx) is a moderator material used in TRIGA and other reactors that may exist in multiple phases with varying stoichiometry, which include the δ phase and the ϵ phase. Current ENDF/B-VIII.0 ZrHx thermal scattering law (TSL) evaluations do not distinguish between phases. These sub-libraries were generated with the LEAPR module of NJOY using historic phonon spectra derived from a central force model and assume incoherent elastic scattering for both bound hydrogen and zirconium, which neglects the effects of crystal structures important for scattering from zirconium bound in ZrHx. In this work, the TSLs for hydrogen and zirconium bound in δ-ZrHx and ϵ-ZrH2 were generated from phonon spectra derived from modern ab initio lattice dynamics methods and ab initio molecular dynamics. Subsequently, TSLs for hydrogen and zirconium in ZrHx and ZrH2 were generated using the Full Law Analysis Scattering System Hub (FLASSH) code. The built-in generalized coherent elastic routine was used to generate the previously neglected elastic contribution from zirconium for this material. The present TSLs provide both a re-evaluation of the current ZrH sub-libraries and expansion of the set of TSLs available for the examination of neutrons in systems with zirconium hydride, permitting explicit treatment of δ and ϵ phases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Role of Coulomb interaction in the phase formation of fcc Ce: Correlation matrix renormalization theory

The effect of electronic Coulomb interaction on the phase formation of fcc Ce lattice is investigated by full ab initio calculations without adjustable Coulomb U and J parameters using the Gutzwiller wavefunction-based correlation matrix renormalization theory (CMRT). Its total energy and pressure as a function of volume agree reasonably well with existing DFT+Gutzwiller calculations and experiments, indicating correct capture of electronic correlation and screening effects within the CMRT formalism. Here, a stable phase is found in line with the experimental α-Ce phase, and a lurking phase is identified supposedly linked with the experimental γ-Ce phase. A criterion based on the local 4f electron charge fluctuation is introduced to confirm the distinct electronic correlation natures of both phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗