Nonlinear electromagnetic response of a uniform electron gas

Density Functional Theory Perspective on the Nonlinear Response of Correlated Electrons across Temperature Regimes

We explore a new formalism to study the nonlinear electronic density response based on Kohn-Sham density functional theory (KS-DFT) at partially and strongly quantum degenerate regimes. It is demonstrated that the KS-DFT calculations are able to accurately reproduce the available path integral Monte Carlo simulation results at temperatures relevant for warm dense matter research. The existing analytical results for the quadratic and cubic response functions are rigorously tested. It is demonstrated that the analytical results for the quadratic response function closely agree with the KS-DFT data. Furthermore, the performed analysis reveals that currently available analytical formulas for the cubic response function are not able to describe simulation results, neither qualitatively nor quantitatively, at small wavenumbers q < 2q_F, with q_F being the Fermi wavenumber. The results show that KS-DFT can be used to describe warm dense matter that is strongly perturbed by an external field with remarkable accuracy. Furthermore, it is demonstrated that KS-DFT constitutes a valuable tool to guide the development of the nonlinear response theory of correlated quantum electrons from ambient to extreme conditions. This opens up new avenues to study nonlinear effects in a gamut of different contexts at conditions that cannot be accessed with previously used path integral Monte Carlo methods.

Nonlinear density response from imaginary-time correlation functions: Ab initio path integral Monte Carlo simulations of the warm dense electron gas

The ab initio path integral Monte Carlo (PIMC) approach is one of the most successful methods in quantum many-body theory. A particular strength of this method is its straightforward access to imaginary-time correlation functions (ITCFs). For example, the well-known density-density ITCF F(q, τ) allows one to estimate the linear response of a given system for all wave vectors q from a single simulation of the unperturbed system. Moreover, it constitutes the basis for the reconstruction of the dynamic structure factor S(q, ω)-a key quantity in state-of-the-art scattering experiments. In this work, we present analogous relations between the nonlinear density response in the quadratic and cubic order of the perturbation strength and generalized ITCFs measuring correlations between up to four imaginary-time arguments. As a practical demonstration of our new approach, we carry out simulations of the warm dense electron gas and find excellent agreement with previous PIMC results that had been obtained with substantially larger computational effort. In addition, we give a relation between a cubic ITCF and the triple dynamic structure factor S(q₁, ω₁; q₂, ω₂), which evokes the enticing possibility to study dynamic three-body effects on an ab initio level.

Nonlinear Light Mixing by Graphene Plasmons

Graphene is known to possess strong optical nonlinearity which turned out to be suitable for creation of efficient saturable absorbers in mode locked fiber lasers. Nonlinear response of graphene can be further enhanced by the presence of graphene plasmons. Here, we report a novel nonlinear effect observed in nanostructured graphene which comes about due to excitation of graphene plasmons. We experimentally detect and theoretically explain enhanced mixing of near-infrared and mid-infrared light in arrays of graphene nanoribbons. Strong compression of light by graphene plasmons implies that the described effect of light mixing is nonlocal in nature and orders of magnitude larger than the conventional local graphene nonlinearity. Both second and third order nonlinear effects were observed in our experiments with the recalculated third-order nonlinearity coefficient reaching values of 4.5 × 10^-6 esu. The suggested effect could be used in variety of applications including nonlinear light modulators, light multiplexers, light logic, and sensing devices.

Frenetic aspects of second order response

A theoretical framework to compute the second order response around equilibrium is provided and the importance of the so called frenetic contribution is illustrated.