Small collection of functions to calculate typical beta-decay observables and other properties.
Log(f) calculation adapted from Juhani Kantele's "Handbook of Nuclear Spectrometry", Academic Press Limited, London.
Neutron penetrability calculation adapted from "Theoretical Nuclear Physics" by Blatt and Weisskopf.
ecoulomb(Z,N)
calculate the coulomb correction for a shell model calculation of isotope (Z,N)
calculateqbetashellmodel(groundstateparent,groundstatechild,zparent,nparent,zchild,nchild)
calculate the beta-decay Q value using shell model calculation of the ground state energies
groundstateparent:absolute energy of the parent ground state in MeV
groundstatechild:absolute energy of the child ground state in MeV
zparent: Z of the parent
nparent: N of the parent
zchild: Z of the child
nparent: N of the child
childActivity(x,A,λ)
A: initial activity
λ: decay probability (ln2/T12)
grandChildActivity(x,A,λ,μ)
A: initial daughter activity
λ: child decay probability
μ: grandchild decay probability
chainActivity(x,A,λ,n)
activity of the n'th member of a decay chain.
A:initial parent nuclei. If initial activity is desired as input use activity/λ[1]
λ:vector containing all decay probabilities in the chain
n:decay curve of n'th member
logf(z,Qᵦ,Eₓ)
Calculate the log10 of the Fermi function for allowed beta decay
z: atomic number of the parent
Qᵦ: β decay Q value in MeV
Eₓ: daughter state energy relative to the ground state energy in MeV
calculateT12(z,Qᵦ,Eₓ,BGT)
calculate halflife of the beta decay of an isotope given feedings to excited states
Qᵦ: β decay Q value in MeV
Eₓ: vector of daughter states relative to the ground state energy in MeV
BGT: vector of BGT values
z: atomic number of the parent
calculateIb(z,Qᵦ,Eₓ,BGT)
calculate branching ratios of the beta decay of an isotope given feedings to excited states
Qᵦ: β decay Q value in MeV
Eₓ: vector of daughter states relative to the ground state energy in MeV
BGT: vector of BGT values
z: atomic number of the parent
logftfromib(z,t₁₂,Qᵦ,Eₓ,Iᵦ)
calculate logft of a given transition to an excitated state
z: atomic number of the parent
t₁₂: decay halflife
Qᵦ: β decay Q value in MeV
Eₓ: energy of daughter state relative to the ground state energy in MeV
Iᵦ: partial branching values
logftfrombgt(bgt)
calculate the logft for a given BGT (not quenched)
nPenetrability(x,mass::Vector,Lorb)
calculates the neutron penetrability p(x,Lorb).
x is the excitation energy above Sₙ, Lorb is the neutron angular momentum
mass[1] is the recoil, mass[2] is the neutron mass.