Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Single-shot high-resolution characterization of optical pulses by spectral phase diversity

Open Access Open Access

Abstract

The concept of spectral phase diversity is proposed and applied to the temporal characterization of optical pulses. The experimental trace is composed of the measured power of a plurality of ancillary optical pulses derived from the pulse under test by adding known amounts of chromatic dispersion. The spectral phase of the pulse under test is retrieved by minimizing the error between the experimental trace and a trace calculated using the known optical spectrum and diagnostic parameters. An assembly composed of splitters and dispersive delay fibers has been used to generate 64 ancillary pulses whose instantaneous power can be detected in a single shot with a high-bandwidth photodiode and oscilloscope. The diagnostic is experimentally shown to accurately characterize pulses from a chirped-pulse–amplification system when its stretcher is detuned from the position for optimal recompression. Pulse-shape reconstruction for pulses shorter than the photodetection impulse response has been demonstrated. Various investigations of the performance with respect to the number of ancillary pulses and the range of chromatic dispersion generated in the diagnostic are presented.

© 2015 Optical Society of America

1. Introduction

Short optical pulses are used in many areas of optics and physics. Measuring their pulse shape became a significant technical challenge after the demonstration of mode-locking because quasi-instantaneous optical nonlinearities can yield much-shorter pulses than can be directly measured via electronic processes. A large variety of strategies and diagnostics have been demonstrated to circumvent the relatively low speed of photodetectors and sampling systems [1], but single-shot characterization of picosecond pulses remains difficult. Techniques that rely on the spatiotemporal duality are an attractive means of measuring the instantaneous power of optical pulses. Because the spectrotemporal variables (ω, t) are Fourier conjugates in a way identical to the one-dimensional (1-D) wave-vector/space coordinates (k, x), temporal imaging systems can be designed to temporally magnify the optical waveform under test and make it directly measurable via photodetection [2,3]. These imaging systems combine spectral and temporal quadratic phase modulations, the former being generally provided by dispersive materials or assemblies including angular dispersion, and the latter by high-bandwidth electro-optic phase modulators or nonlinear wave mixing with ancillary chirped pulses [1]. A photodiode followed by an oscilloscope can characterize the waveform under test without significant loss of details from the convolution by the photodetection impulse response when there is sufficient temporal magnification. Time-to-frequency conversion is an alternative to temporal magnification, where the instantaneous power of the waveform under test is mapped onto the optical spectrum of an optical pulse that can be measured with a spectrometer.

Time magnification and time-to-frequency conversion are technically elegant solutions that are well adapted to the telecommunication environment [4–6] but they have also been applied to the characterization of isolated events with low duty cycle [7,8]. High-energy laser systems are complex large-scale optical systems that operate at a low rate (~1-h shot cycle) and generate short (1- to 100-ps) optical pulses with energy in the kilojoule range [9]. There is significant operational complexity for accurately synchronizing a subsystem that induces quadratic temporal phase modulation to the pulse under test at the end of the laser system. Such lasers require temporal diagnostics for safe operation (e.g., keeping the intensity below the damage threshold of optical components) and interpretation of experimental results (e.g., calculating the on-target intensity). A particular feature of high-energy laser systems is that the output-pulse duration is often increased by detuning the stretcher from its optimal settings. This allows for safe amplification and recompression of pulses at higher energy levels because the fluence damage threshold of optical components, in particular diffraction gratings, depends strongly on the pulse duration. In these conditions, the pulse duration must be monitored accurately to avoid amplifying pulses that are shorter than expected to unsafe energy levels. Because of the large scale of these laser systems, the measurement conditions are often nonideal—e.g., wavefront distortions and pointing variations at the end of the system can make it difficult to implement diagnostics based on free-space nonlinear interactions. Streak cameras are typically used as temporal diagnostics on these systems [10], and optical diagnostics have been demonstrated for improved performance [11–13].

This article describes a temporal diagnostic analogous to the phase-diversity technique used for wavefront metrology [14,15] and demonstrates it in the context of pulse characterization on a high-energy chirped-pulse–amplification (CPA) laser system. The spectral phase-diversity technique allows for characterization of a coherent optical pulse that is shorter than the photodetection impulse response. An optical fiber in the near field of the beam couples the pulse into the fiber-based diagnostic, therefore alleviating the impact of wavefront and pointing variations. Precise synchronization is not required because of the long memory length of commercial oscilloscopes. The diagnostic principle is described in Sec. 2. The experimental implementation and algorithms are presented in Sec. 3; experimental results obtained on the OMEGA EP laser [16] are presented in Sec. 4.

2. Principle

The spectral phase-diversity technique is analogous to the spatial phase-diversity technique used to characterize the wavefront φ(x,y) of a coherent beam with fluence distribution F(x,y) in the near field [Fig. 1(a)]. In this technique, the beam to be characterized is focused and fluence distributions are measured at various longitudinal positions in the far field [15]. The measured fluence distributions Gk are expressed in the far-field coordinates system (x′, y′) by

Gk(x,y)=|FT{F(x,y)exp[iφ(x,y)]exp[iπΔzkλf2(x2+y2)]}|2,
where FT is the Fourier transform, λ is the wavelength, f is the focal length of the focusing element, and Δzk represents the distance offsets in the longitudinal direction. The input spatial phase is reconstructed by processing the set of measured fluence distributions. There are many reconstruction approaches in spatial phase diversity—for example, by parametrization of the unknown phase or by projection on a basis such as Legendre polynomials. The unknown wavefront is obtained by minimizing an error metric quantifying the difference between the measured and calculated sets of fluence distributions.

 figure: Fig. 1

Fig. 1 (a) Principle of spatial phase diversity for wavefront measurements and (b) spectral phase diversity for optical pulse characterization. In (a), the spatially resolved fluence Gk(x,y) is measured at various longitudinal locations Δzk close to the focus to reconstruct the near-field phase φ(x,y). In (b), various amounts of spectral phase ψk are introduced on the pulse under test, yielding the instantaneous power Pk(t), and the spectral phase φ(ω) is reconstructed from the powers Pk(t) measured after photodetection with response R(t).

Download Full Size | PDF

Spectral phase diversity is based on the duality between the spectrotemporal variables (ω, t) and the spatial variables (k, x). The spectral phase φ(ω) of the input field is reconstructed using the instantaneous powers Pk(t) that are measured after known spectral phases ψk(ω) have been introduced [Fig. 1(b)]. The optical spectrum of the pulse S(ω) is assumed to be accurately measured because spectral density measurements pose no significant technical difficulty and are routinely performed. (Equivalently, it is routine to measure the fluence beam profile of the beam to be characterized for wavefront measurements via spatial phase diversity.) One significant practical difference with spatial phase diversity is that temporal measurements are convolutions of the actual instantaneous power with the photodetection impulse response R, leading to the measured data set Pk:

Pk=PkR
with

Pk(t)=|S(ω)exp{i[φ(ω)+ψk(ω)+ωt]}dω|2.

Equation (3) is a 1-D equivalent of Eq. (1) when the induced spectral phases are quadratic functions of the optical frequency ω. Reconstruction of the unknown φ from the known S, ψk, R, and Pk can be addressed as a minimization problem. For our application, we chose the parametrization of the spectral phase with its second-order and third-order terms in the Taylor expansion around the central frequency ω0 of the pulse φ(ω) = φ2 (ωω0)2/2 + φ3 (ωω0)3/6. This adequately represents the pulse-shape variations that must be characterized when tuning the stretcher of the CPA system to adjust the output-pulse duration. We found that the measured data were consistent with essentially zero third-order dispersion because the stretcher and compressor had already been matched accurately; therefore, the results presented here correspond to parametrizing the spectral phase with only its second-order coefficient. This is evidently not a limitation of spectral phase diversity because the technique, like its spatial counterpart, can in principle be applied with any parametrization basis provided that the data (measured set of power versus time {Pk} and optical spectrum) can unambiguously be linked to a set of parameters with sufficient accuracy and precision for the targeted application.

When ψk is a quadratic function inducing relatively small changes to the input pulse shape, Eq. (3) can be developed to yield the temporal transport-of-intensity equation [17]. This equation links the dispersion-induced changes in the time-varying power of an optical field to its temporal phase. It has been used to temporally characterize nonlinear phase shifts [17] and optical pulses after temporal stretching to allow for accurate photodetection [18]. The experimental data presented in this article have a relatively low sampling rate and signal-to-noise ratio compared to the data of Refs. 17 and 18 obtained on repetitive signals with a sampling oscilloscope, and an error minimization strategy appears more suitable, although it is not based on a direct analytic data inversion and requires more computing power.

3. Setup and algorithm

The experimental fiber-based spectral phase-diversity setup is described in Fig. 2. The input pulse propagates in a custom-built fiber assembly composed of multiple sequentially connected 2 × 2 fiber splitters. One output port of splitter j and one input port of splitter j + 1 are connected by a short optical fiber. The other pair of ports is connected by a delay fiber with length l × 2j–1. The fiber assembly (seven splitters, l = 4 m) generates N = 64 output pulses with an interpulse delay τ = 20 ns. The fiber assembly serves two purposes. The first purpose is to delay the output pulses relative to each other so that their instantaneous power can be measured with a single photodiode and single-channel oscilloscope. This property has been used with similar fiber assemblies to increase the single-shot photodetection signal-to-noise ratio for narrowband optical pulses. For narrowband pulses, the output pulse shapes are essentially identical and can be summed up after photodetection to provide a lower-noise determination of the input-pulse shape [19,20]. Chromatic dispersion in the optical fibers induces changes in the spectral phase and, subsequently, the instantaneous power of broadband optical pulses. The second purpose of the fiber assembly is to induce different amounts of chromatic dispersion by propagation in different fiber lengths. The chromatic dispersion of a single-mode fiber at 1053 nm is approximately –40 ps/nm/km (φ2 > 0), leading to a relative dispersion of –0.16 ps/nm between successive output pulses and non-negligible pulse-shape changes for a pulse with a bandwidth of the order of a few nanometers or larger. With these parameters, no interference between adjacent output pulses is observed because the interpulse delay (~20 ns) is significantly longer than the temporal extent of each output pulse (at most a few hundred ps). One of the fiber-assembly outputs is connected to a high-bandwidth InGaAs photodiode (DSC10, Discovery Semiconductors) and a real-time high-bandwidth oscilloscope. The diagnostic is entirely built with single-mode fibers from input to the photodiode to alleviate concerns about polarization alignment. We tested two oscilloscopes for this application—a 45-GHz Teledyne Lecroy Wavemaster 8Zi-A and a 70-GHz Tektronix DPO77002SX, leading to photodetection impulse responses at low power levels with a full width at half maximum (FWHM) of 17 ps and 13.5 ps, respectively. The oscilloscopes nominally have a flat frequency response over their stated bandwidth, and the impulse response is mainly limited by the photodiode. The two oscilloscopes have a record length much longer than the ~1.3-μs temporal extent of the 64 pulses with 20-ns interpulse delay. The oscilloscopes were triggered either on the signal itself or with an external trigger synchronized to the pulse under test. The single-shot characterization of isolated optical pulses generated by a high-energy laser system requires a real-time oscilloscope, but the diagnostic can be operated with a sampling oscilloscope if the source under test delivers a train of identical optical pulses that are separated by more than the temporal extent of the fiber-assembly output (~1.3 μs for 64 pulses). Two optical pulses coupled to the two inputs of the first fiber splitter can be characterized independently in a single shot provided that the 64 pulses that are generated by each input pulse do not temporally overlap, or equivalently, that the two input pulses do not temporally overlap at the first fiber splitter. This condition is easily met in practice, and optical fiber can be added to one input in the unlikely event that the pulses overlap. The possibility of characterizing two independent optical pulses with a single setup has been advantageously put into practice to characterize the outputs of the two OMEGA EP CPA beamlines [16].

 figure: Fig. 2

Fig. 2 Setup for spectral phase diversity, based on a fiber assembly and photodetection. The 2 × 2 splitters combined with dispersive delay fibers generate 2S–1 pulses that are photodetected at one output of the assembly, where S is the number of splitters.

Download Full Size | PDF

An example of photodetection impulse response measured with a subpicosecond mode-locked laser is shown in Fig. 3(a), where the typical noise level (~10% peak-to-valley relative to the full voltage range) and relatively low sampling rate (120 GS/s, i.e., one sample every 8.25 ps) can be seen. Figure 3(b) shows that the duration of the impulse response increases with increasing incident power because of nonlinearities in the photodetector [21,22]. Figure 3(c) shows that the measured FWHM varies only by a few picoseconds over a large range of pulse durations that can be generated by our laser system. Measurement noise, a low sampling rate, and variations in impulse response render the task of directly deconvolving measured data to retrieve the optical pulse shape impractical for pulses shorter than the impulse response.

 figure: Fig. 3

Fig. 3 (a) An example of measured impulse response with the DSC10 photodetector and Lecroy 45-GHz oscilloscope; (b) measured full width at half maximum (FWHM) of the impulse response versus signal level; (c) measured FWHM of an optical pulse system versus calculated FWHM when the OMEGA EP stretcher is detuned.

Download Full Size | PDF

Examples of spectral phase-diversity data measured for an input pulse close to best compression and a chirped pulse are shown in the upper and lower row of Fig. 4, respectively. An input pulse with negligible spectral-phase distortions yields chirped output pulses, with chromatic dispersion and temporal extent increasing with the fiber length in which they have propagated (upper row). An input pulse with negative second-order dispersion experiences recompression in the fiber assembly, yielding some output pulses close to the Fourier-transform limit (lower row). An input pulse with positive second-order dispersion is further stretched in the fiber assembly (not shown). The experimental data intuitively indicate the sign and approximate amount of the input-pulse chromatic dispersion from an understanding of the properties of the fiber assembly.

 figure: Fig. 4

Fig. 4 Examples of experimental data for a pulse (upper row) close to best compression and (lower row) after stretcher detuning. [(a),(d)] Waveform measured by the oscilloscope; [(b),(e)] composite trace composed of the power versus time data for the 64 output pulses after retiming; [(c),(f)] power versus time for the 64 output pulses oversampled 32 times by zero-padding for display purposes. Each line of the experimental traces is peak normalized to 1 for clarity.

Download Full Size | PDF

The input-pulse dispersion was determined by comparing the measured data, i.e., the set of instantaneous powers {Pk,meas}, to the set of instantaneous powers {Pk,calc} calculated with a given set of parameters (e.g., second-order and third-order dispersion), the known optical spectrum, chromatic dispersion of the optical fiber, and photodetection impulse response. A Gaussian impulse response with FWHM equal to the extrapolated FWHM of the measured impulse response at a low power level was consistently used. An error metric ε between the two data sets is defined by

ε=kεk2=kdt[Pk,meas(t)Pk,calc(tτk)]2/2,
where each εk is minimized over the delay τk, which takes into account timing shifts containing no meaningful information on the pulse shape. For each εk, the instantaneous powers are normalized by dtPk,meas(t)2=dtPk,calc(t)2=1. By definition, εk = 0 means that the measured and calculated powers are identical and εk = 1 means that their temporal overlap is zero. With this normalization, the minimal value of εk is 10k,max, where 0k,max is the value of the overlap integral 0k (τ) = ∫dtPk,meas(t)Pk,calc(tτ) maximized over τ. The error ε can be calculated in a computationally efficient manner because each 0k, max is simply the maximum value of the Fourier transform of 0k.

Figure 5 shows ε as a function of the input φ2, assuming φ3 = 0, and as a function of φ2 and φ3 for the lower-row data set of Fig. 4. The error is clearly minimized for a single set of input-pulse parameters, for which the calculated experimental trace is shown on the lower row of Fig. 5. Once the spectral phase parameters are identified, the input pulse shape is calculated by Fourier transformation. As previously stated, the retrieved φ3 was consistently insignificantly close to 0 in our data sets; therefore, all further experimental data presented in this article were obtained by error minimization over φ2 only. Simulations presented in Fig. 6 demonstrate that relatively small amounts of third-order dispersion lead to modifications of the experimental trace that should be identifiable by the reconstruction algorithm provided that the measurement signal-to-noise ratio is high enough. Experimental traces were calculated for a Fourier-transform–limited pulse with the measured front-end spectrum (approximately a 20th-order super-Gaussian spectrum with 7.5-nm FWHM bandwidth) and a pulse having residual third-order dispersion, leading to a 10% peak-power decrease, i.e., φ3 = 0.047 ps3 [Figs. 6(a) and 6(b)]. Figure 6(c) shows the calculated power of output pulse #64 for these two pulses. A change in the slope of the power versus time is observed when φ3 is present. An opposite change was observed for negative φ3, and larger changes were observed for larger magnitudes of the third-order dispersion. Figures 6(d)–6(e) display the calculated experimental traces for the two pulses. The differences in the calculated experimental traces indicate that the diagnostic could be used to characterize a pulse originating from a CPA system with non-optimal parameters and tune the incident angle in a CPA grating compressor.

 figure: Fig. 5

Fig. 5 Examples of phase retrieval for the lower-row data of Fig. 4. (a) Error versus input second-order dispersion; (b) error versus input second-order and third-order dispersion; [(c),(d)] calculated 64-pulse composite trace with the determined φ2 at the sampling rate of the oscilloscope and oversampled by a factor 32 with zero padding.

Download Full Size | PDF

 figure: Fig. 6

Fig. 6 Simulations of the impact of third-order dispersion on the experimental trace using two test pulses, a Fourier-transform–limited pulse A (φ2 = 0, φ3 = 0) and a pulse B with third-order dispersion, leading to a 10% decrease in peak power (φ2 = 0, φ3 = 0.047 ps3). (a) Experimental spectrum (black line) and spectral phase of pulse A and B (blue and red lines, respectively). (b) power of pulse A and B (blue and red lines, respectively). (c) The calculated power of pulse #64 from the fiber assembly for input pulse A and B (blue and red lines, respectively). (d) and (e) The calculated experimental traces for pulses A and B, respectively.

Download Full Size | PDF

4. Experimental results

4.1 Laser description

The phase-diversity diagnostic has been tested on the OMEGA EP Laser System. The two short-pulse beamlines can deliver amplified optical pulses with durations ranging from subpicosecond to 100 ps. The output-pulse duration is adjusted on each system by detuning the single-grating Offner-triplet stretcher [23]. Its spherical mirrors are mounted on a rail and can be translated relative to the stretcher grating while keeping the intermirror distance constant. Each beamline has an optical parametric chirped-pulse–amplification (OPCPA) front end that amplifies the stretched output of a mode-locked laser (~6-nm FWHM). Saturation in the two parametric-amplifier stages leads to a flat-top spectrum (FWHM ~7.5 nm) [24]. Most experimental results presented in this article correspond to the single-shot characterization of the OPCPA pulse after propagation through an entire OMEGA EP beamline and recompression. For statistical purposes, ten acquisitions are performed at each set of experimental settings, e.g., stretcher detuning. The pulse under test is coupled in the fiber assembly by a single-mode fiber at an image plane where the beam has a high-order super-Gaussian square profile. The energy required for accurate characterization is typically 50 pJ in the input fiber, which corresponds to ~5 μJ in the free-space beam. Reducing the beam size or focusing into the fiber could increase the sensitivity of the diagnostic, but this might come at the expense of increased sensitivity to beam pointing and wavefront variations.

4.2 Front-end pulse characterization

Figure 7 displays experimental results obtained for pulse durations ranging from best compression (~500-fs pulse) to 20 ps. The retrieved second-order dispersion and duration agree very well with the stretcher model based on design parameters (groove density = 1740 l/mm, angle of incidence = 72.5°) except for pulses close to best compression (~2 ps and shorter), where the magnitude of the retrieved second-order dispersion and pulse duration are overestimated. This effect appears to be almost negligible very close to best compression, where the retrieved pulse duration (600 fs), is close enough to the Fourier-transform pulse duration (~400 fs) and well within the targeted range allowable for our application. Small modulations of the optical spectrum resulting from self-phase modulation (SPM) were identified for this range of input-pulse durations.

 figure: Fig. 7

Fig. 7 [(a),(b)] Retrieved φ2 and [(c),(d)] FWHM versus calculated stretcher dispersion. [(a),(c)] Data obtained with the Lecroy 45-GHz oscilloscope. [(b),(d)] Data obtained with the Tektronix 70-GHz oscilloscope. Simulated values are plotted with a red line. The FWHM of the impulse response of each oscilloscope is shown in blue on the respective plots.

Download Full Size | PDF

The shape of a pulse propagating in the fiber assembly changes continuously because of chromatic dispersion and its energy changes discretely at each fiber splitter. An approximated SPM-induced phase can be calculated by summing up the temporal phases induced at the peak of the pulse by propagation in each delay fiber. For each fiber, the nonlinear phase is simply calculated as ΓLP(t), where Γ is the fiber nonlinear coefficient (estimated at 2 m–1 ⋅ W–1), L is the fiber length, and P(t) is the instantaneous power of the pulse after propagation in half the fiber length. The sum of these phases at the peak of the pulse is calculated for output pulse #64, which has traveled in all the delay fibers of the seven splitter fiber assembly to quantify the magnitude of the nonlinear phase and identify strategies for reducing SPM, as is done in the next subsection. For an input energy of 50 pJ and Fourier-transform–limited input pulses corresponding to the front end (7.5-nm 20th-order super-Gaussian spectrum) and to on-shot conditions (3.5-nm Gaussian spectrum), the total nonlinear phase is approximately 0.6 rad. Nonlinear propagation has a complex impact on the measured experimental trace and the retrieved optical pulse. The induced error depends on the input-pulse properties and its energy, and the diagnostic should in general be configured to minimize SPM and its impact on the measurement.

Figure 8 presents results taken in the dispersion region around –1.73 ps2. This dispersion value corresponds to an on-shot pulse duration of 10 ps, taking into account the gain narrowing in Nd:glass amplifiers that decreases the optical bandwidth from ~7.5 nm to ~3.5 nm. This regime is operationally highly relevant because the maximum energy for safe operation at this pulse duration is significantly higher than at best compression; therefore, a large fraction of scientific shots are performed at 10 ps. The measured data are very consistent with the system model, and there is little spread of the retrieved dispersion around the calculated value. Over the range of stretcher dispersions shown in Fig. 8, the root mean square (rms) and peak-to-valley precision are 0.05 ps2 and 0.15 ps2, respectively, for the 45-GHz oscilloscope. They are slightly lower with the higher-bandwidth 70-GHz oscilloscope, with rms and peak-to-valley precision of 0.03 ps2 and 0.12 ps2, respectively. The higher precision is attributed to the slightly better intrinsic noise performance of the latter oscilloscope compounded with the shorter impulse response. A 10% change in pulse duration around this stretcher setting corresponds to 0.17 ps2. This indicates that pulse reconstruction with precision better than 10% can be obtained on a single-shot basis with this diagnostic.

 figure: Fig. 8

Fig. 8 Retrieved φ2 versus calculated φ2 for (a) the 45-GHz oscilloscope and (b) the 70-GHz oscilloscope. On each plot, the simulated φ2 is plotted with a continuous red line, a 10% margin is plotted with a dashed red line; the solid circles correspond to the average of the ten single-shot determinations plotted with open circles for each stretcher setting.

Download Full Size | PDF

4.3 Investigation of the fiber-assembly design

The impact of the number of pulses and overall dispersion used for spectral phase diversity has been investigated with experimental data measured with the Tektronix 70-GHz oscilloscope. The experimental traces measured with the seven-splitter fiber assembly are composed of the measured instantaneous power of 64 output pulses. Subsets of these powers have been used in the retrieval algorithm to calculate the rms error expressed by Eq. (4) and simulate an experimental trace obtained with a smaller number of splitters. Broadly speaking, this has been done in two different ways:

  • (1) Reducing the number of splitters and approximately preserving the range of dispersions generated by the fiber assembly by using a coarser dispersion granularity. For example, a 16-pulse trace with pulses {1, 5, …, 61} corresponds to a five-splitter assembly with l = 16 m [Fig. 9(a)].
  • (2) Reducing the number of splitters while preserving the granularity of the dispersions generated by the fiber assembly and decreasing the induced-dispersion range. For example, a 16-pulse trace with pulses {1, 2, …, 16} corresponds to a five-splitter assembly with l = 4 m [Fig. 10(a)], and a 16-pulse trace with pulses {49, 50, …, 64} corresponds to this fiber assembly followed by 192 m of fiber [Fig. 10(e)]. The two cases have the same dispersion range but the range is biased toward low dispersion values in the first case and high dispersion values in the second case.
 figure: Fig. 9

Fig. 9 (a) A five-splitter assembly that generates 16 output pulses while approximately preserving the total range of induced dispersions (N = 16, l = 16 m). The corresponding phase-retrieval performance can be assessed by considering pulses {1, 5, …, 61} from the experimental 64-pulse data. (b) Root mean square (rms) and peak-to-valley precision on the retrieved φ2 for the experimental data measured with the Tektronix 70-GHz oscilloscope when a fraction of the experimental trace is used for processing. [(c)–(f)] The retrieved φ2 versus calculated φ2 for 16, 8, 4, and 2 output pulses, respectively. On each plot, the simulated φ2 is plotted with a continuous red line, a 10% margin is plotted with a dashed red line, and the solid circles correspond to the average of the ten single-shot determinations plotted with open circles for each stretcher setting.

Download Full Size | PDF

 figure: Fig. 10

Fig. 10 (a) A five-splitter assembly that generates 16 low-dispersion output pulses while preserving the dispersion granularity (N = 16, l = 4 m). The corresponding phase-retrieval performance can be assessed by considering pulses {1, 2, …, 16} from the experimental 64-pulse data. [(b),(c)] Retrieved φ2 versus calculated φ2 when using 32 and 16 low-dispersion pulses, respectively. (d) A five-splitter assembly that generates 16 high-dispersion output pulses while preserving the dispersion granularity (N = 16, l = 4 m). The assembly is the same as in (a) and is followed by 192 m of fiber. The corresponding phase-retrieval performance can be assessed by considering pulses {49, 50, …, 64} from the experimental 64-pulse data. [(e),(f)] Retrieved φ2 versus calculated φ2 when using 32 and 16 high-dispersion pulses, respectively. On each plot, the simulated φ2 is plotted with a solid red line, a 10% margin is plotted with a dashed red line, and the full markers correspond to the average of the ten single-shot determinations plotted with open circles for each stretcher setting.

Download Full Size | PDF

Figure 9(b) shows the performance impact of using fewer output pulses with an induced dispersion approximately matching the range of the 64-pulse experimental traces. Decreasing the number of pulses from 64 to 32 has no impact on the diagnostic precision. Processing the experimental trace with 16 or fewer pulses leads to higher values of the peak-to-valley and rms error on the retrieved φ2, but these remain within the ± 10% error margin, as can be seen in Figs. 9(c)–9(f). The noise clearly increases but the average retrieved dispersion is still in very good agreement with the calculated stretcher dispersion. A spectral phase-diversity diagnostic with such fiber assembly accurately reconstructs the input-pulse dispersion on average but its precision decreases as the number of splitters is decreased. This leads to poorer single-shot performance.

Figure 10 shows the performance impact of using fewer output pulses while preserving the dispersion granularity. For Figs. 10(a)–10(c), the error function is calculated with the low-index pulses (e.g, {1, 2, …, 16}), i.e., with a small induced dispersion. For Figs. 10(d)–10(f), the error function is calculated with the high-index pulses (e.g, {49, 50,…, 64}), i.e., with a large induced dispersion. A significant loss of accuracy and precision is clearly observed when reducing the number of splitters in this manner, and the dispersion error is higher than 10% when considering 16 pulses [Figs. 10(c) and 10(f)]. Reducing the total range of dispersion used in the spectral phase-diversity diagnostic negatively impacts its performance.

These simulation results indicate that the accuracy of the spectral phase-diversity diagnostic originates mostly from the total range of dispersion induced by the fiber assembly. Reducing the number of splitters can have a significant impact on the diagnostic accuracy when its performance is limited by self-phase modulation. For a given input energy ε, the energy of each output pulse is ε/2S–1 assuming negligible loss and ideal 3-dB fiber splitters. This indicates that the input energy required for a particular photodetected signal scales like 2S–1. Decreasing the number of splitters while preserving the total dispersion range [e.g., using a fiber assembly similar to that shown in Fig. 9(a)] is a path forward to reducing the nonlinear phase accrued in the diagnostic without significantly impacting the performance. An estimate of the nonlinear phase has been calculated similarly to the estimates given in Sec. 4.2 for an input energy that was decreased by a factor 4 and a fiber assembly having five splitters and l = 16 m. For a Fourier-transform–limited pulse with either the 7.5-nm super-Gaussian spectrum or 3.5-nm Gaussian spectrum, the induced nonlinear phase on pulse #16, which has traveled in all the delay fibers, is approximately 0.2 rad. This strategy will be used in our next-generation diagnostic to make possible an accurate temporal characterization over the full range of achievable OMEGA EP pulse durations.

4.4. On-shot amplified pulse characterization

The data retrieved by spectral phase diversity are compared to results from two currently deployed diagnostics in Fig. 11. Figure 11(a) displays a single-shot autocorrelation measured when the stretcher is set to obtain the shortest pulse achievable on the system. This autocorrelation is compared to an autocorrelation calculated with the measured optical spectrum and retrieved dispersion, demonstrating that the best compression stretcher setting has been identified with subpicosecond precision. Figures 11(b) and 11(c) compare the instantaneous power measured by a streak camera with the power reconstructed with the spectral phase-diversity diagnostic.

 figure: Fig. 11

Fig. 11 Comparison of measurements performed with the spectral phase-diversity diagnostic (red line) and (a) a single-shot autocorrelator (green line) and [(b),(c)] a streak camera (blue line). For (a), an OPCPA pulse is characterized at the output of OMEGA EP when the stretcher was set at the best-compression position. For (b) and (c), large-scale Nd:glass amplifiers were fired and the stretcher was set for a 5-ps or 12-ps amplified output pulse, respectively.

Download Full Size | PDF

This comparison was performed on amplified OMEGA EP shots, for which the stretcher was set to yield either a 5-ps pulse or a 12-ps pulse, taking into account the gain narrowing in the Nd:glass amplifiers to ~3.5 nm. Satisfactory agreement was obtained, considering the significant operational complexity and sources of uncertainty associated with the streak camera [10].

5. Conclusions

In conclusion, we have demonstrated a diagnostic based on the concept of spectral phase diversity to characterize picosecond optical pulses with fast photodetection. The experimental trace consists of the photodetected instantaneous power of optical pulses derived from the pulse under test by adding known amounts of chromatic dispersion. An error minimization algorithm retrieves the spectral phase that best matches the measured data, taking into account the diagnostic parameters. The single-shot, fiber-coupled diagnostic has been demonstrated in the context of stretcher tuning for a high-energy laser system and shows good performance over a large range of pulse durations including pulses several times shorter than the photodetection impulse response.

Improved diagnostic performance is expected from technological advances in high-bandwidth photodiodes, high-bandwidth real-time oscilloscopes, and optimization of the fiber assembly used to implement spectral phase diversity. We are currently using simulations to study the potential performance of spectral phase diversity for more-general pulse shapes than those encountered at the output of a chirped-pulse–amplification system.

Acknowledgment

The authors thank R. Cuffney for manufacturing the fiber assembly used for the experiments and Tektronix for the loan of the DPO77002SX oscilloscope. This material is based upon work supported by the Department of Energy National Nuclear Security Administration under Award Number DE-NA0001944, the University of Rochester, and the New York State Energy Research and Development Authority. The support of DOE does not constitute an endorsement by DOE of the views expressed in this article.

References and links

1. I. A. Walmsley and C. Dorrer, “Characterization of ultrashort electromagnetic pulses,” Adv. Opt. Photonics 1, 308–437 (2009).

2. P. Tournois, J.-L. Vernet, and G. Bienvenu, “Sur l’analogie optique de certains montages électroniques: Formation d’images temporelles de signaux électriques,” C. R. Acad. Sci, Ser. B 267, 375–378 (1968).

3. B. H. Kolner, “Space–time duality and the theory of temporal imaging,” IEEE J. Quantum Electron. 30(8), 1951–1963 (1994). [CrossRef]  

4. M. T. Kauffman, W. C. Banyai, A. A. Godil, and D. M. Bloom, “Time-to-frequency converter for measuring picosecond optical pulses,” Appl. Phys. Lett. 64(3), 270–272 (1994). [CrossRef]  

5. C. V. Bennett and B. H. Kolner, “Upconversion time microscope demonstrating 103× magnification of femtosecond waveforms,” Opt. Lett. 24(11), 783–785 (1999). [CrossRef]   [PubMed]  

6. R. Salem, M. A. Foster, A. C. Turner, D. F. Geraghty, M. Lipson, and A. L. Gaeta, “Optical time lens based on four-wave mixing on a silicon chip,” Opt. Lett. 33(10), 1047–1049 (2008). [CrossRef]   [PubMed]  

7. C. V. Bennett, V. J. Hernandez, B. D. Moran, M. E. Lowry, S. P. Vernon, P. T. Steele, D. Chang, and M. M. Fejer, “Time lens based single-shot ultrafast waveform recording: From high repetition rate to high dynamic range,” in Ultrafast Optics 2011 (2011), paper Th20.

8. S. P. Vernon, M. E. Lowry, K. L. Baker, C. V. Bennett, J. R. Celeste, C. Cerjan, S. Haynes, V. J. Hernandez, W. W. Hsing, G. A. LaCaille, R. A. London, B. Moran, A. S. von Wittenau, P. T. Steele, and R. E. Stewart, “X-ray bang-time and fusion reaction history at picosecond resolution using RadOptic detection,” Rev. Sci. Instrum. 83(10), 10D307 (2012). [CrossRef]   [PubMed]  

9. J. D. Zuegel, S. Borneis, C. Barty, B. LeGarrec, C. Danson, N. Miyanaga, P. K. Rambo, C. LeBlanc, T. J. Kessler, A. W. Schmid, L. J. Waxer, J. H. Kelly, B. Kruschwitz, R. Jungquist, E. Moses, J. Britten, I. Jovanovic, J. Dawson, and N. Blanchot, “Laser challenges for fast ignition,” Fusion Sci. Technol. 49, 453–482 (2006).

10. J. Qiao, P. A. Jaanimagi, R. Boni, J. Bromage, and E. Hill, “Measuring 8-250 ps short pulses using a high-speed streak camera on kilojoule, petawatt-class laser systems,” Rev. Sci. Instrum. 84(7), 073104 (2013). [CrossRef]   [PubMed]  

11. I. Jovanovic, B. C. Haefner, M. Shverdin, M. Taranowski, and C. P. J. Barty, “High-dynamic-range, 200-ps window, single-shot cross-correlator for ultrahigh intensity laser characterization,” in Proceedings of CLEO/QELS (2007), paper JThD137.

12. C. Dorrer, J. Bromage, and J. D. Zuegel, “High-dynamic-range single-shot cross-correlator based on an optical pulse replicator,” Opt. Express 16(18), 13534–13544 (2008). [CrossRef]   [PubMed]  

13. C. H. Sarantos and J. E. Heebner, “Solid-state ultrafast all-optical streak camera enabling high-dynamic-range picosecond recording,” Opt. Lett. 35(9), 1389–1391 (2010). [CrossRef]   [PubMed]  

14. R. A. Gonsalves, “Phase retrieval and diversity in adaptive optics,” Opt. Eng. 21(5), 215829 (1982). [CrossRef]  

15. S.-W. Bahk, J. Bromage, I. A. Begishev, C. Mileham, C. Stoeckl, M. Storm, and J. D. Zuegel, “On-shot focal-spot characterization technique using phase retrieval,” Appl. Opt. 47(25), 4589–4597 (2008). [CrossRef]   [PubMed]  

16. D. N. Maywar, J. H. Kelly, L. J. Waxer, S. F. B. Morse, I. A. Begishev, J. Bromage, C. Dorrer, J. L. Edwards, L. Folnsbee, M. J. Guardalben, S. D. Jacobs, R. Jungquist, T. J. Kessler, R. W. Kidder, B. E. Kruschwitz, S. J. Loucks, J. R. Marciante, R. L. McCrory, D. D. Meyerhofer, A. V. Okishev, J. B. Oliver, G. Pien, J. Qiao, J. Puth, A. L. Rigatti, A. W. Schmid, M. J. Shoup III, C. Stoeckl, K. A. Thorp, and J. D. Zuegel, “OMEGA EP high-energy petawatt laser: progress and prospects,” J. Phys. Conf. Ser. 112(3), 032007 (2008). [CrossRef]  

17. C. Dorrer, “Characterization of nonlinear phase shifts by use of the temporal transport-of-intensity equation,” Opt. Lett. 30(23), 3237–3239 (2005). [CrossRef]   [PubMed]  

18. C. Cuadrado-Laborde, A. Carrascosa, P. Pérez-Millán, A. Díez, J. L. Cruz, and M. V. Andres, “Phase recovery by using optical fiber dispersion,” Opt. Lett. 39(3), 598–601 (2014). [CrossRef]   [PubMed]  

19. J. R. Marciante, W. R. Donaldson, and R. G. Roides, “Averaging of replicated pulses for enhanced-dynamic-range single-shot measurement of nanosecond optical pulses,” IEEE Photonics Technol. Lett. 19(18), 1344–1346 (2007). [CrossRef]  

20. C. Dorrer, “Analysis of pulse replicators for high-bandwidth, high-dynamic-range, single-shot optical characterization,” J. Lightwave Technol. 31(9), 1374–1382 (2013). [CrossRef]  

21. K. J. Williams and R. D. Esman, “Photodiode DC and microwave nonlinearity at high currents due to carrier recombination nonlinearities,” IEEE Photonics Technol. Lett. 10(7), 1015–1017 (1998). [CrossRef]  

22. T. S. Clement, P. D. Hale, K. C. Coakley, and C. M. Wang, “Time-domain measurement of the frequency response of high-speed photoreceivers to 50 GHz,” in Technical Digest:Symposium on Optical Fiber Measurements,2000, P. A. Williams and G. W. Day eds., (National Institute of Standards and Technology, 2000), pp. 121−124.

23. G. Chériaux, P. Rousseau, F. Salin, J. P. Chambaret, B. Walker, and L. F. Dimauro, “Aberration-free stretcher design for ultrashort-pulse amplification,” Opt. Lett. 21(6), 414–416 (1996). [CrossRef]   [PubMed]  

24. C. Dorrer, A. Consentino, D. Irwin, J. Qiao, and J. D. Zuegel, “OPCPA front end and contrast optimization for the OMEGA EP kilojoule, picosecond laser,” J. Opt. A, Pure Appl. Opt. 17, 094007 (2015).

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (11)

Fig. 1
Fig. 1 (a) Principle of spatial phase diversity for wavefront measurements and (b) spectral phase diversity for optical pulse characterization. In (a), the spatially resolved fluence Gk(x,y) is measured at various longitudinal locations Δzk close to the focus to reconstruct the near-field phase φ(x,y). In (b), various amounts of spectral phase ψk are introduced on the pulse under test, yielding the instantaneous power Pk(t), and the spectral phase φ(ω) is reconstructed from the powers P k (t) measured after photodetection with response R(t).
Fig. 2
Fig. 2 Setup for spectral phase diversity, based on a fiber assembly and photodetection. The 2 × 2 splitters combined with dispersive delay fibers generate 2S–1 pulses that are photodetected at one output of the assembly, where S is the number of splitters.
Fig. 3
Fig. 3 (a) An example of measured impulse response with the DSC10 photodetector and Lecroy 45-GHz oscilloscope; (b) measured full width at half maximum (FWHM) of the impulse response versus signal level; (c) measured FWHM of an optical pulse system versus calculated FWHM when the OMEGA EP stretcher is detuned.
Fig. 4
Fig. 4 Examples of experimental data for a pulse (upper row) close to best compression and (lower row) after stretcher detuning. [(a),(d)] Waveform measured by the oscilloscope; [(b),(e)] composite trace composed of the power versus time data for the 64 output pulses after retiming; [(c),(f)] power versus time for the 64 output pulses oversampled 32 times by zero-padding for display purposes. Each line of the experimental traces is peak normalized to 1 for clarity.
Fig. 5
Fig. 5 Examples of phase retrieval for the lower-row data of Fig. 4. (a) Error versus input second-order dispersion; (b) error versus input second-order and third-order dispersion; [(c),(d)] calculated 64-pulse composite trace with the determined φ2 at the sampling rate of the oscilloscope and oversampled by a factor 32 with zero padding.
Fig. 6
Fig. 6 Simulations of the impact of third-order dispersion on the experimental trace using two test pulses, a Fourier-transform–limited pulse A (φ2 = 0, φ3 = 0) and a pulse B with third-order dispersion, leading to a 10% decrease in peak power (φ2 = 0, φ3 = 0.047 ps3). (a) Experimental spectrum (black line) and spectral phase of pulse A and B (blue and red lines, respectively). (b) power of pulse A and B (blue and red lines, respectively). (c) The calculated power of pulse #64 from the fiber assembly for input pulse A and B (blue and red lines, respectively). (d) and (e) The calculated experimental traces for pulses A and B, respectively.
Fig. 7
Fig. 7 [(a),(b)] Retrieved φ2 and [(c),(d)] FWHM versus calculated stretcher dispersion. [(a),(c)] Data obtained with the Lecroy 45-GHz oscilloscope. [(b),(d)] Data obtained with the Tektronix 70-GHz oscilloscope. Simulated values are plotted with a red line. The FWHM of the impulse response of each oscilloscope is shown in blue on the respective plots.
Fig. 8
Fig. 8 Retrieved φ2 versus calculated φ2 for (a) the 45-GHz oscilloscope and (b) the 70-GHz oscilloscope. On each plot, the simulated φ2 is plotted with a continuous red line, a 10% margin is plotted with a dashed red line; the solid circles correspond to the average of the ten single-shot determinations plotted with open circles for each stretcher setting.
Fig. 9
Fig. 9 (a) A five-splitter assembly that generates 16 output pulses while approximately preserving the total range of induced dispersions (N = 16, l = 16 m). The corresponding phase-retrieval performance can be assessed by considering pulses {1, 5, …, 61} from the experimental 64-pulse data. (b) Root mean square (rms) and peak-to-valley precision on the retrieved φ2 for the experimental data measured with the Tektronix 70-GHz oscilloscope when a fraction of the experimental trace is used for processing. [(c)–(f)] The retrieved φ2 versus calculated φ2 for 16, 8, 4, and 2 output pulses, respectively. On each plot, the simulated φ2 is plotted with a continuous red line, a 10% margin is plotted with a dashed red line, and the solid circles correspond to the average of the ten single-shot determinations plotted with open circles for each stretcher setting.
Fig. 10
Fig. 10 (a) A five-splitter assembly that generates 16 low-dispersion output pulses while preserving the dispersion granularity (N = 16, l = 4 m). The corresponding phase-retrieval performance can be assessed by considering pulses {1, 2, …, 16} from the experimental 64-pulse data. [(b),(c)] Retrieved φ2 versus calculated φ2 when using 32 and 16 low-dispersion pulses, respectively. (d) A five-splitter assembly that generates 16 high-dispersion output pulses while preserving the dispersion granularity (N = 16, l = 4 m). The assembly is the same as in (a) and is followed by 192 m of fiber. The corresponding phase-retrieval performance can be assessed by considering pulses {49, 50, …, 64} from the experimental 64-pulse data. [(e),(f)] Retrieved φ2 versus calculated φ2 when using 32 and 16 high-dispersion pulses, respectively. On each plot, the simulated φ2 is plotted with a solid red line, a 10% margin is plotted with a dashed red line, and the full markers correspond to the average of the ten single-shot determinations plotted with open circles for each stretcher setting.
Fig. 11
Fig. 11 Comparison of measurements performed with the spectral phase-diversity diagnostic (red line) and (a) a single-shot autocorrelator (green line) and [(b),(c)] a streak camera (blue line). For (a), an OPCPA pulse is characterized at the output of OMEGA EP when the stretcher was set at the best-compression position. For (b) and (c), large-scale Nd:glass amplifiers were fired and the stretcher was set for a 5-ps or 12-ps amplified output pulse, respectively.

Equations (4)

Equations on this page are rendered with MathJax. Learn more.

G k ( x , y )= | FT{ F( x,y ) exp[ iφ( x,y ) ]exp[ i πΔ z k λ f 2 ( x 2 + y 2 ) ] } | 2 ,
P k = P k R
P k ( t )= | S( ω ) exp{ i[ φ( ω )+ ψ k ( ω )+ωt ] }dω | 2 .
ε= k ε k 2 = k dt [ P k,meas ( t ) P k,calc ( t τ k ) ] 2 / 2,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.