---
title: Realistic Shell Model (RSM) Overview
url: https://www.emergentmind.com/topics/realistic-shell-model-rsm
type: topic
---

# Realistic Shell Model (RSM) Overview

The realistic shell model (RSM) is a shell-model framework in which the valence-space Hamiltonian, and in the fully microscopic formulation also transition operators, are derived from realistic nucleon–nucleon interactions for a truncated model space through renormalization and many-body perturbation theory rather than fixed purely phenomenologically. In the literature considered here, the term covers both the canonical folded-diagram approach based on high-precision or chiral forces and related shell-model constructions that retain realistic monopole and tensor systematics in larger-scale spectroscopic studies. Across these formulations, the central objective is the same: to connect shell-model spectroscopy and electroweak observables to free-space nuclear forces while controlling the effects of degrees of freedom excluded from the chosen valence space [1001.3520] [1101.4460] [1812.04292] [2201.05443].

## 1. Historical development and conceptual scope

The modern RSM descends from the Kuo–Brown program. The 1966 calculation for \(^{18}\)O and \(^{18}\)F, based on an effective interaction derived from the Hamada–Johnston nucleon–nucleon potential, is described as the first successful attempt to provide a description of nuclear structure properties starting from the free nucleon–nucleon potential [1403.5470]. Subsequent work replaced hard-core interactions and limited second-order treatments by chiral forces, \(V_{\rm low-k}\) or related soft interactions, higher-order Goldstone-diagram expansions, and all-order summations of folded diagrams [1403.5470].

In its fully microscopic form, the RSM aims to avoid phenomenological input altogether. One review states explicitly that no phenomenological input is needed because single-particle energies, matrix elements of the two-body interaction, and matrix elements of the electromagnetic multipole operators are derived theoretically within time-dependent degenerate linked-diagram perturbation theory [1101.4460]. Other works adopt a less stringent implementation, for example by replacing the one-body part with experimental single-particle energies or by combining physically motivated pairing, multipole, and monopole terms with limited monopole adjustments [1403.2398] [1312.6191]. This suggests that “realistic shell model” denotes a research program rather than a single fixed Hamiltonian.

A second broadening of scope appears in shell-model studies that emphasize monopole structure and shell evolution. In that literature, the RSM is tied to realistic monopole and tensor components of the interaction, extended Krenciglowa–Kuo decoupling in multishell spaces, and diagnostics such as T-plots for intrinsic-shape analysis [2201.05443]. The common thread remains the attempt to preserve the microscopic content of the underlying nuclear force while working in a practical valence space.

## 2. Microscopic construction of the effective Hamiltonian

The standard starting point is the \(A\)-body Hamiltonian
\[
H = T + V_{NN},
\]
or, when explicit three-body terms are retained,
\[
H_{\rm int} = \sum_i t_i - T_{\rm cm} + \sum_{i<j} V_{NN}^{(i,j)} + \sum_{i<j<k} V_{3N}^{(i,j,k)}.
\]
An auxiliary one-body potential \(U\), usually harmonic-oscillator, defines the partition
\[
H = H_0 + H_1, \qquad H_0 = T + U, \qquad H_1 = V_{NN} - U,
\]
with \(Q=1-P\) the complement of the valence-space projector \(P\) [1001.3520] [1403.2398] [1806.01573]. The goal is an effective Hamiltonian \(H_{\rm eff}\) that reproduces the model-space sector of the exact spectrum:
\[
H_{\rm eff} P \Psi_\alpha = E_\alpha P \Psi_\alpha.
\]

The irreducible vertex function is the \(\hat Q\)-box,
\[
\hat Q(\omega) = P H_1 P + P H_1 Q \frac{1}{\omega - Q H Q} Q H_1 P,
\]
expanded perturbatively in valence-linked, irreducible Goldstone diagrams. In practical calculations, the \(\hat Q\)-box is typically evaluated through second or third order in \(H_1\), including one- and two-body diagrams and, in some implementations, first-order normal-ordered one- and two-body contributions from chiral three-body forces [1403.2398] [1806.01573]. The folded-diagram series is then summed to all orders by Lee–Suzuki or Krenciglowa–Kuo iteration, yielding an energy-independent \(H_{\rm eff}\) [1001.3520]. For GT and \(\beta\beta\) applications in the mass range \(A=48\)–136, a Padé \([2|1]\) approximant is also used to accelerate convergence of the \(\hat Q\)-box expansion [1812.04292].

The one-body part of \(H_{\rm eff}\) defines theoretical single-particle energies, while the two-body part provides the valence-space two-body matrix elements. Some calculations retain the fully theoretical one-body sector; others replace it by experimental single-particle or single-hole energies, explicitly to absorb missing three-nucleon physics or to exploit spectroscopic information in neighboring nuclei [1101.4460] [1403.2398].

Representative model spaces illustrate the breadth of the method.

| Region | Core and valence space | Representative use |
|---|---|---|
| \(p\) shell | \(0p_{3/2},0p_{1/2}\) over a \(^{4}\)He core | Chiral \(3N\) benchmarks against NCSM [1806.01573] |
| \(sd\) shell | \(sd\) shell above a \(^{16}\)O core | Oxygen isotopes and ordinary muon capture [1403.2398] [2508.06221] |
| \(A\approx 48,76\text{–}82,130\text{–}136\) | full \(fp\) over \(^{40}\)Ca; \(f_{5/2}p_{3/2}p_{1/2}g_{9/2}\) over \(^{56}\)Ni; \(g_{7/2}d_{5/2}d_{3/2}s_{1/2}h_{11/2}\) over \(^{100}\)Sn | GT strengths and double-\(\beta\) decay [1812.04292] |
| \(N=82\text{–}126\) neutrons | \(0h_{9/2},1f_{7/2},1f_{5/2},2p_{3/2},2p_{1/2},0i_{13/2}\) over \(^{132}\)Sn | Neutron-rich Sn spectroscopy [1010.4194] |

## 3. Renormalization of transition operators

A defining feature of the fully microscopic RSM is the consistent renormalization of transition operators in the same perturbative framework used for \(H_{\rm eff}\). For an operator \(\Theta\), one introduces a \(\widehat\Theta\)-box analogous to the \(\hat Q\)-box,
\[
\widehat\Theta(\epsilon) = P\Theta P + P\Theta Q \frac{1}{\epsilon - QHQ} Q H_1 P + \cdots,
\]
together with two-energy vertices for folded diagrams. Following Suzuki–Okamoto, the effective operator is assembled from derivatives of \(\widehat\Theta\) and \(\hat Q\). In the GT study, the \(\hat Q\)-box and its derivatives are retained through third order, while the effective operator is truncated at \(\chi_2\); \(\chi_2\) is found to be at the \(\sim 1\%\) level, and \(\chi_3\) and higher are negligibly small [1812.04292].

For Gamow–Teller transitions, the bare one-body operator is
\[
O_{GT}^0 = \sum_{i=1}^{A} \sigma_i \tau_i^\pm .
\]
In the RSM it is renormalized by core polarization, folded diagrams, and blocking corrections induced by the occupation of valence orbitals. In the calculations for \(A=48\)–136, transition-operator diagrams include one-body terms through second order, two-body core-polarization terms through second order, and folded diagrams; blocking corrections are treated via density-dependent one-body operators [1812.04292].

For neutrinoless double-\(\beta\) decay, the closure approximation is used to define the bare two-body operator,
\[
\Theta_\alpha = \sum_{i<j} \tau_i^- \tau_j^- O_\alpha(r_{ij}), \qquad \alpha = F,GT,T,
\]
and the total matrix element is written
\[
M^{0\nu} = M^{0\nu}_{GT} - \Bigl(\frac{g_V}{g_A}\Bigr)^2 M^{0\nu}_F .
\]
In this setting, short-range correlations are incorporated by applying the same \(V_{\rm low-k}\) transformation used for the interaction to the decay operator, and Pauli-principle corrections are approximated by summing connected second-order three-body diagrams over occupied orbitals to produce a density-dependent two-body correction [2001.00890]. The \(^{76}\)Ge-specific study emphasizes that these blocking terms counteract Pauli-violating contributions and reduce the total \(0\nu\beta\beta\) matrix element by about \(10\%\), mostly through the GT channel [1911.12117].

Electromagnetic and other electroweak operators are treated in the same spirit. Fully microscopic calculations include core-polarization corrections to \(E2\) and \(M1\) operators up to third order [1101.4460], and the 2025 ordinary muon capture study extends the same machinery to a process with exchange momenta of approximately \(100\) MeV, with one-body plus leading two-body corrections retained in many-valence-nucleon systems [2508.06221]. A notable limitation is explicit in the GT renormalization study: meson-exchange two-body currents are not included there and are left for future work [1812.04292].

## 4. Truncation schemes, monopole constructions, and large-scale realizations

One strand of the RSM literature addresses computational dimensionality by truncating a realistic shell-model Hamiltonian rather than altering its microscopic origin. In generalized seniority, the truncation parameter \(v\) counts the number of nucleons not participating in a collective \(S\)-pair condensate. For semimagic Ca isotopes in the \(pf\) shell, basis states with \(v\le 3\) are built from
\[
S^\dagger = \sum_a \frac{\alpha_a}{2}\,\hat\jmath_a\,A^\dagger(aa,J=0),
\]
orthonormalized by Gram–Schmidt, and benchmarked against full-shell-model calculations [1203.4856]. The resulting \(v\le 3\) truncation captures binding-energy trends, occupations, and several electromagnetic observables with low computational cost, while static quadrupole moments expose missing non-pairing correlations [1203.4856].

A different realization is the PMMU Hamiltonian,
\[
H = H_0 + H_{PM} + H_m^{MU},
\]
which combines pairing plus quadrupole and octupole multipoles with a monopole Hamiltonian derived from the monopole-based universal force. In this construction, the central Gaussian and \(\pi+\rho\) tensor components of the MU force define the monopole systematics, after which a small number of monopoles are fine-tuned: 10 terms in the \(pf\) shell and 14 terms in the \(pf_{5/2}g_{9/2}\) shell. The resulting fits achieve rms binding-energy errors of \(\sim 0.7\) MeV and \(\sim 0.69\) MeV, respectively [1312.6191]. The PMMU work explicitly presents this as a practical realization of the realistic shell-model ideal for heavier regions.

Large-scale shell-model studies centered on monopole physics provide a third line of development. In that formulation, the interaction is decomposed into monopole and multipole parts, and shell evolution is driven by the monopole contribution to effective single-particle energies,
\[
\epsilon_a(N) = \epsilon_{0;a} + \sum_b V^{\rm mono}(ab)\,n_b .
\]
This perspective underlies the distinction between type-I shell evolution, driven by orbital filling, and type-II shell evolution, driven by particle–hole excitations within a nucleus. It also motivates T-plot visualization in Monte Carlo shell-model calculations and the notion of self-organization of shell structure in collective bands [2201.05443]. These developments remain within the broader RSM literature because the underlying monopole behavior is tied to realistic central and tensor forces.

## 5. Spectroscopic performance across the nuclear chart

For light nuclei, the RSM has been benchmarked directly against no-core methods. In the \(p\) shell, starting from a chiral \(N^3LO\) \(NN\) interaction alone, low-lying spectra of \(^{6}\)Li, \(^{8}\)Li, \(^{8}\)B, \(^{8}\)Be, \(^{10}\)B, \(^{11}\)B, \(^{12}\)C, and \(^{13}\)C agree with NCSM to within a few hundred keV. Without three-body forces, however, the \(^{10}\)B \(3^+–1^+\) inversion remains incorrect; inclusion of the genuine chiral \(N^2LO\) three-body force at first order largely cures the spectroscopic defects and stabilizes the \(p_{1/2}–p_{3/2}\) ESPE gap at \(\sim 4\)–\(5\) MeV [1806.01573].

In the \(sd\) and \(pf\) shells, the method reproduces classical benchmark spectra while also exposing sensitivity to missing three-body contributions. Fully microscopic calculations for \(^{18}\)O give \(E(2_1^+) = 1.905\) MeV versus \(1.982\) MeV experimentally and \(E(4_1^+) = 3.472\) MeV versus \(3.555\) MeV [1101.4460]. In neutron-rich oxygen, a third-order \(N^3LOW\) calculation reproduces the rise of the \(2^+\) energies at \(N=14\) and \(N=16\), but the unshifted one-body spectrum overbinds as neutron number increases and even predicts \(^{26}\)O and \(^{28}\)O to be bound; shifting the single-particle spectrum upward by \(427\) keV restores the binding-energy slope and the drip line at \(N=16\) [1403.2398]. For \(N=40\) isotones, inclusion of the \(\nu d_{5/2}\) orbital is necessary for quantitative agreement with the observed drop in \(E(2^+)\) and the large \(B(E2)\) values in Cr and Fe [1403.2398].

In the \(^{132}\)Sn region, realistic interactions derived from CD-Bonn via \(V_{\rm low-k}\) and folded-diagram methods reproduce the anomalously low \(2^+\) energy of \(^{134}\)Sn and its transition strengths. For \(^{134}\)Sn, the calculated yrast sequence \(2^+_1,4^+_1,6^+_1\) lies at \(0.730\), \(1.270\), and \(1.850\) MeV, close to the experimental \(0.726\), \(1.303\), and \(1.898\) MeV, while \(B(E2;2^+\!\to0^+) = 1.64\) W.u. compares with \(1.42\pm0.20\) W.u. [1010.4194]. The same study predicts almost constant \(2^+\) and \(4^+\) energies from \(^{134}\)Sn to \(^{140}\)Sn and finds no signature of a new shell closure at \(N=90\) [1010.4194].

The pairing analysis around \(^{132}\)Sn identifies the microscopic origin of proton–neutron asymmetry in this region. For two valence protons in \(^{134}\)Te, the diagonal \(J=0\) matrix element evolves from \(V_{\rm low-k}(0^+)=-0.248\) MeV to \(V_{\rm eff}(0^+)=-0.967\) MeV because the one-particle–one-hole core-polarization contribution is \(V_{1p1h}(0^+)=-0.554\) MeV. For two valence neutrons in \(^{134}\)Sn, the corresponding values are \(V_{\rm low-k}(0^+)=-0.403\) MeV, \(V_{1p1h}(0^+)=-0.100\) MeV, and \(V_{\rm eff}(0^+)=-0.654\) MeV. The resulting effective pairing matrix element is therefore roughly \(50\%\) stronger in the proton case, consistent with \(E(2^+;^{134}{\rm Te})=1.28\) MeV versus \(E(2^+;^{134}{\rm Sn})=0.726\) MeV [1211.4758].

## 6. Electroweak observables, misconceptions, and open problems

A recurrent assumption in shell-model phenomenology is that Gamow–Teller observables must be reproduced by imposing an empirical quenching of \(g_A\). In the perturbative RSM treatment of \(A=48\)–136 nuclei, the bare GT operator overestimates \(\Sigma B(p,n)\) and \(2\nu\beta\beta\) matrix elements by factors of \(\sim 2\)–\(4\), whereas the derived \(O_{GT}^{\rm eff}\) yields very good agreement with measured \(\Sigma B(p,n)\), with differences \(\lesssim 10\%\), and reproduces \(2\nu\beta\beta\) matrix elements within experimental errors for \(A=48,76,82,130,136\). No ad hoc quenching \(g_A \to q\cdot g_A\) or renormalization of \(g_l\) and \(g_s\) is used; all renormalization comes from MBPT diagrams [1812.04292].

The same program has been extended to \(0\nu\beta\beta\) decay. For five standard candidates, the fully renormalized operator reduces the bare light-neutrino-exchange matrix element as follows.

| Decay | \(M^{0\nu}\) bare | \(M^{0\nu}(\Theta_{\rm eff})\) |
|---|---:|---:|
| \(^{48}\mathrm{Ca}\to^{48}\mathrm{Ti}\) | 0.53 | 0.30 |
| \(^{76}\mathrm{Ge}\to^{76}\mathrm{Se}\) | 3.35 | 2.66 |
| \(^{82}\mathrm{Se}\to^{82}\mathrm{Kr}\) | 3.30 | 2.72 |
| \(^{130}\mathrm{Te}\to^{130}\mathrm{Xe}\) | 3.27 | 3.16 |
| \(^{136}\mathrm{Xe}\to^{136}\mathrm{Ba}\) | 2.47 | 2.39 |

These calculations indicate a quenching of \(M^{0\nu}\) by \(\sim 10\)–\(20\%\) in most cases, milder than the quenching required for single-\(\beta\) and \(2\nu\beta\beta\) decay, and identify the GT channel as the dominant source of renormalization [2001.00890]. In the dedicated \(^{76}\)Ge study, three different shell-model Hamiltonians with the same bare operator give \(M^{0\nu}=3.40\), \(2.96\), and \(3.25\), differing by \(\lesssim 15\%\), while the second-order renormalized operator plus Pauli blocking reduces the total matrix element from \(2.79\) to \(2.44\) [1911.12117].

Ordinary muon capture provides a complementary high-momentum test. In the \(sd\)-shell RSM study, spectroscopy and electroweak observables are described more successfully by the chiral \(2N+3N\) Hamiltonian than by the pure \(V_{\rm low-k}\) Hamiltonian, and with the fully renormalized ordinary muon-capture operator the chiral calculation reproduces absolute partial rates to better than \(\sim 15\%\) in most channels, whereas the \(V_{\rm low-k}\) version still underestimates by up to \(30\)–\(40\%\) [2508.06221]. This suggests that simultaneous control of spectroscopy and weak processes at \(q\approx 100\) MeV is a stringent discriminator among realistic interactions.

Open problems are identified consistently across the literature. They include explicit three-body forces in medium-mass valence spaces, higher-order \(\hat Q\)-box terms, continuum coupling near driplines, intruder-state and nonperturbative effects, uncertainty estimates from regulator and cutoff variation, full three-body effective operators, and consistent two-body electroweak currents [1001.3520]. The cumulative picture is therefore twofold: the RSM has demonstrated quantitative power for spectroscopy, GT strengths, \(2\nu\beta\beta\), \(0\nu\beta\beta\), and ordinary muon capture, but its long-term completion depends on incorporating the many-body and current operators that several of its own benchmark studies identify as the next required step.

Source: https://www.emergentmind.com/topics/realistic-shell-model-rsm