---
title: Thermal Pseudo Entropy & Analytic Continuation
url: https://www.emergentmind.com/topics/thermal-pseudo-entropy
type: topic
---

# Thermal Pseudo Entropy & Analytic Continuation

Searching arXiv for recent papers on thermal pseudo-entropy and closely related pseudo-entropy work.
Thermal pseudo-entropy is a complex generalization of thermal entropy obtained by analytically continuing the inverse temperature to a complex value and, equivalently, by computing the pseudo-entropy of the transition matrix between thermofield-double states at different times. In the formulation of Caputa et al., the quantity is defined as the von Neumann entropy of a reduced transition matrix that is thermal at complex inverse temperature $\beta+it$, so that ordinary thermal entropy is recovered at $t=0$ [2411.08948]. Closely related constructions appear in earlier subregion pseudo-entropy for time-evolved thermofield-double states in two-dimensional CFTs [2109.00372], and in later work showing that thermal pseudo-entropy is a special case of real-time pseudo-entropy governed at short times by correlations between the physical Hamiltonian and the modular Hamiltonian [2606.14208].

## 1. Definition through thermofield-double transitions

Let a Hamiltonian $H$ have spectrum $\{E_n\}$ and partition function
$$
Z(\beta)=\sum_n e^{-\beta E_n}.
$$
The thermofield-double state on two copies, labeled $L$ and $R$, is
$$
|\Psi_\beta\rangle=\frac{1}{\sqrt{Z(\beta)}}\sum_n e^{-\frac{\beta}{2}E_n}\,|E_n\rangle_L\otimes|E_n\rangle_R.
$$
Evolving only the left copy for real time $t$ gives
$$
|\Psi_\beta(t)\rangle=e^{-iH_L t}|\Psi_\beta\rangle
=\frac{1}{\sqrt{Z(\beta)}}\sum_n e^{-\tfrac{\beta+2it}{2}E_n}\,|E_n\rangle_L\otimes|E_n\rangle_R.
$$
The normalized transition matrix between $|\Psi_\beta\rangle$ and $|\Psi_\beta(t)\rangle$ is
$$
\tau=\frac{|\Psi_\beta(t)\rangle\langle\Psi_\beta|}{\langle\Psi_\beta|\Psi_\beta(t)\rangle},
\qquad
\langle\Psi_\beta|\Psi_\beta(t)\rangle=\frac{Z(\beta+it)}{Z(\beta)}.
$$
Tracing out the right copy yields
$$
\tau_L=\mathrm{Tr}_R\,\tau
=\frac{1}{Z(\beta+it)}\sum_n e^{-(\beta+it)E_n}\,|E_n\rangle_L\langle E_n|.
$$
Thus $\tau_L$ is the thermal density matrix at complex inverse temperature $\beta+it$, and the thermal pseudo-entropy is
$$
S(\beta+it)\equiv-\,\mathrm{Tr}\bigl[\tau_L\log\tau_L\bigr]
= S_{\rm th}(\beta+it).
$$
For ordinary thermal entropy one has
$$
S_{\rm th}(\beta)=(1-\beta\partial_\beta)\log Z(\beta),
$$
so thermal pseudo-entropy is precisely the analytic continuation $\beta\to\beta+it$ of this expression [2411.08948].

At $t=0$, the construction reduces to the usual thermal ensemble:
$$
S(\beta+i0)=S_{\rm th}(\beta).
$$
This establishes thermal pseudo-entropy as a genuine extension of thermal entropy rather than an unrelated complex observable [2411.08948].

## 2. Analytic structure and generalized entropies

The same quantity admits an overlap-based expression,
$$
S(\beta+it)=\log Z(\beta+it)+\langle\theta(t)\rangle,
$$
with
$$
\theta_n(t)\equiv(\beta+it)E_n,
\qquad
\langle\cdots\rangle\equiv\frac{1}{Z(\beta+it)}\sum_n (\cdots)e^{-\theta_n(t)}.
$$
This form makes explicit that the complex phases inherited from the transition amplitudes enter directly into the entropy functional [2411.08948].

A Rényi generalization is defined for integer $n$ by
$$
S^{(n)}(\beta+it)
=\frac{1}{1-n}\log\frac{Z\bigl(n(\beta+it)\bigr)}{\bigl[Z(\beta+it)\bigr]^n},
$$
and the von Neumann limit is recovered as $n\to1$ [2411.08948]. Because the underlying reduced transition matrix is generally non-Hermitian, both the Rényi and von Neumann versions are in general complex-valued.

A further structural result is that, since $S(\beta+it)$ is analytic in the upper half-$t$ plane and vanishes as $\Im\,t\to-\infty$, its real and imaginary parts obey the Kramers-Kronig relations
$$
\Im\,S(\beta+it)=\frac{1}{\pi}\,P\!\!\int_{-\infty}^{\infty}ds\;
\frac{\Re\,S(\beta+is)}{s-t},
$$
and
$$
\Re\,S(\beta+it)= -\frac{1}{\pi}\,P\!\!\int_{-\infty}^{\infty}ds\;
\frac{\Im\,S(\beta+is)}{s-t}.
$$
In explicit examples including two-level systems, oscillators, Schwarzian theory, RMT, and 2D CFTs, direct numerical checks confirm these relations [2411.08948].

These properties distinguish thermal pseudo-entropy from a formal substitution $\beta\mapsto\beta+it$ performed only at the level of notation. The analyticity constraints imply that the real and imaginary parts are not independent data.

## 3. Relation to the spectral form factor and subregion dynamics

A central result is the relation to the finite-temperature spectral form factor
$$
g_\beta(t)=\left|\frac{Z(\beta+it)}{Z(\beta)}\right|^2.
$$
For Rényi index $n$,
$$
\Re\,S^{(n)}(\beta+it)
=\frac{1}{2(1-n)}
\Bigl[\log g_{n\beta}(t)-n\log g_\beta(t)\Bigr]
+S^{(n)}_{\rm th}(\beta).
$$
In particular, for $n=1$, the real part of thermal pseudo-entropy is directly tied to the spectral form factor [2411.08948].

This immediately imports the distinction between integrable and chaotic spectra. Integrable models such as finite-dimensional oscillators and Calogero-Sutherland produce purely periodic-oscillatory $g_\beta(t)$ and hence no ramp or plateau. Chaotic models such as RMT and SYK/Schwarzian exhibit the universal dip-ramp-plateau structure in $g_\beta(t)$, which implies a dip region where $\Re\,S$ decays, a ramp region with $\Re\,S\sim \tfrac12\log t$ when $g(t)\sim t$, and a late plateau at
$$
\Re\,S(\beta+it)\xrightarrow[t\to\infty]{}\tfrac12\,S_{\rm th}(2\beta)
$$
for non-degenerate spectra, with degeneracies only raising this value [2411.08948].

An important precursor appeared in the subregion setting of Nakata et al., who considered a transition matrix between the thermofield-double state and its time-evolved state in two-dimensional CFT and defined the pseudo-entropy of a single interval $A$. For an interval of length $\ell$ on the infinite line,
$$
S_p(A,t)=\frac{c}{6}\,\ln\!\Bigl[\frac{(\beta+it)^2}{\pi^2}\,
\sinh^2\!\Bigl(\frac{\pi\ell}{\beta+it}\Bigr)\Bigr].
$$
Its real part exhibits a four-stage behavior: an initial downward slope from the thermal entropy, a dip at $t\approx t_T\propto\sqrt{\beta\ell}$, a ramp for $t_T\ll t\ll t_H\sim \ell$, and a plateau at the vacuum entanglement entropy for $t\gg t_H$ [2109.00372].

On a spatial circle of circumference $L$, the same subregion construction is modified by torus correlators expressed through Jacobi $\theta$-functions, and one obtains a plateau
$$
S_{\rm plateau}(\ell,L)=\frac{c}{3}\ln\!\Bigl(\frac{L}{\pi}\sin\frac{\pi\ell}{L}\Bigr),
$$
together with revivals at integer multiples of $L$ in integrable theories, or a self-averaging constant in chaotic or holographic CFTs [2109.00372]. This subregion framework does not redefine thermal pseudo-entropy itself, but it shows how the same transition-matrix formalism extends from the whole thermal system to spatial subsystems.

## 4. Continuous spectra, model classes, and universal logarithms

When the density of states near the ground energy $E_0$ behaves as
$$
\rho(E)\sim(E-E_0)^\gamma\qquad (E\gtrsim E_0),
$$
the large-$\beta$ partition function scales as
$$
Z(\beta)\sim e^{-\beta E_0}\,\beta^{-1-\gamma},
$$
and the ordinary thermal entropy behaves as
$$
S_{\rm th}(\beta)\simeq -(1+\gamma)\log\beta+\cdots.
$$
After analytically continuing $\beta\to\beta+it$ and averaging over fast oscillations, one finds for $t\gg\beta$ but below the quantum recurrence time
$$
S(\beta+it)\simeq -(1+\gamma)\log t + O(1).
$$
Hence the coefficient of the logarithmic tail is universal and equals $1+\gamma$ [2411.08948].

Several model classes realize this statement explicitly. In Schwarzian theory, where the edge density implies $\gamma=\tfrac12$, the one-loop exact result of Stanford-Witten gives
$$
S(\beta+it)=S_0+\frac{4\pi^2 C}{\beta+it}-\tfrac32\log(\beta+it)+\cdots,
$$
so that for $t\gg\beta$ one obtains $-\tfrac32\log t+\cdots$. Random Matrix Theory with Wigner semi-circle density also has $\gamma=\tfrac12$ and yields the same $-\tfrac32\log t$ behavior in the real part. For a decompactified free scalar in 2D CFT, $\rho(E)\sim E^{-1/2}$ gives $\gamma=-\tfrac12$, while for $N$ non-compact scalars one has $\gamma=\tfrac N2-1$ and therefore
$$
\Re\,S\sim -\frac{N}{2}\log t
$$
[2411.08948].

By contrast, theories with a discrete gap, including the free fermion CFT, the compact boson at finite radius, and a holographic CFT with sparse primaries, have $\gamma=-1$ by convention and no logarithmic tail [2411.08948]. Symmetric-orbifold CFTs $\mathcal C^N/S_N$ provide large-$N$ integrable examples whose time-averaged spectral form factor exhibits a dip-ramp-plateau structure with steps reflecting contributions of different twisted sectors; their thermal pseudo-entropy mirrors the seed theory at early times and eventually oscillates about the plateau $\tfrac12 S_{\rm th}(2\beta)$ [2411.08948].

The universal coefficient $1+\gamma$ indicates that thermal pseudo-entropy is sensitive to the edge scaling of the density of states rather than only to bulk thermodynamic data. This suggests a use as a spectral probe complementary to standard partition-function asymptotics.

## 5. Modular-Hamiltonian interpretation and short-time response

Later work embeds thermal pseudo-entropy into the general theory of real-time pseudo-entropy. For a bipartite Hilbert space
$$
\mathcal H=\mathcal H_A\otimes\mathcal H_{\bar A},
$$
and two non-orthogonal pure states $\ket{\psi}$ and $\ket{\phi}$, one defines the normalized transition operator
$$
\tau^{\psi|\phi}=\frac{\ket{\psi}\bra{\phi}}{\langle\phi|\psi\rangle},
\qquad
\mathrm{Tr}\,\tau^{\psi|\phi}=1,
$$
its reduced version
$$
\tau^{\psi|\phi}_A=\mathrm{Tr}_{\bar A}\,\tau^{\psi|\phi},
$$
and the pseudo-entropy
$$
S_A^{\psi|\phi}=-\,\mathrm{Tr}_A\Bigl[\tau^{\psi|\phi}_A\,
\log\bigl(\tau^{\psi|\phi}_A\bigr)\Bigr].
$$
For real-time evolution with $\ket{\phi}=\ket{\Psi}$ and $\ket{\psi}=e^{-iHt}\ket{\Psi}$, one finds
$$
S_A(t,0)=S_A(0)-it\,\langle K_A(H-\langle H\rangle)\rangle+\mathcal O(t^2),
$$
where
$$
K_A=-\log\rho_A.
$$
The imaginary part is therefore governed at leading order by the symmetrized covariance of $H$ and $K_A$, while the real part is governed by their commutator [2606.14208].

Thermal pseudo-entropy arises as a special case in a Schmidt-diagonal model with probabilities $p_n=e^{-\beta E_n}/Z(\beta)$. In that setting one obtains
$$
M(\theta)=\frac{Z(\beta-\theta)}{Z(\beta)},
$$
and
$$
S_A(\theta)=\log Z(\beta-\theta)-(\beta-\theta)\frac{Z'(\beta-\theta)}{Z(\beta-\theta)}
=\bigl[1-s\partial_s\bigr]\log Z(s),\qquad s=\beta-\theta.
$$
Analytically continuing $s\to\beta+it$ yields the standard thermal pseudo-entropy
$$
S_{\rm th}(s)=\bigl[1-s\partial_s\bigr]\log Z(s),\qquad s=\beta+it.
$$
Its small-$t$ expansion is
$$
S_{\rm th}(\beta+it)=S_{\rm th}(\beta)-it\,\beta\,\mathrm{Var}_\beta(H)+\mathcal O(t^3),
$$
because
$$
\langle KH\rangle-\langle K\rangle\langle H\rangle
=\beta\,\langle(H-\langle H\rangle)^2\rangle
$$
in the thermal ensemble [2606.14208].

This formulation sharpens the meaning of the imaginary part. It is not merely a branch artifact; rather, in the language of the 2026 analysis, it is a time-oriented modular response generated by the correlation between microscopic time evolution and subsystem coarse graining [2606.14208].

## 6. Finite-dimensional thermal transitions, holographic comparisons, and terminological divergence

A distinct thermal use of pseudo-entropy appears in the two-qubit analysis of Ali-Akbari et al., where the system transitions from an initial thermal state to a final thermal state at fixed temperature under an external magnetic field. The initial Hamiltonian is
$$
H_i=\sigma_1^x\otimes\sigma_2^x+\sigma_1^y\otimes\sigma_2^y+\sigma_1^z\otimes\sigma_2^z,
$$
while the final Hamiltonian adds Zeeman terms,
$$
H_f=H_i-\mu B_1(\sigma_1^z\otimes I_2)-\mu B_2(I_2\otimes\sigma_2^z).
$$
The thermal density matrices are
$$
\rho_i=\frac{e^{-\beta H_i}}{Z_i},\qquad
Z_i=3e^{-\beta}+e^{3\beta},
$$
and
$$
\rho_f=\frac{e^{-\beta H_f}}{Z_f}.
$$
The mixed-to-mixed transition matrix is defined by
$$
\tau\equiv T=\frac{\rho_f\rho_i}{\mathrm{Tr}(\rho_f\rho_i)},
$$
with reduced transition matrix $\tau_A=\mathrm{Tr}_B\,T$, and the subsystem pseudo-entropy
$$
S_{\rm pseudo}\equiv S_A=-\,\mathrm{Tr}_A[\tau_A\ln\tau_A].
$$
Because $\rho_i$ and $\rho_f$ are both block-diagonal in the computational basis, $\tau_A$ is $2\times2$ diagonal with eigenvalues $\lambda_{1A},\lambda_{2A}$, leading to
$$
S_A(B,\beta)=-\bigl[\lambda_{1A}\ln\lambda_{1A}+\lambda_{2A}\ln\lambda_{2A}\bigr].
$$
For the symmetric analytic continuation $B_1=iB$, $B_2=B$, the small-$B$ expansion gives
$$
S_A=\ln2+B^2\cdot F(\beta)+O(B^3),
$$
with real and imaginary parts both quadratic in $B$; at $B=0$ both eigenvalues approach $\tfrac12$, so $S_A=\ln2$ [2602.20843].

The same paper compares this finite-dimensional pseudo-entropy to holographic timelike entanglement entropy (HTEE), computed by replacing the Ryu-Takayanagi surface with a timelike extremal surface in a 5d Einstein-Maxwell background with constant magnetic field. In the exact low-energy solution, both $\Re\,S_{\rm HTEE}$ and $\Im\,S_{\rm HTEE}$ scale as $\sqrt B$ at small $B$, whereas the two-qubit pseudo-entropy scales as $B^2$. Moreover, HTEE diverges as $B\to0$, while the qubit model smoothly tends to $\ln2$ [2602.20843]. The stated physical reason is that the qubit model has only four states and the field appears as a Hamiltonian parameter, whereas in the holographic dual the field backreacts on an infinite-degree-of-freedom geometry.

This comparison is significant because it blocks an overly direct identification between pseudo-entropy in simple thermal quantum systems and timelike holographic entanglement observables. Although both quantities can acquire an imaginary part from timelike analytic continuation, their scaling laws and limiting behavior differ fundamentally [2602.20843].

A separate terminological issue arises in heavy-ion phenomenology, where “pseudo-entropy” has been defined not from a reduced transition matrix but from a normalized hadron transverse-momentum spectrum,
$$
S'_{\rm hadron}=-\sum_j f(p_{T_j})\ln[f(p_{T_j})],
$$
with $p_T$ bins of width $\Delta p_T=0.1$ GeV/$c$ [2112.09473]. In that context it is an information-theoretic measure of the disorder of the hadron $p_T$ distribution, extracted from a Tsallis-Pareto-type fit, and the authors explicitly note that it is not the true thermodynamic entropy of the fireball, but rather a proxy based on single-particle $p_T$ distributions [2112.09473]. This distinction helps prevent confusion between the quantum-information notion of thermal pseudo-entropy based on non-Hermitian transition matrices and an unrelated Shannon-type observable in collision phenomenology.

Source: https://www.emergentmind.com/topics/thermal-pseudo-entropy