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 289 records · Page 16

Prediction of L⁢i 3 ⁢F⁡e 8 ⁢B 8 compound with rapid one-dimensional ion diffusion channels

Using a computational crystal structure search in the Li-Fe-B ternary system, we predict a stable phase of L⁢i 3 ⁢F⁡e 8 ⁢B 8 , featuring 1D channels that enable rapid Li-ion transport. Ab initio molecular dynamics simulations show that the Li-ion diffusion coefficient in L⁢i 3 ⁢F⁡e 8⁢ B 8 surpasses that of common electrode and conductive additive materials by several orders of magnitude. The high diffusion in L⁢i 3 ⁢F⁡e 8 ⁢B 8 can be explained by the Frenkel–Kontorova model, which describes an incommensurate state between the Li diffusion chain and the periodic potential field caused by the FeB backbone structure. The favorable lithium-ion diffusivity and mechanical properties of L⁢i 3 ⁢F⁡e 8 ⁢B 8 make it a promising conductive additive for battery materials. Furthermore, an external magnetic field can further manipulate the properties of this material due to its predicted itinerant ferromagnetism, which also offers a platform for exploring spin-dependent phenomena.

1-dimensional systems↗

First-principles simulation of a shocked H-He mixture along the principal Hugoniot

Recent laser-shock experiments on an H–He mixture containing 11% helium (atomic fraction) have suggested the presence of an immiscibility region inside Jupiter. Reflectivity measurements were used as the primary diagnostic of H–He demixing, with discontinuities in the optical reflectivity proposed as a signature of phase separation under conditions relevant to Jupiter's interior. Here, we investigate shock-compressed H–He using ab initio molecular dynamics simulations with optical properties evaluated within the Kubo–Greenwood formalism. The equation of state and ionic configurations were obtained using the thermal Tr 2 SCANL meta-GGA exchange–correlation (XC) functional, while optical properties were computed using the recently developed RS-KDT0 range-separated thermal hybrid XC, which provides state-of-the-art accuracy for band-gap predictions in the warm dense matter regime. The calculated reflectivity shows overall good agreement with experimental measurements; however, no discontinuity is observed at elevated temperatures. Moreover, the reflectivity predictions for the mixed system are consistent with the experimental measurements in the temperature range where the mixture is inferred to be demixed. Furthermore, these results suggest that reflectivity alone may not provide a unique or sensitive diagnostic of H–He demixing at low helium concentrations under these conditions.

Density functional calculations↗

Glide-resolved photoemission spectroscopy: Measuring topological invariants in nonsymmorphic space groups

The two classes of 3D, time-reversal-invariant insulators are known to subdivide into four classes in the presence of glide symmetry. Here, we extend this classification of insulators to include glide-symmetric Weyl metals and find a finer $\mathbb{Z}_4 ⊕\mathbb {Z}$ classification. We further elucidate the smoking-gun experimental signature of each class in the photoemission spectroscopy of surface states. Measuring the $\mathbb{Z}_4$ topological invariant by photoemission relies on identifying the glide representation of the initial Bloch state before photoexcitation—we show how this is accomplished with relativistic selection rules, combined with standard spectroscopic techniques to resolve both momentum and spin. In this work, our method relies on a spin-momentum locking that is characteristic of all glide-symmetric solids (inclusive of insulators and metals in trivial and topological categories). As an orthogonal application, given a glide-symmetric solid with an ideally symmetric surface, we may utilize this spin-momentum locking to generate a source of fully spin-polarized photoelectrons, which have diverse applications in solid-state spectroscopy. Our ab initio calculations predict Ba 2 Pb, stressed Na 3 Bi, and KHgSb to realize all three, nontrivial insulating phases in the $\mathbb{Z}_4$ classification.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Calculations concerning the reaction C + H3(+) - CH(+) + H2

The ion-molecule reaction of C + H3(+) - CH(+) + H2 has not been studied in the laboratory but is thought to be important in the gas phase synthesis of organic molecules in dense interstellar clouds. This reaction was studied theoretically by performing quantum chemical ab initio calculations on the potential surface. It is found that there is no activation barrier to the reaction and that it proceeds smoothly to the first excited electronic state of CH(+). The rate coefficient as a function of temperature can then be estimated using the proper long-range potentials. The rate coefficient at 10 K is calculated to be 2.9 x 10 to the -9th cu cm/sec.

Talbi, D.↗

Toward an accurate equation of state and B1-B2 phase boundary for magnesium oxide up to terapascal pressures and electron-volt temperatures

By applying auxiliary-field quantum Monte Carlo, we calculate the equation of state (EOS) and B1-B2 phase transition of magnesium oxide (MgO) up to 1 TPa. Here, the results agree with available experimental data at low pressures and are used to benchmark the performance of various exchange-correlation functionals in density functional theory calculations. We determine PBEsol is an optimal choice for the exchange-correlation functional and perform extensive phonon and quantum molecular-dynamics calculations to obtain the thermal EOS. Our results provide a preliminary reference for the EOS and B1-B2 phase boundary of MgO from zero up to 10,500 K.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Electronic, magnetic, and structural properties of CoVMnSb: Ab initio study

We present computational results on electronic, magnetic, and structural properties of CoVMnSb, a quaternary Heusler alloy. Our calculations indicate that this material may crystallize in two energetically close structural phases: inverted and regular cubic. The inverted cubic phase is the ground state, with ferrimagnetic alignment, and around 80% spin polarization. Despite having a relatively large bandgap in the minority-spin channel close to the Fermi level, this phase does not undergo a half-metallic transition under pressure. This is explained by the “pinning” of the Fermi level at the minority-spin states at the Γ point. At the same time, the regular cubic phase is half-metallic and retains its perfect spin polarization under a wide range of mechanical strain. Transition to a regular cubic phase may be attained by applying uniform pressure (but not biaxial strain). In practice, this pressure may be realized by an atomic substitution of non-magnetic atoms (Sb) with another non-magnetic atom (Si) of a smaller radius. Our calculations indicate that 25% substitution of Sb with Si results in a half-metallic regular cubic phase being the ground state. In addition, CoVMnSb0.5Si0.5 retains its half-metallic properties under a considerable range of mechanical pressure, as well as exhibits thermodynamic stability, thus making this alloy attractive for potential spintronic applications. We hope that the presented results will stimulate experimental efforts to synthesize this compound.

Physics↗

Ab Initio Transport Calculations: From Normal to Superconducting Current

Applying the Bogoliubov-de Gennes equations with density-functional theory, it is possible to formulate first principles description of current-phase relationships in superconducting/normal (magnetic)/superconducting tri-layers. Such structures are the basis for the superconducting analog of Magnetoresistive random access memory devices (JMRAM). In a recent paper1 we presented results from the first attempt to formulate such a theory, applied to the Nb/Ni/Nb trilayers. In the present work we provide computational details, explaining how to construct key ingredient (scattering matrices SN ) in a framework of linear muffin-tin orbitals (LMTO).

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗

Ab initio Quantum Chemical Studies of Reactions in Astrophysical Ices. Reactions Involving CH3OH, CO2, CO, HNCO in H2CO/NH3/H2O Ices

While reactions between closed shell molecules generally involve prohibitive barriers in the gas phase, prior experimental and theoretical studies have demonstrated that some of these reactions are significantly enhanced when confined within an icy grain mantle and can occur efficiently at temperatures below 100 K with no additional energy processing. The archetypal case is the reaction of formaldehyde (H2CO) and ammonia (NH3) to yield hydroxymethylamine (NH2CH2OH). In the present work we have characterized reactions involving methanol (CH3OH), carbon dioxide (CO2), carbon monoxide (CO), and isocyanic acid (HNCO) in search of other favorable cases. Most of the emphasis is on CH3OH, which was investigated in the two-body reaction with one H2CO and the three-body reaction with two H2CO molecules. The addition of a second H2CO to the product of the reaction between CH3OH and H2CO was also considered as an alternative route to longer polyoxymethylene polymers of the -CH2O- form. The reaction between HNCO and NH3 was studied to determine if it can compete against the barrierless charge transfer process that yields OCN(-) and NH4(+). Finally, the H2CO + NH3 reaction was revisited with additional benchmark calculations that confirm that little or no barrier is present when it occurs in ice.

Woon, David E.↗

Design of a V–Ti–Ni alloy with superelastic nano-precipitates

Stress-induced martensitic transformations enable metastable alloys to exhibit enhanced strain hardening capacity, leading to improved formability and toughness. As is well-known from transformation-induced plasticity (TRIP) steels, however, the resulting martensite can limit ductility and fatigue life due to its intrinsic brittleness. In this work, we explore an alloy design strategy that utilizes stress-induced martensitic transformations but does not retain the martensite phase. This strategy is based on the introduction of superelastic nano-precipitates, which exhibit reverse transformation after initial stress-induced forward transformation. To this end, utilizing ab-initio simulations and thermodynamic calculations we designed and produced a V 45 Ti 30 Ni 25 (at%) alloy. In this alloy, TiNi is present as nano-precipitates uniformly distributed within a ductile V-rich base-centered cubic (bcc) β matrix, as well as being present as a larger matrix phase. We characterized the microstructure of the produced alloy using various scanning electron microscopy (SEM) and transmission electron microscopy (TEM) methods. Furthermore, the bulk mechanical properties of the alloy are demonstrated through tensile tests, and the reversible transformation in each of the TiNi morphologies were confirmed by in-situ TEM micro-pillar compression experiments, in-situ high-energy diffraction synchrotron cyclic tensile tests, indentation experiments, and differential scanning calorimetry experiments. The observed transformation pathways and variables impacting phase stability are critically discussed.

59 BASIC BIOLOGICAL SCIENCES↗

Anharmonicity and Phase Diagram of Magnesium Oxide in the Megabar Regime

With density functional molecular dynamics simulations, we computed the phase diagram of MgO from 50 to 2000 GPa up to 20000 K. Via thermodynamic integration (TDI), we derive the Gibbs free energies of the B1, B2, and liquid phases and determine their phase boundaries. With TDI and a pseudo-quasi-harmonic approach, we show that anharmonic effects are important and stabilize the B1 phase in particular. As a result, the B1-B2 transition boundary in the pressure-temperature plane exhibits a steep slope. We predict the B1-B2-liquid triple point to occur at approximately T= 10000 K and P = 370 GPa, which is higher in pressure than was inferred with quasi-harmonic methods alone. We predict the principal shock Hugoniot curve to enter the B2 phase stability domain but only over a very small range of parameters. Here, this may render it difficult to observe this phase with shock experiments because of kinetic effects.

3-dimensional systems↗

Multi-reward Reinforcement Learning Based Bond-Order Potential to Study Strain-Assisted Phase Transitions in Phosphorene

Here, we introduce a multi-reward reinforcement learning (RL) approach to train a flexible bond-order potential (BOP) for 2D phosphorene based on ab initio training data sets. Our approach is based on a continuous action space Monte Carlo tree search algorithm that is general and scalable and presents an efficient multiobjective optimization scheme for high-dimensional materials design problems. As a proof-of-concept, we deploy this scheme to parametrize multiple structural and dynamical properties of 2D phosphorene polymorphs. Our RL-trained BOP model adequately captures the structure, energetics, transformation barriers, equation of state, elastic constants, and phonon dispersions of various 2D P polymorphs. We use this model to probe the impact of temperature and strain rate on the phase transition from black (α-P) to blue phosphorene (β-P) through molecular dynamics simulations. A decrease in critical strain for this phase transition with increase in temperature is observed, and the underlying atomistic mechanisms are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bespoke Liquid/Liquid Interfaces (Final Technical Report)

The goal of DE-SC0001815 was to advance the basic science of liquid:liquid interface formation, to develop a deeper understanding of the mechanisms of phase separation and the essential relationships between solution composition, organization and dynamics that underlie the kinetic regime of solvent extraction. This included learning how interfacial organization and dynamics alters the properties of the primary coordination sphere of ions and the free energy of transport of ions complexes across a phase boundary. We relied primarily upon classical molecular dynamics studies to determine the equilibrium ensembles of these complex systems, but also utilized ab-initio MD and cluster-based density functional theory (DFT) calculations when more detailed investigation of the electronic structure was needed. We continued development of graph-theory based analyses to elucidate hierarchical correlations and expanded into geometric topology methods to quantify the collectively organized structures that can organize at a liquid/liquid interface during solute transport. One of the main conclusions was from the observation of two distinct mechanisms for solute transport - those that derive from amplifications of interfacial heterogeneity and surface roughness, and those wherein surface roughness has been dampened and instead collectively organized macrostructures work to bring solutes into the organic phase. It was our aim to create a concrete chemical model of the underlying driving forces behind interfacial primary and secondary structure formation and to map out the energetic features of solute transport so that tailored liquid/liquid can be developed that have characteristic kinetic features associated with mass transport.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Exploring the frontiers of condensed-phase chemistry with a general reactive machine learning potential

Abstract Atomistic simulation has a broad range of applications from drug design to materials discovery. Machine learning interatomic potentials (MLIPs) have become an efficient alternative to computationally expensive ab initio simulations. For this reason, chemistry and materials science would greatly benefit from a general reactive MLIP, that is, an MLIP that is applicable to a broad range of reactive chemistry without the need for refitting. Here we develop a general reactive MLIP (ANI-1xnr) through automated sampling of condensed-phase reactions. ANI-1xnr is then applied to study five distinct systems: carbon solid-phase nucleation, graphene ring formation from acetylene, biofuel additives, combustion of methane and the spontaneous formation of glycine from early earth small molecules. In all studies, ANI-1xnr closely matches experiment (when available) and/or previous studies using traditional model chemistry methods. As such, ANI-1xnr proves to be a highly general reactive MLIP for C, H, N and O elements in the condensed phase, enabling high-throughput in silico reactive chemistry experimentation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Aqueous Phase Aldol Condensation of Formaldehyde and Acetone on Anatase TiO2(101) Surface: A Theoretical Investigation

A mechanistic understanding of catalytic reactions at solid-liquid interface is limited both experimentally and theoretically but attracts much interest. Using density functional theory calculations (DFT) and ab initio molecular dynamics (AIMD) simulations, we investigated the effect of liquid water on a-H abstraction, C-C coupling, and dehydration steps of aldol condensation of formaldehyde and acetone on an anatase TiO2(101) surface. The existence of the aqueous phase lowered the Gibbs energy of activation of dehydration step pronouncedly from 187 to 74 kJ/mol through proton transfer mechanism, making the hydrogenation pathway more favorable in the aqueous phase. In contrast, a mixed route prevails in the vapor phase. This work provides insights into the effect of the bulk water through a proton transfer mechanism on the dehydrogenation, C-C coupling, and dehydration

anatase TiO2, formaldehyde and acetone, aldol cond↗

Nature of molybdenum carbide surfaces for catalytic hydrogen dissociation using machine-learned potentials: an ensemble-averaged perspective

Molybdenum carbides with an electronic structure similar to noble metals have gained attention as a promising low-cost catalyst for biomass valorization and the hydrogen evolution reaction. However, our fundamental understanding of the catalyst surface and how different phases of these catalysts behave at varying reaction conditions is limited to ground state density functional theory calculations as ab initio molecular dynamics (AIMD) is computationally prohibitive at relevant length and time scales. Here, in this work, we train a multi-atomic cluster expansion (MACE) machine-learned interatomic potentials (MLIP) to study hydrogen dissociation and dynamics over Mo, δ-MoC, α-Mo 2 C, and β-Mo 2 C surfaces at varying temperatures and hydrogen partial pressures. Our simulations identify unique and different molecular and atomic hydrogen adsorption sites on different surfaces that do not depend on the temperature. At low hydrogen pressures, the surface coverage is monolayer, which transitions to two-layer adsorption at higher pressures. We find that atomic hydrogen diffusion and recombinations are preferred over molybdenum atom hollow sites, while the diffusion over carbon-terminated facets was negligible, signifying particularly strong C–H interactions. In contrast, molecular hydrogen adsorption occurs mostly atop Mo or the bridging sites. At a comparable hydrogen loading, β-Mo 2 C (001) is the most active surface for hydrogen dissociation reaction. This work provides insights into the dynamic nature of the hydrogen dissociation chemistry and the diversity of hydrogen adsorption sites on molybdenum carbides.

08 HYDROGEN↗

Dynamic Observation of Dendritic Quasicrystal Growth upon Laser-Induced Solid-State Transformation

In this paper, we report the laser-induced solid-state transformation between a periodic “approximant” and quasicrystal in the Al-Cr system during rapid quenching. Dynamic transmission electron microscopy allows us to capture in situ the dendritic growth of the metastable quasicrystals. The formation of dendrites during solid-state transformation is a rare phenomenon, which we attribute to the structural similarity between the two intermetallics. Through ab initio molecular dynamics simulations, we identify the dominant structural motif to be a 13-atom icosahedral cluster transcending the phases of matter.

36 MATERIALS SCIENCE↗

Elucidating the atomistic origin of anharmonicity in tetragonal CH3NH3PbI3 with Raman scattering

Halide perovskite (HP) semiconductors exhibit unique strong coupling between the electronic and structural dynamics. In this work, we use Raman polarization-orientation (PO) measurements and ab initio molecular dynamics (AIMD) to investigate the origin and temperature evolution of the strong structural anharmonicity throughout the tetragonal phase of CH 3 NH 3 PbI 3 . Raman PO measurements reveal a soft modelike spectral feature. This mode shows an unusual continuous increase in damping with temperature which is indicative of an anharmonic potential surface. The analysis of AIMD trajectories identifies two major sources of anharmonicity: the orientational unlocking of the [CH 3 NH 3 ] + ions and large-amplitude octahedral tilting that continuously increases with temperature. Our work suggests that the standard phonon picture cannot describe the structural dynamics of tetragonal CH 3 NH 3 PbI 3 .

36 MATERIALS SCIENCE↗

Gravitational Wave Signal for Quark Matter with Realistic Phase Transition

The cores of neutron stars (NSs) near the maximum mass can realize a transitional change to quark matter (QM). Gravitational waves from binary NS mergers are expected to convey information about the equation of state (EOS) sensitive to the QM transition. Here, in this work, we present the first results of gravitational wave simulation with the realistic EOS that is consistent with ab initio approaches of χ⁢EFT and pQCD and is assumed to go through smooth crossover. We compare them to results obtained with another EOS with a first-order hadron-quark phase transition. Our results suggest that early collapse to a black hole in the post-merger stage after NS merger robustly signifies softening of the EOS associated with the QM onset in the crossover scenario. The nature of the hadron-quark phase transition can be further constrained by the condition that electromagnetic counterparts should be energized by the material left outside the remnant black hole.

79 ASTRONOMY AND ASTROPHYSICS↗