For certain correlated electron–photon systems we construct the exact density-to-potential maps, which are the basic ingredients of a density-functional reformulation of coupled matter-photon problems. We do so for numerically exactly solvable models consisting of up to four fermionic sites coupled to a single photon mode. We show that the recently introduced concept of the intra-system steepening (Dimitrov et al 2016 New J. Phys. 18 083004) can be generalized to coupled fermion-boson systems and that the intra-system steepening indicates strong exchange-correlation effects due to the coupling between electrons and photons. The reliability of the mean-field approximation to the electron–photon interaction is investigated and its failure in the strong coupling regime analyzed. We highlight how the intra-system steepening of the exact density-to-potential maps becomes apparent also in observables such as the photon number or the polarizability of the electronic subsystem. We finally show that a change in functional variables can make these observables behave more smoothly and exemplify that the density-to-potential maps can give us physical insights into the behavior of coupled electron–photon systems by identifying a very large polarizability due to ultra-strong electron–photon coupling.
For certain correlated electron–photon systems we construct the exact density-to-potential maps, which are the basic ingredients of a density-functional reformulation of coupled matter-photon problems. We do so for numerically exactly solvable models consisting of up to four fermionic sites coupled to a single photon mode. We show that the recently introduced concept of the intra-system steepening (Dimitrov et al 2016 New J. Phys. 18 083004) can be generalized to coupled fermion-boson systems and that the intra-system steepening indicates strong exchange-correlation effects due to the coupling between electrons and photons. The reliability of the mean-field approximation to the electron–photon interaction is investigated and its failure in the strong coupling regime analyzed. We highlight how the intra-system steepening of the exact density-to-potential maps becomes apparent also in observables such as the photon number or the polarizability of the electronic subsystem. We finally show that a change in functional variables can make these observables behave more smoothly and exemplify that the density-to-potential maps can give us physical insights into the behavior of coupled electron–photon systems by identifying a very large polarizability due to ultra-strong electron–photon coupling.
Providing sufficient accuracy at low-computational cost, not surprisingly, density-functional theory ( DFT ) is a standard tool in theoretical material science. Similarly, quantum-electrodynamical density-functional theory ( QEDFT ) provides a promising route to describe experiments at the interface of quantum optics and quantum chemistry. Formally, in the Kohn-Sham ( KS ) formulation of DFT and QEDFT the exact many-particle Schrödinger equation is replaced by a set of coupled non-linear one-particle differential equations. The electron-electron, the electron-ion and the electron-photon interactions of the interacting system are mimicked by their corresponding effective KS potential with an, in principle unknown exchange-correlation ( xc ) term. In practice, DFT and QEDFT require approximations for this xc potential. The quality of which fully determines how much the DFT and QEDFT observables deviate from the exact ones. Understanding the shortcomings and improving the quality of approximate functionals makes the study of the exact functional of model systems indispensable. Here, we employ exact diagonalization for the Hamiltonian of different types of real-space lattice models in one spatial dimension to access the ground-state wave function, and thus to explicitly construct the exact functionals of DFT and QEDFT . To study the functional behavior of purely electronic systems, we consider two singlet-state electrons moving in an external potential with two connected atomic wells of variable depth and distance. To study coupled electron-photon systems, we consider the Rabi Hamiltonian with variable electron-photon coupling strength. By varying the model parameters, we systematically construct the functionals for systems from weak to strong coupling between the building blocks of the global system. Using the concept of the effective density-to- potential map, we discuss the intrasystem steepening which is an exact condition on the xc functional that captures the static correlation present in the system and carve out how it relates to the well-known intersystem derivative discontinuity. We further construct arbitrary observables as functionals of the internal variables which are directly affected by such features of the effective density-to-potential map. By constructing the reduced dipole map in the four-particle case of coupled electron-photon system, we show that the density-to-potential map helps to identify and understand recently discovered physical effects in the field of strongly coupled electron-photon systems. This work provides a contribution to the development of approximate function- als required in practical density-functional calculations of arbitrary experimentally accessible observables. Both the intrasystem steepening and the intersystem derivative discontinuity play a crucial role for systems with highly-localized electrons. Including both features in approximate functionals will improve the accu- racy of DFT evaluation of observables for such systems. Those results provide novel insight into the properties of the exact xc potentials in QEDFT . This work opens a new route to construct new functionals using coarse-grained or embed- ding approaches within the general and exact frameworks of DFT and QEDFT.
In the present work, we employ exact diagonalization for model systems on a real-space lattice to explicitly construct the exact density-to-potential and graphically illustrate the complete exact density-to-wavefunction map that underly the Hohenberg–Kohn theorem in density functional theory. Having the explicit wavefunction-to-density map at hand, we are able to construct arbitrary observables as functionals of the ground-state density. We analyze the density-to-potential map as the distance between the fragments of a system increases and the correlation in the system grows. We observe a feature that gradually develops in the density-to-potential map as well as in the density-to-wavefunction map. This feature is inherited by arbitrary expectation values as functional of the ground-state density. We explicitly show the excited-state energies, the excited-state densities, and the correlation entropy as functionals of the ground-state density. All of them show this exact feature that sharpens as the coupling of the fragments decreases and the correlation grows. We denominate this feature as intra-system steepening and discuss how it relates to the well-known inter-system derivative discontinuity. The inter-system derivative discontinuity is an exact concept for coupled subsystems with degenerate ground state. However, the coupling between subsystems as in charge transfer processes can lift the degeneracy. An important conclusion is that for such systems with a near-degenerate ground state, the corresponding cut along the particle number N of the exact density functionals is differentiable with a well-defined gradient near integer particle number.
We present results of a photon-free exchange-correlation functional within the local density approximation (pxcLDA) for quantum electrodynamics density functional theory (QEDFT) that efficiently describes the electron density of many-electron systems across weak to strong light-matter coupling. Building on previous work [I-Te. Lu et al., Phys. Rev. A 109, 052823 (2024)] that captured electron-photon correlations via an exchange-correlation functional derived from the nonrelativistic Pauli-Fierz Hamiltonian and tested on one-electron systems, we use a simple procedure to compute a renormalization factor describing electron-photon correlations and inhomogeneity in the weak-coupling regime by comparing it with quantum electrodynamics coupled-cluster, and previous QEDFT optimized effective potential methods. Across various atoms and molecules, pxcLDA reproduces cavity-modified densities in close agreement with these references. The renormalization factor approaches unity as the system size or collective coupling increases, reflecting an electron-photon exchange-dominated behavior and improved accuracy for larger systems. This approach now offers a practical route to applying QEDFT functionals based on electron density to realistic electron systems.
The exact interaction energy of a many-electron system is determined by the electron pair density, which is not well-approximated in standard Kohn--Sham density functional models. Here we study the (complicated but well-defined) exact universal map from density to pair density. We survey how many common functionals, including the most basic version of the local density approximation (Dirac exchange with no correlation contribution), arise from particular approximations of this map. We develop an algorithm to compute the map numerically and apply it to one-parameter families {αρ(αx)}α>0 of one-dimensional homogeneous and inhomogeneous single-particle densities. We observe that the pair density develops remarkable multiscale patterns which strongly depend on both the particle number and the “width” α−1 of the single-particle density. The simulation results are confirmed by rigorous asymptotic results in the limiting regimes α≫ and α≪1. For one-dimensional homogeneous systems, we show that the whole spectrum of patterns is reproduced surprisingly well by a simple asymptotics-based ansatz which slowly smoothens out the “strictly correlated” α=0 pair density while slowly turning on the α=∞ “exchange” terms as α increases. Our findings lend theoretical support to the celebrated semi-empirical idea [A. D. Becke, J. Chem. Phys., 98 (1993), pp. 5648--5652] to mix in a fractional amount of exchange, albeit not to assuming the mixing to be additive and taking the fraction to be a system-independent constant.
The last two decades have witnessed increasing experimental interest in the study and control of many-electron systems strongly interacting with quantum electromagnetic fields. This includes notable experiments in the areas of cavity and circuit quantum electrodynamics, quantum computing via photon mediated atom entanglement, electromagnetically induced transparency, quantum plasmonics, quantum simulators and chemistry. The description of realistic coupled matter-photon systems requires combining electronic structure methods from material science with quantum-optical methods. In this work, we propose a formally exact and numerically feasible approach, termed quantum electrodynamical time-dependent density functional theory (QED-TDDFT), that generalizes successful TDDFT for electronic structure calculations, to the electron-photon coupling. We develop a general framework for the description of coupled matter-photon systems in all possible realizations of interest, ranging from the fully relativistic case, to well-known quantum-optical model (e.g., the Rabi model). Of particular interest for condensed-matter applications is the density functionalisation of non-relativistic QED. The underlying theory in fact corresponds to standard non-relativistic quantum mechanics for the electrons, but with additional terms, which correct for relativity and electron-photon interactions. However, any application of density-functional theory requires approximations to the exchange-correlation (xc) functional. In this work we construct such an approximation for the description of electron-photon interaction effects in an optical cavity (the theory of reference is QED-TDDFT for many-electron systems coupled to cavity photons). The derived electron-photon optimized effective potential (OEP) reduces in the static limit to the Lamb shift of the ground state. The OEP is then tested on the Rabi model from weak to strong coupling regimes. The results are promising and open a path for describing complex strongly coupled matter-photon systems. Further, we consider the inclusion of relativistic spin-spin effects into the xc functional of standard spin density functional theory. In this way we provide quantum corrections that are essential for the description of magnetic inhomogeneities at the nanoscale. Such an approach improves over the current (classical) micromagnetic treatment of domain walls and skyrmions for spintronic applications.