Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first-principles”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Controlling topology through targeted composite symmetry manipulation in magnetic systems

The possibility of selecting magnetic space groups by orienting the magnetization direction or tuning magnetic orders offers a vast playground for engineering symmetry-protected topological phases in magnetic materials. In this work, we study how selective tuning of symmetry and magnetism can influence and control the resulting topology in a two-dimensional magnetic system, and we illustrate such a procedure in the ferromagnetic monolayer MnPSe 3 . Density functional theory calculations reveal a symmetry-protected accidental semimetallic (SM) phase for out-of-plane magnetization, which becomes an insulator when the magnetization is tilted in-plane, reaching band-gap values close to 100 meV. We identify an order-2 composite antiunitary symmetry and threefold rotational symmetry that induce the band crossing, and we classify the possible topological phases using symmetry analysis, which we support with tight-binding and k · p models. Breaking of inversion symmetry opens a gap in the SM phase, giving rise to a Chern insulator. We demonstrate this explicitly in the isostructural Janus compound Mn 2 ⁢P 2 ⁢S 3 ⁢Se 3 , which naturally exhibits Rashba spin-orbit coupling that breaks inversion symmetry. Our results map out the phase space of topological properties of ferromagnetic transition-metal phosphorus trichalcogenides, and they demonstrate the potential of the magnetization-dependent metal-to-insulator transition as a spin switch in integrated two-dimensional electronics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Thermal disorder and phonon softening in the ferroelectric phase transition of lead titanate

We report a molecular dynamics study of ab initio quality of the ferroelectric phase transition in crystalline PbTi⁢O3. We model anharmonicity accurately in terms of potential energy and polarization surfaces trained on density functional theory data with modern machine learning techniques. Our simulations demonstrate that the transition has a strong order-disorder character, in agreement with diffraction experiments, and provide fresh insight into the approach to equilibrium across the phase transition. We find that the emergence and disappearance of the macroscopic polarization is driven by dipolar switching at the nanometer scale. We also computed the infrared optical absorption spectra in both the ferroelectric and the paraelectric phases, finding good agreement with the experimental Raman frequencies. Often, the almost ideal displacive character of the soft mode detected by Raman scattering in the paraelectric phase has been contrasted with the order-disorder character of the transition suggested by diffraction experiments. We settle this issue by showing that the soft mode coexists with a strong Debye relaxation associated with thermal disordering of the dipoles. The Debye relaxation feature is centered at zero frequency and appears near the transition temperature in both the ferroelectric and the paraelectric phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Melting curves of atomic hydrogen and deuterium calculated using path-integral Monte Carlo

We calculate the melting line of atomic hydrogen and deuterium up to 900 GPa with path-integral Monte Carlo using a machine-learned interatomic potential. We improve upon previous simulations of melting by treating the electrons with reptation quantum Monte Carlo, and by performing solid and liquid simulations using isothermal-isobaric path-integral Monte Carlo. Here, the resulting melting line for atomic hydrogen is higher than previous estimates. There is a small but resolvable decrease in the melting temperature as pressure is increased, which can be attributed to quantum effects.

08 HYDROGEN↗
DOE OSTI2025

Role of magnetic and structural symmetry breaking in forming the Mott insulating gap in Nb 3 ⁢Cl 8

The α-phase of the gapped insulator Nb 3 Cl 8 has recently emerged as the long-sought critical testing bed for examining the importance of strong interelectronic correlation vs symmetry breaking in understanding insulation of such Mott compounds. Structural symmetry breaking detected by density functional theory (DFT) energy lowering (such as dimer formation, disproportionation, or Jahn-Teller distortions) explains insulation in both d-electron Mott-like systems and in non-d-electron cases without recourse to strong correlation. Yet, in Nb 3 Cl 8 , structural symmetry breaking alone (viz. formation of Nb trimers) fails to explain insulation, leading instead to a partially occupied metallic flat band, in contrast with experimental observations. We examine the role of magnetic symmetry breaking, noting that Nb 3 Cl 8 is an observed paramagnet (not an antiferromagnet), thus potentially carrying also short-range ordered magnetic moments. Describing the latter as a polymorphous distribution of nonzero local moments with total zero net magnetization is demonstrated to lower the DFT total energy, while gapping the system without recourse to strong correlation or long-range magnetic order. This suggests that degeneracy removal by symmetry breaking in mean-field-like approaches—either structural, or magnetic, or both—can reduce or eliminate the need for strong correlation, allowing the use of DFT for such Mott systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Weak collisionless shocks mediated by ion gyroviscosity

Collisionless shocks are ubiquitous in space and astrophysical plasmas, and they are essential dynamical features of these systems. Lacking Coulomb collisions, these shocks are mediated by the anomalous dissipation provided by nonlinear plasma instabilities. By numerically resolving the structure of a steady-state, ion gyroviscous shock, we show that ion gyroviscosity, alone, can produce weak (M≲1.1, where M is the sonic Mach number) shocks in a collisionless, magnetized plasma. We emphasize that this effect does not require an appeal to plasma microturbulence. Moreover, while most collisionless systems may be unsuitable to support purely gyroviscous shocks, we argue that gyro-viscous heating may be an overlooked mechanism, generally; and it may be a key driver within magnetohydrodynamic shocks at large. In conclusion, representative examples include the plasma environments produced on the plasma liner experiment and the magnetized liner inertial fusion platforms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Domain Wall Reactions in Multiple-Order Parameter Ferroelectrics

The motion of domain walls is crucial for ferroelectric switching. Conventionally, the switching dynamics is believed to be determined by the motion of one or a few low-energy domain wall types of dominant population. Here, we challenge this conventional idea in multiple-order-parameter ferroelectrics. Using hafnia as example, we show that the multiple-order-parameter nature not only provides various mobile domain walls and defectlike immobile domain walls, but also enables the domain wall reactions. In analogy with chemical reactions where substances react to form new substances, domain walls could also react to form other domain walls during switching. We identify several elementary domain wall reaction types including synthesis, decomposition, and exchange reactions. Domain walls are continually changed by these reactions during switching so that the switching behavior reflects the statistical average of many domain wall types with distinct mobility and stability. These reactions also lead to phenomenon like remanent nuclei and defect-site nucleation that facilitate switching and lower coercive field. Finally, the concept of domain wall reaction is not limited to hafnia but can be generalized to any multiple-order-parameter ferroelectric. As a result, this Letter conceptualizes domain wall reaction, expands the theory of ferroelectric switching, and suggests a practical way for defect engineering to control switching behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nonperturbative Guiding Center Model for Magnetized Plasmas

Perturbative guiding center theory adequately describes the slow drift motion of charged particles in the strongly magnetized regime characteristic of thermal particle populations in various magnetic fusion devices. However, it breaks down for particles with large-enough energy. Here, we report on a data-driven method for learning a nonperturbative guiding center model from full-orbit particle simulation data. We show the data-driven model significantly outperforms traditional asymptotic theory in magnetization regimes appropriate for fusion-born α particles in stellarators, thus opening the door to nonperturbative guiding center calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Weyl semimetal phases and intrinsic spin-Hall conductivity in SbAs ordered alloys

Here, using density functional theory calculations, we investigated possible Weyl semimetal (WSM) phases in antimony arsenide ordered alloys Sb 1-x ⁢As x (x=0, 1/6, 1/3, 1/2, 2/3, 5/6, 1). We find WSM phases for all As compositions of Sb 1-x ⁢As x with broken inversion symmetry, in contrast to Bi 1-x⁢ Sb x where only compositions x=1/2 and 5/6 were predicted to exhibit WSM phases. The WSM phases in Sb 1-x ⁢As x are characterized by the presence of 12 Weyl points, located within 55 meV from the Fermi level in the case of x = 1/2. The robust spin-orbit coupling strength and Berry curvature in these alloys produce large spin-Hall conductivity in the range of 176–602 ($\hslash$/e)(S/cm), comparable to that in the BiSb alloys. Finally, Sb 0.5 ⁢As 0.5 is predicted to be almost lattice matched to GaAs(111), with the Fermi level within the gap of the semiconductor, facilitating growth and characterization, and thus, offering promising integration with conventional semiconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Three-dimensional higher-order saddle-point-induced flatbands in Co-based kagome metals

The saddle point (Van Hove singularity) exhibits a divergent density of states in two-dimensional systems, leading to fascinating phenomena such as strong correlations and unconventional superconductivity, yet it is seldom observed in three-dimensional (3D) systems. In this work we find two types of 3D higher-order saddle points (HOSPs) in emerging 3D kagome metals YbCo 6⁢ Ge 6 and MgCo 6 ⁢Ge 6 . Both HOSPs exhibit a singularity in their density of states, which is significantly enhanced compared to the ordinary saddle point. The HOSP near the Fermi energy generates a flatband extending a large area in the Brillouin zone, potentially amplifying the correlation effect and fostering electronic instabilities. Two types of HOSPs exhibit distinct robustness upon element substitution and lattice distortions in these kagome compounds. Our work paves the way for engineering exotic band structures, such as saddle points and flatbands, and exploring interesting phenomena in Co-based kagome materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Effective many-body interactions in reduced-dimensionality spaces through neural network models

Accurately describing properties of challenging problems in physical sciences often requires complex mathematical models that are unmanageable to tackle head on. Therefore, developing reduced-dimensionality representations that encapsulate complex correlation effects in many-body systems is crucial to advance the understanding of these complicated problems. However, a numerical evaluation of these predictive models can still be associated with a significant computational overhead. To address this challenge, in this paper we discuss a combined framework that integrates recent advances in the development of active-space representations of coupled cluster (CC) downfolded Hamiltonians with neural network approaches. The primary objective of this effort is to train neural networks to eliminate the computationally expensive steps required for evaluating hundreds or thousands of Hugenholtz diagrams, which correspond to multidimensional tensor contractions necessary for evaluating a many-body form of downfolded effective Hamiltonians. Using small molecular systems (the H 2 O and HF molecules) as examples, we demonstrate that training neural networks employing effective Hamiltonians for a few nuclear geometries of molecules can accurately interpolate or extrapolate their forms to other geometrical configurations characterized by different intensities of correlation effects. We also discuss differences between effective interactions that define CC downfolded Hamiltonians with those of bare Hamiltonians defined by Coulomb interactions in the active spaces. Published by the American Physical Society 2024

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Microscopic origin of temperature-dependent magnetism in spin-orbit-coupled transition metal compounds

A few 4 d and 5 d transition metal compounds with various electron fillings were recently found to exhibit magnetic susceptibilities χ and magnetic moments that deviate from the well-established Kotani model. This model has been considered for decades to be the canonical expression for descriing the temperature dependence of magnetism in systems with nonnegligible spin-orbit coupling effects. In this paper, we uncover the origin of such discrepancies and determine the applicability and limitations of the Kotani model by calculating the temperature dependence of the magnetic moments of a series of 4 d (Ru-based) and 5 d (W-based) systems at different electron fillings. For this purpose, we perform exact diagonalization of -derived relativistic multiorbital Hubbard models on finite clusters and compute their magnetic susceptibilities. Comparison with experimentally measured magnetic properties indicates that contributions such as a temperature-independent χ 0 background, crystal field effects, Coulomb and Hund's couplings, and intersite interactions—not included in the Kotani model—are especially crucial for correctly describing the temperature dependence of χ and magnetic moments at various electron fillings in these systems. Based on our results, we propose a generalized approach beyond the Kotani model to accurately describe their magnetism. Published by the American Physical Society 2025

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Measuring the Debye energy in superconductors via two-electron photoemission spectroscopy

In this work, we demonstrate theoretically that double angle resolved photoemission spectroscopy (2eARPES) can directly probe the existence of Cooper pairs away from the Fermi surface, and can thus provide insight into the characteristic energy scale around the Fermi surface, the Debye energy, in which electrons are bound into Cooper pairs. To this end, we compute the photoelectron counting rate 𝑃 (2) in two different types of unconventional superconductors, a 𝑑 𝑥 2 −𝑦 2 -wave superconductor, and a topological superconductor with a broken time-reversal symmetry. We show that 𝑃 (2) provides insight into the relative strength of intra- and interband pairing in multiband systems, as well as into the spin polarization of the bands.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Anisotropy, frustration, and saddle point in the twisted kagome compound ErPdPb

The kagome lattice, with its inherent geometric frustration, provides a rich platform for exploring intriguing magnetic phenomena and topological electronic structures. In reduced-symmetry structures, such as twisted kagome systems involving rare earth elements, additional anisotropy can arise, enabling intriguing properties including spin-ice states, magnetocaloric effects, noncollinear magnetic ordering, and the anomalous Hall effect. Here, we report the synthesis of single crystals of ErPdPb, which features a twisted kagome lattice net of Er atoms within the hexagonal ZrNiAl-type structure, and we investigate its magnetic, electronic, and thermal properties. The material exhibits a highly anisotropic, correlated magnetic state below 2.7 K, as evidenced by magnetic, transport, and heat capacity measurements. Density functional theory (DFT) calculations indicate strong easy-axis anisotropy, consistent with experiment and crystal-field expectations, as well as quasi-one-dimensional bands and a spin-split saddle point at the zone center. Here, the coexistence of competing magnetic interactions along different crystallographic directions suggests an inherent degree of frustration. ErPdPb thus provides a promising platform for exploring the interplay of frustration, anisotropy, and electronic structure in a twisted kagome lattice.

Antiferromagnets↗
DOE OSTI2025

Single-crystalline orthorhombic GdAlGe as a rare-earth magnetic Dirac nodal-line metal

Crystal engineering is a method for discovering new quantum materials and phases, which may be achieved using external pressure or strain. Chemical pressure is unique in that it generates internal pressure perpetually to the lattice. As an example, GdAlSi from the rare-earth (𝑅) 𝑅⁢Al⁢𝑋 (𝑋=Si or Ge) family of Weyl semimetals is considered. Replacing Si with the larger isovalent element Ge creates sufficiently large chemical pressure to induce a structural transition from the tetragonal structure of GdAlSi, compatible with a Weyl semimetallic state, to an orthorhombic phase in GdAlGe, resulting in an inversion-symmetry-protected nodal-line metal. We find that GdAlGe hosts an antiferromagnetic ground state with two successive orderings, at 𝑇 N ⁢1=35K and 𝑇 N ⁢2=30K. In-plane isothermal magnetization shows a magnetic field induced metamagnetic transition at 6.2 T for 2 K. Furthermore, electron-hole compensation gives rise to a large magnetoresistance of ∼100% at 2 K and 14 T. Angle-resolved photoemission spectroscopy measurements and density functional theory calculations reveal a Dirac-like linear band dispersion over an exceptionally large energy range of ∼1.5eV with a high Fermi velocity of ∼10 6 m/s, a rare feature not observed in any magnetic topological materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Ferroelectric phase transition in group-IV monochalcogenides from an equivariant machine learned force field

Group-IV monochalcogenides are a class of layered ferroelectric semiconductors that have demonstrated spontaneous intrinsic polarization above room temperature. Here, in this study, we use the multi-atomic cluster expansion (MACE) machine learning architecture to train and test a force field capable of modeling the structural properties and second-order ferroelectric-to-paraelectric phase transition in a Group-IV monochalcogenide, GeSe. The model captures the double-well potential energy surface associated with the onset of macroscopic polarization in bulk GeSe within 12.5 meV/atom, as well as near-equilibrium properties like the phonon dispersion. The development of this quantitatively accurate force field enables long-time molecular dynamics simulations, which predict the critical temperature of the ferroelectric-to-paraelectric phase transition in bulk GeSe to be T c = 600 K. This study demonstrates the capabilities of equivariant force-fields to accurately describe phenomena associated with structural symmetry breaking.

ferroelectricity↗
DOE OSTI2025

Size-dependent attraction of Cu solutes to clusters formed at Ag grain boundaries

We report a size-dependent solute clustering mechanism at grain boundaries in a sputtered ultrafine-grained Ag-Cu alloy, where large Cu clusters form despite weak individual solute-solute interactions. X-ray diffraction confirms limited Cu solubility in the Ag matrix, while scanning transmission electron microscopy reveals Cu clustering at both ordinary GBs and GB junctions. Density functional theory calculations show that 12-atom Cu clusters are energetically preferred, while smaller three-atom clusters are significantly less stable. Additional calculations demonstrate a marked increase in solute-cluster attraction energy with cluster size. As a result, these findings point to a previously unrecognized pathway for grain-boundary solute clustering in immiscible systems, driven by collective solute-cluster interactions, with implications for segregation behavior and stability in nanocrystalline and ultrafine-grained alloys.

Alloys↗
DOE OSTI2025

Nonresonant two-photon x-ray absorption in Cu

We present a real-space Green's function theory and calculations of two-photon x-ray absorption (TPA). Our focus is on nonresonant 𝐾-shell TPA in metallic Cu, which has been observed experimentally at intense x-ray free electron laser (XFEL) sources. The theory is based on an independent particle Green's function treatment of the Kramers-Heisenberg equation and an approximation for the sum over nonresonant intermediate states in terms of a static quadrupole transition operator. XFEL effects are modeled by a partially depleted 𝑑 band. This approach is shown to give results for 𝐾-shell TPA in quantitative agreement with XFEL experiment and with a Bethe-Salpeter equation approach. Furthermore, we also briefly discuss many-body corrections and TPA sum rules.

Approximation methods for many-body systems↗
DOE OSTI2025

Engineering 2D Square Lattice Hubbard Models in 90° Twisted GeX/SnX (X =S, Se) Moiré Superlattices

Because of the large-period superlattices emerging in moiré two-dimensional (2D) materials, electronic states in such systems exhibit low energy flat bands that can be used to simulate strongly correlated physics in a highly tunable setup. While many investigations have thus far focused on moiré flat bands and emergent correlated electron physics in triangular, honeycomb, and quasi-one-dimensional lattices, tunable moiré realizations of square lattices subject to strong correlations remain elusive. Here, in this work, we propose a feasible scheme to construct moiré square lattice systems by twisting two or more layers of 2D materials in a rectangular lattice by 90°. We demonstrate the concept with twisted GeX/SnX (X =S, Se) moiré superlattices and calculate their electronic structures from first principles. We show that the lowest conduction flat band in these systems can be described by a square lattice Hubbard model with parameters which can be controlled by varying the choice of host materials, number of layers, and external electric fields. In particular, twisted double bilayer GeSe realizes a square lattice Hubbard model with strong frustration due to the next-nearest-neighbor hopping that could host unconventional superconductivity, in close analogy to the Hubbard model for copper-oxygen planes of cuprate high-temperature superconductors. The presented scheme uses 90° twisted 2D materials with rectangular unit cells as a promising platform for realizing the physical phenomena of square lattice Hubbard models, establishing a new route for studying its rich phase diagram of magnetism, charge order, and unconventional superconductivity in a highly tunable setting.

2-dimensional systems↗