Frequency-resolved optical gating

<h2 id="theory">Theory</h2>

<p>FROG and autocorrelation share the idea of combining a pulse with itself in a nonlinear medium.  Since a nonlinear medium will only produce the desired signal when both pulses are present at the same time (i.e. “optical gating”), varying the delay between the pulse copies and measuring the signal at each delay gives a vague estimate of the pulse length.  Autocorrelators measure a pulse by measuring the intensity of the nonlinear signal field.  Estimating the pulse length requires assuming a pulse shape, and the phase of the pulse <a href="/facts/Electric_field/tt4MooBs">electric field</a> cannot be measured at all.  FROG extends this idea by measuring the spectrum of the signal at each delay (hence “frequency-resolved”), instead of just the intensity.  This measurement creates a <a href="/facts/Spectrogram/Jt8zHI0H">spectrogram</a> of the pulse, which can be used to determine the complex electric field as a function of time or frequency as long as the nonlinearity of the medium is known.
</p><p>The FROG spectrogram (usually called a FROG trace) is a graph of intensity as a function of frequency 
  
    
      
        ω
      
    
    {\displaystyle \omega }
  
 and delay 
  
    
      
        τ
      
    
    {\displaystyle \tau }
  
.  The signal field from the nonlinear interaction is easier to express in the <a href="/facts/Time_domain/1w6mP2MU">time domain</a>, however, so the typical expression for the FROG trace includes a <a href="/facts/Fourier_transform/XXhKJtc8">Fourier transform</a>.
</p>

I
          
            FROG
          
        
        (
        ω
        ,
        τ
        )
        =
        
          
            |
            
              
                E
                
                  sig
                
              
              (
              ω
              ,
              τ
              )
            
            |
          
          
            2
          
        
        =
        
          
            |
            
              F
              T
              [
              
                E
                
                  sig
                
              
              (
              t
              ,
              τ
              )
              ]
            
            |
          
          
            2
          
        
        =
        
          
            |
            
              
                ∫
                
                  −
                  ∞
                
                
                  ∞
                
              
              
                E
                
                  s
                  i
                  g
                
              
              (
              t
              ,
              τ
              )
              
                e
                
                  −
                  i
                  ω
                  t
                
              
              
              d
              t
            
            |
          
          
            2
          
        
        .
      
    
    {\displaystyle I_{\text{FROG}}(\omega ,\tau )=\left|E_{\text{sig}}(\omega ,\tau )\right|^{2}=\left|FT[E_{\text{sig}}(t,\tau )]\right|^{2}=\left|\int _{-\infty }^{\infty }E_{sig}(t,\tau )e^{-i\omega t}\,dt\right|^{2}.}

<p>The nonlinear signal field 
  
    
      
        
          E
          
            sig
          
        
        (
        t
        ,
        τ
        )
      
    
    {\displaystyle E_{\text{sig}}(t,\tau )}
  
 depends on the original pulse, 
  
    
      
        E
        (
        t
        )
      
    
    {\displaystyle E(t)}
  
, and the nonlinear process used, which can almost always be expressed as 
  
    
      
        
          E
          
            gate
          
        
        (
        t
        −
        τ
        )
      
    
    {\displaystyle E_{\text{gate}}(t-\tau )}
  
, such that 
  
    
      
        
          E
          
            sig
          
        
        (
        t
        ,
        τ
        )
        =
        E
        (
        t
        )
        
          E
          
            gate
          
        
        (
        t
        −
        τ
        )
      
    
    {\displaystyle E_{\text{sig}}(t,\tau )=E(t)E_{\text{gate}}(t-\tau )}
  
.  The most common nonlinearity is <a href="/facts/Second_harmonic_generation/NzWSrT3J">second harmonic generation</a>, where 
  
    
      
        
          E
          
            gate
          
        
        (
        t
        −
        τ
        )
        =
        E
        (
        t
        −
        τ
        )
      
    
    {\displaystyle E_{\text{gate}}(t-\tau )=E(t-\tau )}
  
.  The expression for the trace in terms of the pulse field is then:
</p>

I
          
            SHG FROG
          
        
        (
        ω
        ,
        τ
        )
        =
        
          
            |
            
              
                ∫
                
                  −
                  ∞
                
                
                  ∞
                
              
              E
              (
              t
              )
              E
              (
              t
              −
              τ
              )
              
                e
                
                  −
                  i
                  ω
                  t
                
              
              
              d
              t
            
            |
          
          
            2
          
        
        .
      
    
    {\displaystyle I_{\text{SHG FROG}}(\omega ,\tau )=\left|\int _{-\infty }^{\infty }E(t)E(t-\tau )e^{-i\omega t}\,dt\right|^{2}.}

<p>There are many possible variations on this basic setup.  If a well-known reference pulse is available, then it may be used as a gating pulse instead of a copy of the unknown pulse.  This is referred to as <a href="/facts/Cross-correlation/Fs1hY3br">cross-correlation</a> FROG or XFROG.  In addition, other non-linear effects besides second harmonic generation may be used, such as third harmonic generation (THG) or polarization gating (PG).  These changes will affect the expression for 
  
    
      
        
          E
          
            gate
          
        
        (
        t
        −
        τ
        )
      
    
    {\displaystyle E_{\text{gate}}(t-\tau )}
  
.
</p>
<h2 id="experiment">Experiment</h2>
<p>In a typical multi-shot FROG setup, the unknown pulse is split into two copies with a beamsplitter.  One copy is delayed by a known amount relative to the other.  Both pulses are focused to the same point in a nonlinear medium, and the spectrum of the nonlinear signal is measured with a spectrometer.  This process is repeated for many delay points.
</p><p>A FROG measurement can be performed on a single shot with some minor adjustments.  The two pulse copies are crossed at an angle and focused to a line instead of a point.  This creates a varying delay between the two pulses along the line focus.  In this configuration, it is common to use a home-made spectrometer, consisting of a <a href="/facts/Diffraction_grating/TUcyt7Do">diffraction grating</a> and a camera, to capture the measurement.
</p>
<h2 id="retrieval-algorithm">Retrieval algorithm</h2>
<p>Although it is theoretically somewhat complex, the method of generalized projections has proven to be an extremely reliable method for retrieving pulses from FROG traces.  Unfortunately, its sophistication is the source of some misunderstanding and mistrust from scientists in the optics community.  Hence, this section will attempt to give some insight into the basic philosophy and implementation of the method, if not its detailed workings.
</p><p>First, imagine a space that contains all possible signal electric fields.  For a given measurement, there is a set of these fields that will satisfy the measured FROG trace.  We refer to these fields as satisfying the data constraint.  There is another set that consists of the signal fields that can be expressed using the form for the nonlinear interaction used in the measurement.  For <a href="/facts/Second-harmonic_generation/NzWSrT3J">second-harmonic generation</a> (SHG), this is the set of fields that can be expressed in the form 
  
    
      
        
          E
          
            sig
          
        
        (
        t
        ,
        τ
        )
        =
        E
        (
        t
        )
        E
        (
        t
        −
        τ
        )
      
    
    {\displaystyle E_{\text{sig}}(t,\tau )=E(t)E(t-\tau )}
  
.  This is referred to as satisfying the mathematical form constraint.
</p><p>These two sets intersect at exactly one point.  There is only one possible signal field that both has the correct intensity to match the data trace and fits the mathematical form dictated by the nonlinear interaction.  To find that point, which will give the pulse we are trying to measure, generalized projections is used.  The generalized projections algorithm operates in this electric field space.  At each step, we find the closest point to the current guess point that will satisfy the constraint for the other set.  That is, the current guess is “projected” onto the other set.  This closest point becomes the new current guess, and the closest point on the first set is found.  By alternating between projecting onto the mathematical constraint set and projecting onto the data constraint set, we eventually end up at the solution.
</p><p>Projecting onto the data constraint set is simple.  To be in that set, the magnitude squared of the signal field has to match the intensity measured by the trace.  The signal field 
  
    
      
        
          E
          
            sig
          
        
        (
        t
        ,
        τ
        )
      
    
    {\displaystyle E_{\text{sig}}(t,\tau )}
  
 is Fourier-transformed to 
  
    
      
        
          E
          
            sig
          
        
        (
        ω
        ,
        τ
        )
      
    
    {\displaystyle E_{\text{sig}}(\omega ,\tau )}
  
. The closest point in the data constraint set is found by replacing the magnitude of 
  
    
      
        
          E
          
            sig
          
        
        (
        ω
        ,
        τ
        )
      
    
    {\displaystyle E_{\text{sig}}(\omega ,\tau )}
  
 by the magnitude of the data, leaving the phase of 
  
    
      
        
          E
          
            sig
          
        
        (
        ω
        ,
        τ
        )
      
    
    {\displaystyle E_{\text{sig}}(\omega ,\tau )}
  
 intact.
</p><p>Projecting onto the mathematical constraint set is not simple.  Unlike the data constraint, there is not an easy way to tell which point in the mathematical constraint set is closest.  A general expression for the distance between the current point and any point in the mathematical constraint set is created, and then that expression is minimized by taking the gradient of that distance with respect the current field guess.  This process is discussed in more detail in <a href="http://frog.gatech.edu/Pubs/DeLong-GenProj-OptLett19-1994.pdf">this paper</a>.
</p><p>This cycle is repeated until the error between the signal guess and the data constraint (after applying the mathematical constraint) reaches some target minimum value.  
  
    
      
        E
        (
        t
        )
      
    
    {\displaystyle E(t)}
  
 can be found by simply integrating 
  
    
      
        
          E
          
            sig
          
        
        (
        t
        ,
        τ
        )
      
    
    {\displaystyle E_{\text{sig}}(t,\tau )}
  
 with respect to delay 
  
    
      
        τ
      
    
    {\displaystyle \tau }
  
.  A second FROG trace is usually constructed mathematically from the solution and compared with the original measurement.
</p>
<h2 id="measurement-confirmation">Measurement confirmation</h2>
<p>One important feature of a FROG measurement is that many more data points are collected than are strictly necessary to find the pulse electric field.  For example, say that the measured trace consists of 128 points in the delay direction and 128 points in the frequency direction.  There are 128×128 total points in the trace.  Using these points, an electric field is retrieved that has 2×128 points (128 for magnitude and another 128 for the phase).  This is a massively <a href="/facts/Overdetermined_system/DMWk4pQR">overdetermined system</a>,<a class="footnote-ref" id="fnref:1" href="#fn:1"><sup>1</sup></a> meaning that the number of equations is much larger than the number of unknowns.  Thus the importance of each individual data point being absolutely correct is greatly reduced.  This is very helpful for real-world measurements that can be affected by detector noise and systematic errors.  Noise is extremely unlikely to affect the measured trace in a way that could be confused with a physical phenomenon in the pulse.  The FROG algorithm tends to “see through” these effects due to the amount of extra information available and the use of a mathematical form constraint in finding a solution.  This means that the error between an experimental FROG trace and a retrieved FROG trace is rarely zero, although it should be quite small for traces without systematic errors.
</p><p>Consequently, significant differences between measured and retrieved FROG traces should be investigated.  The experimental setup may be misaligned, or there may be significant spatio-temporal distortions in the pulse.  If the measurement averages over several or many pulses, then those pulses may vary significantly from each other.
</p>
<h2 id="see-also">See also</h2>
<h3>FROG techniques</h3>
<ul><li><a href="/facts/Grating-eliminated_no-nonsense_observation_of_ultrafast_incident_laser_light_e-fields/dPoKEox8">Grating-eliminated no-nonsense observation of ultrafast incident laser light e-fields</a> (GRENOUILLE), a simplified version of FROG</li>
<li><a href="/facts/Double-Blind_FROG/0TFxy3lw">Double-Blind FROG</a>, for measuring two pulses simultaneously</li></ul>
<h3>Competing techniques</h3>
<ul><li><a href="/facts/Optical_autocorrelation/k6My8cGV">Optical autocorrelation</a>, in its intensity or fringe-resolved (interferometric) version</li>
<li><a href="/facts/Spectral_phase_interferometry_for_direct_electric-field_reconstruction/DRcHIjyU">Spectral phase interferometry for direct electric-field reconstruction</a> (SPIDER)</li>
<li><a href="/facts/Multiphoton_intrapulse_interference_phase_scan/UL11wH0f">Multiphoton intrapulse interference phase scan</a> (MIIPS), a method to characterize and manipulate the ultrashort pulse.</li>
<li>Frequency-resolved electro-absorption gating (FREAG)</li></ul>

<ul><li>Rick Trebino (2002). <i>Frequency-Resolved Optical Gating: The Measurement of Ultrashort Laser Pulses</i>. Springer. <a href="/facts/ISBN_(identifier)/15AdSPa9">ISBN</a> 1-4020-7066-7.</li>
<li>R. Trebino, K. W. DeLong, D. N. Fittinghoff, J. N. Sweetser, M. A. Krumbügel, and D. J. Kane, "<a href="http://frog.gatech.edu/Pubs/Trebino-FROGreview-RSI-1997.pdf">Measuring Ultrashort Laser Pulses in the Time-Frequency Domain Using Frequency-Resolved Optical Gating</a>," <i><a href="/facts/Review_of_Scientific_Instruments/Lwmchn4E">Review of Scientific Instruments</a></i> 68, 3277-3295 (1997).</li></ul>
<h2 id="external-links">External links</h2>
<ul><li><a href="http://frog.gatech.edu">FROG Page by Rick Trebino</a> (co-inventor of FROG)</li></ul>

<h2 id="references">References</h2>

<ol>
<li id="fn:1"><p>Kolesnichenko, Pavel; Zigmantas, Donatas (2023). "Neural-network-powered pulse reconstruction from one-dimensional interferometric correlation traces". Optics Express. 31 (7): 11806–11819. arXiv:2111.01014. doi:10.1364/OE.479638. <a href="/wiki/ArXiv_(identifier)" target="_blank">/wiki/ArXiv_(identifier)</a> <a href="#fnref:1" class="footnote-back-ref">↩</a></p></li>
</ol>

Frequency-resolved optical gating open-in-new

Frequency-resolved optical gating