Search

Using a novel ultrahigh resolution (\Delta E ~ 0.1eV) setup to measure electronic features in x-ray Thomson scattering (XRTS) experiments at the European XFEL in Germany, we have studied the collective plasmon excitation in aluminium at ambient conditions, which we can measure very accurately even at low momentum transfers. As a result, we can resolve previously reported discrepancies between ab initio time-dependent density functional theory simulations and experimental observations. The demonstrated capability for high-resolution XRTS measurements will be a game changer for the diagnosis of experiments with matter under extreme densities, temperatures, and pressures, and unlock the full potential of state-of-the-art x-ray free electron laser (XFEL) facilities to study planetary interior conditions, to understand inertial confinement fusion applications, and for material science and discovery.

The Perdew-Zunger (PZ) self-interaction correction (SIC) is an established tool to correct unphysical behavior in density functional approximations. Yet, PZ-SIC is well-known to sometimes break molecular symmetries. An example of this is the benzene molecule, for which PZ-SIC predicts a symmetry-broken electron density and molecular geometry, since the method does not describe the two possible Kekulé structures on an even footing, leading to local minima [Lehtola et al, J. Chem. Theory Comput. 2016, 12, 3195]. PZ-SIC is often implemented with Fermi-Löwdin orbitals (FLOs), yielding the FLO-SIC method, which likewise has issues with symmetry breaking and local minima [Trepte et al, J. Chem. Phys. 2021, 155, 224109]. In this work, we propose a generalization of PZ-SIC - the ensemble PZ-SIC (E-PZ-SIC) method - which shares the asymptotic computational scaling of PZ-SIC (albeit with an additional prefactor). E-PZ-SIC is straightforwardly applicable to various molecules, merely requiring one to average the self-interaction correction over all possible Kekulé structures, in line with chemical intuition. We showcase the implementation of E-PZ-SIC with FLOs, as the resulting E-FLO-SIC method is easy to realize on top of an existing implementation of FLO-SIC. We show that E-FLO-SIC indeed eliminates symmetry breaking, reproducing a symmetric electron density and molecular geometry for benzene. The ensemble approach suggested herein could also be employed within locally scaled variants of PZ-SIC and their FLO-SIC versions.

This repository contains all PIMC results related to the publication "Ab initio Density Response and Local Field Factor of Warm Dense Hydrogen". Generally, data formats are identical to figures. Exceptions are 3D ITCF data sets for Figs. 2, 8 and 12: #1 k [a_Bohr^{-1}], #2 tau [Ha^{-1}], #3/#4 F(q,tau)x32 and statistical error and the "ITCF" folders with the raw data for F(q,tau): ITCF: #1 tau [Ha^{-1}]; #2/3: F(q,tau) and statistical error The number after "index" in the file names gives the number of the respective q-vector; see "static_structure_factor_key.dat", columns 1 and 2 for the respective index-to-q mapping, with [q]=a_Bohr{-1}

The accurate description of non-ideal quantum many-body systems is of prime importance for a host of applications within physics, quantum chemistry, material science, and related disciplines. At finite temperatures, the gold standard is given by \textit{ab initio} path integral Monte Carlo (PIMC) simulations, which do not require any empirical input, but exhibit an exponential increase in the required compute time for fermionic systems with increasing the system size N. Very recently, it has been suggested to compute fermionic properties without this bottleneck based on PIMC simulations of fictitious identical particles. In the present work, we use this technique to carry out very large (N≤1000) PIMC simulations of the warm dense electron gas and demonstrate that it is capable of providing a highly accurate description of investigated properties, i.e., the static structure factor, the static density response function, and local field correction, over the entire range of length scales.

We present extensive new \emph{ab initio} path integral Monte Carlo (PIMC) results for a variety of structural properties of warm dense hydrogen and beryllium. To deal with the fermion sign problem -- an exponential computational bottleneck due to the antisymmetry of the electronic thermal density matrix -- we employ the recently proposed [\textit{J.~Chem.~Phys.}~\textbf{157}, 094112 (2022); \textbf{159}, 164113 (2023)] -extrapolation method and find excellent agreement with exact direct PIMC reference data where available. This opens up the intriguing possibility to study a gamut of properties of light elements and potentially material mixtures over a substantial part of the warm dense matter regime, with direct relevance for astrophysics, material science, and inertial confinement fusion research.

Hydrogen at extreme temperatures and pressures is ubiquitous throughout our universe and naturally occurs in a variety of astrophysical objects. In addition, it is of key relevance for cutting-edge technological applications, with inertial confinement fusion research being a prime example. In the present work, we present exact \emph{ab initio} path integral Monte Carlo (PIMC) results for the electronic density of warm dense hydrogen along a line of constant degeneracy across a broad range of densities. Using the well-known concept of reduced density gradients, we develop a new framework to identify the breaking of bound states due to pressure ionization in bulk hydrogen. Moreover, we use our PIMC results as a reference to rigorously assess the accuracy of a variety of exchange--correlation (XC) functionals in density functional theory calculations for different density regions. Here a key finding is the importance of thermal XC effects for the accurate description of density gradients in high-energy density systems. Our exact PIMC test set is freely available online and can be used to guide the development of new methodologies for the simulation of warm dense matter and beyond.

This repository contains the PIMC simulation results shown in the publication "Ab Initio Path Integral Monte Carlo Simulations of the Uniform Electron Gas on Large Length Scales"