1. Introduction

Interferometric techniques have always been used to characterize ultrashort laser pulses [1–2]. For a long time it was thought that only pure interferometric characterization techniques such as SPIDER [3–4] were able to directly retrieve the phase and amplitude information from an ultrashort laser pulse. Other characterization techniques, specifically time-frequency measurements such as FROG, have required iterative retrieval algorithms to extract phase information [5, 6]. Recently, there has been significant interest in analyzing information within a frequency resolved interferometric autocorrelation traces and has resulted in a new iterative retrieval algorithm being developed [7] as well as a methodology to convert such a trace into a standard FROG trace [8–9]. Additionally, we have recently experimentally demonstrated a new interferometric-based pulse characterisation technique called MEFISTO (Measuring the electric field by interferometric spectral trace observation) [10]. MEFISTO has allowed for the first time, the analytical characterization of an unknown ultrashort pulse from a time-frequency representation of a pulse.

MEFISTO, as SPIDER, allows the pulse information to be directly extracted from experimental data. However, unlike SPIDER, it is based on recording a time-frequency measurement of the pulse which makes it strongly related to FROG. This ambivalent nature makes it possible to share some of the advantages attributed to both methodologies. On one side, as in SPIDER, the extraction of the phase can be directly performed without the need of an iterative retrieval algorithm. On the other side, MEFISTO possesses the same error checking capabilities and simple experimental arrangement as FROG allowing it to be used in a broader range of applications where collinear geometry is required [9]. Furthermore, as we show in this paper, the time-frequency nature of the method makes possible, as in XFROG [11, 12], the simultaneous characterization of two unknown ultrashort pulses, a technique we refer to as blind-MEFISTO.

The aim of this paper is to outline the theoretical background and discuss the retrieving analytical procedure and practical implications in detail. Our main goal here is to present the generalized theory for Blind-MEFISTO and then go on to discuss the particular properties of the degenerate case where the two pulses are equal (MEFISTO). For this purpose, the paper has been organized as follows: In section 2, we develop the general theory. This is followed by an example of the blind-MEFISTO retrieval. In section 4 some experimental and practical consideration are pointed out. A discussion on the degenerate case is presented in section 5. Finally, general conclusions are given.

2. Theory

MEFISTO is based on the Fourier analysis of a spectrally resolved interferometric correlation trace. The experimental set up to obtain such a trace can be seen in Fig. 1.

 

Fig. 1. Typical experimental set up needed to obtain a Frequency resolved interferometric correlation trace.

Download Full Size | PDF

Two pulses interact collinearly within a nonlinear medium, one of them after passing through a delay arm. The second harmonic generated signal is then directed to a spectrograph to obtain the interferometric time-frequency trace in terms of the time-delay τ and the frequency f. An example of the resulting trace can be seen in Fig. 2(a). In a general case, this trace can be mathematically described as

where E(t) and G(t) are the slowly varying amplitude of the complex electric field centered at the frequencies f ₁ and f ₂. The Fourier transform with respect to the variable t is indicated by F_t and χ is related with the conversion efficiency in the nonlinear process. Note that the two interfering pulses E(t) and G(t) can be either different, what is referred to us as Blind-MEFISTO, or equal, which is the degenerate case that for simplicity is solely named MEFISTO. In Blind-MEFISTO, phase matching must be achieved for all the cross-terms in (1) and, as we will show, the new information carried on the interferometric terms will allow the simultaneous analytical determination of E(t) and G(t). In general, this should restricts the use of Blind-MEFISTO to the case in which the central frequencies of E(t) and G(t) are equal (f ₁=f ₂=f ₀).

We start our analysis by calculating the Fourier transform of Eq. (1) in the τ axis, i.e., Y^SHG (f,κ)=F _τ{I^SHG (f,τ)}. The resulting expression consist of 5 main spectral components [see Fig. 2(b)] at delay-frequencies κ=0, ±f ₀ and ±2f ₀. Since the interferometric trace [Fig. 2(a)] is real, the negative spectral components in Fig. 2(b) are the complex conjugate of the positive ones. We therefore only need to focus upon one of the components (either the positive or the negative) to analyze the information enclosed in the transformed trace. Each of these terms contain information of the pulse phase and intensity and their use will depend on the particular experimental conditions. To highlight the information enclosed in Y^SHG (f,κ), we will separately focus on the mathematical expression for each component and analyze their possibilities.

 

Fig. 2. (a) Frequency resolved interferometric correlation trace and (b) its Fourier transform in the delay-frequency axis, showing the different spectral components.

Download Full Size | PDF

We start analyzing the component at κ=0, which can be written as:

where the second harmonic of the two unknown fields E_SHG (f) and G_SHG (f), are related with the fundamental pulse following F_SHG (f)=${\int }_{-\infty }^{\infty }$ df′F(f′)F(f-f′), with f corresponding to the base-band frequency [13]. The second term in Eq. (2) includes the standard XFROG trace that has been used to characterize ultra short laser pulses [11,12]. It is worth noting that this term, with the use of adequate phase matching conditions, can be obtained even in the case where the spectra of the two pulses do not overlap (f ₁≠f ₂) allowing XFROG to be used to characterise pulses that have differing central frequencies. In a previous paper we outlined a filtering procedure capable to extract the FROG term from the κ=0 component in the degenerate case (E(t)=G(t)). This technique named as CFROG [8] can straightforwardly be extended to obtain the XFROG term. However, when retrieving the pulses, it again relies upon the use of retrieval algorithms since the integral on Eq. (2) does not allow the analytical determination of E(t).

We continue our analysis by focusing upon the spectral component of Y^SHG (f,κ) at κ≈2f ₀. This component can be written as

This expression corresponds to a straight line in the transformed trace [Fig. 2(b)] given by κ=f+2f ₀. This term can be used for different purposes. Firstly, in an experimental trace, the divergence from an ideal delta function gives a measure of the pulse jitter, as suggested in Ref. [7], or in general, the quality of the measurement. In addition, any deviation from a linear dependency between κ (real frequency) and f (base band frequency) will indicate a spectrograph miss-calibration. This component can also determine E(t) when the gating function, G(t), is already known. It can therefore be used as a further verification of the blind-MEFISTO retrieval by checking that the two pulses fulfill

Although the terms at κ=0 and ±2f ₀ are useful to obtain some pulse information, it is the component at κ≈f ₀ of Y^SHG (f,κ) that allows the simultaneous analytical characterization of two unknown ultrashort pulses. The mathematical expression for this component can be written as:

This equation shows the main advantage of Fourier transforming the trace as was done in Fig. 2(b): All the convolutions that appear in the time domain have been transformed into products in the frequency domain. This isolates the individual pulse fields and makes it possible to resolve analytically. To show this, we firstly write all the involved complex magnitudes in polar form, i.e., E(f)=E ₀ U(f)exp(iϕ(f)), G(f)=G ₀ V(f) exp(iγ(f)), and equivalently for the harmonic pulses E_SHG (f)=$E_0^{\mathit{\text{SHG}}}$ U_SHG (f)exp(iϕ_SHG (f)) and G_SHG (f)=$G_0^{\mathit{\text{SHG}}}$ V_SHG (f)exp(iϕ_SHG (f)). Here, for convenience, all spectral profiles U(f) are normalized at the central base-band frequency (U(f=0)=1). Then, taking Y^SHG (f,κ)=R(f,κ) exp(iθ(f,κ)), Eq. (5) can be written as

Here we have defined two parameters χ ₁=${\chi E}_0^{\mathit{\text{SHG}}}$ E ₀ G ₀ and χ ₂=${\chi G}_0^{\mathit{\text{SHG}}}$ E ₀ G ₀, which can be understood as effective conversion efficiency parameters. In section 4, we will go on to describe a procedure to obtain χ ₁ and χ ₂ from the experimental trace. In what follows, and for clarity purposes, we will assume that these parameters are known. Under typical lab conditions the normalized spectra profiles, U(f), V(f), U_SHG (f) and V_SHG (f) can be easily measured. As R(f,κ) and θ(f,κ) are obtained from the interferometric trace, therefore, the only unknowns in Eq. (6) corresponds to the phase of the fundamental and second harmonic pulse, i.e., ϕ(f), γ(f),ϕ_SHG (f) and γ_SHG (f) . In order to fully characterize the pulses these unknowns must be calculated. This can be achieved by first taking two different slices in the transformed space of the interferometric trace, e.g., at κ=f ₀ and κ=f ₀-Δf. Then, by taking the real and imaginary parts in Eq. (6), we can isolate the phase component at κ=f ₀, obtaining

and at κ=f ₀-Δf

where we have defined two functions that relate the different pulse spectra and the trace as

${\Omega }_1(f,\kappa )=\frac{R^2(f,k)+4\left[{\chi }_1^2{U}_{\mathit{SHG}}^2\left(f\right)U^2(f+f_0-\kappa )V^2(\kappa -f_0)-{\chi }_2^2{V}_{\mathit{SHG}}^2\left(f\right)V^2(f+f_0-\kappa )U^2(\kappa -f_0)\right]}{4{\chi }_{1}R(f,\kappa )U_{\mathit{SHG}}\left(f\right)U(f+f_0-\kappa )V(\kappa -f_0)}$

and

${\Omega }_2(f,\kappa )=\frac{R^2(f,k)-4\left[{\chi }_1^2{U}_{\mathit{SHG}}^2\left(f\right)U^2(f+f_0-\kappa )V^2(\kappa -f_0)-{\chi }_2^2{V}_{\mathit{SHG}}^2\left(f\right)V^2(f+f_0-\kappa )U^2(\kappa -f_0)\right]}{4{\chi }_{2}R(f,\kappa )V_{\mathit{SHG}}\left(f\right)V(f+f_0-\kappa )U(\kappa -f_0)}$

Then, by subtracting Eqs. 7(a) and 7(c) we get

From this equation the phase for the spectral component E(f), i.e., ϕ(f) can be determined by taking an arbitrary origin ϕ(0) and adding up Δϕ_k (f) as

In a similar way, by subtracting Eqs. 7(b) and 7(d), we obtain the equation necessary to determine the phase increment for G(f) in terms of the frequency as

and, as in the previous case, by taking an arbitrary origin γ(0) find the spectral phase components as

Notice that the terms γ(0)-γ(-Δf) and ϕ(0)-ϕ(-Δf) in eqs. (8) and (10) are unknown constants. These terms add a linear spectral phase shift that only affects the electric field time origin and therefore they can be decided arbitrarily without affecting the shape of the temporal pulse envelope. Equations (8–11) are the principal result of this work showing that to simultaneously characterize two different pulses in an analytical way is possible.

3. Simultaneous characterization of two unknown pulses: blind-MEFISTO

In this section we will numerically go through an analytical retrieval step by step in order to show details of how the method should be performed in an experimental case.

 

Fig. 3. Experimental data necessary to solve Eqs. (8) and (10). In (a) we show the amplitude and in (b) the phase corresponding to the two slices at κ=f ₀ and κ=f ₀+Δf obtained from the transformed trace in Fig. 2(b). We also show the spectra for the two fundamental pulses (c) and at the SHG frequency (d).

Download Full Size | PDF

First, as commented we need to take the two slices (at κ=f and at κ=f ₀+Δf) of the Fourier transformed interferometric trace shown in Fig. 2(b). The amplitude of the trace R(f,κ) and the phases θ(f,κ) are respectively shown in Fig. 3 (a) and 3(b) at κ=f ₀ and κ=f ₀-Δf. In addition, to evaluate the Ω₁(f,κ) and Ω₂(f,κ) functions, the two fundamental pulses spectra are needed [Fig. 3(c)], which must be obtained experimentally. The two spectra of the second harmonic pulses [Fig. 3(d)] and the parameters χ ₁ and χ ₂ can be either obtained experimentally or directly from the interferometric trace using the method outlined in section 4 (Eqs. 15–16).

In an experimental case Ω₁(f,κ) and Ω₂(f,κ) are therefore evaluated using solely experimental data. The resulting Ω₁(f,κ) and Ω₂(f,κ) function are shown, respectively, in Figures 4(a) and 4(b), at κ=f ₀ and κ=f ₀-Δf.

 

Fig. 4. (a) Function Ω₁(f,κ) and (b) Ω₂(f,κ) at κ=f ₀ and κ=f ₀+Δf. Dashed lines shows the division between domains where the function cos^-1[Ω(f,κ=f ₀)] alternates the sign.

Download Full Size | PDF

Note that the results shown in Figs. 3(b), 4(a), and 4(b) are the necessary data to compute the phase increments in Eqs. (8) and (10). Then, these increments are used in Eqs. (9) and (11), to built up the spectral phases for both pulse. Figures 5(a)–5(b) show the calculated spectral phase obtained with our procedure.

 

Fig. 5. Retrieved pulses (a) E(f) and (b) G(f) using Eqs. (8) and (10) without alternate sign in when changing between the domains shown in Fig. 4.

Download Full Size | PDF

As can be seen, the results only match with the exact phase within a certain region. This originates from not being able to determine the sign of the phase, which is caused by the sign indetermination related with the cos^-1(Ω) function present in Eqs. (8) and (10). Taking into account that cos^-1(Ω) should change continuous along the frequency axis in the [-π,π] range, the problem can be overcome identifying those points where Ω(f,κ) get close to 1 or -1 and there is an abrupt change in the slope. In these positions an alternate sign in the cos^-1(Ω) function should be performed. These sign-switching points (SSP) define the different domains where the sign in the cos^-1(Ω) function changes. These are highlighted in Figures 4(a) and 4(b) with dashed vertical lines. It is worth to mention that finding the right SSP is fundamental to ensure a proper result. In most of the cases, these can be determined by eye inspection or using a simple algorithm. However, for abrupt changes in the Ω(f,κ) function, this determination can be difficult. This can be the case for pulses with spike spectra or large noise. These two situations are commented in the next section.

After applying the criterion described above to Eqs. (8) and (10), and again using Eqs. (9) and (11) to build up the phase, two possible solutions are available. These two solutions originate from the ambiguity in the sign of Eqs. (8) and (10) at the first domain. However, only one of these solutions is correct. The correct one, observed in Figures 6(a)–6(b), results in the spectral phase that best reproduces the original phase. The second solution, when the other sign is chosen at the first donain, greatly differs from the original phase [Figs. 6(c)–6(d)]. This second solution is a spurious solution of our equation system that correctly reproduces the component at term κ≈f ₀ in Fig. 2(b) but it is not solution of the complete interferometric trace. In particular, the pulses in Figs. 6(c)–6(d) do not fulfil Eq. (4), that, as commented before, can be used as error-checking procedure. Experimentally, any source of errors, including the spurious solution, can be detected by numerically generating the interferometric trace using the retrieved pulses and comparing the result with the experimental trace. Equivalently, we can compare the interferometric correlations of the two solutions with the experimental one, which can be directly obtained from the interferometric trace time marginal. We have demonstrated this in Fig. 7. Here, only the interferometric correlation corresponding to the solution in Figs. 6(a)–6(b) is identical to the trace time marginal obtained by integrating the frequency axis of the interferometric trace in Fig. 2(a). Apart from the spurious solution, which can be rejected using the two checking procedures described above, we have not been able to detect other ambiguities as the ones present in XFROG and discussed in Ref. [14].

 

Fig. 6. Obtained pulses alternating the sign at the domains shown in Fig. 4. Result (a) E(f) and (b) G(f) using the sign combination (+,-) in Eqs. (8) and (10). The spurious solution is shown in (c) E(f) and (d) G(f) using the (-,+) sign combination. Here the discontinuity in the phase corresponds to a 2π change.

Download Full Size | PDF

 

Fig. 7. E(f) and G(f) calculated interferometric correlations (blue) and envelope of the interferometric correlation obtained from Fig. 2(a) as the time marginal of the interferometric trace (red). (a) Result using the solution in Figs. 6(a)–6(b). Here the red line coincides with the blue contour and is barely visible. (b) Result using the solution in Figs. 6(c)–6(d).

Download Full Size | PDF

4. Experimental and practical considerations

In this section we discuss important characteristics of MEFISTO that should be considered to ensure easier and successful pulse retrievals.

4.1. Pulse bandwidth

An important issue to consider is the maximum pulse bandwidth that MEFISTO can successfully characterize. This fundamental limit derives from the necessity to prevent overlap between the different components of the transformed trace [see Fig. 2(b)]. Specifically, the two slices at κ=f ₀ and κ=f ₀-Δf must not be affected by the tails of the κ=0 component. In principle, this suggests that the available bandwidth could be equal to the optical carrier, i.e., Δλ _max≈λ ₀. However, in practice, this can be affected by the particular spectral shape.

4.2. Frequency resolution

In our analysis we have implicitly considered that the sampling step in the f and κ axis coincides (Δκ=Δf). Initially, the frequency resolution is in fact given by the time-delay span τ_span , i.e., Δκ=1/τ_span , which experimentally can differ from the spectrograph resolution. However, by using interpolating techniques to fulfil Δκ=Δf, the frequency resolution of the method can be extended to Δf.

In the previous section, we pointed out that large changes in the Ω(f,κ) function can result in errors in determining the SSP. This source of errors can be reduced by increasing the frequency resolution. We must ensure, therefore, that the sampling, and therefore the resolution, results in a reasonable number of points inside the bandwidth.

4.3. Delay resolution.

The delay resolution is determined by the sampling step Δτ. This is elected following Nyquist criterion to resolve the inteferometric fringes: Δτ<1/2f _max<1/4f ₀. From the experimental point of view this results in large data sets and acquisition times (of the order of seconds) [10]. The implementation of a faster procedure can be obtained using undersampling techniques in a similar way to the procedure described in Ref. [8]. A detailed investigation of this situation is beyond the scope of this paper and will be reported elsewhere.

4.4. Determining the conversion efficiency parameters χ₁ and χ₂

In our theoretical development, we have considered all the spectral profiles to be normalized to unity at the central wavelength. The energy of each pulse is therefore associated with the spectral amplitudes E ₀, $E_0^{\mathit{\text{SHG}}}$, G ₀ and $G_0^{\mathit{\text{SHG}}}$. These terms, together with the nonlinear coefficient χ are included in the conversion efficiency parameters χ ₁ and χ ₂. Here we will outline how these parameters can be calculated using the κ=0 and κ=2f ₀ components of the transformed interferometric trace. To achieve this we must first define three separate equations, two that are derived from the κ=0 term and one that is derived from the κ=2f ₀+f term. Firstly, using the same procedure outlined in Ref. [8], the two contributions to the κ=0 term in Eq. (2) can be separated. The first part, involving the delta function, when κ=0 gives

and the second part, involving the CFROG information, when κ=0 and f=0, yields,

where P=${\int }_{-\infty }^{\infty }$ df′U ²(f′)V ²(-f′) can be experimentally determined from the experimental spectra I_E (f)=U ²(f) and I_G (f)=V ²(f). The third equation is obtained by looking at the component at κ=2f ₀+f and taking its modulus;

After some algebra, Eq. (12)–Eq. (14) lead to

where all terms on the right hand side are experimental parameters. The sign before the square root depends on the peak intensities of the two second harmonic pulses and must be chosen accordingly so that when, for example, $I_E^{\mathit{\text{SHG}}}$ >$I_G^{\mathit{\text{SHG}}}$ then χ ₁>χ ₂. By determining Eqs. (15) and (16) we have shown how it is possible to evaluate the functions Ω₁(f,κ) and Ω₂(f,κ) without the need of experimentally measuring the intensities of the second harmonic spectra U_SHG (f) and V_SHG (f).

4.5. Pectral calibration

An experimental measurement, in general, requires two independent measurements: the two fundamental pulse spectra and the interferometric trace. The need for the exact determination of the f ₀ and 2f ₀ places extreme demands and importance on the spectral calibration. These demands can be relaxed significantly by evaluating the second term of Eq. (2) at κ=0, and using it to directly relate f ₀ and 2f ₀ with one another:

By taking a convolution of the two measured fundamental spectra, an accurate spectral registration can be performed that relates the region of the fundamental spectrum with the region of the second harmonic spectrum, where the interferometric trace is obtained. This property allows the MEFISTO technique to be extremely roboust against spectrometer miscalibration.

4.6. Noise.

The presence of noise in the experimental measurement can affect the phase extraction as noise increase the difficulty in determining the SSPs. As an example of the limits that noise impose to the technique we have added white Gaussian noise of 20 dB peak signal-to-noise ratio (SNR) to the slide R(f,κ≅f ₀), shown in Fig. 3(a). Under these conditions it was still possible to detect the SSPs and we were able to obtain the spectral phase inside the spectrum bandwidth [Fig. 8(a)]. A further rising in the noise level, increased the difficulty in distinguishing the SSPs. To overcome this problem, we reject the high-frequency noise components by filtering the Ω₁(f,κ) and Ω₂(f,κ) functions in the f axis. This procedure allows the phase to be obtained at noise levels as high as 13 dB SNR, as shown in Fig. 8(b).

 

Fig. 8. Retrieved phase for a pulse with (a) 20 dB and (b) 13 dB SNR. Solution in (b) has been obtained by low band pass filtering the functions Ω₁(f,κ) and Ω₂(f,κ) in the frequency axis.

Download Full Size | PDF

5. Characterization of an unknown single pulse: degenerate case

This section is devoted to the degenerate case, i.e., when the two interacting pulses are identical. We have experimentally demonstrated MEFISTO for this situation in a previous work [10], here we discuss in more detail some aspect that are particular of this case and are not present in the general case of blind-MEFISTO.

To start, Eqs. (8) and (10) are identical and can be reduced to

In the degenerated case we must take into account that the interferometic trace is symmetric with respect to the delay axis τ thus its Fourier transform must be real. This has two consequences. First, the upper term in the Ω(f,κ) functions is directly Y^SHG (f,κ), where Ω(f,κ) now is written as

$\Omega (f,\kappa )=\frac{Y^{\mathit{SHG}}(f,k)}{4{\chi }_{\mathit{eff}}{U}_{\mathit{SHG}}\left(f\right)U\left(f+f_0-\kappa \right)U\left(\kappa -f_0\right)}$

with χ_eff =${\chi E}_0^{\mathit{\text{SHG}}}$ $E_0^2$. In an experimental situation, any imaginary component in Y^SHG (f,κ) can be attributed to experimental errors. This is because the interferometric nature of this technique requires having a perfectly centered trace, exact delay steps, no laser instabilities, etc. In this situation, when the center of the delay axis is accurately known, the imaginary part can be omitted. This is the equivalent to the symmetrization process performed in the FROG technique. Retaining the absolute value is however preferred in most of the cases since experimental measurement are never performed under an ideal conditions and the delay origin is unknown. The second consequence of the trace being symmetric is that a positive and negative sign in front of Eq. (18) will give the as result E(f) and E*(f). This is equivalent to the intrinsic ambiguity that appears in SHG-FROG measurements.

As with blind-MEFISTO, an important factor that affects the determination of Ω(f,κ) is the effective conversion efficiency parameter χ_eff . In the degenerate case, the procedure to obtain χ_eff is very much simplified. In contrast to blind-MEFISTO where an external measurement of the fundamental spectra U ²(f) is required, the degenerate case allows all the information to be directly obtained from the experimental interferometric trace. This is achieved as follows. First, the amplitude of the second harmonic pulse are obtained from the term at κ=2f ₀+f, resulting in

By normalizing using U_SHG (f=0), the spectral profile U_SHG (f) and the factor √χ·$E_0^{\mathit{\text{SHG}}}$ needed to find χ_eff are obtained. Second, the term at κ=0 in the degenerate case results in

By using the inverse Fourier transform to deconvolve (20) we obtain:

Here again, by normalising U(f=0)=1, the spectral profile of the fundamental pulse U(f) and the factor √χ·$E_0^2$ needed to finally find χ_eff are obtained. By not requiring an extra spectral measurement, all the data used for the pulse measurement is contained within a single set of data. This helps keep everything self-consistent and as a consequence reduces the risk of erroneous errors entering into the final results.

Finally, as in Blind-MEFISTO, the degenerate case can compare the experimental and numerically generated interferometric traces to check the validity of the method. As commented, the experimental interferometric auto/cross correlation trace is obtained by integrating the trace in frequency. This can be compared with the one numerically obtained from the retrieved pulse. With the degenerate case however an even simpler check can be made by making use of the intrinsic 8:1 ratio that intefermetric autocorrelation must have. By integrating the experimental data to obtain the time marginal immediately after acquisition, it is possible to check for the 8:1 ratio. If the ratio is not correct, the data can immediately be discarded and taken again. It should also be noted that in the degenerate case, the use of Eq. (4) as an error checking procedure is not possible, since this component always results in the second harmonic intensity spectrum.

5. Conclusions

In this work we have outlined, the theory and procedure of MEFISTO, a method that allows the complex amplitude of two unknown ultrashort pulses with similar central frequency to be simultaneously obtained. Its time-frequency and interferometric ambivalent nature leads to a technique that possesses advantages and properties characteristic of interferometric techniques as SPIDER, as well as time-frequency measurements as FROG. In this sense, similar to interferometric techniques, the method enables the analytical extraction of pulse information without the need of an iterative retrieval algorithm. Similar to time-frequency techniques, it can be applied to simultaneously characterize two unknown pulses (blind MEFISTO) and enjoys extended error-checking capabilities. In addition, MEFISTO is based on a simpler experimental set up as it is based on a collinear arrangement.

The MEFISTO methodology relies on a Fourier analysis after obtaining a spectrally resolved interferometric correlation trace. In particular, we have described the general theory for blind-MEFISTO, pointing out the use of the different terms obtained during the theoretical development. Practical considerations have been discussed. These include spectral calibration, resolution (temporal and spectral), maximum bandwidth, sign indeterminations, possible problems with pulses whose spectra present abrupt changes, error checking procedures, and effect of noise. To manage these experimental limitations, a set of tools have been developed and explained. The analysis has also demonstrated that the method is not affected by some of the non-trivial ambiguities that are present in other techniques [14]. The particular properties of the degenerate case have also been analyzed to highlight additional error-checking capabilities and self-calibrating techniques that help avoids errors entering the pulse measurement.

Finally, we want to stress that this is, to the best of our knowledge, the first analytical methodology capable to simultaneously characterize two unknown ultrashort laser pulses.