Search NASA⌕ Search

SEARCH · Search NASA

Results for “Phonon Simulation”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Integration of Silica in G4CMP for Phonon Simulations: Framework and Tools for Material Integration

Superconducting detectors with sub-eV energy resolution have demonstrated success setting limits on Beyond the Standard Model (BSM) physics due to their unique sensitivity to low-energy events. G4CMP, a Geant4-based extension for condensed matter physics, provides a comprehensive toolkit for modeling phonon and charge dynamics in cryogenic materials. This paper introduces a technical formalism to support the superconducting qubit and low-threshold detector community in implementing phonon simulations in custom materials into the G4CMP. As a case study, we present the results of a detailed analysis of silica phonon transport properties relevant for simulating substrate backgrounds in Beryllium Electron capture in Superconducting Tunnel junctions (BeEST)-style experiments using G4CMP. Additionally, Python-based tools were developed to aid users in implementing their own materials and are available on the G4CMP repository.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Atomic dynamics in 𝑀⁢Cr⁢𝑋 2 (𝑀=Ag, Cu; 𝑋 = S, Se) across magnetic and superionic transitions

Here, a systematic study of atomic dynamics and thermal properties of the family of layered chalcogenide compounds 𝑀⁢Cr⁢𝑋 2 (𝑀= Ag, Cu; 𝑋 = S, Se) was performed, including neutron and x-ray scattering, thermal characterization, and first-principles simulations. In all compounds, we observe a breakdown of specific phonon modes across the superionic phase transition, for phonons whose eigenvectors exhibit large contributions of mobile ions. In particular, the nondispersive portions of transverse acoustic (TA) branches at short-wavelengths and the low-energy optical phonons with large contributions from Ag + or Cu + become severely damped in the superionic phase. However, well-defined quasiparticles persist in the superionic state for long-wavelength TA phonons. In the case of AgCrS 2 , the coupling of lattice dynamics with its antiferromagnetic transition was also investigated. The magnetic ordering couples with the monoclinic–rhombohedral structural transition, and the Cr 3+ spin arrangement strongly affects the phonon dispersions. We qualitatively reproduce the magnetic and nuclear components of the INS measurement for antiferromagnetic AgCrS 2 by combining models of spin-waves and spin-polarized first-principles phonon simulations. Quasielastic magnetic fluctuations persist in the paramagnetic phase up to high temperature, but are clearly distinguished from the nuclear component through their momentum dependence. Finally, we report measurements of the thermal properties of the selenide compounds and find good agreement with our DFT simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

An Investigation of Thermal Properties of 2D Materials [Dissertation]

Studying the thermal conductivity of 2D materials is important due to the applications of 2D materials in fields such as thermal management, thermoelectricity, renewable energy, and sensors. As such, measurements of the thermal conductivity of these 2D materials become important to measure. Thermal conductivity is often difficult to measure for 2D materials due to their atomically thin nature and many experimental methods for doing so requiring contact with the sample, which can alter the thermal properties. A non-contact method for calculating the thermal conductivity of 2D materials supported on substrates in order to model the thermal conductivity of 2D materials for devices, is proposed and experimentally performed in this dissertation. The optothermal Raman technique is a useful non-contact diagnostic technique useful in determining the thermal conductivity of 2D materials. The optothermal Raman typically does not account for heat losses due to convection or radiation or substrate resistance, which are shown to be important factors to consider when developing an optothermal Raman model. Additionally, the calculation of the interfacial thermal conductance between the bottom surface of the sample and the top surface of the substrate, plays an important role in determining the final value of the thermal conductivity of a supported sample, and will yield differing results based on whether or not the conductance is calculated using an approach such as the Diffuse Mismatch Model (DMM) or calculated directly by varying the laser heating profile (usually done by changing the laser objective). This is shown to be the case for both graphene on Ni, graphene on Cu, and SnSe 2 on Cu. In addition to experimentally calculating the thermal conductivity of a 2D material with the optothermal Raman technique, the thermal conductivity of 2D materials can also be calculated using computational methods. The three-phonon method is a method which can be used to simulate phonon scattering processes and determine the thermal conductivity of semiconductors, wherein phonon scattering is the dominant mechanism which determines the thermal conductivity. The three-phonon method uses relaxation times for phonon scattering with other phonons, electrons, and other material system elements, such as isotopes or material defects, in order to create a single-mode relaxation time approximation (SMRTA), which is used to calculate the final value of the thermal conductivity. An important consideration when determining the thermal conductivity of a 2D material using this method is the device geometry, which is reflected in this work as the phonon-boundary scattering relaxation time. This inclusion is important along with the inclusion of phonon-electron scattering in accurately determining the thermal conductivity of a 2D material. In both the optothermal Raman experiments and the three-phonon method computations, strain is shown to have a demonstrable effect on the thermal conductivity of 2D materials. When a 1.1% strain was applied to the mechanical properties of SnSe, the three-phonon processes yielded a lower thermal conductivity than the no-strain case. For the optothermal Raman experiments, the strain induced in the Cu substrate and transferred to a single-layer graphene (SLG) sample yields a trend where the thermal conductivity of the SLG decreases with respect to strain applied. In the case where the interfacial thermal conductance was calculated directly, the conductance increased with respect to strain applied. This presents strain as a reliable and viable method for tuning the thermal properties of 2D materials for device applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Anharmonic phonons, thermal expansion, and nuclear quantum effects in Zn

The anharmonic behavior of phonons and thermal expansion of hexagonal zinc were studied from 15 to 690 K by inelastic neutron scattering (INS) and ab initio simulations. Phonon spectra were measured for Q-points covering the full Brillouin zone, giving the phonon density of states (DOS) and dispersions along high-symmetry directions. The dispersions were sharp at 15 K, but diffuse intensity was observed at energies above them. The dispersions broadened with temperature T and the diffuse intensity became stronger. This diffuse intensity appeared in all INS measurements and simulations, except for classical molecular dynamics at 15 K. The temperature-dependent effective potential (TDEP) method, which included the nuclear quantum effect from zero-point vibrational dynamics, was used to calculate the free energy and thermal expansion. For T < 100 K nuclear quantum effects were important for obtaining the correct negative thermal expansion, and path integral molecular dynamics (PIMD) was particularly effective for obtaining the negative thermal expansion in the basal plane. A Heisenberg-Langevin model for interacting phonons coupled to a thermal bath was able to reproduce the shape and intensity of the diffuse spectral features.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Modeling athermal phonons in novel materials using the G4CMP simulation toolkit

Understanding phonon and charge propagation in superconducting devices plays an important role in both performing low-threshold dark matter searches and limiting correlated errors in superconducting qubits. The Geant4 Condensed Matter Physics (G4CMP) package, originally developed for the Cryogenic Dark Matter Search (CDMS) experiment, models charge and phonon transport within silicon and germanium detectors and has been validated by experimental measurements of phonon caustics, mean charge-carrier drift velocities, and heat pulse propagation times. Here, in this work, we present a concise framework for expanding the capabilities for phonon transport to a number of other novel substrate materials of interest to the dark matter and quantum computing communities, including sapphire (Al 2 O 3 ), gallium arsenide (GaAs), lithium fluoride (LiF), calcium tungstate (CaWO 4 ), and calcium fluoride (CaF 2 ). We demonstrate the use of this framework in generating phonon transport properties of these materials and compare these properties with experimentally-determined values where available.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Modulation of thermal conductivity of iron-doped ß-Ga2O3 by helium-ion irradiation

This study examines the impact of helium-ion irradiation on the thermal conductivity of ß-Ga2O3. A laser-based spatial domain thermoreflectance technique is used to investigate thermal conductivity map for both un-irradiated and irradiated ß-Ga2O3, which are then validated against simulation results derived from density functional theory-based phonon transport simulations. Since helium bubble evolution was ob- served at the nanoscale using transmission electron microscopy, the simulation study was carried out on eight distinct helium-induced sites in ß-Ga2O3. Our findings indicate a reduction in thermal conductivity for the irradiated samples. Experimental results show a significant reduction in thermal conductivity in irradiated samples, with de- creases of approximately 25% along the [100] direction and 40% along [001] directions. Phonon transport simulations closely replicate these findings, particularly when helium occupying interstitial sites, predicting reductions of ˜53% along [100] and ˜50% along [001] directions. This work underscores the role of irradiation-induced microstructural changes in the heat transport properties of ß-Ga2O3 which is crucial for its application in sensor devices in extreme environments.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Spontaneously formed phonon frequency combs in van der Waals solid CrGeTe 3 and CrSiTe 3

Optical phonon engineering through nonlinear effects has been utilized in ultrafast control of material properties. However, nonlinear optical phonons typically exhibit rapid decay due to strong mode-mode couplings, limiting their effectiveness in temperature or frequency sensitive applications. Here we report the observation of long-lived nonlinear optical phonons through the spontaneous formation of phonon frequency combs in the van der Waals material CrXTe 3 (X=Ge, Si) using high-resolution Raman scattering. Unlike conventional optical phonons, the highest A g mode in CrGeTe 3 splits into equidistant, sharp peaks forming a frequency comb that persists for hundreds of oscillations and survives up to 200K. These modes correspond to localized oscillations of Ge 2 Te 6 clusters, isolated from Cr hexagons, behaving as independent quantum oscillators. Introducing a cubic nonlinear term to the harmonic oscillator model, we simulate the phonon time evolution and successfully replicate the observed comb structure. Similar frequency comb behavior is observed in CrSiTe 3 , demonstrating the generalizability of this phenomenon. Our findings demonstrate that Raman scattering effectively probes high-frequency nonlinear phonon modes, offering insight into the generation of long-lived, tunable phonon frequency combs with potential applications in ultrafast material control and phonon-based technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Advancing simulations of coupled electron and phonon nonequilibrium dynamics using adaptive and multirate time integration

Electronic structure calculations in the time domain provide a deeper understanding of nonequilibrium dynamics in materials. The real-time Boltzmann equation (rt-BTE), used in conjunction with accurate interactions computed from first principles, has enabled reliable predictions of coupled electron and lattice dynamics. However, the timescales and system sizes accessible with this approach are still limited, with two main challenges being the different timescales of electron and phonon interactions and the cost of computing collision integrals. As a result, only a few examples of these calculations exist, mainly for two-dimensional (2D) materials. Here we leverage adaptive and multirate time integration methods to achieve a major step forward in solving the coupled rt-BTEs for electrons and phonons. Relative to conventional (non-adaptive) time-stepping, our approach achieves a 10x speedup for a target accuracy, or greater accuracy by 3–6 orders of magnitude for the same computational cost, enabling efficient calculations in both 2D and bulk materials. This efficiency is showcased by computing the coupled electron and lattice dynamics in graphene up to ~100 ps, as well as modeling ultrafast lattice dynamics and thermal diffuse scattering maps in bulk materials (silicon and gallium arsenide). In addition to improved efficiency, our adaptive method can resolve the characteristic rates of different physical processes, thus naturally bridging different timescales. This enables simulations of longer timescales and provides a framework for modeling multiscale dynamics of coupled degrees of freedom in matter. Our work opens new opportunities for quantitative studies of nonequilibrium physics in materials, including driven lattice dynamics with phonons coupled to electrons, spin, and other degrees of freedom.

Yao, Jia [California Institute of Technology (CalT↗
DOE OSTI2025

Nanoscale Decoupling of Carrier-Phonon Transport in Carbon Nanotube-Halide Perovskite Heterostructures

In conventional semiconductors, electrical and thermal conductivity are typically coupled, posing a challenge in optimizing both simultaneously. Overcoming this inherent trade-off enables strategies for advancing electronic applications. Herein, a strategy is demonstrated to decouple electrical and thermal conductivity trade-off by creating heterostructures of highly conductive single-walled carbon nanotubes (SWCNTs) coated with low conductivity hybrid perovskites. Coating SWCNTs with methylammonium lead iodide perovskite results in an enhancement in electrical conductivity (408-1266 S cm-1) due to p-type doping followed by a threefold decrease of the in-plane thermal conductivity (3.3-1 W m-1 K-1), compared to pristine SWCNTs. Molecular dynamics simulations uncover phonon boundary scattering at the SWCNT/perovskite interface as well as localization of methylammonium-related and softening of the Pb-I-related phonon modes in methylammonium lead iodide perovskite decreasing the thermal conductivity.

14 SOLAR ENERGY↗
DOE OSTI2026

Complex Dynamics in Argyrodite Solid-State Ion Conductors

Argyrodites are a compositionally diverse family of materials that exhibit remarkable ion transport properties. While the average crystal structures of argyrodites have been extensively studied, ion transport in these materials is governed by a confluence of dynamic processes spanning the cation, anion, and polyanionic sublattices. This Perspective synthesizes recent advances in understanding the role of dynamics in structural behavior and ion transport properties. We examine the compositional and structural motifs that govern order−disorder transitions within the argyrodite family and further explore how ion hopping is facilitated by lattice dynamics, from long-range phonons to local rotational dynamics of polyanionic species. Through the lens of dynamics spanning multiple time and length scales, we establish guiding principles that govern transport phenomena and highlight avenues of future study for the argyrodite family of ion conductors.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Ultranarrow bright single-photon emitters in diamond with strong broadband phonon decoupling

Single-photon emitters are fundamental building blocks for quantum information processing, communication and sensing. However, unwanted interactions with bulk phonons in their host environment strongly limit their coherence and controllability. We report single color centers in nanodiamonds that are strongly and comprehensively decoupled from the bulk phononic environment. The color centers feature record-narrow linewidths down to 0.3 nm at room temperature and stable, bright emission, exceeding 10 Mcps in saturation. Notably, the bulk phonon sideband is almost entirely suppressed, revealing the presence of a single localized vibrational mode outside the diamond phonon band. Our observations and simulations point towards a unique mechanism for phonon decoupling in common wide-gap materials, based on a strongly radiative orbital transition coupled to a localized vibrational mode. The new color center enables qualitatively higher performance for applications in quantum networks and nanoscale sensing, and the exploration of new physical resources associated with vibrational states.

Sahoo, Swetapadma [ORNL]↗
DOE OSTI2024

Simulating Superconducting Devices in Novel Materials Using the G4CMP Simulation Toolkit

Understanding the kinematics of phonons and charge propagation in superconducting devices is crucial for conducting low-threshold dark matter searches and minimizing correlated errors in superconducting qubits. For nearly a decade, the Geant4 Condensed Matter Physics (G4CMP) package, originally developed for the Cryogenic Dark Matter Search (CDMS) experiment, has been limited to simulating charge and phonon transport in silicon and germanium materials. In this work, we have expanded the capabilities of G4CMP to include novel substrate materials such as sapphire (Al2O3), Gallium Arsenide (GaAs), Lithium Fluoride (LiF), Calcium Tungstate (CaWO4), and Calcium Fluoride (CaF2). We demonstrate the use of this toolkit in generating phonon transport properties for these materials and compare the results with experimentally determined values where available. Additionally, we explore the performance of superconducting devices utilizing these materials.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Phonon second harmonic generation in NaBr studied by inelastic neutron scattering and computer simulation

The phenomenon of second harmonic generation (SHG) was found for phonons in anharmonic NaBr by inelastic neutron scattering. The temperature dependence of this phonon SHG was measured from 300 K to 650 K. At 300 K the second harmonic (SH) is seen as a high-energy branch around 33 meV, nearly independent of $\overrightarrow{Q}$. The temperature effective potential (TDEP) method and classical molecular dynamics (MD) simulation with machine learning interatomic potential were able to reproduce the SH, and showed that SHG occurs with the flat transverse optical (TO) phonon branch. A classical model of a nonlinear medium explains the intensity and lifetime of the SH, compared to those of the TO modes. Also successful was a quantum model based on the Heisenberg-Langevin equation for interacting phonons coupled to a thermal bath, which also predicts a spectral distribution of the SH. In conclusion, the measured temperature dependence of the intensity of the second harmonic showed that it follows the Planck distribution of a one-phonon quasiparticle, and not two TO phonons.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Coherent phonon flatband generated in GaAs/AlAs superlattices via layer-selective optical pumping

Flatbands, characterized by their dispersionless energy levels in electronic, magnetic, and phononic systems, hold substantial potential for advancements in electronics and quantum information processing. Most flatbands exist in thermal equilibrium and cannot be easily created or annihilated externally, limiting their flexibility as switchable knobs for use in microelectronics and quantum applications. In our work, we demonstrate the generation of a coherent phonon flatband in a GaAs/AlAs superlattice using 800 nm femtosecond laser pulses. This coherent phonon flatband does not correspond to a phonon eigenmode at equilibrium and exhibits strong coupling with two branches of coherently excited longitudinal phonon modes. With molecular dynamics simulations, we show more generally that the coherent phonon flatband can be induced by coherently and spatially modulated optical excitations of superlattice structures. Our results highlight a pathway for coherent phonon flatband creation in the time domain that can be generalized to various superlattice systems, potentially inspiring the realization of coherent flatband generation of other quasiparticles.

Condensed-matter physics↗
DOE OSTI2025

Revealing the Phonon Bottleneck Limit in Negatively Charged CdS Quantum Dots

The capture of photoexcited hot electrons in semiconductors before they lose their excess energy to cooling is a long-standing goal in photon energy conversion. Semiconductor nanocrystals have large electron energy spacings that are expected to slow down electron relaxation by phonon emission, but hot electrons in photoexcited nanocrystals, nevertheless, cool rapidly by energy transfer to holes. This makes the intrinsic phonon-bottleneck-limited hot electron lifetime in nanocrystals elusive. For this work, we used a combination of theory and experiments to probe the hot-electron dynamics of negatively charged cadmium sulfide (CdS) colloidal quantum dots (QDs) in the absence of holes. Experiments found that these hot electrons cooled on a 100 ps time scale. Theoretical simulations predicted that pure phonon-bottleneck-limited hot electron cooling occurs on a similar time scale. This similarity suggests that the experimental measurements reflect the upper limit on the hot-electron lifetimes in these CdS QDs and the lower limit on the rates of processes that can harvest those hot electrons.

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

Synthetic Tuning of Exciton–Phonon Coupling in Janus WS 2(1-x) Se 2x Monolayers Revealed by Resonant Raman Excitation Spectroscopy for Optoelectronic Applications

Janus monolayers, such as WSSe, have broken out-of-plane symmetry and an intrinsic dipole moment, impacting exciton transport, lifetime, and phonon interactions while imbuing piezoelectric, photocatalytic, and Rashba spin-splitting properties to transition metal dichalcogenides (TMDs). The new properties of this atomically thin material can be used for optoelectronic device applications. As TMDs are converted into Janus monolayers, e.g., top selenization of WS2 to WSSe, the bandgap and structure smoothly evolve, impacting not only the formation of excitons but also their complex interactions with different phonon modes. Resonant Raman excitation profiles (REPs) are uniquely well-suited to reveal both excitonic transitions and exciton–phonon coupling. Here, the resonant REPs of $A^{'}_{1}$ WS 2 and A 1 WSSe modes are measured to understand the strength of their coupling with the A, B, and C excitonic bands of a WS2 monolayer throughout its stepwise transformation into Janus WSSe by pulsed laser deposition (PLD) of energetic selenium species. In situ Raman spectroscopy during deposition is used to controllably prepare stable intermediate Janus structures, WS 2(1-x) Se 2x (0 ≤ x ≤ 0.5), for ex situ measurement of their resonant REPs. As x increases, REPs reveal not only pronounced excitonic bands that gradually shift toward lower photon energies but also strong, mode-selective exciton–phonon coupling. First-principles resonant Raman simulations independently predict this spectral behavior and are shown capable of matching the spectrally broadened, experimentally observed REP profiles in this model system, indicating their strong predictive capability for future experiments. The combination of controlled synthesis, REP characterization, and predictive theory employed here demonstrates a powerful pathway to understand and ultimately tune exciton–phonon interactions for future quantum optical devices.

Janus monolayers↗