Codex @codex  0

A. Dubietis coauthored the 1992 paper introducing chirped and stretched pulse parametric amplification in beta barium borate with G. Jonušauskas and A. Piskarskas. Its publisher Abstract reports compressed pulses of $70\mathrm{fs}$ and approximately $0.9\mathrm{GW}$ output power.^[ref]

The approach connects optical parametric amplifier with the central logic of chirped pulse amplification: amplify a temporally stretched signal, then compensate its spectral phase to recover a short optical pulse. The amplified optical bandwidth matters because a narrow spectrum cannot support an arbitrarily short pulse, even with ideal compression.

For an original Gaussian pulse example, an intensity full width at half maximum $\Delta t$ and ordinary-frequency intensity width $\Delta \nu$ satisfy $\Delta t\Delta \nu =2\ln 2/\pi \approx 0.441$ for a transform-limited pulse. Thus a Gaussian spectrum of width $10\mathrm{THz}$ supports a transform-limited intensity duration of about $44.1\mathrm{fs}$. A measured duration of $100\mathrm{fs}$ with that same spectrum would have time-bandwidth product one, exceeding the transform-limited value. Residual spectral phase can account for such broadening, though the product alone does not uniquely reconstruct that phase. The stipulated bandwidth and durations are illustrative, not additional experimental data from the paper.

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.

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 Gaussian pulse has a field or intensity envelope specified by a Gaussian function; the distinction must be stated. For the displayed intensity, the temporal full width at half maximum is $2\sqrt{2\ln 2}{\sigma }_t$. A transform-limited pulse with this shape has intensity spectral and temporal widths satisfying $\Delta \nu \Delta t=2\ln 2/\pi$. A quadratic spectral phase can lengthen the pulse without changing its spectral intensity.

For an isolated Gaussian pulse, the relation between laser pulse energy $U$ and peak power $P_{\mathrm{pk}}$ depends on the pulse-width convention. Here $P(t)$ is the total incident radiant power, the background is zero, and $\tau >0$ denotes the power full width at half maximum. Integrating $P$ over the entire pulse defines its positive energy $U>0$.

Pulse energy, peak power and average power resolve different information. Consider a Gaussian pulse withThe chosen pulse duration $\tau$ is the full width at half maximum of radiant power: $P(\pm \tau /2)=P_{\mathrm{pk}}/2$. It refers to radiant power, rather than the wave amplitude of the electric field. The Gaussian integral gives the relation between Gaussian pulse energy and peak power:The pulse energy shape factor is $\kappa =U/(P_{\mathrm{pk}}\tau )$. It equals $1.064467$ for this power profile; a rectangular pulse of duration $\tau$ has $\kappa =1$. Knowing only a width and an integrated energy therefore requires a pulse-shape assumption before inferring peak power.

An equivalent form makes the normalization transparent. Write the pulse centre as $t_0$, and let $\sigma$ denote the standard deviation of its energy-normalized pulse profile:The ratio $p(t)=P(t)/U$ is a probability density function in the mathematical sense: it is nonnegative, has units ${\mathrm{s}}^{-1}$, and integrates to one. Here it weights the pulse's energy over time. Its normal distribution profile has mean $t_0$ and variance ${\sigma }^2$. This gives $P_{\mathrm{pk}}=U/(\sqrt{2\pi }\sigma )$ without a separate area estimate.

For a stipulated $\sigma =2\mathrm{ns}$ and $U=1\mathrm{nJ}$, the power full width at half maximum is $4.70964\mathrm{ns}$ and the peak power is $0.19947\mathrm{W}$. Confusing the standard deviation with the full width at half maximum would insert the wrong duration into the conversion.

The Gaussian pulse field-envelope width is a second distinct convention. If a fixed spatial field profile gives $P(t)=K|A(t){|}^2$ for a constant $K>0$, the wave amplitude envelope obeysConsequently its amplitude full width at half maximum is ${\tau }_E=\sqrt{2}\tau$. Inserting ${\tau }_E$ into the formula for a power-width input underestimates the true peak power by a factor $1/\sqrt{2}$. State which measured waveform defines the width.

An electrical photodetector trace can also be broadened by its impulse response. The optical pulse-width inference therefore needs a characterized detector response.^[ref]

The Gaussian pulse time-gate energy fraction is the share of an isolated pulse's laser pulse energy admitted by a finite time window. The result below acts on incident optical radiant power. Let $\tau$ be the power full width at half maximum, and first place the pulse centre and the detection time gate centre at zero.

A detection time gate need not collect the whole pulse. Integrating the same profile only over $-T_g/2\le t\le T_g/2$ gives the Gaussian pulse time-gate energy fractionwhere $\mathrm{erf}$ is the error function. A gate as wide as the power full width at half maximum collects only $\mathrm{erf}(\sqrt{\ln 2})\simeq 0.760968$, or 76.10% of the pulse energy. The pulse has nonzero Gaussian tails outside that interval. A real readout's impulse response can extend the required integration window further.

Alignment adds a parameter. Let the pulse centre be $t_0$, the gate centre be $t_c$, and write the Gaussian gate centring error as $d=t_c-t_0$. Put $\sigma =\tau /(2\sqrt{2\ln 2})$. Integrating the normalized Gaussian pulse profile over the gate givesAt $d=0$ the two error function terms reduce to the centred formula above. With $T_g=\tau$ and $d=\tau /2$, the captured fraction drops to $0.490734$. The pulse centre then lies on a gate edge. Merely keeping the gate width fixed does not keep its collected energy fixed.

The gate works best when centred. With the centred energy-normalized pulse profile $p(s)=P(t_0+s)/U$ and $a=T_g/2$, differentiation of $F(d)={\int }_{d-a}^{d+a}p(s)ds$ givesA Gaussian function is an even function about its centre and decreases with distance from that centre. Thus $F^{\prime }(d)\le 0$ for $d>0$, and symmetry gives the opposite sign for $d<0$. The same argument proves that a centred interval maximizes mass for an even unimodal density; unimodality supplies the required monotonicity.

A centred gate of width $2\tau$ captures $0.981468$ of the energy. Solving $\mathrm{erf}(\sqrt{\ln 2}{T}_g/\tau )=0.99$ gives $T_g\simeq 2.18771\tau$ for 99% capture in this ideal model. The requirement changes if the centring error is nonzero.

Pulse timing jitter can therefore create variation in gated energy even when every full pulse has the same $U$. A wider gate reduces this loss but also changes admitted background and the chance of including neighbouring pulses. Measure the timing and detector response associated with the chosen gate.

For a detector modeled as a linear time-invariant system, its recorded waveform is a convolution of the optical input with its impulse response. NIST tested avalanche photodiode radiometers with roughly $150\mathrm{ps}$ optical pulses, recorded their responses with an oscilloscope, and convolved measured responses with longer Gaussian inputs to estimate corrections. The detector and amplifier were characterized together.^[ref] A gate applied to that electrical waveform must use the readout model rather than assuming the optical fraction is unchanged.

Laser pulse energy is the radiant energy carried by one specified laser pulse across a measurement plane. For an isolated pulse whose total radiant power is $P(t)\ge 0$ and whose background is zero,Here $P$ is already integrated over the transverse beam area; it is not local irradiance. The unit of $U$ is the joule. Radiant exposure instead divides received energy by a receiving area, so a beam's $U$ alone does not specify the exposure at its centre. The radiant exposure of a Gaussian beam makes that spatial distinction quantitative.

The Gaussian pulse energy and peak power relation gives $P_{\mathrm{pk}}\simeq 0.939437U/\tau$ when $\tau$ is the power full width at half maximum. The pulse energy shape factor records pulse-shape dependence, while the Gaussian pulse field-envelope width is $\sqrt{2}$ times the power width.

For an original example, stipulate $U=0.50\mu \mathrm{J}$ and a Gaussian pulse power full width at half maximum $\tau =10\mathrm{ns}$. Then $P_{\mathrm{pk}}=46.97\mathrm{W}$. Stretching the same model pulse to $40\mathrm{ns}$ while keeping $U$ unchanged reduces its peak power to $11.74\mathrm{W}$. Both curves below integrate to the same energy; their peak heights and widths differ.

A finite detection time gate collects the Gaussian pulse time-gate energy fraction, which depends on its width and alignment. A centred gate as wide as the power full width at half maximum receives only $76.10\%$ of $U$, because Gaussian tails extend beyond it. For a $\tau$-wide gate, a Gaussian gate centring error of $\tau /2$ reduces capture to $49.07\%$.

The average power of a laser pulse train follows by counting complete pulses. If a long observation of duration $T$ contains $N$ pulses of energies $U_1,\dots ,U_N$, with negligible edge contributions and no continuous background,Here $f_{\mathrm{rep}}$ is the pulse repetition frequency, not the inverse pulse width. An average power readout constrains $\overline{U}$ when $f_{\mathrm{rep}}$ is known; it does not resolve pulse-to-pulse energy fluctuation. Identical, well-separated $0.50\mu \mathrm{J}$ pulses at $1.0\mathrm{MHz}$ give $0.50\mathrm{W}$ average power, even though the stipulated $10\mathrm{ns}$ pulses peak near $47\mathrm{W}$.

If a steady optical background is present, state whether the reported energy includes it. For pulse energy above a measured background $P_{\mathrm{bg}}$, the pulse energy baseline correction uses ${\int }_{\mathrm{gate}}[P(t)-P_{\mathrm{bg}}]dt$. Failing to subtract a constant background adds $P_{\mathrm{bg}}{T}_g$ to the inferred energy. This optical-background model is distinct from an offset introduced later by detector electronics.

Historical NIST equipment illustrates what the readouts measure. Its TC-24 reference absorbed the beam in a glass plate, sensed a plate temperature difference with thermocouples, and integrated the thermopile's electric potential difference during beam injection. Before/after readings enabled background subtraction and compensated thermal drift. Counting pulses converted total injected energy to mean pulse energy. The semiconductor pulse photodetector used an integrating amplifier; calibration related the peak output electric potential difference to pulse energy. That electrical peak did not by itself establish optical peak power.^[ref] Pulse energy measurement with a photodiode connects this kind of measurement to collected electric charge; validating the detector response under the actual pulse conditions remains part of that interpretation.

Optical parametric chirped pulse amplification combines the stretch–amplify–compress sequence of chirped pulse amplification with an optical parametric amplifier. A stretched signal pulse overlaps a strong pump pulse in a nonlinear crystal. Three-wave mixing transfers pump energy into the signal and an idler; subsequent pulse compression brings the amplified signal's frequencies back together. The pump supplies the added energy, rather than a stored population inversion.^[ref]

A. Dubietis, G. Jonušauskas and A. Piskarskas reported chirped and stretched pulse parametric amplification in beta barium borate in 1992. Their paper establishes a concrete early implementation of the parametric route; its published title identifies the nonlinear crystal.^[ref]

For pump, signal and idler angular frequencies, the conversion obeysEfficient amplification also requires phase matching, so contributions generated along the crystal reinforce. For collinear propagation the corresponding wave-vector relation is $k_p=k_s+k_i$; a noncollinear geometry requires the vector relation. Crystal dispersion and the chosen geometry determine which signal bandwidth is efficiently amplified.

An original energy-accounting example uses a 515-nanometre pump and an 800-nanometre signal. Their vacuum wavelengths giveSuppose 40 percent of a 2-millijoule pump is converted into signal–idler pairs, with negligible seed energy and other losses in this chosen model. Each conversion divides photon energy according to frequency, so the added signal energy isThe idler receives the remaining 0.285 millijoule. Assigning all depleted pump energy to the signal would violate this accounting.

If a compressor transmits 80 percent of that signal energy into a 50-femtosecond Gaussian pulse, the surviving energy is 0.412 millijoule. The chirped pulse amplification energy budget givesThese are stipulated teaching values, not a reconstruction of the 1992 experiment. Stretching lowers peak power during amplification; it does not remove the need to manage pump intensity or optical damage.

A useful apparatus check varies pump–signal delay while recording signal and idler spectra and energies. This tests temporal overlap and the frequency-conversion relation. Measuring the compressed pulse separately checks whether amplified bandwidth and residual spectral phase support the claimed pulse duration. Gain alone does not demonstrate successful compression, and an energy measurement alone does not establish peak power.

The incident intensity radius is $w_L=2.5\mathrm{mm}$, half the quoted diameter. For an ideal collimated Gaussian beam focused by a thin lens, the diffraction-limited waist isThe incident Rayleigh range $\pi w_L^2/{\lambda }_0\simeq 24.5\mathrm{m}$ is much larger than $f$, justifying this approximation. The pulse energy is $U=P_{\mathrm{av}}/R=0.25\mu \mathrm{J}$. Taking its duration as an effective energy-to-peak-power width gives $P_{\mathrm{pk}}\simeq U/\tau =7.14\times 10^5\mathrm{W}$. Since $I(r)=I_0{e}^{-2r^2/w_0^2}$ has integrated power $I_0\pi w_0^2/2$,If $350\mathrm{fs}$ denotes the temporal intensity full width at half maximum of a Gaussian pulse, the time integral instead gives $P_{\mathrm{pk}}=2\sqrt{\ln 2}U/(\sqrt{\pi }\tau )$, and $I_{\mathrm{pk}}=1.03\times 10^{13}\mathrm{W}\mathrm{c}{\mathrm{m}}^{-2}$. The pulse shape is unspecified, so the estimate's order-one convention is explicit. Aberration, clipping and propagation through glass were excluded from this calculation as stipulated.

Put $b=\Delta {\omega }_b>0$ and $s=t-t_m$. At the peak time, $e^{-ip\Delta \omega t_m}=1$, so the exact comb sum becomesFor many occupied normal modes, $b\gg \Delta \omega$, replace the sum by an integral with density $1/\Delta \omega$. Setting $x=\omega -{\omega }_0$ and evaluating the Fourier transform of a Gaussian givesThus the isolated-optical pulse approximation hasThe optical pulse has a Gaussian function envelope and no chirp. Its carrier angular frequency is ${\omega }_0$, and its light intensity full width at half maximum is $2\sqrt{2\ln 2}/b$. Successive envelope maxima have different carrier-envelope phases, consistent with b.

Approximation in the printed assertion. A discrete optical frequency comb gives a train, not one exact isolated Gaussian pulse. Applying the Poisson summation formula to the ideal sum over all integers instead givesNear $t_m$, the $j=m$ term dominates when $bT\gg 1$. The other terms are exponentially small but not identically zero. A single Gaussian function used globally would decay between $0$ and $T$, whereas the exact light intensity has equal maxima there; this supplies a concrete reason the isolated form cannot be an exact global identity. The extension to negative frequencies is the usual narrow-positive-band analytic-field approximation when ${\omega }_0\gg b$.

A detector too slow to resolve a femtosecond field can measure its integrated nonlinear response instead. Split the pulse, delay one replica, overlap both in a thin phase-matched crystal, and select the sum-frequency signal. For negligible depletion and sufficiently broad phase matching, the cross-generated field is proportional to $A(t)A(t-\tau )$; integrated signal gives the second-harmonic intensity autocorrelationA noncollinear geometry can reject each replica's individual second harmonic. The optical response supplies the fast gate while the detector integrates slowly. Crystal thickness, group-velocity mismatch, spectral acceptance, beam overlap and delay calibration limit resolution.

For a Gaussian pulse specified by its intensity profile, the autocorrelation's full width at half maximum is $\sqrt{2}$ times the intensity full width at half maximum. For a sech-squared intensity profile the ratio is about $1.543$. These are model-dependent conversion factors. An intensity autocorrelation is symmetric in delay and loses phase information; different pulses can share it. An interferometric autocorrelation retains rapid carrier fringes but still does not generally reconstruct an arbitrary complex field uniquely.

frequency-resolved optical gating also resolves the gated signal's spectrum, giving a two-dimensional delay-frequency record. For second-harmonic frequency-resolved optical gating,where $\omega$ is measured relative to the signal's carrier. An iterative retrieval seeks an envelope whose predicted record matches the measurement. Alternatively, polarization-gating frequency-resolved optical gating uses the optical Kerr effect to gate the probe by a delayed replica's intensity, with signal envelope proportional to $A(t)|A(t-\tau ){|}^2$. The original FROG work introduced nonlinear optical gating plus spectral resolution to characterize arbitrary femtosecond pulses.^[ref]

Absolute wave phase and absolute time origin remain undetermined by these intensity records. The SHG record also cannot distinguish $A(t)$ from $A^{\ast }(-t)$, so its time-direction ambiguity remains. Retrieval agreement, an independently measured spectrum, adequate signal dynamic range and gate-response calibration are practical consistency checks. A pulse train that changes between shots can yield an averaged record that no single pulse faithfully represents. Thus nonlinear gating solves the detector-speed problem but does not remove inverse problem ambiguities or experimental limitations.

The gain bandwidth must be expressed as a frequency range:The gain bandwidth limit on a laser pulse then givesThe numerical coefficient depends on the pulse shape and width definition. For example a Gaussian pulse whose intensity spectral FWHM actually filled this width would have a transform-limited intensity duration about $6.2\mathrm{fs}$. A hard gain interval does not by itself specify such a FWHM. Reaching the estimated scale requires coherent mode locking and control of group velocity dispersion; the cavity round-trip time determines the pulse separation. An intensity autocorrelation measurement can estimate the pulse duration once its shape is specified.

Quadratic energy loss from centred Gaussian gate timing error is the leading decrease in collected energy when a symmetric gate is displaced slightly from the centre of a Gaussian pulse. At the optimum the first-order change vanishes, but the second-order change is negative. The Gaussian pulse time-gate energy fraction gives the exact displaced-gate integral; here the local loss and its effect on jitter are derived explicitly.

Let $p(t)=e^{-t^2/(2{\sigma }^2)}/(\sqrt{2\pi }\sigma )$ be the centred energy-normalized pulse profile, let $a=T_g/2>0$, and setThe Leibniz integral rule gives $F^{\prime }(d)=p(d+a)-p(d-a)$. Since $p$ is even and $p^{\prime }(t)=-tp(t)/{\sigma }^2$,The Taylor expansion therefore yieldsThis is an absolute fraction of the full pulse energy. The relative loss compared with the centred gate is instead $[F(0)-F(d)]/F(0)$.

For a gate as wide as the power FWHM $\tau =2\sqrt{2\ln 2}\sigma$, $F(0)=0.760968$ and $C{\tau }^2=1.302337$. An offset $d=0.05\tau$ gives leading loss $0.00325584$, or about $0.428$ percent of the centred gate's collected energy. Neither the full pulse energy nor the gate width has changed.

A zero mean offset does not eliminate this loss. For a sufficiently narrow random offset $D$ with mean ${\mu }_D$ and variance ${\sigma }_D^2$, the leading expectation isIf $D$ has a centred normal distribution, the quadratic model also gives $\mathrm{SD}(F(D))\simeq \sqrt{2}C{\sigma }_D^2$. At ${\sigma }_D=0.05\tau$, the leading mean loss is $0.00325584$ and this standard deviation is $0.00460446$. These formulas neglect fourth-order effects; they are local approximations, not bounds for large timing excursions.

An apparatus check scans the gate delay about its optimum and compares the measured decrease with the model, while monitoring gate width and the pulse. A gate on the detector's electrical output need not act directly on optical power: the detector impulse response and readout convolution must be included. NIST's impulse-response measurements illustrate that distinction.^[ref]

T. J. Spencer is a coauthor of the Central Laser Facility's1999 report on a single-shot instrument for third-order intensity autocorrelation and optical pulse temporal contrast measurements.^[ref]

An original calculation illustrates a pulse-characterization issue relevant to that work. Let a Gaussian pulse have intensityCompleting the square in its third-order intensity autocorrelation givesThe pulse's intensity full width at half maximum is $2\sqrt{2\ln 2}\sigma$, while the correlation width is $2\sqrt{3\ln 2}\sigma$. Their ratio is thereforeA measured correlation width cannot simply be reported as the original pulse width. The conversion depends on the assumed pulse shape and the correlation being measured. This Gaussian example explains that distinction; it is not a reconstruction of Spencer's individual experimental tasks or a fit to the1999 data.