Search NASA⌕ Search

SEARCH · Search NASA

Results for “Phase 1”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Experimental X-ray Charge-Density Studies–A Suitable Probe for Superconductivity? A Case Study on MgB 2

Case studies of 1T-TiSe 2 and YBa 2 Cu 3 O 7-δ have demonstrated that X-ray diffraction (XRD) studies can be used to trace even subtle structural phase transitions which are inherently connected with the onset of superconductivity in these benchmark systems. However, the utility of XRD in the investigation of superconductors like MgB 2 lacking an additional symmetry-breaking structural phase transition is not immediately evident. Nevertheless, high-resolution powder XRD experiments on MgB 2 in combination with maximum entropy method analyses hinted at differences between the electron density distributions at room temperature and 15 K, that is, below the T c of approx. 39 K. The high-resolution single-crystal XRD experiments in combination with multipolar refinements presented here can reproduce these results but show that the observed temperature-dependent density changes are almost entirely due to a decrease of atomic displacement parameters as a natural consequence of a reduced thermal vibration amplitude with decreasing temperature. Our investigations also shed new light on the presence or absence of magnesium vacancies in MgB 2 samples–a defect type claimed to control the superconducting properties of the compound. Here, we propose that previous reports on the tendency of MgB 2 to form non-stoichiometric Mg 1–x B 2 phases (1 – x ~ 0.95) during high-temperature (HT) synthesis might result from the interpretation of XRD data of insufficient resolution and/or usage of inflexible refinement models. Indeed, advanced refinements based on an Extended Hansen–Coppens multipolar model and high-resolution X-ray data, which consider explicitly the contraction of core and valence shells of the magnesium cations, do not provide any significant evidence for the formation of non-stoichiometric Mg 1–x B 2 phases during HT synthesis.

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

Pressure dependence of antiferromagnetic and superconducting phases in U 2 Rh 1–x Pt x C 2

In this work, we report temperature (T)- and pressure (P)-dependent resistivity measurements on the isostructural compounds U 2 Rh 1–x Pt x C 2 (x = 0, 0.5, and 0.9) from which we construct a T–P–x phase diagram. Compounds with x = 0 and x = 0.5 are antiferromagnets below 22.1 and 9.4 K, respectively, and their Néel temperature (T N ) decreases under applied pressure. For x = 0, the critical pressure P c required to suppress T N to zero temperature is projected to be about 8.8 GPa, but P c for x = 0.5 is between 1.6 and 2.1 GPa. At atmospheric pressure, increasing Pt concentration in U 2 Rh (1–x) Pt x C 2 tunes magnetic transition temperatures to zero at a critical value of x c ≈ 0.7, and, consequently, we surmise the existence of a quantum-phase boundary in the P–x plane at T = 0 K that extends from (P = 0, x = xc) to (P c ≈ 8.8 GPa, x = 0). For x = 0.9, superconductivity appears at T c = 1.09 K , which decreases at a rate of ≈ –1 K/GPa that is nearly twice that found for U 2 PtC 2 whose T c is 1.47 K. Together, these results indicate that domes of magnetism and superconductivity formed with T–P–x variations are detached and that the two broken symmetries are independent of each other. Fluctuations in average composition produce rare regions that play an important role in determining physical properties of materials with noninteger x.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Tracking the phase changes in micelle-based NiGa nanocatalysts for methanol synthesis under activation and working conditions

The hydrogenation of CO 2 into high energy density fuels such as methanol, where the required H 2 is obtained from renewable sources, is of utmost importance for a sustainable society. In recent years, NiGa alloys have attracted attention as promising catalyst material systems for the hydrogenation of CO 2 into methanol at ambient pressures. They thus represent an energy-saving alternative to the Cu-based catalysts employed in today's catalytic industry that require high pressures for the CO 2 hydrogenation. However, the underlying reaction mechanisms for the NiGa system are still under debate. One of the challenges here is to unravel the evolution and coexistence of the different species in the heterogeneous NiGa catalyst system under activation and reaction conditions. To shed light on their evolution under activation in H 2 and their catalytic roles under CO 2 hydrogenation working conditions on well-defined Ni 3 Ga 1 and Ni 5 Ga 3 nanoparticle (NP) catalysts, we employed a multi-probe approach in this study. It included advanced machine learning-based analysis of operando X-ray absorption spectroscopy data combined with operando powder X-ray diffraction and near ambient pressure X-ray photoelectron spectroscopy measurements, as well as reactivity studies using bed-packed mass flow reactors. In addition, we employed atomic force microscopy and scanning transmission electron microscopy for structural characterization. Under H 2 activation at 1 bar total pressure, we concluded the formation of metallic Ni, starting for Ni 3 Ga 1 at 300 °C, and for Ni 5 Ga 3 at 400 °C. At higher temperatures, the formation of NiGa alloys follows. The α'-Ni 3 Ga 1 alloy phase is predominantly formed for the Ni 3 Ga 1 NPs, while the coexistence of α'-Ni 3 Ga 1 , δ-Ni 5 Ga 3 and Ga 2 O 3 phases is observed for the Ni 5 Ga 3 NPs after the H 2 activation. The formation of the Ga 2 O 3 phase also results in the presence of excess metallic Ni. Under CO 2 hydrogenation reaction conditions, Ga partially oxidizes again to form a Ga 2 O 3 -rich particle shell for both NP compositions, yet, to a larger extent for the Ni 3 Ga 1 NPs, which, in turn, feature a higher amount of excess Ni. We reveal that metallic Ni is responsible for the high selectivity of the Ni 3 Ga 1 NPs towards the production of methane in our catalytic tests. In contrast, the Ni 5 Ga 3 NPs display a strong selectivity toward methanol production (>92%), more than one order of magnitude higher than that for the Ni 3 Ga 1 NPs, which we ascribe to the presence of the δ-Ni 5 Ga 3 phase.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Tetragonal structure and uniaxial magnetic anisotropy in the arc-melted (Ce,Zr) 2 (Fe,M) 17 (M = Mo, W, Co) alloys

The tetragonal ThMn 12 -type structure is stabilized in (Ce 1-x Zr x ) 2 T 16 M (x = 0.2 – 0.3; T = Fe or Fe/Co, M = Mo or W) arc-melted alloys. Approximately 5 at% of Mo or W admixture is sufficient to transform the hexagonal Th 2 Ni 17 -type structure of (Ce 1-x Zr x ) 2 Fe 17 into the tetragonal ThMn 12 -type structure of nearly single-phase (Ce 1-x Zr x ) 2 Fe 16 M and/or (Ce 1-x Zr x ) 2 Fe 15 CoM bulk alloys. X-ray Rietveld refinements reveal that Zr and Mo (W) substitute different sites in the tetragonal crystal structure. Zirconium preferentially replaces Ce, whereas Mo (W) substitutes Fe. At room temperature, the tetragonal phases exhibit strong ferromagnetism and a uniaxial magneto-crystalline anisotropy with anisotropy fields of 15 – 17 kOe. These materials possess room-temperature saturation magnetizations of 95 – 105 emu/g and Curie temperature of 420 – 505 K. At some concentrations, concurrent Zr and Mo (W) site occupancy facilitates 1:12 structure formation in bulk alloys with minimal presence of non-magnetic element in the Fe sublattice, thus securing the highest known magnetic moment per Fe atom (∼ 1.37 µ B ) in this type of materials. The intrinsic magnetic characteristics, as well as the absence of critical rare earths (RE), make these compounds interesting for development as low-cost permanent magnets. In conclusion, small Co additions improve the Curie temperature, especially if combined with W.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Pressure-Induced Chemical Bonding Effects on Lattice and Magnetic Instabilities in Antiferromagnetic Insulating CaMn 2 Sb 2

Exotic quantum phenomena often emerge near an electronic delocalization transition (EDT) from an antiferromagnetic insulating phase to a strongly correlated metallic state under pressure. We report the pressure-induced structural and magnetic evolution of the antiferromagnetic insulator CaMn 2 Sb 2 . Single-crystal X-ray diffraction reveals a first-order phase transition near 5.4 GPa from a trigonal P-3m1 structure to a monoclinic P2 1 /m phase accompanied by a ∼7% volume collapse. Residual electron density analysis at intermediate pressures reveals charge localization along Mn–Sb chains, signaling electronic instability preceding the structural transition. Bonding analysis indicates anisotropic Mn–Sb orbital reconfiguration under pressure, driving a distorted square-pyramidal geometry. Neutron scattering confirms the transition and identifies a pressure-induced incommensurate magnetic order, distinct from the ambient antiferromagnetic state. In the monoclinic phase, zigzag Mn chains exhibit antiferromagnetic coupling along the ac-plane, enabled by enhanced orbital overlap. These results establish CaMn 2 Sb 2 as a model system for studying the coupling of structural distortion, charge redistribution, and magnetic order in layered Mn pnictides under pressure.

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

Novel Hydrogen Clathrate Hydrate

We report a new hydrogen clathrate hydrate synthesized at 1.2 GPa and 298 K documented by single-crystal x-ray diffraction, Raman spectroscopy, and first-principles calculations. The oxygen sublattice of the new clathrate hydrate matches that of ice II, while hydrogen molecules are in the ring cavities, which results in the trigonal R3c or R¯3c space group (proton ordered or disordered, respectively) and the composition of ( H 2 O) 6 H 2 . Raman spectroscopy and theoretical calculations reveal a hydrogen disordered nature of the new phase $C$$^{'}_{1}$, distinct from the well-known ordered C 1 clathrate, to which this new structure transforms upon compression and/or cooling. Furthermore, this new clathrate phase can be viewed as a realization of a disordered ice II, unobserved before, in contrast to all other ordered ice structures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Competing Charge/Spin-Stripe and Correlated Metal Phases in Trilayer Nickelates (Pr 1– x La x ) 4 Ni 3 O 8

We report low valent nickelates R n+1 Ni n O 2n+2 (R = rare earth) containing Ni 1+ (d 9 ) with a quasi-two-dimensional (quasi-2D) square-planar coordination geometry possess structural and electronic properties that are similar to those of high T c cuprates, including superconductivity itself in the doped infinite-layer (n = ∞) RNiO 2 system. Within this R n+1 Ni n O 2n+2 nickelate family, the crystallographic isomorphs Pr 4 Ni 3 O 8 and La 4 Ni 3 O 8 exhibit singularly different ground states: Pr 4 Ni 3 O 8 is metallic, and La 4 Ni 3 O 8 is a charge- and spin-stripe-ordered insulator. To explore and understand the ground state evolution from metallic Pr 4 Ni 3 O 8 to stripe-ordered La 4 Ni 3 O 8 in the R 4 Ni 3 O 8 family, we have grown a series of isovalent, substituted single crystals (Pr 1-x La x ) 4 Ni 3 O 8 . Combining thermodynamic, transport, magnetic, and synchrotron X-ray single crystal diffraction measurements, we reveal a T = 0 transition between metallic and stripe-insulator phase regions, with this putative quantum phase transition at x ≈ 0.45. We propose two possible models for (Pr 1-x La x ) 4 Ni 3 O 8 : an electronically inhomogeneous system that could serve as a candidate for exploring quantum Griffiths phase physics or a homogeneous system with a clean quantum critical point at the phase boundary

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

Superconductivity and phase diagrams of CaK(Fe 1- x Mn x ) 4 As 4 single crystals

In order to study the effects of Mn substitution on the superconducting and magnetic ground state of CaKFe 4 As 4 (T c = 35 K), members of the CaK (Fe 1–x Mn x ) 4 As 4 series have been synthesized by high-temperature solution growth in single-crystalline form and characterized by elemental analysis, thermodynamic, and transport measurements. These measurements show that the superconducting transition temperature decreases monotonically and is finally suppressed below 1.8 K as x is increased from 0 to 0.036. For x values greater than 0.016, signatures of a magnetic transition can be detected in both thermodynamic and transport measurements in which kinklike features allow for the determination of the transition temperature T* that increases as Mn substitution increases. Here, a temperature-composition (T–x) phase diagram is constructed, revealing a half-dome of superconductivity with the magnetic transition temperature T* appearing near 26 K for x ~ 0.017 and rising slowly up to 33 K for x ~ 0.036. In addition to the creation of the T–x phase diagram for CaK(Fe 1–x Mn x ) 4 As 4 , specific-heat data are used to track the jump in specific heat at T c ; the CaK(Fe 1–x Mn x ) 4 As 4 data do not follow the scaling of ΔC p with $T$$^{3}_{c}$ as many of the other Fe-based superconducting systems do. These data suggest that, as magnetic pair breaking is present, the jump in C p for a given T c is reduced. Elastoresistivity coefficients 2m 66 and m 11 – m 12 as a function of temperature are also measured. 2m 66 and m 11 – m 12 are qualitatively similar to CaK(Fe 1–x Ni x ) 4 As 4 . This may indicate that the magnetic order in Mn-substituted system may be still the same as CaK(Fe 1–x Mn x ) 4 As 4 . Superconductivity of CaK(Fe 1–x Mn x ) 4 As 4 is also studied as a function of magnetic field. A clear change in H' c2 (T)/T c , where H' c2 (T) is dH' c2 (T)/dT, at x ~ 0.015 is observed and probably is related to change of the Fermi surface due to magnetic order. Coherence lengths and the London penetration depths are also calculated based on H c1 and H c2 data. Coherence lengths as the function of x also show the changes near x = 0.015, again consistent with Fermi-surface changes associated with the magnetic ordering seen for higher-x values.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Pressure-temperature phase diagram of CaK ( Fe 1 – x Mn x ) 4 As 4 ( x = 0.024 )

Resistance measurements on single crystals of CaK(Fe 1–x Mn x ) 4 As 4 (x = 0.024) were performed under hydrostatic pressure up to 5.15 GPa. The pressure dependence of the magnetic and superconducting transition temperatures and that of the superconducting upper critical field are reported. Furthermore, our results show that upon increasing pressure, the magnetic transition temperature T N is suppressed, whereas the superconducting transition temperature T c first increases and then decreases, exhibiting a maximum at a pressure p c corresponding to the intersection of the T N (p) and T c (p) lines. In addition, a minimum in the normalized slope of the superconducting upper critical field as well as a change in the pressure dependence of the inferred superconducting coherence length are observed at p c , suggesting a difference in the Fermi surface of the paramagnetic and antiferromagnetic states. Finally, CaK(Fe 1–x Mn x ) 4 As 4 (x = 0.024) likely goes through a half-collapsed tetragonal phase transition at ~4.3 GPa, further demonstrating that the half-collapsed tetragonal transition pressure in the CaKFe 4 As 4 system is relatively insensitive to transition metal substitution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Distributed phase sensing using two-mode squeezed states in a truncated SU(1,1) interferometer

We demonstrate distributed phase sensing with a truncated SU(1,1) interferometer through the measurement of a linear combination of two phases distributed between the two beams of a two-mode squeezed state. We theoretically analyze the sensitivity enhancement with respect to the corresponding classical strategy and experimentally demonstrate a 2 dB quantum noise reduction in the measurement of a linear combination of two phases.

Hong, Seongjin↗
DOE OSTI2021

Simulations of ENSO Phase-locking in CMIP5 and CMIP6

The characteristics of El-Niño-Southern Oscillation (ENSO) phase-locking in observations and CMIP5 and CMIP6 models are examined in this study. Two metrics based on the peaking month histogram for all El Niño and La Niña events are adopted to delineate the basic features of ENSO phase-locking in terms of the preferred calendar month and strength of this preference. It turns out that most models are poor at simulating the ENSO phase-locking, either showing little peak strengths or peaking at the wrong seasons. By deriving ENSO’s linear dynamics based on the conceptual recharge oscillator (RO) framework through the seasonal linear inverse model (sLIM) approach, various simulated phase-locking behaviors of CMIP models are systematically investigated in comparison with observations. In observations, phase-locking is mainly attributed to the seasonal modulation of ENSO’s SST growth rate. In contrast, in a significant portion of CMIP models, phase-locking is co-determined by the seasonal modulations of both SST growth and phase-transition rates. Further study of the joint effects of SST growth and phase-transition rates suggests that for simulating realistic winter peak ENSO phase-locking with the right dynamics, climate models need to have four key factors in the right combination: (1) correct phase of SST growth rate modulation peaking at the fall; (2) large enough amplitude for the annual cycle in growth rate; (3) amplitude of semi-annual cycle in growth rate needs to be small; and (4) amplitude of seasonal modulation in SST phase-transition rate needs to be small.

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

Quasiparticle random-phase approximation calculations for M 1 transitions with the noniterative finite-amplitude method and application to neutron radiative capture cross sections

Here, we derive the equations of quasiparticle random-phase approximation (QRPA) based on the finite amplitude method with the Hartree-Fock+Bardeen–Cooper–Schrieffer (HF+BCS) single-particle states, and calculate the magnetic dipole (M⁢1) transition for deformed gadolinium isotopes. Our QRPA calculation shows both large spin-flip transitions in the 5 to 10 MeV excitation energy and the low energy orbital transition that would correspond to the M⁢1 scissors mode observed in nuclear experiments. Then, we calculate neutron capture reactions based on the statistical Hauser-Feshbach theory with the photoabsorption cross sections of even-even nuclei given by QRPA. We find that the capture cross section is enhanced due to the contribution from the low energy M⁢1 transition although the calculated capture cross section still underestimates the experimental data. This issue in the calculated capture cross section could be improved by uncertainties of the low energy E⁢1 transition neglected in our QRPA calculation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Phases of Theories with ℤ 𝑁 1-Form Symmetry, and the Roles of Center Vortices and Magnetic Monopoles

We analyze the phases of theories having a microscopic ℤ 𝑁 1-form symmetry, starting with a topological BF theory and deforming it so that only the microscopic symmetry is preserved. These theories have a well-defined notion of confinement, prototypical examples being pure SU⁡(𝑁) and ℤ 𝑁 gauge theories in the continuum and on the lattice. Our analysis shows that the generic phases are in 𝑑 = 2, only the confined phase; in 𝑑 = 3, both the confined phase and the topological BF phase; and in 𝑑 = 4, the confined phase, the topological BF phase, and a Coulomb phase. We construct a ℤ𝑁 lattice gauge theory with a deformation that, surprisingly, produces up to (𝑁−1) photons. We give an interpretation of these findings in terms of the behaviors of two competing drivers of confinement—magnetic monopoles and center vortices—and conclude that proliferation of center vortices is necessary but insufficient for confinement, while proliferation of magnetic monopoles is both necessary and sufficient.

Color confinement↗
DOE OSTI2025

Dissolved Oxygen Redox as the Source of Hydrogen Peroxide and Hydroxyl Radical in Sonicated Emulsive Water Microdroplets

Sonicated emulsive water microdroplets (SEWMs) accelerate and enable a variety of catalyst-free chemical transformations. However, significant unanswered questions remain regarding the chemical intermediates they form and their possible redox origin. In this study, we identified dissolved O 2 as the primary originator of reactive oxygen species (ROS) such as OH• and H 2 O 2. We uncovered the role of dissolved O 2 redox by using a combination of microelectrochemical methods to detect H 2 O 2 , isotopic methods to identify the source of H 2 O 2 , and a combination of electron spin resonance and the DMPO spin trap to detect radicals such as OH•. Notably, we found that H 2 O 2 production is correlated with O 2 content via a reduction pathway enabled by a sufficiently large reducing power that can additionally generate H 2 and even perform Pb electroless deposition on Au and Cu metal substrates. Building on our findings, continuous O 2 bubbling of SEWMs showed accumulation of H 2 O 2 up to ∼88 mM in the aqueous phase within 1 h of sonication, demonstrating the scale-up promise of this method. Distinct to sonochemistry of a single phase, this study advances our understanding of the confluence of redox and chemical reaction mechanisms within SEWMs as a biphasic system. This insight paves the way for improving their reaction kinetics, yield, and selectivity, positioning these attractive redox microreactors as alternatives to traditional electrolyzers.

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

Cation Migration in Physically Paired 2D and 3D Lead Halide Perovskite Films

2D/3D interfaces provide long-term stability for the operation of metal halide perovskite solar cells. However, the cation migration under heat and light can disturb the interface and create a graded cation interface. This study has now probed the cation migration by physically pairing X 2 PbI 4 (X = butylammonium BA, oleylammonium OA, or phenethylammonium PEA) 2D film and (CH 3 NH 3 , MA)PbI 3 3D film at different temperatures and recording changes in the absorption and emission spectra. The migration of the methyl ammonium cation toward the 2D film slowly transforms the n = 1 layered phase into n = 2 and 3 layered phases. The 3D film, on the other hand, exhibits relatively small changes, as the inclusion of spacer cation has little effect on its phase. The apparent activation energy determined from the temperature-dependent cation migration kinetics (E a = 29– 50 kJ mol –1 ) indicates alkyl ammonium cations such as butyl ammonium migrate more readily into the 3D layer than aromatic cations such as PEA. As a result, the ease of migration of spacer cations between physically paired films suggests that the varied composition at the interface should be considered while evaluating the effectiveness of the 2D/3D design strategy for improving the performance of perovskite solar cells.

2D/3D perovskite↗
DOE OSTI2024

Pressure control of magnetic order and excitations in the pyrochlore antiferromagnet MgCr 2 ⁢O 4

MgCr 2 O 4 is one of the best-known realizations of the pyrochlore-lattice Heisenberg antiferromagnet. The strong antiferromagnetic exchange interactions are perturbed by small further-neighbor exchanges such that this compound may in principle realize a spiral spin liquid (SSL) phase in the zero-temperature limit. However, a spin Jahn-Teller transition below T N ≈ 13 K yields a complicated long-range magnetic order with multiple coexisting propagation vectors. We present neutron scattering and thermomagnetic measurements of MgCr 2 O 4 samples under applied hydrostatic pressure up to P = 1.7 GPa demonstrating the existence of multiple close-lying nearly degenerate magnetic ground states. We show that the application of hydrostatic pressure increases the ordering temperature by around 0.8 K per GPa and increases the bandwidth of the magnetic excitations by around 0.5 meV per GPa. We also evidence a strong tendency for the preferential occupation of a subset of magnetic domains under pressure. In particular, we show that the $k$ = (0,0,1) magnetic phase, which is almost negligible at ambient pressure, dramatically increases in spectral weight under pressure. This modifies the spectrum of magnetic excitations, which we interpret unambiguously as spin waves from multiple magnetic domains. Moreover, we report that the application of pressure reveals a feature in the magnetic susceptibility above the magnetostructural transition. We interpret this as the onset of a short-range ordered phase associated with $k$ = (0,0,1), previously not observed in magnetometry measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car–Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson’s equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. Furthermore, this is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500–1000 atoms (e.g., (H₂O)₂₅₆) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.

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

Manipulation of spin orientation via ferroelectric switching in Fe-doped Bi 2 WO 6 from first principles

Atomic-scale control of spins by electric fields is highly desirable for future technological applications. Magnetically doped Aurivillius-phase oxides present one route to achieve this, with magnetic ions substituted into the ferroelectric structure at dilute concentrations, resulting in spin-charge coupling. However, there has been minimal exploration of the ferroelectric switching pathways in this materials class, limiting predictions of the influence of an electric field on magnetic spins in the structure. Here, we determine the ferroelectric switching pathways of the end member of the Aurivillius phase family, Bi 2 WO 6 , using a combination of group theoretic analysis and density functional theory calculations. We find that in the ground state P2 1 ab phase, a two-step switching pathway via C2 and Cm intermediate phases provides the lowest energy barrier. Considering iron substitutions on the W site in Bi 2 WO 6 , we determine the spin easy axis. Finally, by tracking the change in spin directionality during ferroelectric switching, we find that a 90° switch in the polarization direction leads to a 112° reorientation of the spin easy axis. The low-symmetry crystal-field environment of Bi 2 WO 6 and magnetoelastic coupling on the magnetic dopant provide a route to spin control via an applied electric field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗