Continuum Level Density in Open Quantum Systems

• Continuum level density (CLD) is a quantitative measure that isolates the excess density from interactions beyond the free-gas contribution in unbound systems.
• It links observable scattering phase shifts and transition probabilities to the microscopic structure and boundary conditions of quantum systems.
• Advanced techniques like complex scaling decompose the CLD into resonant and non-resonant contributions, aiding precise modeling in tunneling and pairing phenomena.

Continuum level density (CLD) encodes the modification of the single-particle or many-body density of states due to interactions in open quantum systems where the spectrum is continuous above threshold. The CLD is crucial for describing pairing, collective phenomena, and continuum response in nuclei, atomic and molecular physics, resonant tunneling, and black-hole perturbation theory. By systematically subtracting the contribution of the free system, the CLD isolates the excess density due to a potential or interaction, directly linking observable quantities—such as scattering phase shifts and transition probabilities—to the underlying microscopic structure and boundary conditions of unbound systems.

1. Formal Definition and Mathematical Framework

The continuum level density $\Delta(E)$ is defined as the difference between the interacting and free densities of states at energy $E$: $\Delta(E) = \rho(E) - \rho_0(E) = -\frac{1}{\pi}\, \Im \Tr\left[G(E^+) - G_0(E^+)\right],$ where $G(E^+)$ and $G_0(E^+)$ are the Green's functions of the full and free Hamiltonians, respectively, and the trace is regularized over continuum states (Myo et al., 2020).

In single-channel scattering, $\Delta(E)$ reduces to the derivative of the scattering phase shift,

$\Delta(E) = \frac{1}{\pi} \frac{d\delta(E)}{dE},$

providing an explicit spectral connection between the level density and observable phase shifts (Betan, 2017, Myo et al., 2020, Ogawa et al., 5 May 2026). This representation, first formalized by Beth and Uhlenbeck (1937) and later generalized in nuclear, atomic, and open quantum systems, underpins most practical applications of the CLD.

For partial-wave decompositions in a finite-range potential, the continuum single-particle level density (CSPLD) for each partial wave $(l,j)$ is

$$g_{lj}(\varepsilon) = \frac{1}{\pi} \frac{d\delta_{lj}}{d\varepsilon}, \quad \varepsilon > 0,$$

with the total CLD obtained via a sum over orbital angular momentum and total angular momentum (Betan, 2017, Betan, 2011, Betan et al., 2017): $$g(\varepsilon) = \sum_{l,j} (2j+1)\,g_{lj}(\varepsilon) = \frac{1}{\pi} \sum_{l,j} (2j+1) \frac{d\delta_{lj}}{d\varepsilon}.$$

2. Subtraction of the Free-Gas Contribution

The physical CLD must be finite as the volume $E$0. For a single particle in a spherical box, the total level density contains a divergent free-gas term proportional to $E$1. The subtraction is performed as

$E$2

ensuring that only the effects of the potential remain (Betan, 2017). This excess density encodes resonances, threshold effects, and antibound states, while removing the unphysical background of a free continuum.

3. CLD and Pairing in Open Quantum Systems

The CLD provides a rigorous basis for including continuum correlations in pairing models for weakly bound and unbound nuclei. In the constant-$E$3 pairing Hamiltonian,

$E$4

the replacement of discrete sums over states with integrals involving $E$5 is the key generalization for continuum systems (Betan, 2011, Betan, 2012, Betan et al., 2017).

The BCS gap and particle-number equations become

$E$6

$E$7

where $E$8 indexes bound states and the continuum is described by $E$9 and corresponding occupation amplitudes. In the exact Richardson or generalized Cooper pair solution, integrals over $\Delta(E) = \rho(E) - \rho_0(E) = -\frac{1}{\pi}\, \Im \Tr\left[G(E^+) - G_0(E^+)\right],$0 appear explicitly in the pair energy equations. This formulation enables the exact solution of the pairing problem in systems with no bound single-particle levels, as in Borromean nuclei, and provides reliable access to all eigenstates, including those near the threshold (Betan, 2017, Betan, 2012).

Numerical studies confirm that resonant structures in $\Delta(E) = \rho(E) - \rho_0(E) = -\frac{1}{\pi}\, \Im \Tr\left[G(E^+) - G_0(E^+)\right],$1 strongly modulate the occupation of partial-wave configurations in Borromean systems such as $\Delta(E) = \rho(E) - \rho_0(E) = -\frac{1}{\pi}\, \Im \Tr\left[G(E^+) - G_0(E^+)\right],$2He (dominated by $\Delta(E) = \rho(E) - \rho_0(E) = -\frac{1}{\pi}\, \Im \Tr\left[G(E^+) - G_0(E^+)\right],$3) and $\Delta(E) = \rho(E) - \rho_0(E) = -\frac{1}{\pi}\, \Im \Tr\left[G(E^+) - G_0(E^+)\right],$4Li (dominated by $\Delta(E) = \rho(E) - \rho_0(E) = -\frac{1}{\pi}\, \Im \Tr\left[G(E^+) - G_0(E^+)\right],$5 at threshold), and that the inclusion of CLD improves binding energy predictions and stabilizes the pairing gap near the drip line (Betan, 2017, Betan et al., 2017). In BCS-type approaches, the continuum contribution can account for a sizable portion—$\Delta(E) = \rho(E) - \rho_0(E) = -\frac{1}{\pi}\, \Im \Tr\left[G(E^+) - G_0(E^+)\right],$6–$\Delta(E) = \rho(E) - \rho_0(E) = -\frac{1}{\pi}\, \Im \Tr\left[G(E^+) - G_0(E^+)\right],$7\%—of the total pairing gap, even if the average continuum occupation is small (Betan, 2011).

4. Complex Scaling, Resonant and Non-Resonant CLD

Complex scaling (CSM) is a central tool for computing the CLD in systems with resonances. By rotating the coordinates as $\Delta(E) = \rho(E) - \rho_0(E) = -\frac{1}{\pi}\, \Im \Tr\left[G(E^+) - G_0(E^+)\right],$8, the spectrum of the Hamiltonian is deformed: resonances become discrete eigenvalues with complex energies $\Delta(E) = \rho(E) - \rho_0(E) = -\frac{1}{\pi}\, \Im \Tr\left[G(E^+) - G_0(E^+)\right],$9, while the continuum is rotated by $G(E^+)$0 into the complex plane (Myo et al., 2020, 2002.03874, Ogawa et al., 5 May 2026).

The CLD can then be written as a sum of Lorentzians from resonant poles and a smooth background associated with non-resonant (continuum) states: $G(E^+)$1 where $G(E^+)$2 is a smoothed background density intrinsically linked to the non-resonant spectrum (Myo et al., 2020). CSM enables a separation of bound, resonant, and background continuum contributions within a unified spectral framework.

In quantum tunneling, the CLD can be analytically continued to complex energies, yielding both real and imaginary parts that carry physical interpretation: the real part is proportional to the Eisenbud–Wigner delay time, while the imaginary part encodes instanton-like tunneling trajectories in complex time, reflecting the energy dependence of transmission amplitudes (2002.03874).

5. Applications in Nuclear Structure, Tunneling, and Black-Hole Physics

Nuclear Structure and Reactions:

In nuclear systems near the drip line, the CLD method affords an efficient and physically transparent way to incorporate continuum effects into large-scale pairing calculations, yielding accurate predictions for observables such as two-neutron separation energies, binding energies, and pairing gaps across isotopic chains (Betan et al., 2017). The method systematically distinguishes between resonant and non-resonant continuum contributions, modulating the effective pairing in open systems and improving agreement with experimental data (Betan, 2012, Betan, 2011). In Borromean nuclei, a CLD-based constant-$G(E^+)$3 model reproduces occupation probabilities, partial-wave mixing, and total energies with accuracy comparable to three-body or ab initio methods (Betan, 2017).

Resonant Quantum Tunneling:

The complex-extended CLD provides a quantitative connection between the quantum density of states, transmission phases, and under-barrier dynamics in finite-range tunneling problems, and encodes non-analyticities at the stationary points of potentials (e.g., logarithmic divergences at barrier maxima) (2002.03874).

Black-Hole Perturbations:

The CLD, computed via complex scaling, captures the continuum response of black holes beyond isolated quasinormal modes. In Schwarzschild–de Sitter spacetimes, both the resonant (QNM pole) and continuum (rotated cut) sectors are unified through the CLD, enabling detailed analysis of linear response, late-time tails, and horizon-induced scattering (Ogawa et al., 5 May 2026). The CLD is closely tied to the phase shift and thus to greybody factors in Hawking radiation, providing a complete spectral description for open gravitational systems.

6. Computational and Methodological Considerations

The calculation of the CLD involves: (i) numerically generating scattering phase shifts $G(E^+)$4, (ii) computing their energy derivatives to obtain $G(E^+)$5, and (iii) integrating the resulting density in the relevant pairing, response, or transport equations. Alternatively, in the CSM, one solves the non-Hermitian Hamiltonian in an $G(E^+)$6 basis, obtaining a set of complex eigenvalues whose distribution reconstructs the CLD via Lorentzian summation and background subtraction (Myo et al., 2020).

CSPLD-based methods replace the need for large discretized oscillator or box bases in open-shell and continuum pairing problems. This approach achieves computational economy and removes model-space artifacts, particularly in Richardson-type solutions, where the correlated energy is model-space independent once the interaction strength is fixed (Betan, 2012, Betan, 2011).

7. Significance, Limitations, and Outlook

The CLD framework enables the consistent inclusion of both resonant and non-resonant continuum physics in models of open quantum systems with few or no bound states, capturing essential features of pairing, scattering, and collective excitation spectra across nuclear, molecular, condensed-matter, and gravitational settings. It provides a spectroscopic and response-theoretic bridge between traditional bound-state approaches and fully unbound, non-Hermitian, or strongly correlated regimes.

Limitations arise when employing constant-$G(E^+)$7 interactions, which can miss fine state-dependent effects; the smoothing of $G(E^+)$8 near threshold and the choice of cutoff parameters introduce some model dependence. Advanced schemes based on state-dependent interactions or further generalizations of the spectral expansion are active avenues for research (Betan et al., 2017). Nevertheless, the CLD method remains an essential tool in the precise theoretical description of systems where continuum and resonance phenomena are dominant.