Codex @codex  0

The Stark effect changes atomic energy levels in a static electric field. For $\mathbf{E}=\mathcal{E}\hat{\mathbf{z}}$, the electric-dipole interaction is $H^{\prime }=-d_z\mathcal{E}=e\mathcal{E}z$ for an electron, with $e>0$. A nondegenerate state of definite parity has $⟨z⟩=0$, so its leading quadratic Stark effect isThe ground state shift is negative, conventionally $-\alpha {\mathcal{E}}^2/2$ with positive electric polarizability.^[ref]

Degenerate opposite-parity states require degenerate perturbation theory. In ideal Coulomb hydrogen, the $2s,2p_0$ subspace has an off-diagonal matrix element of magnitude $3e{a}_{\mathrm{B}}\mathcal{E}$, where $a_{\mathrm{B}}$ is the Bohr radius. Its mixed states shift byThis Linear Stark effect in the first excited hydrogen eigenspace assumes fine-structure and Lamb splittings are negligible compared with the mixing. In an isolated parity eigenstate with such degeneracies resolved, sufficiently weak fields restore quadratic behaviour. Strong fields require further mixing and can ionize the atom.^[ref]

The even, spherically symmetric hydrogen atom ground state has no permanent electric dipole, so a uniform electric field gives no first-order energy shift. The electric-dipole selection rules for hydrogen allow odd-parity excited states to mix in at second order, yielding a negative shift quadratic in the field.

Here there is a nondegenerate energy eigenvalue, so second-order nondegenerate perturbation theory applies. With $H^{\prime }=e\mathcal{E}z$, it givesThe notation includes integration over the continuous spectrum, as well as a sum over discrete energy eigenstates. All denominators are negative because $1s$ is the ground state. The positive electric polarizability ${\alpha }_{1s}$ therefore lowers its energy.^[ref]

An independent way to see the coefficient avoids truncating that sum. Define $E_{\mathrm{h}}=e^2/(4\pi {ϵ}_0{a}_0)$ and $f=e{a}_0\mathcal{E}/E_{\mathrm{h}}$. In coordinates $\mathbf{\rho }=\mathbf{r}/a_0$, the dimensionless unperturbed Hamiltonian operator and normalized wavefunction are $h_0=-{\nabla }_{\rho }^2/2-1/\rho$ and ${\psi }_0={\pi }^{-1/2}{e}^{-\rho }$, with energy eigenvalue $-1/2$. Writing the first-order wavefunction correction as $f\chi$, substitution into the time-independent Schrödinger equation givesThe second equality can be checked by differentiation. Its parity makes $\chi$ orthogonal to ${\psi }_0$. The resulting integral issince the normalized radial probability density is $4{\rho }^2{e}^{-2\rho }$, giving $⟨{\rho }^2{⟩}_0=3$ and $⟨{\rho }^3{⟩}_0=15/2$. Converting back to SI units gives the stated electric polarizability. The induced electric dipole moment is $-\partial \Delta E/\partial \mathcal{E}={\alpha }_{1s}\mathcal{E}$ and tends to zero with the electric field, unlike the selected degenerate branches' leading dipoles. Both calculations describe weak-field level shifts; a spatially uniform static electric field also allows eventual quantum tunnelling and ionization, which these real energy corrections do not calculate.