Spectral phase interferometry for direct electric-field reconstruction

SPIDER is an interferometric ultrashort pulse measurement technique in the frequency domain based on spectral shearing interferometry. Spectral shearing interferometry is similar in concept to lateral shearing interferometry, except the shearing is performed in the frequency domain. The spectral shear is typically generated by sum-frequency mixing the test pulse with two different quasi-monochromatic frequencies (usually derived by chirping a copy of the pulse itself), although it can also be achieved by spectral filtering or even with linear electro-optic modulators for picosecond pulses. The interference between the two upconverted pulses allows the spectral phase at one frequency to be referenced to the spectral phase at a different frequency, separated by the spectral shear - the difference in frequency of the two monochromatic beams. In order to extract the phase information, a carrier fringe pattern is introduced, typically by delaying the two spectrally sheared copies with respect to one another.

Flow chart describing the SPIDER reconstruction algorithm

The intensity of the interference pattern from two time-delayed spectrally sheared pulses can be written as

${\displaystyle {\begin{aligned}S(\omega )&=|E(\omega )+E(\omega -\Omega )e^{i\omega \tau }|^{2}\\&=I(\omega )+I(\omega -\Omega )+2{\sqrt {I(\omega )I(\omega -\Omega )}}\cos[\phi (\omega )-\phi (\omega -\Omega )-\omega \tau ]\end{aligned}}}$,

where ${\displaystyle E(\omega )}$ is the analytic signal representing the unknown (upconverted) field being measured, ${\displaystyle \Omega }$ is the spectral shear, ${\displaystyle \tau }$ is the time delay, ${\displaystyle I(\omega )=|E(\omega )|^{2}}$ is the spectral intensity and ${\displaystyle \phi (\omega )}$ is the spectral phase. For a sufficiently large delay (from 10 to 1000 times the Fourier transform limited [FTL] pulse duration), the interference of the two time-delayed fields results in a cosine modulation with a nominal spacing of ${\displaystyle \delta \omega \sim 2\pi /\tau }$; and any dispersion of the pulse results in minor deviations in the nominal fringe spacing. Effectively it is these deviations in the nominal phase spacing that yield the dispersion of the test pulse .

The unknown spectral phase of the pulse can be extracted using a simple, direct algebraic algorithm first described by Takeda.[1] The first step involves Fourier transforming the interferogram into the pseudo time domain:

${\displaystyle {\begin{aligned}{\widetilde {S}}({\widetilde {t}})&={\mathfrak {F}}[S(\omega )]\\&={\widetilde {E}}^{dc}({\widetilde {t}})+{\widetilde {E}}^{ac}({\widetilde {t}}-\tau )+{\widetilde {E}}^{-ac}({\widetilde {t}}+\tau )\end{aligned}}}$,

where ${\displaystyle {\widetilde {E}}^{dc}({\widetilde {t}})={\mathfrak {F}}[I(\omega )+I(\omega -\Omega )]}$ is a 'direct current' (dc) term centred at ${\displaystyle {\widetilde {t}}}$ with a width inversely proportional to the spectral bandwidth, and ${\displaystyle {\widetilde {E}}^{\pm ac}({\widetilde {t}}\mp \tau )={\mathfrak {F}}\{{\sqrt {I(\omega )I(\omega -\Omega )}}e^{\pm i[\phi (\omega )-\phi (\omega -\Omega )]}e^{\pm i\omega \tau }\}}$ are two 'alternating current' (ac) sidebands resulting from the interference of the two fields. The dc term contains information about the spectral intensity only, whereas the ac sidebands contain information about the spectral intensity and phase of the pulse (since the ac sidebands are Hermitian conjugates of each other, they contain the same information).

One of the ac sidebands is filtered out and inverse Fourier transformed back into the frequency domain, where the interferometric spectral phase can be extracted:

${\displaystyle {\begin{aligned}D(\omega ,\Omega )&={\mathfrak {F}}^{-1}[{\widetilde {E}}^{ac}({\widetilde {t}}-\tau )]\\&={\sqrt {I(\omega )I(\omega -\Omega )}}e^{i[\phi (\omega )-\phi (\omega -\Omega )]}e^{-i\omega \tau }\end{aligned}}}$.

The final exponential term, resulting from the delay between the two interfering fields, can be obtained and removed from a calibration trace, which is achieved by interfering two unsheared pulses with the same time delay (this is typically performed by measuring the interference pattern of the two fundamental pulses which have the same time-delay as the upconverted pulses). This enables the SPIDER phase to be extracted simply by taking the argument of the calibrated interferometric term:

${\displaystyle {\begin{aligned}\theta (\omega )&=\angle [D_{\text{cal}}(\omega )D^{\ast }(\omega ,\Omega )]\\&=\phi (\omega -\Omega )-\phi (\omega )\end{aligned}}}$.

There are several methods to reconstruct the spectral phase from the SPIDER phase, the simplest, most intuitive and commonly used method is to note that the above equation looks similar to a finite difference of the spectral phase (for small shears) and thus can be integrated using the trapezium rule:

${\displaystyle \phi (\omega _{N}-\Omega /2)\approxeq -\sum _{n=0}^{N}{\frac {\omega _{n+1}-\omega _{n}}{2\Omega }}[\theta (\omega _{n})+\theta (\omega _{n+1})]}$.

This method is exact for reconstructing group delay dispersion (GDD) and third order dispersion (TOD); the accuracy for higher order dispersion depends on the shear: smaller shear results in higher accuracy.

An alternative method us via concatenation of the SPIDER phase:

${\displaystyle {\begin{aligned}\phi (\omega _{0}+N|\Omega |)&={\begin{cases}-\sum _{n=1}^{N}\theta (\omega _{0}+n\Omega )&{\text{if}}\,\Omega >0\\\sum _{n=0}^{N-1}\theta (\omega _{0}+n|\Omega |)&{\text{if}}\,\Omega <0\end{cases}}\end{aligned}}}$

for integer ${\displaystyle N}$ and concatenation grid ${\displaystyle \{\omega _{N}\}=\{\omega _{0}+N|\Omega |\}}$. Note that in the absence of any noise, this would provide an exact reproduction of the spectral phase at the sampled frequencies. However, if ${\displaystyle D(\omega )}$ falls to a sufficiently low value at some point on the concatenation grid, then the extracted phase difference at that point is undefined and the relative phase between adjacent spectral points is lost.

The spectral intensity can be found via a quadratic equation using the intensity of the dc and ac terms (filtered independently via a similar method above) or more commonly from an independent measurement (typically the intensity of the dc term from the calibration trace), since this provides the best signal to noise and no distortion from the upconversion process (e.g. spectral filtering from the phase matching function of a 'thick' crystal).

Spatially encoded arrangement for SPIDER (SEA-SPIDER) is a variant of SPIDER.[2] [3] [4] [5] The spectral phase of an ultrashort laser pulse is encoded into a spatial fringe pattern rather than a spectral fringe pattern.

Other techniques are frequency-resolved optical gating, streak camera with picosecond response times, and multiphoton intrapulse interference phase scan (MIIPS), a method to characterize and manipulate the ultrashort pulse.

Micro-SPIDER is an implementation of SPIDER in which the spectral shear required for a SPIDER measurement is generated in a thick nonlinear crystal with a carefully engineered phase-matching function.[6] [7]

• Spectroscopy

• US patent 6611336, Ian A. Walmsley & Chris Iaconis, "Pulse measurement using frequency shifting techniques", issued 2003-8-26
• Iaconis, C; Walmsley, I. A. (1999), "Self-Referencing Spectral Interferometry for Measuring Ultrashort Optical Pulses", IEEE J. Quantum Electron., 35 (4): 501–509, Bibcode:1999IJQE...35..501I, doi:10.1109/3.753654, S2CID 55097406
• Iaconis, C; Walmsley, I. A. (1998), "Spectral Phase Interferometry for Direct Electric-Field Reconstruction of Ultrashort Optical Pulses", Opt. Lett., 23 (10): 792–794, Bibcode:1998OptL...23..792I, doi:10.1364/OL.23.000792, PMID 18087344
• Walmsley, I. A.; Wong, V. (1996), "Characterization of the Electric Field of Ultrashort Optical Pulses", J. Opt. Soc. Am. B, 13 (11): 2453–2463, Bibcode:1996JOSAB..13.2453W, doi:10.1364/JOSAB.13.002453

• new SPIDER page includes links to example code