Vibrational spectral diffusion

Vibrational spectral diffusion is the evolution of a molecular vibrational transition frequency as its surroundings or configuration fluctuate. It concerns frequency memory within a vibrational response, whereas vibrational relaxation concerns loss or redistribution of vibrational excitation. Their characteristic times need not coincide.

In two-dimensional infrared spectroscopy, an initial frequency label is followed by an independently varied 2D IR waiting time and a later frequency-sensitive readout. The vibrational frequency-frequency correlation function models the relation between those earlier and later frequencies, while a 2D IR center-line slope summarizes part of the measured spectral geometry under specified response assumptions.^[ref]

An original two-state vibrational frequency memory model has equally populated environments with spectroscopic frequencies ${\tilde{\nu }}_0-\Delta$ and ${\tilde{\nu }}_0+\Delta$. Suppose switching in either direction follows a continuous-time Markov chain with rate $k$. With a sign variable $s=\pm 1$, the probability of a flip during $dt$ is $kdt$; each flip reverses its frequency deviation. Hence the normalized sign memory satisfiesThe memory time is $1/(2k)$, not the mean residence time $1/k$ in either state. Returning to the initial state contributes to the measured memory, so these two time scales describe different events.

For stipulated transitions at 1000 and $1010\mathrm{c}{\mathrm{m}}^{-1}$ with $k=0.25\mathrm{p}{\mathrm{s}}^{-1}$, the mean is $1005\mathrm{c}{\mathrm{m}}^{-1}$, the variance is $25\mathrm{c}{\mathrm{m}}^{-2}$ and the memory time is 2 ps. At 2 ps the correlation is $25e^{-1}\approx 9.20\mathrm{c}{\mathrm{m}}^{-2}$. The probability of occupying the initial environment is $(1+e^{-1})/2\approx 0.684$, including trajectories with an even number of switches.

If an independent vibrational population in this illustrative model decays with lifetime 10 ps, its surviving fraction at 2 ps is $e^{-0.2}\approx 0.819$, while normalized frequency memory is only 0.368. This explicitly separates population survival from frequency decorrelation. These are model values, not measured solvent or vibrational constants.

The model uses spectroscopic wavenumber, not angular frequency. Converting to an angular-frequency correlation multiplies by $(2\pi c_0{)}^2$ when the speed of light $c_0$ uses compatible length units. Mixing these conventions changes the fluctuation amplitude without changing the switching time.

To test the interpretation, acquire spectra at several waiting times, characterize pulse timing and phase, and assess overlap, population decay and instrument response. A decaying spectral statistic can constrain frequency fluctuations, but the same observed decay may fit more than one microscopic environment model. The two-state calculation is a reusable mechanism example, not a universal explanation of spectral diffusion.

Two-state vibrational frequency memory models a probe that switches between two frequency environments while retaining a time-dependent relation to its earlier state. It distinguishes state occupancy, a vibrational frequency-state residence time and decay of the vibrational frequency-frequency correlation function.

Let a continuous-time Markov chain switch from state A to B at rate $a>0$ and from B to A at rate $b>0$. Its stationary probabilities areFor spectroscopic wavenumbers ${\nu }_A,{\nu }_B$, the stationary mean is $\overline{\nu }=p_A{\nu }_A+p_B{\nu }_B$ and the variance isThe centered state values form a decaying mode of the two-state transition generator, with rate $a+b$. Consequently,The vibrational frequency-state residence times are instead $1/a$ in A and $1/b$ in B. Correlation decay includes paths that leave and return, so it is not identical to survival without a switch.

An original example uses ${\nu }_A=1000$, ${\nu }_B=1012\mathrm{c}{\mathrm{m}}^{-1}$, $a=0.1$ and $b=0.3\mathrm{p}{\mathrm{s}}^{-1}$. Then $p_A=0.75$, $p_B=0.25$, $\overline{\nu }=1003\mathrm{c}{\mathrm{m}}^{-1}$, $C(0)=27\mathrm{c}{\mathrm{m}}^{-2}$ and ${\tau }_c=2.5\mathrm{ps}$. The residence times are $10$ and approximately $3.33\mathrm{ps}$.

Starting in A, the probability of being in A later isAt $2.5\mathrm{ps}$ it is about $0.842$, while the probability of never leaving A is only $e^{-aT}\simeq 0.779$. Their difference includes return trajectories. At long times, the stationary probability of matching the initial state approaches $p_A^2+p_B^2=0.625$, while normalized frequency covariance tends to zero. A nonzero same-state baseline therefore need not mean persistent covariance.

Primary two-dimensional infrared spectroscopy theory connects waiting-time-dependent spectral geometry with frequency memory under response approximations.^[ref] The rates and wavenumbers above are illustrative. A real assignment also needs spectral resolution, pulse characterization and a population-relaxation model; agreement with one exponential decay does not uniquely establish two microscopic states.

A vibrational frequency-state residence time is the duration a probe remains in a defined frequency environment before leaving it. It differs from the time over which the probe's frequency retains statistical memory, and from a vibrational relaxation lifetime describing excitation loss.

For a state with constant exit rate $k$, let $S(t)$ be the probability of no exit by time $t$. The small-interval exit model givesThe residence-time density is $f(t)=k{e}^{-kt}$, with mean $1/k$ and variance $1/k^2$. The exponential distribution is memoryless: $S(t+s)/S(t)=e^{-ks}$. These conclusions require a constant-rate model, rather than following from the name of a spectral state alone.

For a two-state continuous-time Markov chain with A-to-B rate $a$ and B-to-A rate $b$, the mean residence times are $1/a$ and $1/b$. The normalized two-state vibrational frequency memory instead decays as $e^{-(a+b)T}$. Returning after an exit contributes to occupancy and correlation but not to survival without an exit.

An original symmetric example uses $a=b=k=0.25\mathrm{p}{\mathrm{s}}^{-1}$. Mean residence time is $4\mathrm{ps}$, while frequency-memory time is $2\mathrm{ps}$. Starting in A, the probability of being in A after $T$ iswhereas no-exit probability is $e^{-kT}$. At $T=2\mathrm{ps}$ these are about $0.684$ and $0.607$.

This produces vibrational residence-time sampling bias if a same-state observation is treated as proof that no exit occurred. Using $-\ln P_{AA}(2)/2$ as an exit-rate estimate gives about $0.190\mathrm{p}{\mathrm{s}}^{-1}$, below the stipulated $0.25\mathrm{p}{\mathrm{s}}^{-1}$. Faster time sampling can resolve more excursions, but the estimator must still account for missed transitions and state-assignment uncertainty.

In two-dimensional infrared spectroscopy, the 2D IR waiting time separates initial frequency labeling from later response readout. Primary response theory relates such measurements to frequency memory under explicit approximations.^[ref] A fitted memory decay is therefore not automatically a directly measured single-state residence time. Nonexponential survival can require additional states, time-dependent rates or a different physical model.

Vibrational residence-time sampling bias occurs when finite-time observations miss exits and returns between frequency states, or when imperfect state assignment creates apparent transitions. The resulting apparent dwell behavior need not equal the underlying vibrational frequency-state residence time.

For an original symmetric two-state continuous-time Markov chain with exit rate $k$ in either state, the probability that two samples separated by $\Delta$ differ isIt counts an odd number of switches, not every switch. The expected observed-change count per observation time is thereforewhich approaches $k$ only as the interval becomes small. At long intervals it approaches $1/(2\Delta )$ even if the true process switches much faster.

With $k=0.25\mathrm{p}{\mathrm{s}}^{-1}$ and $\Delta =2\mathrm{ps}$, $p_{\mathrm{flip}}\simeq 0.3161$ and $k_{\mathrm{count}}\simeq 0.1580\mathrm{p}{\mathrm{s}}^{-1}$. Treating same-state occupancy as no-exit survival gives a different naive estimate, $-\ln (1-p_{\mathrm{flip}})/\Delta \simeq 0.190\mathrm{p}{\mathrm{s}}^{-1}$. Both differ from the stipulated true rate.

Under this specific symmetric model, inversion givesUsing the exact model probability recovers the chosen rate. An estimated probability near $1/2$ makes the inverse very sensitive to uncertainty, and the nonlinear inverse is not automatically an unbiased finite-sample estimator.

Misclassification can bias in the opposite direction. If a physically unchanged two-state label is independently reversed with probability $q$ at each observation, the apparent transition probability is $2q(1-q)$. At $q=0.1$ it is $0.18$ despite no physical transition. Faster sampling alone does not remove this labeling error.

A useful measurement comparison varies sampling interval and checks state assignment against suitable references. In two-dimensional infrared spectroscopy, a 2D IR waiting time separates frequency labeling and later response, but spectra are not automatically direct individual-state trajectories.^[ref] Relating a spectral memory statistic to residence times requires a physical observation model as well as a transition model.