Realistic Shell Model (RSM) Overview

• RSM is a framework that derives valence-space Hamiltonians from realistic nucleon–nucleon interactions using many-body perturbation theory, thereby avoiding phenomenological inputs.
• It employs renormalization techniques such as folded-diagram summations and Lee–Suzuki iterations to construct effective operators that accurately reproduce spectroscopic and electroweak observables.
• RSM applications span nuclear spectroscopy, beta decay, and muon capture, with validations against no-core shell model benchmarks while addressing challenges like three-body forces and continuum effects.

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 (Covello et al., 2010, Coraggio et al., 2011, Coraggio et al., 2018, Otsuka, 2022).

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 (Coraggio et al., 2014). 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 (Coraggio et al., 2014).

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 (Coraggio et al., 2011). 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 (Gargano et al., 2014, Kaneko et al., 2013). 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 (Otsuka, 2022). 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$ (Covello et al., 2010, Gargano et al., 2014, Fukui et al., 2018). The goal is an effective Hamiltonian $^{18}$0 that reproduces the model-space sector of the exact spectrum: $^{18}$1

The irreducible vertex function is the $^{18}$2-box,

$^{18}$3

expanded perturbatively in valence-linked, irreducible Goldstone diagrams. In practical calculations, the $^{18}$4-box is typically evaluated through second or third order in $^{18}$5, including one- and two-body diagrams and, in some implementations, first-order normal-ordered one- and two-body contributions from chiral three-body forces (Gargano et al., 2014, Fukui et al., 2018). The folded-diagram series is then summed to all orders by Lee–Suzuki or Krenciglowa–Kuo iteration, yielding an energy-independent $^{18}$6 (Covello et al., 2010). For GT and $^{18}$7 applications in the mass range $^{18}$8–136, a Padé $^{18}$9 approximant is also used to accelerate convergence of the $V_{\rm low-k}$0-box expansion (Coraggio et al., 2018).

The one-body part of $V_{\rm low-k}$1 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 (Coraggio et al., 2011, Gargano et al., 2014).

Representative model spaces illustrate the breadth of the method.

Region	Core and valence space	Representative use
$V_{\rm low-k}$2 shell	$V_{\rm low-k}$3 over a $V_{\rm low-k}$4He core	Chiral $V_{\rm low-k}$5 benchmarks against NCSM (Fukui et al., 2018)
$V_{\rm low-k}$6 shell	$V_{\rm low-k}$7 shell above a $V_{\rm low-k}$8O core	Oxygen isotopes and ordinary muon capture (Gargano et al., 2014, Lyu et al., 8 Aug 2025)
$V_{\rm low-k}$9	full $A$0 over $A$1Ca; $A$2 over $A$3Ni; $A$4 over $A$5Sn	GT strengths and double-$A$6 decay (Coraggio et al., 2018)
$A$7 neutrons	$A$8 over $A$9Sn	Neutron-rich Sn spectroscopy (Covello et al., 2010)

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 = T + V_{NN},$0. For an operator $H = T + V_{NN},$1, one introduces a $H = T + V_{NN},$2-box analogous to the $H = T + V_{NN},$3-box,

$H = T + V_{NN},$4

together with two-energy vertices for folded diagrams. Following Suzuki–Okamoto, the effective operator is assembled from derivatives of $H = T + V_{NN},$5 and $H = T + V_{NN},$6. In the GT study, the $H = T + V_{NN},$7-box and its derivatives are retained through third order, while the effective operator is truncated at $H = T + V_{NN},$8; $H = T + V_{NN},$9 is found to be at the $$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)}.$$0 level, and $$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)}.$$1 and higher are negligibly small (Coraggio et al., 2018).

For Gamow–Teller transitions, the bare one-body operator is

$$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)}.$$2

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 $$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)}.$$3–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 (Coraggio et al., 2018).

For neutrinoless double-$$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)}.$$4 decay, the closure approximation is used to define the bare two-body operator,

$$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)}.$$5

and the total matrix element is written

$$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)}.$$6

In this setting, short-range correlations are incorporated by applying the same $$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)}.$$7 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 (Coraggio et al., 2020). The $$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)}.$$8Ge-specific study emphasizes that these blocking terms counteract Pauli-violating contributions and reduce the total $$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)}.$$9 matrix element by about $U$0, mostly through the GT channel (Itaco et al., 2019).

Electromagnetic and other electroweak operators are treated in the same spirit. Fully microscopic calculations include core-polarization corrections to $U$1 and $U$2 operators up to third order (Coraggio et al., 2011), and the 2025 ordinary muon capture study extends the same machinery to a process with exchange momenta of approximately $U$3 MeV, with one-body plus leading two-body corrections retained in many-valence-nucleon systems (Lyu et al., 8 Aug 2025). 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 (Coraggio et al., 2018).

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 $U$4 counts the number of nucleons not participating in a collective $U$5-pair condensate. For semimagic Ca isotopes in the $U$6 shell, basis states with $U$7 are built from

$U$8

orthonormalized by Gram–Schmidt, and benchmarked against full-shell-model calculations (Caprio et al., 2012). The resulting $U$9 truncation captures binding-energy trends, occupations, and several electromagnetic observables with low computational cost, while static quadrupole moments expose missing non-pairing correlations (Caprio et al., 2012).

A different realization is the PMMU Hamiltonian,

$$H = H_0 + H_1, \qquad H_0 = T + U, \qquad H_1 = V_{NN} - U,$$0

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 $$H = H_0 + H_1, \qquad H_0 = T + U, \qquad H_1 = V_{NN} - U,$$1 tensor components of the MU force define the monopole systematics, after which a small number of monopoles are fine-tuned: 10 terms in the $$H = H_0 + H_1, \qquad H_0 = T + U, \qquad H_1 = V_{NN} - U,$$2 shell and 14 terms in the $$H = H_0 + H_1, \qquad H_0 = T + U, \qquad H_1 = V_{NN} - U,$$3 shell. The resulting fits achieve rms binding-energy errors of $$H = H_0 + H_1, \qquad H_0 = T + U, \qquad H_1 = V_{NN} - U,$$4 MeV and $$H = H_0 + H_1, \qquad H_0 = T + U, \qquad H_1 = V_{NN} - U,$$5 MeV, respectively (Kaneko et al., 2013). 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,

$$H = H_0 + H_1, \qquad H_0 = T + U, \qquad H_1 = V_{NN} - U,$$6

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 (Otsuka, 2022). 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 $$H = H_0 + H_1, \qquad H_0 = T + U, \qquad H_1 = V_{NN} - U,$$7 shell, starting from a chiral $$H = H_0 + H_1, \qquad H_0 = T + U, \qquad H_1 = V_{NN} - U,$$8 $$H = H_0 + H_1, \qquad H_0 = T + U, \qquad H_1 = V_{NN} - U,$$9 interaction alone, low-lying spectra of $Q=1-P$0Li, $Q=1-P$1Li, $Q=1-P$2B, $Q=1-P$3Be, $Q=1-P$4B, $Q=1-P$5B, $Q=1-P$6C, and $Q=1-P$7C agree with NCSM to within a few hundred keV. Without three-body forces, however, the $Q=1-P$8B $Q=1-P$9 inversion remains incorrect; inclusion of the genuine chiral $P$0 three-body force at first order largely cures the spectroscopic defects and stabilizes the $P$1 ESPE gap at $P$2–$P$3 MeV (Fukui et al., 2018).

In the $P$4 and $P$5 shells, the method reproduces classical benchmark spectra while also exposing sensitivity to missing three-body contributions. Fully microscopic calculations for $P$6O give $P$7 MeV versus $P$8 MeV experimentally and $P$9 MeV versus $^{18}$00 MeV (Coraggio et al., 2011). In neutron-rich oxygen, a third-order $^{18}$01 calculation reproduces the rise of the $^{18}$02 energies at $^{18}$03 and $^{18}$04, but the unshifted one-body spectrum overbinds as neutron number increases and even predicts $^{18}$05O and $^{18}$06O to be bound; shifting the single-particle spectrum upward by $^{18}$07 keV restores the binding-energy slope and the drip line at $^{18}$08 (Gargano et al., 2014). For $^{18}$09 isotones, inclusion of the $^{18}$10 orbital is necessary for quantitative agreement with the observed drop in $^{18}$11 and the large $^{18}$12 values in Cr and Fe (Gargano et al., 2014).

In the $^{18}$13Sn region, realistic interactions derived from CD-Bonn via $^{18}$14 and folded-diagram methods reproduce the anomalously low $^{18}$15 energy of $^{18}$16Sn and its transition strengths. For $^{18}$17Sn, the calculated yrast sequence $^{18}$18 lies at $^{18}$19, $^{18}$20, and $^{18}$21 MeV, close to the experimental $^{18}$22, $^{18}$23, and $^{18}$24 MeV, while $^{18}$25 W.u. compares with $^{18}$26 W.u. (Covello et al., 2010). The same study predicts almost constant $^{18}$27 and $^{18}$28 energies from $^{18}$29Sn to $^{18}$30Sn and finds no signature of a new shell closure at $^{18}$31 (Covello et al., 2010).

The pairing analysis around $^{18}$32Sn identifies the microscopic origin of proton–neutron asymmetry in this region. For two valence protons in $^{18}$33Te, the diagonal $^{18}$34 matrix element evolves from $^{18}$35 MeV to $^{18}$36 MeV because the one-particle–one-hole core-polarization contribution is $^{18}$37 MeV. For two valence neutrons in $^{18}$38Sn, the corresponding values are $^{18}$39 MeV, $^{18}$40 MeV, and $^{18}$41 MeV. The resulting effective pairing matrix element is therefore roughly $^{18}$42 stronger in the proton case, consistent with $^{18}$43 MeV versus $^{18}$44 MeV (Covello et al., 2012).

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 $^{18}$45. In the perturbative RSM treatment of $^{18}$46–136 nuclei, the bare GT operator overestimates $^{18}$47 and $^{18}$48 matrix elements by factors of $^{18}$49–$^{18}$50, whereas the derived $^{18}$51 yields very good agreement with measured $^{18}$52, with differences $^{18}$53, and reproduces $^{18}$54 matrix elements within experimental errors for $^{18}$55. No ad hoc quenching $^{18}$56 or renormalization of $^{18}$57 and $^{18}$58 is used; all renormalization comes from MBPT diagrams (Coraggio et al., 2018).

The same program has been extended to $^{18}$59 decay. For five standard candidates, the fully renormalized operator reduces the bare light-neutrino-exchange matrix element as follows.

Decay	$^{18}$60 bare	$^{18}$61
$^{18}$62	0.53	0.30
$^{18}$63	3.35	2.66
$^{18}$64	3.30	2.72
$^{18}$65	3.27	3.16
$^{18}$66	2.47	2.39

These calculations indicate a quenching of $^{18}$67 by $^{18}$68–$^{18}$69 in most cases, milder than the quenching required for single-$^{18}$70 and $^{18}$71 decay, and identify the GT channel as the dominant source of renormalization (Coraggio et al., 2020). In the dedicated $^{18}$72Ge study, three different shell-model Hamiltonians with the same bare operator give $^{18}$73, $^{18}$74, and $^{18}$75, differing by $^{18}$76, while the second-order renormalized operator plus Pauli blocking reduces the total matrix element from $^{18}$77 to $^{18}$78 (Itaco et al., 2019).

Ordinary muon capture provides a complementary high-momentum test. In the $^{18}$79-shell RSM study, spectroscopy and electroweak observables are described more successfully by the chiral $^{18}$80 Hamiltonian than by the pure $^{18}$81 Hamiltonian, and with the fully renormalized ordinary muon-capture operator the chiral calculation reproduces absolute partial rates to better than $^{18}$82 in most channels, whereas the $^{18}$83 version still underestimates by up to $^{18}$84–$^{18}$85 (Lyu et al., 8 Aug 2025). This suggests that simultaneous control of spectroscopy and weak processes at $^{18}$86 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 $^{18}$87-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 (Covello et al., 2010). The cumulative picture is therefore twofold: the RSM has demonstrated quantitative power for spectroscopy, GT strengths, $^{18}$88, $^{18}$89, 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.