🔢 mathematics

Density functionals for the Pauli-Fierz Hamiltonian

This paper establishes a rigorous density-functional theory for the Pauli-Fierz Hamiltonian by proving that the mixed-state constrained-search functional coincides with the convex Lieb functional and that the pure-state Levy-Lieb functional is also convex due to photon purification, thereby unifying current-density-functional theory with quantum electrodynamics in both the zero-coupling and long-wavelength dipole limits.

Original author: Mihály A. Csirik, Andre Laestadius, Markus Penz, Erik I. Tellgren

Published 2026-11-03

📖 5 min read🧠 Deep dive

CC BY 4.0

Original author: Mihály A. Csirik, Andre Laestadius, Markus Penz, Erik I. Tellgren

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the author. For technical accuracy, refer to the original paper. Read full disclaimer

In the microscopic world where atoms and molecules live, matter is never truly alone. Even in a perfect vacuum, particles are constantly bathed in a sea of invisible energy: the electromagnetic field. While we often think of light as something we can see or feel as heat, at the quantum level, this field is made of discrete packets called photons. When electrons move through this field, they exchange these photons, creating a complex dance of attraction and repulsion that determines how matter behaves. For decades, scientists have used a powerful tool called density-functional theory to predict the properties of materials without having to track every single electron. This method simplifies the problem by focusing on the density of electrons—essentially, how crowded they are in different places—rather than their individual paths. However, this standard approach breaks down when the interaction with the electromagnetic field becomes strong, such as inside optical cavities where light is trapped and amplified. In these environments, the electrons and the light become so intertwined that they can no longer be treated as separate entities.

A team of researchers has now built a new mathematical framework to handle this deep connection between matter and light. They focused on a specific model known as the Pauli–Fierz Hamiltonian, which describes a system of electrons coupled to a full, continuous stream of massless photons. Unlike previous models that simplified the light into just a few modes or treated it as a background force, this work treats the light as a dynamic, quantized field that can contain any number of photons, including those with extremely low energy. The researchers discovered that to describe such a system accurately, one cannot rely on electron density alone. Instead, they identified a new, combined quantity that acts as the fundamental variable. This quantity is a sum of two parts: the usual flow of electrons and a second flow generated directly by the interaction with the quantized light field. By tracking this combined current, the team was able to map out exactly which electron and light configurations are physically possible.

One of the most striking findings of this work is how the light field acts as a purifier for the electrons. In quantum mechanics, a system can exist in a "mixed state," which is essentially a statistical blend of different possibilities, or a "pure state," which is a single, definite configuration. Usually, turning a mixed state into a pure one requires a significant amount of energy. However, the researchers proved that in this specific quantum electrodynamics setting, the presence of arbitrarily soft photons—light particles with vanishingly small energy—allows any mixed state of electrons to be represented as a pure state of the combined electron-light system. This happens with almost no energy cost. This result is profound because it means that the mathematical function used to calculate the energy of the system is perfectly smooth and convex, a property that makes the theory much more robust and easier to work with than previous electronic theories. It also implies that the distinction between mixed and pure states, which is often a major hurdle in quantum calculations, effectively disappears when the full light field is taken into account.

The study also explored what happens when the connection between the electrons and the light is turned off or when the light waves become incredibly long compared to the size of the system. In the limit where the light is turned off, the new theory smoothly transitions back to the established methods used for electrons in magnetic fields, confirming that the new framework is consistent with what is already known. When the light waves are stretched out, the theory converges to a well-known approximation used in practical chemistry, known as the dipole limit. This shows that the new, more rigorous approach does not contradict existing practical tools but rather provides a solid mathematical foundation for them. The researchers also demonstrated that while the theory guarantees the existence of a lowest-energy state for the combined system, finding a single, perfect pure state that minimizes the energy is not always possible in the strictest mathematical sense. Instead, the lowest energy is approached by states that contain a few very soft photons, which vanish in the limit.

By establishing these rigorous connections, the paper provides a complete and mathematically sound version of density-functional theory for systems where light and matter are inseparable. This work resolves long-standing questions about whether such a theory could exist and what its basic variables should be. It confirms that the combined current density is the correct key to unlocking the behavior of these systems. The results suggest that for realistic systems involving strong light-matter coupling, such as molecules inside optical cavities, the behavior of the material can be predicted with high confidence using these new functionals. The researchers did not just propose a new idea; they proved that the mathematical objects they defined behave exactly as required for a physical theory, ensuring that the energy calculations are stable and that the predicted states are physically realizable. This lays the groundwork for future simulations of chemical reactions and material properties in environments where light plays an active, quantum mechanical role.

Technical Summary: Density Functionals for the Pauli–Fierz Hamiltonian

Problem Statement
This work addresses the mathematical formulation of Density-Functional Theory (DFT) for nonrelativistic quantum electrodynamics (QED), specifically for systems described by the Pauli–Fierz Hamiltonian. While the Pauli–Fierz Hamiltonian itself is well-understood mathematically, the density-functional side has previously been analyzed only for simplified quantum-optical models (e.g., Rabi and Dicke models) where the system is reduced to two-level systems coupled to a finite number of field modes. The challenge lies in treating realistic systems involving $N$ electrons coupled to the full continuum of massless transverse photons, subject to external scalar and vector potentials, without an infrared cutoff. A central difficulty is that finite-energy states may contain arbitrarily many soft (low-energy) photons, complicating the analysis of representability and the existence of minimizers.

Methodology
The authors adopt the convex-analytic framework established by Lieb for standard DFT, extending it to the Pauli–Fierz setting. The methodology involves:

Variable Selection: The theory is built upon the standard electronic density $\rho$ and a combined current density $j_c$ . This combined current is defined as the sum of the paramagnetic current density ( $j_p$ ) and a field-induced current density ( $j_{ph}$ ) generated by the quantized vector potential. This choice ensures that the external potentials couple to the state only through these two variables.
Functional Definitions: Two universal functionals are defined:
- The Pauli–Fierz Levy–Lieb (PFLL) functional, defined via a pure-state constrained search over normalized wavefunctions in the form domain of the internal Hamiltonian.
- The Pauli–Fierz Lieb (PFL) functional, defined as the convex conjugate of the ground-state energy (reparametrized to absorb the diamagnetic term), which corresponds to a mixed-state constrained search.
Analytical Tools: The proofs rely on rigorous functional analysis, including:
- Photon Purification: Utilizing the infinite-dimensional nature of the photon Fock space to purify mixed electronic states.
- Soft-Photon Removal: Techniques to remove arbitrarily soft photons from minimizing sequences at negligible energy cost while preserving the density and combined current density.
- Localization and Compactness: Using geometric localization in Fock space and partial traces to handle the lack of compactness caused by photon loss.
- Mosco Convergence: Analyzing limits of functionals in the sense of Mosco to establish convergence in zero-coupling and long-wavelength limits.

Key Contributions and Results

Representability: The paper proves that the set of pure-state $N$ -representable density–current pairs coincides exactly with the set of mixed-state representable pairs. Specifically, any finite-energy electronic mixed state can be purified by a pure electron–photon state with finite energy. This establishes that the pure-state and mixed-state $N$ -representability classes are identical and convex.
Equality of Functionals: A primary result is the equality $F^{e,\chi}_{PFLL} = F^{e,\chi}_{PFL}$ . Unlike in standard electronic DFT, where the pure-state Levy–Lieb functional is generally non-convex and its infimum may not be attained, the Pauli–Fierz Levy–Lieb functional is convex and lower semicontinuous. This is attributed to the ability of soft photons to purify mixed states at arbitrarily small energy cost without altering the density or combined current.
Existence and Non-Existence of Optimizers:
- The mixed-state constrained search functional attains its infimum at a density matrix with finite mean photon number.
- Conversely, the pure-state constrained search generally does not attain its infimum, even for a single electron. The authors construct explicit examples where the infimum is approached by a sequence of pure states but not reached by any single pure state.
Soft-Photon Purification: The authors demonstrate that arbitrarily soft photons can purify any finite-energy mixed state at an arbitrarily small energy cost. This mechanism is responsible for the convexity of the pure-state functional and the coincidence of pure and mixed representability.
Zero-Coupling Limit: As the coupling constant $e \to 0$ , the Pauli–Fierz functionals converge to the mixed-state Vignale–Rasolt functional of current-density-functional theory (CDFT). The convergence holds in the sense of Mosco and pointwise.
Long-Wavelength (Dipole) Limit: The paper establishes the convergence of the functionals to a dipole functional in the long-wavelength limit (where the system size is small compared to the photon wavelength). This limiting functional is identified with the constrained search for the standard length-gauge dipole Hamiltonian.

Significance
The paper provides the first rigorous mathematical foundation for QEDFT based on the full Pauli–Fierz Hamiltonian, moving beyond simplified model systems. Its significance lies in:

Clarifying the Mathematical Structure: It resolves the delicate issues of $N$ -representability and the relationship between pure and mixed states in the presence of a quantized radiation field.
Establishing Convexity: It shows that the inclusion of the quantized field fundamentally alters the convexity properties of the universal functional, rendering the pure-state functional convex and lower semicontinuous, a property not shared by standard electronic DFT.
Rigorous Limits: It rigorously connects the full QEDFT to established theories (CDFT) and practical approximations (dipole approximation) through well-defined mathematical limits.
Foundational for Approximations: By establishing the exact properties of the universal functionals, the work provides a necessary benchmark for developing and validating approximate functionals in quantum-electrodynamical DFT.

The authors note that while the pure-state infimum is not generally attained, the existence of mixed-state optimizers and the equality of the functionals ensure that the variational principle remains well-posed. The work also highlights that the "price" for the convexity of the pure-state functional is the potential non-attainment of the minimizer, a phenomenon driven by the infrared freedom of the photon field.

More like this