Diatomic Molecules and Spectroscopy

Roberto Di Remigio, Luca Frediani

Spectroscopy probes the electronic structure of atoms and molecules by measuring their interaction with light. Different portions of the light spectrum can be used to explore different aspects of the electronic structure of a molecule. In this notebook, we will look closely at rotatonial and vibrational spectroscopy of diatomic molecules. Some properties: http://hyperphysics.phy-astr.gsu.edu/hbase/Tables/diatomic.html

Rotational Spectroscopy

Rotational spectroscopy uses microwave radiation (300 MHz to 300 GHz) to investigate the rotational levels of molecule, i.e. the distribution of the nuclear masses in the molecule. In a rotational spectroscopy experiment, transitions between the closely space rotational levels are observed and their energy differences can provide information on the molecular geometry and the masses of the nuclei. For a diatomic molecule, the energy of the transition is: \begin{equation} F(J) = BJ(J+1) \end{equation} with intensity \begin{equation} |\mu_{J+1, J}|^2 = \frac{J+1}{2J+1}\mu_0^2 \end{equation} $B$ and $\mu_0$ are the rotational constant and the dipole moment, respectively: \begin{equation} B = \frac{\hbar}{4\pi c I}; \quad I = \mu R_\mathrm{e}^2 \end{equation} Thus a rotational transition will be deteceted only if the molecule has a nonzero dipole moment in the ground state. Moreover, the selection rule: \begin{equation} \Delta J = \pm 1 \end{equation} applies.

Vibrational Spectroscopy

Vibrational spectoscopy uses infrared radiation (430 THz to 300 GHz) to investigate the relative motions of nuclei in molecules, i.e. molecular vibrations. In the harmonic model of diatomic molecules, the wavenumber for the vibrational levels are obtained as: \begin{equation} G(v) = (v+\frac{1}{2})\tilde{\nu}, \quad \tilde{\nu} = \frac{1}{2\pi c}\sqrt{\frac{k}{\mu}} \end{equation} The intensity of a vibrational transition can be calculated as: \begin{equation} I_{vw} = \left\langle \psi_v|\mu|\psi_w\right\rangle \end{equation} The first-order term in the Taylor expansion with respect to the nuclear displacements determines the selection rules for vibrational spectroscopy in the harmonic approximation. Vibrational motion is never exactly harmonic. In the harmonic model we cannot, in fact, observe bond dissociation which is clearly in contrast with experiment.

Computing the Rovibrational Spectrum

All vibrational levels have a stack of rotational levels associated. Thus, when measuring vibrational spectra in gas phase one can observe transitions between these stacks of rotational levels associated with different vibrational levels. The energy of such transitions is in general expressed as: \begin{equation} S(v, J) = G(v) + F_v(J) \end{equation} where $G(v)$ can contain anharmonic terms and $F_v(J)$ centrifugal distortion terms, i.e. depending on the vibrational quantum number $v$.

In the simplest case, these nonlinear terms can be neglected, obtaining the simplified formula: \begin{equation} S(v, J) = (v+\frac{1}{2})\tilde{\nu} + BJ(J+1) \end{equation} where the $\tilde{\nu}$ and $B$ are the vibrational and rotational constants, respectively. Due to the small energy separation between rotational levels, the rovibrational spectrum will be temperature dependent. Each state, with energy $S(v, J)$, will be represented with a weight according to the Boltzmann distribution: \begin{equation} w(v, J, T) = \frac{\mathrm{deg}(J)\exp(-\frac{S(v, J)}{k_\mathrm{B}T})}{\sum_{v,J}\mathrm{deg}(J)\exp(-\frac{S(v, J)}{k_\mathrm{B}T})} \end{equation} $\mathrm{deg}(J)$ is the degeneration of the rotational state with quantum number $J$: states with higher degeneration will be more represented in the spectrum. The rovibrational spectrum captures the rotational transitions that occur when the vibrational transition $v+1\leftarrow v$ occurs. It consists of three branches:

  • The P branch with all transitions where $\Delta J = -1$.
  • The Q branch with the $\Delta J = 0$ transition. This can only occur with open-shell molecules.
  • The R branch with all transitions where $\Delta J = +1$.

Tasks

  1. Derive a general formula for the energy difference (in wavenumbers) of transitions in the P, Q and R branches of the spectrum
  2. Set up general script(s) to calculate and plot the rovibrational spectrum of a diatomic molecule. More in detail, the script(s) should:
    • Calculate the vibrational and rotational constants, given the masses, bond length and bond stifness constant.
    • Calculate the energies of the rotational levels for the P and R branches. We assume the Q branch to be forbidden.
    • Calculate the relative intensities of each transition in the P and R branches in terms of the Boltzmann distribution weight.
    • Plot the spectrum as wavenumber vs relative intensity.
  3. Compute the spectrum for the HCl molecule. An experimentally measured spectrum is reported below.

In [ ]: