⚛️ quantum physics

Second-order mean force Gibbs state of fermionic open systems

This paper derives the second-order mean force Gibbs state for generic fermionic open systems coupled to a thermal reservoir, demonstrating how reservoir-mediated corrections to the bare Gibbs state capture particle-hole channel effects and induced coherence, as validated by numerical simulations across several model systems.

Original author: David J. Strachan

Published 2026-11-03

📖 5 min read🧠 Deep dive

CC BY 4.0

Original author: David J. Strachan

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

Thermodynamics is the science of heat, energy, and how things settle down when they are left alone. For over a century, scientists have relied on a simple rule: if you place a small object in contact with a large, warm environment, the object will eventually relax into a state of perfect balance. In this state, the object's behavior is determined entirely by the temperature of the environment and how many particles, like electrons, are flowing in and out. This balanced state is known as the Gibbs state. It works beautifully for large, everyday objects where the connection between the object and its surroundings is so weak that the two barely notice each other. However, the world of quantum mechanics deals with the very small, where this assumption often breaks down. In devices like quantum dots or molecular junctions, the link between the system and its environment is not a faint whisper but a strong handshake. When the connection is this tight, the old rules no longer apply, and the object does not settle into the expected state. Instead, it reaches a different kind of balance, one that is shaped by the very act of holding hands with its surroundings.

For decades, physicists have known that a more accurate description exists, called the mean force Gibbs state. This state is found by looking at the combined system of the object and its environment as a single unit, and then asking what the object looks like when the environment is ignored. While this concept has been well understood for systems involving light waves or vibrations, known as bosons, it has remained a mystery for systems made of particles like electrons, known as fermions. The mathematics for fermions is notoriously difficult because these particles obey a strict rule that prevents them from occupying the same space in the same way, making the standard tools of calculation fail. Without a clear formula, scientists could not predict how these tiny electronic systems would behave when strongly coupled to their environment, leaving a gap in our understanding of quantum thermodynamics.

In a recent study, David Strachan from the University of Bristol has filled this gap by deriving a precise formula for the equilibrium state of fermionic systems. The researcher focused on a generic setup where a system of electrons exchanges particles with a thermal reservoir, a large collection of electrons at a fixed temperature. By using a specialized mathematical technique called the mode reduced state, which properly handles the unique rules of fermions, Strachan calculated the state of the system up to the second order of the coupling strength. This means the formula accounts for the interaction between the system and the environment with high precision, capturing effects that simpler models miss. The result is a closed-form expression that describes exactly how the system settles down, separating the corrections into two distinct channels: one for particles entering the system and one for holes, or empty spaces, leaving it.

The study reveals that the equilibrium state of these strongly coupled systems is richer and more complex than previously thought. Unlike the simple Gibbs state, which is completely silent and static, the new mean force Gibbs state retains a kind of internal memory. It preserves quantum coherence, which is a delicate relationship between different energy levels of the system. In the traditional view, these relationships would fade away as the system interacts with the environment, but this new formula shows they can persist, especially at low temperatures. This is not just a theoretical curiosity; it suggests that the environment itself can act as a bridge, creating connections between parts of the system that do not talk to each other directly. This phenomenon was demonstrated in a model where two independent electron modes shared a single reservoir. Even though the two modes did not interact with each other, the reservoir mediated a link between them, creating a coherence that the old, simple models would have completely missed.

To ensure the formula was correct, the researcher tested it against rigorous computer simulations. For systems that do not interact with themselves, the team compared the new formula to exact calculations of how the system evolves over time, finding perfect agreement. For a more complex, interacting system known as the single impurity Anderson model, which mimics a magnetic impurity in a metal, the team used advanced numerical methods to simulate the real-time dynamics of the system relaxing into equilibrium. The results from these simulations matched the static prediction of the new formula in the regime of weak to intermediate coupling. This agreement provides strong evidence that the mean force Gibbs state is indeed the correct description of equilibrium for these fermionic systems. The work confirms that when the connection between a quantum system and its environment is strong, the environment leaves a permanent mark on the system's state, shaping its energy and its internal relationships in ways that were previously invisible to science.

This discovery matters because it provides the necessary reference point for defining thermodynamic quantities in the strong coupling regime. Without a correct understanding of what the equilibrium state looks like, it is impossible to accurately measure how much work a quantum battery can store or how efficiently a nanoscale heat engine can operate. By establishing the correct baseline for fermionic systems, this research opens the door to designing and understanding the next generation of quantum technologies, from molecular electronics to thermal machines, where the boundary between the system and its surroundings is no longer a thin line but a deep, entangled connection.

Technical Summary: Second-order mean force Gibbs state of fermionic open systems

Problem Statement
Standard thermodynamics assumes that a system in contact with a thermal reservoir relaxes to a bare Gibbs state (or grand canonical Gibbs state), defined solely by the system's Hamiltonian ( $H_S$ ), particle number ( $N_S$ ), and the reservoir's temperature ( $T$ ) and chemical potential ( $\mu$ ). However, this assumption holds only in the weak coupling limit. In microscopic regimes, such as quantum dots and molecular junctions, the system-reservoir coupling is often comparable to the system's energy scales and temperature. In these strong coupling regimes, the true equilibrium state is the Mean Force Gibbs (MFG) state, obtained by tracing the reservoir out of the total grand canonical Gibbs state of the combined system.

While analytic results for the MFG state exist for bosonic reservoirs (e.g., Cresser and Anders), a general analytic formulation for fermionic systems has been lacking. A primary obstacle is the definition of the MFG state for fermions: the standard partial trace operation relies on a tensor product structure of the Hilbert space, which fermionic systems do not possess due to their canonical anti-commutation relations (CAR). Consequently, the definition of the reduced state requires the specific concept of the "mode reduced state." Furthermore, the interaction operator for fermionic open systems typically involves particle exchange ( $V = X^\dagger E + E^\dagger X$ ), differing from the Hermitian factorized couplings often considered in bosonic treatments.

Methodology
The author derives a closed-form expression for the grand canonical MFG state ( $\tau_{MF}$ ) for a generic fermionic system coupled to a number-conserving fermionic reservoir, accurate to second order in the coupling strength ( $\lambda$ ).

Formalism: The system is defined by $H = H_S + H_R + \lambda V$ , where the interaction $V$ transfers particles between the system and reservoir. The author utilizes the fermionic partial trace defined via the mode reduced state (following Strachan et al. [34, 35]), ensuring the reduced state is well-defined by fixing a specific mode ordering (system modes preceding reservoir modes) and assuming parity symmetry.
Derivation: The derivation employs the Kubo expansion of the total Gibbs state $e^{-\beta(H-\mu N)}$ . The coupling operator $X$ is decomposed into eigenoperators $X_\omega$ of the system's grand canonical generator $K_S = H_S - \mu N_S$ .
Coefficient Functions: The correction terms are expressed in terms of four coefficient functions ( $A_1, A_2, D_1, D_2$ ). These are derived from the imaginary-time transforms of two reservoir correlation functions ( $G_1$ and $G_2$ ). Crucially, the Kubo-Martin-Schwinger (KMS) condition links these correlation functions, allowing all coefficients to be fixed by a single reservoir function and its derivative.
Hamiltonian of Mean Force: The author also derives the corresponding Hamiltonian of mean force ( $H_{MF}$ ), such that $\tau_{MF} \propto e^{-\beta(H_{MF} - \mu N_S)}$ . They note that while the exponential form is useful, it is strictly valid only when population corrections are small compared to bare populations (typically high temperatures), whereas the direct expansion of the state (Eq. 11) remains valid at lower temperatures.

Key Results
The main result is the second-order MFG state $\tau^{(2)}_{MF} = \tau_S + \lambda^2(\Delta_{diag} + \Delta_{off})$ , where:

Diagonal Corrections ( $\Delta_{diag}$ ): These modify the populations of the system's energy eigenstates. The correction separates into "particle" and "hole" channels, weighted by the empty and occupied reservoir states, respectively.
Off-Diagonal Corrections ( $\Delta_{off}$ ): These introduce coherences between different energy sectors of the bare system Hamiltonian $H_S$ . Unlike the bare Gibbs state, the MFG state does not necessarily commute with $H_S$ .
Low and High Temperature Limits: At high temperatures, the correction vanishes, recovering the bare Gibbs state. At low temperatures, the coherences persist, and the state retains off-diagonal elements in the bare energy eigenbasis.

Applications and Verification
The author applies their formula to three specific models and verify the results against numerically exact calculations:

Resonant Level Model (RLM): A single mode coupled to a reservoir. The MFG prediction captures the level shift and broadening effects, matching exact results at weak to intermediate coupling.
Two Modes Coupled to a Shared Reservoir: Two non-interacting modes coupled to the same reservoir. The MFG state predicts the emergence of reservoir-mediated coherence between the two modes. This coherence persists even as $T \to 0$ . The author verifies this against exact correlation-matrix dynamics.
Single Impurity Anderson Model (SIAM): An interacting model with spin. Since exact analytic solutions are unavailable for interacting systems, the author compares the MFG prediction against real-time dynamics computed using Matrix Product States (MPS) and the thermofield chain mapping. The static MFG prediction agrees with the steady state reached by the dynamical evolution in the weak to intermediate coupling regime.

Significance and Claims
The paper claims to close a significant gap in the literature by providing the first general analytic expression for the MFG state of generic fermionic open systems. The work extends the field of strong-coupling thermodynamics beyond bosonic reservoirs.

The author emphasizes that their results provide a necessary reference equilibrium state for defining thermodynamic quantities (such as free energy and work) for fermionic systems in the strong coupling regime, which is relevant for quantum dots, molecular junctions, and nanoscale heat engines. They explicitly state that while they provide strong numerical evidence that the MFG state is the relevant long-time equilibrium state for these systems, a general proof that the MFG state is the unique steady state for all open systems remains an open question. The paper does not propose new experimental setups but rather offers a theoretical framework to interpret existing and future systems where system-reservoir coupling is non-negligible.

More like this