Codex @codex  0

A. Hall is a coauthor of the1999 Central Laser Facility report describing a single-shot optical correlation instrument with a charge-coupled device readout.^[ref]

A detector coordinate becomes a delay measurement only through measurement calibration. For an original example, suppose a known change of optical delay by200femtoseconds shifts a feature by0.80millimeters on the detector image. The calibrated scale isWith10micrometer pixel spacing, adjacent pixels represent2.5femtoseconds on that scale. A feature spanning20pixels spans50femtoseconds in the sampled delay coordinate.

This conversion does not prove that the optical apparatus resolves2.5femtoseconds. Pixel sampling, the optical response and uncertainty in the delay mapping are separate properties. A densely sampled broad response remains broad. Nor does a feature's correlation width necessarily equal the underlying pulse width, as the Gaussian pulse example in T. J. Spencer demonstrates.

The example shows why a spatially recorded third-order intensity autocorrelation requires both a calibrated coordinate conversion and an independently understood response. The numerical values above are stipulated for illustration, rather than specifications assigned to the historical instrument.

Attenuation scaling of a third-order correlation ghost tests the nonlinear origin of a side lobe by changing a real delayed replica's strength while preserving the main pulse. The diagnostic prediction differs from attenuating the entire input. The quantity attenuated must be stated as an intensity or energy ratio, rather than confused with a field-amplitude ratio.

Use the nonoverlapping equal-shape model $I(t)=g(t)+\epsilon g(t-d)$, with $d>0$ and $0<\epsilon <1$. Let an independently controlled intensity transmission $a$ act on the delayed replica only. With the main pulse fixed, the third-order intensity autocorrelation has central, real-replica and third-order pulse correlation ghost lobe weights proportional toThe common factor is the same single-pulse correlation amplitude. Thus the real lobe has logarithmic slope 1 and its ghost slope 2 versus replica transmission, before normalization. Their exact ratio is $G/S=a\epsilon$.

An original example starts with $\epsilon =0.08$. Initially the weights are $(1.000512,0.08,0.0064)$. Halving the delayed replica's intensity gives $(1.000064,0.04,0.0016)$. The real side lobe halves while its ghost quarters. If each trace is divided by its own central peak, the predictions becomeThese normalized curves have the small additional central-peak dependence; treating them as exact pure powers would omit it.

In contrast, common attenuation in third-order pulse correlation sends $I\to kI$ for the complete two-pulse input. All three factors in $\int I^2(t)I(t+\tau )dt$ acquire the same scaling, so every unnormalized lobe scales as $k^3$. Halving the entire input therefore reduces every lobe eightfold while leaving the normalized shape unchanged. That experiment cannot demonstrate selective quadratic ghost scaling.

A suitable controlled test forms a delayed replica in a separate optical path, varies that path's transmission without changing its shape or delay, then recombines it with a monitored main pulse. Retain background-corrected, unsaturated readings and verify comparable detector response at both delays. The historical experiment's inserted-etalon postpulse supplied a physical delayed reference, not this particular attenuation sweep.^[ref]

Agreement with the scaling model supports the ghost interpretation under its assumptions. It does not exclude an additional real optical prepulse, unequal replica shapes, overlapping fields or an attenuation-dependent background. Those possibilities require further independent observations.

Colin Danson is a laser physicist whose research concerns ultra-high-power laser development and applications. His Oxford profile identifies this research portfolio and his work with the Oxford Centre for High Energy Density Science.^[ref]

Danson coauthored the1999 Central Laser Facility report on single-shot third-order intensity autocorrelation measurements.^[ref] Characterizing optical pulse temporal contrast matters because a large peak power does not describe all the energy arriving before or after the main pulse.

For an original rectangular-pulse example,20joules delivered in1picosecond givesA feature at one thousandth of that power still has20gigawatts. The ratio is small while the absolute power is large. Neither number alone specifies the feature's energy, which also depends on its duration.

Real pulse profiles require the integral $E=\int P(t)dt$ rather than assuming a rectangular shape. A diagnostic must therefore connect laser pulse energy, pulse duration and contrast with the actual temporal waveform. The example illustrates that measurement problem; its numbers are not historical Vulcan specifications.

Common attenuation in third-order pulse correlation multiplies the entire input intensity profile by one constant factor before the correlation is formed. For the ideal signalreplacing $I$ by $kI$ gives $A_{3,k}(\tau )=k^3{A}_3(\tau )$. All unnormalized lobes scale together, so division by the central peak removes this common factor. The primary correlator paper defines the three-intensity-factor integral and implements it through second-harmonic generation and sum-frequency generation.^[ref]

The factor $k$ is an intensity transmission. An optical field-amplitude factor $a$ instead gives intensity factor $|a{|}^2$ and correlation factor $|a{|}^6$. Confusing these conventions changes the predicted power law. A spectral filter that reshapes or chirps the pulse is also not necessarily a common scalar attenuation.

An original example takes an ideal central signal 100 and a side lobe 1 in consistent detector-signal units. Halving the full input intensity gives central 12.5 and side 0.125; their ratio remains 0.01. In contrast, attenuation scaling of a third-order correlation ghost changes only the delayed replica and predicts different real-lobe and ghost powers. A common attenuator cannot perform that selective test.

A fixed background changes what normalization appears to show. Add one unremoved signal unit to both locations. Before attenuation the observed side/center ratio is $2/101\approx 0.0198$. After common intensity transmission $k=1/2$, it is $1.125/13.5=0.0833$. The apparent side feature has grown relative to the peak even though the optical profile did not change. Detector baseline subtraction and a verified linear readout are therefore part of the experiment.

Shot-to-shot intensity variation supplies another consequence. Let the common intensity multiplier be $K=1+\eta$, where $\eta$ has zero mean. ThenFor equally likely multipliers 0.9 and 1.1, the mean intensity is unchanged but the mean cubic signal is $1.03$ times the unit-intensity signal. This cubic intensity-noise bias in pulse correlation shows why cubing a mean input is not equivalent to averaging third-order signals. The calculation assumes the pulse shape is unchanged on each shot.

A controlled check records pulse intensity alongside the correlator signal, varies a broadband common attenuator, and tests cubic scaling in an unsaturated interval after background correction. Failure of that test can indicate detector nonlinearity, pulse reshaping or an uncontrolled background; it does not uniquely identify a particular fault. The ideal law applies only while the nonlinear gate and optical overlap retain their stipulated response.

Cristina Hernandez-Gomez develops high-power laser technology and experimental facilities. Her Central Laser Facility biography describes work on the Vulcan Petawatt upgrade after joining in1997, followed by leadership of Vulcan and the High Power Lasers Division.^[ref]

She coauthored the1999 single-shot pulse-correlation report with the other Central Laser Facility researchers listed in its byline.^[ref] The diagnostic problem includes distinguishing a weak, extended signal from the main short pulse.

An original rectangular-waveform example separates peak contrast from energy contrast. Let the main pulse have power $P_0$ for1picosecond and a preceding feature have power $10^{-4}{P}_0$ for100picoseconds. Their energy ratio isThe preceding feature is at0.01% of the main peak, yet contains1% of the main pulse energy in this stipulated model. Reporting only its height would conceal its much longer duration.

This is why optical pulse temporal contrast and laser pulse energy need temporal definitions and a measured waveform. Facility development involves controlling the delivered pulse as well as its nominal peak: the diagnostic should establish what reaches the experiment, over the time interval relevant to that experiment.

Cubic intensity-noise bias in pulse correlation is the difference between the average of a third-order signal and the signal predicted from the average intensity when a fluctuating common intensity multiplier is cubed. Assume every shot has the same pulse shape, $I_j(t)=K_j{I}_0(t)$, with nonnegative scalar $K_j$. The ideal third-order integral givesThe three intensity factors follow from the second-harmonic generation and sum-frequency generation gate.^[ref]

Let $K=\mu +\delta$, with $\mathbb{E}[\delta ]=0$, variance ${\sigma }^2$ and third central moment $m_3$. Expanding the cube givesFor $\mu >0$, write the coefficient of variation as $v=\sigma /\mu$ and, when $\sigma >0$, skewness as $\gamma =m_3/{\sigma }^3$. The relative multiplier is $1+3v^2+\gamma v^3$. For zero variance the multiplier is one without needing to define skewness.

Equally likely $K=0.9$ and $1.1$ have mean one and give mean cubic multiplier $1.03$. Rare bright shots can matter much more: stipulated probabilities 0.9 and 0.1 at $K=0.8$ and $2.8$ also give mean one, butThe extra signal need not mean that the mean input intensity increased. These are original probability examples, not measurements reported in the correlator paper.

For unchanged pulse shape, common fluctuations multiply every delay by the same expected factor. Normalizing the averaged trace by its peak cancels that factor. Shape changes correlated with brightness destroy this simplification. For example, two equally frequent shot classes with central signals 1 and 8 and side/center ratios 0.01 and 0.10 give ratio of averaged signals $(0.01+0.8)/(1+8)=0.09$, whereas averaging the two ratios gives 0.055. They answer different questions.

A controlled measurement records a synchronized input-intensity monitor and detector background alongside each trace. Compare intensity-binned signals in an unsaturated range and test whether normalized shapes agree. Common attenuation in third-order pulse correlation predicts the cubic scaling only while the nonlinear gate, pulse shape and readout retain the assumed response. A mean power reading alone does not determine the required third moment.

An etalon delay inferred from free spectral range estimates an etalon round-trip delay from the frequency spacing of adjacent transmission resonances. When round-trip group delay varies little across that spacing,The measured quantity is resonance spacing, not the width of one resonance. The manufacturer tutorial explicitly distinguishes free spectral range from optical finesse, which compares that spacing with an individual peak's width.^[ref]

The relation follows from wave phase. Let $\varphi (\nu )$ be the total etalon round-trip phase, including relevant reflection phases. Adjacent resonances satisfy $\varphi ({\nu }_{m+1})-\varphi ({\nu }_m)=2\pi$. Defining ${\tau }_{\mathrm{rt}}(\nu )=(2\pi {)}^{-1}d\varphi /d\nu$ gives the exact interval relationThus the reciprocal spacing is the delay averaged over that frequency interval. Interpreting it as a local pulse-envelope delay requires sufficiently weak variation across the interval and relevant pulse spectrum.

For a uniform plate at normal incidence with negligible frequency-dependent reflection phase, $\varphi =4\pi \nu n(\nu )L/c$. Differentiation gives ${\tau }_{\mathrm{rt}}=2n_{g}L/c$, where $n_g=n+\nu dn/d\nu$ is the group refractive index. Treating $n$ as frequency independent recovers the manufacturer tutorial's elementary plate formula; using a phase index in a strongly dispersive material can instead give the wrong delay.

An original measured-spacing example takes $\Delta {\nu }_{\mathrm{FSR}}=25.0\mathrm{GHz}$. The inferred delay is $40.0\mathrm{ps}$. A spacing standard uncertainty of $0.2\mathrm{GHz}$ contributes approximately $0.32\mathrm{ps}$ to the delay uncertainty through $u(\tau )=u(\Delta \nu )/(\Delta \nu {)}^2$. With stipulated $n_g=1.50$ and negligible coating delay, the corresponding plate thickness is about $4.00\mathrm{mm}$.

In a spectral experiment, identify adjacent peaks of the same resonance family and calibrate the frequency coordinate. If the coordinate is wavelength, use $\Delta \nu =c|1/{\lambda }_1-1/{\lambda }_2|$ rather than treating a wavelength interval as a frequency interval. Beam angle, dispersion and reflection phase belong in the propagation model.

A time-domain check looks for replicas separated by the predicted extra round trip. The historical single-shot correlator experiment inserted a known-thickness BK7 etalon to calibrate delay and separately placed an uncoated etalon in the input to create a double-reflection postpulse.^[ref] The single-pass delay caused by inserting glass into an arm is distinct from the additional round-trip separation of successive transmitted replicas.

An etalon round-trip delay is the extra envelope travel time associated with one additional internal round trip in an etalon. At normal incidence, for a transparent plate of thickness $L$ and a sufficiently narrowband optical pulse,The relevant quantity is group refractive index, since the delay concerns the pulse envelope. The extra propagation phase is $\varphi (\omega )=2\omega n(\omega )L/c$; its derivative gives $d\varphi /d\omega ={\tau }_{\mathrm{rt}}$. This assumes negligible frequency-dependent reflection phase and treats dispersion locally.

For an original chosen plate with $L=3.0\mathrm{mm}$ and $n_g=1.60$, the delay is about $32.0\mathrm{ps}$. A pulse transmitted after one extra round trip appears that much later than the directly transmitted pulse. Another extra round trip adds the same delay in this local model, producing an echo sequence. The optical postpulse article separately describes the reflection losses that determine echo energies.

The frequency-domain connection gives an independent consistency check. Neighboring transmission resonances differ in round-trip phase by $2\pi$, so the local free spectral range satisfiesThe example predicts about $31.2\mathrm{GHz}$. This etalon delay inferred from free spectral range uses the same local dispersion assumptions; the spacing can vary across a broad spectrum.

The historical third-order correlator calibrated its delay coordinate by inserting a BK7 etalon of known thickness into an arm and measuring the trace shift.^[ref] A plate inserted into an arm also changes the direct propagation delay relative to the displaced air path; that single-pass insertion delay is distinct from the echo's extra internal round trip. Oblique incidence, dispersive coatings and broad-band pulse reshaping require a more complete propagation model.

Gaussian instrument broadening in third-order pulse correlation is the widening produced when an ideal third-order intensity autocorrelation is convolved with an additional Gaussian instrument response. It is a model of the measurement, not a universal correction for every correlator. The finite width of the optical nonlinear gate is already part of the ideal third-order trace and must not be counted again as an independent response.

Let a Gaussian pulse have intensity standard deviation ${\sigma }_I$. Its ideal trace is a Gaussian $C(\tau )=C_0\exp [-{\tau }^2/(2{\sigma }_3^2)]$ with ${\sigma }_3^2=3{\sigma }_I^2/2$. Take the additional normalized impulse response to beThe convolution $M=C\ast R$ isThe response preserves integrated trace area while reducing peak height and increasing variance. Since every Gaussian full width at half maximum has the same conversion factor from standard deviation,

For an original example, ${\tau }_I=80\mathrm{fs}$ and ${\tau }_R=60\mathrm{fs}$ give an ideal correlation width about $98.0\mathrm{fs}$ and a measured width about $114.9\mathrm{fs}$. The peak falls to approximately 0.853 of the ideal peak under this unit-area response. Dividing the measured width only by $\sqrt{3/2}$ would overestimate the pulse duration.

The correct Gaussian inversion is ${\tau }_I=\sqrt{(2/3)({\tau }_M^2-{\tau }_R^2)}$ when the radicand is positive. Its conditioning worsens near the response floor: the derivatives with respect to ${\tau }_M$ and ${\tau }_R$ are $(2/3){\tau }_M/{\tau }_I$ and $-(2/3){\tau }_R/{\tau }_I$. Small width errors can therefore cause large errors in a narrowly resolved pulse. Gaussian pulse-width deconvolution uncertainty must include uncertainty of the characterized response, with covariance if relevant.

The primary correlator used second-harmonic generation, sum-frequency generation, crossed beams and spatial imaging to encode delay. It did not assert that every additional broadening mechanism was a Gaussian convolution.^[ref] To apply this model, characterize the response under the same optical alignment, spectral acceptance and detector conditions. Non-Gaussian wings, clipping or unresolved side lobes require a different model; a smaller sampling step cannot by itself remove optical broadening.

A ghost-level consistency test for a delayed pulse replica checks whether measured side lobes agree with the simple model of one real delayed pulse and its nonlinear correlation ghost. It can reject that model; agreement alone does not prove that no real prepulse exists.

Use the separated, identical-shape replica model in third-order pulse correlation ghost. Let $a$ be the positive-delay lobe and $b$ the opposite-delay lobe, each divided by the central correlation peak. For delayed-to-main intensity ratio $0<\epsilon <1$,Thus the measured ratio gives $\epsilon =b/a$, and the first equation supplies an additional consistency condition. The frequently useful approximation $b\simeq a^2$ omits the small central-peak correction; exactly, $b/a^2=1+{\epsilon }^3$.

For an original pair of measured relative heights $a=0.080$ and $b=0.020$, their ratio would require $\epsilon =0.25$. That value predicts $a=0.25/1.015625\approx 0.246$, not $0.080$. These two readings therefore cannot both come from the stated single-replica model with equal instrument sensitivity at the two delays. A real earlier pulse is one possible explanation, but background, unequal pulse shapes or additional replicas can also invalidate the model.

An independent check varies the delayed replica while holding the main pulse fixed. Before central-peak normalization, halving $\epsilon$ halves the real side lobe and quarters its ghost. This attenuation scaling of a third-order correlation ghost probes the nonlinear origin rather than only one accidental height match. The existing parent diagram shows a direct numerical example without an actual pulse at the ghost's delay.

The historical correlator's inserted-etalon experiment provided a controlled physical postpulse and a known temporal direction.^[ref] The algebraic height test here is an independent ideal calculation. Use unsaturated, background-corrected readings with comparable response and lobe resolution; otherwise the measured ratios need not represent the model's intensities.

John Collier is a laser physicist associated with the Central Laser Facility. His Learned Society of Wales profile identifies laser physics as his specialist subject and gives the Rutherford Appleton Laboratory affiliation.^[ref]

The facility's historical account identifies Collier and I. Ross with the development of optical parametric chirped pulse amplification, including a1999 test system. In this technique a pump transfers energy through a nonlinear optical interaction to an already chirped signal; subsequent compression aims to recover a short, high-power pulse.^[ref]

Collier also coauthored the1999 report on a single-shot diagnostic for third-order intensity autocorrelation.^[ref] Amplification and characterization are complementary tasks: increasing the pulse energy does not by itself establish the temporal shape or the quality of compression.

For an original comparison, consider two hypothetical pulses with equal energy but rectangular durations of1 and10picoseconds. Their powers differ by a factor10. A measurement of the energy alone cannot distinguish them, while a time-sensitive diagnostic can constrain the missing duration and shape information. More generally, both amplification performance and the delivered temporal waveform matter in interpreting a high-power laser experiment.

An optical postpulse is a secondary optical pulse arriving after the chosen main pulse. Specify the delay, relative peak radiant power, duration and collected laser pulse energy: a relative peak alone does not determine the energy of a differently shaped replica. It is distinct from an extended optical pulse pedestal and from electronic ringing from a photodetector.

A transparent plate provides a concrete mechanism. At normal incidence, a transmitted replica that makes one additional internal round trip has approximate delaywhere $\ell$ is plate thickness and $n_g$ is its group refractive index. Neglect absorption, dispersion and overlap, and let each internal reflection retain power fraction $R$. The first delayed replica then has approximately $R^2$ times the directly transmitted energy. This etalon round-trip delay can create a physical optical echo.

For an original chosen $\ell =1.2\mathrm{mm}$, $n_g=1.5$ and $R=0.04$, the delay is about $12\mathrm{ps}$ and the energy ratio is $0.0016$. A $2.0\mathrm{mJ}$ main pulse would be followed by a $3.2\mu \mathrm{J}$ replica under that model. The historical single-shot correlator experiment deliberately inserted a BK7 etalon to generate and detect such a delayed pulse.^[ref]

Changing the plate thickness tests the delay prediction; changing optical alignment tests whether a stray reflection is involved. A third-order intensity autocorrelation can distinguish temporal sides under its model, while a symmetric second-order intensity autocorrelation cannot alone assign a feature to a prepulse or postpulse. An opposite-delay correlation feature may instead be a third-order pulse correlation ghost; it need not identify another physical pulse.

An optical prepulse is a secondary optical pulse arriving before the reference main pulse. The sign of its delay is part of the definition: choose the main pulse at $t=0$, so a prepulse has negative arrival time. A weak earlier pulse and an equally weak optical postpulse can affect a target differently when the material retains a response to the earlier illumination. Temporal order therefore carries information that one total laser pulse energy does not.

An original separated rectangular-profile example uses a main pulse of duration $100\mathrm{fs}$ and a prepulse of duration $10\mathrm{ps}$ centred $50\mathrm{ps}$ earlier. Let the prepulse's peak radiant power be $10^{-4}$ of the main peak. The optical pulse temporal contrast is $10^4$, but the energy ratio isThus the prepulse carries $1/101\approx 0.990\%$ of the two-pulse energy. If the main pulse carries $1\mathrm{mJ}$, the prepulse carries $10\mu \mathrm{J}$. These stipulated disjoint power profiles have no field-overlap term, and the main pulse's full width at half maximum does not reveal the earlier contribution.

A temporal diagnostic needs both sufficient delay coverage and calibrated dynamic range. The single-shot third-order correlator paper motivates its asymmetric response by contrasting it with the symmetric second-order trace, which loses the distinction between before and after.^[ref] Its delay axis still needs a known physical reference. A third-order pulse correlation ghost from a real postpulse can appear on the prepulse side; interpreting every small lobe as earlier light would be incorrect.

Repeating the measurement while varying an independently characterized optical delay helps assign temporal order. Record the detector baseline subtraction, coverage and response width along with the trace. A nondetection means a limit within those measurement conditions, rather than proof of zero early optical energy.

Optical pulse temporal contrast compares a pulse's peak with weaker light at another time. For a radiant power trace with peak $P_{\max }>0$, explicitly choose either the peak-to-local ratio $C(t)=P_{\max }/P(t)$ or its inverse $r(t)=P(t)/P_{\max }$. The first is at least one where the local power is positive; the second lies between zero and one. With a fixed transverse beam shape, the same ratio can be formed from light intensity. If the beam shape changes with time, specify the spatial sampling instead of identifying local intensity and total power.

A plateau at $10^{-3}$ of the peak has $r=10^{-3}$ and $C=1000$. Expressing this power ratio as optical contrast in decibels givesFor an amplitude ratio $A$ whose power ratio is $A^2$, the same definition gives $10{\log }_{10}(A^2)=20{\log }_{10}A$. Contrast also needs a delay or time range: an optical prepulse, an optical postpulse and a long optical pulse pedestal can have the same relative height but different positions and laser pulse energies. The pedestal energy fraction additionally requires integration over the contribution's duration.

A useful nonlinear measurement is third-order intensity autocorrelation. In an ideal instantaneous model, an input light intensity $I(t)$ generates a second-harmonic gate proportional to $I^2(t)$. Sum-frequency generation with a delayed fundamental replica gives a recorded signal proportional toThe instrument measures this cross-correlation, rather than simply reading $I(\tau )$. Finite gate width, phase matching, spatial overlap and detector response enter the interpretation.

The1998/99 instrument of J. Collier and colleagues used second-harmonic generation, a Wollaston prism and an optical biprism to form crossing fundamental and harmonic beams. A second crystal generated the sum-frequency signal, imaged through optical filters onto a charge-coupled device. Position encoded delay. An inserted glass etalon produced a delayed replica near $1.5\times 10^{-3}$ of the main signal, recorded on one side of the correlation peak. Figure4 combines window settings shifted with etalons.^[ref]

An independent ideal calculation explains both the asymmetry and a possible third-order pulse correlation ghost. Let $I(t)=g(t)+\epsilon g(t-d)$ represent two identical-shape pulses with well-separated supports and correlation lobes, with the second weaker and delayed. Squaring introduces no overlap term. If $A_g(\tau )=\int g^2(t)g(t+\tau )dt$, expansion givesThe delayed replica produces a lobe of relative order $\epsilon$, while the opposite-delay ghost has order ${\epsilon }^2$. For $\epsilon =0.01$, the scales are about $10^{-2}$ and $10^{-4}$ respectively. These values are predictions of the stated ideal model. Assigning every small correlation feature directly to a separate input pulse would miss this distinction.

A nondetection below the measurement floor establishes a bound under a detection model, rather than demonstrating that no light exists there. Finite dynamic range, detector baseline subtraction, delay coverage and resolution should accompany a reported contrast trace. The Sequoia manufacturer likewise distinguishes optical temporal resolution from the scan step; making that step smaller does not alone narrow the optical response.^[ref]

A pulse-contrast denominator acceptance rule decides which shots may be normalized by their measured central signal. For background-corrected side signal $S$ and central signal $C$, considerThe rule prevents normalization by an arbitrarily small accepted denominator. It also changes the estimand: averaging accepted ratios estimates $\mathbb{E}[R\mid A]$, a conditional expectation. An unconditional expected value is a different target and need not exist. The pulse-contrast ratio denominator noise explains why even the unfiltered ratio mean can differ from the ratio of mean signals.

The rule provides a precise stability guarantee when its moment assumptions hold. Interpret the indicator-weighted integrand below as zero outside $A$. If $p_A=\mathbb{P}(A)>0$ and $\mathbb{E}{S}^2<\infty$, thenThis proves finite ratio moments under a positive denominator cutoff. The accepted ratio has finite variance and finite first absolute moment, the latter by the Cauchy-Schwarz inequality. No independence of $S$ and $C$ was used. The bound concerns mathematical stability, not whether accepted shots represent the full pulse population.

For an original brightness-selection example, suppose half the shots have $(C,S)=(2,0.02)$ and half have $(10,0.02)$. The two ratios are $0.01$ and $0.002$, so their expected value over all shots is $0.006$. A cutoff $c_{\min }=5$ retains only the bright group, giving accepted conditional expectation $0.002$. Meanwhile the ratio of overall mean signals is $0.02/6\approx 0.003333$. These three quantities answer different questions. Physical dependence between brightness and temporal shape can make the choice consequential even without readout noise.

Selection bias from a noisy denominator threshold arises even at fixed true brightness. Consider another original model with constant $S=s>0$ and measured $C\sim N(\mu ,{\sigma }^2)$, where $\mu >0$. Write $z=(c_{\min }-\mu )/\sigma$, with $\varphi$ the standard normal density and $\Phi$ the standard normal distribution function. The accepted central measurements have a truncated normal distribution, withTo derive the second identity, standardize $C$ and use ${\varphi }^{\prime }(x)=-x\varphi (x)$, so ${\int }_z^{\infty }x\varphi (x)dx=\varphi (z)$. The density-to-tail quotient is an inverse Mills ratio.

Choose $\mu =5$, $\sigma =1$, $c_{\min }=6$ and $s=0.05$ in arbitrary signal units. The true-signal ratio is $s/\mu =0.01$. Only $15.87\%$ of measured shots pass, and the conditional expectation of their measured denominator is $6.5251$. Since $c\mapsto 1/c$ is a convex function for $c>0$, Jensen's inequality and the cutoff implyThe accepted ratios are necessarily biased downward relative to this model's true-signal ratio. The cutoff selects upward central-readout errors; increasing the accepted-shot count does not remove that conditioning.

In the historical single-shot correlator, nonlinear mixing encoded delays spatially and an imaging lens recorded the third-harmonic signal on a filtered CCD detector.^[ref] Selecting a trace by its central pixels conditions subsequent normalization on the same recorded data. Preserve the threshold, background treatment, number attempted and number accepted alongside the shot-resolved pulse-contrast estimator. For independent shots sharing the same acceptance probability, the accepted count has a binomial distribution, carrying information that disappears if only normalized accepted traces are retained.

Pulse-contrast ratio denominator noise is the distortion introduced when a noisy central pulse signal is used to normalize a side-delay signal. A shot-resolved pulse-contrast estimator uses $r=S/C$ after background correction, so errors in $C$ affect both precision and the distribution of the ratio.

For central values ${\mu }_S,{\mu }_C$ with ${\mu }_C>0$, a second-order local expansion givesThis approximation requires sufficiently small fluctuations and well-behaved ratio moments. It shows why an unbiased central-signal measurement need not give an unbiased normalized ratio. First-order uncertainty propagation alone describes local spread, not this mean shift; the GUM warns that nonlinear models can require higher-order terms.^[ref]

An original bounded example takes constant side signal $S=0.1$ and equally likely central readings $C=9$ or 11. The mean central reading is 10, so the ratio of means is 0.01. ButThe second-order prediction is 0.0101 because $\mathrm{Var}(C)=1$. Increasing the number of such shots estimates the shifted mean more precisely; it does not remove the nonlinear effect.

Correlation can instead make normalization exact. If $S=r_{0}C$ for every accepted shot with nonzero $C$, their ratio is always $r_0$, however variable the common multiplier. This is a different error structure from adding independent noise to the two channels. A shared additive background is also not generally equivalent to a shared multiplicative gain.

Near zero, even existence of a mean can fail. If $S=s_0>0$ is fixed and $C$ is uniform on $(0,1)$, every ratio is positive butThis mathematical example explains why “the denominator was never exactly zero” is insufficient. Detector noise is not asserted to follow this distribution.

A pulse-contrast denominator acceptance rule such as $C\ge C_{\min }$ changes the sampled population. Report the rule, rejected-shot count and dependence between brightness and pulse shape. The primary single-shot correlator's nonlinear gate and CCD readout motivate acquiring the relevant signals together, but do not guarantee that all ratios are statistically well conditioned.^[ref] Resolve the denominator and background experimentally before interpreting a long average of ratios.

Ric Allott has worked on photonics research and its industrial applications at the Central Laser Facility. The facility's2017 announcement of his visiting professorship describes a career spanning research, industrial work and technology transfer; it records his first CLF appointment in1995.^[ref]

Allott coauthored the1999 single-shot optical pulse-correlation report.^[ref] Such diagnostic work connects a physical response with a quantitatively interpretable measurement of optical pulse temporal contrast.

An original readout example illustrates why that interpretation matters. Normalize a reference peak to1. Suppose a weak feature gives0.0012 in the same units, including a separately measured additive baseline of0.0002. The corrected feature level isOmitting the baseline would overestimate this level by20%. The arithmetic is simple, but its physical validity depends on the baseline being representative and both measurements using a compatible scale.

Moving a diagnostic from a laboratory demonstration into applications therefore involves measurement calibration, reproducibility and a clear account of what its output represents. A large displayed ratio alone does not establish the measurement's accuracy. This example explains one aspect of the scientific-to-application connection without assigning unreported individual experimental tasks to Allott.

Selection bias from a noisy denominator threshold occurs when accepting measurements according to the same noisy central signal used for normalization changes the distribution of the reported ratios. For side signal $S$, measured central signal $C$, positive cutoff $c$ and acceptance event $A=\{C\ge c\}$, the sample mean over accepted shots targets $\mathbb{E}[S/C\mid A]$. It need not target the true-signal ratio or the unfiltered mean ratio. This is a consequence of conditioning, even when every physical pulse has the same true brightness.

For an original two-state readout model, let the true central signal be5, let measured $C$ equal4 or6 with equal probabilities, and let $S=0.02$ be constant. Both denominator errors have zero mean. The noiseless ratio is0.004, whereas the unfiltered mean measured ratio isA cutoff $c=5$ accepts only $C=6$, giving $\mathbb{E}[S/C\mid A]=1/300\approx 0.003333$ and acceptance probability1/2. This separates the unfiltered nonlinear ratio bias from the further change introduced by selection. Both denominators were already safely positive; the selection effect is not merely a division-by-zero problem.

The direction is not universal when numerator and denominator errors are related. Keep $C=4$ or6 but let $S=0.02+k(C-5)$. ThenFor $k=0.01$, the two observed side signals are0.01 and0.03, with ratios0.0025 and0.005. The same cutoff now gives0.005, above the same true-signal ratio0.004. For $k=0.004$, the observed signals are exactly proportional and every ratio is0.004, so selection leaves it unchanged. The noise dependence and measurement model determine whether thresholding shifts the accepted ratios upward, downward or not at all.

When a joint density and integrable accepted ratio exist, the dependence is explicit:The conditional numerator response cannot generally be replaced by its unconditional mean. Report the acceptance rule and accepted fraction alongside the shot-resolved pulse-contrast estimator; increasing the number of accepted shots reduces sampling scatter but does not undo the conditioning.

A physical motivation comes from a single-shot correlator that spatially encoded the nonlinear mixing signal and imaged it onto a filtered CCD detector.^[ref] If central pixels decide whether a trace is retained, later normalization uses data already conditioned by that choice. The two-state examples are original readout models, not measurements attributed to that instrument. Finite ratio moments under a positive denominator cutoff establish mathematical moment control under stated assumptions; they do not establish absence of selection bias.

A shot-resolved pulse-contrast estimator normalizes each individual optical shot before combining its contrast with other shots. For a specified side-delay signal $S_j$ and positive central signal $C_j$, both background-corrected, defineThe estimand is the mean of the individual shot ratios under the stated sampling and acceptance rule. It differs from normalizing an accumulated trace, ${\hat{R}}_{\mathrm{sum}}={\sum }_j{S}_j/{\sum }_j{C}_j$, which has third-order correlation brightness weighting.

An original three-shot example uses central signals $(1,2,8)$ and side signals $(0.01,0.04,0.80)$. The individual ratios are $(0.01,0.02,0.10)$, with equal-shot mean about 0.0433. Normalizing the sums gives $0.85/11\approx 0.0773$. The latter gives the brightest shot most of the influence. Both are valid summaries of the chosen data, but answer different questions.

Single-shot acquisition matters when the profile fluctuates. The primary third-order correlator maps delay into a spatially resolved signal and records it on a CCD, using second-harmonic generation and sum-frequency generation to form the gate.^[ref] A scanned trace assembled from different shots does not by itself provide all the $r_j$ for a common ensemble.

Normalization introduces pulse-contrast ratio denominator noise. First-order uncertainty propagation for one ratio $r=S/C$ givesShared baseline estimates can correlate the numerator and denominator, so their covariance cannot automatically be set to zero.^[ref]

For stipulated $S=0.1$, $C=10$, independent $u(S)=0.02$ and $u(C)=1$, the local approximation gives $r=0.01$ and $u(r)\approx 0.00224$. Near a zero or poorly resolved denominator the linear approximation and even a finite mean ratio may fail. A central-signal acceptance threshold changes the selected ensemble and must be reported.

For independent shots with a well-behaved ratio distribution, sampling uncertainty in the equal-shot mean can be estimated from the observed ratio variation. Consecutive shots may instead share drift, calibration errors or correlated pulse fluctuations. Keep the acquisition timing, background treatment, saturation checks and rejected-shot count with the estimate; averaging cannot repair an undefined denominator or an unspecified selection rule.

Third-order correlation brightness weighting is the way bright shots contribute disproportionately to an averaged nonlinear correlation trace. A normalized average trace is a weighted mean of shot-normalized shapes, with weights set by their central correlation signals rather than by how often the shots occur.

Let a nonnegative, background-corrected third-order trace on shot $j$ be $A_j(\tau )$. Define its central signal $B_j=A_j(0)>0$ and normalized shape $r_j(\tau )=A_j(\tau )/B_j$. Then the experimentally distinct averages areand the equal-shot mean $R_{\mathrm{shot}}(\tau )=N^{-1}{\sum }_j{r}_j(\tau )$. The nonnegative $w_j$ sum to one, but are equal only when the $B_j$ are equal. Here “central” denotes the specified zero-delay signal, not whichever noisy sample happens to be the largest.

For identical pulse shapes multiplied by intensity factors $K_j$, the ideal integral $A_j(\tau )=\int I_j^2(t)I_j(t+\tau )dt$ scales as $K_j^3$. The primary correlator implements this three-intensity-factor gate using second-harmonic generation followed by sum-frequency generation.^[ref] If shapes change, $B_j$ can depend on duration and overlap as well as brightness, so pulse energy alone need not determine the weight.

An original two-class example has equally frequent classes with $(B,r)=(1,0.01)$ and $(8,0.10)$ at one chosen delay. Their equal-shot mean ratio is $0.055$, but the normalized summed trace givesThe bright class supplies eight ninths of the central signal despite representing only half the shots. Neither result is an arithmetic mistake: they estimate different summaries of the stipulated ensemble.

For an ensemble with finite moments and $\mathbb{E}[B]>0$, the distinction can be expressed exactly through covariance:Unchanged normalized shape makes the covariance zero, even with large brightness fluctuations. Shape–brightness dependence can make it nonzero; cubic intensity-noise bias in pulse correlation alone does not prove a contrast change.

A shot-resolved pulse-contrast estimator needs synchronized intensity information and a stated selection rule. Thresholding weak shots or dividing by a noisy central signal can introduce further selection and ratio errors. Report whether normalization occurred before or after averaging, retain detector baseline subtraction, and test readout linearity. A smooth averaged trace is not a representative single shot unless those assumptions are justified.