CamiXon.jl

A package for image analysis of backscattered light

CamiXon.jl

Introduction

CamiXon is a package for the numerical solution of the radial Schrödinger equation allowing for screening.

In this package the solution is obtained for a single electron, acting as a spectator in the central field of the atomic nucleus screened by 'the other' electrons (the screening electrons).

The starting point is the 1D Schrödinger equation,

\[\tilde{χ}_{l}^{′′}+2[Z_{\mathrm{eff}}(ρ)/ρ-l(l+1)/2ρ^{2}+ε_{l}]\tilde{χ}_{l}=0,\]

where $\tilde{\chi}_{l}(\rho)=\rho\tilde{R}_{l}(\rho)$ is the reduced radial wavefunction and $ε_{l}=-κ _{l}^{2}$ is the corresponding binding energy in Hartree atomic units (a.u.). As compared to the hydrogenic case, the atomic number $Z$ has been replaced by $Z_{\mathrm{eff}}(ρ)$, the effective nuclear charge at radial distance $ρ$ from the atomic center (in a.u.). In other words, the energy of the electron in the Coulomb field of the nucleus is replaced by an effective, central field potential of the form

\[U_{\mathrm{CF}}(ρ)=-Z_{\mathrm{eff}}(ρ)/ρ=-Z/ρ+U_{\mathrm{scr}}(ρ),\]

consisting of the bare Coulomb contribution, $-Z/ρ$, and the screening field $U_{\mathrm{scr}}(ρ)$, subject to the boundary conditions $U_{\mathrm{scr}}(0)=Z$ and $\mathrm{lim}_{ρ→\infty}U_{\mathrm{scr}}(ρ)=Z_{c}/ρ$. Here $Z_{c}$ is the Rydberg charge; i.e. the effective nuclear charge for a spectator electron in the far field $(ρ→\infty)$. Within these assumptions we can optimize $U_{\mathrm{scr}}(ρ)$, while preserving the bare Coulomb field close to the nucleus as well as the Rydberg potential in the far field. The price we pay is that the radial Schrödinger equation has to be solved numerically by radial integration. Our strategy is to use both inward and outward integration and match the two branches by equating the two solutions for the wavefunction, $χ(ρ)$, and its derivative, $χ^′(ρ)$, at a point near the classical turning point of the radial motion of the electron. The basics of the solution can be found in the book Atomic Structure Theory by Walter R. Johnson.

Illustration: the hydrogen 3d orbital

Shown below are the reduced radial wavefunction $(χ)$ and its derivative $(χ^′)$ in the near field (left), the far field (right), and in the region near the classical turning point (center).

Codata

CamiXon.Codata — Type

Codata

Object to hold the natural constants from CODATA. It is best created with the function castCodata

The fields are:

.∆νCs: Cs hyperfine transition frequency (::Value)
.c: speed of light in vacuum (::Value)
.h: Planck constant (::Value)
.ħ: Planck constant - reduced (::Value)
.e: elementary charge (::Value)
.kB: Boltzmann constant (::Value)
.NA: Avogadro constant (::Value)
.Kcd: Luminous efficacy (::Value)
.me: electron rest mass (::Value)
.R∞: Rydberg constant (::Value)
.Ry: Rydberg frequency (::Value)
.Eh: Hartree a.u. (::Value)
.α: fine-structure constant (::Value)
.μ0: magnetic permitivity of vacuum (::Value)
.ε0: electric permitivity of vacuum (::Value)
.KJ: Josephson constant (::Value)
.RK: Von Klitzing constant (::Value)
.R: Molar gas constant (::Value)
.u: unified atomic mass unit (::Value)
.matE: unit conversion matrix (Matrix{Float64})

Example:

codata = castCodata(2018)
codata.μ0
  Value(1.2566370621250601e-6, "N A⁻²")

codata.μ0.val
  1.2566370621250601e-6