# Presentation: 947707

## Nuclear Structure

**Erich Ormand**
**N-Division, Physics and Adv. Technologies Directorate**
**Lawrence Livermore National Laboratory**

**This work was carried out under the auspices of the U.S. Department of Energy by the University of California, Lawrence Livermore National Laboratory under contract No. W-7405-Eng-48. This work was supported in part from LDRD contract 00-ERD-028 **

- This document was prepared as an account of work sponsored by an agency of  the United States Government. Neither the United States Government nor the University of California nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise, does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or the University of California. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or the University of California, and shall not be used for advertising or product endorsement purposes.

- 15thNational Nuclear Physics
- Summer School
- June 15-27, 2003

- 26-Jun-01

## Nuclear Structure

**At least we speak a common language. It could be worse:**

**Nuclear Structure**

- 26-Jun-01

## Nuclear Structure 

**The nuclear many-body problem is one of the hardest problems in all of physics!**
**Do we really know how they are put together?**
**This is a fundamental question in nuclear physics and we are now getting some interesting answers – for example, three-nucleon forces are important for structure**
**Atomic nuclei make up the vast majority of matter that we can see (and touch). How did they (and we) get there?**
**One of the key questions in science**
**Nucleosynthesis **
**Structure plays an important role**

**Nuclear Structure **

- 26-Jun-01

## Nuclear Structure 

**Computationally exact methods for A up to 16**
**Approximate many-body methods for A up to 60**
**Mostly mean-field pictures for A greater than 60 or so**

**Nuclear Structure **

- 26-Jun-01

## Nuclear Structure 

**So, there will be formulae**

**Nuclear Structure **

- 26-Jun-01

## Before we get started: Some useful reference materials I

***Structure of the Nucleus*****, M.A. Preston and R.K. Bhaduri, (Addison-Wesley,Reading, MA, 1975)**
***Nuclear Models*****, W. Greiner and J.A. Maruhn, (Springer Verlag, Berlin, 1996)**
***Basic Ideas and Concepts in Nuclear Physics*****, K. Heyde (IoP Publishing, Bristol, 1999)**
***Nuclear Structure vols. I & II*****, A. Bohr and B. Mottelson, (W.A. Benjamin, New York, 1969) & **
***Nuclear Theory, vols. I-III*****, J.M. Eisenberg and W. Greiner, (North Holland, Amsterdam, 1987)**

**Before we get started: ****Some useful reference materials I**

- General references for nuclear-structure physics:

- 26-Jun-01

## Before we get started: Some useful reference materials II

***The Nuclear Many-Body problem, *****P. Ring and P. Schuck, (Springer Verlag, Berlin, 1980)**
***Theory of the Nuclear Shell Model*****, R.D. Lawson, (Clarendon Press, Oxford, 1980)**
***A Shell Model Description of Light Nuclei, *****I.S. Towner, (Clarendon Press, Oxford, 1977)**
***The Nuclear Shell Model Towards the Drip Line*****, B.A. Brown, Progress in Particle and Nuclear Physics 47, 517 (2001)**

**Before we get started: ****Some useful reference materials II**

- References for many-body problem:

- Review for applying the shell model near the drip lines:

- 26-Jun-01

## Nuclear masses, what nuclei exist?

**There are global Volume, Surface,  Symmetry, and Coulomb terms**
**And specific corrections for each nucleus due to pairing and shell structure**

**Nuclear masses, what nuclei exist?**

**One goal in theory is to accurately describe the binding energy**

- Values for the parameters,
- A.H. Wapstra and N.B. Gove, Nulc. Data Tables 9, 267 (1971)

- 26-Jun-01

## Nuclear masses, what nuclei exist?

**Liquid-drop isn’t too bad!**
**There are notable problems though.**
**Can we do better and what about the microscopic structure?**

- Bethe-Weizsäcker formula

- From A. Bohr and B.R. Mottleson,
*Nuclear Structure*, vol. 1, p. 168
- Benjamin, 1969, New York

- 26-Jun-01

## Nuclear masses, what nuclei exist?

**How about deformation?**
**For each energy term, there are also shape factors dependent on the quadrupole deformation parameters ********** and ********

**Nuclear masses, what nuclei exist?**

**Note that the liquid drop always has a minimum for a spherical shape!**
**So, where does deformation come from?**

- 26-Jun-01

## Shell structure - evidence in atoms

**Atomic ionization potentials show sharp discontinuities at shell boundaries**

- From A. Bohr and B.R. Mottleson,
*Nuclear Structure*, vol. 1, p. 191
- Benjamin, 1969, New York

- 26-Jun-01

## Shell structure - neutron separation energies

**So do neutron separation energies**

- From A. Bohr and B.R. Mottleson,
*Nuclear Structure*, vol. 1, p. 193
- Benjamin, 1969, New York

- 26-Jun-01

## More evidence of shell structure

**2, 8, 20, 28, 50, 82, and 126**

**More evidence of shell structure**

**Experimental Single-particle Effect (MeV)**

- From W.D. Myers and W.J. Swiatecki, Nucl. Phys. **81**, 1 (1966)

- Deviation from average binding energy

- 26-Jun-01

## Origin of the shell model

**Goeppert-Mayer and Haxel, Jensen, and Suess proposed the independent-particle shell model to explain the magic numbers**

**Origin of the shell model**

- M.G. Mayer and J.H.D. Jensen, *Elementary Theory of Nuclear *
*Shell Structure*,p. 58, Wiley, New York, 1955

**Harmonic oscillator with spin-orbit is a reasonable approximation to the nuclear mean field**

- 26-Jun-01

## Nilsson Hamiltonian - Poor man’s Hartree-Fock

**Nilsson Hamiltonian - Poor man’s Hartree-Fock**

- Prolate, =0 ̊

- Oblate, =60 ̊

- From J.M. Eisenberg and W. Greiner, *Nuclear Models*,
-  p 542, North Holland, Amsterdam, 1987

- 26-Jun-01

## Nuclear massesShell corrections to the liquid drop

**In general, the liquid drop does a good job on the bulk properties**
**The oscillator doesn’t!**
**But we need to put in corrections due to shell structure **
**Strutinsky averaging; difference between the energy of the discrete spectrum and the averaged, smoothed spectrum **

**Nuclear masses****Shell corrections to the liquid drop**

- Discrete spectrum

- Smoothed spectrum

- 26-Jun-01

## Nilsson-Strutinsky and deformation

**Nilsson-Strutinsky and deformation**

**Nilsson-Strutinsky is a mean-field type approach that allows for a comprehensive study of nuclear deformation under rotation and at high temperature**

- From C.G. Andersson,et al., Nucl. Phys. **A268**, 205 (1976)

- 26-Jun-01

## Nuclear massesPairing corrections to the liquid drop

**Nuclear masses****Pairing corrections to the liquid drop**

- From A. Bohr and B.R. Mottleson,
*Nuclear Structure*, vol. 1, p. 170
- Benjamin, 1969, New York

- 26-Jun-01

## How well do mass formulae work?

**How well do mass formulae work?**

**Predictions towards the neutron-drip line tend to diverge!**
**Why?**

**Mass predictions for Sn isotopes **
**and experiment**

- 26-Jun-01

## How well do mass formulae work?

**How well do mass formulae work?**

**The location of the neutron drip-line in rather uncertain!**

- 26-Jun-01

## Mass Formulae

**Show that the most stable Z****0**** value is**

**Mass Formulae**

**Homework:**

**More advanced homework:**

**Assume symmetric fission and calculate the energy released. Approximately at what A value is the energy released > 0?**

- 26-Jun-01

## What are the forces in nuclei?

**NN-scattering**
**The deuteron**
**From these we infer the form of the nucleon-nucleon interaction**
**The starting point is, of course, the Yukawa hypothesis of meson exchange**
**Pion, rho, sigma, two pion, etc.**
**However, it is also largely phenomenological**
**Deuteron binding energy: 2.224 MeV**
**Deuteron quadrupole moment: 0.282 fm****2**
**Scattering lengths and ranges for pp, nn, and analog pn channels**
**Unbound!**
**Note that V****pp ******** V****nn**** ******** V****pn**
**Some of the most salient features are the Tensor force and a strong repulsive core at short distances**

**What are the forces in nuclei?**

- 26-Jun-01

## NN-interactions 

**R.B. Wiringa, V.G.J. Stoks, R. Schiavilla, PRC51, 38 (1995)**
**Coulomb + One pion exchange + intermediate- and short-range**
**Bonn potential**
**R. Machleidt, PRC63, 024001 (2001)**
**Based on meson-exchange**
**Non-local**
**Effective field theory**
**C. Ordóñez, L. Ray, U. van Kolck, PRC53, 2086 (1996); E. Epelbaoum, W. Glöckle, Ulf-G. Meißner, NPA637, 107 (1998)**
**Based on Chiral Lagrangians**
**Expansion in momentum relative to a cutoff parameter (~ 1 GeV)**
**Generally has a soft core**
**All are designed to reproduce the deuteron and NN-scattering**

**NN-interactions **

- 26-Jun-01

## NN-interactions

**Elastic scattering in momentum space**
**Or, through a Fourier transform, coordinate space (                  )**
**Off-shell component  present in the Bonn potentials**
**Non-local (depends on the energies of the initial and final states)**
**Plays a role in many-body applications and provides more binding**

**NN-interactions**

- Tensor operator

- 26-Jun-01

## Three-Nucleon interactions

**Two-nucleon interactions under bind**
**Note CD-Bonn has a little more binding due to non-local terms**
**Further evidence is provided by *****ab initio***** calculations for ****10****B**
**NN-interactions give the wrong ground-state spin!**
**More on this later**

**Three-Nucleon interactions**

**Tucson-Melbourne**
**S.A. Coon and M.T. Peña, PRC48, 2559 (1993)**
**Based on two-pion exchange and intermediate ********’s**
**The exact form of NNN is not known**

**There is mounting evidence that three-body forces are very important**

- 26-Jun-01

## Isospin

**Nearly the same mass, *****q******p*****=1, *****q******n*****=0**
**Heisenberg introduced a spin-like quantity with the z-component defining the electric charge**
**Protons and neutrons from an isospin doublet**
**Add isospin using angular momentum algebra, e.g., two particles:**
**T=1**
**T=0**

**Isospin**

**With T=0, symmetry under p ******** n**

- 26-Jun-01

## Isospin

**Even-even *****N=Z*****: *****T=*****0**
***N************ Z*****: *****T=T******z***
**Odd-odd *****N=Z*****: *****T=*****0 ****or**** T=1 ****(Above A=22, essentially degenerate)**
**If V****pp****=V****nn****=V****pn****, isospin-multiplets have the same energy and isospin is a good quantum number**

**Isospin**

- 26-Jun-01

## Isospin

**NN-interaction has scalar, vector, and tensor components in isospin space**
**Note that the Coulomb interaction contributes to each component and is the largest!!!**

**Isospin**

- 26-Jun-01

## Coulomb-displacement energies

**We apply the Wigner-Eckart theorem**

**Coulomb-displacement energies**

- 26-Jun-01

## Isobaric-Mass-Multiplet Equation

**Isobaric-Mass-Multiplet Equation**

- 26-Jun-01

## Coulomb-displacement energies

**Can we use the IMME to predict the proton drip-line?**
**Binding energy difference between mirror nuclei**
**In a T-multiplet, often the binding energy of the neutron-rich mirror is measured**

**Coulomb-displacement energies**

- Calculated with theory
- 30 keV uncertainty

- Simple estimate from a charged sphere with radius r=1.2A1/2

- 26-Jun-01

## Coulomb-displacement energies

**B.A. Brown, PRC42, 1513 (1991); W.E. Ormand, PRC53, 214 (1996); B.J. Cole, PRC54, 1240 (1996); W.E. Ormand, PRC55, 2407 (1997); B.A. Brown et al., PRC65, 045802 (2002)**

**Coulomb-displacement energies**

**Yes! **
**Coulomb displacement energies provide an accurate method to map the proton drip line up to A=71**

- Use nature nature to give us the strong interaction part, i.e., the *a*-coefficient by adding the Coulomb-displacement to the experimental binding energy of the neutron-rich mirror

- 26-Jun-01

## Now, lets start to get at the structure of nuclei

**For three and four, we turn to Faddeev and Faddeev-Yakubovsky formulations**

**Now, lets start to get at the structure of nuclei**

**Exact methods exist for A ******** 4**

**Jacobi coordinates**

**Identical particles**

- W. Glöcle in *Computational Nuclear Physics*,
- Springer-Verlag, Berlin, 1991

- 26-Jun-01

## Three-particle harmonic oscillator

**Three-particle harmonic oscillator**

**Harmonic oscillators are always a convenient way to start solving the problem**

**Jacobi coordinates**

**Show that we can rewrite *****H***** as:**

**Finally, we have to couple L to the spin S and anti-symmetrize**

- 26-Jun-01

## More than four particles

**More than four particles**

**This is where life starts to get very hard!**
**Why?**
**Because there are so many degrees of freedom.**
**What do we do?**
**Green’s Function Monte Carlo**
**Coupled-cluster**
**Shell model**
**I’ll tell you about this approach**

- 26-Jun-01

## Many-body Hamiltonian

**Start with the many-body Hamiltonian**
**Introduce a mean-field *****U*****  ****to yield basis**
**The mean field determines the shell structure**
**In effect, nuclear-structure calculations rely on perturbation theory**

- Residual interaction

**The success of any nuclear structure calculation depends on the choice of the mean-field basis and the residual interaction!**** **

- 26-Jun-01

## Single-particle wave functions

**With the mean-field, we have the basis for building many-body states**
**This starts with the single-particle, radial wave functions, defined by the radial quantum number *****n*****, orbital angular momentum *****l*****, and z-projection *****m***
**Now include the spin wave function:**
**Two choices, jj-coupling or ls-coupling**
**Ls-coupling**
**jj-coupling is very convenient when we have a spin-orbit (*****l************s*****) force**

- 26-Jun-01

## Multiple-particle wave functions

**Total angular momentum, and isospin; **
**Anti-symmetrized, two particle, jj-coupled wave function**
**Note *****J+T*****=odd if the particles occupy the same orbits**
**Anti-symmetrized, two particle, LS-coupled wave function**

- 26-Jun-01

## Two-particle wave functions

**Of course, the two pictures describe the same physics, so there is a way to connect them**
**Recoupling coefficients**
**Note that the wave functions have been defined in terms of     and     , but often we need them in terms of the relative coordinate **
**We can do this in two ways**
**Transform the operator**

**Quadrupole, l=2, component is large and very important**

- 26-Jun-01

## Two-particle wave functions in relative coordinate

**Use Harmonic-oscillator wave functions and decompose in terms of the relative and center-of-mass coordinates, i.e.,**
**Harmonic oscillator wave functions are a very good approximation to the single-particle wave functions**
**We have the useful transformation**
**2*****n*****1*****+l*****2*****+2n*****2*****+l*****2*****=*****2*****n+l+2N+L***
**Where the *****M(nlNL;n*****1*****l*****1*****n*****2*****l*****2*****)***** is known as the Moshinksy bracket **
**Note this is where we use the jj to LS coupling transformation**
**For some detailed applications look in *****Theory of the Nuclear Shell Model*****, R.D. Lawson, (Clarendon Press, Oxford, 1980)**

- 26-Jun-01

## Many-particle wave function

**To add more particles, we just continue along the same lines**
**To build states with good angular momentum, we can bootstrap up from the two-particle case, being careful to denote the distinct states**
**This method uses Coefficients of Fractional Parentage (CFP)**
**Or we can make a many-body Slater determinant that has only a specified *****J******z***** and *****T******z***** and project *****J***** and *****T***

**In general Slater determinants are more convenient**

- 26-Jun-01

## Second Quantization

**Second quantization is one of the most useful representations in many-body theory**
**Creation and annhilation operators**
**Denote |0******** as the state with no particles (the vacuum)**
***a*****+*****i******  *****creates a particle in state *****i;***
***a******i****** *****annhilates a particle in state *****i;***
**Anticommuntation relations:**
**Many-body Slater determinant **

- 26-Jun-01

## Second Quantization

**Operators in second-qunatization formalism**
**Take any one-body operator *****O*****, say quadrupole E2 transition operator *****er******2******Y******2***********, the operator is represented as:**
**where *********j|O|i********* is the single-particle matrix element of the operator *****O***
**The same formalism exists for any *****n*****-body operator, e.g., for the NN-interaction**
**Here, I’ve written the two-body matrix element with an implicit assumption that it is anti-symmetrized, i.e.,**

- 26-Jun-01

## Second Quantization

**Matrix elements for Slater determinants**

**Second quantization makes the computation of expectation values for the many-body system simpler**

- 26-Jun-01

## Second Quantization

**Angular momentum tensors**
**Creation operators rotate as tensors of rank *****j***
**Not so for annihilation operators**
**Anti-symmetrized, two-body state**

- 26-Jun-01

## The mean field

**One place to start for the mean field is the harmonic oscillator**
**Specifically, we add the center-of-mass potential**
**The Good:**
**Provides a convenient basis to build the many-body Slater determinants**
**Does not affect the intrinsic motion**
**Exact separation between intrinsic and center-of-mass motion**
**The Bad:**
**Radial behavior is not right for large *****r ***
**Provides a confining potential, so all states are effectively bound**

- 26-Jun-01