Coherence Thermofield Dynamics Overview

• Coherence thermofield dynamics is a finite-temperature framework that represents mixed thermal states as pure states in a doubled Hilbert space, enabling pure-state propagation.
• It enables efficient evaluation of vibrationally resolved electronic spectra and offers fidelity and survival amplitude diagnostics in many-body quantum and chaotic systems.
• The approach leverages advanced propagation schemes, such as split-operator and semiclassical methods, to overcome computational challenges in low-dimensional settings.

Searching arXiv for the specified papers to ground the article in current records. Coherence thermofield dynamics is a finite-temperature formulation in which a mixed thermal state is represented as a pure state in a doubled Hilbert space, so that thermal correlation functions can be recast as wavepacket autocorrelation functions and propagated by a standard zero-temperature Schrödinger equation on an augmented configuration space. In vibronic spectroscopy this construction is used to evaluate vibrationally resolved electronic spectra at nonzero temperatures, while in many-body and quantum-chaotic settings the same thermofield structure supports fidelity- and survival-amplitude diagnostics of coherence, decoherence, and critical behavior (Zhang et al., 2023, Kröninger et al., 29 Jul 2025, Xu et al., 2020).

1. Doubled-space purification and thermofield states

The formal basis of coherence thermofield dynamics is the thermofield doubling of Hilbert space. For a Hamiltonian $H$ with eigenpairs $\{E_n,|n\rangle\}$ and partition function $Z(\beta)=\mathrm{Tr}\,e^{-\beta H}$, the thermofield double state is

$$|\mathrm{TFD}(0)\rangle=\frac{1}{\sqrt{Z(\beta)}}\sum_n e^{-\beta E_n/2}\,|n\rangle\otimes|n\rangle,$$

and, equivalently in ancilla notation,

$$|I\rangle=\sum_n |n\rangle\otimes|\tilde n\rangle,\qquad |\psi(\beta)\rangle=Z^{-1/2}\sum_n e^{-\beta E_n/2}|n\rangle\otimes|\tilde n\rangle .$$

Tracing out the ancillary or tilde factor reproduces the Gibbs state,

$$\rho(\beta)=e^{-\beta \hat H}/Z=\mathrm{Tr}_{\mathrm{ancilla}}\bigl[|\psi(\beta)\rangle\langle\psi(\beta)|\bigr],$$

so the thermofield state is a purification rather than a modification of the thermal ensemble (Zhang et al., 2023, Xu et al., 2020).

In coherence thermofield dynamics for spectroscopy, the initial doubled-space wavefunction is written in the coordinate basis as

$$\Psi_{\mathrm{tf}}(q,q';0)=\langle q|\sqrt{\hat\rho}|q'\rangle,$$

with doubled space $$\mathcal H_{\mathrm{TF}}=\mathcal H_q\otimes\mathcal H_{q'}$$. The associated thermofield state may also be expressed as

$$|\Psi_{\rm tf}(0)\rangle=(\hat\rho^{1/2}\otimes I)\,|I\rangle.$$

This construction converts finite-temperature density-operator data into a pure wavepacket in an augmented space of doubled dimension (Kröninger et al., 29 Jul 2025).

A broader geometric version of thermofield doubling appears in coherent-state and coadjoint-orbit formulations, where $$\mathcal H_{\rm TFD}=\mathcal H\otimes\tilde{\mathcal H}$$ and the thermofield vacuum $\{E_n,|n\rangle\}$0 is an entangled state satisfying $\{E_n,|n\rangle\}$1. In field-theoretic language this doubling can be organized on $\{E_n,|n\rangle\}$2, with the second copy assigned opposite orientation (Nair, 2015).

2. Correlation functions and coherence observables

The specific content of coherence thermofield dynamics is the rewriting of thermal coherence dynamics as pure-state propagation. In the Condon approximation, the finite-temperature dipole correlation is

$\{E_n,|n\rangle\}$3

In the doubled space, the coherence-form thermofield Hamiltonian is

$\{E_n,|n\rangle\}$4

or, in operator form for vibronic problems,

$\{E_n,|n\rangle\}$5

The thermal dipole autocorrelation then becomes a zero-temperature wavepacket autocorrelation: $\{E_n,|n\rangle\}$6 This avoids solving the von Neumann equation for the coherence and replaces it with a Schrödinger problem in the augmented space (Zhang et al., 2023, Kröninger et al., 29 Jul 2025).

In many-body thermofield dynamics, coherence is often monitored through survival amplitudes and fidelities. If one evolves one copy,

$\{E_n,|n\rangle\}$7

and

$\{E_n,|n\rangle\}$8

If one evolves both copies with

$\{E_n,|n\rangle\}$9

then

$Z(\beta)=\mathrm{Tr}\,e^{-\beta H}$0

The difference between $Z(\beta)=\mathrm{Tr}\,e^{-\beta H}$1 and $Z(\beta)=\mathrm{Tr}\,e^{-\beta H}$2 is therefore a matter of the evolution convention, not of the thermofield principle itself (Liu et al., 2024, Xu et al., 2020).

Under nonunitary evolution the overlap with the initial thermofield state remains the relevant coherence measure. If $Z(\beta)=\mathrm{Tr}\,e^{-\beta H}$3 and $Z(\beta)=\mathrm{Tr}\,e^{-\beta H}$4, then

$Z(\beta)=\mathrm{Tr}\,e^{-\beta H}$5

quantifies the remaining coherence under decoherence (Xu et al., 2020).

3. Propagation schemes for finite-temperature spectroscopy

A central use of coherence thermofield dynamics is the computation of finite-temperature vibronic spectra. In the split-operator implementation, the exact propagator is approximated by the second-order factorization

$Z(\beta)=\mathrm{Tr}\,e^{-\beta H}$6

with augmented potential

$Z(\beta)=\mathrm{Tr}\,e^{-\beta H}$7

and total kinetic energy defined on the doubled coordinate set. On a tensor-product grid, potential factors act by pointwise multiplication and kinetic factors by Fourier transformation in each coordinate. The method is numerically exact but the augmented grid doubles the number of coordinates, so the cost scales like $Z(\beta)=\mathrm{Tr}\,e^{-\beta H}$8 for $Z(\beta)=\mathrm{Tr}\,e^{-\beta H}$9 physical degrees of freedom and is limited to very low $$|\mathrm{TFD}(0)\rangle=\frac{1}{\sqrt{Z(\beta)}}\sum_n e^{-\beta E_n/2}\,|n\rangle\otimes|n\rangle,$$0, stated in practice as $$|\mathrm{TFD}(0)\rangle=\frac{1}{\sqrt{Z(\beta)}}\sum_n e^{-\beta E_n/2}\,|n\rangle\otimes|n\rangle,$$1–$$|\mathrm{TFD}(0)\rangle=\frac{1}{\sqrt{Z(\beta)}}\sum_n e^{-\beta E_n/2}\,|n\rangle\otimes|n\rangle,$$2 (Zhang et al., 2023).

The benchmark reported for a one-dimensional Morse test case uses a harmonic ground surface,

$$|\mathrm{TFD}(0)\rangle=\frac{1}{\sqrt{Z(\beta)}}\sum_n e^{-\beta E_n/2}\,|n\rangle\otimes|n\rangle,$$3

and a Morse excited surface,

$$|\mathrm{TFD}(0)\rangle=\frac{1}{\sqrt{Z(\beta)}}\sum_n e^{-\beta E_n/2}\,|n\rangle\otimes|n\rangle,$$4

At scaled temperatures $$|\mathrm{TFD}(0)\rangle=\frac{1}{\sqrt{Z(\beta)}}\sum_n e^{-\beta E_n/2}\,|n\rangle\otimes|n\rangle,$$5, the thermofield split-operator propagation agrees perfectly with Boltzmann averaging over individual Franck–Condon wavepacket propagations, including the appearance of hot bands at higher temperature (Zhang et al., 2023).

A semiclassical alternative combines coherence thermofield dynamics with the Herman–Kluk initial-value representation. With $$|\mathrm{TFD}(0)\rangle=\frac{1}{\sqrt{Z(\beta)}}\sum_n e^{-\beta E_n/2}\,|n\rangle\otimes|n\rangle,$$6 and $$|\mathrm{TFD}(0)\rangle=\frac{1}{\sqrt{Z(\beta)}}\sum_n e^{-\beta E_n/2}\,|n\rangle\otimes|n\rangle,$$7,

$$|\mathrm{TFD}(0)\rangle=\frac{1}{\sqrt{Z(\beta)}}\sum_n e^{-\beta E_n/2}\,|n\rangle\otimes|n\rangle,$$8

The Monte Carlo implementation samples initial conditions from the Husimi density and propagates two independent classical trajectories, one on the excited surface and one on the negative-ground surface. For Morse potentials of increasing anharmonicity evaluated at various temperatures, both Herman–Kluk coherence thermofield dynamics and single-trajectory thawed Gaussian coherence thermofield dynamics are accurate at low anharmonicity, but at higher anharmonicity the thawed Gaussian method fails to capture emerging hot bands, whereas the Herman–Kluk thermofield approach reproduces them. A direct comparison with a Boltzmann-averaged Herman–Kluk calculation shows perfect agreement, while the thermofield route obtains the result in a single simulation at roughly the cost of two zero-temperature Herman–Kluk runs (Kröninger et al., 29 Jul 2025).

The computational limitation of full-grid thermofield propagation motivates hybridization with techniques already developed for zero-temperature split-operator Fourier dynamics, including MCTDH, multilayer MCTDH, TT-SOFT, matching-pursuit coherent-state bases, adaptive or moving real-space grids, and high-order geometric compositions (Zhang et al., 2023).

4. Quantum chaos, decoherence, and characteristic times

In isolated many-body systems, thermofield coherence displays the same universal structures familiar from spectral statistics. For a thermofield state prepared from a single-copy Hamiltonian, the unitary fidelity $$|\mathrm{TFD}(0)\rangle=\frac{1}{\sqrt{Z(\beta)}}\sum_n e^{-\beta E_n/2}\,|n\rangle\otimes|n\rangle,$$9 starts at $$|I\rangle=\sum_n |n\rangle\otimes|\tilde n\rangle,\qquad |\psi(\beta)\rangle=Z^{-1/2}\sum_n e^{-\beta E_n/2}|n\rangle\otimes|\tilde n\rangle .$$0 and exhibits the characteristic decay, dip, ramp, and plateau of quantum-chaotic systems. In the stochastic Sachdev–Ye–Kitaev model, the four regimes are described by decay from unity to a first minimum at $$|I\rangle=\sum_n |n\rangle\otimes|\tilde n\rangle,\qquad |\psi(\beta)\rangle=Z^{-1/2}\sum_n e^{-\beta E_n/2}|n\rangle\otimes|\tilde n\rangle .$$1, a linear rise or ramp due to level-repulsion correlations, and saturation to the plateau

$$|I\rangle=\sum_n |n\rangle\otimes|\tilde n\rangle,\qquad |\psi(\beta)\rangle=Z^{-1/2}\sum_n e^{-\beta E_n/2}|n\rangle\otimes|\tilde n\rangle .$$2

at the Heisenberg time $$|I\rangle=\sum_n |n\rangle\otimes|\tilde n\rangle,\qquad |\psi(\beta)\rangle=Z^{-1/2}\sum_n e^{-\beta E_n/2}|n\rangle\otimes|\tilde n\rangle .$$3 (Xu et al., 2020).

Independent energy dephasing on each copy is modeled by the stochastic perturbation

$$|I\rangle=\sum_n |n\rangle\otimes|\tilde n\rangle,\qquad |\psi(\beta)\rangle=Z^{-1/2}\sum_n e^{-\beta E_n/2}|n\rangle\otimes|\tilde n\rangle .$$4

with real Gaussian white noise $$|I\rangle=\sum_n |n\rangle\otimes|\tilde n\rangle,\qquad |\psi(\beta)\rangle=Z^{-1/2}\sum_n e^{-\beta E_n/2}|n\rangle\otimes|\tilde n\rangle .$$5. Noise averaging yields a Lindblad equation for the doubled density matrix,

$$|I\rangle=\sum_n |n\rangle\otimes|\tilde n\rangle,\qquad |\psi(\beta)\rangle=Z^{-1/2}\sum_n e^{-\beta E_n/2}|n\rangle\otimes|\tilde n\rangle .$$6

In this setting decoherence introduces an information-loss timescale $$|I\rangle=\sum_n |n\rangle\otimes|\tilde n\rangle,\qquad |\psi(\beta)\rangle=Z^{-1/2}\sum_n e^{-\beta E_n/2}|n\rangle\otimes|\tilde n\rangle .$$7 through the short-time decay

$$|I\rangle=\sum_n |n\rangle\otimes|\tilde n\rangle,\qquad |\psi(\beta)\rangle=Z^{-1/2}\sum_n e^{-\beta E_n/2}|n\rangle\otimes|\tilde n\rangle .$$8

For the SYK model with large $$|I\rangle=\sum_n |n\rangle\otimes|\tilde n\rangle,\qquad |\psi(\beta)\rangle=Z^{-1/2}\sum_n e^{-\beta E_n/2}|n\rangle\otimes|\tilde n\rangle .$$9 and approximate Gaussian density of states,

$$\rho(\beta)=e^{-\beta \hat H}/Z=\mathrm{Tr}_{\mathrm{ancilla}}\bigl[|\psi(\beta)\rangle\langle\psi(\beta)|\bigr],$$0

with $$\rho(\beta)=e^{-\beta \hat H}/Z=\mathrm{Tr}_{\mathrm{ancilla}}\bigl[|\psi(\beta)\rangle\langle\psi(\beta)|\bigr],$$1 (Xu et al., 2020).

The effect of decoherence can be expressed as a temporal coarse-graining of the spectral form factor: $$\rho(\beta)=e^{-\beta \hat H}/Z=\mathrm{Tr}_{\mathrm{ancilla}}\bigl[|\psi(\beta)\rangle\langle\psi(\beta)|\bigr],$$2 As $$\rho(\beta)=e^{-\beta \hat H}/Z=\mathrm{Tr}_{\mathrm{ancilla}}\bigl[|\psi(\beta)\rangle\langle\psi(\beta)|\bigr],$$3 increases, the dip becomes shallower, its minimum shifts to later times, the pre-dip quantum noise is suppressed, and the span of the linear ramp is shortened. If $$\rho(\beta)=e^{-\beta \hat H}/Z=\mathrm{Tr}_{\mathrm{ancilla}}\bigl[|\psi(\beta)\rangle\langle\psi(\beta)|\bigr],$$4, quantum-noise fluctuations around the plateau are quenched while the plateau height $$\rho(\beta)=e^{-\beta \hat H}/Z=\mathrm{Tr}_{\mathrm{ancilla}}\bigl[|\psi(\beta)\rangle\langle\psi(\beta)|\bigr],$$5 remains unchanged. If $$\rho(\beta)=e^{-\beta \hat H}/Z=\mathrm{Tr}_{\mathrm{ancilla}}\bigl[|\psi(\beta)\rangle\langle\psi(\beta)|\bigr],$$6, the fidelity decays monotonically to the plateau with no visible dip or ramp (Xu et al., 2020).

A common misconception is that thermofield purification removes decoherence from the physical subsystem. The qubit example in thermofield dynamics shows the opposite distinction clearly: the full thermofield state remains pure in the doubled space, but after tracing out the tilde degrees of freedom the physical qubit is diagonal in the energy basis and therefore fully decohered in that basis. In that construction the loss of physical coherence is exactly encoded as entanglement with the tilde copy (Petronilo et al., 2021).

5. Fisher zeros, criticality, and thermofield coherence

Analytic continuation of the inverse temperature, $$\rho(\beta)=e^{-\beta \hat H}/Z=\mathrm{Tr}_{\mathrm{ancilla}}\bigl[|\psi(\beta)\rangle\langle\psi(\beta)|\bigr],$$7, links thermofield coherence to the complex partition function $$\rho(\beta)=e^{-\beta \hat H}/Z=\mathrm{Tr}_{\mathrm{ancilla}}\bigl[|\psi(\beta)\rangle\langle\psi(\beta)|\bigr],$$8. In the one-dimensional transverse-field Ising model,

$$\rho(\beta)=e^{-\beta \hat H}/Z=\mathrm{Tr}_{\mathrm{ancilla}}\bigl[|\psi(\beta)\rangle\langle\psi(\beta)|\bigr],$$9

the Fisher zeros are the points $$\Psi_{\mathrm{tf}}(q,q';0)=\langle q|\sqrt{\hat\rho}|q'\rangle,$$0 in the complex-$$\Psi_{\mathrm{tf}}(q,q';0)=\langle q|\sqrt{\hat\rho}|q'\rangle,$$1 plane at which $$\Psi_{\mathrm{tf}}(q,q';0)=\langle q|\sqrt{\hat\rho}|q'\rangle,$$2. In the thermodynamic limit these zeros condense onto smooth curves rather than remaining isolated points (Liu et al., 2024).

Two patterns are identified. For $$\Psi_{\mathrm{tf}}(q,q';0)=\langle q|\sqrt{\hat\rho}|q'\rangle,$$3, there are open lines parallel to the $$\Psi_{\mathrm{tf}}(q,q';0)=\langle q|\sqrt{\hat\rho}|q'\rangle,$$4 axis, which drift to $$\Psi_{\mathrm{tf}}(q,q';0)=\langle q|\sqrt{\hat\rho}|q'\rangle,$$5 as $$\Psi_{\mathrm{tf}}(q,q';0)=\langle q|\sqrt{\hat\rho}|q'\rangle,$$6. For all $$\Psi_{\mathrm{tf}}(q,q';0)=\langle q|\sqrt{\hat\rho}|q'\rangle,$$7, there are closed ovals encircling the imaginary-$$\Psi_{\mathrm{tf}}(q,q';0)=\langle q|\sqrt{\hat\rho}|q'\rangle,$$8 axis, which shrink toward $$\Psi_{\mathrm{tf}}(q,q';0)=\langle q|\sqrt{\hat\rho}|q'\rangle,$$9 as $$\mathcal H_{\mathrm{TF}}=\mathcal H_q\otimes\mathcal H_{q'}$$0 increases. Exactly at $$\mathcal H_{\mathrm{TF}}=\mathcal H_q\otimes\mathcal H_{q'}$$1, the topology changes: the open curves disappear and only the closed ovals remain. The paper identifies this qualitative change in Fisher-line topology as the signal of the quantum critical point (Liu et al., 2024).

The same complexified partition function governs thermofield survival dynamics: $$\mathcal H_{\mathrm{TF}}=\mathcal H_q\otimes\mathcal H_{q'}$$2 The short-time regime is Gaussian,

$$\mathcal H_{\mathrm{TF}}=\mathcal H_q\otimes\mathcal H_{q'}$$3

so the initial coherence decay is set by thermal fluctuations through the specific heat. At low temperature, the long-time behavior contains oscillations with the gap $$\mathcal H_{\mathrm{TF}}=\mathcal H_q\otimes\mathcal H_{q'}$$4, and for the one-dimensional transverse-field Ising model this reduces to

$$\mathcal H_{\mathrm{TF}}=\mathcal H_q\otimes\mathcal H_{q'}$$5

At the critical point $$\mathcal H_{\mathrm{TF}}=\mathcal H_q\otimes\mathcal H_{q'}$$6, the spectrum becomes equally spaced for large $$\mathcal H_{\mathrm{TF}}=\mathcal H_q\otimes\mathcal H_{q'}$$7, producing a revival at

$$\mathcal H_{\mathrm{TF}}=\mathcal H_q\otimes\mathcal H_{q'}$$8

The reported self-similarity

$$\mathcal H_{\mathrm{TF}}=\mathcal H_q\otimes\mathcal H_{q'}$$9

for $$|\Psi_{\rm tf}(0)\rangle=(\hat\rho^{1/2}\otimes I)\,|I\rangle.$$0 is presented as a manifestation of critical scale invariance (Liu et al., 2024).

This suggests a useful synthesis: in thermofield language, Fisher zeros are not merely a thermodynamic diagnostic but also a map of real-time coherence structures, because nonanalytic features of the Loschmidt echo occur when the contour $$|\Psi_{\rm tf}(0)\rangle=(\hat\rho^{1/2}\otimes I)\,|I\rangle.$$1 crosses the Fisher-zero manifold. The same work notes that $$|\Psi_{\rm tf}(0)\rangle=(\hat\rho^{1/2}\otimes I)\,|I\rangle.$$2 can be realized and probed in monitored quantum circuits (Liu et al., 2024).

6. Quantum simulation and neighboring formalisms

Gate-based preparation of thermofield states gives coherence thermofield dynamics a direct quantum-computing realization. For a single spin-$$|\Psi_{\rm tf}(0)\rangle=(\hat\rho^{1/2}\otimes I)\,|I\rangle.$$3, the finite-temperature thermofield vacuum is

$$|\Psi_{\rm tf}(0)\rangle=(\hat\rho^{1/2}\otimes I)\,|I\rangle.$$4

with reduced physical density matrix

$$|\Psi_{\rm tf}(0)\rangle=(\hat\rho^{1/2}\otimes I)\,|I\rangle.$$5

The magnetization is therefore

$$|\Psi_{\rm tf}(0)\rangle=(\hat\rho^{1/2}\otimes I)\,|I\rangle.$$6

A gate decomposition uses $$|\Psi_{\rm tf}(0)\rangle=(\hat\rho^{1/2}\otimes I)\,|I\rangle.$$7, nearest-neighbor CNOT gates, and a controlled phase, and the reported algorithm has circuit depth linear in system size (Petronilo et al., 1 Apr 2026).

The same thermofield logic appears in didactic qubit constructions. In a two-qubit realization of the thermal vacuum,

$$|\Psi_{\rm tf}(0)\rangle=(\hat\rho^{1/2}\otimes I)\,|I\rangle.$$8

the full doubled state is pure, with off-diagonal amplitude

$$|\Psi_{\rm tf}(0)\rangle=(\hat\rho^{1/2}\otimes I)\,|I\rangle.$$9

while the reduced physical qubit has Bloch vector $$\mathcal H_{\rm TFD}=\mathcal H\otimes\tilde{\mathcal H}$$0 and purity

$$\mathcal H_{\rm TFD}=\mathcal H\otimes\tilde{\mathcal H}$$1

This sharply separates coherence in the enlarged Hilbert space from coherence accessible on the physical subsystem alone (Petronilo et al., 2021).

The term “coherence thermofield dynamics” should also be distinguished from other TFD constructions that use the word “coherent.” Thermal coherent states defined via the Lie–Trotter product formula treat the thermalizing operator and the displacement operator symmetrically and are equivalent to more conventional thermal coherent states up to parameterization and phase; they were analyzed through uncertainty relations, quasiprobability distributions, and an optical parametric oscillator realization (Azuma et al., 2013). Generalized coherent states for deformed bosons in TFD extend this logic to Barut–Girardello and Klauder–Perelomov families using the Diagonal Operator Ordering Technique, with explicit links to $$\mathcal H_{\rm TFD}=\mathcal H\otimes\tilde{\mathcal H}$$2, internal energy, and free energy (Popov, 22 Dec 2025). Coherent-state path-integral TFD, finally, reformulates thermofield dynamics in terms of coadjoint-orbit actions and fields on $$\mathcal H_{\rm TFD}=\mathcal H\otimes\tilde{\mathcal H}$$3; in $$\mathcal H_{\rm TFD}=\mathcal H\otimes\tilde{\mathcal H}$$4 dimensions the opposite orientations of the two components lead, in the commutative limit, to the Einstein–Hilbert action as the difference of two Chern–Simons actions (Nair, 2015).

Taken together, these strands delimit the scope of coherence thermofield dynamics in the strict sense. In spectroscopy it denotes the rewriting of finite-temperature coherence propagation as pure-state dynamics in doubled space. In many-body dynamics it provides fidelity- and partition-function-based probes of chaos, decoherence, and criticality. In quantum simulation it offers a unitary encoding of thermal states on doubled registers. The unifying statement across these settings is that thermal mixed-state structure is transferred to entanglement structure in an enlarged Hilbert space, while the observable consequences depend on whether one retains or discards access to the auxiliary copy.