This website uses the Flot javascript library to plot [r,rho(r)] data and model the charge density of various nuclei. Data for these models is taken from Atomic and Nuclear Data Tables, Volumes 14, 36 and 60, which are provided on the Downloads page.
This website also provides the source code that generated the distributions, and are provided for anyone to edit as they see fit.
ρ=ρo(1+α(r/a)2)exp(-(r/a)2)
With ρo calculated such that 4π ∫ ρ(r)r2dr = Ze
Same as HO but with α as a free parameter
ρ=ρo/(1+exp((r - c)/z))
With ρo calculated such that 4π ∫ ρ(r)r2dr = Ze
ρ=ρo(1+wr2/c2)/(1+exp((r-c/z))
With ρo calculated such that 4π ∫ ρ(r)r2dr = Ze
ρ=ρo/(1+exp((r2-c2)/a2))
With ρo calculated such that 4π ∫ ρ(r)r2dr = Ze
ρ=ρo(1+ wr2/c2)/(1+exp((r2-c2/z2))
With ρo calculated such that 4π ∫ ρ(r)r2dr = Ze
Σ R(n)*sin(x)/x, with x = n * π * r/RR
The Data File is set up as follows:
| Atom | Z | A | R(n) | R(n) | .... | RR |
Where R(n) is the list of the Fourier-Bessel coefficients, up to R[17].
ρ= Σ Ai{exp(-[(r-Ri)/γ]2)+exp(-[(r+Ri)/ γ]2)}, with Ai = ZeQi/[2 π2/3γ3(1+2Ri2/ γ2)]
The Data File is set up as follows:
| Z | A | Atom | rms | RP | Qi | Ri | Qi | Ri | ... |
Where Qi and Ri are the amplitude and position of the Gaussians.
The values of Qi indicate the fraction of the charge contained in the ith Gaussian normalized such that ΣiQi=1
RP is the rms radius of the Gaussians. RP
= γ
√ 3/2