# Presentation: 350413

## 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 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

**Lecture #3**

- 26-Jun-01

## Electro-magnetic transitions

**Ok, now what do I do with the states?**
**Well, for one excited states decay!**
**Electromagnetic decays**
**Electric multipole (EL)**
** ********J given by |J****i****-L| ******** J****f ******** J****i****+L**
**Parity change (-1)****L**
**Magnetic multipole (ML)**
** ********J given by |J****i****-L| ******** J****f ******** J****i****+L **
**Parity change (-1)****L+1**

- 26-Jun-01

## Transition life-times

**Define the “B”-values**
**The transition rate is**
**Note the important phase-space factor, *****E************2L******+******1***

- 26-Jun-01

## Electromagnetic transitions

**How well does the shell-model work?**
**Not well at all with free electric charges!**
**Ok, with free g-factors**
**So, where did we do wrong?????!!!!!**
**Remember we renormalized the interaction**
**This accounts for excitations not included in the active valence space**
**What about the operators?**
**We also have to renormalize the transition operators!**
***e******p*****~1.3 and *****e******n*****~0.5**
**Free g-factors for M1 transitions are not bad (but some renormalization is needed like adding ****[********xY****2****]****1****)**
**Only with these renormalized (effective) operators, we can get excellent agreement with experiment**

- 26-Jun-01

## Estimates for electromagnetic transitions

**Weisskopf estimates**
**Assume ******** constant over nuclear volume, zero outside**

- 26-Jun-01

## Estimates for electromagnetic transitions

**Use the Weisskopf estimates to determine**

- 26-Jun-01

## How big are nuclei?

**Electron scattering**
**Current-current interaction**
**Charge form factor**

- 26-Jun-01

## How big are nuclei?

**Electron scattering**

- 26-Jun-01

## Is there anyway to probe the neutrons?

**Yes, again with electron scattering**
**But we must look to the parity violating part**
**Neutral current - Z-boson!**
**Parity-violating electron scattering also provides a test of the Standard Model**

- 26-Jun-01

## The Signal

**PV elastic electron scattering:**
**Charges: **
**Include intrinsic electric and magnetic form factors**
**Neutral current couples more strongly to the neutron distribution**

- 26-Jun-01

## Nuclear Structure effects

**With isospin symmetry**
**No nuclear structure effects **
**low *****q***** ******** 1% measurement of sin****2****q****W**
**Deviations are a signature for “new” Physics**
**Exotic neutral currents**
**Strangeness form factor, *****F******C*****(s)**** (higher q)**

**Deviations from *****q******2***** behavior could signal new Physics**
**or be due to Nuclear Physics**

- 26-Jun-01

## Nuclear Structure effects

**Isospin-symmetry is broken**
**Coulomb interaction (larger)**
**strong interaction (smaller): isotensor or charge-dependent interaction**
**v****(pn) ****- (v****(pp)**** + v****(nn)****)/2 **
**Mix states with ********T****max****=2**
**Effect on observables:**

- 26-Jun-01

## Nuclear Structure effects

**Correction due to breaking of isospin symmetry**
**Overall agreement with recent *****ab initio***** calculation (Navratil and Barrett)**
** ****G(*****q*****) < 1% for *****q***** < 0.9 fm****-1**** (1% measurement for 0.3 < *****q***** < 1.1 fm****-1****,  Musolf & Donnelly, NPA546, 509 (1992).**

******(q) < 1% for q < 0.9 fm****-1**** and q=2.4± 0.1 fm****-1**

- 26-Jun-01

## Neutron Radii

**What good is parity violation?**
**Assume Standard Model correct to 1% level - infer neutron distribution**
**Very little precision data regarding the distribution of neutrons**
**Useful for mean-field models - improve extrapolation to the drip line**
**Hadron scattering - strong in-medium effects**

**Experiments planned for ****208****Pb at TJNAF**

- 26-Jun-01

## The weak interaction in the shell model

** ********-decay and neutrino absorption**
** ********-decay**
**Partial half-life**
**Fermi (F): **
**Gamow-Teller (GT): **
**g****A****=1.2606 ******** 0.0075**
**GT is very dependent on model space and shell-model interaction**
**Spin-orbit and quasi SU(4) symmetry**
**Meson-exchange currents modify B(GT)**
**For an effective operator, GT must be renormalized: multiply by ~ 0.75**
**Total half life:                                        Branching ratio:**

- 26-Jun-01

## The weak interaction in the shell modelIsospin-symmetry violation

**Isospin is approximately conserved (~ 1% level)**
**For transitions, isospin violation enters in two places**
**One-body transition density as ******** no longer has good isospin**
**One-body matrix element**
**Note we have a proton(neutron) converted to a neutron(proton)**
**Due to the Coulomb interaction protons and neutrons have different radial wave functions, so we need the overlap: **
**Important for high-precision tests of the vector current (~0.4%)**
**For GT this effect can be large, and is essentially contained in the global factor of 0.75 obtained empirically**
**Mirror transitions are no longer the same!!!**

- 26-Jun-01

## Superallowed Fermi -decay

**Test of the Standard Model**
**Cabibbo-Kobayashi-Maskawa matrix**

- 26-Jun-01

## Superallowed Fermi -decay

- C for 10C: 0.15(9)%

- Ab initio calculation

**Unitarity condition: 0.9956(8)****stat****(7)****sys**

- Current Situation

- 26-Jun-01

## Neutrino absorption

**As for ********-decay, neutrino absorption requires Fermi and Gamow-Teller matrix elements**
**We carried out calculations for ****23****Na and ****40****Ar**
**Ormand et al., PLB308, 207 (1993), Ormand et al., PLB345, 343 (1995)**
**ICARUS and proposed bolometric detectors**
**For ****40****Ar Gamow-Teller is very important**
**Total is twice as large as Fermi contribution**
**Counter to original design assumption**
**Can we trust the calculation?**
** ********-decay of analog**

- 26-Jun-01

## Checking the calculation

**For ****40****Ar, look at ********-decay of ****40****Ti**
** ********-delayed proton emitter**
**Calculated half-life: 55 ******** 5 ms**
**Exp: 52.7 ******** 1.5 ms and 54 ******** 2 ms**

- 26-Jun-01

## Checking the calculation

**But there are problems with B(GT) strength**

- Theory:=11.50.710-43 cm

**There is no substitute for experiment when available**

- 26-Jun-01

## What about heavier nuclei?

**Above A ~ 60 or so the number of configurations just gets to bed too large ~ 10****10****!**
**Here, we need to think of more approximate methods**
**The easiest place to start is the mean-field of Hartree-Fock**
**But, once again we have the problem of the interaction**
**Repulsive core causes us no end of grief!!**
**We will, at some point use effective interactions like the Skyrme force**

- 26-Jun-01

## Hartree-Fock

**There are many choices for the mean field, and Hartree-Fock is one optimal choice**
**We want to find the best single Slater determinant ********0**** ****so that**
**Thouless’ theorem**
**Any other Slater determinant ******** not orthogonal to ********0**** may be written as**
**Where *****i***** is a state occupied in ********0**** and *****m***** is unoccupied**
**Then**

- 26-Jun-01

## Hartree-Fock

**Let *****i,j,k,l***** denote occupied states and *****m,n,o,p***** unoccupied states**
**After substituting back we get**
**This leads directly to the Hartree-Fock single-particle Hamiltonian *****h***** with matrix elements between any two states ********** and ********

- 26-Jun-01

## Hartree-Fock

**We now have a mechanism for defining a mean-field**
**It does depend on the occupied states**
**Also the matrix elements with unoccupied states are zero, so the first order 1p-1h corrections do not contribute**
**We obtain an eigenvaule equation (more on this later)**
**Energies of A+1 and A-1 nuclei relative to A**

- 26-Jun-01

## Hartree-Fock – Eigenvalue equation

**Two ways to approach the eigenvalue problem**
**Coordinate space where we solve a Schrödinger-like equation**
**Expand in terms of a basis, e.g., harmonic-oscillator wave function**
**Expansion**
**Denote basis states by Greek letters, e.g., ********
**From the variational principle, we obtain the eigenvalue equation**

- 26-Jun-01

## Hartree-Fock – Solving the eigenvalue equation

**As I have written the eigenvalue equation, it doesn’t look to useful because we need to know what states are occupied**
**We use three steps**
**Make an initial guess of the occupied states and the expansion coefficients *****C******i*********
**For example the lowest Harmonic-oscillator states, or a Woods-Saxon and *****C******i***********=***********i*********
**With this ansatz, set up the eigenvalue equations and solve them**
**Use the eigenstates ****|*****i********* from step 2 to make the Slater determinant ********0****, go back to step 2 until the coefficients *****C******i************ *****are unchanged**

**The Hartree-Fock equations are solved self-consistently**

- 26-Jun-01

## Hartree-Fock – Coordinate space

**Here, we denote the single-particle wave functions as ***********i*****(r) **
**These equations are solved the same way as the matrix eigenvalue problem before**
**Make a guess for the wave functions ***********i*****(r) and Slater determinant ********0**
**Solve the Hartree-Fock differential equation to obtain new states ***********i*****(r) **
**With these go back to step 2 and repeat until ***********i*****(r) are unchanged**

- Direct or Hartree term: UH

- Exchange or Fock term: UF

**Again the Hartree-Fock equations are solved self-consistently**

- 26-Jun-01

## Hartree-Fock

**M. Moshinsky, Am. J. Phys. 36, 52 (1968). Erratum, Am. J. Phys. 36, 763 (1968).**
**Two identical spin-1/2 particles in a spin singlet interact via the Hamiltonian**
**Use the coordinates                      and                        to show the exact energy and wave function are**
**Note that since the spin wave function (S=0) is anti-symmetric, the spatial wave function is symmetric**

**Hard homework problem:**

- 26-Jun-01

## Hartree-Fock

**The Hartree-Fock solution for the spatial part is the same as the Hartree  solution for the S-state. Show the Hartree energy and radial wave function are: **

**Hard homework problem:**

- 26-Jun-01

## Hartree-Fock with the Skyrme interaction

**In general, there are serious problems trying to apply Hartree-Fock with realistic NN-interactions (for one the saturation of nuclear matter is incorrect)**
**Use an effective interaction, in particular a force proposed by Skyrme**
***P*********** is the spin-exchange operator**
**The three-nucleon interaction is actually a density dependent two-body, so replace with a more general form, where ********** determines the incompressibility of nuclear matter**

- 26-Jun-01

## Hartree-Fock with the Skyrme interaction

**One of the first references: D. Vautherin and D.M. Brink, PRC5, 626 (1972)**
**Solve a Shrödinger-like equation**
**Note the effective mass *****m************* ***
**Typically, *****m************* < m*****, although it doesn’t have to, and is determined by the parameters *****t******1***** and *****t******2***
**The effective mass influences the spacing of the single-particle states**
**The bias in the past was for *****m*************/m***** ~ 0.7 because of earlier calculations with realistic interactions**

***z* labels protons or neutrons

- 26-Jun-01

## Hartree-Fock calculations

**The nice thing about the Skyrme interaction is that it leads to a computationally tractable problem**
**Spherical (one-dimension)**
**Deformed**
**Axial symmetry (two-dimensions)**
**No symmetries (full three-dimensional)**
**There are also many different choices for the Skyrme parameters**
**They all do some things right, and some things wrong, and to a large degree it depends on what you want to do with them**
**Some of the leading (or modern) choices are:**
**M*********, M. Bartel *****et al*****., NPA386, 79 (1982)**
**SkP [includes pairing], J. Dobaczewski and H. Flocard, NPA422, 103 (1984)**
**SkX, B.A. Brown, W.A. Richter, and R. Lindsay, PLB483, 49 (2000)**
**Apologies to those not mentioned!**
**There is also a finite-range potential based on Gaussians due to D. Gogny, D1S, J. Dechargé and D. Gogny, PRC21, 1568 (1980).**
**Take a look at J. Dobaczewski *****et al*****., PRC53, 2809 (1996) for a nice study near the neutron drip-line and the effects of unbound states**

- 26-Jun-01

## Nuclear structure

**Remember what our goal is:**
**To obtain a quantitative description of all nuclei within a microscopic frame work**
**Namely, to solve the many-body Hamiltonian:**

- Residual interaction

- Perturbation Theory

- 26-Jun-01

## Nuclear structure

**Hartree-Fock is the optimal choice for the mean-field potential *****U(r)*****!**
**The Skyrme interaction is an “effective” interaction that permits a wide range of studies, e.g., masses, halo-nuclei, etc.**
**Traditionally the Skyrme parameters are fitted to binding energies of doubly magic nuclei, rms charge-radii, the incompressibility, and a few spin-orbit splittings**
**One goal would be to calculate masses for all nuclei**
** ****By fixing the Skyrme force to known nuclei, maybe we can get 500 keV accuracy that CAN be extrapolated into the unknown region**
**This will require some input about neutron densities – parity-violating electron scattering can determine <r****2****>****p****-<r****2****>****n****. **
**This could have an important impact**

- 26-Jun-01

## Hartree-Fock calculations

**Hartree-Fock calculations**

**Permits a study of a wide-range of nuclei, in particular, those far from stability and with exotic properties, halo nuclei**

- H. Sagawa, PRC65, 064314 (2002)

- The tail of the radial density depends on the separation energy
- S. Mizutori et al. PRC61, 044326 (2000)

- Drip-line studies
- J. Dobaczewski *et al*., PRC53, 2809 (1996)

- 26-Jun-01

## What can Hartree-Fock calculations tell us about shell structure?

**Shell structure**
**Because of the self-consistency, the shell structure can change from nucleus to nucleus**

**J. Dobaczewski *****et al*****., PRC53, 2809 (1996)**

**As we add neutrons, traditional shell closures are changed, and may even disappear! **
**This is **_**THE**_** challenge in trying to predict the structure of nuclei at the drip lines!**

- 26-Jun-01

## Beyond mean field

**Hartee-Fock is a good starting approximation**
**There are no particle-hole corrections to the HF ground state**
**The first correction is**
**However, this doesn’t make a lot of sense for Skyrme potentials**
**They are fit to closed-shell nuclei, so they effectively have all these higher-order corrections in them!**
**We can try to estimate the excitation spectrum of one-particle-one-hole states – Giant resonances**
**Tamm-Dancoff approximation (TDA)**
**Random-Phase approximation (RPA)**

- You should look these up!
- A Shell Model Description of Light Nuclei, I.S. Towner
- The Nuclear Many-Body Problem, Ring & Schuck

- 26-Jun-01

## Nuclear structure in the future

**Nuclear structure in the future**

**With newer methods and powerful computers, the future of nuclear structure theory is bright!**

- 26-Jun-01