Class nucmass_dz_fit_33 (o2scl)¶
-
class
o2scl
::
nucmass_dz_fit_33
: public o2scl::nucmass_fit_base¶ The 33-parameter Duflo-Zuker mass formula.
This class is designed to provide essentially identical results to the original Duflo-Zuker code (see Duflo95) at
http://amdc.in2p3.fr/theory/dz31.f
(Two parameters were added by Duflo and Zuker after the fact to the original 31-parameter code, and still referred to as
dz31.f
.)The default values of nucmass::m_neut and nucmass::m_prot are adjusted to make sure that the mass excesses match the values given in the original.
- Todo:
Document each field.
Some explanations about the individual terms come from [MendozaTemis10] and the work by G. Bertsch at http://www.phys.washington.edu/users/bertsch/duflo2.ps
a[0]
: “Full master term”. Density sqaured divided by cube root of A. This is the master term which includes the bulk energy from the liquid droplet model and the harmonic oscillator effectsa[2]
: “Full spin-orbit term +”a[4]
: “Full spin-orbit term -“a[6]
: “Full cross term”a[8]
: “Partial master term”a[10]
: “Partial spin-orbit term +”a[12]
: “Partial spin-orbit term -“a[14]
: “S3”, polarizability of the valence spin-orbit shella[16]
: “SQ”, “QQM”, a neutron-proton interactiona[18], a[19]
: “D3” balance of monopole effectsa[20], a[21], a[24], a[25]
: “QQ+”/”QQ-“, Quadrupole terms, corresponding to filling equidistant Nilsson orbitsa[22], a[23]
: “D0”, Loss of monopole and gain of quadrupole energy for intrudersa[26]
: “TT”, Symmetry energya[28]
: “SS”a[30]
: “C”, Coulomb energya[31]
: “P0”, First pairing energy terma[32]
: “P1”, Second pairing energy term
For odd parameters up to
a[29]
, the odd parameter is just the preceeding even term divided by the cube root of A.Note that the original code states that,
.“for i even
a(i,program) =a(i-1,paper)*a(i,paper)”
Workspace vectors used internally
These are allocated in the class constructor.
-
tensor3
op
¶
-
tensor3
onp
¶
-
tensor3
ot
¶
-
ubvector_int
jup
¶
-
ubvector_int
jud
¶
-
ubmatrix_int
noc
¶
-
nucmass_dz_fit_33
()¶
-
~nucmass_dz_fit_33
()¶
-
const char *
type
()¶ Return the type,
"nucmass_dz_fit_33"
.
-
bool
is_included
(int Z, int N)¶ Return false if the mass formula does not include specified nucleus.
-
double
binding_energy
(int Z, int N)¶ Return the binding energy in MeV.
This function reproduces the original function called
EMASSDZ()
, except that, to be consistent with the other nuclear mass classes, it returns the binding energy with the opposite sign from the original.
-
double
binding_energy_d
(double Z, double N)¶ Return the binding energy in MeV.
-
double
mass_excess
(int Z, int N)¶ Given
Z
andN
, return the mass excess in MeV.
-
double
mass_excess_d
(double Z, double N)¶ Given
Z
andN
, return the mass excess in MeV.