Hawking Radiation in non-Hermitian Microscopic Analogues

AI summary

Overview

Research area: Analogue gravity in open quantum systems — specifically, how black-hole and white-hole horizons, and fermionic Hawking-pair correlations, emerge from a microscopic Lindblad master equation rather than from a hand-written non-Hermitian Hamiltonian.

Technical level: Advanced. The paper assumes familiarity with Lindblad master equations, third quantization of quadratic fermionic Liouvillians, biorthogonal quantum mechanics, non-Hermitian exceptional points, and Bogoliubov scattering. The prose is readable, but the derivations are not beginner material.

Scope (one sentence): The paper derives an exactly solvable quadratic fermionic Lindbladian whose rapidity spectrum contains black-hole and white-hole exceptional cones, extracts an effective Painlevé-Gullstrand geometry from it, and formulates many-body witnesses of analogue Hawking-pair correlations for the resulting non-Hermitian single-particle problem.

What This Paper Is About

Most analogue-horizon proposals in non-Hermitian systems start from an effective gain/loss Hamiltonian with non-reciprocal hopping, which describes only the no-jump (conditional) short-time dynamics of an open system. This paper instead starts from a genuinely Hermitian fermionic tight-binding chain coupled to Markovian reservoirs, derives the non-Hermitian horizon model as a limit of that microscopic Lindblad dynamics, and then asks which horizon features survive when quantum jumps, relaxation, and the full long-time steady state are included. It then builds a fermionic many-body (Nambu) extension to test whether the outgoing channels actually carry Hawking-partner correlations.

Key Contributions

1. A microscopic Lindblad realization of an analogue horizon. A Hermitian tight-binding chain with $N$ unit cells and two fermionic modes $a_j, b_j$ per cell, coupled to four linear jump operators with rates $\kappa$ and $\gamma$ ($\kappa > |\gamma|$), produces — in the no-jump limit — an effective non-Hermitian Hamiltonian with a gain/loss term $i\gamma\sum_j(a_j^\dagger a_j - b_j^\dagger b_j)$ and a non-reciprocal next-nearest-neighbor hopping term of strength $\kappa$, plus an overall shift $-2i\kappa N$ that pushes energies into the lower-half complex plane.

2. Exact solution of the full Lindbladian and identification of a black-hole/white-hole pair. Using third quantization (doubling to $4N$ Majorana fermions from $2N$ complex fermions), the damping matrix $X$ block-diagonalizes into two sectors $\tilde X_c$ and $\tilde X_d$ whose spectra are identical except for the sign of the non-reciprocal next-nearest-neighbor hopping. The full Lindblad spectrum therefore contains both a black-hole cone $\beta_c(q)$ and a white-hole cone $\beta_d(q)$, whereas the effective non-Hermitian Hamiltonian contains only one.

3. An emergent Painlevé-Gullstrand geometry with four momentum branches. Near $q=\pi$, the exceptional points form cones $\gamma = \mp f(r)k \pm \sqrt{\tau^2 k^2}$ (with $\tau=1$), matched to the metric $ds^2 = [1-f^2(r)],dt^2 \pm 2f(r),dt,dr + (dr)^2$, where $f(r)$ is a two-interface $\tanh$ profile with asymptotic values $\kappa_1$, $\kappa_2$, $\kappa_1$ and interface value $\kappa_3 = (\kappa_1+\kappa_2)/2$. Solving for $k$ yields black-hole branches $k_{\rm BH}^{\pm}$ and white-hole branches $k_{\rm WH}^{\pm}$, with the $+$ ($-$) sign denoting particle (antiparticle) channels at asymptotic infinity.

4. Single-particle channel identification and many-body Hawking-pair witnesses. A biorthogonal flux built from asymptotic left and right Bloch eigenvectors selects the outgoing exterior and interior channels. A fermionic Gaussian Nambu extension with an antisymmetric pairing matrix turns the stationary many-body problem into frequency-resolved Bogoliubov scattering, yielding Hawking and partner occupations, their anomalous correlation, and the real-space covariance, checked against covariance-positivity and Pauli constraints.

Main Findings

• The effective non-Hermitian limit reproduces the known gain/loss non-reciprocal model. Taking the no-jump conditional evolution of the Lindblad equation gives $H^{\rm eff}{\rm nH}[\kappa] = H{\rm TB} + i\gamma\sum_j(a_j^\dagger a_j - b_j^\dagger b_j) + \frac{\kappa}{2}\sum_j(-a_{j+1}^\dagger a_j + a_j^\dagger a_{j+1} + b_{j+1}^\dagger b_j - b_j^\dagger b_{j+1})$, i.e. exactly the gain/loss lattice with non-reciprocal next-nearest-neighbor hopping previously used to mimic horizons.

• The horizon structure survives the full Lindblad dynamics. The exact rapidity spectrum retains a black-hole sector ($\tilde X_c$) and a white-hole sector ($\tilde X_d$), differing only by the sign of the non-reciprocal hopping. This is the paper's central point: the analogue geometry is not an artifact of discarding quantum jumps.

• Steady-state particle densities show sharp transitions at both horizons. Using the momentum-space correlation matrix $C(q)$ from the Lyapunov equation $X^\dagger(q)C(q) + C(q)X(q) = iY(q)$, the steady-state densities $|\langle a^\dagger a\rangle_{\rm SS}|$ and $|\langle b^\dagger b\rangle_{\rm SS}|$ display sharp transitions at $r = r_1$ and $r = r_2$, for both the black-hole cone (Fig. 1) and the white-hole cone (Fig. 2), and for both particle and antiparticle momentum channels. These figures use $r_1 = 100$, $r_2 = 200$, $N = 300$, $l = 1$, $\gamma = 0.3$, with zoom windows $r \in {r_1-5, r_1+5}$ and $r \in {r_2-5, r_2+5}$.

• The horizons do not cancel in the steady state. The authors state explicitly that the equal-momentum density plots alone do not establish whether the black-hole and white-hole contributions cancel, and that the steady-state current (constructed from the correlation matrix) is the observable that addresses this.

• Validity conditions are stated, not assumed. The momentum-space derivation requires the $f(r)$-commutator corrections to be subleading, controlled by $|(\kappa_2-\kappa_1)q/(2l)| \ll |\gamma|$ and $|(\kappa_2-\kappa_1)/(2l)| \ll |\tau|$, with Schwarzschild-like parameter bounds $(\kappa_1-\kappa_2)/l = 1/r_1$, $\kappa_1 > 1$, $\kappa_2 < 1$, $\kappa_3 = 1$.

• Real-space covariance shows cross-horizon ridges. The Hawking-partner correlations produce cross-horizon ridges in the covariance that follow the expected equal-travel-time trajectory and are strongly suppressed in the no-crossing and uniform control cases.

• Fermionic positivity is satisfied. The frequency-resolved Hawking-partner correlations obey the fermionic covariance-positivity constraint; the normalized anomalous correlation stays below unity, as required by Pauli exclusion. The authors are careful to note that this constraint demonstrates compatibility with a physical fermionic state and is not, by itself, a nonclassicality witness — a fermionic entanglement criterion is needed for that.

Methodology in Plain English

1. Build the open system from honest ingredients. Write down an ordinary, Hermitian tight-binding chain with two sites per unit cell. Couple it to four engineered dissipative channels (two local particle-gain processes, two bond-type processes). All rates are real and non-negative, enforced by $\kappa > |\gamma|$.

2. Ask what happens with no jumps detected. The standard Lindblad formalism gives a conditional evolution governed by an effective non-Hermitian Hamiltonian. Expanding it reveals a gain/loss imbalance between the two sublattices and a directional next-nearest-neighbor hopping — the ingredients associated with analogue horizons.

3. Make the hopping space-dependent. Replace the constant $\kappa$ by a two-interface $\tanh$ profile $f(r)$ that equals $\kappa_1$ far to the left, $\kappa_2$ in the middle, and $\kappa_1$ again far to the right, with interfaces at $r_1$ and $r_2$. This gives a black-hole-like interface at one end and a white-hole-like interface at the other.

4. Solve the full dissipative problem exactly. Because the Liouvillian is quadratic in fermions, it can be diagonalized by doubling the number of modes (Majorana doubling) and computing the spectrum of a single non-Hermitian damping matrix. This is exact — no truncation in time, no neglect of jumps.

5. Go to momentum space and look for exceptional points. Under periodic boundary conditions, the damping matrix becomes a $2\times2$ Bloch problem in each sector. The exceptional points (where eigenvalues coalesce) form cones near $q=\pi$; identifying these cones with light-cone kinematics gives the effective metric.

6. Compute steady-state observables. Solve the Lyapunov equation for the momentum-space correlation matrix, then convert back to Majorana and fermionic bases to get number densities and currents. Plot densities at the momenta that sit on the exceptional cones.

7. Add a pairing sector for Hawking pairs. Introduce a fermionic Gaussian Nambu extension with an antisymmetric pairing matrix, treat the stationary state as a frequency-resolved Bogoliubov scattering problem, and read off Hawking occupations, partner occupations, anomalous correlations, and the real-space covariance. Run no-crossing and uniform control cases to confirm the ridges are horizon-specific.

Why This Matters

Impact on research. The paper closes a gap that is often glossed over in analogue-gravity work: non-Hermitian models are usually introduced as effective descriptions, leaving open whether the horizon survives the unconditional, trace-preserving dynamics. By deriving the horizon model from a Lindblad equation and then solving that equation exactly, the authors show that a black-hole and a white-hole sector coexist in the full rapidity spectrum — a structural result that effective-Hamiltonian treatments cannot see. It also links three normally separate literatures: microscopic open-system engineering, non-Hermitian exceptional-point geometry, and fermionic many-body probes of Hawking radiation.

Potential real-world applications (as motivated by this line of work; the paper itself reports no device demonstration):

• Cold-atom and optical-lattice quantum simulators with engineered dissipation, where loss and gain terms can be tuned site-by-site and bond-by-bond.
• Photonic and circuit-QED lattices supporting non-reciprocal hopping and gain/loss patterning, which are natural hosts for the $\kappa$- and $\gamma$-type terms in the model.
• Quantum simulation of curved-spacetime kinematics in tabletop fermionic platforms, as a testbed for horizon physics where gravitational experiments are inaccessible.
• Fermionic correlation metrology, since cross-horizon covariance ridges and frequency-resolved partner correlations are the kind of nonlocal quantities that dissipation-engineered lattices can in principle measure.

Industry relevance. Dissipation is usually treated as an adversary in quantum hardware; this work treats it as a design resource. That framing is directly relevant to groups engineering controlled loss in superconducting and photonic processors, to benchmarking of non-Hermitian mode structure in hardware, and to quantum-sensing schemes that exploit nonlocal correlations rather than single-site occupancies.

Future Directions

• Experimental realization and readout. The paper is theoretical; the natural next step is a concrete platform proposal specifying which lattice, which loss/gain mechanism, and which measurement protocol would reveal the steady-state density transitions at $r_1$ and $r_2$ and the cross-horizon covariance ridges.

• A genuine fermionic entanglement criterion. The authors explicitly note that covariance positivity only confirms compatibility with a physical fermionic state. Applying an entanglement criterion adapted to the fermionic algebra would upgrade the correlations from "consistent with a physical state" to "genuinely quantum."

• Beyond quadratic dynamics. The exact solvability rests on the Lindbladian being quadratic. Whether the black-hole/white-hole sector structure, and the Hawking-partner signatures, survive interactions or number-nonconserving terms is open.

• Which vacuum? The gravitational discussion in the introduction emphasizes that the physically relevant state for a collapse-formed black hole is the Unruh vacuum, empty at past null infinity and regular on the future horizon. Establishing which state the steady state of this Lindbladian corresponds to — and whether a late-time thermal flux follows — is not settled by the results reported here.

• Scale and dimensionality. The reported steady-state calculations use $N = 300$ with $r_1 = 100$, $r_2 = 200$, $l = 1$, $\gamma = 0.3$. Extending to two-dimensional horizon geometries and to parameter regimes where the $f(r)$-commutator corrections are larger would test how robust the cone structure is.

Target Audience

Advanced graduate students and researchers in quantum open systems, non-Hermitian physics, analogue gravity, and cold-atom or photonic quantum simulation. It will be most valuable to readers who already know the Lindblad formalism and third quantization and want to see how horizon geometry and Hawking-pair correlations can be derived from — rather than postulated alongside — a microscopic dissipative model. Readers seeking an introduction to analogue Hawking radiation should start elsewhere first; readers seeking a rigorous test of whether non-Hermitian analogue horizons survive decoherence and jumps will find this paper directly aimed at them.

Authors’ abstract

Non-Hermitian systems have recently emerged as a convenient platform to realize black-hole physics. Here, we propose a microscopic open quantum system, namely a fermionic tight-binding chain coupled to Markovian reservoirs, which emulates black-hole properties. The resulting quadratic Lindbladian admits an effective non-Hermitian limit that produces a gain and loss lattice model with a non-reciprocal next-nearest-neighbor hopping. The rapidity spectrum of the full Lindbladian forms tilted exceptional cones that define an effective Painlevé-Gullstrand geometry and separates into black-hole and white-hole sectors, while steady-state particle densities and currents retain clear signatures at the horizon positions. We then analyze the effective non-Hermitian Hamiltonian that arise from the microscopic quantum system. The biorthogonal flux identifies the outgoing exterior and interior channels. Finally, we introduce a fermionic Gaussian Nambu extension and formulate Hawking-pair witnesses through frequency-resolved scattering, nonlocal density-density correlations, and anomalous Hawking-partner amplitudes. The resulting correlations display analogue Hawking radiation signatures, while the frequency-resolved Hawking-partner correlations satisfy the fermionic covariance-positivity constraint. These results establish a connection between microscopic open-system dynamics, emergent non-Hermitian geometry, and fermionic many-body probes of analogue Hawking radiation.