SEARCH · Search NASA

Results for “first-principles phonon calculations”

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.

SourceResearch connectionFrom yearThrough year

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

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Machine learning approach for vibronically renormalized electronic band structures

Here, we present a machine learning (ML) method for efficient computation of vibrational thermal expectation values of physical properties from first principles. Our approach is based on the nonperturbative frozen phonon formulation in which stochastic Monte Carlo algorithm is employed to sample configurations of nuclei in a supercell at finite temperatures based on a first-principles phonon model. A deep-learning neural network is trained to accurately predict physical properties associated with sampled phonon configurations, thus bypassing the time-consuming ab initio calculations. To incorporate the point-group symmetry of the electronic system into the ML model, group-theoretical methods are used to develop a symmetry-invariant descriptor for phonon configurations in the supercell. We apply our ML approach to compute the temperature dependent electronic energy gap of silicon based on density functional theory (DFT). We show that, with less than a hundred DFT calculations for training the neural network model, an order of magnitude larger number of sampling can be achieved for the computation of the vibrational thermal expectation values. Our work highlights the promising potential of ML techniques for finite temperature first-principles electronic structure methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Density functional theory-based surrogate kinetic models for heterogeneous reactions of hydrocarbon intermediates on silicon carbide

The increasing demand for high-performance materials in advanced technologies highlights the importance of achieving a fundamental understanding and potential control of silicon carbide (SiC) deposition processes. However, existing models often lack sufficient theoretical detail, relying heavily on empirical data and offering limited predictive capability. In particular, the complex surface chemistry governing SiC growth remains poorly understood. This study addresses these challenges by employing density functional theory (DFT) to investigate key heterogeneous reactions involving hydrocarbon intermediates on SiC surfaces, including dehydrogenation, hydrogenation, and carbon deposition. Transition state searches were conducted to identify reaction pathways and energy barriers. While first-principles calculations offer high accuracy, they are computationally intensive. To extend the utility of these first-principles results, vibrational analyses were performed using phonon-based statistical thermochemistry to compute temperature-dependent reaction rates which were used to develop Arrhenius-type surrogate kinetic models. Furthermore, the resulting framework provides a more rigorous, physically grounded basis for integrating atomistic insights into continuum-scale modeling, ultimately enabling improved prediction and optimization of SiC film growth in high-performance material systems.

Density Functional Theory↗

DOE OSTI2026

Machine Learning an Ab-Initio Based Bond-Order Potential for Bismuthene

Bismuthene is a heavy 2D material whose strong spin–orbit coupling and recently observed single-element ferroelectricity have intensified interest in its structural, vibrational, and transport properties. Accurate modeling of these behaviors requires a short-range interatomic potential that can reproduce the underlying bonding physics at a fraction of the computational cost of first-principles methods. However, such a potential is currently unavailable. Here, in this work, we construct a Tersoff bond-order potential for β-bismuthene using a reinforcement-learning framework that integrates a continuous Monte Carlo Tree Search with a simplex-based local optimizer. The optimized parameter sets reproduce first-principles lattice constants, cohesive energy, the equation of state, elastic constants, and phonon dispersion. We validate the models by performing thermal-conductivity calculations and uniaxial fracture simulations our findings confirm the reliability of the resulting models across multiple thermomechanical regimes. Comparison of the three best solutions reveals how differences in pairwise interactions, angular terms, and bond-order behavior govern phonon features and mechanical responses. We demonstrate an interpretable and computationally efficient potential for bismuthene and demonstrate a general reinforcement-learning strategy for developing bond-order models in emerging 2D materials.

deformation↗

DOE OSTI2025

Fermi-surface driven frustrations in charge ordered kagome metal LuNb 6 ⁢Sn 6

The charge density wave (CDW), a translational symmetry breaking electronic liquid, plays a pivotal role in correlated quantum materials, such as high-𝑇 𝑐 superconductors and topological semimetals. Recently, CDWs that possibly intertwine with superconductivity and magnetism are observed in various kagome metals. However, the nature of CDWs and the role of the Fermi-surface (FS) topology in these materials remain an unresolved challenge. In this paper, we reveal the formation of CDWs in the newly discovered kagome metal LuNb 6 ⁢Sn 6 . Here, we observe a ‘‘yield sign’’-like hollow triangular diffuse scattering pattern and nearly complete softening of a flat optical phonon band near 𝐐 𝐻 = (1/3,1/3,1/2). The scattering intensity of the diffuse scattering displays divergent behavior as decreasing temperature until 𝑇 CDW = 70 K, where a competing CDW at 𝑄 CDW = (1/3,1/3,1/3) emerges. Combined with scanning tunneling microscopy/spectroscopy and first-principles calculations, our observations support a Fermi-surface driven frustration that suppresses the leading CDW instability at 𝐐 𝐻 . This frustration is relieved by the emergence of a subleading CDW at 𝑄 CDW . These results provide insights into the interplay between the Fermi surface, strong electron-phonon coupling, and charge density wave formation in quantum materials.

Yang, F. Z. [Oak Ridge National Laboratory (ORNL),↗

DOE OSTI2026

Atomically Resolved Acoustic Dynamics Coupled with Magnetic Order in a Van der Waals Antiferromagnet

Magnetoelastic coupling in van der Waals (vdW) magnetic materials enables a unique interplay between the spin and lattice degrees of freedom. Characterizing the elastic responses with atomic and femtosecond resolution across the magnetic transition is essential for guiding the design of magnetically tunable actuators and strain-mediated spintronic devices. Here, ultrafast X-ray diffraction employed at a free-electron laser reveals that the atomic displacements, wave vectors, and dispersion relations of acoustic phonon modes in a vdW antiferromagnet FePS3 are coupled with the magnetic order, by tracking both in-plane and out-of-plane Bragg peaks upon optical excitation across the Néel temperature (TN). One transverse mode shows that a quasi-out-of-plane atomic displacement undergoes a significant directional change across TN. Its quasi-in-plane wave vector is derived by comparing the measured sound velocity and the first-principles calculations. The other transverse mode is an interlayer shear acoustic mode whose amplitude is strongly enhanced in the antiferromagnetic phase, exhibiting eight times stronger amplitude than the longitudinal acoustic mode below TN. The atomically resolved characterization of acoustic phonon dynamics that couple with magnetic ordering opens opportunities for harnessing unique magnetoelastic coupling in vdW magnets on ultrafast timescales.

coherent acoustic phonons↗

DOE OSTI2025

Time-domain theory of transient heat conduction in the local limit

Ultrafast and nanoscale heat conduction demands a unified theoretical framework that rigorously bridges macroscopic transport equations with microscopic material properties derived from statistical physics. Existing empirical generalizations of Fourier's law often lack a solid microscopic foundation, failing to connect observed non-Fourier behavior with underlying atomic-scale mechanisms. In this work, we present a time-domain theory of transient heat conduction rooted in Zwanzig's statistical theory of irreversible processes. Central to this framework is the time-domain transport function $\overleftrightarrow{𝑍}$⁡(𝑡) defined through equilibrium time-correlation functions of heat fluxes. This function generalizes the conventional concept of steady-state thermal conductivity, governing the transition of conduction dynamics from onset second sound type wave propagation at finite speeds to diffusion-dominated behavior across broad temporal and spatial scales. Unlike phonon hydrodynamic models that rely on mesoscopic constructs such as phonon drift velocity, our approach provides a quantitative and microscopic description of intrinsic memory effects in transient heat fluxes and applies universally to bulk materials at any temperature or length scale. By integrating atomistic-scale first-principles calculations with continuum-level macroscopic equations, this framework offers a robust foundation for numerical simulations of transient temperature fields. Furthermore, it facilitates the interpretation and design of transient thermal grating experiments using nanometer-scale heat sources and ultrafast laser systems in the extreme ultraviolet and x-ray wavelength ranges, advancing our understanding of heat dissipation dynamics.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2025

Temperature-Dependent Mechanical Properties of Ni-Based Concentrated Alloys: Insights from First-Principles Calculations

The present work focuses on predicting temperature-dependent mechanical properties of Ni-based concentrated alloys Ni 18 Cr 10 Co 10 Fe 6 M 4 (abbreviated by X 44 M 4 , with M = Al, V, Mn, Fe, Nb, Mo, and W) using density functional theory (DFT). These predictions are based on shear (plastic) and elastic deformations, utilizing the special quasirandom structure (SQS), the phonon-based quasiharmonic approach (QHA), and the quasistatic approach. The resulting properties include coefficient of thermal expansion via QHA, ideal shear strength (τ IS ), and stable and unstable stacking fault energies (γ SF and γ US ) through pure alias shear deformation, and elastic constants (c ij ), bulk modulus (B 0 ), and shear modules (G 0 ) via elastic deformation. Notably, predicting accurate γ SF is challenging due to uncertainties that can exceed the γ SF values. τ IS and γ US exhibit a strong linear relationship, enabling the accurate prediction of γ US based on the precisely determined τ IS . All mechanical properties of X 44 M 4 decrease with increasing temperature, except for some γ SF cases such as X 44 M 4 with M = V, Mn, Fe, Mo, and W. Among the X 44 M 4 alloys, X 44 Nb 4 exhibits the lowest τ IS , γ US , and G 0 values, and the highest B 0 /G 0 ratio, while X 44 Mn 4 has the lowest B 0 and B 0 /G 0 ratio. We found that volume is a crucial descriptor for understanding and modeling mechanical properties (except B0 and maybe also γ SF ) affected by alloying elements and temperature. Ni-based dilute alloys (e.g., Ni 11 M 1 and Ni 31 M 1 ) and concentrated alloys (e.g., X 44 M 4 ) show similar trends in mechanical properties influenced by alloying elements and temperature, simplifying the analysis and design of Ni-based alloys.

Elastic properties↗

DOE OSTI2026

Understanding polaronic transport in complex oxides by combining precise synthesis and first-principles many-body theory

In complex oxides, charge carriers often couple strongly with lattice vibrations to form polarons–entangled electron–phonon quasiparticles whose transport properties remain difficult to characterize. Experimental access to intrinsic polaronic transport requires ultraclean samples, while theoretical description demands methods beyond low-order perturbation theory. Here, we show a predictive theory–experiment workflow to study polaron transport in complex oxides. Focusing on a prototypical polaronic oxide, anatase TiO 2 , we combine growth of high-quality oxygen-vacancy-doped films using hybrid molecular beam epitaxy with a first-principles electron–phonon diagrammatic Monte-Carlo (FEP-DMC) framework recently developed for accurate polaron predictions. Our films exhibit record-high electron mobility for anatase TiO 2 , in excellent agreement with FEP-DMC calculations conducted prior to experiment, which predict a room-temperature mobility of 45 ± 15 cm −2 V −1 s −1 and a mobility-temperature scaling of μ ∝ T −1.9 ± 0.077 . Microscopic analysis using scanning transmission electron microscopy and x-ray photoelectron spectroscopy reveals the role of oxygen vacancies in modulating transport at lower temperatures. FEP-DMC further provides quantitative insight into polaron formation energy, phonon cloud distribution, lattice distortion around the polaron, and the polaronic contribution to mobility. Together, these results provide a deeper microscopic understanding of large-polaron transport in a complex oxide and provide the blueprint to characterize other polaronic materials.

anatase TiO2↗

DOE OSTI2026

ALD-Derived WO 3– x Leads to Nearly Wake-Up-Free Ferroelectric Hf 0.5 Zr 0.5 O 2 at Elevated Temperatures

Breaking the memory wall in advanced computing architectures will require complex 3D integration of emerging memory materials such as ferroelectrics─either within the back-end-of-line (BEOL) of CMOS front-end processes or through advanced 3D packaging technologies. Achieving this integration demands that memory materials exhibit high thermal resilience, with the capability to operate reliably at elevated temperatures, such as 125°C, due to the substantial heat generated by front-end transistors. However, silicon-compatible HfO 2 -based ferroelectrics tend to exhibit antiferroelectric-like behavior in this temperature range, accompanied by a more pronounced wake-up effect, posing significant challenges to their thermal reliability. Here, we report that by introducing a thin tungsten oxide (WO 3–x ) layer─known as an oxygen reservoir─and carefully tuning its oxygen content, ultrathin Hf 0.5 Zr 0.5 O 2 (5 nm) films can be made robust against the ferroelectric-to-antiferroelectric transition at elevated temperatures. This approach not only minimizes polarization loss in the pristine state but also effectively suppresses the wake-up effect, reducing the required wake-up cycles from 10 5 to only 10 at 125°C, a qualifying temperature for back-end memory integrated with front-end logic, as defined by the JEDEC standard. First-principles density functional theory (DFT) calculations reveal that WO 3 enhances the stability of the ferroelectric orthorhombic phase (o-phase) at elevated temperatures by increasing the tetragonal-to-orthorhombic phase energy gap and promoting favorable phonon mode evolution, thereby supporting o-phase formation under both thermodynamic and kinetic constraints.

36 MATERIALS SCIENCE↗

DOE OSTI2026

High-Performance Near-Infrared Quantum Emission from Color Centers in hBN

Color centers hosted in hexagonal boron nitride (hBN) have emerged as a highly promising platform for single-photon emission and spin-photon technologies relevant to quantum communication and quantum networking. As a wide bandgap van der Waals material, hBN can host optically active quantum defects across a broad spectral range. Here, we demonstrate a simple and scalable oxygen-plasma process that reproducibly creates single quantum emitters in hBN with blinking-free zero-phonon lines (ZPLs) spanning near-infrared (NIR) from 700 up to 971 nm. These emitters combine MHz-level brightness, single-photon purity up to 99.9%, and ultranarrow cryogenic line widths down to 2.7 GHz under quasi-resonant excitation, placing them in a particularly attractive regime for quantum photonics. Photostability measurements further reveal resistance to photobleaching, subnanometer spectral stability over long time scales, and near-shot-noise-limited intensity fluctuations. Analysis of the phonon sidebands shows weak vibronic coupling and ZPL-dominated emission, with Debye–Waller factors approaching 50%. Control experiments together with elemental mapping support oxygen incorporation as a necessary ingredient in activating the NIR emitter population, while first-principles calculations identify O N V N and O N V N H as the leading defect candidates. These results establish a high-performance NIR quantum-emitter platform in hBN for free-space quantum networking and future integrated quantum-photonic architectures.

2D materials↗

DOE OSTI2025

First-Principles Studies on Sc 2 RuZ (Z = Si, Ge, Sn) Inverse Heusler Alloys: Structural, Electronic, and Transport Properties

The continuous demand for efficient, nontoxic, and thermally stable materials for room-temperature energy conversion motivates the exploration of novel thermoelectric systems beyond the traditional magnetic Heusler alloys. While full and half-Heusler compounds, especially Co-, Ni-, and Mn-based systems, have demonstrated promising thermoelectric properties, their typically high operating temperatures and magnetic complexities limit their applicability in ambient thermal management. In this context, we investigate whether Sc-based inverse Heusler alloys can offer a viable nonmagnetic alternative with competitive thermoelectric performance. In this work, we perform a systematic first-principles study of the inverse Heusler compounds Sc 2 RuZ (Z = Si, Ge, Sn), focusing on their structural, electronic, mechanical, and thermodynamic-thermoelectric properties. Density Functional Theory (DFT) was employed to compute optimized lattice structures and band dispersion, while dynamical stability was assessed via phonon calculations. Thermoelectric transport coefficients, including Seebeck coefficient, electrical conductivity, and thermal conductivity, were estimated using the semiclassical Boltzmann transport theory within the constant relaxation time approximation. Our results show that all Sc 2 RuZ compounds are thermodynamically stable semiconductors with indirect band gaps of 0.12–0.16 eV and exhibit high elastic moduli, especially Sc 2 RuSn, which demonstrates superior stiffness and incompressibility. Importantly, all compounds display promising room-temperature thermoelectric characteristics, including high Seebeck coefficients and power factors. These findings reveal that Sc 2 RuZ alloys represent a rare class of stable, nonmagnetic inverse Heusler semiconductors with intrinsic thermoelectric potential at room temperature, unlike many existing Heusler systems optimized for spintronics or high-temperature operation. This work expands the known design space for Heusler-based thermoelectrics and offers a theoretical basis for experimental realization of efficient, low-temperature, nonmagnetic thermoelectric materials.

alloys↗

DOE OSTI2026

Pressure-tuned plethora of ferroelectric phases in CuInP 2 S 6

Layered chalcogenides like CuInP 2 S 6 are room temperature ferroelectrics. Modest compression even increases the electric polarization, raising questions about the origin of this unusual trend and other properties under pressure. In this work, we combine synchrotron-based infrared and Raman scattering spectroscopies, x-ray diffraction, and resistivity measurements with first-principles calculations of the lattice dynamics and energy landscape to unravel the influence of pressure on symmetry, polarization, and metallicity. We uncover a remarkable series of phase transitions across a series of polar space groups: monoclinic Cc → trigonal P31c (prismatic sulfur) → P31c (octahedral sulfur). True metallicity develops above 63 GPa, significantly higher than in related MPS 3 materials (M = Mn, Co, Fe, Ni), offering a picture of competing states of matter that is different than previously supposed. Detailed examination of pressure trends within the Cc phase also reveals phonon lifetime changes and streaking of satellite x-ray peaks that correlate with the maximum polarization. We discuss these tendencies in terms of Cu + ion migration, phase formation, and the overall energy landscape. Our findings place the high pressure behavior of CuInP 2 S 6 on a firm foundation and pave the way for the development of structure-property relations in this family of complex chalcogenides.

ferroelectric phases↗

DOE OSTI2025

Intercalated chalcogenides Fe 1/4 TaS 2 and Fe 1/3 TaS 2 under extreme pressure—temperature conditions

Intercalation is an important strategy for enhancing functionality in complex chalcogenides. This is because layered materials can be endowed with intriguing properties by filling the van der Waals gap with various ions and molecules, which, in addition to their unique chemistry, break symmetry in new ways. In order to explore the properties of intercalated metal dichalcogenides under extreme pressure-temperature conditions, we combine diamond anvil cell techniques, Raman scattering spectroscopy, magnetic susceptibility, and first-principles calculations to reveal the pressure–temperature phase diagram of Fe 1/3 TaS 2 . This system hosts a compression-driven structural phase transition to a polar state as well as remnant charge density wave signatures deriving from the host lamella of the 2H-parent compound. We also explore the role of different A-site patterns and determine that, by comparison, Fe 1/4 TaS 2 is soft and flexible due to the lower metal density inside the van der Waals gap. These findings open the door to entirely new states of matter with exciting property combinations, including metallicity, polarity, chirality, and altermagnetism – depending upon the conditions – that can support a wide range of spintronics and phononics applications.

Materials science↗

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 OSTI2026

Carbon-nitrogen complex as an alternative to the 𝑇 center in Si

A first-principles study of a carbon-nitrogen (CN) impurity complex in silicon is presented as an isoelectronic alternative to the T center [(CCH) Si ]. The latter has been pursued for applications in quantum information science, yet its sensitivity to the presence of hydrogen is still problematic. Our proposed complex has no hydrogen, thereby eliminating this issue. First, we show that the CN complex is stable against decomposition into substitutional and interstitial defects. Next, we show that due to being isoelectronic to the T center, the CN complex has a similar electronic structure, and therefore could be used in similar applications. We assess several low-energy configurations of the CN complex, finding (CN) Si to be stable and have the largest Debye-Waller factor. We predict a zero-phonon line (ZPL) of 828 meV (in the telecom S-band) and a radiative lifetime of 4.2 μs, comparable to the T center. Due to the presence of a bound exciton, choice of the exchange-correlation functional and also supercell-size scaling of the ZPL and transition dipole moment require special scrutiny; we rigorously justify our extrapolation schemes that allow computing values in the dilute limit.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Leveraging Polymorphism in YbCuBi to Map Transport and Elastic Properties

AMX Zintl compounds with the hexagonal ZrBeSi structure have gained significant attention for their remarkable vacancy tolerance and low thermal conductivity. Their 2D honeycomb sublattice, composed of M–X covalent bonds, is believed to contribute to high anharmonicity and unusual thermal transport properties. In this study, we explore the temperature-dependent polymorphism of YbCuBi as a model system to investigate the relationship between the structure and elastic and thermal transport properties in AMX Zintls. YbCuBi undergoes a structural transition from the “flat” Cu–Bi layers in the ZrBeSi structure to corrugated layers in the LiGaGe structure below 410 K, resulting in a distortion of its centrosymmetric structure. To probe the effects of this crystallographic transition, we employ inelastic neutron scattering and temperature-dependent resonant ultrasound spectroscopy. These experimental findings, coupled with first-principles calculations and thermal conductivity measurements, allow us to elucidate a direct relationship between corrugation of the honeycomb lattice and the observed changes in elastic and thermal transport properties. These insights can be extended to other Zintl phases with similar structure types, providing a platform for the rational design of functional materials with tailored thermal properties.

Chemical structure↗

DOE OSTI2025

First-principles computations of the Stark shift of a defect-bound exciton: The case of the T center in silicon

The T center in silicon has recently drawn a lot of attention for its potential in quantum information science. The sensitivity of the zero-phonon line (ZPL) to electrical field was recently investigated by a combination of different experimental methods but there are still few first principles studies on the Stark shift of the T center. Dealing with the defect-bound exciton nature of the excited state is particularly challenging using density functional theory because of the large spatial delocalization associated with the wavefunction. Here, in this work, we tackle this issue by performing a convergence study over the supercell size. We obtain an exciton binding energy of 28.5 meV, in good agreement with experimental results. We then calculate the Stark shift through the dipole moment change of the ZPL transition of the T center using the modern theory of polarization formalism and find a modest linear coefficient of Δ⁢𝜇=0.79⁢D along X and Δ⁢𝜇=0.03⁢𝐷 along Y. We discuss our results in light of the recent experimental measurements of the Stark shift. Our analysis suggests that bound-exciton defects could be particularly sensitive to local field effect as a result of their large spatial extent.

color centers↗

DOE OSTI2025

𝐴𝑏 initio density-matrix approach to exciton coherence: Phonon scattering, Coulomb interactions, and radiative recombination

Relaxation processes following light excitation in semiconductors are key in materials-based quantum technology applications. These processes are broadly studied in atomically thin transition-metal dichalcogenides, quasi-two-dimensional excitonic semiconductors in which atomistic design allows for tunable excited-state properties, such as relaxation lifetimes and photoinduced coherence. In this work, we present a density-matrix-based approach to compute exciton relaxation within a many-body ab initio perspective. We expand our previously developed Lindblad density-matrix formalism to capture multichannel electron-hole pair relaxation processes, including phonon and Coulomb scattering as well as radiative recombination, and we study their effect on the time-resolved excited-state propagation. Using monolayer MoSe 2 as a prototypical example, we examine many-body effects on the time-dependent dynamics of photoactive excitations, exploring how the electron-hole pair interactions are reflected in variations of the excitation energy, spectral signature, and state coherence. In conclusion, our method supplies a detailed understanding of exciton relaxation mechanisms in realistic materials, offering a previously unexplored pathway to study excited-state dynamics in semiconductors from first principles.

Band structure methods↗