---
title: Broeckhove-Deumens Theorem in Shell Models
url: https://www.emergentmind.com/topics/broeckhove-deumens-theorem
type: topic
---

# Broeckhove-Deumens Theorem in Shell Models

Searching arXiv for the specified paper and closely related background on non-orthogonal shell-model / symmetry-restored methods.
The Broeckhove–Deumens theorem, as stated in recent shell-model work, is an existence theorem about the discretization of continuous non-orthogonal manifolds in Hilbert space. For a separable Hilbert space $\mathscr L^2$ and a closed subspace $\mathscr H = \overline{\mathrm{span}\,\Gamma} \subseteq \mathscr L^2$, where $\Gamma$ is a continuous, dense set of non-orthogonal states, the theorem asserts the existence of a discrete countable non-orthogonal subset $\Gamma_0 \subset \Gamma$ such that $\mathscr H = \overline{\mathrm{span}\,\Gamma_0}$ [2507.09073]. In nuclear-structure applications, this statement underpins the use of discrete sets of symmetry-restored non-orthogonal Slater determinants in place of very large orthonormal shell-model configuration spaces. The 2025 work “Exact solutions of the nuclear shell-model secular problem: Discrete Non-Orthogonal Shell Model within a Variation After Projection approach” formulates this connection explicitly and presents numerical realizations in sd-shell nuclei, in $^{48}$Cr, and in $^{78}$Ni [2507.09073].

## 1. Formal statement and Hilbert-space setting

The paper gives the theorem in the following form: for a separable Hilbert space $\mathscr L^2$ and a closed subspace $\mathscr H = \overline{\mathrm{span}\,\Gamma}$, with $\Gamma$ a continuous, dense set of non-orthogonal states, there exists a discrete countable non-orthogonal subset $\Gamma_0 \subset \Gamma$ such that
$$
\mathscr H = \overline{\mathrm{span}\,\Gamma_0}.
$$
Equivalently, any $\psi \in \mathscr H$ can be represented as
$$
\psi = \sum_k c_k \phi_k,
$$
with $\{\phi_k\} \subset \Gamma_0$ and convergence in the Hilbert-space norm [2507.09073].

In the shell-model context, the relevant space is the valence-space many-body Hilbert space $\mathscr H_{SM}$. The standard orthonormal basis is the set of spherical Slater determinants, identified in the paper with the usual $M$-scheme basis. Since this space is finite dimensional in the applications considered, the countable subset implied by the theorem becomes finite in practice. The significance of the theorem in this setting is therefore precise: it guarantees the existence of a finite family of non-orthogonal Slater determinants spanning the same space as the standard shell-model basis [2507.09073].

The paper also emphasizes the theorem’s limitations. It does not specify the nature of the continuous family $\Gamma$, does not characterize the discrete subset $\Gamma_0$, and does not provide a constructive procedure for obtaining it. Its content is existential rather than algorithmic. A central contribution of the 2025 work is therefore not a new theorem, but an explicit numerical construction showing how such a subset can be realized for realistic shell-model Hamiltonians [2507.09073].

## 2. Shell-model secular problem in a non-orthogonal representation

The shell-model Hamiltonian is written as a standard one- plus two-body operator in a valence space of spherical harmonic-oscillator orbitals,
$$
\hat H = \sum_{ij} t_{ij} c^\dagger_i c_j + \frac 1 4 \sum_{ijkl} V_{ijkl} \, c^\dagger_i c^\dagger_j c_l c_k ,
$$
with antisymmetry conditions
$$
V_{ijkl} = -V_{jikl} = -V_{ijlk} = V_{jilk}.
$$
Here the creation and annihilation operators obey the canonical anticommutation relations, and the indices label spherical single-particle orbitals in the chosen valence shell [2507.09073].

Rather than expanding eigenstates in an orthonormal configuration basis, the paper uses non-orthogonal Slater determinants $\{|\phi_q\rangle\}$, optionally supplemented by many-particle–many-hole excitations. Two frameworks are distinguished: projection after variation (PAV), based on constrained Hartree–Fock plus $NpNh$ excitations, and variation after projection (VAP), based on Thouless-parameterized determinants [2507.09073].

States of good $J^\pi$ are constructed as
$$
|\psi^{\pi JM}_n\rangle = \sum_{q,K} C^{\pi J}_{n,qK}\, \mathcal P^J_{MK} P^\pi |\phi_q\rangle,
$$
and in the extended PAV treatment as
$$
\begin{aligned}
|\psi^{\pi JM}_n\rangle &= \sum_{q,K} C^{\pi J}_{n,qK}\, \mathcal P^J_{MK} P^\pi |\phi_q\rangle \\
&\quad + \sum_{q,K,N} C^{\pi J}_{n,qK,N}\, \mathcal P^J_{MK} P^\pi |\phi_q(N\mathrm p N\mathrm h)\rangle.
\end{aligned}
$$
Here $P^\pi$ is the parity projector and $\mathcal P^J_{MK}$ is the angular-momentum projector from intrinsic projection $K$ to laboratory projection $M$ [2507.09073].

The basic objects are the norm and Hamiltonian kernels,
$$
\mathcal N^{\pi J}_{qK, q'K'} = \langle \phi_q | \mathcal P^J_{KK'} P^\pi | \phi_{q'} \rangle,
\qquad
\mathcal H^{\pi J}_{qK, q'K'} = \langle \phi_q | \hat H \, \mathcal P^J_{KK'} P^\pi | \phi_{q'} \rangle,
$$
which define the Hill–Wheeler–Griffin generalized eigenvalue problem
$$
\mathcal H^{\pi J} \, C^{\pi J}_n = E^{\pi J}_n \, \mathcal N^{\pi J} \, C^{\pi J}_n.
$$
Written in standard linear-algebra form, this is the generalized eigenproblem $Hc = ESc$ in a non-orthogonal basis [2507.09073].

For the VAP construction, the determinants are generated through the Thouless parameterization
$$
|\phi_q\rangle = \mathcal N_0
\exp\!\left(\sum_{ij} Z^{(q)}_{ij} a^\dagger_i a_j\right) |\phi^{(q)}_0\rangle,
$$
with $Z^{(q)}_{ij}$ a skew-symmetric complex matrix and $|\phi^{(q)}_0\rangle$ a fixed reference Slater determinant for that branch. This parameterization provides the variational degrees of freedom used to optimize projected states directly [2507.09073].

## 3. Numerical realization in the discrete non-orthogonal shell model

The paper interprets “realization” of the theorem as the explicit construction, for a given realistic shell-model Hamiltonian, of a finite discrete set $\Gamma_0$ of non-orthogonal Slater determinants such that the generalized eigenvalue problem in that space reproduces the exact shell-model eigenvalues and eigenstates for low-lying states of interest [2507.09073].

In the PAV implementation, the starting point is a continuous generator-coordinate set of triaxial Hartree–Fock states labeled by quadrupole deformation coordinates $(\beta,\gamma)$. A basis-selection procedure identified as the Caurier basis-selection technique is then used to extract a minimal discrete subset. The procedure begins with an HF minimum, selects subsequent states by minimizing the lowest eigenvalue of the generalized eigenproblem in the enlarged space, and continues iteratively until convergence. After this HF stage, the basis is enlarged by selected $NpNh$ excitations that further improve the variational energy [2507.09073].

This yields the PAV realization
$$
\Gamma_0^{(\text{PAV})} = \{ |\phi_q\rangle, |\phi_q(N\mathrm p N\mathrm h)\rangle\}.
$$
The importance of this construction lies in its diagnostic role. The paper shows that a manifold restricted to $(\beta,\gamma)$-constrained HF determinants is generally insufficient: it spans a proper subspace of the shell-model space and leaves residual correlation energy unrecovered even when the number of HF generator points becomes very large. By contrast, the enlarged HF plus $NpNh$ set can reach the exact shell-model ground-state energies in the test cases considered [2507.09073].

The VAP implementation is more compact. It constructs generic non-orthogonal determinants through direct minimization of the projected energy, thereby generating states that effectively encode the physics of many $NpNh$ configurations. The paper explicitly states that VAP is “strictly equivalent to performing $NpNh$ excitations” in the full space, and presents VAP as the practical device that finds a discrete subset satisfying the theorem’s completeness requirement in concrete shell-model calculations [2507.09073].

This suggests a useful conceptual distinction between theorem and realization. The theorem asserts only the existence of a discrete spanning subset. The discrete non-orthogonal shell model provides an explicit algorithmic route to such a subset, with the quality of the realization judged by whether the resulting Ritz eigenvalues coincide with exact shell-model energies to numerical precision.

## 4. Variation after projection as the constructive mechanism

In the VAP framework, the projected energy functional is
$$
E^{\pi J}_n =
\frac{\langle \psi^{\pi JM}_n|\hat H|\psi^{\pi JM}_n\rangle}
     {\langle \psi^{\pi JM}_n|\psi^{\pi JM}_n\rangle},
\qquad
\langle \psi^{\pi JM}_n|\psi^{\pi JM}_n\rangle = 1,
$$
with
$$
|\psi^{\pi JM}_n\rangle = \sum_{q,K} C^{\pi J}_{n,qK} \mathcal P^J_{MK} P^\pi |\phi_q\rangle.
$$
This is the Ritz variational principle applied to symmetry-restored states [2507.09073].

Variation with respect to the mixing amplitudes $C^{\pi J}_{n,qK}$ produces the generalized eigenvalue problem. Variation with respect to the intrinsic Slater determinants, using the Thouless parameterization, yields the projected Brillouin condition
$$
\frac{\displaystyle\sum_{qK,q'K'}
 C^{*\pi J}_{n,qK}C^{\pi J}_{n,q'K'}
 \langle \phi_q | a^\dagger_j a_i (\hat H - E^{\pi J}_n) \mathcal P^J_{KK'}P^\pi |\phi_{q'}\rangle}
{\displaystyle\sum_{qK,q'K'}
C^{*\pi J}_{n,qK}C^{\pi J}_{n,q'K'}
\langle \phi_q | \mathcal P^J_{KK'}P^\pi |\phi_{q'}\rangle } = 0,
\qquad \forall q,i,j.
$$
The paper identifies this as the resonating Hartree–Fock condition of Fukutome, generalized to symmetry projection [2507.09073].

The implementation uses a hybrid optimization strategy. The first $q_0$ Slater determinants are kept fixed once optimized, the full Brillouin condition is applied only to the most recently added determinants, and a quasi-Newton algorithm of L-BFGS type is used to minimize the projected energy [2507.09073]. The practical consequence is that VAP searches directly in the manifold of symmetry-restored non-orthogonal determinants rather than relying on an a priori deformation mesh plus explicit combinatorial excitation generation.

Within the logic of the theorem, this role is decisive. PAV with deformation coordinates alone fails because the generator manifold is not complete for the shell-model space under study. PAV plus $NpNh$ corrections can restore completeness, but the enlargement can become combinatorially large. VAP bypasses this by optimizing the determinants in the projected space itself, so that each determinant functions as a highly compressed carrier of many-shell-model configurations [2507.09073].

## 5. Empirical demonstrations in sd-shell nuclei, $^{48}$Cr, and $^{78}$Ni

The paper tests the theorem’s realization in several nuclear systems. In sd-shell nuclei with the USDB interaction, it examines $^{20}$Ne, $^{24}$Mg, $^{28}$Si, and $^{26}$Al. For $^{24}$Mg, DNO-SM(PAV) with only $(\beta,\gamma)$-constrained HF determinants gives a ground-state energy of $-86.73278$ MeV, whereas the exact shell-model result is $-87.10445$ MeV, a difference of approximately $0.37$ MeV. The paper reports that this discrepancy does not disappear even when approximately 6000 HF points are included across the $(\beta,\gamma)$ plane, indicating that the pure deformation manifold spans only a proper subspace of the full shell-model space [2507.09073].

When $NpNh$ excitations are added in PAV, the agreement becomes essentially exact. For the same nucleus, the enlarged calculation yields $-87.10428$ MeV versus the exact $-87.10445$ MeV. Similar exact agreement is reported for $^{20}$Ne and $^{28}$Si. The dominant excitations are described as $N = 2,4,6,8$ particle–hole pairs coupled to $K = 0$, that is, pair-condensate-like configurations [2507.09073].

The VAP results are more compressed still. For $^{20}$Ne, $^{24}$Mg, $^{28}$Si, and $^{26}$Al, the differences between DNO-SM(VAP) and exact shell-model energies are at the level of a few $10^{-4}$ MeV. Two representative basis-size comparisons are highlighted: $^{28}$Si is reproduced with 45 VAP Slater determinants versus 93,710 standard shell-model configurations, and $^{24}$Mg with 16 determinants versus 28,503 shell-model configurations [2507.09073]. The paper presents these calculations as numerical evidence that a finite set of projected non-orthogonal Slater determinants can span the full valence space for the low-lying states considered.

In $^{48}$Cr with the KB3 interaction, the VAP approach is applied state by state along the yrast band. The resulting energies for $0^+_1,2^+_1,\dots,16^+_1$ agree with exact shell-model values to within $\leq 0.013$ MeV, often exactly within the quoted digits, with basis sizes ranging from 50 determinants for $0^+_1$ down to 12 determinants for $16^+_1$ [2507.09073]. The $\gamma$-ray transitions,
$$
E_\gamma(J) = E(J) - E(J-2),
$$
reproduce the backbending pattern exactly, including the impact of proton–neutron pairing correlations. In the paper’s interpretation, this shows that number-conserving non-orthogonal Slater determinants under VAP fully capture pairing correlations in this rotational band [2507.09073].

For $^{78}$Ni in the pf-sdg space with the PFSDG-U interaction, the shell-model dimension is reported as approximately $2\times 10^{11}$ in the $M$-scheme, beyond direct full diagonalization. Truncated shell-model calculations up to 10p–10h yield a ground-state energy of $-372.71668$ MeV and an extrapolated value of $-372.72850$ MeV. DNO-SM(VAP), starting from a spherical HF reference and adding optimized non-orthogonal determinants, converges to $E_{gs} = -372.73275$ MeV, lower than the extrapolated shell-model energy and presented as essentially the best variational bound in the full valence space [2507.09073]. The paper further reports exponential-like convergence with the number of determinants and saturation at the quoted value.

## 6. Conceptual scope, limitations, and broader connections

The theorem’s conceptual content is straightforward in functional-analytic terms: a dense continuous manifold of non-orthogonal vectors admits a discrete countable subset with the same closed span. In finite-dimensional shell-model spaces, this reduces to the statement that the full valence-space Hilbert space can be spanned by a finite family of non-orthogonal Slater determinants, provided that the family is sufficiently rich [2507.09073].

The limitations arise at the level of realization rather than theorem. The paper notes that the completeness of a particular continuous generator manifold is not automatic. For the triaxial deformation manifold $\Gamma(\beta,\gamma)$, separability and completeness are not trivial to prove in general, and the explicit calculations show that the subspace generated by $(\beta,\gamma)$-HF states is not the full shell-model space. Persistent energy differences, even with a dense deformation mesh, provide the practical evidence for this restriction [2507.09073]. A common misconception is therefore that any continuous generator-coordinate manifold is automatically complete; the paper explicitly argues against this.

The non-orthogonal basis must also be numerically well behaved. Because overlap matrices $\mathcal N^{\pi J}$ can become ill-conditioned when determinants are nearly linearly dependent, basis selection and stabilization are required. The use of Caurier’s selection method in PAV and controlled determinant addition in VAP addresses this issue operationally [2507.09073].

The results are also scoped to low-lying states. The paper presents exact or near-exact reproduction for ground states and yrast bands, and describes this as a realization of the theorem for “low-lying states of interest” [2507.09073]. This suggests that the practical algorithm is optimized for the physically relevant low-energy sector, even though the theorem itself is formulated at the level of the entire closed span.

The broader theoretical placement given in the paper links the theorem to generator coordinate methods, the Hill–Wheeler–Griffin framework, Peierls–Yoccoz symmetry restoration, non-orthogonal configuration interaction, resonating Hartree–Fock, and VAMPIR-like strategies [2507.09073]. In this perspective, the Broeckhove–Deumens theorem serves as the mathematical foundation for replacing continuous generator-coordinate superpositions by discrete non-orthogonal expansions without loss of completeness, while the discrete non-orthogonal shell model shows that this replacement can be made numerically exact for realistic shell-model Hamiltonians.

A plausible implication is that the theorem’s practical value is greatest not as a standalone formal result, but as a justification for aggressively compressed, symmetry-restored non-orthogonal representations of shell-model eigenstates. In the examples studied, the compression is extreme, yet the resulting energies coincide with exact diagonalization where available, or improve upon the best truncated-extrapolated shell-model estimates in spaces near current computational limits [2507.09073].

Source: https://www.emergentmind.com/topics/broeckhove-deumens-theorem