Optical Pump-Probe Reflectance (OPPR)

• Optical Pump-Probe Reflectance (OPPR) is an ultrafast spectroscopy technique that measures pump-induced, time-dependent reflectance changes to probe transient dielectric properties.
• It utilizes a model-independent inversion of differential reflectance data to reconstruct changes in the complex dielectric function and refractive indices without relying solely on empirical fitting.
• OPPR enables detailed studies of carrier dynamics, coherent phonons, and spin-light coupling in various materials with sub-femtosecond time resolution and enhanced spectral sensitivity.

Searching arXiv for recent and foundational OPPR-related papers to ground the article. Optical pump–probe reflectance (OPPR) is an ultrafast optical spectroscopy method in which a pump pulse generates a transient change in a material and a delayed probe pulse measures the corresponding relative change in reflectance, typically written as $\Delta R/R$ or $\Delta R/R_{0}$. Across the literature, OPPR is used to interrogate non-equilibrium electronic dynamics, coherent collective modes, carrier relaxation, spin-light coupling, transient conductivity, and selective-reflection nonlinearities in solids, semiconductors, correlated materials, and atomic vapors (Ashoka et al., 2021, Baldini et al., 2016, Wang et al., 2011, Mondal et al., 2017, Lopez et al., 2024, Sautenkov et al., 2023). A central methodological development is the model-independent inversion of differential reflectance into transient changes of the complex dielectric function, $$\Delta \epsilon(\omega,t)=\Delta\epsilon_{1}(\omega,t)+i\Delta\epsilon_{2}(\omega,t)$$, which is intended to remove ambiguities caused by the simultaneous influence of reflection, transmission, and scattering on conventional pump–probe observables (Ashoka et al., 2021).

1. Definition and observable structure

In its most common form, OPPR records the pump-induced, time-dependent change in reflectance after photoexcitation by a femtosecond pump pulse, for example

$$\Delta R/R(\omega,t)=\frac{R_{\mathrm{pumped}}(\omega,t)-R_{\mathrm{unpumped}}(\omega)}{R_{\mathrm{unpumped}}(\omega)}$$

or equivalently

$$\Delta R/R_{0}\equiv (R_{\mathrm{on}}-R_{\mathrm{off}})/R_{\mathrm{off}}\,,$$

depending on the notation adopted in a given experiment (Baldini et al., 2016, Wang et al., 2011). In reflection geometry, the measured signal is often interpreted as arising from pump-induced changes in the real and imaginary parts of the dielectric function, $\epsilon_{1}$ and $\epsilon_{2}$, or equivalently in the refractive-index components $n$ and $k$ (Ashoka et al., 2021, Talbayev et al., 2014).

For a planar interface, the equilibrium optical response is expressed through Fresnel relations. In the semi-infinite normal-incidence limit from air into a medium,

$R(\omega)=\frac{(n-1)^2+k^2}{(n+1)^2+k^2}\,,$

with

$\Delta R/R_{0}$0

For thin films on substrates, transfer-matrix or analytic expressions $\Delta R/R_{0}$1 are used instead (Ashoka et al., 2021). A linearized treatment of small photoinduced perturbations gives

$\Delta R/R_{0}$2

and therefore

$\Delta R/R_{0}$3

where $\Delta R/R_{0}$4 and $\Delta R/R_{0}$5 are computed from the static dielectric function (Ashoka et al., 2021).

This formulation establishes the basic interpretive problem of OPPR: $\Delta R/R_{0}$6 is not itself a material parameter, but a projection of the transient optical state through equilibrium optics and geometry. The motivation for quantitative OPPR is therefore to reconstruct $\Delta R/R_{0}$7, $\Delta R/R_{0}$8, and $\Delta R/R_{0}$9 rather than to interpret raw reflectance changes phenomenologically (Ashoka et al., 2021).

2. Dielectric-function inversion and quantitative OPPR

A model-independent route to transform differential transmission or reflection spectra into changes in the transient dielectric function is given by a variational Kramers–Kronig procedure (Ashoka et al., 2021). The protocol begins by measuring the equilibrium reflectance $$\Delta \epsilon(\omega,t)=\Delta\epsilon_{1}(\omega,t)+i\Delta\epsilon_{2}(\omega,t)$$0, and if available transmission $$\Delta \epsilon(\omega,t)=\Delta\epsilon_{1}(\omega,t)+i\Delta\epsilon_{2}(\omega,t)$$1, then retrieving $$\Delta \epsilon(\omega,t)=\Delta\epsilon_{1}(\omega,t)+i\Delta\epsilon_{2}(\omega,t)$$2 and $$\Delta \epsilon(\omega,t)=\Delta\epsilon_{1}(\omega,t)+i\Delta\epsilon_{2}(\omega,t)$$3 with a variational dielectric function composed of a small set of Lorentz oscillators plus a dense grid of weak KK-constrained oscillators (Ashoka et al., 2021). This supplies analytic $$\Delta \epsilon(\omega,t)=\Delta\epsilon_{1}(\omega,t)+i\Delta\epsilon_{2}(\omega,t)$$4 and $$\Delta \epsilon(\omega,t)=\Delta\epsilon_{1}(\omega,t)+i\Delta\epsilon_{2}(\omega,t)$$5 and permits analytical differentiation of $$\Delta \epsilon(\omega,t)=\Delta\epsilon_{1}(\omega,t)+i\Delta\epsilon_{2}(\omega,t)$$6 to obtain $$\Delta \epsilon(\omega,t)=\Delta\epsilon_{1}(\omega,t)+i\Delta\epsilon_{2}(\omega,t)$$7 and $$\Delta \epsilon(\omega,t)=\Delta\epsilon_{1}(\omega,t)+i\Delta\epsilon_{2}(\omega,t)$$8 across the measured spectral window.

The transient step is then formulated as a constrained inverse problem. One records $$\Delta \epsilon(\omega,t)=\Delta\epsilon_{1}(\omega,t)+i\Delta\epsilon_{2}(\omega,t)$$9 at selected delays, represents $$\Delta R/R(\omega,t)=\frac{R_{\mathrm{pumped}}(\omega,t)-R_{\mathrm{unpumped}}(\omega)}{R_{\mathrm{unpumped}}(\omega)}$$0 on a dense energy grid of KK-paired basis functions, enforces KK consistency so that $$\Delta R/R(\omega,t)=\frac{R_{\mathrm{pumped}}(\omega,t)-R_{\mathrm{unpumped}}(\omega)}{R_{\mathrm{unpumped}}(\omega)}$$1 is the Hilbert transform of $$\Delta R/R(\omega,t)=\frac{R_{\mathrm{pumped}}(\omega,t)-R_{\mathrm{unpumped}}(\omega)}{R_{\mathrm{unpumped}}(\omega)}$$2, and solves a least-squares problem of the form

$$\Delta R/R(\omega,t)=\frac{R_{\mathrm{pumped}}(\omega,t)-R_{\mathrm{unpumped}}(\omega)}{R_{\mathrm{unpumped}}(\omega)}$$3

with respect to the variational amplitudes (Ashoka et al., 2021). The result is a delay-dependent reconstruction of $$\Delta R/R(\omega,t)=\frac{R_{\mathrm{pumped}}(\omega,t)-R_{\mathrm{unpumped}}(\omega)}{R_{\mathrm{unpumped}}(\omega)}$$4 and $$\Delta R/R(\omega,t)=\frac{R_{\mathrm{pumped}}(\omega,t)-R_{\mathrm{unpumped}}(\omega)}{R_{\mathrm{unpumped}}(\omega)}$$5 without invoking an assumed bandshape or kinetic model.

The same work states that the dielectric perturbation can then be converted to refractive-index changes using

$$\Delta R/R(\omega,t)=\frac{R_{\mathrm{pumped}}(\omega,t)-R_{\mathrm{unpumped}}(\omega)}{R_{\mathrm{unpumped}}(\omega)}$$6

The detailed expression printed in the supplied technical guide for $$\Delta R/R(\omega,t)=\frac{R_{\mathrm{pumped}}(\omega,t)-R_{\mathrm{unpumped}}(\omega)}{R_{\mathrm{unpumped}}(\omega)}$$7 appears in truncated form; however, the text explicitly states that the inversion yields $$\Delta R/R(\omega,t)=\frac{R_{\mathrm{pumped}}(\omega,t)-R_{\mathrm{unpumped}}(\omega)}{R_{\mathrm{unpumped}}(\omega)}$$8 and $$\Delta R/R(\omega,t)=\frac{R_{\mathrm{pumped}}(\omega,t)-R_{\mathrm{unpumped}}(\omega)}{R_{\mathrm{unpumped}}(\omega)}$$9 from the recovered dielectric response (Ashoka et al., 2021). Benchmarking against time-resolved Frequency Domain Interferometry (FDI), which directly measures $$\Delta R/R_{0}\equiv (R_{\mathrm{on}}-R_{\mathrm{off}})/R_{\mathrm{off}}\,,$$0 with sub-femtosecond resolution, shows quantitative agreement within experimental uncertainty (Ashoka et al., 2021). This is the principal basis for treating OPPR as a quantitative dielectric probe rather than only a differential observable.

A worked example is provided for a CsPbBr$$\Delta R/R_{0}\equiv (R_{\mathrm{on}}-R_{\mathrm{off}})/R_{\mathrm{off}}\,,$$1 thin film at $$\Delta R/R_{0}\equiv (R_{\mathrm{on}}-R_{\mathrm{off}})/R_{\mathrm{off}}\,,$$2 fs. At $$\Delta R/R_{0}\equiv (R_{\mathrm{on}}-R_{\mathrm{off}})/R_{\mathrm{off}}\,,$$3 eV, static ellipsometry gives $$\Delta R/R_{0}\equiv (R_{\mathrm{on}}-R_{\mathrm{off}})/R_{\mathrm{off}}\,,$$4 and $$\Delta R/R_{0}\equiv (R_{\mathrm{on}}-R_{\mathrm{off}})/R_{\mathrm{off}}\,,$$5, corresponding to $$\Delta R/R_{0}\equiv (R_{\mathrm{on}}-R_{\mathrm{off}})/R_{\mathrm{off}}\,,$$6 and $$\Delta R/R_{0}\equiv (R_{\mathrm{on}}-R_{\mathrm{off}})/R_{\mathrm{off}}\,,$$7. The static Fresnel reflectance is $$\Delta R/R_{0}\equiv (R_{\mathrm{on}}-R_{\mathrm{off}})/R_{\mathrm{off}}\,,$$8, the derivatives are $$\Delta R/R_{0}\equiv (R_{\mathrm{on}}-R_{\mathrm{off}})/R_{\mathrm{off}}\,,$$9 and $\epsilon_{1}$0, and a measured $\epsilon_{1}$1 is inverted to $\epsilon_{1}$2 and $\epsilon_{1}$3. The derived refractive-index changes are $\epsilon_{1}$4 and $\epsilon_{1}$5, consistent with FDI-measured $\epsilon_{1}$6 at the same delay (Ashoka et al., 2021).

3. Experimental architectures and detection strategies

OPPR instrumentation spans broadband visible reflectivity, spectrally filtered detection, and reflection-geometry THz probing. A representative femtosecond broadband transient reflectivity instrument for solids combines a Ti:Sapphire oscillator at 800 nm with 45 fs pulses, an amplified 3–10 kHz output up to 3 mJ/pulse, pump–probe beam splitting, broadband white-light generation in 3 mm CaF$\epsilon_{1}$7, and collection into a spectrometer coupled to a single-shot Hamamatsu CMOS line array with 1024 pixels (Baldini et al., 2016). In that configuration, the probe spans 1.72–2.92 eV, the effective time resolution is approximately 45 fs, and a single-shot readout synchronized up to 10 kHz resolves coherent oscillations with amplitudes below $\epsilon_{1}$8 in $\epsilon_{1}$9 (Baldini et al., 2016). Group-velocity-dispersion correction is performed by fitting the wavelength dependence of time zero to a polynomial up to third order and shifting each trace by its extracted $\epsilon_{2}$0 (Baldini et al., 2016).

The same apparatus illustrates several standard OPPR design choices: a motorized retroreflector delay line with step size 1 $\epsilon_{2}$1m corresponding to 6.7 fs, pump spot dimensions of about $\epsilon_{2}$2, probe spot dimensions of about $\epsilon_{2}$3, and cryogenic sample control from 8 to 340 K in a closed-cycle He cryostat (Baldini et al., 2016). The raw observables are the pumped and unpumped reflectances, from which the differential reflectivity is computed directly (Baldini et al., 2016).

A distinct implementation is reflection-geometry pump–THz probe spectroscopy. In the system reported for tunable near- to mid-infrared pumping, a Ti:sapphire amplifier at 800 nm and about 35 fs drives a two-output optical parametric amplifier and difference-frequency generation in GaSe to generate pump wavelengths from 3 to 15 $\epsilon_{2}$4m and beyond, while a THz probe is generated in ZnTe and reflected from the sample at about $\epsilon_{2}$5 incidence (Zhang et al., 2017). Detection uses electro-optic sampling in a 1 mm ZnTe crystal. The measured quantity is the transient change in reflected THz electric field, $\epsilon_{2}$6, with normalized response $\epsilon_{2}$7 in the small-change approximation (Zhang et al., 2017). The stated frequency coverage is about 0.25 to 2.5 THz, time resolution is set by the optical gate at about 35 fs, and a $\epsilon_{2}$8 transient change in reflected THz electric field is resolved in EuSbTe$\epsilon_{2}$9 under 15 $n$0m excitation (Zhang et al., 2017).

Spectrally resolved detection can substantially alter OPPR sensitivity. In diamond, transient reflectivity from coherent optical phonons measured with optical band-pass filters shows an approximately 35-fold enhancement in coherent oscillation sensitivity relative to spectrally integrated detection (Nakamura et al., 2016). The theoretical basis is that phonon-induced sidebands at $n$1 and $n$2 enter the reflectivity modulation with opposite sign, producing substantial cancellation upon full spectral integration, whereas filtering isolates one sideband and prevents cancellation (Nakamura et al., 2016). In the reported experiment, the unfiltered oscillation amplitude is of order $n$3, while filtered detection with FB900 yields a much larger oscillatory signal at the same 40.0 THz phonon frequency and about 6 ps damping time (Nakamura et al., 2016).

4. Dynamical regimes accessed by OPPR

Carrier and excitonic dynamics

Few-layer MoS$n$4 provides a canonical example of exciton-sensitive OPPR. Using a 780 nm, approximately 100 fs Ti:Sapphire pump and a tunable 645–700 nm probe derived from an OPO and BBO doubling, differential reflection was measured near an excitonic resonance with pump fluence from about 15 to 100 $n$5 (Wang et al., 2011). The $n$6 spectrum closely resembles the derivative of the resonance, and the best-fit perturbation of a Gaussian resonance involves an approximate red shift of $n$7 nm, a peak amplitude increase of about $n$8, and a linewidth narrowing of about $n$9 (Wang et al., 2011). At fixed wavelength the dynamics are well fitted by a single exponential,

$k$0

with decay times increasing monotonically from about 150 ps on the blue side of resonance to about 400 ps on the red side (Wang et al., 2011). The peak signal is proportional to pump fluence, while the decay time is fluence independent between 15 and 100 $k$1, which the authors interpret as linear carrier dynamics without significant many-body effects at those densities (Wang et al., 2011).

Polaron formation and spin-correlated relaxation

In Eu$k$2Y$k$3MnO$k$4, OPPR at 1.55 and 3.10 eV was used to track Jahn–Teller relaxation, polaron trapping, and localized spin excitation after photoexcitation (Talbayev et al., 2014). The experiment employed 70 fs pulses from a Ti:Sapphire regenerative amplifier at 250 kHz, pump fluence of approximately $k$5 with linearity confirmed down to $k$6, and temperatures from 5 to 300 K (Talbayev et al., 2014). The transient response at low temperature is described by a double-exponential form

$k$7

where the first term parameterizes the rise of a slow component associated with localized spin excitation and the second describes nanosecond-scale recombination (Talbayev et al., 2014). Representative values at 5 K are $k$8 and $k$9 ps at 1.55 eV, and $R(\omega)=\frac{(n-1)^2+k^2}{(n+1)^2+k^2}\,,$0 and $R(\omega)=\frac{(n-1)^2+k^2}{(n+1)^2+k^2}\,,$1 ps at 3.10 eV, with $R(\omega)=\frac{(n-1)^2+k^2}{(n+1)^2+k^2}\,,$2–5 ns (Talbayev et al., 2014). The temperature dependence shows that $R(\omega)=\frac{(n-1)^2+k^2}{(n+1)^2+k^2}\,,$3 and $R(\omega)=\frac{(n-1)^2+k^2}{(n+1)^2+k^2}\,,$4 are essentially zero above about 150 K, begin to grow below about 100–120 K, and do not display sharp anomalies at $R(\omega)=\frac{(n-1)^2+k^2}{(n+1)^2+k^2}\,,$5 K or $R(\omega)=\frac{(n-1)^2+k^2}{(n+1)^2+k^2}\,,$6 K, instead tracking the onset of short-range spin correlations near 100 K (Talbayev et al., 2014).

Coherent phonons and collective modes

Broadband OPPR on La$R(\omega)=\frac{(n-1)^2+k^2}{(n+1)^2+k^2}\,,$7CuO$R(\omega)=\frac{(n-1)^2+k^2}{(n+1)^2+k^2}\,,$8 demonstrates the ability of transient reflectivity to resolve coherent collective modes in strongly correlated quantum systems (Baldini et al., 2016). With a 3.10 eV, 45 fs pump at absorbed fluence 4.4 mJ/cm$R(\omega)=\frac{(n-1)^2+k^2}{(n+1)^2+k^2}\,,$9 and probe continuum from 1.72 to 2.92 eV, the transient map shows a resolution-limited rise of about 50 fs with superimposed oscillations (Baldini et al., 2016). Temporal traces are fitted with multicomponent oscillators containing five frequencies, 3.68, 4.55, 6.93, 8.23, and 13.0 THz with uncertainty $\Delta R/R_{0}$00 THz, and Fourier analysis identifies $\Delta R/R_{0}$01 phonon modes up to about 13 THz in agreement with spontaneous Raman data (Baldini et al., 2016). The extracted amplitude versus probe photon energy is described as a “Raman matrix element” and is used for quantitative electron-phonon coupling analysis (Baldini et al., 2016).

In diamond, the coherent optical phonon at 40.0 THz with period $\Delta R/R_{0}$02 fs and damping time about 6 ps is detected by transient reflectivity, and the experiment explicitly shows that the phase of the oscillatory signal depends on whether the high-energy or low-energy sideband of the probe spectrum is selected (Nakamura et al., 2016). This establishes that not only the magnitude but also the sign and phase of OPPR signals can depend strongly on spectral detection window (Nakamura et al., 2016).

5. Coherence, polarization selection, and magneto-optical effects

OPPR in reflection geometry can detect perturbed free-induction decay (PFID) at negative pump–probe delay. For a two-level ensemble described by optical Bloch equations, a weak probe first generates a free-induction-decay polarization $\Delta R/R_{0}$03, and a later pump perturbs that residual polarization, generating a third-order component $\Delta R/R_{0}$04 that radiates into the probe-reflection direction (Mondal et al., 2018). After monochromator filtering, double chopping, and lock-in detection, the measured spectrally resolved signal takes the form

$\Delta R/R_{0}$05

for negative delay $\Delta R/R_{0}$06, with $\Delta R/R_{0}$07 (Mondal et al., 2018). At resonance, $\Delta R/R_{0}$08, the negative-delay signal becomes a pure exponential $\Delta R/R_{0}$09, permitting direct extraction of the coherence time. In GaAs/AlGaAs multi-quantum wells at 4 K, degenerate 100 fs pulses tuned to the heavy-hole exciton yielded $\Delta R/R_{0}$10 ps, in good agreement with the value inferred from the continuous-wave photoluminescence-excitation linewidth $\Delta R/R_{0}$11 meV via $\Delta R/R_{0}$12 ps (Mondal et al., 2018).

Magneto-optical OPPR introduces an additional polarization-selective mechanism through relativistic spin–light coupling. Starting from the Foldy–Wouthuysen-transformed Dirac–Kohn–Sham Hamiltonian, Mondal et al. identify an Angular Magneto-Electric Hamiltonian

$\Delta R/R_{0}$13

which can be recast as a Zeeman-like coupling to an induced opto-magnetic field $\Delta R/R_{0}$14 (Mondal et al., 2017). In a pump–probe experiment the total field contains pump–pump, probe–probe, and pump–probe cross terms, and the cross term is transiently active during pulse overlap (Mondal et al., 2017). For parallel linear pump and probe polarizations, $\Delta R/R_{0}$15 and $\Delta R/R_{0}$16. For perpendicular linear polarizations, the theory predicts a nonzero opto-magnetic field along $\Delta R/R_{0}$17 proportional to $\Delta R/R_{0}$18, so that at short delays $\Delta R/R_{0}$19 the coherent field grows linearly in $\Delta R/R_{0}$20 (Mondal et al., 2017). The paper states that this explains why a sizable coherent Kerr signal appears in the perpendicular configuration but not in the parallel one, and gives an estimated field of approximately $\Delta R/R_{0}$21T for the experiment of Bigot et al. (Mondal et al., 2017).

A common misconception is that negative-delay or polarization-dependent reflectivity signals necessarily require a more complex nonlinear wave-mixing geometry. The cited works show, more narrowly, that reflection-geometry pump–probe experiments can access PFID and polarization-dependent coherent opto-magnetic fields within pump–probe architectures, provided spectral filtering, polarization control, and lock-in selection are handled explicitly (Mondal et al., 2018, Mondal et al., 2017).

6. Nonlinear depth profiles, selective reflection, and theoretical generalizations

When the photoinduced response depends nonlinearly on incident pump intensity, the depth profile of photoconductivity can no longer be assumed to retain a fixed shape as fluence changes. In the formalism of “Nonlinear photoconductivity in pump-probe spectroscopy. I. Optical coefficients,” the inhomogeneous medium is written for $\Delta R/R_{0}$22 with

$\Delta R/R_{0}$23

and the probe obeys the scalar wave equation

$\Delta R/R_{0}$24

with reflection amplitude determined by the surface logarithmic derivative $\Delta R/R_{0}$25 (Lopez et al., 2024). The paper gives analytic profiles for linear exponential photoconductivity, saturable profiles, two-photon-absorption (TPA) profiles, and combined saturable+TPA profiles (Lopez et al., 2024). Depending on the profile, the wave equation reduces to the Bessel equation, Gauss hypergeometric equation, or Heun equation, yielding exact analytic expressions for the optical coefficients (Lopez et al., 2024). In the thin-film limit, a saturable profile produces $\Delta R/R_{0}$26, while the TPA thin-film form is stated to be indistinguishable from linear behavior in $\Delta R/R_{0}$27 at low fluence (Lopez et al., 2024). This provides a systematic framework for constraining systematic uncertainty associated with nonlinear photoconductivity profile distortion (Lopez et al., 2024).

A different nonlinear reflectance regime is selective reflection from a YAG window–high-density rubidium vapor interface. There, the reflection coefficient is written as

$\Delta R/R_{0}$28

with $\Delta R/R_{0}$29 for YAG at 780 nm and $\Delta R/R_{0}$30 for the vapor (Sautenkov et al., 2023). For $\Delta R/R_{0}$31, the resonant change in reflectance is proportional to $\Delta R/R_{0}$32, making the signal dispersion-like (Sautenkov et al., 2023). The linear susceptibility near the Rb D$\Delta R/R_{0}$33 resonance is approximated by a Lorentzian with linewidth dominated by dipole–dipole broadening,

$\Delta R/R_{0}$34

so that for $\Delta R/R_{0}$35–$\Delta R/R_{0}$36 the broadening spans 13.2–39.6 GHz (Sautenkov et al., 2023). Pump intensities from 1.2 to 8.8 kW$\Delta R/R_{0}$37cm$\Delta R/R_{0}$38 reduce both resonance magnitude and width through resonant saturation, with the measured selective-reflection width collapsing from about 16 GHz to about 8 GHz at the lowest density and from about 42 GHz to about 20 GHz at the highest density over the reported intensity range (Sautenkov et al., 2023). For $\Delta R/R_{0}$39, the resonance amplitude decreases from about $\Delta R/R_{0}$40 at zero pump to about $\Delta R/R_{0}$41 at 8.8 kW$\Delta R/R_{0}$42cm$\Delta R/R_{0}$43 (Sautenkov et al., 2023). The paper suggests all-optical modulation as an application, with estimated modulation depth 50–60% and response times 4–15 ps inferred from $\Delta R/R_{0}$44 (Sautenkov et al., 2023).

On the theoretical side, a generalized two-time linear-response theory for pumped systems extends OPPR analysis to generic pump–probe delays and probe frequencies (Eskandari-asl et al., 2024). In that framework, the total Hamiltonian is split into a pumped part $\Delta R/R_{0}$45 and a weak probe coupling $\Delta R/R_{0}$46, leading to a two-time susceptibility $\Delta R/R_{0}$47 and a two-time optical conductivity $\Delta R/R_{0}$48 (Eskandari-asl et al., 2024). The out-of-equilibrium dielectric function is then

$\Delta R/R_{0}$49

and for $\Delta R/R_{0}$50-polarized light at incidence angle $\Delta R/R_{0}$51 the reflectance is approximated in the XUV by the usual Fresnel formula involving $\Delta R/R_{0}$52 (Eskandari-asl et al., 2024). The work emphasizes that inter-band transitions dominate sharp resonance features, intra-band motion contributes smoother background modulations, and inclusion of a local dipole term introduces odd-order pump-frequency oscillations in $\Delta R/R_{0}$53, whereas pure Peierls coupling yields only even harmonics (Eskandari-asl et al., 2024). The associated Dynamical Projective Operatorial Approach permits one-shot pump-only time propagation followed by postprocessing for all pump–probe delays and probe spectra (Eskandari-asl et al., 2024).

7. Assumptions, limitations, and methodological practice

Quantitative OPPR depends on a sequence of approximations and calibration choices that are explicit in the cited works. The dielectric inversion of $\Delta R/R_{0}$54 assumes the perturbative limit, namely that $\Delta R/R_{0}$55 remains small enough for first-order Taylor expansion to hold (Ashoka et al., 2021). It also assumes that KK reconstruction is sufficiently accurate within a finite measured spectral window, which requires that evolving features be local and well bounded or, as a practical recommendation, that the spectral range include the full feature and its tails beyond band edges (Ashoka et al., 2021). Errors in the static dielectric baseline propagate directly into $\Delta R/R_{0}$56 and $\Delta R/R_{0}$57, and neglected Fabry–Pérot fringes, film-thickness uncertainty, or inhomogeneities can distort the inversion if not built into $\Delta R/R_{0}$58 (Ashoka et al., 2021). Strong spectral shifts or high-fluence effects such as bandgap renormalization may require repeated expansion around a shifted baseline or inclusion of higher-order terms (Ashoka et al., 2021).

Instrumental practices likewise determine reliability. The broadband OPPR guide recommends high static-spectrum signal-to-noise ratio, stated as $\Delta R/R_{0}$59 counts, careful chirp calibration, homogeneous excitation densities, and fluence-linearity tests (Ashoka et al., 2021). The 45 fs broadband apparatus reports a noise floor of $\Delta R/R_{0}$60 RMS after 1000 shots per integration, at least 50 repeats per matrix, and binning of about 5–20 pixels in photon energy (Baldini et al., 2016). In THz reflection geometry, the principal error sources are laser intensity and timing jitter, beam-pointing drift, residual water-vapor absorption, detector-crystal alignment, and pump-induced sample heating or damage at high fluence (Zhang et al., 2017). In PFID measurements, the signal may vanish under spectral integration, so monochromator selection is not incidental but essential to observing the negative-delay coherence term (Mondal et al., 2018). In coherent-phonon reflectivity, spectral integration can suppress the oscillatory signal by destructive interference between sidebands, whereas band-pass filtering can enhance the detected oscillation amplitude by a factor of about 35 (Nakamura et al., 2016).

Several broader implications follow from the literature but should be stated cautiously. The collected results suggest that OPPR is not a single fixed observable but a family of reflectance-based measurements whose interpretation depends on geometry, spectral resolution, polarization configuration, depth profile, and inversion model. They also suggest that reflectance-only data can be elevated from qualitative kinetics to quantitative dielectric reconstruction when equilibrium optical constants, derivative optics, KK constraints, and independent benchmarks such as FDI are available (Ashoka et al., 2021). A plausible implication is that future OPPR work will continue to converge with transfer-matrix modeling, two-time nonequilibrium response theory, and spectrally selective detection rather than relying only on empirical fitting of raw $\Delta R/R_{0}$61 traces (Ashoka et al., 2021, Eskandari-asl et al., 2024, Lopez et al., 2024).