Search NASA⌕ Search

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.

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

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

First-Principles Understanding of Optical Excitations within Low-Dimensional Materials

The objective of this project is to utilize first-principles computational approaches to understand optical excitations within two-dimensional monolayers and heterostructures, with the ultimate goal of designing new materials by modifying the chemical and physical structure on the nanoscale. In 2D heterostructures, the nature and movement of the exciton is determined to a large extent by the electron-electron and electron-phonon strengths within individual monolayers, as well as the coupling between layers, all of which in turn depend on chemistry, solid-state screening properties, and quantum confinement. Controlling the nature and migration properties of excitons requires an understanding of the complex relationship between electrons, phonons, defects, disorder, and dynamics. The aim of this project is to determine, by analysis of highly accurate density functional theory (DFT) and many-body perturbation theory (MBPT) calculations, how these complexities can be decomposed into simple, tunable parameters. The specific goals of this project are to 1) understand electron conductivity, absorption and transparency of monolayer semiconductors and metals; 2) develop a theory of excitons near defects in low dimensions; and 3) understand the influence of inter-layer interactions in heterostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Exciton thermalization dynamics in monolayer MoS2: A first-principles Boltzmann equation study

Understanding exciton thermalization is critical for optimizing optoelectronic and photocatalytic processes in many materials. However, it is hard to access the dynamics of such processes experimentally, especially on systems such as monolayer transition metal dichalcogenides, where various low-energy excitations pathways can compete for exciton thermalization. Here, we study exciton dynamics due to exciton-phonon scattering in monolayer MoS2 from a first-principles, interacting Green's function approach, to obtain the relaxation and thermalization of low-energy excitons following different initial excitations at different temperatures. We find that the thermalization occurs on a picosecond time scale at 300 K but can increase by an order of magnitude at 100 K. The long total thermalization time, owing to the nature of its excitonic band structure, is dominated by slow spin-flip scattering processes in monolayer MoS2. In contrast, thermalization of excitons in individual spin-aligned and spin-anti-aligned channels can be achieved within a few hundred fs when exciting higher-energy excitons. We further simulate the intensity spectrum of time-resolved angle-resolved photoemission spectroscopy experiments and anticipate that such calculations may serve as a map to correlate spectroscopic signatures with microscopic exciton dynamics.

Chan, Yang-hao↗
DOE OSTI2024

E(n)-Equivariant cartesian tensor message passing interatomic potential

Machine learning potential (MLP) has been a popular topic in recent years for its capability to replace expensive first-principles calculations in some large systems. Meanwhile, message passing networks have gained significant attention due to their remarkable accuracy, and a wave of message passing networks based on Cartesian coordinates has emerged. However, the information of the node in these models is usually limited to scalars, and vectors. In this work, we propose High-order Tensor message Passing interatomic Potential (HotPP), an E(n) equivariant message passing neural network that extends the node embedding and message to an arbitrary order tensor. By performing some basic equivariant operations, high order tensors can be coupled very simply and thus the model can make direct predictions of high-order tensors such as dipole moments and polarizabilities without any modifications. The tests in several datasets show that HotPP not only achieves high accuracy in predicting target properties, but also successfully performs tasks such as calculating phonon spectra, infrared spectra, and Raman spectra, demonstrating its potential as a tool for future research.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Scattering of phonons by edge dislocation and thermal conductivity of nanocrystalline silicon

For this work, the temperature dependent thermal conductivity of nanocrystalline silicon with the grain size of 500 nm, 144 nm, and 76 nm, was calculated based on the intrinsic phonon properties using a first-principles-based method and recent analytical theory of phonon-dislocation scattering. The results demonstrate, in the nanocrystalline silicon, that the low temperature (T) limit of thermal conductivity is directly proportional to T 2 instead of the conventional Casimir limit of T 3 . Calculated frequency-dependent scattering rates also indicate regimes of dominant scattering mechanisms: frequency-independent boundary scattering dominates in the very low phonon frequency regime of less than 0.6 THz; frequency-dependent boundary scattering dominates in the intermediate phonon frequency regime between 0.6 and 6 THz; and Umklapp scattering dominates in the high phonon frequency regime beyond 6 THz. The thermal conductivity is linear proportional to nanograin size in the region of d < 100nm and it shows independence at d > 1μm. This work proves the importance of dislocation-phonon scattering to thermal conductivity when the size of nanograins is small or at low temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Doubling the superconducting transition temperature of ultraclean wafer-scale aluminum nanofilms

Here, we studied the role of reduced dimensionality and disorder in the superconducting properties of wafer-scale aluminum (Al) nanofilms. This new generation of ultrathin films were grown using molecular beam epitaxy and depict normal-state sheet resistance at least 20 times lower than the quantum resistance ℎ⁡/(4⁢𝑒 2 ). Defying general expectations, the superconducting transition temperature of our films increases with decreasing Al film thickness, reaching 2.4 K for a 3.5-nm-thick Al film grown on GaAs: twice that of bulk Al (1.2 K). Surface phonon softening is shown to impact superconductivity in pure ultrathin films, offering a route for materials engineering in two dimensions.

2-dimensional systems↗
DOE OSTI2025

Thermally induced structural competitiveness and metastability of body-centered-cubic iron under nonequilibrium conditions

The structure and stability of iron near melting at multi-megabar pressures are of significant interest in high pressure physics and earth and planetary sciences. While the body-centered cubic (BCC) phase is generally recognized as unstable at lower temperatures, its stability relative to the hexagonal close-packed (HCP) phase at high temperatures (approximately 0.5 eV) in the Earth’s inner core (IC) remains a topic of ongoing theoretical and experimental debate. Motivated by a series of studies attempting to understand the state of the IC and resulting in contradictory findings, we aim to explore the physical factors that may contribute to this confusion in a broader context. Our ab initio calculations show a significant drop in energy, the emergence of a plateau and a local minimum in the potential energy surface, and stabilization of all phonon modes of BCC iron at elevated electron temperatures (>1–1.5 eV). These effects increase the competition among the BCC, HCP, and the face-centered cubic (FCC) phases and lead to the metastability of the BCC structure. Furthermore, at pressure-temperature conditions relevant to the IC, we can expect substantial coexistence of FCC and HCP iron due to their similar free energies and low transition barrier (∼0.2 eV/atom), and the thermodynamic stability of BCC iron is enhanced by its substantial lattice vibration entropy. In conclusion, these findings provide a clear theoretical framework for understanding iron phase relations and solidification processes, both experimentally and in the IC.

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↗
NASA NTRSDate not supplied

Machine Learning Approaches for Rare-Earth Silicate Environmental Barrier Coating Thermochemical and Thermomechanical Property Predictions

Environmental barrier coatings (EBCs) are a necessary enabling technology for the transition from superalloys to silicon carbide (SiC) ceramic matrix composites (CMCs) in gas turbine engines for increased efficiency and decreased fuel costs. SiC-based CMCs are prone to oxidation-based degradation in the engine hot section, and rare-earth (RE) silicates are promising candidates for EBCs due to their close thermal expansion match to the composite substrate and oxidation resistance. However, the design of EBCs is hindered by the large chemical space of candidate materials and the difficulty in obtaining material properties for engineering optimization. This is especially difficult as research continues into mixed-cation or “high-entropy” RE silicates. First-principles computational methods such as density functional theory (DFT) are highly effective at calculating material properties to guide coating design but are limited by their computational cost. Atomistic simulations have the potential to both accelerate property calculations and expand the properties able to be calculated due to their lower computational compared to DFT. However, they require interatomic potentials (IAPs) specific to the material system of interest, and, to our knowledge, there are no suitable IAPs for RE silicates. Machine learning (ML) is a promising technique to accelerate material property predictions indirectly by generating IAPs for atomistic simulations or via direct prediction. In this work, we present two ML approaches to accelerate the calculation of RE silicate properties relevant to EBC design: 1) a ML-derived interatomic potential (IAP) for atomistic simulations of yttrium disilicate (Y2Si2O7) from DFT training data, and 2) a neural network (NN) model to directly predict thermochemical properties of RE silicates and oxides directly from easily obtainable unit cell parameters. Classical MD simulations using the IAP yield lattice properties and bond lengths in good agreement with both DFT and experimental results from x-ray diffraction. Thermodynamic properties calculated using the finite-displacement phonon method and quasi-harmonic approximation were orders of magnitude faster than DFT with good agreement to the DFT results. The IAP was also used to calculate properties such as coefficient of thermal expansion (CTE) that require large simulation supercells and are therefore difficult with DFT. The IAP correctly predicted the anisotropic nature of the CTE in three different phases of Y2Si2O7. The NN model predicts constant pressure heat capacity, Cp, orders of magnitude faster than DFT calculations, which can enable its use as a surrogate model for multiscale simulations. The two methods presented in this work demonstrate the utility of ML for accelerating the prediction of RE silicate properties, which can in turn accelerate EBC design and optimization.

machine learning↗
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

Anomalous thermal transport in Eshelby twisted van der Waals nanowires

Dislocations in van der Waals (vdW) layered nanomaterials induce strain and structural changes that substantially impact thermal transport. Understanding these effects could enable the manipulation of dislocations for improved thermoelectric and optoelectronic applications, but experimental insights remain limited. Here, in this study, we use synthetic Eshelby twisted vdW GeS nanowires (NWs) with single screw dislocations as a model system to explore the interplay between dislocation-induced structural modifications and lattice thermal conductivity. Our measurements reveal a monoclinic structure stabilized by the dislocation, leading to a substantial drop in thermal conductivity for larger-diameter NWs (70% at room temperature), supported by first-principles calculations. Interestingly, we also find an anomalous enhancement of thermal conductivity with decreasing diameter in twisted NWs, contrary to typical trends in non-twisted GeS NWs. This is attributed to increased conductivity near the NW cores due to compressive strain around the central dislocations, and aligns with a density-functional-theory-informed core–shell model. Our results highlight the critical role of dislocations in thermal conduction, providing fundamental insights for defect and strain engineering in advanced thermal applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Direct Observation of Enhanced Electron-Phonon Coupling in Copper Nanoparticles in the Warm-Dense Matter Regime

Warm-dense matter (WDM) is a highly-excited state that lies at the confluence of solids, plasmas, and liquids and that cannot be described by equilibrium theories. The transient nature of this state when created in a laboratory, as well as the difficulties in probing the strongly-coupled interactions between the electrons and the ions, make it challenging to develop a complete understanding of matter in this regime. Here in this work, by exciting isolated ~8 nm nanoparticles with a femtosecond laser below the ablation threshold, we create uniformly-excited WDM. We then use photoelectron spectroscopy to track the instantaneous electron temperature and directly extract the strongest electron-ion coupling observed experimentally to date. By directly comparing with state-of-the-art theories, we confirm that the superheated nanoparticles lie at the boundary between hot solids and plasmas, with associated strong electron-ion coupling. This is evidenced both by the fast energy loss of electrons to ions, as well as a strong modulation of the electron temperature by acoustic oscillations in the nanoparticle. This work demonstrates a new route for experimental exploration and theoretical validation of the exotic properties of WDM.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2024

First-Principles Calculations of the Electrical Conductivity of Carbon Nanotubes Functionalized with Copper and Nitrogen: Implications for Electronics, Energy Storage, and Nanodevices

In this work, we investigate the electrical conductivity of carbon nanotubes (CNTs), with a particular focus on the effects of doping. Using first-principles-based approaches, we study the electronic structure, phonon dispersion, and electron–phonon scattering to understand the finite-temperature electrical transport properties in CNTs. Our study covers both prototypical metallic and semiconducting CNTs, with special emphasis on the influence of typical defects such as vacancies and the incorporation of copper or nitrogen, such as pyridinic N, pyrrolic N, graphitic N, and oxidized N. Our theoretical study shows significant improvements in the electrical conduction properties of copper-CNT composites, especially when semiconducting CNTs are functionalized with nitrogen. Doping is found to cause significant changes in the electronic density of states near the Fermi level, which affects the electrical conductivity. Calculations show that certain types of functional groups, such as N-pyrrolic, result in more than 30-fold increase in the conductivity of semiconducting CNTs compared to Cu-incorporated CNTs alone. For metallic CNTs, the conductivity is in agreement with existing experimental data, and our prediction of significant increases in conductivity with N-pyrrolic functional group is consistent with recent experimental results, demonstrating the effectiveness of doping in modifying conductivity. In conclusion, our study provides valuable insight into the electronic properties of doped CNTs and contributes to the development of ultrahigh conductivity CNT composites.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Dynamic mode decomposition of nonequilibrium electron-phonon dynamics: accelerating the first-principles real-time Boltzmann equation

Abstract Nonequilibrium dynamics governed by electron–phonon ( e -ph) interactions plays a key role in electronic devices and spectroscopies and is central to understanding electronic excitations in materials. The real-time Boltzmann transport equation (rt-BTE) with collision processes computed from first principles can describe the coupled dynamics of electrons and atomic vibrations (phonons). Yet, a bottleneck of these simulations is the calculation of e –ph scattering integrals on dense momentum grids at each time step. Here we show a data-driven approach based on dynamic mode decomposition (DMD) that can accelerate the time propagation of the rt-BTE and identify dominant electronic processes. We apply this approach to two case studies, high-field charge transport and ultrafast excited electron relaxation. In both cases, simulating only a short time window of ~10% of the dynamics suffices to predict the dynamics from initial excitation to steady state using DMD extrapolation. Analysis of the momentum-space modes extracted from DMD sheds light on the microscopic mechanisms governing electron relaxation to a steady state or equilibrium. The combination of accuracy and efficiency makes our DMD-based method a valuable tool for investigating ultrafast dynamics in a wide range of materials.

36 MATERIALS SCIENCE↗
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↗
DOE OSTI2024

High-Temperature Observation of Intralayer, Interlayer, and Rydberg Excitons in Bulk Van Der Waals Alloy Single Crystals

Transition metal dichalcogenides exhibit remarkable optical properties due to the diverse number of strongly bound excitons, which can be fine-tuned by alloying. Despite a flurry of research activity in characterizing these excitons, a comprehensive and profound understanding of their behavior with temperature is lacking. Here, we report the rich spectrum of excitonic features within bulk van der Waals alloy Mo0.5 W0.5 S2 and Mo0.5 W0.5 Se2 single crystals through temperature-dependent reflectance spectroscopy and first-principles calculations. We observed Rydberg excitons and interlayer excitons in both the single crystals. Notably, we provide the first experimental evidence of highly energetic A' and B' excitons in Mo0.5 W0.5 S2 at room temperature. The strong carrier-phonon scattering significantly broadens the A', B', and interlayer excitons at room temperature in bulk Mo0.5 W0.5 S2 single crystal compared to its selenide. Our findings, supported by density functional theory and Bethe-Salpeter equation calculations, signify the crucial role of carrier-phonon interactions. These results open pathways for next-generation optoelectronic devices and quantum technologies operating at high temperature.

excitons↗