Jellium Research Articles

Revealing quasi-excitations in the low-density homogeneous electron gas with model exchange–correlation kernels

Time-dependent density functional theory within the linear response regime provides a solid mathematical framework to capture excitations. The accuracy of the theory, however, largely depends on the approximations for the exchange–correlation (xc) kernels. Away from the long-wavelength (or q = 0 short wave-vector) and zero-frequency (ω = 0) limit, the correlation contribution to the kernel becomes more relevant and dominant over exchange. The dielectric function, in principle, can encompass xc effects relevant to describe low-density physics. Furthermore, besides collective plasmon excitations, the dielectric function can reveal collective electron–hole excitations, often dubbed “ghost excitons.” Besides collective excitons, the physics of the low-density regime is rich, as exemplified by a static charge-density wave that was recently found for rs > 69, and was shown to be associated with softening of the plasmon mode. These excitations are seen to be present in much higher density 2D homogeneous electron gases of rs ≳ 4. In this work, we perform a thorough analysis with xc model kernels for excitations of various nature. The uniform electron gas, as a useful model of real metallic systems, is used as a platform for our analysis. We highlight the relevance of exact constraints as we display and explain screening and excitations in the low-density region.

Static Self-Energy and Effective Mass of the Homogeneous Electron Gas from Quantum Monte Carlo Calculations

Static Self-Energy and Effective Mass of the Homogeneous Electron Gas from Quantum Monte Carlo Calculations

Fermionic physics from abinitio path integral Monte Carlo simulations of fictitious identical particles.

The abinitio path integral Monte Carlo (PIMC) method is one of the most successful methods in statistical physics, quantum chemistry and related fields, but its application to quantum degenerate Fermi systems is severely hampered by an exponential computational bottleneck: the notorious fermion signproblem. Very recently, Xiong and Xiong [J. Chem. Phys. 157, 094112 (2022)] have suggested to partially circumvent the signproblem by carrying out simulations of fictitious systems guided by an interpolating continuous variable ξ ∈ [-1, 1], with the physical Fermi- and Bose-statistics corresponding to ξ = -1 and ξ = 1. It has been proposed that information about the fermionic limit might be obtained by calculations within the bosonic sector ξ > 0 combined with an extrapolation throughout the fermionic sector ξ < 0, essentially bypassing the signproblem. Here, we show how the inclusion of the artificial parameter ξ can be interpreted as an effective penalty on the formation of permutation cycles in the PIMC simulation. We demonstrate that the proposed extrapolation method breaks down for moderate to high quantum degeneracy. Instead, the method constitutes a valuable tool for the description of large Fermi-systems of weak quantum degeneracy. This is demonstrated for electrons in a 2D harmonic trap and for the uniform electron gas (UEG), where we find excellent agreement (∼0.5%) with exact configuration PIMC results in the high-density regime while attaining a speed-up exceeding 11 orders of magnitude. Finally, we extend the idea beyond the energy and analyze the radial density distribution (2D trap), as well as the static structure factor and imaginary-time density-density correlation function (UEG).

The electron-centric approach to the exchange-correlation energy.

The Kohn-Sham theory addresses the challenge of representing the kinetic energy by re-quantizing density functional theory at a level of non-interacting electrons. It transforms the many-electron problem into a fictitious non-interacting electron problem, with the many-electron effects concealed within the exchange-correlation (XC) energy, which is expressed in terms of the electron density ρ(r). Unlike the wave function, ρ(r) can be viewed as a classical quantity, and expressing the XC energy in terms of it circumvents the need for correlated wave functions. In this work, we once again employ the re-quantization strategy and determine the XC energy using a local one-particle Schrödinger equation. The ground-state eigenfunction of the corresponding Hamiltonian is a reference point (r) dependent orbital φr,σ(u, σ') which is subsequently used to generate the XC hole and the XC energy. The spin coordinate is denoted by σ and u is the electron-electron separation. The one-particle equation for φr,σ(u, σ') includes a local potential vr,σ(u, σ') that we approximate using two simple physical constraints. We assess the approximation by applying it to the helium iso-electronic series, the homogeneous electron gas, and the dissociation of the hydrogen molecule.

Quasiparticle approach to collective quantum dielectric response

In current research, we use a generalized quantum multistream model to develop an effective quasiparticle theory for quantum many-body effects. The N-electron Schrödinger–Poisson stream model is reduced to a system of coupled differential equations with new wavefunction representation for collective quantum excitations in the many electron system. The current theory is then applied to the collective quantum statistical behavior of homogenous electron gas. Moreover, the generalized energy dispersion relation, which incorporates the quasiparticle band structure, is used to calculate the linear dielectric response of collective quantum excitations in the electron gas with arbitrary degree of degeneracy beyond many-body theories, limiting assumptions such as the independent electron and the random phase approximations. Important parameters of electron gas such as the dynamic structure factor, the loss function, the static charge screening, optical reflectivity, and the electronic stopping power are investigated as applications of current theory. The quasiparticle theory incorporates effects both due to single-electron excitations as well as the electrostatic interaction among electrons in a single picture. Existence of Van-Hove-like singularity at the plasmon wavenumber leads to distinct features of quasiparticle response to electromagnetic perturbations in the electron gas. It is shown that collective quantum excitations in high density electron gas below a given critical electron temperature are blocked due to existence of a large quasiparticle energy bandgap above the Fermi level. A new equation of states is given for the quasiparticle excitation in the electron gas, based on the transition probability of electrons to the quasiparticle level. It is found that, the screening potential of a static charge in quasiparticle model has an oscillatory Lennard–Jones-type attractive form.

Electronic density response of warm dense hydrogen on the nanoscale.

The properties of hydrogen at warm dense matter (WDM) conditions are of high importance for the understanding of astrophysical objects and technological applications such as inertial confinement fusion. In this work, we present extensive ab initio path integral Monte Carlo results for the electronic properties in the Coulomb potential of a fixed ionic configuration. This gives us unique insights into the complex interplay between the electronic localization around the protons with their density response to an external harmonic perturbation. We find qualitative agreement between our simulation data and a heuristic model based on the assumption of a local uniform electron gas model, but important trends are not captured by this simplification. In addition to being interesting in their own right, we are convinced that our results will be of high value for future projects, such as the rigorous benchmarking of approximate theories for the simulation of WDM, most notably density functional theory.

Nonlocal Functionals Inspired by the Strongly Interacting Limit of DFT: Exact Constraints and Implementation.

Capturing strong correlation effects remains a key challenge for the development of improved exchange-correlation (XC) functionals in density functional theory. The recently proposed multiple radii functional (MRF) [J. Phys. Chem. Lett. 2017, 8, 2799; J. Chem. Theory Comput. 2019, 15, 3580] was designed to capture strong correlation effects seamlessly, as its mathematical structure draws from that of the exact XC functional in the limit of infinite correlations. The MRF functional provides a framework for building approximations along the density-fixed adiabatic connection, delivers accurate XC energy densities in the standard DFT gauge (same as that of the exact exchange energy density), and is free of one-electron self-interaction errors. To facilitate the development of XC functionals based on the MRF, we examine the behavior of the MRF functional when applied to uniform and scaled densities and consider how it can be made exact for the uniform electron gas. These theoretical insights are then used to build improved forms for the fluctuation function, an object that defines XC energy densities within the MRF framework. We also show how the MRF fluctuation function for physical correlation can be easily readjusted to accurately capture the XC functional in the limit of infinite correlations, demonstrating the versatility of MRF for building approximations for different correlation regimes. We describe the implementation of MRF using densities expanded on Gaussian basis sets, which improves the efficiency of previous grid-based MRF implementations.

Physical insights from imaginary-time density–density correlation functions

An accurate theoretical description of the dynamic properties of correlated quantum many-body systems, such as the dynamic structure factor S(q, ω), is important in many fields. Unfortunately, highly accurate quantum Monte Carlo methods are usually restricted to the imaginary time domain, and the analytic continuation of the imaginary-time density–density correlation function F(q, τ) to real frequencies is a notoriously hard problem. Here, it is argued that often no such analytic continuation is required because by definition, F(q, τ) contains the same physical information as does S(q, ω), only represented unfamiliarly. Specifically, it is shown how one can directly extract key information such as the temperature or quasi-particle excitation energies from the τ domain, which is highly relevant for equation-of-state measurements of matter under extreme conditions [T. Dornheim et al., Nat. Commun. 13, 7911 (2022)]. As a practical example, ab initio path-integral Monte Carlo results for the uniform electron gas (UEG) are considered, and it is shown that even nontrivial processes such as the roton feature of the UEG at low density [T. Dornheim et al., Commun. Phys. 5, 304 (2022)] are manifested straightforwardly in F(q, τ). A comprehensive overview is given of various useful properties of F(q, τ) and how it relates to the usual dynamic structure factor. In fact, working directly in the τ domain is advantageous for many reasons and opens up multiple avenues for future applications.

Analysing the dynamic structure of warm dense matter in the imaginary-time domain: theoretical models and simulations.

Rigorous diagnostics of experiments with warm dense matter are notoriously difficult. A key method is X-ray Thomson scattering (XRTS), but the interpretation of XRTS measurements is usually based on theoretical models that entail various approximations. Recently, Dornheim et al. [Nat. Commun. 13, 7911 (2022)] introduced a new framework for temperature diagnostics of XRTS experiments that is based on imaginary-time correlation functions. On the one hand, switching from the frequency to the imaginary-time domain gives one direct access to a number of physical properties, which facilitates the extraction of the temperature of arbitrarily complex materials without relying on any models or approximations. On the other hand, the bulk of theoretical work in dynamic quantum many-body theory is devoted to the frequency domain, and, to the best of our knowledge, the manifestation of physics properties within the imaginary-time density-density correlation function (ITCF) remains poorly understood. In the present work, we aim to fill this gap by introducing a simple, semi-analytical model for the imaginary-time dependence of two-body correlations within the framework of imaginary-time path integrals. As a practical example, we compare our new model to extensive ab initio path integral Monte Carlo results for the ITCF of a uniform electron gas, and find excellent agreement over a broad range of wavenumbers, densities and temperatures. This article is part of the theme issue 'Dynamic and transient processes in warm dense matter'.

Gaussian expansion of Yukawa non‐local kinetic energy functionals: Application to metal clusters

AbstractThe development of kinetic energy (KE) functionals is one of the current challenges in density functional theory (DFT). The Yukawa non‐local KE functionals [Phys. Rev. B 103, 155127 (2021)] have been shown to describe accurately the Lindhard response of the homogeneous electron gas (HEG) directly in the real space, without any step in the reciprocal space. However, the Yukawa kernel employs an exponential function which cannot be efficiently represented in conventional Gaussian‐based quantum chemistry codes. Here, we present an expansion of the Yukawa kernel in Gaussian functions. We show that for the HEG this expansion is independent of the electronic density, and that for general finite systems the accuracy can be easily tuned. Finally, we present results for atomistic sodium clusters of different sizes, showing that simple Yukawa functionals can give superior accuracy as compared to semilocal functionals.

Neural network ansatz for periodic wave functions and the homogeneous electron gas

The authors introduce here WAPNET, a periodic neural network (NN) variational ansatz for solving the ground state of a homogeneous electron gas with high accuracy over a broad range of the density coupling constant ${r}_{s}$. Going beyond recent work for molecules, this contribution establishes NN models as flexible and powerful ansatz for electronic structure calculations in extended systems. In all density regimes, WAPNET-based variational Monte Carlo results are comparable to or better than state-of-the-art benchmarks obtained by diffusion Monte Carlo with iterative backflow.

$m^*$ of two-dimensional electron gas: A neural canonical transformation study

The quasiparticle effective mass m^\astm* of interacting electrons is a fundamental quantity in the Fermi liquid theory. However, the precise value of the effective mass of uniform electron gas is still elusive after decades of research. The newly developed neural canonical transformation approach [Xie et al., J. Mach. Learn. 1, (2022)] offers a principled way to extract the effective mass of electron gas by directly calculating the thermal entropy at low temperature. The approach models a variational many-electron density matrix using two generative neural networks: an autoregressive model for momentum occupation and a normalizing flow for electron coordinates. Our calculation reveals a suppression of effective mass in the two-dimensional spin-polarized electron gas, which is more pronounced than previous reports in the low-density strong-coupling region. This prediction calls for verification in two-dimensional electron gas experiments.

QMC-consistent static spin and density local field factors for the uniform electron gas

Analytic mathematical models for the static spin (${G}_{\ensuremath{-}}$) and density (${G}_{+}$) local field factors for the uniform electron gas (UEG) as functions of wave vector and density are presented. These models closely fit recent quantum Monte Carlo (QMC) data and satisfy exact asymptotic limits. A simple functional form for ${G}_{\ensuremath{-}}$ is developed; the same functional form parametrized for ${G}_{+}$ yields an improvement over previous work. The QMC-computed ${G}_{\ifmmode\pm\else\textpm\fi{}}$ are consistent with a rapid crossover between theoretically derived small-$q$ and large-$q$ expansions of ${G}_{\ifmmode\pm\else\textpm\fi{}}$. These expansions are completely determined by ${r}_{\mathrm{s}}$, the UEG correlation energy per electron, and the UEG on-top pair distribution function. We demonstrate their utility by computing uniform electron gas correlation energies over a range of densities. These models, which hold over an extremely wide range of densities, are recommended for use in practical time-dependent density functional theory calculations of simple metallic systems. A revised model of the spin susceptibility enhancement is developed that fits QMC data, and does not show a ferromagnetic instability at low density.

New features of the TFmix code: Thermodynamic properties of electrons in mixtures

New features of the TFmix code: Thermodynamic properties of electrons in mixtures

Complexity of Implementing Trotter Steps

Quantum dynamics can be simulated on a quantum computer by exponentiating elementary terms from the Hamiltonian in a sequential manner. However, such an implementation of Trotter steps has gate complexity depending on the total Hamiltonian term number, comparing unfavorably to algorithms using more advanced techniques. We develop methods to perform faster Trotter steps with complexity sublinear in the number of terms. We achieve this for a class of Hamiltonians whose interaction strength decays with distance according to power law. Our methods include one based on a recursive block encoding and one based on an average-cost simulation, overcoming the normalization-factor barrier of these advanced quantum simulation techniques. We also realize faster Trotter steps when certain blocks of Hamiltonian coefficients have low rank. Combining with a tighter error analysis, we show that it suffices to use $\left(\eta^{1/3}n^{1/3}+\frac{n^{2/3}}{\eta^{2/3}}\right)n^{1+o(1)}$ gates to simulate uniform electron gas with $n$ spin orbitals and $\eta$ electrons in second quantization in real space, asymptotically improving over the best previous work. We obtain an analogous result when the external potential of nuclei is introduced under the Born-Oppenheimer approximation. We prove a circuit lower bound when the Hamiltonian coefficients take a continuum range of values, showing that generic $n$-qubit $2$-local Hamiltonians with commuting terms require at least $\Omega(n^2)$ gates to evolve with accuracy $\epsilon=\Omega(1/poly(n))$ for time $t=\Omega(\epsilon)$. Our proof is based on a gate-efficient reduction from the approximate synthesis of diagonal unitaries within the Hamming weight-$2$ subspace, which may be of independent interest. Our result thus suggests the use of Hamiltonian structural properties as both necessary and sufficient to implement Trotter steps with lower gate complexity.

Unravelling the nonlinear ideal density response of many-body systems

Nonlinear density response theory is revisited focusing on the harmonically perturbed finite temperature uniform electron gas. Within the non-interacting limit, brute force quantum kinetic theory calculations for the quadratic, cubic, quartic and quintic responses reveal a deep connection with the linear response. Careful analysis of the static long wavelength limit led us to conjecture a canonical non-interacting form that expresses arbitrary order nonlinear responses as the weighted sum of the linear responses evaluated at all multiple harmonics. This harmonic expansion is successfully validated against ab initio path integral Monte Carlo simulations.

Energy response and spatial alignment of the perturbed electron gas.

We study the linear energy response of the uniform electron gas to an external harmonic perturbation with a focus on resolving different contributions to the total energy. This has been achieved by carrying out highly accurate abinitio path integral Monte Carlo (PIMC) calculations for a variety of densities and temperatures. We report a number of physical insights into effects such as screening and the relative importance of kinetic and potential energies for different wave numbers. A particularly interesting finding is obtained from the observed non-monotonic behavior of the induced change in the interaction energy, which becomes negative for intermediate wave numbers. This effect is strongly dependent on the coupling strength and constitutes further direct evidence for the spatial alignment of electrons introduced in earlier works [T. Dornheim et al., Commun. Phys. 5, 304 (2022)]. The observed quadratic dependence on the perturbation amplitude in the limit of weak perturbations and the quartic dependence of perturbation amplitude corrections are consistent with linear and nonlinear versions of the density stiffness theorem. All PIMC simulation results are freely available online and can be used to benchmark new methods or as input for other calculations.

Electronic interaction of slow hydrogen, helium, nitrogen, and neon ions with silicon

We investigate the electronic excitation of silicon by hydrogen, helium, nitrogen, and neon ions for ion energies ranging from several tens to a few hundred kiloelectronvolts. Experiments are carried out in transmission geometry using a time-of-flight medium energy ion scattering system. The targets are self-supporting, single-crystalline Si (100) foils with nominal thicknesses of 50 and 200 nm. Stopping cross-sections (SCSs) are derived and compared with datasets available from the literature and predictions from theory. The results for H projectiles reveal good agreement with literature datasets, within quoted uncertainties. For He projectiles, the results show good agreement with most of the literature data. For N ions, higher values than reported in the literature are measured. For Ne, where literature data are scarce, we extend the velocity regime for which data exist by a factor of two toward higher velocities. The electronic SCS is found to be proportional to ion velocity for all impinging ions, for velocities below the Bohr velocity $(v&lt;{v}_{0})$. Comparison with theoretical predictions for a homogeneous free electron gas indicates strong contributions of local energy loss processes different from those expected for electron-hole pair excitation in binary collisions for all ions different from protons.

Time-dependent exchange-correlation hole and potential of the electron gas

The exchange-correlation hole and potential of the homogeneous electron gas have been investigated within the random-phase approximation, employing the plasmon-pole approximation for the linear density response function. The angular dependence as well as the time dependence of the exchange-correlation hole are illustrated for a Wigner-Seitz radius ${r}_{s}=4$ (atomic unit). It is found that there is a substantial cancellation between exchange and correlation potentials in space and time, analogous to the cancellation of exchange and correlation self-energies. Analysis of the sum rule explains why it is more advantageous to use a noninteracting Green function than a renormalized one when calculating the response function within the random-phase approximation and consequently the self-energy within the well-established $GW$ approximation. The present study provides a starting point for more accurate and comprehensive calculations of the exchange-correlation hole and potential of the electron gas with the aim of constructing a model based on the local density approximation as in density functional theory.

Correlation energy of the paramagnetic electron gas at the thermodynamic limit

The variational and diffusion quantum Monte Carlo methods are used to calculate the correlation energy of the paramagnetic three-dimensional homogeneous electron gas at intermediate to high density. Ground-state energies in finite cells are determined using Slater-Jastrow-backflow trial wave functions, and finite-size errors are removed using twist-averaged boundary conditions and extrapolation of the energy per particle to the thermodynamic limit of infinite system size. Our correlation energies in the thermodynamic limit are more accurate than previous results. The present diffusion quantum Monte Carlo energies together with our recently reported [Phys. Rev. B 105, 245135 (2022)] results at low density, are used to parametrize the correlation energy of the electron gas using a functional form that satisfies the exact asymptotic behavior at high density.