The Woods-Saxon potential is

The LS coefficient is

So, this is the neutron single particle levels.

An adventure in Nuclear Physics!

August 15, 2019

Basic single particle energy, Woods-Saxon Leave a comment

The Woods-Saxon potential is

The LS coefficient is

So, this is the neutron single particle levels.

March 21, 2019

Basic mean field, Woods-Saxon Leave a comment

In the last post, I hope I explained how to find the Woods-Saxon energy levels for a given parameters. I just searched the best fit Woods-Saxon parameters to best fit the neutron single-particle energy for 209Pb. This is a double magic nucleus, and the there is no large fragmentation for the neutron excitation energy, thus, the outermost neutron in this nucleus can be considered as a good single-particle.

The best fit parameters are

The rms difference is just 78 keV!! For each level, the difference is not more than 50 keV!

So, how can we interpret this fitting result? It means, the energy levels can be well explained by WS mean field.

Recall that, the mean field actually included a lot things, it is the effective single particle potential that a nucleon is feeling. When we pushing the rms value to be minium. We are actually finding the best mean field. And the difference is due to the residual interaction.

But what contained in the mean field and residual interaction? The total Hamiltonian is

the mean field approach is to add an artificial potential

such that $latex H_{R} $ is minimum.

And by fitting energy levels using Woods-Saxon energy levels, we are basically doing the same thing! And the mean field is explicitly containing 2-body force, 3N-force and so on. So, can we say “because the Woods-Saxon can explain the energy level very well, the tensor force or other force is insignificant” ? the answer is no. Because the Woods-Saxon potential explicitly contains tensor force and other force.

March 21, 2019

Basic Runge-Kutta, Woods-Saxon Leave a comment

Hello everyone, I have been very busy lately, working on many things at the same time. Writing experiment proposal, testing a new detector array, flighting with the noise. Coding a program to read digitizer data for real time. Update the analysis code for coming experiments. Update/making some calculation codes, one of them is the calculation of Woods-Saxon energy levels.

Using RK4 method, a 2ndary 1-D differential potential can be solved easily (with some effort).

The Woods-Saxon potential with spin-orbital coupling is:

There are 6 parameters, 3 for the central term, 3 for the spin-orbital term.

And the equation is

Since the radial wavefunction has to be zero when it is far away, we can search for energy, such that the wavefunction is close to zero at a certain range. ( this is the only method I know how to practically implemented. I can imagine that, there may be some other method, transform the whole equation with a fixed r and directly solve the energy. Something like least-square method for linear system can be transformed into Matrix formalism and be directly solved. But I don’t know how)

The code is available at github (currently only for neutron).

I also have an excel to solve a single energy. Feel free to leave in comment and I will email it to you.

March 28, 2016

Basic rms, Woods-Saxon Leave a comment

The Woods-Saxon potential is

where is the half-maximum radius that is the reduced radius and is the nuclear mass number, and is the diffusiveness parameter.

The root-mean-squared (rms) radius is

where the “extract” is because of spherical coordinate.

The integration is a polynomial

In mathematica, the sum is notated by,

Thus, the rms radius for is

For , ,

February 16, 2011

Basic, Porperties 16O, 28Pb, Fermi, Woods-Saxon Leave a comment

this is a collective model of the nuclei density vs radius. it has another name Fermi-shape.

where is central density, or density at r = 0. is the radius of half density and is the diffuseness. when a is large, the “tail” of the shape will be longer.

the radius is measured in unit of fm .

can be vary from 1 fm to 7 or 8 fm. for 16O, it is about 2 fm. and for 208Pb, it is about 6fm .

is more or less the same for different nuclei.

the density of nuclei can be mass density or charge density, the Woods-Saxon also gives a good approximation.

for mass density, a . for comparison, water density is 1 kg per meter cube.

and charge density is about .