Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Structural phase transition”

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 1,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Topological phase transition and quantum spin Hall state in TlBiS{sub 2}

We have investigated the bulk and surface electronic structures and band topology of TlBiS{sub 2} as a function of strain and electric field using ab-initio calculations. In its pristine form, TlBiS{sub 2} is a normal insulator, which does not support any non-trivial surface states. We show however that a compressive strain along the (111) direction induces a single band inversion with Z{sub 2} = (1;000), resulting in a Dirac cone surface state with a large in-plane spin polarization. Our analysis shows that a critical point lies between the normal and topological phases where the dispersion of the 3D bulk Dirac cone at the Γ-point becomes nearly linear. The band gap in thin films of TlBiS{sub 2} can be tuned through an out-of-the-plane electric field to realize a topological phase transition from a trivial insulator to a quantum spin Hall state. An effective k·p model Hamiltonian is presented to simulate our first-principles results on TlBiS{sub 2}.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2021

Semiconductor to Metal Transition in the Solids/Melts of Te and Pseudo-Binary of Hg1-xCdxTe for x = 0, 0.1 and 0.2

This paper reports the experimentally measured thermal and electrical conductivities from which the Lorenz numbers as functions of temperature were derived for the solids and melts of elemental Te and samples of Hg1-xCdxTe of x = 0, 0.1 and 0.2. The structural changes in the vicinity of the solid-liquid phase transition of elemental Te and various Hg1-xCdxTe solid solutions were assessed by closely examining these experimental properties and the derived Lorenz numbers. The values of Lorenz numbers for the solids of these samples were higher than the value of L0, i.e., the Lorenz number for the free electron gas derived by Sommerfeld, implying the characteristics of non-degenerate semiconductor. With the increasing value of x, or the larger bandgap, the materials become more non-degenerate with larger deviation from L0. As the solids started to melt, the measured values of Lorenz number for all samples started to decrease and eventually reached and maintained at values close to L0 at higher temperatures. The trend of the Lorenz number indicates that the transition from non-degenerate semiconductor of the solid samples to metal or degenerate semiconductor when the samples temperature reached above their liquidus temperatures.

thermal conductivity↗
DOE OSTI2016

Crystalline structure and magnetic behavior of the Ni 41 Mn 39 In 12 Co 8 alloy demonstrating giant magnetocaloric effect

Magnetic cooling is a green cooling technology, which is more energy efficient than existing fluid-compression cooling machines. Ni 41 Mn 39 In 12 Co 8 alloy, which demonstrates promising magnetocaloric performances, was investigated using neutron diffraction and thermomagnetic measurements. The austenite structure is cubic L2 1 $(Fm\bar{3}m),$ while that of the martensite is a mix of 8 and 6 M modulated monoclinic structures $(P\quad 1\quad 2/m\quad 1).$ The austenitic site occupancy refinements reveal that all substituting Co atoms occupy Ni-sites. Most Mn atoms (65%) are in the Mn-sites and the rest go to In-sites (about 35%) and Ni-sites (less than 5%). This disorder of the magnetic atoms (Mn, Ni and Co) in the austenitic phase remains unchanged during the martensitic transition. The distortions of the interatomic distances due to the modulation of the martensitic structures further enhance the disorder in the magnetic interactions. Thermomagnetic measurements indicate that the austenitic phase is ferromagnetic. Cooling to below 250 K, where the alloy loses its ferromagnetic nature, and down to 50 K, the lack of any antiferromagnetic Bragg peaks suggests no antiferromagnetic ordering in the martensitic phase. Here, at very low temperatures in the martensitic phase, spin glass magnetic nature is identified by magnetic measurements, and the spin-glass transition temperature is ~19 K.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Time-resolved x-ray diffraction and Raman studies of the phase transition mechanisms of methane hydrate

The mechanisms by which methane hydrate transforms from an sI to sH structure and from an sH to filled-ice Ih structure were examined using time-resolved X-ray diffractometry (XRD) and Raman spectroscopy in conjunction with charge-coupled device camera observation under fixed pressure conditions. The XRD data obtained for the sI–sH transition at 0.8 GPa revealed an inverse correlation between sI and sH, suggesting that the sI structure is replaced by sH. Meanwhile, the Raman analysis demonstrated that although the 12-hedra of sI are retained, the 14-hedra are replaced sequentially by additional 12-hedra, modified 12-hedra, and 20-hedra cages of sH. With the sH to filled-ice Ih transition at 1.8 GPa, both the XRD and Raman data showed that this occurs through a sudden collapse of the sH structure and subsequent release of solid and fluid methane that is gradually incorporated into the filled-ice Ih to complete its structure. This therefore represents a typical reconstructive transition mechanism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Grazing-Angle Neutron Diffraction Study of the Water Distribution in Membrane Hemifusion: From the Lamellar to Rhombohedral Phase

The water distribution between lipid bilayers is important in understanding the role of the hydration force at different steps of the membrane fusion pathway. Here in this study, we used grazing-angle neutron diffraction to map out the water distribution in lipid bilayers transiting from a lamellar structure to the hemifusion “stalk” structure in a rhombohedral phase. Under osmotic pressure exerted by different levels of relative humidity, the lipid membrane sample was maintained in equilibrium at different lattices suitable for neutron diffraction. The D 2 O used to hydrate the lipid membrane sample stood out from the lipid in the reconstructed structure because of its much higher coherent neutron scattering length density. The density map indicates that water dissociated from the headgroup in the lamellar phase. In the rhombohedral phase, water was significantly reduced and was squeezed into pockets around the stalk. This study complements earlier structural studies by grazing-angle X-ray diffraction, which is sensitive to only the parts of the structure with high electron density (such as phosphors). The experiment also demonstrated that the recently developed time-of-flight small-angle neutron scattering beamline at the Spallation Neutron Source is suitable for grazing-angle neutron diffraction to provide the structures of large unit cells on the order of a few nanometers, such as biomembrane structures.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1997

Peptide Folding and Translocation Across the Water-Membrane Interface

The ability of small peptides to organize at aqueous interfaces was examined by performing a series of large-scale, molecular dynamics computer simulations of several peptides composed of two amino acids, nonpolar leucine (L) and polar glutamine (Q). The peptides differed in size and sequence of the amino acids. Studies on dipeptides LL, LQ, QL and QQ were extended to two heptamers, LQQLLQL and LQLQLQL, designed to maximize interfacial stability of an alpha-helix and a beta-strand, respectively, by exposing polar side chains to water and nonpolar side chains to a nonpolar phase. Finally, a transition of an undecamer, composed entirely of leucine residues, from a disordered structure in water to an alpha-helix in a nonpolar phase representing the interior of the membrane was investigated. Complete folding of a peptide in solution was accomplished for the first time in computer simulations. The simulations revealed several basic principles governing the sequence-dependent organization of peptides at interfaces. Short peptides tend to accumulate at interfaces and acquire ordered structures, providing that they have a proper sequence of polar and nonpolar amino acids. The dominant factor determining the interfacial structure of peptides is the hydrophobic effect, which is manifested at aqueous interfaces as a tendency for polar and nonpolar groups of the solute to segregate into the aqueous and nonpolar phases, respectively. If peptides consist of nonpolar residue's only, they become inserted into the nonpolar phase. As demonstrated by the example of the leucine undecamer, such peptides fold into an alpha-helix as they partition into the nonpolar medium. The folding proceeds through an intermediate, called 3-10-helix, which remains in equilibrium with the alpha-helix. Once in the nonpolar environment, the peptides can readily change their orientation with respect to the interface from parallel to perpendicular, especially in response to local electric fields. The ability of nonpolar peptides to modify both the structure and orientation with respect to the interface from parallel to perpendicular, especially in response to local electric fields. The ability of nonpolar peptides to modify both the structure and orientation with changing external conditions may have provided a simple mechanism of transmitting signals from the environment to the interior of a cell.

Pohorille, Andrew↗
DOE OSTI2023

Metal–insulator transition in composition-tuned nickel oxide films

Abstract Thin films of the solid solution Nd 1 − x La x NiO 3 are grown in order to study the expected 0 K phase transitions at a specific composition. We experimentally map out the structural, electronic and magnetic properties as a function of x and a discontinuous, possibly first order, insulator–metal transition is observed at low temperature when x = 0.2. Raman spectroscopy and scanning transmission electron microscopy show that this is not associated with a correspondingly discontinuous global structural change. On the other hand, results from density functional theory (DFT) and combined DFT and dynamical mean field theory calculations produce a 0 K first order transition at around this composition. We further estimate the temperature-dependence of the transition from thermodynamic considerations and find that a discontinuous insulator–metal transition can be reproduced theoretically and implies a narrow insulator–metal phase coexistence with x . Finally, muon spin rotation ( µ SR) measurements suggest that there are non-static magnetic moments in the system that may be understood in the context of the first order nature of the 0 K transition and its associated phase coexistence regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Room-temperature photo-induced martensitic transformation in a protein crystal

Martensitic transformations are the first-order crystal-to-crystal phase transitions that occur mostly in materials such as steel, alloys and ceramics, thus having many technological applications. These phase transitions are rarely observed in molecular crystals and have not been detected in protein crystals. Reversibly switchable fluorescent proteins are widely used in biotechnology, including super-resolution molecular imaging, and hold promise as candidate biomaterials for future high-tech applications. Here, we report on a reversibly switchable fluorescent protein, Tetdron, whose crystals undergo a photo-induced martensitic transformation at room temperature. Room-temperature X-ray crystallography demonstrates that at equilibrium Tetdron chromophores are all in the trans configuration, with an ∼1:1 mixture of their protonated and deprotonated forms. Irradiation of a Tetdron crystal with 400 nm light induces a martensitic transformation, which results in Tetdron tetramerization at room temperature revealed by X-ray photocrystallography. Crystal and solution spectroscopic measurements provide evidence that the photo-induced martensitic phase transition is coupled with the chromophore deprotonation, but no trans–cis isomerization is detected in the structure of an irradiated crystal. It is hypothesized that protein dynamics assists in the light-induced proton transfer from the chromophore to the bulk solvent and in the ensuing martensitic phase transition. The unique properties of Tetdron may be useful in developing novel biomaterials for optogenetics, data storage and nanotechnology.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ultrafast Martensitic Phase Transition Driven by Intense Terahertz Pulses

We report on an ultrafast nonequilibrium phase transition with a strikingly long-lived martensitic anomaly driven by above-threshold single-cycle terahertz pulses with a peak field of more than 1 MV/cm. A nonthermal, terahertz-induced depletion of low-frequency conductivity in Nb 3 Sn indicates increased gap splitting of high-energy 12 bands by removal of their degeneracies, which induces the martensitic phase above their equilibrium transition temperature. In contrast, optical pumping leads to a Γ 12 gap thermal melting. Such light-induced nonequilibrium martensitic phase exhibits a substantially enhanced critical temperature up to ~100 K, i.e., more than twice the equilibrium temperature, and can be stabilized beyond technologically relevant, nanosecond time scales. Together with first-principle simulations, we identify a compelling terahertz tuning mechanism of structural order via Γ 12 phonons to achieve the ultrafast phase transition to a metastable electronic state out of equilibrium at high temperatures far exceeding those for equilibrium states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Quantum Computation of Dynamical Quantum Phase Transitions and Entanglement Tomography in a Lattice Gauge Theory

Strongly coupled gauge theories far from equilibrium may exhibit unique features that could illuminate the physics of the early universe and of hadron and ion colliders. Studying real-time phenomena has proven challenging with classical-simulation methods but is a natural application of quantum simulation. To demonstrate this prospect, we quantum compute nonequal-time correlation functions and perform entanglement tomography of nonequilibrium states of a simple lattice gauge theory, the Schwinger model, using a trapped-ion quantum computer by IonQ Inc. As an ideal target for near-term devices, a recently predicted dynamical quantum phase transition in this model is studied by preparing, quenching, and tracking the subsequent nonequilibrium dynamics in three ways: (i) overlap echos signaling dynamical transitions, (ii) nonequal-time correlation functions with an underlying topological nature, and (iii) the entanglement structure of nonequilibrium states, including entanglement Hamiltonians. These results constitute the first observation of a dynamical quantum phase transition in a lattice gauge theory on a quantum computer and are a first step toward investigating topological phenomena in nuclear and high-energy physics using quantum technologies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Pair mobility functions for rigid spheres in concentrated colloidal dispersions: Force, torque, translation, and rotation

The formulation of detailed models for the dynamics of condensed soft matter including colloidal suspensions and other complex fluids requires accurate description of the physical forces between microstructural constituents. In dilute suspensions, pair-level interactions are sufficient to capture hydrodynamic, interparticle, and thermodynamic forces. In dense suspensions, many-body interactions must be considered. Prior analytical approaches to capturing such interactions such as mean-field approaches replace detailed interactions with averaged approximations. However, long-range coupling and effects of concentration on local structure, which may play an important role in, e.g., phase transitions, are smeared out in such approaches. An alternative to such approximations is the detailed modeling of hydrodynamic interactions utilizing precise couplings between moments of the hydrodynamic traction on a suspended particle and the motion of that or other suspended particles. For two isolated spheres, a set of these functions was calculated by Jeffrey and Onishi [J. Fluid Mech. 139, 261–290 (1984)] and Jeffrey [J. Phys. Fluids 4, 16–29 (1992)]. Along with pioneering work by Batchelor, these are the touchstone for low-Reynolds-number hydrodynamic interactions and have been applied directly in the solution of many important problems related to the dynamics of dilute colloidal dispersions [G. K. Batchelor and J. T. Green, J. Fluid Mech. 56, 375–400 (1972) and G. K. Batchelor, J. Fluid Mech. 74, 1–29 (1976)]. Toward extension of these functions to concentrated systems, here we present a new stochastic sampling technique to rapidly calculate an analogous set of mobility functions describing the hydrodynamic interactions between two hard spheres immersed in a suspension of arbitrary concentration, utilizing accelerated Stokesian dynamics simulations. These mobility functions provide precise, radially dependent couplings of hydrodynamic force and torque to particle translation and rotation, for arbitrary colloid volume fraction ϕ. The pair mobilities (describing entrainment of one particle by the disturbance flow created by another) decay slowly with separation distance: as 1/r, for volume fractions 0.05 ≤ ϕ ≤ 0.5. For the relative mobility, we find an initially rapid growth as a pair separates, followed by a slow, 1/r growth. Up to ϕ ≤ 0.4, the relative mobility does not reached the far-field value even beyond separations of many particle sizes. In the case of ϕ = 0.5, the far-field asymptote is reached but only at a separation of eight radii and after a slow 1/r growth. At these higher concentrations, the coefficients also reveal liquid-like structural effects on pair mobility at close separations. These results confirm that long-range many-body hydrodynamic interactions are an essential part of the dynamics of concentrated systems and that care must be taken when applying renormalization schemes.

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

The Ho–Ni–Ge system: Isothermal section and new rare-earth nickel germanides

The Ho–Ni–Ge system has been investigated at 1070 K and up to ~60 at% Ho by X-ray diffraction and microprobe analyses. Besides the eight known compounds, HoNi{sub 5}Ge{sub 3} (YNi{sub 5}Si{sub 3}-type), HoNi{sub 2}Ge{sub 2} (CeAl{sub 2}Ga{sub 2}-type), Ho{sub 2}NiGe{sub 6} (Ce{sub 2}CuGe{sub 6}-type), HoNiGe{sub 3} (SmNiGe{sub 3}-type), HoNi{sub 0.2÷0.6}Ge{sub 2} (CeNiSi{sub 2}-type), Ho{sub 37÷34}Ni{sub 6÷24}Ge{sub 57÷42} (AlB{sub 2}-type), HoNiGe (TiNiSi-type), Ho{sub 3}NiGe{sub 2} (La{sub 3}NiGe{sub 2}-type), the ternary system contains four new compounds: Ho{sub 3}Ni{sub 11}Ge{sub 4} (Sc{sub 3}Ni{sub 11}Ge{sub 4}-type), HoNi{sub 3}Ge{sub 2} (ErNi{sub 3}Ge{sub 2}-type), Ho{sub 3}Ni{sub 2}Ge{sub 3} (Hf{sub 3}Ni{sub 2}Si{sub 3}-type) and ~Ho{sub 5}Ni{sub 2}Ge{sub 3} (unknown structure). Quasi-binary solid solutions were observed at 1070 K for Ho{sub 2}Ni{sub 17}, HoNi{sub 5}, HoNi{sub 7}, HoNi{sub 3}, HoNi{sub 2}, HoNi and Ho{sub 2}Ge{sub 3}, but no detectable solubility was found for the other binary compounds in the Ho–Ni–Ge system. Based on the magnetization measurements, the HoNi{sub 5}Ge{sub 3}, HoNi{sub 3}Ge{sub 2} and Ho{sub 3}Ni{sub 11}Ge{sub 4} (and isostructural (Tb, Dy){sub 3}Ni{sub 11}Ge{sub 4}) compounds have been found to show paramagnetic behavior down to 5 K, whereas Ho{sub 3}Ni{sub 2}Ge{sub 3} exhibits an antiferromagnetic transition at ~7 K. Additionally, the crystal structure of the new isostructural phases (Y, Yb)Ni{sub 3}Ge{sub 2} (ErNi{sub 3}Ge{sub 2}-type), Er{sub 3}Ni{sub 11}Ge{sub 4} (Sc{sub 3}Ni{sub 11}Ge{sub 4}-type) and (Y, Tb, Dy, Er, Tm){sub 3}Ni{sub 2}Ge{sub 3} (Hf{sub 3}Ni{sub 2}Si{sub 3}-type) has been also investigated. - Graphical abstract: The Ho–Ni–Ge system has been investigated at 1070 K and up to ~60 at.% Ho by X-ray and microprobe analyses. Besides the eight known compounds, i.e. HoNi{sub 5}Ge{sub 3} (YNi{sub 5}Si{sub 3}-type), HoNi{sub 2}Ge{sub 2} (CeAl{sub 2}Ga{sub 2}-type), Ho{sub 2}NiGe{sub 6} (Ce{sub 2}CuGe{sub 6}-type), HoNiGe{sub 3} (SmNiGe{sub 3}-type), HoNi{sub 0.2÷0.6}Ge{sub 2} (CeNiSi{sub 2}-type), Ho{sub 37÷34}Ni{sub 6÷24}Ge{sub 57÷42} (AlB{sub 2}-type), HoNiGe (TiNiSi-type), Ho{sub 3}NiGe{sub 2} (La{sub 3}NiGe{sub 2}-type), the ternary system contains four new compounds: Ho{sub 3}Ni{sub 11}Ge{sub 4} (Sc{sub 3}Ni{sub 11}Ge{sub 4}-type), HoNi{sub 3}Ge{sub 2} (ErNi{sub 3}Ge{sub 2}-type), Ho{sub 3}Ni{sub 2}Ge{sub 3} (Hf{sub 3}Ni{sub 2}Si{sub 3}-type) and ~Ho{sub 5}Ni{sub 2}Ge{sub 3} (unknown structure). Quasi-binary solid solutions were found to form at 1070 K from the binary Ho{sub 2}Ni{sub 17}, HoNi{sub 5}, HoNi{sub 7}, HoNi{sub 3}, HoNi{sub 2}, HoNi and Ho{sub 2}Ge{sub 3} compounds, while no detectable solubility was observed for the other binary compounds in the Ho–Ni–Ge system. Based on the magnetization measurements, the HoNi{sub 5}Ge{sub 3}, HoNi{sub 3}Ge{sub 2} and Ho{sub 3}Ni{sub 11}Ge{sub 4} (and isostructural (Tb, Dy){sub 3}Ni{sub 11}Ge{sub 4}) compounds have been found to show paramagnetic behavior down to 5 K, whereas Ho{sub 3}Ni{sub 2}Ge{sub 3} exhibits an antiferromagnetic transition at ~7 K. Additionally, the crystal structure of the new isostructural phases (Y, Yb)Ni{sub 3}Ge{sub 2} (ErNi{sub 3}Ge{sub 2}-type), Er{sub 3}Ni{sub 11}Ge{sub 4} (Sc{sub 3}Ni{sub 11}Ge{sub 4}-type) and (Y, Tb, Dy, Er, Tm){sub 3}Ni{sub 2}Ge{sub 3} (Hf{sub 3}Ni{sub 2}Si{sub 3}-type) has been also investigated. - Highlights: • Ho–Ni–Ge system has been investigated at 1070 K and up to ~60 at% Ho. • Eight known ternary holmium nickel germanides were confirmed in Ho–Ni–Ge. • Four new holmium nickel germanides were detected in Ho–Ni–Ge. • Eight new rare earth nickel germanides were detected in (Y, Tb, Dy, Er–Yb)–Ni–Ge. • HoNi{sub 5}Ge{sub 3}, HoNi{sub 3}Ge{sub 2},(Tb, Dy, Ho){sub 3}Ni{sub 11}Ge{sub 4} and Ho{sub 3}Ni{sub 2}Ge{sub 3} are paramagnet down to 5–7 K.

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

A new structural relaxation pathway of low-density amorphous ice

Low-density amorphous (LDA) ice is involved in critical cosmological processes and has gained prominence as one of the at least two distinct amorphous forms of ice. Despite these accolades, we still have an incomplete understanding of the structural diversity that is encompassed within the LDA state and the dynamic processes that take place upon heating LDA. Heating the high-pressure ice VIII phase at ambient pressure is a remarkable example of temperature-induced amorphisation yielding LDA. We investigate this process in detail using X-ray diffraction and Raman spectroscopy and show that the LDA obtained from ice VIII is structurally different from the more “traditional” states of LDA which are approached upon thermal annealing. This new structural relaxation pathway involves an increase of structural order on the intermediate range length scale. In contrast with other LDA materials the local structure is more ordered initially and becomes slightly more disordered upon annealing. We also show that the cascade of phase transitions upon heating ice VIII at ambient pressure includes the formation of ice IX which may be connected with the structural peculiarities of LDA from ice VIII. Overall, this study shows that LDA is a structurally more diverse material than previously appreciated.

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

Synthesis and physical properties of the 10.6 K ferromagnet NdIr3

The magnetic, transport, and thermodynamic properties of $PuNi_3$ are reported here. Its $PuN_3$-type crystal structure (space group $R-3m$), with lattice parameters $ a = 5.3262(1) Å$ and $c = 26.2218(3) Å_1$ was confirmed by powder x-ray diffraction. Our measurements indicate that $NdIr_3$ reflects a previously unreported paramagnetic to ferromagnetic phase transition below $T_C = 10.6K$ (determined by using the Arrott plot). The magnetic susceptibility obeys the Curie-Weiss law with an effective magnetic moment of μ eff = 3.63 ( 2 ) μ B / Nd and a paramagnetic Curie temperature θ P = 5.4 ( 6 ) K . The heat-capacity anomaly at T C confirms a bulk nature of the transition, though Δ C p = 11.7 J mo l - 1 K - 1 is lower than expected for $J = 9/2$ and instead close to the $J = 1/2$ system. This indicates that Nd ions are subject to the crystalline electrical field that removes spin degeneracy and leaves the Nd ions in a doublet ground state. Resistivity and heat-capacity measurements reveal an anomaly at around 70 K. The origin of this anomaly is unknown, but the lack of thermal hysteresis suggests that it is a second-order phase transition and may be related to electronic instabilities. The calculated electronic structure is observed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Predictions of Boron Phase Stability Using an Efficient Bayesian Machine Learning Interatomic Potential

Thermodynamic phase stability of three elemental boron allotropes, i.e., α-B, β-B, and γ-B, was investigated using a Bayesian interatomic potential trained via a sparse Gaussian process (SGP). SGP potentials trained with datasets from on-the-fly active learning achieve quantum mechanical level accuracy when employed in molecular dynamics simulations to predict wide-ranging thermodynamic, structural, and vibrational properties. The simulated phase diagram (500~1400 K and 0~16 GPa) agrees with experimental measurements. The SGP-based MD simulations also successfully predicted that the B13 defect is critical in stabilizing β-B below 700 K. At higher temperatures, the entropy becomes the dominant factor, making β-B the more stable phase over α-B. Furthermore, this Letter demonstrates that SGP potentials based on a training set consisting of defect-free-only systems could make correct predictions of defect-related phenomena in solid-state crystals, paving the path to investigate crystal phase stability and transitions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗