Quantum Mechanical Calculation of Transition Probabilities and Stark Shifts
In the realm of atomic physics and quantum dynamics, understanding how atoms interact with electromagnetic fields is fundamental to mastering light-matter interaction. Two of the most critical parameters in this study are transition probabilities and Stark shifts. While the former quantifies the likelihood of an atom jumping between energy levels, the latter describes how an external electric field perturbs and splits those very levels.
Accurate calculations of these quantities are not merely theoretical exercises; they are indispensable for practical applications in high-resolution spectroscopy, laser physics, plasma diagnostics, and the development of quantum information processing technologies. This article explores the quantum mechanical frameworks used to derive these values, moving from fundamental perturbation theory to numerical computational methods.
Fermi's Golden Rule and Transition Probabilities
The probability of an atomic transition from an initial state $|i\rangle$ to a final state $|f\rangle$ is a measure of the coupling strength between the atom and an external field. In most spectroscopic contexts, this interaction is treated using time-dependent perturbation theory.
The Perturbation Hamiltonian
When an atom is subjected to an oscillating electromagnetic field, the interaction is typically modeled within the Electric Dipole Approximation. The perturbation Hamiltonian, $\hat{H}'$, is expressed as:
$$\hat{H}' = -\hat{\mu} \cdot \mathbf{E}(t)$$
where $\hat{\mu}$ is the electric dipole moment operator and $\mathbf{E}(t)$ represents the time-varying electric field.
Deriving the Transition Rate
For a transition driven by a monochromatic field, the rate of transition $W_{i \to f}$ (the probability per unit time) is governed by Fermi's Golden Rule:
$$ W_{i \to f} = \frac{2\pi}{\hbar} |\langle f | \hat{H}' | i \rangle|^2 \rho(E_f) $$
Here, $\langle f | \hat{H}' | i \rangle$ is the transition matrix element, which encapsulates the spatial overlap and symmetry between the two states, and $\rho(E_f)$ is the density of final states.
In practical atomic modeling:
- For Hydrogen-like systems: The wavefunctions are analytically known, allowing for direct calculation of the matrix elements.
- For Multi-electron atoms: The complexity increases significantly. One must account for electron-electron correlations and spin-orbit coupling. For instance, calculating the transition probability for the Sodium D-lines requires precise consideration of the fine structure to correctly predict the splitting and intensities.
Selection Rules and Line Shapes
The transition matrix element $\langle f | \hat{\mu} | i \rangle$ is subject to strict selection rules derived from the symmetry of the atomic states. For electric dipole transitions, these typically include:
- $\Delta l = \pm 1$ (change in orbital angular momentum)
- $\Delta m = 0, \pm 1$ (change in magnetic quantum number)
Furthermore, a transition is never a discrete delta function in reality. The Line Shape Function describes the distribution of the transition intensity. In the case of spontaneous emission, the line shape is typically Lorentzian, with a Full Width at Half Maximum (FWHM) denoted by $\Gamma$, which is inversely proportional to the excited state lifetime $\tau$:
$$\Gamma = \frac{\hbar}{\tau}$$
Integrating this line shape allows physicists to determine the total absorption or emission cross-section and the Einstein $A$ coefficients, which characterize the spontaneous emission rates.
The Quantum Mechanical Treatment of the Stark Effect
The Stark effect refers to the shifting and splitting of atomic energy levels in the presence of an external static electric field. This phenomenon is a direct consequence of the field breaking the spatial symmetry of the atomic potential.
Perturbation Theory in Non-Degenerate Systems
For many atoms, such as alkali metals, the energy levels are non-degenerate due to the quantum defect (the deviation from the pure hydrogenic model caused by the presence of a core). In these cases, the first-order Stark shift $\Delta E^{(1)}$ is often zero because the expectation value of the electric dipole operator vanishes due to parity conservation: $\langle n | \hat{\mu} | n \rangle = 0$.
Consequently, we must look to second-order perturbation theory to describe the energy shift:
$$ \Delta E^{(2)} = \sum_{n' \neq n} \frac{|\langle n' | \hat{\mu} \cdot \mathbf{E} | n \rangle|^2}{E_n - E_{n'}} $$
This result demonstrates that the shift is proportional to the square of the electric field strength ($E^2$), a phenomenon known as the quadratic Stark effect. The magnitude of the shift depends on the coupling between the state $|n\rangle$ and all other accessible states $|n'\rangle$, weighted by their energy proximity.
The Linear Stark Effect in Hydrogen
Hydrogen atoms present a unique case due to their high degree of degeneracy. In hydrogen, states with the same principal quantum number $n$ but different orbital angular momentum $l$ can have the same energy. This degeneracy allows for a non-zero first-order correction, leading to the linear Stark effect, where the energy shift is directly proportional to the field strength $E$.
The splitting is characterized by the formula:
$$\Delta E \propto n(k^2 - l^2/2)E$$
where $k$ represents the parabolic quantum numbers.
The Strong-Field Regime
When the external electric field is sufficiently strong, the perturbation becomes comparable to the internal atomic Coulomb field, and standard perturbation theory fails. In this regime, one must employ degenerate perturbation theory or, more commonly, numerical diagonalization. The new eigenstates, known as Stark states, become linear combinations of the original basis states. While they may lose their definite $l$ character, they retain certain symmetries, such as the projection of total angular momentum $m_j$ along the field axis.
Numerical Implementation and Practical Application
In modern computational physics, analytical solutions are often supplemented or replaced by numerical methods to handle complex multi-electron systems.
Matrix Diagonalization Approach
To calculate the Stark splitting of a specific manifold (e.g., the $n=2$ level of Hydrogen), the problem is transformed into an eigenvalue problem. The steps are as follows:
- Basis Construction: Define a basis set of degenerate or near-degenerate states (e.g., $|2,0,0\rangle, |2,1,0\rangle, |2,1,1\rangle, |2,1,-1\rangle$).
- Hamiltonian Matrix Assembly: Construct a matrix where the diagonal elements represent the unperturbed energies and the off-diagonal elements are the matrix elements of the interaction $\langle \psi_a | \hat{\mu} \cdot \mathbf{E} | \psi_b \rangle$.
- Eigenvalue Solution: Solve the secular equation $\det(H - E I) = 0$.
For the $n=2$ level of Hydrogen in a field $E_z$, the $4 \times 4$ matrix will yield three distinct energy levels: one shifted downwards, one shifted upwards, and one remaining at the original energy (due to the $m=0$ symmetry). This numerical approach provides a direct bridge between the abstract Schrödinger equation and observable spectral splitting.
Conclusion
The calculation of transition probabilities and Stark shifts represents a cornerstone of quantum spectroscopy. While Fermi's Golden Rule provides the framework for understanding how atoms absorb and emit radiation, perturbation theory allows us to predict how external environments reshape atomic structures. As we move toward more complex systems—such as highly charged ions or Rydberg atoms—the synergy between rigorous theoretical derivation and robust numerical computation remains the primary driver of discovery in atomic and molecular physics.