---
title: Photonic Doppler Velocimetry Fundamentals
url: https://www.emergentmind.com/topics/photonic-doppler-velocimetry-pdv
type: topic
---

# Photonic Doppler Velocimetry Fundamentals

Searching arXiv for recent and foundational PDV papers to support the article.
arxiv_search(query="Photonic Doppler Velocimetry PDV time lens ejecta Bayesian", max_results=10, sort_by="relevance")
Photon or Photonic Doppler Velocimetry (PDV) is a fiber-based heterodyne diagnostic that measures the frequency shift of light backscattered from moving targets and converts that shift to an apparent velocity through the Doppler relation. In its standard form, a local oscillator or stable reference field is mixed with the Doppler-shifted return from a moving surface or particulate ejecta, and the resulting beat signal is analyzed in time–frequency space, typically by a Short-Time Fourier Transform (STFT), to infer a velocity history. Across recent work, PDV appears in three closely connected forms: as a conventional single-channel fiber interferometer for fast free-surface or shock-front measurements, as a transport-sensitive diagnostic for dense ejecta where multiple scattering invalidates naïve single-scattering interpretations, and as a platform for methodological extensions such as time-lens temporal magnification and Bayesian time-domain inference [2101.02119] [2312.04256] [2508.13695].

## 1. Measurement principle and Doppler mapping

In a single-channel fiber PDV, the probe delivers light to the target and also returns a small, non-Doppler-shifted reference reflection; the returning reference and Doppler-shifted light are innately mixed, routed by a circulator to a photodetector. The detector measures
\[
I(t) = I_0 + I_d + \sqrt{I_0 I_d}\,\sin\!\big(2\pi f_b(t)t + \phi\big),
\]
where \(I_0\) and \(I_d\) are the intensities of the reference and Doppler return, respectively, and \(\phi\) is their relative phase. The beat frequency is \(f_b(t)=|f_d(t)-f_0|\), with \(f_0\) the laser frequency and \(f_d(t)\) the Doppler-shifted return [2101.02119].

For a moving mirror, the beat frequency–velocity relation is
\[
f_b(t)=\pm 2\,[v(t)/c]\,f_0.
\]
For normal incidence in air, using \(f_0=c/\lambda_0\), this reduces to
\[
f_b(t)=\frac{2v(t)}{\lambda_0}.
\]
For non-normal incidence at angle \(\theta\),
\[
f_b(t)=\frac{2v(t)\cos\theta}{\lambda_0},
\]
and in a medium of refractive index \(n\),
\[
f_b(t)=\frac{2nv(t)}{\lambda_0}.
\]
At \(\lambda_0=1550\) nm, \(v=1\) km/s gives \(f_b\approx 1.29\) GHz [2101.02119].

A formally equivalent backscatter relation is used in ejecta work:
\[
v=\frac{\lambda_0}{2}f_b,\qquad \Delta f=f_b=\frac{2v}{\lambda_0},
\]
or, in angular-frequency form,
\[
v=\frac{\omega\lambda_0}{4\pi}.
\]
In the single-scattering limit, spectrogram lines therefore map directly to particle or surface velocities along the probe axis [2406.14578].

A more general vector Doppler formula is required when geometry departs from monostatic backscatter or when multiple scattering is present:
\[
\Delta\omega=(\mathbf{k}_s-\mathbf{k}_i)\cdot\mathbf{v},
\]
with \(|\mathbf{k}|=n\,2\pi/\lambda_0\). In the monostatic backscatter case, \(\mathbf{k}_s=-\mathbf{k}_i\), giving \(\Delta f=2nv/\lambda_0\). This relation remains useful as an apparent-velocity mapping even when the detected photons have undergone multiple scattering, but in that regime it no longer has the direct one-photon/one-velocity meaning of the single-scattering case [2312.04256].

## 2. Interferometric architectures and analysis workflows

The standard instrumentation described across the cited work is a fiber-interferometric heterodyne system operating near \(\lambda_0=1.55~\mu\mathrm{m}\). Typical implementations use a continuous-wave or narrow-line laser, a fiber probe, a circulator, a fast photodetector, and a high-bandwidth oscilloscope or digitizer. In the time-domain Bayesian study, typical systems use near-infrared lasers at \(\lambda \approx 1550\) nm and GHz-bandwidth digitizers such as 25 GHz, 50 GS/s instruments; in the TL-PDV study, the PDV operates in the C-band around \(\lambda_0=1550\) nm \((f_0\approx 194~\mathrm{THz})\), which is well matched to four-wave-mixing time-lens technology [2508.13695] [2101.02119].

The established analysis workflow forms a spectrogram by applying an STFT to the oscilloscope trace. The time–frequency content is represented as
\[
S(t,\omega)=\left|\int \mathcal{I}(\tau)\,w(\tau-t)\,e^{i\omega\tau}\,d\tau\right|^2,
\]
with gate window \(w(t)\) of width \(T_w\). In practical PDV, the window form, duration, and separation determine the trade-off between temporal and frequency resolution. In the validation reported for the Bayesian method, an STFT with a 5 ns window yielded a velocity uncertainty of \(0.1275\ \mathrm{km\,s^{-1}}\); the same work emphasizes that STFT windowing broadens discontinuities and limits recovery at very early times [2406.14578] [2508.13695].

A distinct line of work replaces ridge-picking on a spectrogram by direct inference on the oscilloscope trace. The forward model used in the Bayesian time-domain method writes the mean signal as
\[
\mu(t)=\sin\!\left[\phi_0+\int_{t_0}^{t}f(t')\,dt'\right],
\]
with
\[
f(t)=\left|\nu_\mathrm{T}-\nu_\mathrm{R}+\frac{2\nu_\mathrm{T}n_0}{c}v(t)\right|.
\]
Velocity is parameterized as a one-dimensional piecewise-linear time series in acceleration segments, and inference is performed with the Bilby library and the Dynesty nested sampling engine using 1500 live points. The method recovered injected velocities from synthetic data within the 95% credible interval, and on quartz shock data it produced velocities consistent with STFT-derived values while interpolating across low-signal regions [2508.13695].

This division between STFT-based and time-domain inference is methodological rather than conceptual. Both are built on the same heterodyne beat physics. A plausible implication is that PDV should be regarded less as a single analysis protocol than as a measurement class in which the optical hardware, transport model, and inverse method can vary substantially while the Doppler encoding remains the common core.

## 3. From single scattering to radiative transport in ejecta

In ejecta applications, PDV has traditionally been used in the single-scattering regime to read out velocity distributions. That approximation breaks down when the optical thickness \(b\) becomes large:
\[
b=\int \frac{dz}{\ell_s(z)},
\]
where \(\ell_s\) is the scattering mean free path. The 2024 ejecta study states that velocity-only readout is reliable when \(b\ll 1\), with \(b\lesssim 0.3\) given as a practical confidence level; the 2023 multiple-scattering study shows that a distinct peak at the free-surface frequency appears for \(b\le 4\) and is not visible for larger \(b\). These are complementary statements about different analysis criteria, not contradictory thresholds [2406.14578] [2312.04256].

The rigorous formulation connects the measured spectrogram to the specific intensity \(I(\mathbf{r},\mathbf{u},t,\omega)\), the radiative flux density at position \(\mathbf{r}\), direction \(\mathbf{u}\), time \(t\), and frequency \(\omega\). The spectrogram–intensity relation is
\[
\delta(k-k_R)\,S(t,\omega)=T_w|\mathcal{A}_0|^2\int_G\big[I(\mathbf{r},\mathbf{u},t,\omega_0+\omega)+I(\mathbf{r},\mathbf{u},t,\omega_0-\omega)\big]\;\mathbf{u}\cdot\mathbf{n}\;d\mathbf{u}\,d\mathbf{r},
\]
where \(G\) is the detector etendue and \(\mathbf{n}\) is the detector normal. The measured PDV spectrogram is therefore an etendue-integrated specific intensity at frequencies \(\omega_0\pm\omega\), scaled by the STFT window and reference amplitude [2312.04256].

The generalized radiative transfer equation (RTE) used for dynamic ejecta is
\[
\left[\frac{1}{v_E(\mathbf{r},t,\omega)}\frac{\partial}{\partial t}+\mathbf{u}\cdot\nabla_{\mathbf{r}}+\frac{1}{\ell_e(\mathbf{r},t,\omega)}\right]I(\mathbf{r},\mathbf{u},t,\omega)
=
\frac{1}{\ell_s(\mathbf{r},t,\omega)}
\int p(\mathbf{r},\mathbf{u},\mathbf{u}',t,\omega,\omega')\,I(\mathbf{r},\mathbf{u}',t,\omega')\,d\mathbf{u}'\,\frac{d\omega'}{2\pi},
\]
with energy velocity \(v_E\), extinction mean-free path \(\ell_e\), scattering mean-free path \(\ell_s\), and phase function \(p\). The dependence on \((\omega,\omega')\) captures inelastic scattering, specifically Doppler coupling [2312.04256] [2406.14578].

The size distribution enters through the optical coefficients:
\[
\frac{1}{\ell_e(\mathbf{r},t,\omega)}=\int \rho(\mathbf{r},t)\,\sigma_e(a,\omega)\,h(\mathbf{r},t,a)\;da,\qquad
\frac{1}{\ell_s(\mathbf{r},t,\omega)}=\int \rho(\mathbf{r},t)\,\sigma_s(a,\omega)\,h(\mathbf{r},t,a)\;da,
\]
where \(\rho\) is particle number density and \(\sigma_e,\sigma_s\) are Mie-theory extinction and scattering cross sections. The phase function couples angle and frequency through
\[
\delta\!\Big[\omega'-\omega-k_R(\mathbf{u}'-\mathbf{u})\cdot\mathbf{v}\Big],
\]
so the spectrogram reflects a convolved size–velocity field rather than a direct backscatter velocity map [2406.14578].

This multiple-scattering literature corrects a common misconception: PDV is not intrinsically a velocity-only diagnostic. In optically thick ejecta, the spectrogram morphology depends on particle size statistics through \(\ell_e\), \(\ell_s\), and \(p\), and particle-size information can be recovered by forward modeling when the transport physics is included [2406.14578].

## 4. Vacuum ejecta spectra and closed-form parameter extraction

For vacuum ejecta, one line of work combines full Monte Carlo light transport with a closed-form single-scattering interpretation. The GPU-accelerated Monte Carlo reconstruction models photons propagating through a granular ejecta layer, with scattering and absorption probabilities computed from Mie theory. The algorithm samples Beer–Lambert penetration, scattering versus absorption, directional scattering from the Mie phase function, and Doppler frequency updates at scattering events; outgoing photons are histogrammed to obtain a synthetic PDV spectrum. Using \(10^7\) photons, GPU acceleration up to \(\times 20\) on an Nvidia GTX960 reduced runtime to \(\sim 500\) s per case [2012.00963].

Within the single-scattering interpretation introduced in the same work, the dominant PDV peak can be related directly to ejecta parameters. The peak position satisfies
\[
\frac{\tilde v_{\text{peak}}}{v_{fs}} = 1 + \frac{\ln(\tau_0)}{\beta},
\]
and the peak curvature satisfies
\[
\frac{I''(\tilde v_{\text{peak}})}{I(\tilde v_{\text{peak}})} = -\left(\frac{\beta}{v_{fs}}\right)^2.
\]
These relations enable extraction of the velocity-profile coefficient \(\beta\) from curvature and optical thickness \(\tau_0\) from peak position, provided \(v_{fs}\) is known. With independent size information, \(\tau_0\) can then be converted to areal mass \(m_0\) [2012.00963].

The paper validates this procedure on vacuum ejecta from shock-loaded grooved Sn at \(P_{SB}\approx 28\) GPa with \(v_{fs}\approx 2013\) m/s. From the processed PDV spectrum over a \(0.2\)–\(0.8~\mu\)s window, the main peak occurs at \(\tilde v/v_{fs}\approx 1.32\) with curvature \(\approx -110\) in velocity-normalized units, yielding
\[
\beta \approx 10.5,\qquad \tau_0 \approx 28.79.
\]
Independent Mie-scattering measurements gave \(d_m=0.6~\mu\mathrm{m}\) and \(\sigma=0.5\), from which the areal mass was inferred as \(m_0\approx 7.5~\mathrm{mg/cm^2}\). The PDV-derived velocity profile and the piezoelectric probe mass–velocity function were reported to be in good quantitative agreement [2012.00963].

The same study also makes clear where the closed-form model ceases to be adequate. It reproduces the dominant peak well for moderate-to-large optical thickness in vacuum, but it does not reproduce the free-surface peak that emerges when optical thickness is small. Full multiple-scattering Monte Carlo is therefore preferred when the spectrum exhibits multiple peaks, pronounced low-velocity structure, or strong sensitivity to transport anisotropy [2012.00963].

## 5. Ejecta in gas: joint recovery of velocity and size information

The 2024 ejecta work extends PDV from vacuum single-scattering analysis to ejecta transported in gas, where drag, breakup, and multiple scattering alter the spectrogram substantially. The experimental context is a launch tube with barrel inner diameter \(\Phi_{\text{barrel}}=98\) mm, a grooved tin disk with surface grooves of \(60~\mu\mathrm{m}\times 8~\mu\mathrm{m}\), and a copper flyer impacting tin at \(1650~\mathrm{m/s}\) to produce \(P_{\text{shock}}=29.5\) GPa. Ejecta transport is studied in vacuum \((10^{-5}\,\mathrm{bar})\), helium \((5\,\mathrm{bar})\), and air \((1\,\mathrm{bar})\). The PDV is on-axis and in reflection, with a single probe that both illuminates and collects, using \(\lambda_0=1.55~\mu\mathrm{m}\) in the analysis [2406.14578].

The initial ejecta is modeled as spherical particles with joint size–velocity distribution
\[
g(a,\mathbf{v})=h(a)\,j(\mathbf{v}),
\]
with the velocity distribution derived from the integrated ejected mass–velocity curve
\[
M(v)=M_s\exp\!\left[-\beta\left(\frac{v}{v_s}-1\right)\right],\qquad
M_s=12\,\mathrm{mg/cm^2},\ \beta=11.7,\ v_s=v_{\min}=2060\,\mathrm{m/s},\ v_{\max}=3350\,\mathrm{m/s}.
\]
In gas, particles experience drag and possibly breakup. The simulations use the KIVA-II quadratic drag law and a breakup criterion based on the Weber number
\[
\mathrm{We}_p=\frac{2a_p\,\rho_g\,|\mathbf{v}_p-\mathbf{v}_g|^2}{\sigma_p},
\]
with breakup when \(\mathrm{We}_p>15\) [2406.14578].

The forward model couples hydrodynamics to optics. Ejecta transport is simulated with the CEA Phenix code, typically with \(\sim 2500\) numerical particles and 180 time steps of \(\delta t=0.16~\mu\mathrm{s}\), requiring \(\sim 2\) h on 1 AMD EPYC 7763 with 64 cores. For optical modeling, the launch tube is discretized into \(\sim 100\) layers, Mie routines provide \(\sigma_e,\sigma_s\), and the generalized RTE is solved by Monte Carlo random walks. A typical spectrogram involves \(\sim 5.12\times 10^9\) draws over 180 times and 2500 particles, with compute time \(\sim 1\) h 20 min on 80 AMD EPYC 7763 CPUs [2406.14578].

The inversion strategy is iterative forward modeling rather than a closed-form inversion. In vacuum, a power-law \(h(a)\) with \(\alpha=5\), \(a_{\min}=1~\mu\mathrm{m}\), \(a_{\max}=6~\mu\mathrm{m}\) yielded excessive optical thickness \((b\approx 16)\) and hid the free surface; replacing it by a lognormal
\[
h(a)=\frac{1}{a\,\sigma\sqrt{2\pi}}\exp\!\left[-\frac{\ln^2(a/a_0)}{2\sigma^2}\right],\qquad
\sigma=0.5,\ a_0=2.25~\mu\mathrm{m},
\]
truncated to \([1,6]~\mu\mathrm{m}\), reduced \(b\) to \(\sim 6\) and restored the free-surface return. In helium, increasing drag coefficients \(C_d\) within the KIVA-II model aligned the simulated upper envelope with experiment. In air, drag plus breakup produced an early plateau near \(2800~\mathrm{m/s}\) at \(4\)–\(6~\mu\)s, suppression of the free-surface return by \(\sim 8~\mu\)s, and re-acceleration bands around \(1800~\mathrm{m/s}\) between \(10\) and \(15~\mu\)s; simulations captured these qualitative features but retained a broader long-term velocity spread than experiment, indicating limitations of the independence assumption \(g(a,\mathbf{v})=h(a)j(\mathbf{v})\) and suggesting correlated size–velocity initial distributions [2406.14578].

A second common misconception is therefore that multiple scattering merely adds nuisance broadening. The gas-transport analysis shows that spectrogram features such as free-surface disappearance, dynamic-range changes, slowing-down envelopes, and re-acceleration bands carry information about size evolution and breakup as well as velocity.

## 6. Bandwidth extension, dynamic range, and emerging PDV modalities

Conventional PDV is limited by the maximum detectable beat frequency set by the electronics. If \(f_{b,\max}\approx f_s/2\) by Nyquist, then the maximum measurable velocity is approximately
\[
v_{\max}\approx \left(\frac{\lambda_0}{2}\right)f_{b,\max}.
\]
At \(\lambda_0=1550\) nm, \(f_{b,\max}=64\) GHz corresponds to \(v_{\max}\approx 50\) km/s, described as consistent with current PDV limits in the time-lens study [2101.02119].

Time-lens Photon Doppler Velocimetry (TL-PDV) addresses this bandwidth bottleneck by inserting a temporal imaging system into the optical path between the PDV mixing stage and the photodetector. The time lens is a four-wave-mixing (FWM) device that applies a quadratic temporal phase
\[
\phi(t)=\frac{1}{2}Kt^2,
\]
or, in the notation used in the paper,
\[
\phi_f(\tau)=\frac{\tau^2}{2\phi_f''}=-\frac{\tau^2}{\phi_p''},
\]
with the phase realized by FWM between a linearly chirped pump and the input signal in a nonlinear medium such as a highly nonlinear fiber or integrated photonic waveguide [2101.02119].

Temporal imaging comprises three stages: pre-dispersion of the input waveform, the FWM time lens, and post-dispersion of the idler. The imaging condition is
\[
\frac{1}{\phi_1''}+\frac{1}{\phi_2''}=-\frac{1}{\phi_f''},
\]
with temporal magnification
\[
M=-\frac{\phi_2''}{\phi_1''}.
\]
Because temporal magnification stretches time by \(M\), instantaneous frequency scales inversely:
\[
f_{\text{out}}(t)=\frac{f_{\text{in}}(t)}{M}.
\]
The PDV beat frequency at the output is therefore
\[
f_{b,\mathrm{TL}}(t)=\pm \frac{2}{M}\,[v(t)/c]\,f_0,
\]
and velocity is recovered through
\[
v(t)=\left(\frac{\lambda_0}{2}\right)M f_{b,\mathrm{TL}}(t)
\]
for normal incidence in air [2101.02119].

The simulated TL-PDV implementation uses a mode-locked laser pump with 0.5 ps pulse width at 100 MHz repetition rate, dispersion parameter
\[
D= - \left(\frac{2\pi c}{\lambda^2}\right)k_2 = 17~\mathrm{ps\,nm^{-1}\,km^{-1}}
\]
for all fibers, pre-dispersion \(\xi_1=1.1\xi\), post-dispersion \(\xi_2=11\xi\), and a pump chirped by a TOD-limited dispersive fiber of length \(\xi_p=\xi/2\) with \(\xi=4163\) m. These parameters satisfy the imaging condition and yield \(M=10\) [2101.02119].

Under that magnification, a PDV beat spanning \(\sim 120\) GHz for velocities sweeping from \(+100\) km/s to \(-100\) km/s at 1550 nm is reduced to \(\sim 12\) GHz, bringing the signal within the detection band of “10s of GHz” digitizers. The effective velocity range extension scales approximately as
\[
v_{\max,\mathrm{TL}}\approx M\left(\frac{\lambda_0}{2}\right)f_{b,\max},
\]
so a PDV limited to \(12\) GHz can, with \(M=10\), measure original beat frequencies up to \(120\) GHz [2101.02119].

The time-lens study compares TL-PDV with heterodyne electrical downconversion, leapfrog PDV, and time-stretched PDV. The stated distinctions are that electrical downconversion may introduce spectral artifacts and can reduce sensitivity at low velocities, leapfrog PDV extends range at the cost of complexity and expense, and time-stretched PDV creates replicas and requires complex optical hardware. TL-PDV instead magnifies the optically encoded PDV signal directly. This suggests that bandwidth extension in PDV can be pursued either electronically, optically, or inferentially: TL-PDV changes the signal before detection, whereas the Bayesian time-domain method changes the inverse problem after detection [2101.02119] [2508.13695].

Across these modalities, the principal limitations remain explicit. For TL-PDV they include pump depletion, phase mismatch, dispersion errors, nonlinear phase noise, third-order dispersion, and overlap of adjacent magnified windows. For Bayesian inference they include amplitude-model misspecification, multi-surface or multi-velocity reflections, down-shift branch ambiguity, and computational cost, with runs often taking many hours to converge. In both cases, the limiting factor is not the Doppler relation itself but the degree to which the full signal-generation model captures the experiment [2101.02119] [2508.13695].

Source: https://www.emergentmind.com/topics/photonic-doppler-velocimetry-pdv