---
title: Monopole Interaction Truncation
url: https://www.emergentmind.com/topics/monopole-interaction-based-truncation
type: topic
---

# Monopole Interaction Truncation

Monopole-interaction-based truncation is a configuration-interaction shell-model strategy in which the monopole part of the effective Hamiltonian is used to rank occupation-defined subspaces and to retain only the energetically favored sectors of Hilbert space. In the modern shell-model literature, the method rests on the decomposition \(H_A=H_m+H_M\), where \(H_m\) is the monopole, or mean-field, part and \(H_M\) is the multipole residual interaction; low-lying states are then approximated by configurations with the lowest monopole energies or centroids, while angular momentum conservation and rotational symmetry are preserved by construction [2507.11796][2511.03161].

## 1. Formal definition in the shell-model Hamiltonian

In the shell model with a frozen core, the effective Hamiltonian is written as
\[
H_A=\sum_{\alpha} \epsilon_{\alpha} a^{\dagger}_{\alpha} a_{\alpha} + \frac{1}{4} \sum_{\alpha\beta\gamma\delta JT}
\langle j_{\alpha} j_{\beta} | V | j_{\gamma} j_{\delta} \rangle_{JT}
A^{\dagger}_{JT; j_{\alpha} j_{\beta}}A_{JT; j_{\delta} j_{\gamma}},
\]
with \(\alpha\) labeling single-particle orbitals, \(\epsilon_\alpha\) the single-particle energies, and \(\langle j_\alpha j_\beta|V|j_\gamma j_\delta\rangle_{JT}\) the antisymmetrized two-body matrix elements. Following the standard decomposition used in shell-model monopole analyses, one writes
\[
H_A=H_m+H_M,
\]
where \(H_m\) is the monopole Hamiltonian and \(H_M\) the multipole Hamiltonian [2507.11796][1207.5056].

The monopole matrix elements are the angular-momentum-weighted averages of the diagonal two-body matrix elements:
\[
V_{m;\alpha\beta}=
\frac{\displaystyle\sum_J (2J+1)\,\langle j_{\alpha} j_{\beta} | V | j_{\alpha} j_{\beta} \rangle_{J,T}}
{\displaystyle\sum_J (2J+1)}.
\]
The monopole expectation value in a state \(|\Psi\rangle\) is
\[
E_m=\langle \Psi|H_m|\Psi\rangle
=\sum_{\alpha}\epsilon_{\alpha}\langle \hat N_{\alpha}\rangle
+\sum_{\alpha\le\beta}V_{m;\alpha\beta}
\left\langle
\frac{\hat N_{\alpha}(\hat N_{\beta}-\delta_{\alpha\beta})}{1+\delta_{\alpha\beta}}
\right\rangle.
\]
Because this term depends only on orbital occupancies and not on detailed angular-momentum coupling, it defines the effective mean field, controls shell evolution, and sets the gross ordering of configurations; the multipole part is then responsible for pairing, quadrupole correlations, seniority splittings, and other spectroscopic fine structure [1207.5056][1510.05340].

This separation is the formal basis of monopole-interaction-based truncation. If \(H_m\) dominates the smooth part of the spectrum, then configurations with low monopole energy are natural candidates for retention, while high-monopole-energy configurations can be discarded with comparatively small impact on low-lying states. This suggests a mean-field-guided importance ordering rather than an a priori \(n\hbar\omega\) or particle-hole counting.

## 2. Partitions, centroids, and occupation-space ordering

The practical object in monopole truncation is the partition, or configuration, defined by a fixed set of orbital occupation numbers. For a partition \(P\), with occupations \(N_{P;\alpha}\), the monopole energy is
\[
E_m^P
= \sum_{\alpha} \epsilon_{\alpha} N_{P;\alpha}
+ \sum_{\alpha \le \beta} V_{m;\alpha\beta}
\frac{N_{P;\alpha}(N_{P;\beta} - \delta_{\alpha\beta})}{1+\delta_{\alpha\beta}}.
\]
Within a given partition, all basis states have the same monopole energy, because the monopole term depends only on occupations, not on the detailed coupling of single-particle angular momenta [2507.11796].

The same logic appears in centroid-based formulations. If a configuration subspace is labeled by occupations \(\alpha=\{n_a\}\), its dimension is
\[
d_\alpha=\prod_a \binom{N_a}{n_a},
\qquad N_a=2j_a+1,
\]
and its energy centroid is
\[
E_\alpha=\sum_a n_a \epsilon_a
+\frac{1}{2}\sum_{ab} n_a(n_b-\delta_{ab})U(ab),
\]
where \(U(ab)\) is the monopole potential obtained from the angular-momentum-weighted average of the diagonal two-body matrix elements. These centroids are traces over occupation-defined subspaces and can be computed without constructing the many-body Hamiltonian matrix [2511.03161].

This centroid formulation makes the ordering principle explicit. One ranks configurations by \(E_\alpha\) or \(E_m^P\), identifies the minimum monopole energy,
\[
E_{\text{mono\_min}}=\min_P E_m^P,
\]
and then retains only the subset of partitions lying sufficiently close to it. The resulting truncation is not arbitrary: it is tied to the average interaction energy of occupation patterns. In this sense, monopole-interaction-based truncation is an occupation-space implementation of a mean-field energy hierarchy.

A related development is the TRACER code, which computes monopole potentials, configuration dimensions, and configuration centroids rapidly from shell-model interaction files, and then uses those centroids to optimize energy-based truncation schemes for codes such as BIGSTICK [2511.03161].

## 3. Truncation schemes and algorithmic realizations

The simplest scheme is the sharp cutoff at the partition level. A partition is retained if
\[
E_m^P-E_{\text{mono\_min}} \le E_{\rm cutoff}.
\]
All basis states belonging to the retained partitions are kept; all others are discarded. As \(E_{\rm cutoff}\) increases, the truncated space approaches the full space monotonically [2507.11796].

A second scheme is the distribution-type cutoff truncation. Instead of a hard boundary, one defines
\[
X(E_{\rm mono})=
\frac{1}{1+\exp\left[\frac{E_{\rm mono}-E_{\rm mono\_min}-E_{\rm cutoff}}{a}\right]},
\]
with \(a\) a smoothing parameter. For each partition, one generates a random number \(r\in[0,1]\); if \(r<X(E_m^P)\), that partition is retained. Low-monopole-energy partitions therefore have \(X\approx 1\), whereas higher-energy partitions are sampled only with some probability [2507.11796].

A third variant applies the same distribution function inside each retained partition at the \(J\)-scheme level. If a partition contains \(nJT\) many \(J\)-scheme basis states, the allowed number is set to
\[
n_{\text{allowed}} = X(E_m^P)\times nJT.
\]
Low-energy partitions keep nearly all their basis states, while high-energy partitions are heavily thinned. This is a more aggressive reduction that directly targets the projected basis [2507.11796].

A different but closely related construction is the Approximately Cutoff-Energy-dependent truncation. Here one seeks a weight-based truncation of the type
\[
W_\alpha=\sum_a n_{a,\alpha}w_a \le W_{\rm max},
\]
where the orbital weights \(w_a\) are integers chosen so that the weight-based truncated space approximates a centroid cutoff \(E_\alpha\le E_{\rm cut}\). The optimization is driven by the fitness function
\[
\Phi(E_{\rm cut},\vec{w},W_{\rm max},{\cal F})
=
\sum_{\alpha \in {\cal P}(\vec{w}, W_{\rm max})}
d_\alpha (E_\alpha - E_{\rm cut}),
\]
with \({\cal F}\) the enveloping configuration space and \({\cal P}\) the subset satisfying the weight criterion. TRACER performs a Monte Carlo search over \(\vec w\) and scans \(W_{\rm max}\) to minimize \(\Phi\), then writes a new single-particle-space file with the optimized weights for use in BIGSTICK [2511.03161].

In the KTH version of NuShellX, monopole truncation is implemented by generating partitions, computing \(E_m^P\), selecting important partitions, constructing the truncated \(M\)-scheme basis, and then restoring good total \(J\) by projection. Because truncation is imposed at the level of partitions and projected basis states rather than by arbitrary removal of \(M\)-scheme determinants, angular momentum conservation and rotational symmetry are preserved [2507.11796].

## 4. Physical basis: shell evolution, ESPEs, and configuration hierarchies

The physical justification for the method comes from the empirical hierarchy between monopole and multipole effects. In the tin calculations based on a CD-Bonn-derived interaction, the monopole contribution dominates the absolute scale of energies and their smooth evolution with neutron number, whereas the multipole contribution is relatively small and peaks near mid-shell. The same work shows that light tin isotopes are dominated by \(0g_{7/2}\)-\(1d_{5/2}\) configurations, that the \(0h_{11/2}\) orbit becomes half-filled around \(N=22\)–\(24\), corresponding to \(^{122}\)–\(^{124}\)Sn, and that a truncation to the \(0g_{7/2}1d_{5/2}\) space is natural near \(N=50\), with explicit \(h_{11/2}\) inclusion required toward mid-shell [1207.5056].

This occupation-driven hierarchy is often expressed through effective single-particle energies. Standard practice is to write
\[
\epsilon_\alpha^{\rm eff}
=
\varepsilon_\alpha + \sum_\beta V_{m;\alpha\beta}\,\langle \hat N_\beta\rangle,
\]
so that monopole filling shifts shell gaps and determines which orbitals are near the Fermi surface. Although several of the cited works use this logic without always printing the formula explicitly, the practical consequence is consistent: orbits with low \(\epsilon_\alpha^{\rm eff}\) and sizable occupancy define the active space, while orbits remaining high in energy can be removed or treated approximately through effective monopole shifts [1207.5056][1510.05340].

The PMMU framework in the \(pf_{5/2}g_{9/2}\) shell makes the same point in a broader region. Its Hamiltonian,
\[
H=H_0+H_{PM}+H_m,
\]
combines single-particle energies, pairing-plus-multipole terms, and a monopole interaction obtained from the monopole-based universal force and empirical fits. In Cu isotopes, the proton \(f_{5/2}\) effective single-particle energy is pulled down relative to \(p_{3/2}\) as neutrons fill \(g_{9/2}\), and for \(N>46\), \(f_{5/2}\) becomes the lowest proton orbit. In Ge at \(N=40\), attractive \(T=0\) monopole terms increase the \(p_{1/2}\)–\(g_{9/2}\) neutron gap and suppress \(g_{9/2}\) excitations in the \(0_2^+\) state, whereas yrast states require those excitations. This suggests a state-dependent truncation strategy in which monopole-favored intruders are retained for collective yrast structures but can be more strongly restricted for near-closed-shell excitations [1510.05340].

The same lesson appears in the Ni region through the monopole tensor interaction. In \(^{64}\)Ni, the Monte Carlo shell model with the A3DA-m interaction places the prolate \(0^+_4\) state at \(3463.1\) keV, with a collective \(2^+\) state at about \(3749\) keV. The prolate minimum is driven by neutron \(g_{9/2}\) and \(d_{5/2}\) occupancy together with correlated proton cross-shell excitations; when the proton-neutron monopole part is frozen, the prolate minimum nearly vanishes. A plausible implication is that monopole truncation is not equivalent to simply minimizing basis size: it must retain precisely those intruder occupations that receive large monopole-tensor energy gains [2008.04857].

Cross-shell applications based on \(V_{MU}\) extend this logic to truncated multi-\(\hbar\omega\) spaces. In boron, carbon, nitrogen, and oxygen isotopes, a single Hamiltonian is used in \(0\hbar\omega\), \((0+2)\hbar\omega\), and \((0\text{–}3)\hbar\omega\) spaces without deriving a separate effective interaction for each truncation, and the inclusion of \(2\hbar\omega\) excitations is important for drip lines, spectra, and electromagnetic properties. In the PMMU program for the \(pf\) and \(pf_{5/2}g_{9/2}\) regions, limited excitations out of deeply bound orbitals are introduced only where needed computationally, while the monopole sector preserves the shell evolution and effective single-particle-energy hierarchy that justifies the restricted valence spaces [1209.5587][1312.6191].

## 5. Quantitative performance in Sn, Xe, and Pb

Recent large-scale tests establish the method as a practical truncation strategy rather than only a qualitative guide. In the Pb, Sn, and Xe studies, convergence is assessed by increasing \(E_{\rm cutoff}\) and monitoring energies, excitation energies, and selected \(B(E2)\) values [2507.11796].

| System | Full-space scale | Representative truncated result |
|---|---:|---|
| \(^{202}\)Pb | \(M^\pi=0^+\) dimension \(411{,}184\); full spaces up to \(\sim 3.4\times10^9\) | At \(E_{\rm cutoff}=10\) MeV, about \(65\%\) of the total basis gives \(0_1^+=0.147\) MeV versus full-space \(0.056\) MeV |
| \(^{106}\)Sn | large full space in the \(0g_{7/2},1d_{5/2},1d_{3/2},2s_{1/2},0h_{11/2}\) space | At \(E_{\rm cutoff}=10\) MeV, about \(46\%\) of the full dimension gives ground-state energy \(-6.505\) MeV versus full \(-6.768\) MeV |
| \(^{108}\)Xe | ground-state \(M\)-scheme dimension \(\sim 7.4\times10^7\) | At \(E_{\rm cutoff}=8\) MeV, about \(50\%\) of the model space gives \(E(0_1^+)=-15.426\), \(E(2_1^+)=-14.991\), \(E(4_1^+)=-14.323\) MeV |
| \(^{110}\)Xe | ground-state \(M\)-scheme dimension \(\sim 1.7\times10^9\) | At \(E_{\rm cutoff}=10\) MeV, about \(37\%\) of the full dimension gives \(E(0_1^+)=-19.391\), \(E(2_1^+)=-19.098\), \(E(4_1^+)=-18.536\) MeV |

For \(^{202}\)Pb, the distribution-type cutoff with \(a=0.5\) MeV converges faster than the sharp cutoff at low cutoff energies. At \(E_{\rm cutoff}=6\) MeV, the sharp cutoff gives a \(0_1^+\) energy of about \(0.617\) MeV, while the distribution-type cutoff gives a minimum \(0_1^+\) energy of about \(0.518\) MeV. The basis-level distribution scheme improves convergence further for excited states, and combining monopole and seniority truncation yields energies closer to the full-space result at a given basis fraction [2507.11796].

For \(^{106}\)Sn, the full-space calculation gives \(B(E2;2_1^+\rightarrow 0_1^+)=0.027\,e^2b^2\), smaller than the experimental \(0.048\,e^2b^2\), likely because of missing contributions from the \(g_{9/2}\) orbital below the \(Z=50\) core. Monopole truncation converges rapidly to that full-space shell-model value, which indicates that the truncation error is smaller than the model-space error in this case [2507.11796].

For \(^{108}\)Xe, energies are essentially converged with less than \(20\%\) of the full \(J\)-scheme basis, and the \(B(E2;2_1^+\rightarrow 0_1^+)\) value also stabilizes rapidly. For \(^{110}\)Xe, the excitation energies \(E_x(2_1^+)\) and \(E_x(4_1^+)\) converge faster than the absolute binding energy, and around \(37\%\) of the basis is sufficient for an essentially converged low-lying spectrum. These examples show that the method remains effective for systems with both protons and neutrons and for nuclei with enhanced collectivity [2507.11796].

## 6. Relation to seniority truncation, limitations, and current scope

Monopole truncation is frequently paired with seniority truncation, but the two are conceptually distinct. Monopole truncation ranks partitions by mean-field energy, whereas seniority truncation emphasizes the dominance of monopole pairing by favoring basis states with low quasi-seniority. In the recent Pb and Sn tests, adding seniority information on top of monopole selection improves convergence for ground and low-lying yrast states, although high-spin states require higher-seniority components and therefore converge more slowly under an overly restrictive seniority bias [2507.11796].

The method also has explicit limitations. Centroid-based ranking does not encode all multipole correlations. Collective intruder states can be built from configurations whose monopole centroids are not especially low, and group-theoretical truncations based on seniority, SU(3), or Sp(3,R) can capture collectivity more directly. The centroid paper notes that computing analogous centroids for such symmetry-based partitions is not currently practical; existing symmetry-guided schemes therefore rely on different selection criteria [2511.03161].

Implementation details matter in mixed proton-neutron systems. In the current NuShellX realization, truncation is applied separately in the proton and neutron spaces, whereas KSHELL can use the full proton-neutron monopole Hamiltonian for truncation. The resulting differences in \(^{108}\)Xe and \(^{110}\)Xe show that monopole truncation is not purely abstract: the bookkeeping of partitions and the treatment of proton-neutron cross terms affect numerical performance [2507.11796].

At a broader methodological level, the PMMU and \(V_{MU}\) programs suggest that monopole-interaction-based truncation is most effective when the monopole sector has already been calibrated to reproduce shell evolution and effective single-particle-energy trends. Realistic interactions often get the multipole part approximately right but require monopole corrections; once those corrections are imposed, the valence space becomes naturally stratified by monopole energy, and truncation by partitions or optimized weights becomes physically meaningful [1510.05340][1312.6191].

In current usage, monopole-interaction-based truncation is therefore best understood as a hierarchy-building tool. It defines which configurations are mean-field favored, which orbitals are active near the Fermi surface, and which sectors of a very large CI space should be retained exactly before more selective pairing-, seniority-, or collectivity-based refinements are applied.

Source: https://www.emergentmind.com/topics/monopole-interaction-based-truncation