Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Physical Review

Eigenvalues of the Nuclear Pairing Hamiltonian

Richard R. Chasman

Published Oct 1, 1963• 17 citations
Publisher

Abstract

A new technique is developed in this paper for solving Hamiltonians of the form $H=\ensuremath{\Sigma}\stackrel{}{k>0}{\ensuremath{\epsilon}}_{k}({{a}_{k}}^{+}{a}_{k}+{{a}_{\ensuremath{-}k}}^{+}{a}_{\ensuremath{-}k})\ensuremath{-}\ensuremath{\Sigma}\stackrel{}{k,l>0}G(k,l){{a}_{k}}^{+}{{a}_{\ensuremath{-}k}}^{+}{a}_{\ensuremath{-}l}{a}_{l},$ where $G(k, l)$ may have the forms: $G(k, l)=G$ (independent of $k$ and $l$), $G(k, l)={G}_{k}{G}_{l}$, or $G(k, l)=G(k, l)$ (nonseparable in $k$ and $l$). The eigenvalues obtained with this technique are compared with the exact solutions of Kerman, Lawson, and Macfarlane. The agreement with the exact solutions is fairly good.

DOI
10.1103/physrev.132.343

Related papers

"Seniority-Zero" Excited States of the Pairing Hamiltonian

Richard R. Chasman

The seniority-zero excited-state spectrum of the Hamiltonian $H=\ensuremath{\Sigma}\stackrel{}{k>0}{\ensuremath{\epsilon}}_{k}({{a}_{k}}^{\ifmmode\dagger\else\textdagger\fi{}}{a}_{k}+{{a}_{\ensuremath{-}k}}^{\ifmmode\dagger\else\textdagger\fi{}}{a}_{\ensuremath{-}k})\ensuremath{-}G \ensuremath{\Sigma}\stackrel{}{k,l>0}{{a}_{k}}^{\ifmmode\dagger\else\textdagger\fi{}}{{a}_{\ensuremath{-}k}}^{\ifmmode\dagger\else\textdagger\fi{}}{a}_{\ensuremath{-}l}{a}_{l}$ is studied in detail. Approximate excited-state wave functions of the form ${\ensuremath{\psi}}^{\ensuremath{\lambda}}=\ensuremath{\Sigma}\stackrel{}{\ensuremath{\alpha}}F(\ensuremath{\lambda},\ensuremath{\alpha})\ensuremath{\Sigma}\stackrel{}{i}D(\ensuremath{\alpha},i){{a}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}{{a}_{\ensuremath{-}i}}^{\ifmmode\dagger\else\textdagger\fi{}}\ensuremath{\Sigma}\stackrel{}{j>i}D(\ensuremath{\alpha},j){{a}_{j}}^{\ifmmode\dagger\else\textdagger\fi{}}{{a}_{\ensuremath{-}j}}^{\ifmmode\dagger\else\textdagger\fi{}}\ensuremath{\cdots}|0〉$ are developed. These solutions are compared with exact solutions of a small system and the agreement is quite good. A somewhat larger system is studied in order to see how the excited-state spectrum changes as the number of separable functions is increased. The method is applied to heavy nuclei and is in good agreement with observed ${0}^{+}$ excited states for nuclei having 144 to 150 neutrons. A theoretical ${0}^{+}$ spectrum is displayed for each even neutron system from 144 to 152 neutrons.

Physical Review•Apr 26, 1965•20 citations

Solutions of nuclear pairing

A. B. Balantekin, Yamaç Pehlivan

We give the exact solution of orbit dependent nuclear pairing problem between two nondegenerate energy levels using the Bethe ansatz technique. Our solution reduces to previously solved cases in the appropriate limits including Richardson's treatment of reduced pairing in terms of rational Gaudin algebra operators.

arXiv:0710.3941•Physical Review C•Nov 14, 2007•18 citations

Study of Approximations in the Nuclear Pairing-Force Problem

Mordechai Bishari, I. Unna, A. Mann

The solvable two-level pairing model is used to test two new number-conserving approximations. These approximations are compared with the exact results, with perturbation theory, with the BCS method, and with the projected BCS method of Kerman, Lawson, and Macfarlane. It is found that second-order perturbation theory and one of the tested approximations work rather well for not too high values of the pairing-force strength, while for higher values the other approximation turns out to be quite good.

Physical Review C•Apr 1, 1971•13 citations

Approximation Techniques for Large Finite Quantum Many-body Systems

Shen Yong Ho

In this thesis, we will show how certain classes of quantum many-body Hamiltonians with $\\su{2}_1 \\oplus \\su{2}_2 \\oplus \\ldots \\oplus \\su{2}_k$ spectrum generating algebras can be approximated by multi-dimensional shifted harmonic oscillator Hamiltonians. The dimensions of the Hilbert spaces of such Hamiltonians usually depend exponentially on $k$. This can make obtaining eigenvalues by diagonalization computationally challenging. The Shifted Harmonic Approximation (SHA) developed here gives good predictions of properties such as ground state energies, excitation energies and the most probable states in the lowest eigenstates. This is achieved by solving only a system of $k$ equations and diagonalizing $k\\times k$ matrices. The SHA gives accurate approximations over wide domains of parameters and in many cases even across phase transitions.\n\nThe SHA is first illustrated using the Lipkin-Meshkov-Glick (LMG) model and the Canonical Josephson Hamiltonian (CJH) which have $\\su{2}$ spectrum generating algebras. Next, we extend the technique to the non-compact $\\su{1,1}$ algebra, using the five-dimensional quartic oscillator (5DQO) as an example. Finally, the SHA is applied to a $k$-level Bardeen-Cooper-Shrieffer (BCS) pairing Hamiltonian with fixed particle number. The BCS model has a $\\su{2}_1 \\oplus \\su{2}_2 \\oplus \\ldots \\oplus \\su{2}_k$ spectrum generating algebra.\n\nAn attractive feature of the SHA is that it also provides information to construct basis states which yield very accurate eigenvalues for low-lying states by diagonalizing Hamiltonians in small subspaces of huge Hilbert spaces. For Hamiltonians that involve a smaller number of operators, accurate eigenvalues can be obtained using another technique developed in this thesis: the generalized Rowe-Rosensteel-Kerman-Klein equations-of-motion method (RRKK). The RRKK is illustrated using the LMG and the 5DQO. In RRKK, solving unknowns in a set of $10\\times 10$ matrices typically gives estimates of the lowest few eigenvalues to an accuracy of at least eight significant figures. The RRKK involves optimization routines which require initial guesses of the matrix representations of the operators. In many cases, very good initial guesses can be obtained using the SHA. \n\nThe thesis concludes by exploring possible future developments of the SHA.

•Mar 3, 2010•0 citations

An iterative approach for the exact solution of the pairing Hamiltonian

Xin Guan, Chong Qi

Computer Physics Communications•Feb 17, 2022•6 citations

Approximate solutions for general pairing Hamiltonian by variational method

A. Kamimura

Lettere al nuovo cimento della societa italiana di fisica/Lettere al nuovo cimento•Mar 1, 1978•3 citations
A new technique is developed in this paper for solving Hamiltonians of the form H=\ensuremathΣk>0\ensuremathϵk(ak+ak+a\ensuremath−k+a\ensuremath−k)\ensuremath−\ensuremathΣk,l>0G(k,l)ak+a\ensuremath−k+a\ensuremath−lal,H=\ensuremath{\Sigma}\stackrel{}{k>0}{\ensuremath{\epsilon}}_{k}({{a}_{k}}^{+}{a}_{k}+{{a}_{\ensuremath{-}k}}^{+}{a}_{\ensuremath{-}k})\ensuremath{-}\ensuremath{\Sigma}\stackrel{}{k,l>0}G(k,l){{a}_{k}}^{+}{{a}_{\ensuremath{-}k}}^{+}{a}_{\ensuremath{-}l}{a}_{l},H=\ensuremathΣk>0​\ensuremathϵk​(ak​+ak​+a\ensuremath−k​+a\ensuremath−k​)\ensuremath−\ensuremathΣk,l>0​G(k,l)ak​+a\ensuremath−k​+a\ensuremath−l​al​, where G(k,l)G(k, l)G(k,l) may have the forms: G(k,l)=GG(k, l)=GG(k,l)=G (independent of kkk and lll), G(k,l)=GkGlG(k, l)={G}_{k}{G}_{l}G(k,l)=Gk​Gl​, or G(k,l)=G(k,l)G(k, l)=G(k, l)G(k,l)=G(k,l) (nonseparable in kkk and lll). The eigenvalues obtained with this technique are compared with the exact solutions of Kerman, Lawson, and Macfarlane. The agreement with the exact solutions is fairly good.
The seniority-zero excited-state spectrum of the Hamiltonian $H=\ensuremath{\Sigma}\stackrel{}{k>0}{\ensuremath{\epsilon}}_{k}({{a}_{k}}^{\ifmmode\dagger\else\textdagger\fi{}}{a}_{k}+{{a}_{\ensuremath{-}k}}^{\ifmmode\dagger\else\textdagger\fi{}}{a}_{\ensuremath{-}k})\ensuremath{-}G \ensuremath{\Sigma}\stackrel{}{k,l>0}{{a}_{k}}^{\ifmmode\dagger\else\textdagger\fi{}}{{a}_{\ensuremath{-}k}}^{\ifmmode\dagger\else\textdagger\fi{}}{a}_{\ensuremath{-}l}{a}_{l}$ is studied in detail. Approximate excited-state wave functions of the form ${\ensuremath{\psi}}^{\ensuremath{\lambda}}=\ensuremath{\Sigma}\stackrel{}{\ensuremath{\alpha}}F(\ensuremath{\lambda},\ensuremath{\alpha})\ensuremath{\Sigma}\stackrel{}{i}D(\ensuremath{\alpha},i){{a}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}{{a}_{\ensuremath{-}i}}^{\ifmmode\dagger\else\textdagger\fi{}}\ensuremath{\Sigma}\stackrel{}{j>i}D(\ensuremath{\alpha},j){{a}_{j}}^{\ifmmode\dagger\else\textdagger\fi{}}{{a}_{\ensuremath{-}j}}^{\ifmmode\dagger\else\textdagger\fi{}}\ensuremath{\cdots}|0〉$ are developed. These solutions are compared with exact solutions of a small system and the agreement is quite good. A somewhat larger system is studied in order to see how the excited-state spectrum changes as the number of separable functions is increased. The method is applied to heavy nuclei and is in good agreement with observed 0+{0}^{+}0+ excited states for nuclei having 144 to 150 neutrons. A theoretical 0+{0}^{+}0+ spectrum is displayed for each even neutron system from 144 to 152 neutrons.
The solvable two-level pairing model is used to test two new number-conserving approximations. These approximations are compared with the exact results, with perturbation theory, with the BCS method, and with the projected BCS method of Kerman, Lawson, and Macfarlane. It is found that second-order perturbation theory and one of the tested approximations work rather well for not too high values of the pairing-force strength, while for higher values the other approximation turns out to be quite good.
In this thesis, we will show how certain classes of quantum many-body Hamiltonians with su21oplussu22oplusldotsoplussu2k\\su{2}_1 \\oplus \\su{2}_2 \\oplus \\ldots \\oplus \\su{2}_ksu21​oplussu22​oplusldotsoplussu2k​ spectrum generating algebras can be approximated by multi-dimensional shifted harmonic oscillator Hamiltonians. The dimensions of the Hilbert spaces of such Hamiltonians usually depend exponentially on kkk. This can make obtaining eigenvalues by diagonalization computationally challenging. The Shifted Harmonic Approximation (SHA) developed here gives good predictions of properties such as ground state energies, excitation energies and the most probable states in the lowest eigenstates. This is achieved by solving only a system of kkk equations and diagonalizing ktimeskk\\times kktimesk matrices. The SHA gives accurate approximations over wide domains of parameters and in many cases even across phase transitions.\n\nThe SHA is first illustrated using the Lipkin-Meshkov-Glick (LMG) model and the Canonical Josephson Hamiltonian (CJH) which have su2\\su{2}su2 spectrum generating algebras. Next, we extend the technique to the non-compact su1,1\\su{1,1}su1,1 algebra, using the five-dimensional quartic oscillator (5DQO) as an example. Finally, the SHA is applied to a kkk-level Bardeen-Cooper-Shrieffer (BCS) pairing Hamiltonian with fixed particle number. The BCS model has a su21oplussu22oplusldotsoplussu2k\\su{2}_1 \\oplus \\su{2}_2 \\oplus \\ldots \\oplus \\su{2}_ksu21​oplussu22​oplusldotsoplussu2k​ spectrum generating algebra.\n\nAn attractive feature of the SHA is that it also provides information to construct basis states which yield very accurate eigenvalues for low-lying states by diagonalizing Hamiltonians in small subspaces of huge Hilbert spaces. For Hamiltonians that involve a smaller number of operators, accurate eigenvalues can be obtained using another technique developed in this thesis: the generalized Rowe-Rosensteel-Kerman-Klein equations-of-motion method (RRKK). The RRKK is illustrated using the LMG and the 5DQO. In RRKK, solving unknowns in a set of 10times1010\\times 1010times10 matrices typically gives estimates of the lowest few eigenvalues to an accuracy of at least eight significant figures. The RRKK involves optimization routines which require initial guesses of the matrix representations of the operators. In many cases, very good initial guesses can be obtained using the SHA. \n\nThe thesis concludes by exploring possible future developments of the SHA.