Finite Temperature Green's Function Approach for Excited State and Thermodynamic Properties of Cool to Warm Dense Matter

Electronic Free Energy Surface of the Nitrogen Dimer Using First-Principles Finite Temperature Electronic Structure Methods

We use full configuration interaction and density matrix quantum Monte Carlo methods to calculate the electronic free energy surface of the nitrogen dimer within the free-energy Born-Oppenheimer approximation. As the temperature is raised from T = 0, we find a temperature regime in which the internal energy causes bond strengthening. At these temperatures, adding in the entropy contributions is required to cause the bond to gradually weaken with increasing temperature. We predict a thermally driven dissociation for the nitrogen dimer between 22,000 to 63,200 K depending on symmetries and basis set. Inclusion of more spatial and spin symmetries reduces the temperature required. The origin of these observations is explored using the structure of the density matrix at various temperatures and bond lengths.

Ab initio calculation of X-ray and related core-level spectroscopies: Green's function approaches

X-Ray and related spectroscopies are powerful probes of atomic, vibrational, and electronic structure. In order to unlock the full potential of such experimental techniques, accurate and efficient theoretical and computational approaches are essential. Here we review the status of a variety of first-principles and nearly first principles techniques for X-ray spectroscopies such as X-ray absorption, X-ray emission, and X-ray photoemission, with a focus on Green's function based methods. In particular, we describe the current state of multiple scattering Green's function techniques available in the FEFF10 code and cumulant Green's function techniques for including the effects of many-body electronic excitations. Illustrative examples are shown for a variety of materials and compared with other theoretical and experimental results.

Piecewise Interaction Picture Density Matrix Quantum Monte Carlo

The density matrix quantum Monte Carlo (DMQMC) set of methods stochastically samples the exact $N$-body density matrix for interacting electrons at finite temperature. We introduce a simple modification to the interaction picture DMQMC method (IP-DMQMC) which overcomes the limitation of only sampling one inverse temperature point at a time, instead allowing for the sampling of a temperature range within a single calculation thereby reducing the computational cost. At the target inverse temperature, instead of ending the simulation, we incorporate a change of picture away from the interaction picture. The resulting equations of motion have piecewise functions and use the interaction picture in the first phase of a simulation, followed by the application of the Bloch equation once the target inverse temperature is reached. We find that the performance of this method is similar to or better than the DMQMC and IP-DMQMC algorithms in a variety of molecular test systems.

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.

Equation-of-Motion Coupled-Cluster Cumulant Green's Function for Excited States and X-Ray Spectra

Green’s function methods provide a robust, general framework within many-body theory for treating electron correlation in both excited states and x-ray spectra. Conventional methods using the Dyson equation or the cumulant expansion are typically based on the GW self-energy approximation. In order to extend this approximation in molecular systems, a non-perturbative real-time coupled-cluster cumulant Green’s function approach has been introduced, where the cumulant is obtained as the solution to a system of coupled first order, non-linear differential equations. This approach naturally includes non-linear corrections to conventional cumulant Green’s function techniques where the cumulant is linear in the GW self-energy. The method yields the spectral function for the core Green’s function, which is directly related to the x-ray photoemission spectra (XPS) of molecular systems. The approach also yields very good results for binding energies and satellite excitations. The x-ray absorption spectrum (XAS) is then calculated using a convolution of the core spectral function and an effective, one-body XAS. Here this approach is extended to include the full coupled-cluster-singles (CCS) core Green’s function by including the complete form of the non-linear contributions to the cumulant as well as all single, double, and triple cluster excitations in the CC amplitude equations. This approach naturally builds in orthogonality and shake-up effects analogous to those in the Mahan-Noizeres-de Dominicis edge singularity corrections that enhance the XAS near the edge. The method is illustrated for the XPS and XAS of NH₃.

Studying the lifetime of charge and heat carriers due to intrinsic scattering mechanisms in FeVSb half-Heusler thermoelectric

This work, presents a study of lifetime of carriers due to intrinsic scattering mechanisms viz. electron-electron interaction (EEI), electron-phonon interaction (EPI) and phonon-phonon interaction (PPI) in a promising half-Heusler thermoelectric FeVSb. Using the full-GWmethod, the effect of EEI and temperature on the valence and conduction band extrema and band gap are studied. The lifetime of carriers with temperature are estimated at these band extrema. At 300 K, estimated value of lifetime at VBM (CBM) is ∼1.91 × 10^-14 s (∼2.05 × 10^-14 s). The estimated ground state band gap considering EEI is ∼378 meV. Next, the effect of EPI on the lifetime of electrons and phonons with temperature are discussed. The comparison of two electron lifetimes suggests that EEI should be considered in transport calculations along with EPI. The average acoustic, optical and overall phonon lifetimes due to EPI are studied with temperature. Further, the effect of PPI is studied by computing average phonon lifetime for acoustic and optical phonon branches. The lifetime of the acoustic phonons are higher compared to optical phonons which indicates acoustic phonons contribute more to lattice thermal conductivity (κ_ph). The comparison of phonon lifetime due to EPI and PPI suggests that, above 500 K EPI is the dominant phonon scattering mechanism and cannot be ignored inκ_phcalculations. Lastly, a prediction of the power factor and figure of merit of n-type and p-type FeVSb is made by considering the temperature dependent carrier lifetime for the electronic transport terms. This study shows the importance of considering EEI in electronic transport calculations and EPI in phonon transport calculations in FeVSb. Our study is expected to provide results to further explore the thermoelectric transport in this material.

Investigating the effect of temperature dependent many-body interactions on electronic structures of SnTe in the Matsubara-time domain

Recently, SnTe has gained attention due to its non-trivial topological nature and eco-friendly thermoelectric applications. We report a detailed temperature dependent electronic structure of this compound using DFT andGWmethods. The calculated values of bandgaps by using PBEsol andG₀W₀methods are found to be in good agreement with the experiment, whereas mBJ underestimates the bandgap. The averaged value of diagonal matrix elements of fully screened Coulomb interaction (W̄) atω= 0 eV for Sn (Te) 5porbitals is ∼1.39 (∼1.70) eV. The nature of frequency dependentW̄(ω)reveals that the correlation strength of this compound is relatively weaker and hence the excited electronic state can be properly studied by full-GWmany-body technique. The plasmon excitation is found to be important in understanding this frequency dependentW̄(ω). The temperature dependent electron-electron interactions (EEI) reduces the bandgaps with increasing temperature. The value of bandgap at 300 K is obtained to be ∼161 meV. The temperature dependent lifetimes of electronic state alongW-L-Γ direction are also estimated. This work suggests that EEI is important to explain the high temperature transport behaviour of SnTe.

Finite-size effects in the reconstruction of dynamic properties from ab initio path integral Monte Carlo simulations

We systematically investigate finite-size effects in the dynamic structure factor S(q,ω) of the uniform electron gas obtained via the analytic continuation of ab initio path integral Monte Carlo data for the imaginary-time density-density correlation function F(q,τ). Using the recent scheme by Dornheim et al. [Phys. Rev. Lett. 121, 255001 (2018)PRLTAO0031-900710.1103/PhysRevLett.121.255001], we find that the reconstructed spectra are not afflicted with any finite-size effects for as few as N=14 electrons both at warm dense matter (WDM) conditions and at the margins of the strongly correlated electron liquid regime. Our results further corroborate the high quality of our current description of the dynamic density response of correlated electrons, which is of high importance for many applications in WDM theory and beyond.

Equation of motion coupled-cluster cumulant approach for intrinsic losses in x-ray spectra

We present a combined equation of motion coupled-cluster cumulant Green's function approach for calculating and understanding intrinsic inelastic losses in core level x-ray absorption spectra (XAS) and x-ray photoemission spectra. The method is based on a factorization of the transition amplitude in the time domain, which leads to a convolution of an effective one-body absorption spectrum and the core-hole spectral function. The spectral function characterizes intrinsic losses in terms of shake-up excitations and satellites using a cumulant representation of the core-hole Green's function that simplifies the interpretation. The one-body spectrum also includes orthogonality corrections that enhance the XAS at the edge.

Ab initio simulations and measurements of the free-free opacity in aluminum

The free-free opacity in dense systems is a property that both tests our fundamental understanding of correlated many-body systems, and is needed to understand the radiative properties of high energy-density plasmas. Despite its importance, predictive calculations of the free-free opacity remain challenging even in the condensed matter phase for simple metals. Here we show how the free-free opacity can be modelled at finite-temperatures via time-dependent density functional theory, and illustrate the importance of including local field corrections, core polarization, and self-energy corrections. Our calculations for ground-state Al are shown to agree well with experimental opacity measurements performed on the Artemis laser facility across a wide range of extreme ultraviolet wavelengths. We extend our calculations across the melt to the warm-dense matter regime, finding good agreement with advanced plasma models based on inverse bremsstrahlung at temperatures above 10 eV.

Ab initio Path Integral Monte Carlo Results for the Dynamic Structure Factor of Correlated Electrons: From the Electron Liquid to Warm Dense Matter

The accurate description of electrons at extreme density and temperature is of paramount importance for, e.g., the understanding of astrophysical objects and inertial confinement fusion. In this context, the dynamic structure factor S(q,ω) constitutes a key quantity as it is directly measured in x-ray Thomson scattering experiments and governs transport properties like the dynamic conductivity. In this work, we present the first ab initio results for S(q,ω) by carrying out extensive path integral Monte Carlo simulations and developing a new method for the required analytic continuation, which is based on the stochastic sampling of the dynamic local field correction G(q,ω). In addition, we find that the so-called static approximation constitutes a promising opportunity to obtain high-quality data for S(q,ω) over substantial parts of the warm dense matter regime.