PotentialMarinescu1994Cesium
Class Methods
|
Initialize the potential. |
Calculate the effective centrifugal potential V_l(x) in atomic units. |
|
Calculate the effective potential V_sqrt(x) from the sqrt transformation in atomic units. |
|
Calculate the classical turning point z_i of the state if it would be a hydrogen atom. |
|
Calculate the model potential by Marinescu et al. (1994) in atomic units. |
|
Calculate the Coulomb potential V_Col(x) in atomic units. |
|
Calculate the total effective potential V_eff(x) in atomic units. |
|
|
Calculate the classical inner turning point z_i for the given state. |
Class Attributes and Properties
Static dipole polarizability in atomic units (a.u.), used for the parametric model potential. |
|
Whether this potential is the default potential for the species. |
|
(a_1, a_2, a_3, a_4)} for the parametric model potential. |
|
r_c} to truncate the unphysical short-range contribution of the polarization potential. |
|
The short name of the atomic species. |
|
The tag for these potential parameters. |
- class rydstate.species.cesium.PotentialMarinescu1994Cesium(l_r)[source]
Initialize the potential.
- Parameters:
l_r (
int) – Orbital angular momentum of the Rydberg electron.args (object)
kwargs (object)
- Return type:
CachedT
- species: ClassVar[str] = 'Cs'
The short name of the atomic species.
- is_default: ClassVar[bool] = True
Whether this potential is the default potential for the species.
- alpha_c_marinescu_1994: ClassVar[float] = 15.644
Static dipole polarizability in atomic units (a.u.), used for the parametric model potential.
- r_c_dict_marinescu_1994: ClassVar = {0: 1.9204693, 1: 2.13383095, 2: 0.93007296, 3: 1.99969677}
r_c} to truncate the unphysical short-range contribution of the polarization potential.
- Type:
Cutoff radius {l
- model_potential_parameter_marinescu_1994: ClassVar = {0: (3.49546309, 1.475338, -9.72143084, 0.02629242), 1: (4.69366096, 1.71398344, -24.6562428, -0.09543125), 2: (4.32466196, 1.61365288, -6.7012885, -0.74095193), 3: (3.01048361, 1.40000001, -3.20036138, 0.00034538)}
(a_1, a_2, a_3, a_4)} for the parametric model potential.
- Type:
Parameters {l
- calc_effective_potential_centrifugal(x)
Calculate the effective centrifugal potential V_l(x) in atomic units.
The effective centrifugal potential is given as
\[V_{l_r}(x) = \frac{l_r(l_r+1)}{2x^2}\]where x = r / a_0 and l_r is the orbital angular momentum quantum number of the Rydberg electron.
- Parameters:
x (
TypeVar(XType, NDArray,float)) – The dimensionless radial coordinate x = r / a_0, for which to calculate the potential.- Returns:
The effective centrifugal potential V_{l_r}(x) in atomic units.
- Return type:
TypeVar(XType, NDArray,float)
- calc_effective_potential_sqrt(x)
Calculate the effective potential V_sqrt(x) from the sqrt transformation in atomic units.
The sqrt transformation potential arises from the transformation from the wavefunction u(x) to w(z), where x = r / a_0 and w(z) = z^{-1/2} u(x=z^2) = (r/a_0)^{-1/4} sqrt(a_0) r R(r). Due to the transformation, an additional term is added to the radial Schrödinger equation, which can be written as effective potential V_{sqrt}(x) and is given by
\[V_{sqrt}(x) = \frac{3}{32x^2}\]- Parameters:
x (
TypeVar(XType, NDArray,float)) – The dimensionless radial coordinate x = r / a_0, for which to calculate the potential.- Returns:
The sqrt transformation potential V_sqrt(x) in atomic units.
- Return type:
TypeVar(XType, NDArray,float)
- calc_hydrogen_turning_point_z(n)
Calculate the classical turning point z_i of the state if it would be a hydrogen atom.
The hydrogen turning point is defined as the point, where for the idealized hydrogen atom the potential equals the energy, i.e. V_Col(r_i) + V_l(r_i) = E. This is exactly the case at
\[r_i = n^2 - n \sqrt{n^2 - l_r(l_r + 1)}\]and z_i = sqrt{r_i / a_0}.
- Parameters:
n (
int) – Principal quantum number of the state.- Returns:
The inner hydrogen turning point z_i in the scaled dimensionless coordinate z_i = sqrt{r_i / a_0}.
- Return type:
float
- calc_model_potential(x)
Calculate the model potential by Marinescu et al. (1994) in atomic units.
The model potential, see
reference, is given by\[V_{mp,marinescu}(x) = - \frac{Z_{l}}{x} - \frac{\alpha_c}{2x^4} (1 - e^{-x^6/x_c**6})\]where Z_{l} is the effective nuclear charge, \(\alpha_c\) is the static core dipole polarizability, and x_c is the effective core size.
\[Z_{l} = 1 + (Z - 1) \exp(-a_1 x) - x (a_3 + a_4 x) \exp(-a_2 x)\]with the nuclear charge Z.
- Parameters:
x (
TypeVar(XType, NDArray,float)) – The dimensionless radial coordinate x = r / a_0, for which to calculate potential.- Returns:
The four parameter potential V_{mp,marinescu}(x) in atomic units.
- Return type:
TypeVar(XType, NDArray,float)
- calc_potential_coulomb(x)
Calculate the Coulomb potential V_Col(x) in atomic units.
The Coulomb potential is given as
\[V_{Col}(x) = -1 / x\]where x = r / a_0.
- Parameters:
x (
TypeVar(XType, NDArray,float)) – The dimensionless radial coordinate x = r / a_0, for which to calculate the potential.- Returns:
The Coulomb potential V_Col(x) in atomic units.
- Return type:
TypeVar(XType, NDArray,float)
- calc_total_effective_potential(x)
Calculate the total effective potential V_eff(x) in atomic units.
The total effective potential includes all physical and effective potentials:
\[V_{eff}(x) = V(x) + V_{l_r}(x) + V_{sqrt}(x)\]where V(x) is the physical potential (either Coulomb or a model potential), V_{l_r}(x) is the effective centrifugal potential, and V_{sqrt}(x) is the effective potential from the sqrt transformation.
Note that we on purpose do not include the spin-orbit potential for several reasons:
The fine structure corrections are important for the energies of the states. This includes a) spin-orbit coupling, b) Darwin term, and c) relativistic corrections to the kinetic energy. Since we (obviously) can not include the latter two in the potential, it is only consistent to not include the spin-orbit term either.
The model potentials are generated without the spin-orbit term, since their accuracy is not sufficient to resolve the fine structure corrections at small distances. (This can also be seen by running Numerov for low lying states with an energy changed by e.g. 1%, which will lead to almost no change in the wavefunction.)
- Parameters:
x (
TypeVar(XType, NDArray,float)) – The dimensionless radial coordinate x = r / a_0, for which to calculate potential.- Returns:
The total potential V_eff(x) in atomic units.
- Return type:
TypeVar(XType, NDArray,float)
- calc_turning_point_z(energy_au, dz=0.001)
Calculate the classical inner turning point z_i for the given state.
The classical turning point is defined as the point, where the total effective potential of the Rydberg model equals the energy, i.e. V_eff(r_i) = E.
Note: Because we use the total effective potential, even for l=0 the turning point is not at r=0. The advantage of this is, that this definition of the turning point should correspond to where w(z) should have its last change of sign in the second derivative.
- Parameters:
energy_au (
float) – The energy, for which to calculate the classical turning point in atomic units.dz (
float) – The precision of the turning point calculation.
- Returns:
The inner turning point z_i in the scaled dimensionless coordinate z_i = sqrt{r_i / a_0}.
- Return type:
float
- reference: ClassVar[str] = 'M. Marinescu et al., Phys. Rev. A 49, 982 (1994), https://journals.aps.org/pra/abstract/10.1103/PhysRevA.49.982'
- tag: ClassVar[str] = 'marinescu_1994'
The tag for these potential parameters.