pith. sign in

arxiv: 2606.04055 · v1 · pith:ZNWKV6PKnew · submitted 2026-06-02 · 🪐 quant-ph · nucl-th

Boson Models with Interactions of Arbitrary Order

P. Van Isacker This is my paper

Pith reviewed 2026-06-28 09:39 UTC · model grok-4.3

classification 🪐 quant-ph nucl-th
keywords boson modelsmatrix elementsk-body interactionsangular momentummany-body systemstensor operatorsrecursion
0
0 comments X

The pith

Closed formulas are given for matrix elements between N-boson states for interactions of any order k when there is one or two boson types.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops explicit ways to evaluate matrix elements for quantal systems of N bosons of p types, possibly including scalar s bosons, under a rotationally invariant Hamiltonian that conserves boson number. For p equal to one or two it supplies closed expressions valid for any interaction order k; for arbitrary p it supplies a recursive construction. These expressions let any Hamiltonian matrix element be written symbolically as a linear combination of the underlying k-body interaction matrix elements. The same machinery also covers matrix elements of general tensor operators that need not be scalars or conserve particle number. The resulting formulas support both symbolic manipulation and numerical diagonalization of the many-boson problem.

Core claim

Closed formulas are given for matrix elements between N-boson states for any k if p=1 and p=2. A recursive procedure is defined for arbitrary k and p. With the expressions derived it is possible to express symbolically a Hamiltonian matrix element between N-boson states as a linear combination of k-body interaction matrix elements. More generally, the formulas allow the evaluation of matrix elements of tensor operators that are not necessarily scalar nor boson-number conserving.

What carries the argument

Closed expressions and recursive procedure for matrix elements of k-body interactions (and general tensor operators) between states built from bosons carrying definite angular momenta l1 to lp.

If this is right

  • Any Hamiltonian matrix can be assembled directly from its k-body building blocks without expanding every term by hand.
  • The same algebraic machinery yields matrix elements for non-scalar and particle-number-nonconserving operators.
  • Numerical codes can now treat interactions of arbitrarily high order once the k-body reduced matrix elements are supplied.
  • The recursive construction for general p reduces the problem of many boson types to repeated application of the p=2 case.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The approach may simplify large-scale diagonalizations in algebraic boson models used for nuclear spectra.
  • It supplies a route to embed higher-order interactions into existing computer programs that previously handled only two-body terms.
  • Extension to mixed boson-fermion systems could follow by treating the fermions as an additional species with their own angular momenta.

Load-bearing premise

The Hamiltonian is rotationally invariant, conserves boson number, and acts on states built from bosons with definite angular momenta.

What would settle it

Explicit enumeration of all matrix elements for a three-body interaction (k=3) in a three-type boson system (p=3) with small N, compared against the output of the recursive procedure, would confirm or refute the formulas.

Figures

Figures reproduced from arXiv: 2606.04055 by P. Van Isacker.

Figure 1
Figure 1. Figure 1: FIG. 1: The expectation value of the fourth-order operators (68) in the ground state [PITH_FULL_IMAGE:figures/full_fig_p026_1.png] view at source ↗
read the original abstract

The paper considers quantal many-boson systems that are described by a rotationally invariant and boson-number conserving Hamiltonian. The properties of a generic model are studied which treats N bosons of p different kinds with non-zero angular momenta l_1,l_2,...,l_p, possibly augmented with a (number of) scalar s boson(s). The order k of the interaction between the bosons is arbitrary and closed formulas are given for matrix elements between N-boson states for any k if p=1 and p=2. A recursive procedure is defined for arbitrary k and p. With the expressions derived in the paper it is possible to express symbolically a Hamiltonian matrix element between N-boson states as a linear combination of k-body interaction matrix elements. More generally, the formulas allow the evaluation of matrix elements of tensor operators that are not necessarily scalar nor boson-number conserving. The numerical implementation of the formalism is discussed and illustrated with a few examples.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 2 minor

Summary. The paper considers quantal many-boson systems described by a rotationally invariant and boson-number conserving Hamiltonian acting on states built from N bosons of p different kinds with angular momenta l_1 to l_p, possibly augmented by scalar s bosons. It provides closed formulas for matrix elements of k-body tensor operators between N-boson states when p=1 and p=2 (for arbitrary k), together with a recursive procedure for general p and k. These allow any N-boson Hamiltonian matrix element to be expressed symbolically as a linear combination of elementary k-body interaction matrix elements; the formalism also covers non-scalar and non-number-conserving tensor operators. Numerical implementation is discussed and illustrated with examples.

Significance. If the derivations hold, the work supplies a practical algebraic toolkit for boson models with arbitrary-order interactions, reducing the evaluation of large N-body matrix elements to k-body building blocks via standard angular-momentum recoupling and boson commutation relations. The closed forms for small p and the recursion for general p, together with the explicit numerical examples, constitute a clear advance over purely numerical approaches and support both symbolic manipulation and reproducible computations.

minor comments (2)
  1. Abstract: the phrase 'a few examples' for numerical implementation would be more informative if the specific values of p, k, N, and l_i were stated, allowing readers to gauge the scope immediately.
  2. The manuscript would benefit from an explicit statement (near the model definition) of whether the s bosons are optional or always present, and how their number is treated in the recursion.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive summary, significance assessment, and recommendation of minor revision. No major comments were provided in the report, so there are no specific points requiring point-by-point response or rebuttal.

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained in angular-momentum algebra

full rationale

The paper derives closed formulas and a recursive procedure for N-boson matrix elements of k-body tensor operators directly from rotational invariance, boson-number conservation, angular-momentum recoupling coefficients, and standard commutation relations applied to states built from p species of l_i bosons (plus optional s bosons). No step reduces a claimed prediction or result to a fitted parameter, a self-citation chain, an imported uniqueness theorem, or an ansatz smuggled from prior work by the same authors; the central expressions are presented as algebraic identities obtained from first-principles recoupling rather than from data or prior self-referential results. The derivation is therefore independent of the paper's own inputs and qualifies as non-circular.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work is a derivation within standard quantum mechanics; no free parameters, ad-hoc axioms, or new entities are introduced in the abstract.

axioms (2)
  • standard math Boson states are constructed in the usual way from creation operators with definite angular momentum.
    Invoked implicitly when defining N-boson states with angular momenta l1,...,lp.
  • domain assumption The Hamiltonian is rotationally invariant and conserves boson number.
    Stated at the opening of the abstract as the class of models considered.

pith-pipeline@v0.9.1-grok · 5682 in / 1352 out tokens · 21207 ms · 2026-06-28T09:39:17.038843+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

42 extracted references · 1 canonical work pages

  1. [1]

    The first one equals ⟨sN−n |(s†)k−m(˜s)k−m′ |sN−n ′ ⟩p (k−m)!(k−m ′)! ,(33) withn≡n 1 +n 2 andn ′ ≡n ′ 1 +n ′

    With the result B† Λ1 ×B † Λ2 (L) M † = (−)L−M ˜BΛ1 × ˜BΛ2 (L) −M ,(29) the interaction can be rewritten as X L X Λ1Λ2 X Λ′ 1Λ′ 2 vkL Λ1Λ2,Λ′ 1Λ′ 2 (s†)k−m(˜s)k−m′ p (k−m)!(k−m ′)! (−)LB† Λ1Λ2L · ˜BΛ′ 1Λ′ 2L ,(30) where theBoperators are defined as B† Λ1Λ2LM ≡ B† Λ1 ×B † Λ2 (L) M , ˜BΛ1Λ2LM ≡ ˜BΛ1 × ˜BΛ2 (L) M .(31) Matrix elements of a single interaction...

  2. [2]

    This factor is known from Eq. (22) while the second is given by ⟨NΓ 1Γ2;J|(−) LB† Λ1Λ2L · ˜BΛ′ 1Λ′ 2L|NΓ ′ 1Γ′ 2;J⟩ = (−)J ′ 1+J2+L′ 1+L2+L+J(2L+ 1) X R L1 L2 L L′ 2 L′ 1 R J1 J2 J J ′ 2 J ′ 1 R × ⟨Γ1∥ B† Λ1 × ˜BΛ′ 1 (R) ∥Γ′ 1⟩⟨Γ2∥ B† Λ2 × ˜BΛ′ 2 (R) ∥Γ′ 2⟩.(34) The reduced matrix elements on the right-hand side of Eq. (34) are known from Eq. (25). Matrix...

  3. [3]

    With the current algorithm calculations insdg-IBM can be carried out by restricting the number ofgbosons to a lower value than the total boson numberN

    and therefore an alternative, hitherto unknown, algorithm for the construction of CFPs must be devised. With the current algorithm calculations insdg-IBM can be carried out by restricting the number ofgbosons to a lower value than the total boson numberN. IX. EXAMPLES A. The classical limit of thesd-IBM Geometric insight into algebraic models can be obtai...

  4. [4]

    This proves that a triaxial shape associated with (cos 3γ)2 requires at least cubic interactions between thedbosons [32]

    If mod 3(υ+υ ′) = 0, thent max = (υ+υ ′)/3. This proves that a triaxial shape associated with (cos 3γ)2 requires at least cubic interactions between thedbosons [32]. Since the 3d-boson states inv 2 ddd·ddd have seniorityυ= 1, this interaction does not contribute to (cos 3γ)2, in contrast to the interactionsv L ddd·ddd withL̸= 2, which haveυ+υ ′ = 6

  5. [5]

    As a consequence, quadratic interactions between thedbosons are independent ofγin the classical limit

    If mod 3(υ+υ ′) = 1, thent max = (υ+υ ′−4)/3. As a consequence, quadratic interactions between thedbosons are independent ofγin the classical limit. Similarly, since the sum of the seniorities inv 2 dddd2·dddd2 orv 4 dddd2·dddd2 isυ+υ ′ = 4, these four-body interactions are also independent ofγin the classical limit

  6. [6]

    Quartic interactions between thed bosons can have at mostυ+υ ′ = 8, in which caset max = 2, yielding a term in (cos 3γ)2 (t= 2) and one independent ofγ(t= 0)

    If mod 3(υ+υ ′) = 2, thent max = (υ+υ ′ −2)/3. Quartic interactions between thed bosons can have at mostυ+υ ′ = 8, in which caset max = 2, yielding a term in (cos 3γ)2 (t= 2) and one independent ofγ(t= 0). 23 B. The ˆQµ + ˆLµ expansion of thesd-Hamiltonian Given the importance of quadrupole deformation in nuclei, several authors have explored the results ...

  7. [7]

    Iachello, R.D

    F. Iachello, R.D. Levine,Algebraic Theory of Molecules(Oxford University Press, Oxford, 1995)

    1995

  8. [8]

    Iachello, A

    F. Iachello, A. Arima,The Interacting Boson Model(Cambridge University Press, Cambridge, 1987)

    1987

  9. [9]

    Arima, T

    A. Arima, T. Ohtsuka, F. Iachello, I. Talmi,Collective nuclear states as symmetric couplings of proton and neutron excitations, Phys. Lett. B66(1977) 205

    1977

  10. [10]

    Elliott, A.P

    J.P. Elliott, A.P. White,An isospin invariant form of the interacting boson model, Phys. Lett. B97(1980) 169. 29

    1980

  11. [11]

    de-Shalit, I

    A. de-Shalit, I. Talmi,Nuclear Shell Theory(Academic Press, New York, 1963)

    1963

  12. [12]

    Talmi,Simple Models of Complex Nuclei

    I. Talmi,Simple Models of Complex Nuclei. The Shell Model and Interacting Boson Model (Harwood, Chur, 1993)

    1993

  13. [13]

    Skouras, S

    L.D. Skouras, S. Kossionides,Generalized fractional parentage coefficients for shell-model cal- culations, Comp. Phys. Comm.39(1986) 197

    1986

  14. [14]

    Racah,Theory of complex spectra, Phys

    G. Racah,Theory of complex spectra, Phys. Rev.63(1943) 367

    1943

  15. [15]

    Zerguine, P

    S. Zerguine, P. Van Isacker,Spin-aligned neutron-proton pairs inN=Znuclei, Phys. Rev. C 83(2011) 064314

    2011

  16. [16]

    Wybourne,Symmetry Principles and Atomic Spectroscopy(Wiley-Interscience, New York, 1970)

    B.G. Wybourne,Symmetry Principles and Atomic Spectroscopy(Wiley-Interscience, New York, 1970)

    1970

  17. [17]

    Weyl,The Classical Groups(Princeton University Press, Princeton, 1939)

    H. Weyl,The Classical Groups(Princeton University Press, Princeton, 1939)

    1939

  18. [18]

    Gheorghe, A.A

    A. Gheorghe, A.A. Raduta,New results for the missing quantum numbers labelling the quadrupole and octupole boson basisJ. Phys. A: Math. Gen.37(2004) 10951

    2004

  19. [19]

    Devi, V.K.B

    Y.D. Devi, V.K.B. Kota, Pramana J. Phys.39(1992) 413

    1992

  20. [20]

    Arima, F

    A. Arima, F. Iachello,Interacting boson model of collective nuclear states II. The rotational limit, Ann. Phys. (NY)111(1978) 201

    1978

  21. [21]

    Scholten, KVI Internal report no

    O. Scholten, KVI Internal report no. 63, 1979

    1979

  22. [22]

    Scholten,The Interacting Boson Approximation Model and Applications, Ph.D

    O. Scholten,The Interacting Boson Approximation Model and Applications, Ph.D. Thesis, University of Groningen, The Netherlands, 1980

    1980

  23. [23]

    Van Isacker, Fortran code ibm.f

    P. Van Isacker, Fortran code ibm.f

  24. [24]

    Garc´ ıa-Ramos, K

    J.E. Garc´ ıa-Ramos, K. Heyde,Nuclear shape coexistence: A study of the even–even Hg isotopes using the interacting boson model with configuration mixing, Phys. Rev. C89(2014) 014306

    2014

  25. [25]

    Heinze,Eine Methode zur L¨ osung beliebiger bosonischer und fermionischer Vielteilchensys- teme, Ph.D

    S. Heinze,Eine Methode zur L¨ osung beliebiger bosonischer und fermionischer Vielteilchensys- teme, Ph.D. Thesis, University of Cologne, Germany, 2008

    2008

  26. [26]

    Brown,The nuclear shell model towards the drip lines, Prog

    B.A. Brown,The nuclear shell model towards the drip lines, Prog. Part. Nucl. Phys.47(2001) 517

    2001

  27. [27]

    Brown, W.D.M

    B.A. Brown, W.D.M. Rae,The shell-model code NuShellX@MSU, Nucl. Data Sheets120 (2014) 115

    2014

  28. [28]

    Caurier, G

    E. Caurier, G. Mart´ ınez-Pinedo, F. Nowacki, A. Poves, A.P. Zuker,The shell model as a unified view of nuclear structure, Rev. Mod. Phys.77(2005) 427

    2005

  29. [29]

    Shimizu, T

    N. Shimizu, T. Abe, Y. Tsunoda, Y. Utsuno, T. Yoshida, T. Mizusaki, M. Honma, T. Otsuka, 30 New-generation Monte Carlo shell model for the K computer era, Prog. Theor. Exp. Phys. 2012, 1, 01A205

    2012

  30. [30]

    Thick-restart block Lanczos method for large-scale shell-model calculations.Comp

    Shimizu, N.; Mizusaki, T.; Utsuno, Y.; Tsunoda, Y. Thick-restart block Lanczos method for large-scale shell-model calculations.Comp. Phys. Comm.244(2019) 372

    2019

  31. [31]

    Johnson, W.E

    C.W. Johnson, W.E. Ormand, P.G. Krastev,Factorization in large-scale many-body calcula- tions, Comp. Phys. Comm.184(2013) 2761

    2013

  32. [32]

    Johnson, W.E

    C.W. Johnson, W.E. Ormand, K.S. McElvain, H. Shan,BIGSTICK: A flexible configuration- interaction shell-model code, arXiv:1801.08432

    work page arXiv

  33. [33]

    Van Isacker, Mathematica codeibm.m

    P. Van Isacker, Mathematica codeibm.m

  34. [34]

    Frank, P

    A. Frank, P. Van Isacker,Algebraic Methods in Molecular and Nuclear Structure Physics (Wiley-Interscience, New York, 1994)

    1994

  35. [35]

    Ginocchio, M.W

    J.N. Ginocchio, M.W. Kirson,Relationship between the Bohr collective Hamiltonian and the interacting-boson model, Phys. Rev. Lett.44(1980) 1744

    1980

  36. [36]

    Dieperink, O

    A.E.L. Dieperink, O. Scholten, F. Iachello,Classical limit of the interacting-boson model, Phys. Rev. Lett.44(1980) 1747

    1980

  37. [37]

    Bohr, B.R

    A. Bohr, B.R. Mottelson,Features of nuclear deformations produced by the alignment of in- dividual particles or pairs, Phys. Scripta22(1980) 468

    1980

  38. [38]

    Van Isacker, J.-Q

    P. Van Isacker, J.-Q. Chen,Classical limit of the interacting boson Hamiltonian, Phys. Rev. C24(1981) 684

    1981

  39. [39]

    Fortunato, C.E

    L. Fortunato, C.E. Alonso, J.M. Arias, J.E. Garc´ ıa-Ramos, A. Vitturi,Phase diagram for a cubic-Qinteracting boson model Hamiltonian: Signs of triaxiality, Phys. Rev. C84(2011) 014326

    2011

  40. [40]

    Poves, F

    A. Poves, F. Nowacki, Y. Alhassid,Limits on assigning a shape to a nucleus, Phys. Rev. C 101(2020) 054307

    2020

  41. [41]

    Warner, R.F

    D.D. Warner, R.F. Casten,Predictions of the interacting boson approximation in a consistent Qframework, Phys. Rev. C28(1983) 1798

    1983

  42. [42]

    Sorgunlu, P

    B. Sorgunlu, P. Van Isacker,Triaxiality in the interacting boson model, Nucl. Phys. A808 (2008) 27. 31

    2008