Machine-Learned, Finite Temperature Fermi-Operator Expansions Suitable for GPUs and AI-Hardware
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Herein we study the nonlinear evolution of very oblique small-scale Alfvén waves with k ⊥ d i ≳1. At these scales, the waves become significantly compressive, unlike in magnetohydrodynamics, due to the Hall term in the equations. We demonstrate that when frequencies are small compared to the ion gyrofrequency and amplitudes are small compared to unity, no new nonlinear interaction appears due to the Hall term alone at the lowest non-trivial order, even when k ⊥ d i ~1. However, at the second non-trivial order, we discover that the Hall physics leads to a slow but resonant nonlinear interaction between co-propagating Alfvén waves, an inherently three-dimensional effect. Including the effects of finite temperature, finite frequency, and electron inertia, the two-fluid Alfvén wave also becomes dispersive once one or more of k ⊥ ρ s , k ⊥ d e , or k ∥ d i becomes significant: for oblique waves at low β as studied here, this can be at a much smaller scale than d i . We show that the timescale for one-dimensional steepening of two-fluid Alfvén waves is only significant at these smaller dispersive scales, and also derive an expression for the amplitude of driven harmonics of a primary wave. Importantly, both new effects are absent in gyrokinetics and other commonly used reduced two-fluid models. Our calculations have relevance for the interpretation of laboratory Alfvén wave experiments, as well as shedding light on the physics of turbulence in the solar corona and inner solar wind, where the dominant nonlinear interaction between counter-propagating waves is suppressed, allowing these new effects to become important.
Abstract Symmetry-protected topological phases 1–4 cannot be described by any local order parameter and are beyond the conventional symmetry-breaking model 5 . They are characterized by topological boundary modes that remain stable under symmetry respecting perturbations 1–4,6–8 . In clean, gapped systems without disorder, the stability of these edge modes is restricted to the zero-temperature manifold; at finite temperatures, interactions with mobile thermal excitations lead to their decay 9–11 . Here we report the observation of a distinct type of topological edge mode 12–14 , which is protected by emergent symmetries and persists across the entire spectrum, in an array of 100 programmable superconducting qubits. Through digital quantum simulation of a one-dimensional disorder-free stabilizer Hamiltonian, we observe robust long-lived topological edge modes over up to 30 cycles for a wide range of initial states. We show that the interaction between these edge modes and bulk excitations can be suppressed by dimerizing the stabilizer strength, leading to an emergent U(1) × U(1) symmetry in the prethermal regime of the system. Furthermore, we exploit these topological edge modes as logical qubits and prepare a logical Bell state, which exhibits persistent coherence, despite the system being disorder-free and at finite temperature. Our results establish a viable digital simulation approach 15–18 to experimentally study topological matter at finite temperature and demonstrate a potential route to construct long-lived, robust boundary qubits in disorder-free systems.
Thermal density functional theory is commonly used in simulations of warm dense matter, a highly energetic phase characterized by substantial thermal effects and by correlated electrons demanding quantum mechanical treatment. Methods that account for temperature dependence, such as Mermin–Kohn–Sham finite-temperature density functional theory and free energy density functional theory, are now employed with more regularity and available in many standard code packages. However, approximations from zero-temperature density functional theory are still often used in temperature-dependent simulations using thermally weighted electronic densities as an input to exchange–correlation functional approximations, a practice known to miss temperature-dependent effects in the exchange–correlation free energy of these systems. Here, in this work, the temperature-dependent adiabatic connection is demonstrated and analyzed using a well-known parameterization of the uniform electron gas free energy. Useful tools based on this formalism for analyzing and constraining approximations of the exchange–correlation at zero temperature are leveraged for the finite-temperature case. Inspired by the Lieb–Oxford inequality, which provides a lower bound for the ground-state exchange–correlation energy, bounds for the exchange–correlation at finite temperatures are approximated for various degrees of electronic correlation.
We analyze, in perturbation theory, a theory of weakly interacting fractons and non-relativistic fermions in a 2+1 dimensional Quantum Field Theory. In particular we compute the 1-loop corrections to the self energies and interaction vertex, and calculate the associated 1-loop Renormalization Group flows of the coupling constants. Surprisingly, we find that the fracton-fermion coupling does not flow due to an emergent coordinate-dependent symmetry of the effective Lagrangian, making this model a well-defined quantum field theory. We provide additional discussions on the regularization and renormalization of interacting fractonic theories, as well as both qualitative and quantitative remarks regarding the theory at finite temperature and finite chemical potential.
The widely used thermal Hartree–Fock (HF) theory is generalized to include the effect of electron correlation while maintaining its quasi-independent-particle framework. An electron-correlated internal energy (or grand potential) is postulated in consultation with the second-order finite-temperature many-body perturbation theory (MBPT), which then dictates the corresponding thermal orbital (quasiparticle) energies in such a way that all fundamental thermodynamic relations are obeyed. The associated density matrix is of a one-electron type, whose diagonal elements take the form of the Fermi–Dirac distribution functions, when the grand potential is minimized. The formulas for the entropy and chemical potential are unchanged from those of Fermi–Dirac or thermal HF theory. The theory thus stipulates a finite-temperature extension of the second-order Dyson self-energy of one-particle many-body Green’s function theory and can be viewed as a second-order, diagonal, frequency-independent, thermal inverse Dyson equation. At low temperatures, the theory approaches finite-temperature MBPT of the same order, but it may outperform the latter at intermediate temperatures by including additional electron-correlation effects through orbital energies. Here, a physical meaning of these thermal orbital energies is proposed (encompassing that of thermal HF orbital energies, which has been elusive) as a finite-temperature version of Janak’s theorem.
In this review, we provide an up-to-date account of quantitative holographic descriptions of the strongly coupled quark-gluon plasma (QGP) produced in heavy-ion collisions, based on the class of gauge-gravity Einstein-Maxwell-Dilaton (EMD) models. Holography is employed to tentatively map the QCD phase diagram at finite temperature onto a dual theory of charged, asymptotically AdS black holes in 5D. With a quantitative focus on the hot QCD phase diagram, the EMD models reviewed are adjusted to describe lattice results for the finite-temperature QCD equation of state, with 2+1 flavors and physical quark masses, at zero chemical potential and vanishing electromagnetic fields. The predictive power of EMD models is tested by quantitatively comparing their predictions for the hot QCD equation of state at nonzero baryon density and the corresponding state-of-the-art lattice QCD results. The shear and bulk viscosities predicted by these EMD models are also compared to the corresponding profiles favored by the latest phenomenological multistage models describing different heavy-ion data. We report preliminary results from a Bayesian analysis which provide systematic evidence that lattice results at finite temperature and zero baryon density strongly constrains the free parameters of EMD models. Remarkably, the set of parameters constrained by lattice results at zero chemical potential produces EMD models in quantitative agreement with lattice QCD results also at finite baryon density. We also review results for equilibrium and transport properties from magnetic EMD models, describing the QGP at finite temperatures and magnetic fields. Lastly, we provide a critical assessment of the main limitations and drawbacks of the holographic models reviewed in the present work, and point out some perspectives we believe are of fundamental importance for future developments.
We propose a systematic coset construction of non-equilibrium effective field theories (EFTs) governing the long-distance and late-time dynamics of relativistic, finite-temperature condensed matter systems. Our non-equilibrium coset construction makes significant advances beyond more standard coset constructions in that it takes advantage of recently-developed techniques, which allow the formulation of non-equilibrium effective actions that account for quantum and thermal fluctuations as well as dissipation. Because these systems exist at finite temperature, the EFTs live on the closed-time-path of the Schwinger-Keldysh contour. Since the coset construction and the non-equilibrium effective actions may be unfamiliar to many readers, we include brief introductions to these topics in an effort to make this paper self-contained. To demonstrate the legitimacy of this coset construction, we successfully reproduce the known EFTs for fluids and superfluids at finite temperature. Then, to demonstrate its utility, we construct novel EFTs for solids, supersolids, and four phases of liquid crystals, all at finite temperature. We thereby combine the non-equilibrium effective action and the coset construction to create a powerful tool that can be used to study many-body systems out of thermal equilibrium.
Ternary ABX 3 perovskites made of corner-sharing BX 6 octahedra have long featured prominently in solid-state chemistry and condensed matter physics. Still, the joint understanding of their two main subgroups—halides and oxides—has not been fully developed. Indeed, unlike the case in simpler compounds having a single, robust repeated motif (“monomorphous”), certain cubic perovskites can manifest a non-thermal (= intrinsic) distribution of local motifs (“polymorphous networks”). Such static deformations can include positional degrees of freedom (e.g., atomic displacements and octahedral tilting) or magnetic moment degrees of freedom in paramagnets. Unlike thermal motion, such static distortions do not time-average to zero, being an expression of the intrinsic symmetry breaking preference of the chemical bonding. Here, the present study compares electronic structure features of oxide and halide perovskites starting from the static polymorphous distribution of motifs described by Density Functional Theory (DFT) minimization of the internal energy, continuing to finite temperature thermal disorder modeled via finite temperature DFT molecular dynamics. We find that (i) different oxide vs. halide ABX 3 compounds adopt different energy-lowering symmetry-breaking modes. The calculated pair distribution function (PDF) of SrTiO 3 from the first-principles agrees with recently measured PDF. (ii) In both oxides and halides, such static distortions lead to band gap blueshifts with respect to undistorted cubic Pm-3m structure. (iii) For oxide perovskites, high-temperature molecular dynamics simulations initiated from the statically distorted polymorphous structures reveal that the thermally-induced distortions can lead to a band gap redshift. (iv) In contrast, for cubic halide perovskite CsPbI 3 , both the intrinsic distortions and the thermal distortions contribute in tandem to band gap blueshift, the former, intrinsic effect being dominant. (v) In the oxide SrTiO 3 and CaTiO 3 (but not in halide) perovskites, octahedral tilting leads to the emergence of a distinct Γ–Γ direct band gap component as a secondary valley minimum to the well-known indirect R–Γ gap. Understanding such intrinsic vs. thermal effects on oxide vs. halide perovskites holds the potential for designing target electronic properties.
We calculate the magnetic dipole $\gamma$-ray strength functions in a chain of even-mass neodymium isotopes $^{144-152}$Nd in the framework of the configuration-interaction (CI) shell model. We infer the strength function by applying the maximum entropy method (MEM) to the exact imaginary-time response function calculated with the shell-model Monte Carlo (SMMC) method. The success of the MEM depends on the choice of a good strength function as a prior distribution. We investigate two choices for the prior strength function: the static path approximation (SPA) and the quasiparticle random-phase approximation (QRPA). We find that the QRPA is a better approximation at low temperatures (i.e., near the ground state), while the SPA is a better choice at finite temperatures. We identify a low-energy enhancement (LEE) in the MEM deexcitation $M1$ strength functions of the even-mass neodymium isotopes and compare with recent experimental results for the total deexcitation $\gamma$-ray strength functions. The LEE is already seen in the SPA strength function but not in the QRPA strength function, indicating the importance of large-amplitude static fluctuations around the mean field in reproducing the LEE. Our method is currently the only one which can reproduce LEE in heavy open-shell nuclei where conventional CI shell model calculations are prohibited. With the onset of deformation as number of neutrons increases along the chain of neodymium isotopes, we observe that some of the LEE strength transfers to a low-energy excitation, which we interpret as a finite-temperature ``scissors'' mode. Here, we also observe a finite-temperature spin-flip mode.
Moiré systems offer an exciting playground to study the many-body effects of strongly correlated electrons in regimes that are not easily accessible in conventional material settings. Motivated by a recent experiment on WSe 2 /WS 2 moiré bilayers [Y. Tang et al., Nature (London) 579, 353 (2020)], which realized a triangular superlattice with a small hopping (of approximately 10 Kelvin), with tunable density of holes, we explore the Hubbard model on the triangular lattice for intermediate temperatures t ≲ T < U. Employing finite-temperature Lanczos calculations and closely following the fitting protocols used in the experiment, we recover the observed trends in the reported Curie-Weiss temperature Θ with filling, using the reported interaction strength U/t=20. Here, we focus on the large increase of |Θ| on decreasing the density below half filling and the sign change of Θ at higher fillings, which signals the onset of ferromagnetism. The increase in |Θ| is also seen in the t–J model (the low-energy limit of the Hubbard model) in the intermediate temperature range, which we clarify is opposite to the trend in its high temperature limit. Differences between the low, intermediate, and high temperature regimes are discussed. Our numerical calculations also capture the crossover between short-range antiferromagnetic to ferromagnetic order in the intermediate temperature regime, a result broadly consistent with the experimental findings. We find that this behavior is a finite-temperature remnant of the underlying zero-temperature phase transition, which we explore with ground-state density matrix renormalization group calculations. We provide evidence of ferromagnetism characterized by weak (but robust) correlations that explain the small Θ seen in the experiment.
We report characterizing the temperature-dependent electronic structures of optoelectronic materials is vital for their application in high-temperature processes. Allen–Heine–Cardona (AHC) theory is an established theory to determine the impact of electron–phonon coupling on the temperature renormalization of the electronic structure. A question of fundamental interest is how applicable this approach can be in a broad range of materials. Herein, we establish the consensus on the temperature dependence of band gaps between AHC predictions and experimental measurements, in a versatile set of nine MO x metal oxides and six ABO 3 perovskites that are utilized as gas sensors. Both calculations and measurements point to the appreciable band gap downshifts within the temperature range relevant to typical high-temperature applications in all examined materials. The correlation between the band gaps and temperatures was quantified at zero- and finite-temperature limits using the analytical parameters of the O’Donnell model. We found fairly comparative agreements between AHC- and experiment-derived parameters, especially for the parameters (S and $\langle\hbar$ω$\rangle$) defining the temperature dependence behavior of band gaps at the finite-temperature regime. Our findings provide important implications for the first-principles modeling of the electronic structures and their temperature dependency, which should be the focus toward high-temperature optoelectronic applications.
In this paper, the effects of warm ions on electrostatic, kinetic instabilities observed in partially magnetized plasmas, such as the electron cyclotron drift instability (ECDI) and modified two-stream instability (MTSI), are investigated. The kinetic dispersion relation that accounts for warm, non-magnetized ions and warm, magnetized electrons is solved using a complex root-finding technique [A. C. Denig and K. Hara, Phys. Plasmas 30, 032108 (2023)]. It is observed that finite ion temperature introduces ion Landau damping effects, decreasing the growth rate of most unstable modes. The high-frequency, short-wavelength ECDI modes are more strongly affected by ion Landau damping than the low-frequency, long-wavelength MTSI modes. The solutions of the kinetic dispersion relation for MTSI are compared with the fluid dispersion relation of MTSI, which was originally proposed by McBride et al. [Phys. Fluids 15 2367–2383 (1972)]. The differences in the dispersion characteristics between the kinetic and fluid MTSI limits are assessed, illustrating the effects of electron and ion Landau damping on kinetic instabilities in partially magnetized plasmas.
While significant magnetic interactions exist in lithium transition metal oxides, commonly used as Li-ion cathodes, the interplay between magnetic couplings, disorder, and redox processes remains poorly understood. In this work, we focus on the high-voltage spinel LiNi 0.5 Mn 1.5 O 4 (LNMO) cathode as a model system on which to apply a computational framework that uses first principles-based statistical mechanics methods to predict the finite temperature magnetic properties of materials and provide insights into the complex interplay between magnetic and chemical degrees of freedom. Density functional theory calculations on multiple distinct Ni–Mn orderings within the LNMO system, including the ordered ground-state structure (space group P4332), reveal a preference for a ferrimagnetic arrangement of the Ni and Mn sublattices due to strong antiferromagnetic superexchange interactions between neighboring Mn 4+ and Ni 2+ ions and ferromagnetic Mn–Mn and Ni–Ni couplings, as revealed by magnetic cluster expansions. These results are consistent with qualitative predictions using the Goodenough-Kanamori-Anderson rules. Simulations of the finite temperature magnetic properties of LNMO are conducted using Metropolis Monte Carlo. We find that a “semiclassical” Monte Carlo sampling method based on the Heisenberg Hamiltonian accurately predicts experimental magnetic transition temperatures observed in magnetometry measurements. This study highlights the importance of a robust computational toolkit that accurately captures the complex chemomagnetic interactions and predicts finite temperature magnetic behavior to help analyze experimental magnetic and magnetic resonance spectroscopy data acquired ex situ and operando.
Effects of finite-temperature quasiparticle self-energy corrections on x-ray absorption spectra are investigated within the finite-temperature quasiparticle local density GW approximation up to temperatures T of order the Fermi temperature. To facilitate the calculations, we parametrize the quasiparticle self-energy using low-order polynomial fits. Here we show that temperature-driven decrease in the electron lifetime substantially broadens the spectra in the near-edge region with increasing T. However, the quasiparticle shift is most strongly modified near the onset of plasmon excitations.
Density functional theory (DFT) calculations, within the virtual crystal alloy approximation, are performed, along with the development of a Landau-type model employing a symmetry-allowed analytical expression of the internal energy and having parameters determined from first principles, to investigate properties and energetics of Al1-xScxN ferroelectric nitrides in their hexagonal forms. These DFT computations and this model predict the existence of two different types of minima, namely, the fourfold-coordinated wurtzite (WZ) polar structure and a five-fold coordinated paraelectric hexagonal phase (denoted as H5), for any Sc composition up to 40%. The H5 minimum progressively becomes the lowest-energy state within hexagonal symmetry as the Sc concentration increases from 0 to 0.4. Furthermore, the model points to several key findings. Examples include the crucial role of the coupling between polarization and strains to create the WZ minimum, in addition to polar and elastic energies, and that the origin of the H5 state overcoming the WZ phase as the global minimum within hexagonal symmetry when increasing the Sc composition mostly lies in the compositional dependency of only two parameters-one linked to the polarization and another one being purely elastic in nature. Other examples are that forcing Al1-xScxN systems to have no or a weak change in lattice parameters when heating them allows us to reproduce their finite-temperature polar properties well and that a value of the axial ratio close to that of the ideal WZ structure implies a large polarization at low temperatures but not necessarily at high temperatures because of the ordered-disordered character of the temperature-induced formation of the WZ state. Such findings should allow for a better fundamental understanding of (Al,Sc)N ferroelectric nitrides, which may be used to design efficient devices having, e.g., low operating voltages.