This page was generated from the Jupyter notebook dipole_matrix_elements.ipynb.

Calculating dipole matrix elements

from ryd_numerov.rydberg import RydbergState

state_i = RydbergState("Rb", 60, 2, j=3 / 2, m=1 / 2)
state_f = RydbergState("Rb", 60, 3, j=5 / 2, m=1 / 2)

radial = state_i.calc_radial_matrix_element(state_f, k_radial=1)
angular = state_i.calc_angular_matrix_element(state_f, "ELECTRIC", k_angular=1, q=0)
dipole = state_i.calc_matrix_element(state_f, "ELECTRIC", k_radial=1, k_angular=1, q=0)

print(f"Numerov radial matrix element: {radial}")
print(f"Numerov angular matrix element: {angular}")
print(f"Numerov dipole matrix element: {dipole}")
assert dipole.magnitude == (radial * angular).magnitude

Numerov radial matrix element: 137.7830757530064 bohr
Numerov angular matrix element: 0.4898979485566356
Numerov dipole matrix element: 67.49964615722135 bohr * elementary_charge