Understanding warm dense matter: from theory to experiment
Explore the source record for details and available documents.
Physics - Plasma physics↗
Engineering topics
Dornheim, Tobias
Publications and source records attributed to Dornheim, Tobias.
Towards a new era for warm dense matter diagnostics
Explore the source record for details and available documents.
Physics - Plasma physics↗
Ab initio path integral Monte Carlo simulation of warm dense matter
Explore the source record for details and available documents.
Physics - Plasma physics↗
Roadmap for warm dense matter physics
This roadmap presents the state-of-the-art, current challenges and near future developments anticipated in the thriving field of warm dense matter physics. Originating from strongly coupled plasma physics, high pressure physics and high energy density science, the warm dense matter physics community has recently taken a giant leap forward. This is due to spectacular developments in laser technology, diagnostic capabilities, and computer simulation techniques. Only in the last decade has it become possible to perform accurate enough simulations \& experiments to truly verify theoretical results as well as to reliably design experiments based on predictions. Consequently, this roadmap discusses recent developments and contemporary challenges that are faced by theoretical methods, and experimental techniques needed to create and diagnose warm dense matter. A large part of this roadmap is dedicated to specific warm dense matter systems and applications in astrophysics, inertial confinement fusion and novel material synthesis.
FOS: Physical sciences↗
Toward first principles-based simulations of dense hydrogen
Accurate knowledge of the properties of hydrogen at high compression is crucial for astrophysics (e.g., planetary and stellar interiors, brown dwarfs, atmosphere of compact stars) and laboratory experiments, including inertial confinement fusion. There exists experimental data for the equation of state, conductivity, and Thomson scattering spectra. However, the analysis of the measurements at extreme pressures and temperatures typically involves additional model assumptions, which makes it difficult to assess the accuracy of the experimental data rigorously. On the other hand, theory and modeling have produced extensive collections of data. They originate from a very large variety of models and simulations including path integral Monte Carlo (PIMC) simulations, density functional theory (DFT), chemical models, machine-learned models, and combinations thereof. At the same time, each of these methods has fundamental limitations (fermion sign problem in PIMC, approximate exchange–correlation functionals of DFT, inconsistent interaction energy contributions in chemical models, etc.), so for some parameter ranges accurate predictions are difficult. Recently, a number of breakthroughs in first principles PIMC as well as in DFT simulations were achieved which are discussed in this review. Here we use these results to benchmark different simulation methods. We present an update of the hydrogen phase diagram at high pressures, the expected phase transitions, and thermodynamic properties including the equation of state and momentum distribution. Furthermore, we discuss available dynamic results for warm dense hydrogen, including the conductivity, dynamic structure factor, plasmon dispersion, imaginary-time structure, and density response functions. We conclude by outlining strategies to combine different simulations to achieve accurate theoretical predictions that are based on first principles.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Effects of mosaic crystal instrument functions on x-ray Thomson scattering diagnostics
Mosaic crystals, with their high integrated reflectivities, are widely employed in spectrometers used to diagnose high energy density systems. X-ray Thomson scattering (XRTS) has emerged as a powerful diagnostic tool of these systems, providing in principle direct access to important properties such as the temperature via detailed balance. However, the measured XRTS spectrum is broadened by the spectrometer instrument function (IF), and without careful consideration of the IF one risks misdiagnosing system conditions. Here, we consider in detail the IF of 40 and 100 μm mosaic Highly Annealed Pyrolytic Graphite crystals, and how the broadening varies across the spectrometer in an energy range of 6.7–8.6 keV. Notably, we find a strong asymmetry in the shape of the IF toward higher energies. As an example, we consider the effect of the asymmetry in the IF on the temperature inferred via XRTS for simulated 80 eV CH plasmas and find that the temperature can be overestimated if an approximate symmetric IF is used. We, therefore, expect a detailed consideration of the full IF will have an important impact on system properties inferred via XRTS in both forward modeling and model-free approaches.
47 OTHER INSTRUMENTATION↗
Imposing correct jellium response is key to predict the density response by orbital-free DFT
Orbital-free density functional theory constitutes a computationally highly effective tool for modeling electronic structures of systems ranging from room-temperature materials to warm dense matter. Its accuracy critically depends on the employed kinetic energy (KE) density functional, which has to be supplied as an external input. In this work we consider several nonlocal and Laplacian-level KE functionals and use an external harmonic perturbation to compute the static density response at T = 0 K in the linear and beyond-linear response regimes. We test for the satisfaction of exact conditions in the limit of uniform densities and for how approximate KE functionals reproduce the density response of realistic materials (e.g., Al and Si) against the Kohn-Sham DFT reference, which employs the exact KE. The results illustrate that several functionals violate exact conditions in the uniform electron gas (UEG) limit. We find a strong correlation between the accuracy of the KE functionals in the UEG limit and in the strongly inhomogeneous case. This empirically demonstrates the importance of imposing the limit of UEG response for uniform densities and validates the use of the Lindhard function in the formulation of kernels for nonlocal functionals. This conclusion is substantiated by additional calculations for bulk aluminum (Al) with a face-centered cubic (fcc) lattice and silicon (Si) with an fcc lattice, body-centered cubic (bcc) lattice, and semiconducting crystal diamond state. The analysis of fcc Al, and fcc as well as bcc Si data follows closely the conclusions drawn for the UEG, allowing us to extend our conclusions to realistic systems that are subject to density inhomogeneities induced by ions.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Imaginary-time correlation function thermometry: A new, high-accuracy and model-free temperature analysis technique for x-ray Thomson scattering data
The accurate interpretation of experiments with matter at extreme densities and pressures is a notoriously difficult challenge. In a recent work [Dornheim et al., Nat. Commun. 13, 7911 (2022)], we have introduced a formally exact methodology that allows extracting the temperature of arbitrary complex materials without any model assumptions or simulations. Here, we provide a more detailed introduction to this approach and analyze the impact of experimental noise on the extracted temperatures. In particular, we extensively apply our method both to synthetic scattering data and to previous experimental measurements over a broad range of temperatures and wave numbers. We expect that our approach will be of high interest to a gamut of applications, including inertial confinement fusion, laboratory astrophysics, and the compilation of highly accurate equation-of-state databases.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Electronic density response of warm dense matter
Matter at extreme temperatures and pressures—commonly known as warm dense matter (WDM)—is ubiquitous throughout our Universe and occurs in astrophysical objects such as giant planet interiors and brown dwarfs. Moreover, WDM is very important for technological applications such as inertial confinement fusion and is realized in the laboratory using different techniques. A particularly important property for the understanding of WDM is given by its electronic density response to an external perturbation. Such response properties are probed in x-ray Thomson scattering (XRTS) experiments and are central for the theoretical description of WDM. In this work, we give an overview of a number of recent developments in this field. To this end, we summarize the relevant theoretical background, covering the regime of linear response theory and nonlinear effects, the fully dynamic response and its static, time-independent limit, and the connection between density response properties and imaginary-time correlation functions (ITCF). In addition, we introduce the most important numerical simulation techniques, including path-integral Monte Carlo simulations and different thermal density functional theory (DFT) approaches. From a practical perspective, we present a variety of simulation results for different density response properties, covering the archetypal model of the uniform electron gas and realistic WDM systems such as hydrogen. Moreover, we show how the concept of ITCFs can be used to infer the temperature from XRTS measurements of arbitrary complex systems without the need for any models or approximations. Finally, we outline a strategy for future developments based on the close interplay between simulations and experiments.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Accurate temperature diagnostics for matter under extreme conditions
The experimental investigation of matter under extreme densities and temperatures, as in astrophysical objects and nuclear fusion applications, constitutes one of the most active frontiers at the interface of material science, plasma physics, and engineering. The central obstacle is given by the rigorous interpretation of the experimental results, as even the diagnosis of basic parameters like the temperature T is rendered difficult at these extreme conditions. Here, we present a simple, approximation-free method to extract the temperature of arbitrarily complex materials in thermal equilibrium from X-ray Thomson scattering experiments, without the need for any simulations or an explicit deconvolution. Our paradigm can be readily implemented at modern facilities and corresponding experiments will have a profound impact on our understanding of warm dense matter and beyond, and open up a variety of appealing possibilities in the context of thermonuclear fusion, laboratory astrophysics, and related disciplines.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Towards a quantum fluid theory of correlated many-fermion systems from first principles
Correlated many-fermion systems emerge in a broad range of phenomena in warm dense matter, plasmonics, and ultracold atoms. Quantum hydrodynamics (QHD) complements first-principles methods for many-fermion systems at larger scales. We illustrate the failure of the standard Bohm potential central to QHD for strong perturbations when the density perturbation is larger than about 10^{-3} 10 − 3 of the mean density. We then extend QHD to this regime via the many-fermion Bohm potential from first-principles. This may lead to more accurate QHD simulations beyond their common application domain in the presence of strong perturbations at scales unattainable with first-principles methods.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Path integral Monte Carlo approach to the structural properties and collective excitations of liquid $$^3{\text {He}}$$ without fixed nodes
Abstract Due to its nature as a strongly correlated quantum liquid, ultracold helium is characterized by the nontrivial interplay of different physical effects. Bosonic $$^4{\text {He}}$$ 4 He exhibits superfluidity and Bose-Einstein condensation. Its physical properties have been accurately determined on the basis of ab initio path integral Monte Carlo (PIMC) simulations. In contrast, the corresponding theoretical description of fermionic $$^3{\text {He}}$$ 3 He is severely hampered by the notorious fermion sign problem, and previous PIMC results have been derived by introducing the uncontrolled fixed-node approximation. In this work, we present extensive new PIMC simulations of normal liquid $$^3{\text {He}}$$ 3 He without any nodal constraints. This allows us to to unambiguously quantify the impact of Fermi statistics and to study the effects of temperature on different physical properties like the static structure factor $$S({\mathbf {q}})$$ S ( q ) , the momentum distribution $$n({\mathbf {q}})$$ n ( q ) , and the static density response function $$\chi ({\mathbf {q}})$$ χ ( q ) . In addition, the dynamic structure factor $$S({\mathbf {q}},\omega )$$ S ( q , ω ) is rigorously reconstructed from imaginary-time PIMC data. From simulations of $$^3{\text {He}}$$ 3 He , we derived the familiar phonon–maxon–roton dispersion function that is well-known for $$^4{\text {He}}$$ 4 He and has been reported previously for two-dimensional $$^3{\text {He}}$$ 3 He films (Nature 483:576–579 (2012)). The comparison of our new results for both $$S({\mathbf {q}})$$ S ( q ) and $$S({\mathbf {q}},\omega )$$ S ( q , ω ) with neutron scattering measurements reveals an excellent agreement between theory and experiment.
97 MATHEMATICS AND COMPUTING↗
Momentum distribution of the uniform electron gas at finite temperature: Effects of spin polarization
We carry out extensive direct path integral Monte Carlo (PIMC) simulations of the uniform electron gas (UEG) at finite temperature for different values of the spin-polarization ξ. This allows us to unambiguously quantify the impact of spin effects on the momentum distribution function n(k) and related properties. We find that interesting physical effects like the interaction-induced increase in the occupation of the zero-momentum state n(0) substantially depend on ξ. Our results further advance the current understanding of the UEG as a fundamental model system, and are of practical relevance for the description of transport properties of warm dense matter in an external magnetic field. All PIMC results are freely available online and can be used as a benchmark for the development of methods and applications.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
First-principles modeling of plasmons in aluminum under ambient and extreme conditions
The theoretical understanding of plasmon behavior is crucial for an accurate interpretation of inelastic scattering diagnostics in many experiments. In this work, we highlight the utility of linear response time-dependent density functional theory (LR-TDDFT) as a first-principles framework for consistently modeling plasmon properties. We provide a comprehensive analysis of plasmons in aluminum from ambient to warm dense matter conditions and assess typical properties such as the dynamical structure factor, the plasmon dispersion, and the plasmon lifetime. We compare our results with scattering measurements and with other TDDFT results as well as models such as the random phase approximation, the Mermin approach, and the dielectric function obtained using static local field corrections of the uniform electron gas parametrized from path-integral Monte Carlo simulations. We conclude that results for the plasmon dispersion and lifetime are inconsistent between experiment and theories and that the common practice of extracting and studying plasmon dispersion relations is an insufficient procedure to capture the complicated physics contained in the dynamic structure factor in its full breadth.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗