Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite-temperature dynamics”

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.

53 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Suppression of Composition g -modes in Chemically Equilibrating Warm Neutron Stars

We investigate the impact of chemical equilibration and the resulting bulk viscosity on nonradial oscillation modes of warm neutron stars at temperatures up to T ≈ 5 MeV, relevant for protoneutron stars and neutron star postmerger remnants. In this regime, the relaxation rate of weak interactions becomes comparable to the characteristic frequencies of composition g-modes in the core, resulting in resonant damping. To capture this effect, we introduce the dynamical sound speed, a complex, frequency-dependent generalization of the adiabatic sound speed that encodes both the restoring force and the dissipative effects of bulk compression. Using realistic weak reaction rates and three representative equations of state, we compute the complex frequencies of composition g-modes with finite-temperature profiles. We find that bulk viscous damping becomes increasingly significant with temperature and can completely suppress composition g-modes. In contrast, the f-mode remains largely unaffected by bulk viscosity due to its nearly divergence-free character. Our results highlight the sensitivity of g-mode behavior to thermal structure, weak reaction rates, and the equation of state, and establish the dynamical sound speed as a valuable descriptor characterizing oscillation properties in dissipative neutron star matter.

Neutron stars↗
DOE OSTI2026

Thermal exchange-correlation functionals: Capturing quantum electron behavior in warm, dense plasmas

We summarize and give perspective upon recent progress in developing non-empirical constraint-based thermal (i.e., free energy) exchange-correlation (XC) density functionals essential for accurate description of the quantum behavior of electrons in warm, dense plasmas. After delineating the critical role of ground-state functionals for zero-temperature, time-dependent DFT, we outline the underpinnings of local density approximation, generalized gradient approximation (GGA), and meta-GGA XC free-energy functionals. Two basic thermalization principles for upgrading ground-state XC functionals to successful thermal ones are emphasized. Then, we turn to a long-standing challenge, assessment of the accuracy of well-founded functionals. Unlike the ground state, there are a few exact results for large T and P. An exception is path integral Monte Carlo (PIMC) data for dense H/D and He plasmas. For those, we did ab initio molecular dynamics simulations under selected thermodynamic conditions employing five thermal XC functionals: two approximate thermal GGAs, fully thermal GGA, an approximate meta-GGA, and fully thermal meta-GGA. Comparisons with the PIMC data show that functionals thermalized by augmenting a non-thermal functional with a lower-level thermal contribution are inferior to functionals with thermal XC and spatial inhomogeneity effects taken into account at the same level of refinement. We believe this and similar evidence should be convincing to the high-energy density physics community of the necessity of use of proper thermal XC functionals in simulation studies of finite-temperature quantum effects in warm, dense plasmas.

Ab-initio molecular dynamics↗
DOE OSTI2022

Finite-temperature electron-capture rates for neutron-rich nuclei near N = 50 and effects on core-collapse supernova simulations

In this paper the temperature dependence of stellar electron-capture (EC) rates is investigated, with a focus on nuclei near N = 50, just above Z = 28, which play an important role during the collapse phase of core-collapse supernovae (CCSN). Two new microscopic calculations of stellar EC rates are obtained from relativistic and nonrelativistic finite-temperature quasiparticle random-phase approximation approaches, for a conventional grid of temperatures and densities. In both approaches, EC rates due to Gamow-Teller transitions are included. In the relativistic calculation, contributions from first-forbidden transitions are also included and add strongly to the EC rates. The new EC rates are compared with large-scale shell-model calculations for the specific case of 86 Kr, providing insight into the finite-temperature effects on the EC rates. At relevant thermodynamic conditions for core collapse, the discrepancies between the different calculations of this paper are within about one order of magnitude. Numerical simulations of CCSN are performed with the spherically symmetric gr1d simulation code to quantify the impact of such differences on the dynamics of the collapse. These simulations also include EC rates based on two parametrized approximations. A comparison of the neutrino luminosities and enclosed mass at core bounce shows that differences between simulations with different sets of EC rates are relatively small (≈ 5 %), suggesting that the EC rates used as inputs for these simulations have become well constrained.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

I ntera C hem : Exploring Excited States in Virtual Reality with Ab Initio Interactive Molecular Dynamics

InteraChem is an ab initio interactive molecular dynamics (AI-IMD) visualizer that leverages recent advances in virtual reality hardware and software, as well as the graphical processing unit (GPU)-accelerated TeraChem electronic structure package, in order to render quantum chemistry in real time. We introduce the exploration of electronically excited states via AI-IMD using the floating occupation molecular orbital-complete active space configuration interaction method. The optimization tools in InteraChem enable identification of excited state minima as well as minimum energy conical intersections for further characterization of excited state chemistry in small- to medium-sized systems. We demonstrate that finite-temperature Hartree–Fock theory is an efficient method to perform ground state AI-IMD. InteraChem allows users to track electronic properties such as molecular orbitals and bond order in real time, resulting in an interactive visualization tool that aids in the interpretation of excited state chemistry data and makes quantum chemistry more accessible for both research and educational purposes.

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

Adaptive variational quantum minimally entangled typical thermal states for finite temperature simulations

Scalable quantum algorithms for the simulation of quantum many-body systems in thermal equilibrium are important for predicting properties of quantum matter at finite temperatures. Here we describe and benchmark a quantum computing version of the minimally entangled typical thermal states (METTS) algorithm for which we adopt an adaptive variational approach to perform the required quantum imaginary time evolution. The algorithm, which we name AVQMETTS, dynamically generates compact and problem-specific quantum circuits, which are suitable for noisy intermediate-scale quantum (NISQ) hardware. We benchmark AVQMETTS on statevector simulators and perform thermal energy calculations of integrable and nonintegrable quantum spin models in one and two dimensions and demonstrate an approximately linear system-size scaling of the circuit complexity. We further map out the finite-temperature phase transition line of the two-dimensional transverse field Ising model. Finally, we study the impact of noise on AVQMETTS calculations using a phenomenological noise model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Electronic Structure of the Au Nanoparticle-TiO 2 Heterojunction: Influence of Nanoparticle Size, Shape, Oxygen Vacancies, and Temperature

The electronic structure of the gold nanoparticle-titanium dioxide (AuNP-TiO 2 ) heterojunction plays a critical role for charge transfer and recombination dynamics that underpin its photocatalytic function. However, building a representative model to capture the key physics remains a significant challenge. Here, we investigate the influence of the AuNP size and shape, as well as oxygen vacancy (V O ) defects at the anatase-phase TiO 2 (101) surface and the temperature of the heterojunction, on its interfacial electronic properties. Using density functional theory (DFT), we compare the closed-shell Au 20 and open-shell Au 19 clusters interfaced with pristine and V O defect TiO 2 surfaces. We find that the presence of a V O defect transforms pure TiO 2 from a p-type to an n-type semiconductor, reversing the interfacial band bending from downward to upward. For the heterosystem, density of states (DOS) analysis shows that V O minimally affects the open-shell Au 19 -TiO 2 system, but it significantly alters the closed-shell Au 20 -TiO 2 system, converting it from p-type to n-type at the Γ-point at 0 K. Furthermore, ab initio molecular dynamics (AIMD) simulations at 300 K reveal significant thermal fluctuations in AuNP positions relative to the TiO 2 surface. These fluctuations result in dynamic variations in the gap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) in all systems studied: Au 19 NP-Pristine TiO 2 , Au 20 NP-Pristine TiO 2 , Au 19 NP-V O defect TiO 2 , and Au 20 NP-V O defect TiO 2 . Complementary Bader charge analysis performed at both 0 K and finite-temperature AIMD snapshots supports the emergence of an upward band bending and the formation of a Schottky barrier at the V O containing heterojunctions. Notably, we find that at finite temperature, an Au atom can dynamically passivate the V O , leading to a pronounced widening of the HOMO-LUMO gap in the AuNP-V O defect TiO 2 heterojunctions. Our computational findings underscore the pivotal role of the V O defect and thermal effects in modulating interfacial band alignment, electronic states, and HOMO-LUMO gaps, providing insights for designing Au-TiO 2 heterojunctions with tailored electronic properties.

Mora Perez, Carlos [Lawrence Berkeley National Lab↗
DOE OSTI2020

Finite spatial-grid effects in energy-conserving particle-in-cell algorithms

Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. We also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length. Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. Finally, we also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Signatures of deconfined quark phases in binary neutron star mergers

We investigate the quark deconfinement phase transition in the context of binary neutron star (BNS) mergers. We employ a new finite-temperature composition-dependent equation of state (EOS) with a first order phase transition between hadrons and deconfined quarks to perform numerical relativity simulations of BNS mergers. Here, the softening of the EOS due to the phase transition causes the merger remnants to be more compact and to collapse to a black hole (BH) at earlier times. The phase transition is imprinted on the postmerger gravitational wave (GW) signal duration, amplitude, and peak frequency. However, this imprint is only detectable for binaries with sufficiently long-lived remnants. Moreover, the phase transition does not result in significant deviations from quasi-universal relations for the postmerger GW peak frequency. We also study the impact of the phase transition on dynamical ejecta, remnant accretion disk masses, r-process nucleosynthetic yields and associated electromagnetic (EM) counterparts. While there are differences in the EM counterparts and nucleosynthesis yields between the purely hadronic models and the models with phase transitions, these can be primarily ascribed to the difference in remnant collapse time between the two. An exception is the non-thermal afterglow caused by the interaction of the fastest component of the dynamical ejecta and the interstellar medium, which is systematically boosted in the binaries with phase transition as a consequence of the more violent merger they experience.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Thermodynamic stability of a spin microemulsion in Rashba spin-orbit-coupled bosons

Recent finite-temperature numerical simulations have unveiled a quantum “spin” microemulsion analog, found by raising the temperature of a stripe supersolid phase in a Rashba spin-orbit coupled Bose gas. This microemulsion state is a highly correlated, isotropic normal fluid where atoms self-arrange based on their internal pseudospin into patterns that resemble bicontinuous microemulsions. This finding leaves several open questions regarding the broader accessibility of this phase in experiments. Here, we use equilibrium finite-temperature numerical simulations based on a coherent-state path integral representation to perform a computational investigation into the thermodynamic stability of the spin microemulsion state. Numerical simulations emphasize the requirement of a nearly, but not perfectly, isotropic spin-orbit coupling in order to achieve the microemulsion phase in cold-atom experiments. Moreover, the microemulsion state exists independent of miscibility of the pseudospin components and for a wide range of pseudospin population imbalance, suggesting a high degree of flexibility in choosing the atom and hyperfine states in an experimental realization. Lastly, we demonstrate this feasibility by mimicking a Rashba spin-orbit-coupled 87 Rb experiment in an isotropic harmonic trap, where we confirm the microemulsion's existence via its density profile and equilibrium quasimomentum distribution.

Complex Langevin dynamics↗
DOE OSTI2020

Identification of non-Fermi liquid fermionic self-energy from quantum Monte Carlo data

Quantum Monte Carlo (QMC) simulations of correlated electron systems provide unbiased information about system behavior at a quantum critical point (QCP) and can verify or disprove the existing theories of non-Fermi liquid (NFL) behavior at a QCP. However, simulations are carried out at a finite temperature, where quantum critical features are masked by finite-temperature effects. Here, we present a theoretical framework within which it is possible to separate thermal and quantum effects and extract the information about NFL physics at T = 0. We demonstrate our method for a specific example of 2D fermions near an Ising ferromagnetic QCP. We show that one can extract from QMC data the zero-temperature form of fermionic self-energy Σ(ω) even though the leading contribution to the self-energy comes from thermal effects. We find that the frequency dependence of Σ(ω) agrees well with the analytic form obtained within the Eliashberg theory of dynamical quantum criticality, and obeys ω 2/3 scaling at low frequencies. Our results open up an avenue for QMC studies of quantum critical metals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Theory of finite-temperature two-dimensional structural transformations in group-IV monochalcogenide monolayers

An account of two-dimensional (2D) structural transformations in 2D ferroelectrics predicts an evolution from a structure with Pnm 2 1 symmetry into a structure with square P4/nmm symmetry and is consistent with experimental evidence, while another argues for a transformation into a structure with rectangular Pnmm symmetry. An analysis of the assumptions made in these models is provided here, and six fundamental results concerning these transformations are contributed as follows: (i) Softened phonon modes produce rotational modes in these materials. (ii) The transformation to a structure with P4/nmm symmetry occurs at the lowest critical temperature T c . (iii) The hypothesis that one unidirectional optical vibrational mode underpins the 2D transformation is unwarranted. (iv) Being successively more constrained, a succession of critical temperatures ( T c < T c ' < T c " ) occurs in going from molecular dynamics calculations with the NPT and NVT ensembles onto models with unidirectional oscillations. (v) The choice of exchange-correlation functional impacts the estimate of the critical temperature. (vi) Importantly, the correct physical picture of these transformations is one in which rotational modes confer a topological character to the 2D transformation via the proliferation of vortices.

2d ferroelectrics↗
DOE OSTI2021

Anomalous thermoelectricity at the two-dimensional structural transition of SnSe monolayers

The thermoelectric figure of merit ZT comprises electronic and vibrational contributions that drastically change across phase transitions, and the most common theoretical ab initio approach to thermoelectricity fails to describe the evolution of ZT across finite-temperature structural transitions in its entirety. Furthermore, while the thermoelectric behavior of bulk SnSe has been extensively studied, SnSe monolayers have been experimentally realized only recently, and the existent prediction of thermoelectricity on this two-dimensional material is unreliable because it misses its structural transition altogether. SnSe monolayers (and similar GeS, GeSe, GeTe, SnS, and SnTe monolayers) experience a temperature-induced two-dimensional Pnm2 1 → P4/nmm structural transition precipitated by the softening of vibrational modes, and we describe their thermoelectric properties across the phase transition, using molecular dynamics data to inform both electronic and vibrational coefficients directly and within the same footing. Similar to recent experimental observations pointing to an overestimated ZT past the transition temperature in bulk SnSe, we find a smaller ZT on SnSe monolayers when compared to its value predicted by the standard paradigm, due to the dramatic changes in the electrical conductivity and lattice thermal conductivity as the structural transition ensues. Here, the process described here lends a strong focus to both the vibrational and electronic evolutions throughout the structural transition, and it applies to thermoelectric materials undergoing thermally driven solid-to-solid structural phase transitions in one, two, and three dimensions.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Spectral-partitioned Kohn-Sham density functional theory

Here we introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Finite-Time Teleportation Phase Transition in Random Quantum Circuits

How long does it take to entangle two distant qubits in a quantum circuit evolved by generic unitary dynamics? Here, we show that if the time evolution is followed by measurements of all but two infinitely separated test qubits, then the entanglement between them can undergo a phase transition and become nonzero at a finite critical time t c . The fidelity of teleporting a quantum state from an input qubit to an infinitely distant output qubit shows the same critical onset. Specifically, these finite-time transitions occur in short-range interacting two-dimensional random unitary circuits and in sufficiently long-range interacting one-dimensional circuits. The phase transition is understood by mapping the random continuous-time evolution to a finite-temperature thermal state of an effective spin Hamiltonian, where the inverse temperature equals the evolution time in the circuit. In this framework, the entanglement between two distant qubits at times t > t c corresponds to the emergence of long-range ferromagnetic spin correlations below the critical temperature. We verify these predictions using numerical simulation of Clifford circuits and propose potential realizations in existing platforms for quantum simulation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Phase diagram of the three-dimensional subsystem toric code

Subsystem quantum error-correcting codes typically involve measuring a sequence of noncommuting parity check operators. They can sometimes exhibit greater fault tolerance than conventional codes, which use commuting checks. However, unlike subspace codes, it is unclear if subsystem codes—in particular their advantages—can be understood in terms of ground-state properties of a physical Hamiltonian. In this paper, we address this question for the three-dimensional subsystem toric code (3D STC), as recently constructed by Kubica and Vasmer [], which exhibits single-shot error correction. Motivated by a conjectured relation between single-shot properties and thermal stability, we study the zero- and finite-temperature phases of an associated noncommuting Hamiltonian. By mapping the Hamiltonian model to a pair of 3D Z 2 gauge theories coupled by a kinetic constraint, we find various phases at zero temperature, all separated by first-order transitions: There are 3D toric code-like phases with deconfined point-like excitations in the bulk, and there are phases with a confined bulk supporting a 2D toric code on the surface when appropriate boundary conditions are chosen. The latter is similar to the surface topological order present in 3D STC. However, the similarities between the single-shot correction in 3D STC and the confined phases are only partial: they share the same sets of degrees of freedom, but they are governed by different dynamical rules. Instead, we argue that the process of single-shot error correction can more suitably be associated with a path (rather than a point) in the zero-temperature phase diagram, a perspective, which inspires alternative measurement sequences enabling single-shot error correction. Moreover, since none of the above-mentioned phases survives at nonzero temperature, the single-shot error-correction property of the code does not imply thermal stability of the associated Hamiltonian phase. Published by the American Physical Society 2024

Li, Yaodong (ORCID:0000000337421944)↗
DOE OSTI2023

Recent progress on in-medium properties of heavy mesons from finite-temperature EFTs

Mesons with heavy flavor content are an exceptional probe of the hot QCD medium produced in heavy-ion collisions. In the past few years, significant progress has been made toward describing the modification of the properties of heavy mesons in the hadronic phase at finite temperature. Ground-state and excited-state thermal spectral properties can be computed within a self-consistent many-body approach that employs appropriate hadron-hadron effective interactions, providing a unique opportunity to confront hadronic Effective Field Theory predictions with recent and forthcoming lattice QCD simulations and experimental data. In this article, we revisit the application of the imaginary-time formalism to extend the calculation of unitarized scattering amplitudes from the vacuum to finite temperature. These methods allow us to obtain the ground-state thermal spectral functions. The thermal properties of the excited states that are dynamically generated within the molecular picture are also directly accessible. We present here the results of this approach for the open-charm and open-bottom sectors. We also analyze how the heavy-flavor transport properties, which are strongly correlated to experimental observables in heavy-ion collisions, are modified in hot matter. In particular, transport coefficients can be computed using an off-shell kinetic theory that is fully consistent with the effective theory describing the scattering processes. The results of this procedure for both charm and bottom transport coefficients are briefly discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1989

Conserving approximations for strongly correlated electron systems - Bethe-Salpeter equation and dynamics for the two-dimensional Hubbard model

A semianalytical approach is described for strongly correlated electronic systems which satisfies microscopic conservation laws, treats strong frequency and momentum dependences, and provides information on both static and dynamic properties. This approach may be used to treat large systems and temperatures lower than those currently accessible to finite-temperature quantum Monte Carlo techniques. Examples of such systems include heavy-electron compounds, organic Bechegaard salts, bis-(ethylenedithiolo)-TTF superconductors, and the oxide superconductors. The technique is based on the derivation and self-consistent solution of infinite-order conserving approximations. The technique is used to derive a low-temperature phase diagram and dynamic correlation functions for the two-dimensional Hubbard lattice model.

Bickers, N. E.↗