Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userBypassing the malfunction junction in warm dense matter simulations
473
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Simulation of warm dense matter requires computational methods that capture both quantum and classical behavior efficiently under high-temperature, high-density conditions. Currently, density functional theory molecular dynamics is used to model electrons and ions, but this methods computational cost skyrockets as temperatures and densities increase. We propose finite-temperature potential functional theory as an in-principle-exact alternative that suffers no such drawback. We derive an orbital-free free energy approximation through a coupling-constant formalism. Our density approximation and its associated free energy approximation demonstrate the methods accuracy and efficiency.
rate research
Read More
In a recent Letter [T.~Dornheim emph{et al.}, Phys.~Rev.~Lett.~textbf{125}, 085001 (2020)], we have presented the first emph{ab initio} results for the nonlinear density response of electrons in the warm dense matter regime. In the present work, we extend these efforts by carrying out extensive new path integral Monte Carlo (PIMC) simulations of a emph{ferromagnetic} electron gas that is subject to an external harmonic perturbation. This allows us to unambiguously quantify the impact of spin-effects on the nonlinear density response of the warm dense electron gas. In addition to their utility for the description of warm dense matter in an external magnetic field, our results further advance our current understanding of the uniform electron gas as a fundamental model system, which is important in its own right.
The calculation of chemical potential has traditionally been a challenge in atomistic simulations. One of the most used approaches is Widoms insertion method in which the chemical potential is calculated by periodically attempting to insert an extra particle in the system. In dense systems this method fails since the insertion probability is very low. In this paper we show that in a homogeneous fluid the insertion probability can be increased using metadynamics. We test our method on a supercooled high density binary Lennard-Jones fluid. We find that we can obtain efficiently converged results even when Widoms method fails.
Warm dense matter (WDM) -- an exotic state of highly compressed matter -- has attracted high interest in recent years in astrophysics and for dense laboratory systems. At the same time, this state is extremely difficult to treat theoretically. This is due to the simultaneous appearance of quantum degeneracy, Coulomb correlations and thermal effects, as well as the overlap of plasma and condensed phases. Recent breakthroughs are due to the successful application of density functional theory (DFT) methods which, however, often lack the necessary accuracy and predictive capability for WDM applications. The situation has changed with the availability of the first textit{ab initio} data for the exchange-correlation free energy of the warm dense uniform electron gas (UEG) that were obtained by quantum Monte Carlo (QMC) simulations, for recent reviews, see Dornheim textit{et al.}, Phys. Plasmas textbf{24}, 056303 (2017) and Phys. Rep. textbf{744}, 1-86 (2018). In the present article we review recent further progress in QMC simulations of the warm dense UEG: namely, textit{ab initio} results for the static local field correction $G(q)$ and for the dynamic structure factor $S(q,omega)$. These data are of key relevance for the comparison with x-ray scattering experiments at free electron laser facilities and for the improvement of theoretical models. In the second part of this paper we discuss simulations of WDM out of equilibrium. The theoretical approaches include Born-Oppenheimer molecular dynamics, quantum kinetic theory, time-dependent DFT and hydrodynamics. Here we analyze strengths and limitations of these methods and argue that progress in WDM simulations will require a suitable combination of all methods. A particular role might be played by quantum hydrodynamics, and we concentrate on problems, recent progress, and possible improvements of this method.
We study nonideal mixing effects in the regime of warm dense matter (WDM) by computing the shock Hugoniot curves of BN, MgO, and MgSiO_3. First, we derive these curves from the equations of state (EOS) of the fully interacting systems, which were obtained using a combination of path integral Monte Carlo calculations at high temperature and density functional molecular dynamics simulations at lower temperatures. We then use the ideal mixing approximation at constant pressure and temperature to rederive these Hugoniot curves from the EOS tables of the individual elements. We find that the linear mixing approximation works remarkably well at temperatures above ~2*10^5 K, where the shock compression ratio exceeds ~3.2. The shape of the Hugoniot curve of each compound is well reproduced. Regions of increased shock compression, that emerge because of the ionization of L and K shell electrons, are well represented and the maximum compression ratio on the Hugoniot curves is reproduced with high precision. Some deviations are seen near the onset of the L shell ionization regime, where ionization equilibrium in the fully interacting system cannot be well reproduced by the ideal mixing approximation. This approximation also breaks down at lower temperatures, where chemical bonds play an increasingly import role. However, the results imply that equilibrium properties of binary and ternary mixtures in the regime of WDM can be derived from the EOS tables of the individual elements. This significantly simplifies the characterization of binary and ternary mixtures in the WDM and plasma phases, which otherwise requires large numbers of more computationally expensive first-principles computer simulations.
Exploring and understanding ultrafast processes at the atomic level is a scientific challenge. Femtosecond X-ray Absorption Spectroscopy (XAS) is an essential experimental probing technic, as it can simultaneously reveal both electronic and atomic structures, and thus unravel their non-equilibrium dynamic interplay which is at the origin of most of the ultrafast mechanisms. However, despite considerable efforts, there is still no femtosecond X-ray source suitable for routine experiments. Here we show that betatron radiation from relativistic laser-plasma interaction combines ideal features for femtosecond XAS. It has been used to investigate the non-equilibrium transition of a copper sample brought at extreme conditions of temperature and pressure by a femtosecond laser pulse. We measured a rise time of the electron temperature below 100 fs. This first experiment demonstrates the great potential of the betatron source and paves the way to a new class of ultrafast experiments.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
