---
title: Quasi-Classical Doorway-Window Approximation
url: https://www.emergentmind.com/topics/quasi-classical-doorway-window-approximation
type: topic
---

# Quasi-Classical Doorway-Window Approximation

Searching arXiv for the cited topic and papers to ground the article in current literature.
The quasi-classical doorway-window approximation is a factorized representation of time-resolved nonlinear spectroscopic signals in which the preparation of a nonstationary molecular ensemble by one set of pulses is encoded in a “doorway” function and its later interrogation by another set of pulses is encoded in a “window” function. In the formulation summarized in "WaveMixings.jl: a Julia package for performing on-the-fly time-resolved nonlinear electronic spectra from quasi-classical trajectories" [2509.03917] and "Ab Initio Simulation of Femtosecond Time-Resolved Multi-Pulse Spectroscopies applied to the Heptazine$\cdots$H$_2$O Complex" [2603.22671], the exact quantum response is reduced to averages over quasi-classical trajectories by combining a Wigner representation of the initial nuclear density with classical or surface-hopping propagation of nuclei. The resulting framework is designed for pump-probe, dispersed transient absorption, two-dimensional, 2D-FLEX, strong-field, pump-push-probe, and pump-induced two-dimensional spectroscopies, and is explicitly intended for on-the-fly electronic-structure simulations [2509.03917].

## 1. Response-function formulation and the doorway-window factorization

In the third-order response-function formalism, the heterodyne-detected four-wave mixing signal can be written, under the rotating-wave and impulsive approximations and assuming well-separated non-overlapping pulses, as a sum over ground-state bleach (GSB), stimulated emission (SE), and excited-state absorption (ESA) pathways [2509.03917]. For pump-probe spectroscopy this takes the form
$$
S^{(3)}(T_1,T_2,T_3)\propto \mathrm{Re}\sum_{k=\mathrm{GSB,SE,ESA}} a_k\,\mathrm{Tr}\!\left[W_k(T_3)\,e^{iLT_2}\,D(T_1)\right],
$$
where \(T_1\) is the coherence time between the first two pump pulses, \(T_2\equiv T\) is the population delay, \(T_3\equiv T'\) is the detection-time coherence interval, \(L\rho\equiv[H,\rho]/\hbar\) is the Liouvillian, \(D(T_1)\) is the doorway operator, and \(W_k(T_3)\) is the window operator for pathway \(k\). The coefficients are \(a_k=+1\) for GSB and SE and \(a_{\mathrm{ESA}}=-1\) for ESA [2509.03917].

A conceptually equivalent formulation appears in the multi-pulse generalization, where the detected signal is reorganized so that the preparation stage and the detection stage enter as separate phase-space factors [2603.22671]. In that representation, the pump-probe signal is written schematically as
$$
S(T)\approx \sum_{k=\mathrm{GSB,SE,ESA}} \int dX_0\, D(X_0)\,W_k[X(T;X_0)],
$$
with nuclear phase-space point \(X=(R,P)\), doorway \(D(X_0)\), and window \(W_k[X(T)]\) [2603.22671]. This factorized structure is the defining feature of the approximation.

The terminology is operational. The doorway encodes how the initial pulse sequence prepares an excited-state population or coherence, while the window encodes how the later pulse sequence detects the system through bleach, emission, or absorption channels [2603.22671]. In this sense, the approximation separates state preparation from readout without discarding the dependence of either object on nuclear geometry, transition dipoles, or pulse spectra.

## 2. Quasi-classical reduction from quantum response to trajectory averages

The quasi-classical approximation is obtained by reorganizing the trace so that nuclei evolve classically on adiabatic surfaces while electronic coherence contributions are evaluated at fixed nuclear geometries [2509.03917]. After a Wigner representation of the ground-state nuclear density and Monte Carlo sampling of initial conditions \(\{R_g(0),P_g(0)\}\), the pump-probe signal becomes
$$
S_{\mathrm{int}}(T,\omega_{\mathrm{pr}})
=
\mathrm{Re}\Big\langle
\sum_{k=\mathrm{GSB,SE,ESA}} a_k\,
W_k[\omega_{\mathrm{pr}};R_e(T),P_e(T)]\,
D[\omega_{\mathrm{pu}};R_g(0),P_g(0)]
\Big\rangle_{\mathrm{MC}},
$$
with explicit doorway and window functions built from adiabatic energy gaps, transition-dipole magnitudes, spectral envelopes, and the ground-state Wigner distribution [2509.03917].

For the integral pump-probe signal, the doorway is
$$
D(\omega_{\mathrm{pu}};R_g,P_g)=
|\mu_{ge}(R_g)|^2\,
E_{\mathrm{pu}}[\omega_{\mathrm{pu}}-\Delta E_{ge}(R_g)]\,
\mathrm{Wig}_g(R_g,P_g),
$$
and the three windows are
$$
W_{\mathrm{GSB}}(\omega_{\mathrm{pr}};R_g(T),P_g(T))
=
|\mu_{ge}(R_g(T))|^2
E_{\mathrm{pr}}[\omega_{\mathrm{pr}}-\Delta E_{ge}(R_g(T))],
$$
$$
W_{\mathrm{SE}}(\omega_{\mathrm{pr}};R_e(T),P_e(T))
=
|\mu_{ge}(R_e(T))|^2
E_{\mathrm{pr}}[\omega_{\mathrm{pr}}-\Delta E_{ge}(R_e(T))],
$$
$$
W_{\mathrm{ESA}}(\omega_{\mathrm{pr}};R_e(T),P_e(T))
=
\sum_{f\in\{II\}}
|\mu_{ef}(R_e(T))|^2
E_{\mathrm{pr}}[\omega_{\mathrm{pr}}-\Delta E_{fe}(R_e(T))].
$$
Here \(\mu_{ab}(R)\) are transition-dipole magnitudes, \(\Delta E_{ab}(R)=E_b(R)-E_a(R)\) are adiabatic energy gaps, \(E_{\mathrm{pu}}(\omega)\) and \(E_{\mathrm{pr}}(\omega)\) are pump and probe spectral envelopes, and \(\mathrm{Wig}_g\) is the ground-state Wigner distribution [2509.03917].

The same reduction is described in the multi-pulse treatment as the replacement of the exact quantum trace by a Wigner phase-space integral,
$$
\mathrm{Tr}[\cdot]\rightarrow \int dX\,\mathrm{Wigner}\{\cdot\}(X),
$$
together with classical phase-space evolution \(X(t)=\Phi_t(X_0)\) [2603.22671]. This places the doorway-window approximation at the interface of nonlinear response theory, semiclassical initial-state sampling, and mixed quantum-classical dynamics.

## 3. Core assumptions and their domain of validity

The derivation in the quasi-classical doorway-window framework invokes four explicit approximations [2509.03917]. First, pump and probe are assumed non-overlapping and separated by \(T\gg\) pulse durations, which allows mixed-time integrals to factorize into separate doorway and window operators. Second, nuclear motion during each ultrashort pulse is neglected, so nuclear coordinates are frozen at the time of each light-matter interaction. Third, electronic coherences during the short coherence intervals \(T_1\) and \(T_3\) evolve at fixed \(R\), contributing phase factors \(\exp[\pm i\,\Delta E(R)\,t/\hbar]\) but no nuclear dynamics. Fourth, the population evolution during \(T_2\) is replaced by classical or surface-hopping trajectories, so the quantum trace becomes a Monte Carlo average over sampled initial conditions and stochastic hops [2509.03917].

The range of validity is stated directly. Assumptions (i)–(iii) hold when pulses are \(\lesssim 10\) fs and \(T_1,T_3\) are controlled by pulse durations in pump-probe conditions [2509.03917]. Assumption (iv) inherits the accuracy of the chosen mixed quantum-classical propagation scheme, with Tully’s surface hopping given as an example [2509.03917].

The multi-pulse treatment emphasizes analogous constraints. It states that the framework is based on the usual “short, non-overlapping pulses” assumption and identifies limitations that follow from the approximation: it neglects nuclear and electronic coherences beyond classical transport, cannot describe interference between pathways, requires non-overlapping short pulses, and produces signal artifacts at overlap times \((<20\ \mathrm{fs})\) that are not quantitative [2603.22671]. Accuracy also inherits that of the classical-trajectory and electronic-structure level; the example given is that ADC(2) may mis-locate high states [2603.22671].

These statements delimit a common misconception. The method is not a full quantum treatment of nonlinear response, even though it starts from quantum response functions. It is instead a quasi-classical reduction that preserves explicit pulse, energy-gap, and transition-dipole dependence while replacing exact nuclear quantum propagation by classical or mixed quantum-classical dynamics [2509.03917].

## 4. Signal classes and generalizations beyond pump-probe

Within the doorway-window construction, transient absorption pump-probe is the basic third-order case. The multi-pulse summary gives the pump-probe signal as
$$
S_{\mathrm{PP}}(T;\omega_{\mathrm{pr}})
=
\sum_{k=\mathrm{GSB,SE,ESA}}
\int dX_0\,
D(X_0;\omega_{\mathrm{pu}})
\,W_k[\Phi_T(X_0);\omega_{\mathrm{pr}}],
$$
with
$$
D(X_0;\omega_{\mathrm{pu}})=
\sum_e |\mu_{e0}(R_0)|^2 |E_{\mathrm{pu}}[\omega_{e0}(R_0)]|^2 \rho_0(X_0),
$$
and windows
$$
W_{\mathrm{GSB}}=\sum_e |\mu_{0e}(R)|^2 |E_{\mathrm{pr}}[\omega_{e0}(R)]|^2,
\quad
W_{\mathrm{SE}}=\sum_e |\mu_{ge}(R)|^2 |E_{\mathrm{pr}}[\omega_{eg}(R)]|^2,
$$
$$
W_{\mathrm{ESA}}=-\sum_{e,f} |\mu_{fe}(R)|^2 |E_{\mathrm{pr}}[\omega_{fe}(R)]|^2.
$$
The minus sign of ESA is explicitly identified as reflecting its opposite phase contribution in heterodyne detection [2603.22671].

The framework extends to dispersed-probe, 2D, 2D-FLEX, and strong-field expressions by introducing electronic dephasing rates \((\nu)\), Fourier transforms in detection times, or non-linear Rabi phase factors [2509.03917]. In the WaveMixings.jl presentation, the same doorway-window logic is stated to apply to dispersed pump-probe, 2D spectroscopy, 2D-FLEX, and strong-field formulas by changing the doorway and window definitions in a controlled way [2509.03917].

The 2026 generalization broadens the method to multi-pulse spectroscopies [2603.22671]. For pump-push-probe (PPP), the preparation consists of a pump step followed by a push step, so the doorway becomes
$$
D_{\mathrm{ppp}}(X_0)
=
\sum_{e,f}
|\mu_{e0}(R_0)|^2
|E_{\mathrm{pu}}[\omega_{e0}(R_0)]|^2
|\mu_{fe}(R_1)|^2
|E_{\mathrm{push}}[\omega_{fe}(R_1)]|^2
\rho_0(X_0),
$$
with \(R_1=R(t=T_1;X_0)\), and the probe windows are evaluated at \(R_2=R(t=T_1+T_2;X_0)\). The signal is then
$$
S_{\mathrm{PPP}}(T_2;\omega_{\mathrm{pr}})
=
\sum_k
\int dX_0\,
D_{\mathrm{ppp}}(X_0)\,
W_k[\Phi_{T_2}(\Phi_{T_1}(X_0));\omega_{\mathrm{pr}}]
$$
[2603.22671].

For pump-induced two-dimensional spectroscopy (P-2D), the quasi-classical doorway-window form is
$$
S_{\mathrm{P\mbox{-}2D}}(\omega_1,\omega_3;T_2)
=
\sum_k
\int dX_0\,
D_{\mathrm{2D}}(X_0;\omega_1)\,
W^{\mathrm{2D}}_k[\Phi_{T_2}(X_0);\omega_3],
$$
where \(D_{\mathrm{2D}}\) and \(W^{\mathrm{2D}}_k\) are analytic expressions almost identical to those in standard 2D spectroscopy, but with the reference ensemble \(X_0\) taken from the pump-prepared distribution at \(t=T_1\) [2603.22671]. This suggests that the doorway-window architecture is not tied to a specific pulse count; rather, it provides a modular factorization pattern whose complexity grows through nested preparation stages and modified detection functionals.

## 5. Computational realization in on-the-fly simulations

The approximation is explicitly intended for on-the-fly simulations driven by ab initio electronic-structure data [2509.03917]. In WaveMixings.jl, the required input for each trajectory \(j\) is the time series of adiabatic energies \(E_n(R_j(t))\) and transition dipoles \(\mu_{mn}(R_j(t))\), together with the sampling weights [2509.03917]. The high-level workflow for an integral pump-probe signal consists of seven steps: generation of \(N_{\mathrm{MC}}\) initial conditions from the Wigner distribution, sampling of the initial excited state \(e\) from the pump-weighted transition probability, doorway evaluation, ground-state propagation, GSB-window assembly, excited-state propagation on manifold \(\{I\}\) using Tully’s surface hopping, and SE/ESA window evaluation, followed by Monte Carlo averaging,
$$
S_{\mathrm{int}}(T,\omega_{\mathrm{pr}})
=
\mathrm{Re}\Big[
\frac{1}{N_{\mathrm{MC}}}
\sum_{i=1}^{N_{\mathrm{MC}}}
D_i\,
\big(
W_{\mathrm{GSB}}^{(i)}(T)
+
W_{\mathrm{SE}}^{(i)}(T)
-
W_{\mathrm{ESA}}^{(i)}(T)
\big)
\Big]
$$
[2509.03917].

The Julia implementation is organized around a main composite type,
```julia
struct WaveMixing{T<:Real}
    pump_freq::Vector{T}
    probe_freq::Vector{T}
    pulse_dur::T
    trajectories::Vector{Traj}
    …
end
```
with modules for data import, constructors and unit conversion, post-processing, and core signal evaluation [2509.03917]. The package contains parsing of surface-hopping outputs, filtering of completed trajectories, direct import of preprocessed arrays of energies and dipoles, unit converters, pulse-shape generators for Gaussian, \(\mathrm{sech}^2\), and Lorentzian profiles, normalization, interpolation and smoothing, difference maps, orientational averages, and plotting routines [2509.03917]. Core functions include pathway-resolved and modality-resolved methods such as `gsb!`, `se!`, `esa!`, `gsb_dispersion!`, `gsb_strong!`, `gsb_2D!`, `se_2D!`, `esa_2D!`, `flex_2D!`, and `trf!` [2509.03917].

The implementation details are paired with explicit performance choices: all inner loops are written in type-stable Julia code; hot functions are annotated with `@inbounds` and `@simd`; heavy linear algebra uses BLAS from `LinearAlgebra.jl`; multiple dispatch on `WaveMixing` and signal types eliminates `if/else` at runtime; and external dependencies are minimized to `LinearAlgebra` and `DelimitedFiles` from the standard library [2509.03917]. The package is open-source on Codeberg.org, and its docstrings are accessible via Julia help mode [2509.03917].

The on-the-fly protocol described for the heptazine\(\cdots\)H\(_2\)O complex provides a complementary ab initio realization [2603.22671]. It starts by sampling the \(T=0\) harmonic Wigner distribution of the electronic ground state using an MP2 Hessian, computing vertical excitation energies and transition dipoles to manifold I at the ADC(2) level with basis cc-pVDZ(+aug), selecting pump-resonant trajectories weighted by \(|\mu_{e0}|^2|E_{\mathrm{pu}}[\omega_{e0}]|^2\), and propagating them with velocity-Verlet plus Landau-Zener surface hopping using the Belyaev-Lebedev formula [2603.22671]. For PPP, higher-state energies and transition dipoles to manifold II are evaluated at the push time, push-resonant trajectories are selected, and propagation continues on \(\{I\cup II\}\) before probe windows are accumulated into the signal average [2603.22671]. Practical details are explicitly specified: MP2/cc-pVDZ for the ground state, ADC(2)/aug-cc-pVDZ for excitations, TURBOMOLE for electronic structure, and ZagHop for surface hopping [2603.22671].

## 6. Benchmarks, applications, and interpretive scope

An illustrative benchmark in WaveMixings.jl reproduces the transient absorption pump-probe spectrum of pyrazine from Gelin et al., *J. Chem. Theory Comput.* **17**, 2394 (2021) [2509.03917]. The setup uses 591 completed surface-hopping trajectories, obtained from 600 launched trajectories filtered to 591, with manifold \(\{I\}\) containing 4 states up to 5.2 eV and manifold \(\{II\}\) containing 30 states up to 10 eV [2509.03917]. Gaussian pump and probe pulses with \(\tau=5\) fs are employed; the pump is centered at \(\omega_{\mathrm{pu}}=5.2\) eV with 200 steps over \(\pm 0.02\) eV, and the probe spans 2.0–6.0 eV in 0.02 eV steps [2509.03917].

The reported performance for the full GSB+SE+ESA integral pump-probe signal on 591 trajectories and a \(200\times 200\) \((T\times \omega_{\mathrm{pr}})\) grid is a total wall-clock time of approximately 128 s on a 6-core i5-12500H laptop with Julia 1.13.0-DEV, with a typical memory footprint below 2 GB [2509.03917]. Convergence is quantified: signals converge within 5% RMS error for \(N_{\mathrm{MC}}\gtrsim 500\) trajectories, and doubling to 1,000 trajectories reduces noise to \(\lesssim 2\%\) [2509.03917]. Accuracy is described in pathway-resolved and modality-resolved terms: the package reproduces peak positions and oscillatory GSB and ESA trace features with period \(\sim 33\) fs in perfect agreement with published results; dispersed pump-probe and TRF spectra with dephasing \(\nu=0.001\) eV fully match reference line shapes; 2D nonrephasing and rephasing signals at \(T=60\) fs show correct cross-peak patterns; and 2D-FLEX reveals vibronic doublets at \((\omega_{\mathrm{ex}},\omega_{\mathrm{det}})\approx(4.5\ \mathrm{eV},3.2\ \mathrm{eV})\) in line with theory [2509.03917]. Parameter studies further state that narrowing the pulse duration to \(\tau\to 2\) fs broadens spectral features by \(\sim 50\) meV, that dephasing \(\nu=0.005\) eV yields \(\sim 100\) fs \(T_2\) coherence lifetimes in dispersed signals, and that alternative state-sampling distributions affect early-time GSB amplitude by \(\lesssim 10\%\) [2509.03917].

The heptazine\(\cdots\)H\(_2\)O application demonstrates the extension to multi-pulse spectroscopy [2603.22671]. The protocol yields approximately 300 trajectories on manifold I after pump selection and approximately 287 trajectories on \(\{I\cup II\}\) after push selection, with \(T_1=100\) fs, \(T_2\) up to 100 fs, and \(\Delta t\approx 0.5\) fs [2603.22671]. The reported spectral interpretation is that pump-stimulated experiments provide much richer information on ultrafast radiationless relaxation dynamics than conventional pump-probe and 2D experiments [2603.22671]. Specifically, pump-probe SE and ESA reproduce ultrafast \((<20\ \mathrm{fs})\) \(\pi\pi^\*\to n\pi^\*\) internal conversion and the appearance of a vibronically hot \(S_1(\pi\pi^\*)\) state at \(\sim 2\) eV; PPP shows a push-induced re-excitation of the hot \(S_1\) ensemble followed by rapid cascade back to \(S_1\), with an intense oscillating SE band near 2.8 eV that red-shifts and decays over \(\sim 20\) fs and a broad ESA from 0–4 eV that builds up over \(\sim 30\) fs; and P-2D maps show a \(\sim 3\) eV-wide emission band below the diagonal even at \(T_2=20\) fs, together with broad elongated ESA attributed to population cycling among bright \(S_4\)–\(S_6\) and dark charge-transfer states [2603.22671].

These examples define the practical scope of the approximation. It is presented as numerically efficient because only a single ensemble of classical trajectories is needed and no multi-dimensional integrals over Liouville pathways are required; it is described as conceptually transparent because \(D\) specifies preparation and \(W\) specifies emission or absorption at the probe; and it is stated to be fully compatible with on-the-fly electronic structure for systems with many modes and many electronic states [2603.22671]. A plausible implication is that the method is particularly well suited to studies where electronic-structure cost, rather than formal response-theory algebra, is the limiting factor.

## 7. Relation to neighboring methodologies and common points of confusion

The quasi-classical doorway-window approximation belongs to the broader class of mixed quantum-classical methods for nonlinear spectroscopy, but its distinctive feature is the explicit factorization of the signal into preparation and detection objects evaluated on trajectory ensembles [2509.03917]. In the WaveMixings.jl formulation, this structure is preserved across integral and dispersed pump-probe, 2D, 2D-FLEX, TRF, and strong-field signals by changing the definitions of the relevant doorway and window functions rather than replacing the overall framework [2509.03917]. In the multi-pulse extension, the same logic is carried over to PPP and P-2D by introducing additional preparation stages before the final window evaluation [2603.22671].

A frequent source of confusion is the role of coherence. The approximation does include electronic coherence phase factors during short coherence intervals \(T_1\) and \(T_3\), evaluated at fixed nuclear geometry [2509.03917]. However, it does not retain nuclear and electronic coherences beyond classical transport and therefore cannot describe interference between pathways in the fully quantum sense [2603.22671]. Another common confusion concerns pulse overlap. The factorization relies on short, non-overlapping pulses; signal artifacts at overlap times below 20 fs are explicitly stated to be non-quantitative [2603.22671]. A further point is that the method is not independent of the chosen dynamics engine: its accuracy inherits that of the mixed quantum-classical propagation and the electronic-structure level employed, whether exemplified by Tully’s surface hopping [2509.03917] or by Landau-Zener surface hopping at the ADC(2) level [2603.22671].

Taken together, these features place the quasi-classical doorway-window approximation in a specific methodological niche. It is a reduction of nonlinear optical response theory to phase-space averages that preserves pathway resolution, pulse dependence, and direct compatibility with on-the-fly simulations, while accepting the standard semiclassical compromises associated with frozen-coordinate interactions, classical population transport, and non-overlapping ultrashort pulses [2509.03917].

Source: https://www.emergentmind.com/topics/quasi-classical-doorway-window-approximation