Two-state vibrational frequency memory

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.