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

Chemical short-range order in derivative Cr–Ta–Ti–V–W high entropy alloys from the first-principles thermodynamic study

The development of high-entropy alloys (HEAs) focuses on exploring compositional regions in multi-component systems with all alloy elements in equal or near-equal atomic concentrations. Initially it was based on the main idea that high mixing configurational entropy contributions to the alloy free energy could promote the formation of a single solid solution phase. By using the ab-initio based Cluster Expansion (CE) Hamiltonian model constructed for the quinary bcc Cr–Ta–Ti–V–W system in combination with Monte Carlo (MC) simulations, we show that the phase stability and chemical short-range order (SRO) of the equiatomic quinary and five sub-quaternary systems, as well as their derivative alloys, can dramatically change the order–disorder transition temperatures (ODTT) as a function of alloy compositions. In particular, it has been found, that the equiatomic quaternary Ta–Ti–V–W and Cr–Ta–Ti–W alloys had the lowest order–disorder transition temperature (500 K) among all the analysed equiatomic compositions. In all investigated alloy systems, the strongest chemical ordering has been observed between Cr and V, which led to the conclusion that decreasing the concentration of either Cr or V might be beneficial in terms of decreasing the ODTT. It also predicts that increasing concentration of Ti significantly decreases the ODTT. Our analysis of chemical SRO as a function of alloy composition allows to understand the microstructure evolution of HEAs as a function of temperature in excellent agreement with available experimental observations. Importantly, our free energy of mixing and SRO calculations predict that the origin of precipitates formed by Cr- and V-rich in the sub-quaternary Cr–Ta–V–W system is driven by the thermodynamics. The modelling results are in an excellent agreement with experimental observation of Cr and V segregation in the W 0.38 Ta 0.36 Cr 0.15 V 0.11 alloy which in turns shows an exceptional radiation resistance.

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

Floquet band engineering and topological phase transitions in 1T’ transition metal dichalcogenides

Using ab initio tight-binding approaches, we investigate Floquet band engineering of the 1T' phase of transition metal dichalcogenides (MX 2 , M = W, Mo and X = Te, Se, S) monolayers under the irradiation with circularly polarized light. Furthermore, our first principles calculations demonstrate that light can induce important transitions in the topological phases of this emerging materials family. For example, upon irradiation, Te-based MX2 undergoes a phase transition from quantum spin Hall (QSH) semimetal to time-reversal symmetry broken QSH insulator with a nontrivial band gap of up to 92.5 meV. On the other hand, Se- and S-based MX 2 undergoes the topological phase transition from the QSH effect to the quantum anomalous Hall effect and into trivial phases with increasing light intensity. From a general perspective, this theoretical work brings further insight into non-equilibrium topological systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Phase transitions and potential ferroelectricity in noncentrosymmetric KNaNbOF 5

Increasing the number of known ferroelectrics requires expanding the chemical design space and mechanisms producing ferroelectricity. To that end, we examine the displacive, noncentrosymmetric-to-centrosymmetric phase transition in the oxyfluoride KNaNbOF 5 using ab initio calculations and Landau theory. Here we predict an intermediate Pnma phase occurs in the transition from the known low- and high-temperature phases, Pna2 1 and Cmcm, respectively, which we subsequently confirm using high-resolution and in situ 19 F magic-angle spinning nuclear magnetic resonance. Using the intermediate structure, we then calculate the monodomain polarization switching barrier in KNaNbOF 5 to be ~ 93 meV per formula unit, comparable to existing ferroelectrics. The reversal of the polarization is due to uncompensated antipolar displacements of oxygen and fluorine and does not require oxygen and fluorine site exchange as in other nonswitchable heteroanionic materials, which makes KNaNbOF 5 a viable oxyfluoride for a switchable electric polarization.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Structural and magnetic transitions of CoFeMnNiAl high-entropy alloys caused by composition and annealing

Tailoring magnetic behavior by composition and annealing is an effective way. The phase transition from the face-centered-cubic (FCC) to B2 structure influenced by the Al content on the CoFeMnNiAl x alloys leads to the enhanced magnetization. However, the excessive Al addition in the B2-structured alloys impairs the magnetization. The phase transformation caused by annealing also indicates that the saturation magnetization of the alloy depends on the volume fraction of the B2 phase, and a higher content of the B2 phase is beneficial to the large saturation magnetization. Ab initio calculations are used to explain the magnetic behavior of the present HEAs and the effects of phase structures on magnetic characteristics. The phase stability of the CoFeMnNiAl x alloy is thoroughly studied by both the existing phase-formation empirical criteria and annealing. Here, the existing empirical criteria used to predict the phase formation are not universal. However, combining the ΔH mix -δ relation and the φ value, whether the single disordered solid-solution phase is formed can be preliminarily predicted.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Liquid-Vapor Coexistence and Critical Point of Mg 2 SiO 4 From Ab Initio Simulations

Hypervelocity impact-driven vaporization is characteristic of late-stage planet formation. Yet the behavior and properties of liquid-vapor mixtures of planetary materials of interest are typically unknown. Multiphase equations of state used in hydrodynamic simulations of planet impacts therefore lack reliable data for this important phenomenon. Here, we present the first constraints on the liquid-vapor critical point and coexistence phase boundary of Mg2SiO4 computed from ab initio molecular dynamics simulations. We found that the vapor is depleted in magnesium and enriched in silica and oxygen, while the coexisting liquid is enriched in magnesium and depleted in oxygen, from which we infer vaporization is incongruent. The critical point was estimated from an equation of state fit to the data. The results are in line with recent calculations of MgSiO3 and together confirm that extant multiphase equation of state (EOS) models used in planetary accretion modeling significantly underestimate the amount of supercritical material postimpact.

58 GEOSCIENCES↗
DOE OSTI2023

Ab Initio Simulations of Tritium Diffusion in Al-Rich Iron Aluminide Coating Phases

Surface coatings of steels used in extreme conditions and corrosive environments generally aim to provide protection and increased durability. In the case of tritium-producing burnable absorber rods (TPBARs) used in nuclear reactors, a 316 stainless steel has been coated with Al. Scanning transmission electron microscopy (STEM) characterization of the coating found three Al-rich (>60 atom % Al) iron aluminide alloys identified as hexagonal FeNiAl 5 , monoclinic Fe 4 Al 13 , and orthorhombic Fe 2 Al 5 . Density functional theory simulations using nudged elastic band have been performed to investigate the diffusion of interstitial tritium in each Al-rich iron aluminide phase. While FeNiAl 5 and Fe 4 Al 13 can be viewed as the stacking of two layers, the structural peculiarity of Fe 2 Al 5 is that channels of variable Al vacancy content are present along the c-axis. Therefore, three stoichiometries for Fe 2 Al x phase, namely, Fe 2 Al 4 , Fe 2 Al 5 , and Fe 2 Al 6 , have been considered to evaluate the impact of Al vacancy concentration on tritium diffusion behavior. Altogether, we found that at 600 K, tritium diffusion decreases from a faster rate in the channels of Fe 2 Al x phases (D T ≤ 10 –11 m 2 ·s –1 ) to Fe 4 Al 13 (D T ≈ 10–12 m2·s–1), and finally in FeNiAl5 (D T ≈ 10 –13 m 2 ·s –1 ). Here, we also find that interstitial tritium generally diffuses faster in Fe–Al coating phases than in the tritium breeding material γ-LiAlO 2 (D T ≈ 10 –14 m 2 ·s –1 ) but slightly slower than in 316 stainless steel (D T ≈ 10 –10 m 2 ·s –1 ).

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

Condensed Phase Water Molecular Multipole Moments from Deep Neural Network Models Trained on Ab Initio Simulation Data

Ionic solvation phenomena in liquids involve intense interactions in the inner solvation shell. For interactions beyond the first shell, the ion–solvent interaction energies result from the sum of many smaller-magnitude contributions that can still include polarization effects. Deep neural network (DNN) methods have recently found wide application in developing efficient molecular models that maintain near-quantum accuracy. In this work, we extend the DeePMD-kit code to produce accurate molecular multipole moments in the bulk and near interfaces. The new method is validated by comparing the DNN moments with those generated by ab initio simulations. The moments are used to compute the electrostatic potential at the center of a molecular-sized hydrophobic cavity in water. The results show that the fields produced by the DNN models are in quantitative agreement with the AIMD-derived values. These efficient methods will open the door to more accurate solvation models for large solutes such as proteins.

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

Structures of elemental potassium at terapascal pressures

Here we investigate the pressure-temperature phase diagram of elemental potassium (K) up to multiterapascal (TPa) pressures using ab initio random structure searching (AIRSS), discovering eight structural phase transitions beyond the already known double hexagonal close-packed (dhcp) structure. Starting at 1.15 TPa, K transitions from the dhcp structure and passes through a variety of close-packed hexagonal and trigonal phases (dhcp → P6 3 /mmc → P3m1 → P$\bar{6}m2$), which differ only in their stacking sequence along the c axis. At 2.55 TPa, K adopts the hexagonal close-packed (hcp) structure, and at 7.21 TPa transitions into a complex orthorhombic phase of Fddd symmetry with 64 atoms in its (conventional) unit cell, before assuming an Ibam structure at 8.60 TPa, and, eventually, the face centered cubic (fcc) structure at 23.6 TPa, which persists well into the petapascal (PPa) regime. Further, we calculate the full pressure-temperature dependence of the melting line of K through extensive molecular dynamics simulations. We study the evolution of the bonding topology of K with pressure, finding that K passes through two quasi-molecular phases featuring diatomic pairs (the Fddd and Ibam structures), before ultimately becoming a high pressure electride (HPE) in the fcc phase. The electron-phonon and superconducting properties of K at these extreme pressures are also investigated, where we find the critical temperature T c rises to a maximum of between 7.65 K and 15.70 K in the hcp phase, before eventually decreasing to essentially zero in the fcc phase. Our results fully elucidate the structural and electronic behavior of K under the most extreme conditions, and provide new case studies for multi-TPa dynamic compression experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Ab Initio Calculations on the Ground and Excited Electronic States of Thorium–Ammonia, Thorium–Aza-Crown, and Thorium–Crown Ether Complexes

Positively charged metal–ammonia complexes are known to host peripheral, diffuse electrons around their molecular skeleton. The resulting neutral species form materials known as expanded or liquid metals. Alkali, alkaline earth, and transition metals have been investigated previously in experimental and theoretical studies of both the gas and condensed phase. This work is the first ab initio exploration of an f-block metal–ammonia complex. The ground and excited states are calculated for Th0–3+ complexes with ammonia, crown ethers, and aza-crown ethers. For Th3+ complexes, the one valence electron Th populates the metal’s 6d or 7f orbitals. For Th0–2+, the additional electrons prefer occupation of the outer s- and p-type orbitals of the complex, except Th(NH3)10, which uniquely places all four electrons in outer orbitals of the complex. Although thorium coordinates up to ten ammonia ligands, octa-coordinated complexes are more stable. Crown ether complexes have a similar electronic spectrum to ammonia complexes, but excitations of electrons in the outer orbitals of the complex are higher in energy. Aza-crown ethers disfavor the orbitals perpendicular to the crowns, attributed to the N-H bonds pointing along the plane of the crowns.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

High-Temperature Thermodynamics of Uranium from Ab Initio Modeling

We present high-temperature thermodynamic properties for uranium in its γ phase (γ-U) from first-principles, relativistic, and anharmonic theory. The results are compared to CALPHAD modeling. The ab initio electronic structure is obtained from density-functional theory (DFT) that includes spin–orbit coupling and an added self-consistent orbital-polarization (OP) mechanism for more accurate treatment of magnetism. The first-principles method is coupled to a lattice dynamics scheme that is used to model anharmonic lattice vibrations, namely, Self-Consistent Ab Initio Lattice Dynamics (SCAILD). The methodology can be summarized in the acronym DFT + OP + SCAILD. Upon thermal expansion, γ-U develops non-negligible magnetic moments that are included for the first time in thermodynamic theory. The all-electron DFT approach is shown to model γ-U better than the commonly used pseudopotential method. In addition to CALPHAD, DFT + OP + SCAILD thermodynamic properties are compared with other ab initio and semiempirical modeling and experiments. Our first-principles approach produces Gibbs free energy that is essentially identical to CALPHAD. The DFT + OP + SCAILD heat capacity is close to CALPHAD and most experimental data and is predicted to have a significant thermal dependence due to the electronic contribution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Ab initio simulations on the pure Cr lattice stability at 0K: Verification with the Fe-Cr and Ni-Cr binary systems

Significant discrepancies have been observed and discussed on the lattice stability of Cr between the predictions from the ab initio calculations and the CALPHAD approach. In the current work, we carefully examined the possible structures for pure Cr and reviewed the history back from how Kaufman originally determined the Gibbs energy of FCC-Cr in the 1970s. The reliability of Cr lattice stability derived by the CALPHAD and ab initio approaches was systematically discussed. It is concluded that the Cr lattice stability based on the CALPHAD approach has large uncertainty. Meanwhile, we cannot claim that the ab initio H FCC-Cr is error-free as FCC-Cr is an unstable phase under ambient conditions. The present work shows that the ab initio HFCC-Cr can be a viable scientific approach. As both approaches have their limitations, the present work propose to integrate the ab initio results into the CALPHAD platform for the development of the next generation CALPHAD database. As a result, the Fe-Cr and Ni-Cr binary systems were chosen as two case studies demonstrating the capability to adopt the ab initio Cr lattice stability directly into the current CALPHAD database framework.

36 MATERIALS SCIENCE↗
DOE OSTI2020

First-principles study of the low-temperature charge density wave phase in the quasi-one-dimensional Weyl chiral compound ( TaSe 4 ) 2 I

Using ab initio density functional theory, we study herein the lattice phase transition of quasi-one-dimensional (TaSe 4 ) 2 I. In the undistorted state, the strongly anisotropic semimetal band structure presents two nonequivalent Weyl points. In previous efforts, two possible Ta-tetramerization patterns were proposed to be associated with the low-temperature structure. Our phonon calculations indicate that the orthorhombic F222 CDW-I phase is the most likely ground state for this quasi-one-dimensional system. In addition, the monoclinic C2 CDW-II phase may also be stable according to the phonon dispersion spectrum. Since these two phases have very similar energies in our density functional theory calculations, both these Ta-tetramerization distortions likely compete or coexist at low temperatures. The semimetal-to-insulator transition is induced by a Fermi-surface-driven instability that supports the Peierls scenario, which affects the Weyl physics developed above T CDW . Furthermore, the spin-orbit coupling generates Rashba-like band splittings in the insulating charge density wave phases.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Visualizing interfacial collective reaction behaviour of Li–S batteries

Benefiting from high energy density (2,600 Wh kg -1 )) and low cost, lithium-sulfur (Li-S) batteries are considered promising candidates for advanced energy-storage systems. Despite tremendous efforts in suppressing the long-standing shuttle effect of lithium polysulfides, understanding of the interfacial reactions of lithium polysulfides at the nanoscale remains elusive. This is mainly because of the limitations of in situ characterization tools in tracing the liquid-solid conversion of unstable lithium polysulfides at high temporal-spatial resolution. There is an urgent need to understand the coupled phenomena inside Li-S batteries, specifically, the dynamic distribution, aggregation, deposition and dissolution of lithium polysulfides. Here, by using in situ liquid-cell electrochemical transmission electron microscopy, we directly visualized the transformation of lithium polysulfides over electrode surfaces at the atomic scale. Notably, an unexpected gathering-induced collective charge transfer of lithium polysulfides was captured on the nanocluster active-centre-immobilized surface. It further induced an instantaneous deposition of nonequilibrium Li2S nanocrystals from the dense liquid phase of lithium polysulfides. Without mediation of active centres, the reactions followed a classical single-molecule pathway, lithium polysulfides transforming into Li 2 S 2 and Li 2 S step by step. Molecular dynamics simulations indicated that the long-range electrostatic interaction between active centres and lithium polysulfides promoted the formation of a dense phase consisting of Li + and S n 2- (2 < n ≤ 6), and the collective charge transfer in the dense phase was further verified by ab initio molecular dynamics simulations. Finally, the collective interfacial reaction pathway unveils a new transformation mechanism and deepens the fundamental understanding of Li-S batteries.

25 ENERGY STORAGE↗
DOE OSTI2021

Ground State Properties of the Wide Band Gap Semiconductor Beryllium Sulfide (BeS)

We report the results from self-consistent calculations of electronic, transport, and bulk properties of beryllium sulfide (BeS) in the zinc-blende phase, and employed an ab-initio local density approximation (LDA) potential and the linear combination of atomic orbitals (LCAO). We obtained the ground state properties of zb-BeS with the Bagayoko, Zhao, and Williams (BZW) computational method, as enhanced by Ekuma and Franklin (BZW-EF). Our findings include the electronic energy bands, the total (DOS) and partial (pDOS) densities of states, electron and hole effective masses, the equilibrium lattice constant, and the bulk modulus. The calculated band structure clearly shows that zb-BeS has an indirect energy band gap of 5.436 eV, from Γ to a point between Γ and X, for an experimental lattice constant of 4.863 Å. This is in excellent agreement with the experiment, unlike the findings of more than 15 previous density functional theory (DFT) calculations that did not perform the generalized minimization of the energy functional, required by the second DFT theorem, which is inherent to the implementation of our BZW-EF method.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Superconducting properties of MoTe 2 from ab initio anisotropic Migdal-Eliashberg theory

Molybdenum ditelluride (MoTe 2 ) is attracting considerable interest since it is the archetypal type-II Weyl semimetal and a candidate for topological superconductivity. In this work, we investigate the superconducting phase diagram of two MoTe 2 polymorphs using the ab initio anisotropic Migdal-Eliashberg theory, and we show that the superconducting dome originates from the synergistic contribution of the density of states at the Fermi level and the transverse acoustic Te modes in the 1$\textit{T}’$ phase. We find that the electron and hole pockets carry trivial $\textit{s}$-wave order parameters of slightly different magnitude, reminiscent of a two-gap structure as suggested by recent experiments. We suggest that a possible route for enhancing the superconducting critical temperature, and realizing $s_{+ –}$ pairing, in the $T_d$ phase is to exploit its nontrivial band topology via electron doping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Reversible hydrogen storage in multilayer graphane: Lattice dynamics, compressibility, and heat capacity studies

Multilayer graphane (hydride of graphite) is a crystalline hydrocarbon of composition CH, which can be synthesized from graphite and molecular hydrogen at pressures above 2GPa [V.E. Antonov et al. Carbon 100 (2016) 465]. Using X-ray diffraction, this compound was tentatively identified as the “graphane II” phase of 3D-graphane predicted by ab initio calculations [X.-D. Wen et al. PNAS 108 (2011) 6833] and consisting of layers of 2D-graphane in the “chair” conformation. When heated in a vacuum, the compound does not form any intermediate hydrocarbons and reversibly decomposes back into graphite and hydrogen at 770–920 K. In the present work, almost single-phase samples of graphite hydride and deuteride were synthesized at 7.4 GPa and 870 K. Their investigation by inelastic neutron scattering supplemented by ab initio calculations gave spectra g(E) of the phonon density of states with a gap of about 15 meV at approx. 100 meV, which is a unique identifier for the chair form of graphane. The equation of state V(P) of the hydride was studied at room temperature and hydrogen pressures up to 53 GPa by synchrotron X-ray diffraction in a diamond anvil cell. Further, the graphane II phase did not react with the surrounding hydrogen and did not undergo any phase transformations upon the compression and after heating to 1500 K at 53GPa. The high thermal and pressure stability of this exotic phase makes it an important part of the C–H system. The obtained g(E) spectra of graphite hydride and deuteride were used to calculate temperature dependences of their heat capacity. Measurements of the heat capacity at temperatures 120–673 K confirmed the good accuracy of these calculations.

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

Observed metallization of hydrogen interpreted as a band structure effect

A recent experimental study of the metallization of hydrogen tracked the direct band gap and vibron frequency via infrared measurements up to ~425 GPa (Loubeyre et al (2020 Nature 577 631). Above this pressure, the direct gap has a discontinuous drop to below the minimum experimentally accessible energy (~0.1 eV). The authors suggested that this observation is caused by a structural phase transition between the C2/c-24 molecular phase to another molecular phase such as Cmca-12. Here, through ab initio calculations of pressure dependent vibron frequency and direct band gap, we find that the experimental data is consistent with the C2/c-24 phase up to 425 GPa, and suggest that this consistency extends beyond that pressure. Specifically, we find that qualitative changes in the band structure of the C2/c-24 phase lead to a discontinuous drop of the direct band gap, which can explain the observed drop without a structural transition. In conclusion, this alternative scenario, which naturally explains the absence of hysteresis in the measurements, will hopefully motivate further experimental studies to ascertain the structure of the phase above the high pressure 'phase transition'.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Atomic-Scale Dynamic Interaction of H 2 O Molecules with Cu Surface

Atomic-scale interaction of water vapor with metal surfaces beyond surface adsorption under technologically relevant conditions remains mostly unexplored. Using aberration-corrected environmental transmission electron microscopy, we reveal the dynamic surface activation of Cu by H 2 O at elevated temperature and pressure. In this study, we find a structural transition from flat to corrugated surface for the Cu(011) under low water-vapor pressure. Increasing the water-vapor pressure leads to the surface reaction of Cu with dissociated H 2 O, resulting in the formation of a metastable “bilayer” Cu-O-H phase. Corroborated by density functional theory and ab initio molecular dynamics calculations, the cooperative O and OH interaction with Cu is responsible for the formation and subsurface propagation of this phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗