# 14.3: Gaussian-Stochastic Model for Spectral Diffusion

- Page ID
- 107302

We will begin with a *classical *description of how random fluctuations in frequency influence the absorption lineshape, by calculating the dipole correlation function for the resonant transition. This is a Gaussian stochastic model for fluctuations, meaning that we will describe the time-dependence of the transition energy as random fluctuations about an average value through a Gaussian distribution.

\[ \begin{align} \omega (t) &= \langle \omega \rangle + \delta \omega (t) \label{13.9} \\[4pt] \langle \delta \omega (t) \rangle &= 0 \label{13.10} \end{align}\]

The fluctuations in \(\omega\) allow the system to explore a Gaussian distribution of transitions frequencies characterized by a variance:

\[\Delta = \sqrt {\left\langle \omega^{2} \right\rangle - \langle \omega \rangle^{2}} = \sqrt {\left\langle \delta \omega^{2} \right\rangle} \label{13.11}\]

In many figures the width of the Gaussian distribution is labeled with the standard deviation (here \(\Delta\)). This is meant to symbolize that \(\Delta\) is the parameter that determines the width, and not that it is the line width. For Gaussian distributions, the full line width at half maximum amplitude (FWHM) is \(2.35 \Delta\).

The time scales for the frequency shifts will be described in terms of a frequency correlation function

\[C _ {\delta \omega s \omega} (t) = \langle \delta \omega (t) \delta \omega ( 0 ) \rangle \label{13.12}\]

Furthermore, we will describe the time scale of the random fluctuations through a correlation time \(\tau_c\).

The absorption lineshape is described with a dipole time-correlation function. Let’s treat the dipole moment as an internal variable to the system, whose value depends on that of \(\omega\). Qualitatively, it is possible to write an equation of motion for \(\mu\) by associating the dipole moment with the displacement of a bound particle (\(x\)) times its charge, and using our intuition regarding how the system behaves. For a unperturbed state, we expect that \(x\) will oscillate at a frequency \(\omega\), but with perturbations, it will vary through the distribution of available frequencies. One function that has this behavior is

\[x (t) = x _ {0} e^{- i \omega (t) t} \label{13.13}\]

If we differentiate this equation with respect to time and multiply by charge we have

\[\frac {\partial \mu} {\partial t} = - i \omega (t) \mu (t) \label{13.14}\]

Although it is a classical equation, note the similarity to the quantum Heisenberg equation for the dipole operator:

\[\partial \mu / \partial t = i H (t) \mu / \hbar + h . c .\]

The correspondence of \(\omega (t)\) with \(H (t) / \hbar\) offers some insight into how the quantum version of this problem will look.

The solution to Equation \ref{13.14} is

\[\mu (t) = \mu ( 0 ) \exp \left[ - i \int _ {0}^{t} d \tau \, \omega ( \tau ) \right] \label{13.15}\]

Substituting this expression and Equation \ref{13.9} into the dipole correlation function gives or

\[C _ {\mu \mu} (t) = | \mu |^{2} e^{- i \langle \omega \rangle t} F (t) \label{13.16}\]

where

\[F (t) = \left\langle \exp \left[ - i \int _ {0}^{t} d \tau \, \delta \omega ( \tau ) \right] \right\rangle \label{13.17}\]

The dephasing function (\(F(t)\)) is obtained by performing an equilibrium average of the exponential argument over fluctuating trajectories. For ergodic systems, this is equivalent to averaging long enough over a single trajectory.

The dephasing function is a bit complicated to work with as written. However, for the case of Gaussian statistics for the fluctuations, it is possible to simplify \(F(t)\) by expanding it as a **cumulant expansion of averages*** *(see appendix below for details).

\[F (t) = \exp \left[ - i \int _ {0}^{t} d \tau^{\prime} \, \left\langle \delta \omega \left( \tau^{\prime} \right) \right\rangle + \frac {i^{2}} {2 !} \int _ {0}^{t} d \tau^{\prime} \, \int _ {0}^{t} d \tau^{\prime \prime} \, \left\{\left\langle \delta \omega \left( \tau^{\prime} \right) \delta \omega \left( \tau^{\prime \prime} \right) \right\rangle - \left\langle \delta \omega \left( \tau^{\prime} \right) \right\rangle \left\langle \delta \omega \left( \tau^{\prime \prime} \right) \right\rangle \right\} \right] \label{13.18}\]

In this expression, the first term is zero, since \(\langle \delta \omega \rangle = 0\). Only the second term survives for a system with Gaussian statistics. Now recognizing that we have a stationary system, we have

\[F (t) = \exp \left[ - \frac {1} {2} \int _ {0}^{t} d \tau^{\prime} \, \int _ {0}^{t} d \tau^{\prime \prime} \, \left\langle \delta \omega \left( \tau^{\prime} - \tau^{\prime \prime} \right) \delta \omega ( 0 ) \right\rangle \right] \label{13.19}\]

We have rewritten the dephasing function n terms of a correlation function that describes the fluctuating energy gap. Note that this is a classical exception, so there is no time-ordering to the exponential. \(F(t)\) can be rewritten through a change of variables (\(\tau = \tau^{\prime} - \tau^{\prime \prime}\)):

\[F (t) = \exp \left[ - \int _ {0}^{t} d \tau ( t - \tau ) \langle \delta \omega ( \tau ) \delta \omega ( 0 ) \rangle \right] \label{13.20}\]

So the Gaussian stochastic model allows the influence of the frequency fluctuations on the lineshape to be described by \(C _ {\delta \omega \delta v} (t)\) a frequency correlation function that follows Gaussian statistics. Note, we are now dealing with two different correlation functions \(C _ {\delta \omega \delta \omega}\) and \(C _ {\mu \mu}\). The frequency correlation function encodes the dynamics that result from molecules interacting with the surroundings, whereas the dipole correlation function describes how the system interacts with a light field and thereby the absorption spectrum.

Now, we will calculate the lineshape assuming that \(C _ {\delta \omega \delta \omega}\) decays with a correlation time \(\tau_c\) and takes on an exponential form

\[C _ {\delta \omega \delta \omega} (t) = \Delta^{2} \exp \left[ - t / \tau _ {c} \right] \label{13.21}\]

Then Equation \ref{13.20} gives

\[F (t) = \exp \left[ - \Delta^{2} \tau _ {c}^{2} \left( \exp \left( - t / \tau _ {c} \right) + t / \tau _ {c} - 1 \right) \right] \label{13.22}\]

which is in the form we have seen earlier \(F (t) = \exp ( - g (t) )\)

\[g (t) = \Delta^{2} \tau _ {c}^{2} \left( \exp \left( - t / \tau _ {c} \right) + t / \tau _ {c} - 1 \right) \label{13.23}\]

to interpret this lineshape function, let’s look at its limiting forms.

Long correlation times (\( t \ll \tau_c\))

This corresponds to the inhomogeneous case where \(C _ {\delta \omega \delta \omega} (t) = \Delta^{2}\), a constant. For \(t \ll \tau _ {c}\), we can perform a short time expansion of exponential

\[e^{- t / \tau _ {c}} \approx 1 - \frac {t} {\tau _ {c}} + \frac {t^{2}} {2 \tau _ {c}^{2}} + \ldots \label{13.24}\]

and from Equation \ref{13.23} we obtain

\[g (t) = \Delta^{2} t^{2} / 2 \label{13.25}\]

At short times, the dipole correlation function will have a Gaussian decay with a rate given by \(\Delta^{2}\):

\[F (t) = \exp \left( - \Delta^{2} t^{2} / 2 \right)\]

This has the proper behavior for a classical correlation function, i.e., even in time

\[C _ {\mu \mu} (t) = C _ {\mu \mu} ( - t ).\]

In this limit, the absorption lineshape is:

\[\begin{align} \sigma ( \omega ) & = | \mu |^{2} \int _ {- \infty}^{+ \infty} d t \, e^{i \omega t} e^{- i ( \omega ) t - g (t)} \\ & = | \mu |^{2} \int _ {- \infty}^{+ \infty} d t \, e^{i ( \omega - ( \omega ) ) t} e^{- \Delta^{2} t^{2} / 2} \\ & = \sqrt {\frac {2 \pi} {\Delta^{2}}} | \mu |^{2} \exp \left( - \frac {( \omega - \langle \omega \rangle )^{2}} {2 \Delta^{2}} \right) \label{13.26} \end{align}\]

We obtain a Gaussian inhomogeneous lineshape centered at the mean frequency with a width dictated by the frequency distribution.

Short Correlation Times (\( t \gg \tau_c\))

This corresponds to the homogeneous limit in which you can approximate

\[C _ {\delta \omega \delta \omega} (t) = \Delta^{2} \delta (t).\]

For \(t \gg \tau _ {c}\) we set \(e^{- t / \tau _ {c}} \approx 0\), \(t / \tau _ {c} > > 1\), and Equation \ref{13.23} gives

\[g (t) = - \Delta^{2} \tau _ {c} t \label{13.27}\]

If we define the constant

\[\Delta^{2} \tau _ {c} \equiv \Gamma \label{13.28}\]

we see that the dephasing function has an exponential decay:

\[F (t) = \exp [ - \Gamma t ] \label{13.29}\]

The lineshape for short correlation times (or fast fluctuations) takes on a Lorentzian shape

\[ \begin{array} {c} {\sigma ( \omega ) = | \mu |^{2} \int _ {- \infty}^{+ \infty} d t \, e^{i ( \omega - \langle \omega ) t} e^{- \Gamma t}} \\ {\operatorname {Re} \sigma ( \omega ) = | \mu |^{2} \frac {\Gamma} {( \omega - \langle \omega \rangle )^{2} + \Gamma^{2}}} \end{array} \label{13.30}\]

This represents the homogeneous limit. Even with a broad distribution of accessible frequencies, if the system explores all of these frequencies on a time scale fast compared to the inverse of the distribution (\(\Delta \tau _ {\mathrm {c}} > 1\), then the resonance will be “motionally narrowed” into a Lorentzian line.

More generally, the envelope of the dipole correlation function will look Gaussian at short times and exponential at long times.

The correlation time is the separation between these regimes. The behavior for varying time scales of the dynamics (\(\tau_c\)) are best characterized with respect to the distribution of accessible frequencies (\(\Delta\)). So we can define a factor

\[\kappa = \Delta \cdot \tau _ {c} \label{13.31}\]

\(\kappa \ll 1\) is the *fast modulation limit *and \(\kappa \gg 1\) is the *slow modulation limit*. Let’s look at how \(C _ {\delta o \delta \omega}\), \(F (t)\), and \(\sigma _ {a b s} ( \omega )\) change as a function of \(\kappa\).

We see that for a fixed distribution of frequencies \(\Delta\) the effect of increasing the time scale of fluctuations through this distribution (decreasing \(\tau_c\)) is to gradually narrow the observed lineshape from a Gaussian distribution of static frequencies with width (FWHM) of \(2.35\Delta\) to a motionally narrowed Lorentzian lineshape with width (FWHM) of

\[\Delta^{2} \tau _ {c} / \pi = \Delta \cdot \kappa / \pi.\]

This is analogous to the motional narrowing effect first described in the case of temperature dependent NMR spectra of two exchanging species. Assume we have two resonances at \(\omega_A\) and \(\omega_B\) associated with two chemical species that are exchanging at a rate \(k_{AB}\)

\[\ce{A <=> B}\]

If the rate of exchange is slow relative to the frequency splitting, \(k _ {A B} < < \omega _ {A} - \omega B\) then we expect two resonances, each with a linewidth dictated by the molecular relaxation processes (\(T_2\)) and transfer rate of each species. On the other hand, when the rate of exchange between the two species becomes faster than the energy splitting, then the two resonances narrow together to form one resonance at the mean frequency.

## Appendix: The Cumulant Expansion

For a statistical description of the random variable \(x\), we wish to characterize the moments of \(x\): \(\langle x \rangle\), \(\langle x^2 \rangle\), .... Then the average of an exponential of \(x\) can be expressed as an expansion in moments

\[\underbrace{\left\langle e^{i k x} \right\rangle = \sum _ {n = 0}^{\infty} \frac {( i k )^n} {n !} \left\langle x^n \right\rangle}_{\text{expansion in moments}} \label{13.31A}\]

An alternate way of expressing this expansion is in terms of cumulants

\[\underbrace{\left\langle e^{i k x} \right\rangle = \exp \left( \sum _ {n = 1}^{\infty} \frac {( i k )^{n}} {n !} c _ {n} (x) \right)}_{\text{expansion in cumulants}} \label{13.32}\]

where the first few cumulants are:

\[ \begin{align} c _ {1} (x) &= \langle x \rangle \tag{mean} \label{13.33} \\[4pt] c _ {2} (x) &= \left\langle x^{2} \right\rangle - \langle x \rangle^{2} \label{13.34} \tag{variance} \\[4pt] c _ {3} (x) &= \left\langle x^{3} \right\rangle - 3 \langle x \rangle \left\langle x^{2} \right\rangle + 2 \langle x \rangle^{3} \tag{skewness} \label{13.35} \end{align}\]

An expansion in cumulants converges much more rapidly than an expansion in moments, particularly when you consider that \(x\) may be a time-dependent variable. Particularly useful is the observation that all cumulants with \(n > 2\) vanish for a system that obeys Gaussian statistics.

We obtain the cumulants above by expanding Equation \ref{13.31} and \ref{13.32}, and comparing terms in powers of \(x\). We start by postulating that, instead of expanding the exponential directly, we can instead expand the exponential argument in powers of an operator or variable \(H\)

\[F = \exp [ c ] = 1 + c + \frac {1} {2} c^{2} + \cdots \label{13.36}\]

\[c = c _ {1} H + \frac {1} {2} c _ {2} H^{2} + \cdots \label{13.37}\]

Inserting Equation \ref{13.37} into Equation \ref{13.36} and collecting terms in orders of \(H\) gives

\[\begin{aligned} F & = 1 + \left( c _ {1} H + \frac {1} {2} c _ {2} H^{2} + \cdots \right) + \frac {1} {2} \left( c _ {1} H + \frac {1} {2} c _ {2} H^{2} + \cdots \right)^{2} + \cdots \\ & = 1 + \left( c _ {1} \right) H + \frac {1} {2} \left( c _ {2} + c _ {1}^{2} \right) H^{2} + \cdots \end{aligned} \label{13.38}\]

Now comparing this with the expansion of the exponential operator (of \(H\))

\[\begin{align} F & = \exp [ f H ] \\ & = 1 + f _ {1} H + \frac {1} {2} f _ {2} H^{2} + \cdots \label{13.39} \end{align}\]

allows one to see that

\[\begin{array} {l} {c _ {1} = f _ {1}} \\ {c _ {2} = f _ {2} - f _ {1}^{2}} \end{array} \label{13.40}\]

The cumulant expansion can also be applied to time-correlations. Applying this to the time-ordered exponential operator we obtain:

\[\begin{align} F (t) & = \left\langle \exp _ {+} \left[ - i \int _ {0}^{t} d t \omega (t) \right] \right\rangle \\ & \approx \exp \left[ c _ {1} (t) + c _ {2} (t) \right] \label{13.42} \end{align} \]

\[\begin{aligned} c _ {1} & = - i \int _ {0}^{t} d \tau \langle \omega ( \tau ) \rangle \\ c _ {2} & = - \int _ {0}^{t} d \tau _ {2} \int _ {0}^{\tau _ {2}} d \tau _ {1} \left\{\left\langle \omega \left( \tau _ {2} \right) \omega \left( \tau _ {1} \right) \right\rangle - \left\langle \omega \left( \tau _ {2} \right) \right\rangle \left\langle \omega \left( \tau _ {1} \right) \right\rangle \right\} \\ & = - \int _ {0}^{t} d \tau _ {2} \int _ {0}^{\tau _ {2}} d \tau _ {1} \left\langle \delta \omega \left( \tau _ {2} \right) \delta \omega \left( \tau _ {1} \right) \right\rangle \end{aligned} \label{13.43}\]

For Gaussian statistics, all higher cumulants vanish.

## Readings

- Kubo, R., A Stochastic Theory of Line-Shape and Relaxation. In Fluctuation, Relaxation and Resonance in Magnetic Systems, Ter Haar, D., Ed. Oliver and Boyd: Edinburgh, 1962; pp 23- 68.
- McHale, J. L., Molecular Spectroscopy. 1st ed.; Prentice Hall: Upper Saddle River, NJ, 1999.
- Mukamel, S., Principles of Nonlinear Optical Spectroscopy. Oxford University Press: New York, 1995.
- Schatz, G. C.; Ratner, M. A., Quantum Mechanics in Chemistry. Dover Publications: Mineola, NY, 2002; Sections 7.4 and 7.5.
- W. Anderson, P., A Mathematical Model for the Narrowing of Spectral Lines by Exchange or Motion. Journal of the Physical Society of Japan 1954, 9, 316-339.
- Wang, C. H., Spectroscopy of Condensed Media: Dynamics of Molecular Interactions. Academic Press: Orlando, 1985.