---
title: Swanson Hamiltonian Overview
url: https://www.emergentmind.com/topics/swanson-hamiltonian
type: topic
---

# Swanson Hamiltonian Overview

Searching arXiv for recent and foundational papers on the Swanson Hamiltonian to ground the article in the literature.
The Swanson Hamiltonian is a non-Hermitian quadratic oscillator Hamiltonian built from bosonic creation and annihilation operators and characterized by number-nonconserving quadratic terms. In its standard form,
\[
H_S=\omega\,a^\dagger a+\alpha\,a^2+\beta\,(a^\dagger)^2
\]
or, with the zero-point term retained,
\[
H_S=\omega\Bigl(a^\dagger a+\tfrac12\Bigr)+\alpha\,a^2+\beta\,(a^\dagger)^2,
\]
with \([a,a^\dagger]=1\) and typically real parameters \(\omega,\alpha,\beta\) [1501.03922]. For \(\alpha\neq\beta\), the Hamiltonian is non-Hitian because the anomalous terms \(a^2\) and \((a^\dagger)^2\) are not adjoints of each other with equal coefficients, yet it remains a central example of \(\mathcal{PT}\)-symmetric and pseudo-Hermitian quantum mechanics. Its importance derives from three connected facts: it is exactly diagonalizable in broad parameter regimes, it admits an equivalent Hermitian representation under a similarity transformation, and it has become a template for studying exceptional points, pseudo-bosonic structures, supersymmetric extensions, noncommutative deformations, time-dependent metrics, and multi-mode generalizations [2401.17189].

## 1. Definition and operator structure

The defining feature of the Swanson Hamiltonian is the coexistence of the standard oscillator term with pair-annihilation and pair-creation terms. In the operator language used throughout the literature,
\[
H_S=\omega\Bigl(a^\dagger a+\tfrac12\Bigr)+\alpha a^2+\beta(a^\dagger)^2,
\]
where \(\omega\) is the oscillator frequency and \(\alpha,\beta\) control the anomalous quadratic couplings [1501.03922]. The term \(\omega(a^\dagger a+\tfrac12)\) is the usual harmonic-oscillator Hamiltonian, while \(\alpha a^2\) and \(\beta(a^\dagger)^2\) are pair-annihilation and pair-creation contributions [1501.03922].

For real parameters, Hermiticity is lost whenever \(\alpha\neq\beta\), since
\[
H_S^\dagger=\omega\,a^\dagger a+\alpha\,(a^\dagger)^2+\beta\,a^2
\]
and therefore \(H_S\neq H_S^\dagger\) unless the two anomalous coefficients match in the appropriate way [2401.17189]. Nevertheless, the model is widely treated as \(\mathcal{PT}\)-symmetric in the standard bosonic realization with parity acting as \(a\mapsto -a\), \(a^\dagger\mapsto -a^\dagger\), and time reversal sending \(i\mapsto -i\) [2401.17189]. In more general complex-parameter treatments, the \(\mathcal{PT}\)-symmetry condition is stated as \(\beta=\alpha^*\) [2606.04062].

This elementary-looking quadratic form is sufficiently flexible to generate several qualitatively different physical regimes. Depending on the signs of the effective mass parameter and of the squared effective frequency, the model can represent an ordinary harmonic oscillator, a negative-mass oscillator, an inverted oscillator or parabolic barrier, and a negative-mass parabolic barrier [2108.01004]. This classification is one reason the Swanson Hamiltonian functions as a canonical laboratory for non-Hermitian spectral theory.

## 2. Pseudo-Hermiticity and Hermitian equivalence

A basic structural result is that the Swanson Hamiltonian is pseudo-Hermitian. One introduces a positive-definite metric operator \(\zeta=\rho^2\) such that
\[
H_S^\dagger=\zeta\,H_S\,\zeta^{-1},
\]
and then defines the equivalent Hermitian Hamiltonian
\[
h=\rho\,H_S\,\rho^{-1},
\]
which is manifestly Hermitian in the reference Hilbert space [1501.03922]. This is the standard quasi-Hermitian construction for the model.

In Bogoliubov or Dyson form, one seeks a squeeze operator
\[
S(\theta)=\exp\!\Bigl[\frac{\theta}{2}(a^2-(a^\dagger)^2)\Bigr]
\]
that implements
\[
b=S(\theta)\,a\,S(\theta)^{-1}=\cosh\theta\,a+\sinh\theta\,a^\dagger,
\]
with a corresponding expression for \(b^\dagger\), so that the transformed Hamiltonian becomes
\[
h=S(\theta)\,H_S\,S(\theta)^{-1}=\Omega\Bigl(b^\dagger b+\tfrac12\Bigr),
\qquad
\Omega=\sqrt{\omega^2-4\alpha\beta}.
\]
The off-diagonal terms cancel when
\[
\tanh(2\theta)=\frac{\alpha-\beta}{\omega},
\]
and the reality of the transformed oscillator frequency requires
\[
\omega^2-4\alpha\beta\ge 0
\]
[2401.17189].

In coordinate space, the equivalent Hermitian operator can also be written as a Schrödinger Hamiltonian. One explicit realization uses
\[
\Upsilon=\exp\!\Bigl[-\,\frac{\beta-\alpha}{2(\omega-\alpha-\beta)}\,\frac{x^2}{b_0^2}\Bigr],
\]
which maps \(H_{SW}\) to
\[
\mathfrak h_{SW}
=
\Upsilon\,H_{SW}\,\Upsilon^{-1}
=
\frac{p^2}{2m}+\frac12 m\Omega^2 x^2,
\]
where
\[
m=\frac{\hbar}{(\omega-\alpha-\beta)b_0^2},
\qquad
\Omega=\sqrt{\omega^2-4\alpha\beta}
\]
[2212.13173]. This form makes explicit that the Swanson model can be interpreted as a harmonic or inverted oscillator depending on the sign of \(\Omega^2\).

The positivity of the metric is central. In one formulation,
\[
\eta=\rho^\dagger\rho=\exp[\theta(a^{\dagger2}-a^2)]
\]
is positive-definite for \(|\tanh(2\theta)|<1\), equivalently in the real-spectrum regime [2606.04062]. This metric underlies the biorthogonal inner product and the associated physical interpretation of observables.

## 3. Spectrum, phases, and exceptional points

The spectral behavior of the Swanson Hamiltonian is controlled by the effective frequency
\[
\Omega=\sqrt{\omega^2-4\alpha\beta}.
\]
When \(\Omega\in\mathbb R\), the model has a real discrete spectrum
\[
E_n=\Omega\Bigl(n+\tfrac12\Bigr),\qquad n=0,1,2,\dots
\]
in the Hermitian representation, and therefore the original non-Hermitian Hamiltonian shares the same discrete spectrum in the unbroken regime [2401.17189].

A more refined parameter-space classification uses
\[
m=\frac{\hbar}{(\omega-\alpha-\beta)b_0^2},
\qquad
\Omega^2=\omega^2-4\alpha\beta.
\]
With these quantities, four regions arise in the \((\alpha/\omega,\beta/\omega)\)-plane. Region I, with \(m>0\) and \(\Omega^2>0\), is an ordinary harmonic oscillator with discrete real spectrum
\[
E_n=\hbar|\Omega|(n+\tfrac12).
\]
Region III, with \(m<0\) and \(\Omega^2>0\), is a negative effective mass oscillator with discrete real spectrum of opposite sign,
\[
E_n=-\hbar|\Omega|(n+\tfrac12).
\]
Regions II and IV, where \(\Omega^2<0\), correspond to parabolic-barrier regimes and support resonant ladders together with a real continuous band [2108.01004].

In the broken-\(\mathcal{PT}\) treatment based on complex scaling, the same transition is described by taking \(\Omega=i\Gamma\) when \(\omega^2<4\alpha\beta\). Then the discrete energies are written as
\[
E_n^\pm=\pm i\hbar|\Omega|\Bigl(n+\tfrac12\Bigr),
\]
so the spectrum appears in complex-conjugate form [2405.04599]. This formulation emphasizes resonance behavior and non-unitary dynamics.

Exceptional points occur on the critical hypersurface
\[
\omega^2-4\alpha\beta=0.
\]
At this boundary, the effective frequency vanishes and the discrete levels collapse [2108.01004]. One description states that all discrete levels \(E_n=\pm\hbar|\Omega|(n+\tfrac12)\) coalesce at \(E=0\), and the Hamiltonian develops an infinite-order Jordan block at that energy [2108.01004]. A separate formulation identifies the same locus as the transition where the squeezing angle diverges, the similarity map becomes singular, and the Hamiltonian ceases to be diagonalizable [2401.17189]. These are not contradictory descriptions; they reflect different operator realizations and spectral embeddings.

One subtle boundary is \(\omega=\alpha+\beta\). In that limit the effective mass diverges in the coordinate representation, and the Hamiltonian reduces, after a different similarity transformation, to
\[
h=\frac{\hbar(\alpha+\beta)}{2}\bigl(2x\,d/dx+1\bigr),
\]
with discrete generalized eigenfunctions \(x^n/\sqrt{n!}\) and \(((-1)^n/\sqrt{n!})\delta^{(n)}(x)\) carrying positive and negative ladders of real eigenvalues [2108.01004]. This boundary is distinguished by anti-pseudo-Hermiticity rather than by exceptional points.

## 4. Biorthogonal, rigged-Hilbert-space, and complex-scaling formulations

Because the Swanson Hamiltonian is non-Hermitian, the conventional Hilbert-space eigenbasis is not always sufficient. In the non-\(\mathcal{PT}\)-symmetric phase, the eigenfunctions are treated as generalized eigenfunctionals in a Gel'fand triplet
\[
\Phi\subset\mathcal H\subset\Phi^\times,
\]
with right and left eigenfunctionals satisfying
\[
H^\times|\phi_E\rangle=E|\phi_E\rangle,
\qquad
\langle\psi_E|H^\times=E\langle\psi_E|
\]
[2108.01004]. A similarity map \(\Upsilon\) sends the problem to an oscillator-type Hermitian operator, after which one constructs either Hermite-function states or parabolic-cylinder-function states and then transforms them back [2108.01004].

In this setting, biorthogonality is imposed through a positive operator \(U=\Upsilon^2\), and the right and left generalized eigenfunctions provide completeness in the distributional sense [2108.01004]. This is particularly useful in the parabolic-barrier regime, where resonant states and continuous-spectrum states coexist.

An alternative formulation uses the Complex Scaling Method. One introduces the non-unitary dilatation
\[
U(\theta)=\exp\!\Bigl(\frac{i\theta}{2}(xp+px)\Bigr),
\qquad
x\mapsto e^{i\theta}x,\quad p\mapsto e^{-i\theta}p,
\]
which transforms the coordinate-space Hamiltonian into a complex-scaled operator \(H(\theta)\) [2405.04599]. Further Bogoliubov-type maps reduce \(H(\theta)\) to a rotated harmonic or inverted oscillator \(\mathfrak h(\theta)\) [2405.04599].

The complex-scaled Hamiltonian admits right and left eigenfunctions \(\phi_n^\pm(\theta,x)\) and \(\psi_n^\pm(\theta,x)\) satisfying
\[
H(\theta)\phi_n^\pm=E_n^\pm\phi_n^\pm,
\qquad
H_c(\theta)\psi_n^\pm=E_n^{\pm*}\psi_n^\pm,
\]
together with biorthonormality and completeness relations [2405.04599]. A central advantage of the complex-scaled representation is that it yields square-integrable eigenfunctions for suitable \(\theta\)-sectors, whereas the rigged-Hilbert-space approach works directly with distributional Gamow-type states [2405.04599]. Both approaches are stated to reproduce the same physical decay rates and scattering quantities once the scaling sector is chosen appropriately [2405.04599].

A related but distinct development is the fully pseudo-bosonic reformulation, where one replaces the canonical pair \((a,a^\dagger)\) by a \(\mathcal D\)-pseudo-bosonic pair \((a,b)\) satisfying \([a,b]f=f\) on a common dense invariant domain. In that setting the Hamiltonian
\[
H=\omega\,ba+\lambda(b^2+a^2),\qquad \omega>2\lambda,
\]
is diagonalized by a Bogoliubov-type transformation to
\[
H=\Omega\,BA+\gamma\,\mathbf 1,
\]
with a complete biorthogonal eigensystem and associated bi-coherent states [2210.11326]. This formulation generalizes the Swanson construction beyond the standard adjoint relation \(b=a^\dagger\).

## 5. Supersymmetry, pseudo-supersymmetry, and generalized differential realizations

A substantial branch of the literature embeds the Swanson Hamiltonian into supersymmetric and pseudo-supersymmetric frameworks. One starting point is the generalized differential realization
\[
\tilde a=A(x)\frac{d}{dx}+B(x),
\qquad
\tilde a^\dagger=-A(x)\frac{d}{dx}+B(x)-A'(x),
\]
leading to the generalized Swanson Hamiltonian
\[
H_{GS}
=
\omega\bigl(\tilde a^\dagger \tilde a+\tfrac12\bigr)
+\alpha\tilde a^2+\beta\tilde a^{\dagger2}
\]
[1101.1040]. When \([\tilde a,\tilde a^\dagger]=\) constant, the Hermitian equivalent can be transformed locally into a harmonic-oscillator equation, but the model need not be globally isospectral to the harmonic oscillator because the coordinate map \(z(x)=\int dx/A(x)\) may cover a finite interval rather than the full real line [1101.1040]. This non-isospectrality is one of the recurring subtleties of generalized Swanson systems.

Pseudo-supersymmetric factorization begins from a first-order differential operator \(\xi=a(x)d/dx+b(x)\) and the non-Hermitian Hamiltonian
\[
H_-=\omega(\xi^\dagger\xi+\tfrac12)+\alpha\xi^2+\beta\xi^{\dagger2},
\qquad \alpha\neq\beta.
\]
A similarity transformation \(h_-=\rho H_-\rho^{-1}\) yields a Hermitian Hamiltonian factorized as \(h_-=\mathcal A^\dagger\mathcal A\), with partner \(h_+=\mathcal A\mathcal A^\dagger\), intertwining operator \(\eta_1=\rho^{-1}\mathcal A\rho\), and pseudo-supersymmetric relation \(\eta_1H_-=H_+\eta_1\) [1108.0106]. For a specific choice of \(a(x)\), \(b(x)\), and \(\tilde b(x)\), this construction produces isotonic and nonlinear isotonic oscillator potentials \(V_+(x)\) and \(V_-(x)\), along with exact eigenvalues and confluent-hypergeometric or Laguerre eigenfunctions [1108.0106].

The second-derivative pseudo-supersymmetric generalization pushes the factorization further. One defines a diagonal SSUSY Hamiltonian
\[
H_{SS}=\operatorname{diag}(h_+,h_-),
\qquad
h_\pm=-\,d\,\tilde a^2(x)\,d+V_\pm(x),
\]
together with second-order differential operators
\[
A^+=\xi_2^\dagger\xi_1^\dagger,
\qquad
A^-=\xi_1\xi_2,
\]
where \(\xi_1=\tilde a(x)d+b_1(x)\) and \(\xi_2=\tilde a(x)d+b_2(x)\) [1501.03922]. These satisfy the intertwining relations
\[
A^-h_+=h_-A^-,
\qquad
A^+h_-=h_+A^+.
\]
The associated quasi-Hamiltonian
\[
K=Q^+Q^-+Q^-Q^+=H_{SS}^2-2\lambda H_{SS}+\mu
\]
is a fourth-order differential operator [1501.03922]. For the special case \(\mu=\lambda^2\), it factorizes as
\[
K=(H_{SS}-\lambda)^2
\]
subject to the perfect-square constraint
\[
\xi_2\xi_2^\dagger=\xi_1^\dagger\xi_1
\]
[1501.03922]. This construction explicitly embeds the non-Hermitian Swanson oscillator into a higher-order supersymmetric algebra.

A related one-step extension uses a pseudo-Hermite-polynomial seed \(\phi_m(z)=H_m(iz)e^{z^2/2}\), \(m=2,4,6,\dots\), to generate a rationally extended Hermitian partner \(h^{(-)}\), and then pulls this back by similarity to a non-Hermitian extended Swanson Hamiltonian \(H_{\rm ext}\) [1502.03165]. The extended model has higher-order ladder operators \(K\), \(K^\dagger\), a polynomial Heisenberg algebra, and a spectrum
\[
\sigma(H_{\rm ext})=\{-2m-1\}\cup\{2n+1\mid n\neq m\},
\]
which differs from the equally spaced standard Swanson spectrum by the presence of an extra level [1502.03165].

## 6. Extensions, deformations, and dynamical developments

The Swanson Hamiltonian has served as the seed for numerous extensions. In two dimensions, one construction defines
\[
H_{2D}=H_{S,1}+H_{S,2}
=
\sum_{i=1}^2\Bigl[\omega_i a_i^\dagger a_i+\alpha_i a_i^2+\beta_i(a_i^\dagger)^2+\tfrac{\omega_i}{2}\Bigr],
\]
which, after similarity transformation, becomes the sum of two extended one-dimensional oscillators. Under the resonance condition
\[
n_1(\omega_1^2-4\alpha_1\beta_1)=n_2(\omega_2^2-4\alpha_2\beta_2),
\]
one obtains algebraically independent integrals of motion proving superintegrability [1502.03165].

A different two-dimensional generalization is formulated on a noncommutative plane. There the Hamiltonian \(H_{\nu,\theta}\) is diagonalized by pairs of \(\mathcal D\)-pseudo-bosons \(A_j,B_j\), leading to explicit biorthogonal eigenstates, eigenvalues
\[
E_{n_1,n_2}=\frac{n_1+n_2+\tfrac12}{\cos(2\nu)},
\]
and bi-coherent states with a resolution of the identity on a dense domain [1711.10181]. This indicates that the algebraic Swanson mechanism persists under simultaneous non-Hermiticity and spatial noncommutativity.

Time dependence introduces an additional layer of structure. For the generalized time-dependent Hamiltonian
\[
H(t)=\omega(t)\Bigl(a^\dagger a+\tfrac12\Bigr)+\alpha(t)a^2+\beta(t)a^{\dagger2},
\]
one seeks a time-dependent Dyson map \(\eta(t)\) such that
\[
h(t)=\eta H\eta^{-1}+i\dot\eta\,\eta^{-1}
\]
is Hermitian, equivalently
\[
H^\dagger(t)\rho(t)-\rho(t)H(t)=i\dot\rho(t),
\qquad
\rho(t)=\eta^\dagger(t)\eta(t)
\]
[1606.04807]. With an \(\mathrm{SU}(1,1)\)-type ansatz for \(\eta(t)\), the Hermiticity of \(h(t)\) reduces to coupled nonlinear differential equations for time-dependent parameters, after which the Lewis–Riesenfeld invariant method yields exact solutions [1606.04807].

The dynamical behavior of the Swanson oscillator has also been analyzed semiclassically. For the \(\mathcal{PT}\)-symmetric quadratic Hamiltonian
\[
\hat{\mathcal H}
=
\frac12(\hat p^2+\hat q^2)-\frac{i\delta}{2}(\hat p\hat q+\hat q\hat p),
\]
metriplectic flow equations produce closed-form expressions for the phase-space metric and trajectories [1409.6456]. Although the spectrum remains real and discrete after similarity transformation to a harmonic oscillator with frequency \(\omega=\sqrt{1+\delta^2}\), the classical metric, the real phase-space trajectories, and Gaussian wave packets can diverge in finite time for \(|\delta|\ge1\) [1409.6456]. This is frequently cited as a caution against identifying real spectrum with completely benign dynamics.

Finite-temperature pseudo-Hermitian thermodynamics has likewise been worked out. In the \(\mathcal{PT}\)-unbroken regime, the partition function of the equivalent harmonic oscillator is
\[
Z(\beta)=\frac{1}{2\sinh(\beta\hbar\Omega/2)},
\]
while in the \(\mathcal{PT}\)-broken inverted-oscillator regime one obtains a resonance-sector partition function
\[
Z_r(\beta)=\frac{1}{2i\sin(\beta|\Omega|/2)}
\]
with real internal energy
\[
U_r(\beta)=\tfrac12\hbar|\Omega|\cot(\beta|\Omega|/2)
\]
[2212.13173]. The same framework supports double-time Green’s functions and linear response under periodic perturbations [2212.13173].

More recently, out-of-time-ordered correlators have been computed for the exactly solvable quadratic Swanson Hamiltonian and for Kerr and driven extensions. In the unbroken phase,
\[
[x(t),p(0)]=i\cos(\Omega t),
\qquad
C(t)=\cos^2(\Omega t),
\]
while in the broken phase with \(\Omega=i\Gamma\),
\[
[x(t),p(0)]=i\cosh(\Gamma t),
\qquad
C(t)=\cosh^2(\Gamma t).
\]
The exponential growth in the broken phase is explicitly interpreted as linear instability rather than quantum chaos [2606.04062].

## 7. Generalizations to coupled and fermionic models

The Swanson construction has also been exported to coupled and fermionic settings. A two-mode interacting system is built from two Swanson-like quadratic Hamiltonians plus a \(\mathcal{PT}\)-invariant bilinear interaction,
\[
\mathcal H^{(\mathrm{int})}
=
\alpha(a_1^\dagger a_2+a_2^\dagger a_1)
+\beta(a_1 a_2^\dagger+a_2 a_1^\dagger)
+2\delta(i a_1^\dagger a_2^\dagger-i a_2 a_1),
\]
with the full problem reduced to a \(4\times4\) non-Hermitian matrix [2110.05865]. In the simplifying case \(\beta=-\alpha\), the eigenvalues are available in closed form, and exceptional points arise when the nested square roots collapse. Under the additional conditions
\[
\alpha=\pm\lambda,
\qquad
\gamma^2=\delta^2,
\]
all four roots coalesce at \(E=\omega\), producing fourth-order exceptional points [2110.05865].

A fermionic extension with two fermionic oscillators exhibits exceptional points and a quantum phase transition. The model is described as a fermionic extension of the Swanson scheme with bilinear couplings that do not conserve particle number, and its exceptional points are characterized by coalescing eigenstates and self-orthogonality with respect to the biorthogonal inner product [2401.17189]. The same work emphasizes that the Swanson oscillator remains a standard benchmark for non-Hermitian quantum mechanics [2401.17189].

A three-site fermionic Swanson-like model provides an exactly solvable example of a triple exceptional point. After reduction to a traceless \(3\times3\) block
\[
H(a,A,B,d,u)=
\begin{pmatrix}
-2 & A/a & d\\
a & 1-u & 0\\
B/d & 0 & 1+u
\end{pmatrix},
\]
the secular equation is cubic, and the EP3 locus is given in closed form by
\[
A^{(EP3)}(u)=\frac{(1-u)^3}{2u},
\qquad
B^{(EP3)}(u)=-\frac{(1+u)^3}{2u},
\qquad u\neq 0
\]
[2606.03789]. The same model exhibits a non-empty “corridor of unitarity” near the EP3 when perturbations are scaled by fractional powers of a small parameter [2606.03789]. This suggests a broader role for Swanson-type constructions in higher-order non-Hermitian degeneracies.

A final strand of generalization is interpretive rather than algebraic. The Swanson oscillator has been used as the non-Hermitian structure underlying a quantization of the Schwarzschild horizon area. In that framework, a similarity transformation maps the Swanson Hamiltonian to a scaled harmonic oscillator and yields the area spectrum
\[
A_n=8\pi G\hbar\Bigl(n+\tfrac12\Bigr),
\]
together with a logarithmic entropy correction proportional to \(-\tfrac12\ln A\) [2407.08358]. This does not alter the core quantum-mechanical definition of the model, but it illustrates the portability of the Swanson formalism across mathematical physics contexts.

The Swanson Hamiltonian therefore occupies a distinctive place in non-Hermitian quantum theory: it is simple enough to be solved exactly, but rich enough to support pseudo-Hermitian metrics, biorthogonal spectral theory, supersymmetric factorizations, exceptional-point phenomena of several orders, and a wide range of deformations and applications. Its enduring role in the literature reflects this combination of algebraic tractability and structural depth.

Source: https://www.emergentmind.com/topics/swanson-hamiltonian