Abstract

In a recent Letter, Dornheim et al. [PRL 125, 085001 (2020)] have investigated the nonlinear density response of the uniform electron gas in the warm dense matter regime. More specifically, they have studied the cubic response function at the first harmonic, which cannot be neglected in many situations of experimental relevance. In this work, we go one step further and study the full spectrum of excitations at the higher harmonics of the original perturbation based on extensive new ab initio path integral Monte Carlo (PIMC) simulations. We find that the dominant contribution to the density response beyond linear response theory is given by the quadratic response function at the second harmonic in the moderately nonlinear regime. Furthermore, we show that the nonlinear density response is highly sensitive to exchange-correlation effects, which makes it a potentially valuable new tool of diagnostics. To this end, we present a new theoretical description of the nonlinear electronic density response based on the recent effective static approximation to the local field correction [PRL 125, 235001 (2020)], which accurately reproduces our PIMC data with negligible computational cost.

Highlights

  • Linear response theory (LRT), i.e., the assumption of a linear response of a system of interest to a sufficiently small perturbation, is ubiquitous throughout physics and related disciplines

  • We have obtained extensive path integral Monte Carlo (PIMC) results carrying out simulations of a harmonically perturbed electron gas, which has allowed us to obtain the full spectrum of excitations at the integer harmonics of the original perturbation

  • The second potentially important nonlinear term is given by the cubic response at the first harmonic studied in Ref. [56], whereas the cubic response at the third harmonic is practically negligible in most realistic situations

Read more

Summary

INTRODUCTION

Linear response theory (LRT), i.e., the assumption of a linear response of a system of interest to a sufficiently small perturbation, is ubiquitous throughout physics and related disciplines. In a recent letter [56], three of us carried out extensive ab initio path integral Monte Carlo (PIMC) simulations of the harmonically perturbed electron gas in the WDM regime without any assumptions about the response being linear First and foremost, this has allowed us to check the validity range of LRT in WDM systems. We find that the nonlinear response functions strongly depend on electronic exchange-correlation effects, making them a potentially valuable new tool of diagnostics To this end, we introduce a theory of the nonlinear density response based on available representations of the LFC [58,59], which is capable to accurately reproduce our PIMC data with negligible computational cost.

THEORY
Path integral Monte Carlo
Simulation idea
Theory of nonlinear density response
Noninteracting case
Random phase approximation
Going beyond the RPA description by using LFC
PIMC approach to nonlinear density response functions
RESULTS
Temperature dependence
SUMMARY AND DISCUSSION
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call