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

Rapid programmable pulse shaping of femtosecond pulses at the MHz repetition rate

Open Access Open Access

Abstract

We have shown experimentally the successful engineering of femtosecond pulse shaping at a 76 MHz repetition rate input pulse with an acousto-optic modulator (AOM). High repetition rate (HRR) femtosecond laser pulse shaping using an AOM in the Fourier plane was incomprehensible because of its intrinsic 100 kHz acoustic update limit. We demonstrate an effective way of pulse selection and a calibration routine (Fourier shift theorem), which enables generation of ∼10 MHz shaped output pulses from the HRR input pulse train. We have generated a temporally shifted rectangular shaped pulse profile by applying modulation on both the phase and amplitude of the ‘sinc’ RF modulation function.

© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Numerous laser spectroscopic applications depend on amplitude, phase, and frequency of the applied electromagnetic field. In a number of diverse research areas, such as quantum coherent control, imaging, optical communications [13], faithful generation of laser pulses with control in phase and amplitude is of paramount importance. Theoretically, it has been demonstrated that shaped laser pulses can selectively break bonds [4,5], populate selective high vibrational energy levels in anharmonic ladder [6,7], can manipulate curve-crossings between states [8,9], and can suppress intramolecular vibrational relaxation [10]. These studies mainly prompted the beginning of femtosecond optical pulse shaping. Since its first demonstration, the advancement and reliability of femtosecond pulse shaping to generate arbitrarily shaped (amplitude/phase) optical pulses has been used in several complex experiments. Tracing of electronic coherence via multi-dimensional spectroscopy [1113], coherent control of vibrational and/or electronic excitations to manipulate chemical reactions [2,1416], coherent control of dark and bright states in atomic and more complex quantum systems [1718] are few of the well-known experiments performed using pulse shaping technology. In recent times the usefulness of pulse shaping setup has been greatly enhanced by its use in multidimensional electronic spectrometers [11,19] where pulse shaping is commonly used to create two time-delayed pulses to excite the sample. Inherent phase stability of the time delayed pulses created in the pulse shaper simplifies the multi-dimensional experiments.

Optical pulses are defined by their amplitude, phase, and polarization variation with time (and/or frequency). Ultrafast pulse shaping is a powerful tool to precisely manipulate electric field into desired temporal and spectral forms [20]. Fourier spectral filtering in a 4-F geometrical setup [21] is the most commonly used approach [22]. As linear time-invariant filters, liquid crystal modulator (LCM), and acousto-optic modulator (AOM) are extensively used to shape pulses by controlled management of phase and intensity of the laser pulses in Fourier plane of the 4-F setup [23]. Beside 4-F pulse shaper, there are other pulse shaping technologies, which are routinely used for femtosecond pulse shaping. Acousto-optic programmable dispersive filter (AOPDF) [2,24,25] is used for ultrafast pulse shaping, which is based on a collinear acousto-optic interaction that maximizes the interaction length. In the AOPDF, the acoustic wave creates a longitudinal transient grating which maximizes this interaction length. In this technology, an acousto-optic crystal acts as a highly birefringent material for both the acoustic as well as for the optical waves. In the crystal, grating generated by the acoustic wave shapes the optical pulse by diffracting the ordinary incoming optical wave to the diffracted extraordinary optical wave. Thus, optical pulse shaping both in phase and amplitude is achieved through anisotropic interaction. For this reason, the AOPDF does not have to be positioned in the Fourier plane of a 4-F dispersive line, and the technique involves a long material interaction length, either a chirped pulse input or a compressor after the AOPDF is necessary for dispersion compensation. Digital micro-mirror device (DMD) based pulse shaper is also used for ultrafast pulse shaping [26,27]. In this technology, the frequency spectrum of the input pulsed laser is first spread horizontally, and then mapped to a thin strip on the DMD programmed with phase modulation patterns.

Among these approaches, AOM has advantages over LCM in the form of continuous modulation of intensity and phase with the ability of rapid waveform update rate. The AOM pulse shaper acts as diffracting modulator, which generates background free shaped optical pulses [28]. Above mentioned abilities of AOM gives it an edge over LCM as a modulator, when modulation of shot to shot pulse is required (i.e., molecular dynamics measurements). Driving AOM with tailored radiofrequency (RF) generates the acoustic waveform in the AOM, which acts as traveling wave grating. This traveling wave grating acts as a spatial filter and transforms the shape and phase of the input laser pulse [29]. Due to this traveling grating nature of the spatial filter, the refresh rate of a typical AOM is ∼100 kHz (depends on the length of the crystal). However, to function as a time-invariant filter (traveling acoustic waveform in the AOM), the repetition rate of the input laser has to be lower than the refresh-rate mentioned above. The calibration of the AOM is also necessary to design the tailored RF input. Previously, calibration was carried out by reducing the repetition rate of the laser, below the acoustic update rate of the AOM, using pulse picker [2831]. This process ensures the concept of time-invariance of the traveling acoustic filter, and for this reason, pulse shaping at MHz repetition rate input pulse using AOM has been incomprehensible. In last few years, it has been shown, in pump-probe experiments and also in 2D electronic spectroscopy that, the technique involving pulses with high repetition rate (100 kHz) yields much faster data acquisition and enhancement in signal to noise (S/N) ratio [32,33] than the ones using low repetition rate laser pulses. Keeping that trend of innovative advancement in mind, we demonstrate here, femtosecond pulse shaping of MHz repetition rate pulses using AOM. In recent years, DMD-enabled ultrafast pulse shaper has been used to achieve continuous tuning of group velocity dispersion modulation at a rate of 2 MHz [26,27] but the efficiency of the shaper is relatively low (2%) because of the binary hologram setup of the DMD. Here we demonstrate a different technique to calibrate MHz AOM pulse-shaper using principles of Shift Theorem of Fourier Transform. We have characterized the shaped output pulses using this calibrated result. We have successfully generated and reproduced shaped pulses with ∼10 MHz repetition rate, keeping all the advantages provided by AOM pulse shaper (effective efficiency ∼12%). These kinds of pulses are quite advantageous for those experiments where the need of shaped pulse repetition rate is much higher (MHz).

2. Experimental setup

Schematic of the experimental setup is shown in Fig. 1, which utilizes a linear spectral filtering technique on the frequency components of the input pulse to generate a desired shaped output pulse. This AOM based pulse shaper consists of a holographic grating pair (1200 grooves/mm), a pair of concave mirrors (175 mm focal length), arranged in 4-F configuration, to create zero dispersion line on the focal plane (Fourier plane). The AOM (∼4.1 cm TeO2 crystal) is then placed in the Fourier plane as a linear time-invariant Fourier filter, which can modulate both amplitude and phase simultaneously. The first grating disperses incoming laser pulse to its frequency components. The first-order diffracted beam is thereafter focused to diffraction-limited spots on the AOM by the first curved mirror (this forward spatial Fourier transform converts angular dispersion (ω) into spatial dispersion on the Fourier plane). The second mirror converts spatial dispersion into angular dispersion (inverse Fourier transform), which is then re-collimated by the second grating. In the absence of AOM, the pulse shaper operates as a unit magnification device.

 figure: Fig. 1.

Fig. 1. Schematic layout of the optical set-up contains zero dispersion 4-F geometry with AOM in the Fourier plane. The generated RF wave is converted to an acoustic wave by a piezoelectric transducer to the AOM. This traveling acoustic wave inside AOM acts as a diffractive grating mask. In our particular method, only a fraction of pulses is diffracted. The need and usefulness of this fractionalized diffraction are discussed later.

Download Full Size | PDF

We have used an oscillator (MIRA 900F, Coherent Inc.) as our laser pulse source, which produces bandwidth limited ∼150fs full width at half maxima (FWHM) pulses, centered at 780 nm, with a repetition rate of 76 MHz. The tailored RF pulses are generated by an arbitrary function generator (LeCroy, LW420A) with a resolution of 400 Msamples/s. We have used the quadrature circuit [34] for mixing shaped RF with a carrier RF (200 MHz). This synthesized RF has been used as input to the transducer of the AOM. The output optical pulse is then cross-correlated with the unmodulated zero-order pulse, which is the reflected beam from the first grating. A thin BBO (1 mm) crystal is used to generate a cross-correlated second harmonic signal, which is detected by a slow photomultiplier tube (1P28, Hamamatsu). The signal is then fed to a lock-in amplifier (SRS 830 DSP Lock-in Amplifier) to further amplify and increase the signal to noise by order of 200.

3. Calibration of the pulse shaper

In AOM pulse shaping, first laser pulse (temporal profile) Ein(tl) is diffracted by a grating into individual frequency components (ξinl)). The ξinl) is then modulated in the Fourier plane with a transfer function M(ωl) [2], which is applied through the AOM. The output pulse ξoutl)can be written as (ωl is the optical frequency, here we have used subscript ‘l’ for components of light and subscript ‘RF’ for components of applied radio signal):

$${\xi _{out}}({\omega _l}) = {\xi _{in}}({\omega _l})M({\omega _l})\,$$
The transfer function M(ωl) is directly related to the applied RF pulse (f(tRF)). To predict and/or routinely generate desired pulse shapes, one needs to know this one to one relationship. Previously [29,35], this relationship has been measured via propagating a small radio pulse through AOM with a series of triggered delays, and by measuring the corresponding diffracted spectrum. This also ensures the fact that ‘appropriate synchronization with the laser pulse to RF pulse’ has been achieved. The technique, however, required that the laser repetition rate must be lower than the acoustic update rate (100 kHz) of AOM. In this framework, the traveling acoustic wave acts as a linear time-invariant filter. In our case, the repetition rate of the laser pulse is much higher (76 MHz). This prevents precise selection of the spectrum from a single pulse, as several light pulses will be deflected even for the smallest RF pulse window. The concept of ‘triggering and synchronization’ between RF pulse and laser pulse, shot to shot, is pointless in case of repetition rate higher than 100 kHz.

3.1 Optimized selection of RF window

Though selection and shaping of individual laser pulse is impossible, MHz repetition rate pulse shaping via AOM is achievable. Here, we call this MHz pulse shaping via AOM, an ‘averaged method of pulse shaping’. This is due to the fact that there is no possible way to deliver a single transfer function to individual laser pulses which are appearing at 76 MHz repetition rate to the AOM. The applied RF pulse creates a travelling grating in the AOM crystal and the neighboring optical pulses interact with different parts of the travelling grating. This can be summarized in the cartoon as shown in Fig. 2.

 figure: Fig. 2.

Fig. 2. Cartoon representation of the selection process involved in the MHz pulse shaping using AOM. The applied RF window creates traveling grating, which deflects several pulses (1-9) in sequence. The fact is that all deflected pulses do not contain the entire spectrum, as they are selected via a traveling grating inside the AOM.

Download Full Size | PDF

Here, the acoustic update rate is much slower than the repetition rate of the laser pulse, leading to the point that a traveling grating in the AOM interacts and deflects several pulses, (Fig. 2) when laser repetition rate is in MHz range. Not only that, whenever an RF pulse arrives at the AOM crystal, it will always encounter the presence of the laser spectrum. Thus, whenever we try to modulate MHz repetition rate pulse using AOM, the need of previously mentioned method, “appropriate synchronization with the laser pulse to RF pulse”, becomes meaningless. When an RF pulse travels through AOM, it not only deflects several light pulses, but it slices different portions of the spectrum of the selected pulses (Fig. 2). Depending on the width of the RF pulse window, the number of deflected pulses will vary. The pulses 1-3 and 7-9 (Fig. 2) are partially diffracted. But pulses 4-6 are the copy of the original pulse (here ‘copy’ means pulses without any spectral slices, though, of course, in the strictest sense these copies will not be exactly identical to the original input as the dispersion from AOM and grating would have some influence). The deflected pulses are then cross-correlated with the pulses from the original source. Cross-correlation signal is generated as an average effect of all the pulses in correlation, which in this case, is an average measurement. This leads to deformation from the original pulse (here original pulse means output pulse from the shaper when the pulse is only diffracted by 200 MHz carrier radio frequency and not modulated at all by conventional intensity or phase modulation). If the numbers of these copies are increased, the averaged signal will also resemble the original pulse.

An RF pulse segment is created with 10µs duration and with a variable window (window duration of = 5µs, 4µs, 3µs, etc.).Here, a train of RF square pulse of different pulse width (variable RF window) has been generated and each pulse is separated by 10-ΔRF µs (see Fig. 2). It is then applied as an input to the AOM transducer, and the deflected pulses are then cross-correlated with the original pulse and cross-correlation signal is then measured (Fig. 3). This shows us, for our case, a 150fs FWHM Gaussian pulses, an RF segment of 10µs with the 2µs window is good enough to produce an optimized ratio between completely diffracted pulses to those of partially diffracted pulses (which are generating the average signal from the original pulses). Putting it simply, since the applied spectral phase is constant and equals to zero, it does not matter that the RF pulse is sliding through the spectrum for each laser shot. This is because, “constant zero is zero”, no matter at which frequency the input laser pulse is applied. This is the heart of this novel technique, generation of the optimized number of copied pulses.

 figure: Fig. 3.

Fig. 3. Cross-correlation signal from pulses deflected using 10 µs RF pulse segment with variable duration of radio pulse window: (a) When radio pulse window is too narrow (i.e., 0.1 µs and 0.05 µs), deformation in cross-correlation signal is significant as the number of completely diffracted pulses are rather low. (b) The enhanced cross-correlation features show the fact that pulse width decreases as the duration of radio window increases. After reaching 2 µs and beyond, the pulse width does not change and satisfactorily overlaps with the cross-correlation trace of the original pulse (output from the pulse shaper when no modulation is present).

Download Full Size | PDF

We have captured the deflected pulse-train signals, generated by applying different RF windows (Fig. 4) in AOM. From these figures, we can easily comprehend the facts described in cartoon representation in Fig. 2. As the moving grating passes through the AOM, it deflects 76 MHz pulse train to a top-hat like waveform structure. The waveform has a rise time, an equilibrium time and a fall time. In this selection processes, inside the rising and falling part of the waveform, spectrally sliced pulses (Fig. 2— slices 1-3 and 7-9) reside, whereas the pulses within the equilibrium portion are completely diffracted pulses (slices 4-6). The total waveform created by a 2µs windowed RF pulse is 2.614µs, and the equilibrium portion is ∼1.69µs. This means that the percent of the completely diffracted pulses are 65% in this average selection and measurement. The copied pulses for the 1µs case are only 35% of the total deflected pulses and are ∼1% for the 0.5µs case. However, it is to be noted that, in this kind of selection procedure also, there is a reduction in the effective repetition rate of the shaped laser pulses. The use of, a 2 µs RF pulse window with a 10 µs RF pulse segment duration, means that ∼20% of the pulse from the 76 MHz pulse train are diffracted. Furthermore, only 65% of the diffracted pulses are “completely diffracted” pulses, which are correctly modulated. Combining these reduction factors, the effective repetition rate of the diffracted pulse is ∼ 10 MHz. This is several folds higher than the repetition rate achievable with other previous pulse shaping approaches using AOM. One more point to be noted that the shaped pulses are still separated ∼13 ns apart from each other.

 figure: Fig. 4.

Fig. 4. Measured Photodiode signals of the deflected pulses are presented here. (a) The deflected signal due to different RF window function. (b) 2µs windowed Photodiode signal and some of its zoomed in parts. (c) The envelope of the waveform signals in (a).

Download Full Size | PDF

In our case, it is prominent from the pictorial representation that the spectrum of the input pulse is not spread all over the available aperture of the AOM. This has been done deliberately to minimize the amount of spectrally sliced pulses. The available resolution of the shaper thus has been reduced, but it is the necessary requirement of this technique to work. This method of pulse selection is also the main obstacle to previously mentioned triggering and synchronization procedure to calibrate pulse shaper. In regard to the points discussed above, we can find the available modulation pixel of the AOM. AOM is 4.1 cm in length, and the acoustic velocity in TeO2 crystal is 4.2 mm/µs, which translate into a transit time Δτ∼9.76 µs. In a TeO2 crystal, the minimum spatial feature size is of 21 µm [2,36]. The RF bandwidth of the AOM is 100 MHz, which translates ∼976 available AOM pixels. The effective modulation window of 2µs creates 8.4 mm spatial feature, and the spectral spread of the laser pulse on the AOM crystal is about 5.5 mm. Thus,

$$\lambda = d(\sin \phi - \sin \theta )\,\,\,\,\,\,\,\,\,\textrm{and}\,\,\,\,\,\,\,x^{\prime} = f \times \tan ({\theta _\lambda } \pm {\theta _{\lambda 0}})$$
Here ϕ is the incident angle, θ is diffracted angle from the grating, d is grating spacing, x’ is spatial dispersion on the AOM, f is the focal length of the spherical mirror. The total angle created after diffraction is ∼1.8o which results in a spatial dispersion of 5.5 mm on the AOM. For the 5.5 mm spatial dispersion of the pulse spectrum on the AOM, this provides (5.5 mm)*(1 us/4.2 mm)*(100 MHz) ∼ 130 effective pixels. This is a significant reduction of the resolution of the AOM shaper, but it is necessary for this technique to work.

In our system, the diffraction efficiency of the grating pair is 85%. The efficiency of the curved mirrors is 95%, and the first order deflection efficiency of the AOM is ∼20%. The selection scheme further reduces the average output of the shaper by 13% (∼10 MHz output from 76 MHz input). The overall efficiency of the average power output from the shaper is thus 1.3% of the average input power. However, the peak output does not get reduced by the pulse selection routine, and therefore the effective efficiency of the peak power output is 12.35%.

3.2 Standardization between transfer function and RF pulse

In order to generate an identical output pulse as compared to that of the original pulse, a well-defined RF pulse window (the width of the window depends on the pulse-width of the laser pulse) has to be applied. In this framework, if a modulation (intensity and/or phase) is applied in the windowed part, it will modify the optical pulse at the Fourier plane of the 4-F. The question is, how much the output pulse shape resembles the expected modulation created by the RF, and if it is so, how a desired shaped pulse can be created where the optical pulse is in MHz repetition rate? To formulate the exact transformation between transfer function (M(ωl) and RF signal (f(tRF)), calibration of the pulse shaper is an absolute necessity. As mentioned previously, calibration via ‘triggering and synchronization’ is impossible in this case. Here, this technical barrier has been overcome by using Fourier Shift Theorem principle to calibrate our MHz AOM pulse shaper.

The Fourier shift theorem says that a delay in the time domain corresponds to a linear phase sweep term in the frequency domain. This can be represented as:

$$f(t + \alpha )\buildrel {\textrm{FourierTransform}} \over \longleftrightarrow {e^{2\pi i\omega \alpha }}F(\omega )$$
In practice, delayed pulse generation by applying linear phase sweep in the frequency domain has been generated earlier [36]. There, authors have demonstrated how to generate, ‘Rapid ultrafast tunable delay line’ by using this principle. The authors have also established one-to-one correspondence between a single temporarily shifted laser pulse to the applied RF signal. The broadening in the shifted output pulse is also explained, which emerges from the higher-order effects of the zero-dispersion line [36]. The question was, whether the same principles would hold or not in our case of average selection scheme.

Applying linear phase sweep as a transfer function, we have generated several instances of the 2µs window, in which we have applied sinusoidal radio modulation of frequency 1 MHz, 5 MHz, and 10 MHz, etc. Then cross-correlation signals from the linear phase swept (spectral domain) output pulses were measured (Fig. 5). To our satisfaction, it was found that the peak shift occurs in the measured cross-correlation signal (temporal domain). The broadening of the output pulses has also been observed in the cross-correlation measurements, which is coming from higher order effects. In our 4-F alignment (see Fig. 1), the input optical pulse is non-normal to that of AOM, which invokes an inherent optical path length difference for each frequency component. This optical path-length difference for different frequency components induces dispersion and results in pulse broadening [36]. The interesting fact is that, when we prepared the sinusoidal modulation within in 1µs RF window, it shows the same peak shift effect. Not only that, both 2µs and 1µs window results in the exact same amount of peak shift in corresponding sinusoidal modulation frequency (Fig. 6). This proves another crucial point that, though the modulation experienced by first pulse and the last pulse (selected by the windowed RF) is significantly different, the average effect (cross-correlation measurement) merge to give a single output peak shift. A linear slope shifted in space (experienced by first and last pulses) is still a linear slope throughout, and thus the temporal shift is not affected (only an offset result, i.e., a zero-order spectral phase affecting the carrier-envelope phase that is not detected here). This experimental finding proves another curtail point, which is, the modulation function M(ωl) is a slowly and smoothly varying function in comparison to the spectrally dispersed optical pulse (ξinl)), which requires the condition, that the smallest feature of M(ωl) is larger than the effective Gaussian beam waist (which, for our case 6.3 µm) at the central Fourier plane. This is because each of these pulses is temporarily very close to the other (only ∼13 ns apart), and the modulation window (2µs) is much larger than this pulse to pulse separation.

 figure: Fig. 5.

Fig. 5. Time segmented RF pulses applied to AOM, which creates the delayed optical pulses in the temporal domain. Applied RF pulse is a 2µs window within a 10µs time segment with 1 MHz to 40 MHz oscillation frequencies. In (a), the real part of the RF signal generated using 1 MHz oscillation frequency has been pictorially represented and (b) is the corresponding imaginary part of the sinusoidal modulation, and (c) shows the cross-correlation traces of delayed optical pulses created by different oscillation frequencies.

Download Full Size | PDF

 figure: Fig. 6.

Fig. 6. Peak shift due to different linear phase sweep for both the case of 1 µs and 2 µs window is represented here. Though 2 µs window creates pulses whose average effect resembles that of original pulse and 1 µs window does not do so, the linear sweep results in the same shift for both the cases.

Download Full Size | PDF

This explains why one can time-delay the pulses using a linear spectral phase even when the modulation RF wave has preceded through the crystal from laser shot to laser shot (Fig. 6). In other words, a straight line is translationally invariant. Thus, one can indeed obtain a calibration of the AOM pulse shaper, even with 76 MHz pulses as input. This linear relationship depends on several factors, such as response time of the AOM (depends on the material of AOM), the bandwidth of the light pulse, the repetition rate of the femtosecond laser pulse. To do MHz pulse shaping using AOM, one has to find out the RF window, as well as has to establish a linear relationship between the laser spectra to the RF pulse within the RF window (acoustic waveform in the AOM).

The above-described results are the facts that lead us to correlate the transfer function (M(ωl)) to the applied RF modulation frequency using the mathematical methods, described in Fourier Shift theorem. Using Matlab script, we have found out the exact relationship by back calculating the experimental data. First, we transform spectral and temporal representations of the pulse to the corresponding pixel representations (henceforth pixelation, which we denote as the dimensionless discretization of the relative variables for the numerical procedure). The smallest measured temporal resolution is dtl ( = 6.666 fs). Temporal domain pulse Ein(tl) is transformed to discretized version, Ein(pt) via pixel transform (PT), by replacing tl to tl/dtl. This generates Nt number of temporal pixel points, with pt (=tl/dtl) being the temporal pixel variable. Similarly, spectral domain pulse ξinl) is also remodeled to ξin(pω), with Nω spectral points and pω (=ωl/dωl) being spectral pixel variable (Fig. 7).

 figure: Fig. 7.

Fig. 7. Pixelated representation of spectral domain pulse and radio pulse: The goal is to establish one to one correlation between the two-pixel forms.

Download Full Size | PDF

The input RF function, f(tRF) is also transformed into f(pRF) via PT, and it contains NRF radio pixel points. In the relative pixelated representation, the number of pixels is set to equal, that is, Nt=Nω=NRF=N. We know that the transfer function, M(ωl), is directly related to f(tRF). The goal is to establish a relation between pixelated transfer function (M(pω)) and pixelated radio pulse f(pRF). Thus, pω is a scaled version of pRF (${p_\omega } = c \times {p_{RF}}$, where c is constant). We would like to extract this constant scaling value.

In our experiments, we have measured the amount of peak shift in temporal domain caused by application of smoothly varying phase sweep in the frequency domain. Fourier Shift Theorem says that linear phase sweep in frequency corresponds to peak shift in time and vice versa. In this scenario, if the spectral pixel imposes linear phase sweep, it will generate a peak shift of 1 in a temporal pixel. In the experiment, we have measured the shift in peak position to that of the zero-delay position. A resultant peak shift of ‘α’ in the temporal profile will generate a peak shift of ‘a’ (=α/dtl) in pixel representation (Fig. 8). This means an applied radio-pulse, as a whole, has delivered a 2πa linear phase sweep to the spectral pixel. As the radio pulse is generated by sinusoidal radio frequency modulation, we can easily characterize it and can draw a conclusive correlation between radio pulse and transfer function in pixelated representation (Fig. 8).

 figure: Fig. 8.

Fig. 8. Flow chart to find out the exact correspondence between pixelated transfer function (M(pω)) and pixelated radio pulse f(p­RF) from experimental measurements of peak shift via application of different sinusoidal modulation. Bold lines are main flow chart, and dotted lines are for comparison and checking.

Download Full Size | PDF

4. Comparison among desired and measured shaped output pulses

4.1 Intensity modulation

Taking advantage of the mathematical relationship mentioned in the earlier section, we have designed some arbitrary RF pulses for intensity modulation and predicted cross-correlation of the shaped output. RF pulses are generated through a simple combination of ‘sin’ functions with 10µs time-segment and 1µs square window. A series of these linearly combined RF pulses are used to intensity modulate a Gaussian pulse. Two of them are represented here in Fig. 9(a) and 9(b) (using equations 2a and 2b) with appropriately tailored RF pulses:

$$\begin{aligned}{f_1}({t_{RF}}) &= \sin (2\pi \nu _{RF}^1{t_{RF}}) + 0.35\sin (2\pi \nu _{RF}^2{t_{RF}}) + 0.25\sin (2\pi \nu _{RF}^4{t_{RF}})\\ &\quad + 0.1\sin (2\pi \nu _{RF}^6{t_{RF}}) \end{aligned}$$
$$\begin{aligned}&{f_2}({t_{RF}}) = \sin (2\pi \nu _{RF}^2{t_{RF}}) + 0.2\sin (2\pi \nu _{RF}^4{t_{RF}}) + 0.4\sin (2\pi \nu _{RF}^6{t_{RF}})\\ & \textrm{where}\,\,\,\,\,\nu _{RF}^1 = 1\textrm{MHz,}\ \nu _{RF}^2 = 2\textrm{MHz},\nu _{RF}^4 = 4\textrm{MHz},\nu _{RF}^6 = 6\textrm{MHz} \end{aligned}$$
Cross-correlations of output pulses match reliably with the predicted traces, calculated using experimentally found α value. And the modulation, represented in Eq. (2a), changes to
$${f_1}({p_{RF}}) = \sin (2\pi .1.{p_{RF}}) + 0.35\sin (2\pi .2.{p_{RF}}) + 0.25\sin (2\pi .4.{p_{RF}})\, + 0.1\sin (2\pi .6.{p_{RF}})$$
$$\downarrow $$
$${f_1}({p_\omega }) = \sin (2\pi .1.\frac{a}{N}{p_\omega }) + 0.35\sin (2\pi .2.\frac{a}{N}{p_\omega }) + 0.25\sin (2\pi .4.\frac{a}{N}{p_\omega })\,\, + 0.1\sin (2\pi .6.\frac{a}{N}{p_\omega })$$
Presence of many ‘sin’ modulations divides the input pulse into many smaller segments with their respective intensities proportional to their RF counterparts. Since 1µs RF window does not generate enough number of copied pulses to give cross-correlation trace resembling the original pulse, there is an obvious broadening of the fragmented sections. This effect is observed when we tried to create temporal square pulses by using RF ‘sinc’ function with 10µs time segment and 2µs window.

 figure: Fig. 9.

Fig. 9. Intensity modulation through a combination of ‘sin’ functions and cross-correlation traces of shaped pulse is represented. Input RF pulse corresponding to equations (2a, 2b) are depicted in (a) and (b). Experimental and expected values corresponding to RF (a) and (b) is presented in (c) and (d). Using the previously mentioned method, we get the expected cross-correlation, which is in good agreement with experiment.

Download Full Size | PDF

4.2 Intensity and phase modulation: square pulse

Due to a balancing problem in our RF-circuity, there is a constant DC leakage of the RF carrier with the tailored RF pulse, which in combination with the MHz repetition rate laser pulse creates a constant background. Though contribution from this is negligible, in the ‘sin’ amplitude modulation case, the effect is paramount when we tried to create temporal square pulses by using RF ‘sinc’ function with 10µs time segment and 2µs window.

The background signal from unmodulated pulse is comparable to the signal from the original square pulse. Still, by including a fixed background, the cross-correlation trace of the square pulse again matches the expected shape as calculated by using the pixel transform technique (Fig. 10(a)).

$$\begin{aligned}& {f_{{\textrm{sinc}} }}({t_{RF}}) = \frac{{\sin (\pi \nu _{RF}^{10}{t_{RF}})}}{{\pi \nu _{RF}^{10}{t_{RF}}}}.\exp \left[{i\left({2\pi \nu_{RF}^{10}{t_{RF}} - 20\pi \nu_{RF}^{10}t_{RF}^3} \right)} \right]\,\\& \textrm{where}\,\,\,\,\nu _{RF}^{10} = 10\textrm{MHz} \end{aligned}$$
To confirm that a proper square pulse is generated via ‘sinc’ modulation, we repeated our experiments using an additional linear sweep to produce a time-shifted square pulse from the background. We found out that the time-shifted square pulse to be slightly deformed (Fig. 10(b), black dot trace). This is due to the deformation caused by the dispersion in the AOM which is considerable in case of square pulse generation as the RF ‘sinc’ modulation has much more complexity than the previously applied modulations (combination of ‘sin’ modulations).

 figure: Fig. 10.

Fig. 10. Intensity modulation through ‘sinc’ RF function generates a square pulse in the temporal domain. Experimental cross-correlation of the generated square pulse is represented: (a) Cross-correlation of the output from only ‘sinc’ RF modulation; (b) same ‘sinc’ RF with a linear sweep and pre-compensated ‘sinc’ RF modulation with a linear sweep.

Download Full Size | PDF

Interestingly, however, both the deformed time-shifted square pulses, generated from positive and negative linear phase swept ‘sinc’ modulation, creates cross-correlation traces that are mirror images of each other with respect to the zero position. This leads to the key information about this deformation, which is, that the linear sweep adds a higher order phase in the traveling acoustic wave in the AOM. This inherent higher order phase will invoke chirp to the input pulse. We can compensate for this chirp by using a pre-compensated RF pulse, which adds a positive or negative first order chirp based on the frequency sweep to the input pulse. The imaginary and real part of the Eq. (3) is presented in Fig. 11(a). Surprisingly, only a first-order correction in the input pulse compensates the deformation effectively, and the cross-correlation trace matches with the expected trace (Fig. 11(b)). This is done by adding only a cubic phase in the tailored RF pulse as shown in Eq. (3) above.

 figure: Fig. 11.

Fig. 11. (a) The plot of the imaginary and real part of the pre-compensated RF signal that is generated by adding a third-order phase with the linear frequency sweep. Experimental cross-correlation signal and expected square pulse (using αRF) are presented in (b) that are quite good agreement with each other.

Download Full Size | PDF

5. Conclusion

In conclusion, we have used AOM pulse-shaping setup as our experimental device to successfully alter the temporal shape of the pulse and get the desired outcome by using high repetition rate laser (76 MHz) as input. We have achieved much higher repetition rate pulse shaping (∼10 MHz) via AOM pulse shaper technique. The essential part of this pulse shaping technique is the establishment of the square windows RF signal by ensuring the linear peak shift (linear relationship between the RF pulse and spectrum of the light pulse). We have taken full advantage of the principle of Shift-Theorem of Fourier transform to calibrate our high repetition rate laser pulse-shaping setup in a new way (by transforming into a pixelated form). Using this calibration, we have successfully demonstrated several shaped pulses applying intensity modulation, which are in good agreement with predicted shapes. Finally, we have successfully prepared a shifted square pulse that is free of deformation using a pre-compensation scheme of simultaneous intensity and phase modulation. Using our pixelated calibration scheme, we have characterized all the output shaped pulses. In this present first demonstration experiments for high-repetition experiments with AOM pulse shaping, a reduced effective resolution AOM pulse shaper was used that somewhat limits the possible freedom and capabilities of the AOM pulse-shaping. The limitation comes in the RF modulation, where the modulation function has to be smooth and slowly varying, compared to the time gap between laser shot to laser shot.

This kind of pulse shaping is beneficial for measuring optical free induction decay, where it is advantageous to have the full MHz shaped pulse train produced. The other point is that the shaped pulse sequences are created using an acousto-optic pulse shaper, which means that the experimental setup is the same regardless of whether a sequence of two or of hundreds of pulses is created. This technique is especially useful in samples where phase-matching is not possible such as single-molecule spectroscopy or in samples where it is impossible to detect the transmission, for example in highly scattering samples. To our success, we have shown the high fidelity of these shaped pulses in terms of reproducibility, i.e., in one case of not generating the shaped pulse, we recover a bunch of deformed shapes while in the other case when we do shape the pulse, we detect the shape faithfully.

Funding

Science and Engineering Research Board (SERB) (Individual Intramural Funds); Indian Space Research Organisation (ISRO) (STC).

Acknowledgment

SD and SNB are thankful for Senior Research Fellowship of the Council of Scientific and Industrial Research, India. All the authors acknowledge S. Goswami for extensive language editing of this article.

References

1. W. S. Warren, H. Rabitz, and M. Dahleh, “Coherent Control of Quantum Dynamics: The Dream is Alive,” Science 259(5101), 1581–1589 (1993). [CrossRef]  

2. D. Goswami, “Optical pulse shaping approaches to coherent control,” Phys. Rep. 374(6), 385–481 (2003). [CrossRef]  

3. H. Kawashima, M. M. Wefers, and K. A. Nelson, “Femtosecond pulse shaping, multiple-pulse spectroscopy, and optical control,” Annu. Rev. Phys. Chem. 46(1), 627–656 (1995). [CrossRef]  

4. S. Shi and H. Rabitz, “Optimal control of selective vibrational excitation of harmonic linear chain molecules: analytic solution and restricted forms for the optimal fields,” J. Chem. Phys. 92(5), 2927–2937 (1990). [CrossRef]  

5. S. Rice, B. Amstrup, R. Carlson, and A. Matro, “The use of pulse shaping to control the photodissociation of a diatomic molecule: preventing the best from being the enemy of the good,” J. Chem. Phys. 95(21), 8019–8027 (1991). [CrossRef]  

6. J. S. Melinger, S. Gandhi, A. Hariharan, J. X. Tull, and W. S. Warren, “Generation of narrowband inversion with broadband laser pulses,” Phys. Rev. Lett. 68(13), 2000–2003 (1992). [CrossRef]  

7. S. Chelkowski, A. Bandrauk, and P. B. Corkum, “Efficient molecular dissociation by a chirped ultrashort infrared laser pulse,” Phys. Rev. Lett. 65(19), 2355–2358 (1990). [CrossRef]  

8. P. Gross, D. Neuhauser, and H. Rabitz, “Optimal control of curve-crossing systems,” J. Chem. Phys. 96(4), 2834–2845 (1992). [CrossRef]  

9. E. E. Aubanel, J. Gauthier, and A. D. Bandrauk, “Molecular stabilization and angular distribution in photodissociation of H2+ in intense laser fields,” Phys. Rev. A 48(3), 2145–2152 (1993). [CrossRef]  

10. D. Goswami and W. S. Warren, “Control of chemical dynamics by restricting intramolecular vibrational relaxation,” J. Chem. Phys. 99(6), 4509–4517 (1993). [CrossRef]  

11. S. H. Shim and M. T. Zanni, “How to turn your pump–probe instrument into a multidimensional spectrometer: 2D IR and Vis spectroscopies via pulse shaping,” Phys. Chem. Chem. Phys. 11(5), 748–761 (2009). [CrossRef]  

12. K. Gundogdu, K. W. Stone, D. B. Turner, and K. A. Nelson, “Multidimensional coherent spectroscopy made easy,” Chem. Phys. 341(1-3), 89–94 (2007). [CrossRef]  

13. K. W. Stone, K. Gundogdu, D. B. Turner, X. Q. Li, S. T. Cundiff, and K. A. Nelson, “Two-quantum 2D FT electronic spectroscopy of biexcitons in GaAs quantum wells,” Science 324(5931), 1169–1173 (2009). [CrossRef]  

14. K. Ohmori, “Wave-Packet and coherent control dynamics,” Annu. Rev. Phys. Chem. 60(1), 487–511 (2009). [CrossRef]  

15. M. Dantus and V. V. Lozovoy, “Experimental coherent laser control of physicochemical processes,” Chem. Rev. 104(4), 1813–1860 (2004). [CrossRef]  

16. A. Assion, T. Baumert, M. Bergt, T. Brixner, B. Kiefer, V. Syfried, M. Strehle, and G. Gerber, “Control of Chemical Reactions by Feedback-Optimized Phase-Shaped Femtosecond Laser Pulses,” Science 282(5390), 919–922 (1998). [CrossRef]  

17. D. Meshulach and Y. Silberberg, “Coherent quantum control of two-photon transitions by a femtosecond laser pulse,” Nature 396(6708), 239–242 (1998). [CrossRef]  

18. S. Chen, A. Jaron-Becker, and A. Becker, “Time-dependent analysis of few-photon coherent control schemes,” Phys. Rev. A 82(1), 013414 (2010). [CrossRef]  

19. S. H. Shim, D. B. Strasfeld, and M. T. Zanni, “Generation and characterization of phase and amplitude shaper femtosecond mid-IR pulses,” Opt. Express 14(26), 13120–13130 (2006). [CrossRef]  

20. A. M. Weiner and J. P. Heritage, “Picosecond and femtosecond Fourier pulse shape synthesis,” Rev. Phys. Appl. 22(12), 1619–1628 (1987). [CrossRef]  

21. O. E. Martinez, “Matrix formalism for pulse compressors,” IEEE J. Quantum Electron. 24(12), 2530–2536 (1988). [CrossRef]  

22. A. M. Weiner, “Femtosecond optical pulse shaping and processing,” Prog. Quantum Electron. 19(3), 161–237 (1995). [CrossRef]  

23. A. M. Weiner, “Femtosecond pulse shaping using spatial light modulators,” Rev. Sci. Instrum. 71(5), 1929–1960 (2000). [CrossRef]  

24. F. Verluise, V. Laude, Z. Cheng, C. Spielmann, and P. Tournois, “Amplitude and phase control of ultrashortpulses by use of an acousto-optic programmable dispersive filter: pulse compression and shaping,” Opt. Lett. 25(8), 575–577 (2000). [CrossRef]  

25. S. Coudreau, D. Kaplan, and P. Tournois, “Ultraviolet acousto-optic programmable dispersive filter laser pulse shaping in KDP,” Opt. Lett. 31(12), 1899–1901 (2006). [CrossRef]  

26. C. Gu, Y. Chang, D. Zhang, J. Cheng, and S. Chen, “Femtosecond Laser Pulse Shaping at Megahertz Rate via a Digital Micromirror Device,” Opt. Lett. 40(17), 4018–4021 (2015). [CrossRef]  

27. C. Gu, D. Zhang, Y. Chang, and S. Chen, “Digital Micromirror Device-based Ultrafast Pulse Shaping for Femtosecond Laser,” Opt. Lett. 40(12), 2870–2873 (2015). [CrossRef]  

28. W. Hillegas, J. X. Tull, D. Goswami, D. Strickland, and W. S. Warren, “Femtosecond laser pulse shaping by use of microsecond radio-frequency pulses,” Opt. Lett. 19(10), 737–739 (1994). [CrossRef]  

29. M. R. Fetterman, D. Goswami, D. Keusters, W. Yang, J.-K. Rhee, and W. S. Warren, “Ultrafast pulse shaping: amplification and characterization,” Opt. Express 3(10), 366–375 (1998). [CrossRef]  

30. M. Roth, M. Mehendale, A. Bartelt, and H. Rabitz, “Acousto-optical shaping of ultraviolet femtosecond pulses,” Appl. Phys. B 80(4-5), 441–444 (2005). [CrossRef]  

31. M. A. Dugan, J. X. Tull, and W. S. Warren, “High-resolution acousto-optic shaping of unamplified and amplified femtosecond laser pulses,” J. Opt. Soc. Am. B 14(9), 2348–2358 (1997). [CrossRef]  

32. F. Kanal, S. Keiber, R. Eck, and T. Brixner, “100-kHz shot-to-shot broadband data acquisition for high-repetition-rate pump–probe spectroscopy,” Opt. Express 22(14), 16965–16975 (2014). [CrossRef]  

33. N. M. Kearns, R. D. Mehlenbacher, A. C. Jones, and M. T. Zanni, “Broadband 2D electronic spectrometer using white light and pulse shaping: noise and signal evaluation at 1 and 100 kHz,” Opt. Express 25(7), 7869 (2017). [CrossRef]  

34. J. X. Tull, M. A. Dugan, and W. S. Warren, “High-Resolution, Ultrafast Laser Pulse Shaping and Its Applications,” Adv. Magn. Opt. Reson. 20, 1–65 (1997). [CrossRef]  

35. A. Nag, P. A. Chaphekar, and D. Goswami, “Applying genetic algorithm optimization to a folded geometry acousto-optic modulated spatial pulse shaper,” Rev. Sci. Instrum. 81(1), 013101 (2010). [CrossRef]  

36. W. Yang, D. Keusters, D. Goswami, and W. S. Warren, “Rapid ultrafine-tunable optical delay line at the 1.55-µm wavelength,” Opt. Lett. 23(23), 1843–1845 (1998). [CrossRef]  

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. Schematic layout of the optical set-up contains zero dispersion 4-F geometry with AOM in the Fourier plane. The generated RF wave is converted to an acoustic wave by a piezoelectric transducer to the AOM. This traveling acoustic wave inside AOM acts as a diffractive grating mask. In our particular method, only a fraction of pulses is diffracted. The need and usefulness of this fractionalized diffraction are discussed later.
Fig. 2.
Fig. 2. Cartoon representation of the selection process involved in the MHz pulse shaping using AOM. The applied RF window creates traveling grating, which deflects several pulses (1-9) in sequence. The fact is that all deflected pulses do not contain the entire spectrum, as they are selected via a traveling grating inside the AOM.
Fig. 3.
Fig. 3. Cross-correlation signal from pulses deflected using 10 µs RF pulse segment with variable duration of radio pulse window: (a) When radio pulse window is too narrow (i.e., 0.1 µs and 0.05 µs), deformation in cross-correlation signal is significant as the number of completely diffracted pulses are rather low. (b) The enhanced cross-correlation features show the fact that pulse width decreases as the duration of radio window increases. After reaching 2 µs and beyond, the pulse width does not change and satisfactorily overlaps with the cross-correlation trace of the original pulse (output from the pulse shaper when no modulation is present).
Fig. 4.
Fig. 4. Measured Photodiode signals of the deflected pulses are presented here. (a) The deflected signal due to different RF window function. (b) 2µs windowed Photodiode signal and some of its zoomed in parts. (c) The envelope of the waveform signals in (a).
Fig. 5.
Fig. 5. Time segmented RF pulses applied to AOM, which creates the delayed optical pulses in the temporal domain. Applied RF pulse is a 2µs window within a 10µs time segment with 1 MHz to 40 MHz oscillation frequencies. In (a), the real part of the RF signal generated using 1 MHz oscillation frequency has been pictorially represented and (b) is the corresponding imaginary part of the sinusoidal modulation, and (c) shows the cross-correlation traces of delayed optical pulses created by different oscillation frequencies.
Fig. 6.
Fig. 6. Peak shift due to different linear phase sweep for both the case of 1 µs and 2 µs window is represented here. Though 2 µs window creates pulses whose average effect resembles that of original pulse and 1 µs window does not do so, the linear sweep results in the same shift for both the cases.
Fig. 7.
Fig. 7. Pixelated representation of spectral domain pulse and radio pulse: The goal is to establish one to one correlation between the two-pixel forms.
Fig. 8.
Fig. 8. Flow chart to find out the exact correspondence between pixelated transfer function (M(pω)) and pixelated radio pulse f(p­RF) from experimental measurements of peak shift via application of different sinusoidal modulation. Bold lines are main flow chart, and dotted lines are for comparison and checking.
Fig. 9.
Fig. 9. Intensity modulation through a combination of ‘sin’ functions and cross-correlation traces of shaped pulse is represented. Input RF pulse corresponding to equations (2a, 2b) are depicted in (a) and (b). Experimental and expected values corresponding to RF (a) and (b) is presented in (c) and (d). Using the previously mentioned method, we get the expected cross-correlation, which is in good agreement with experiment.
Fig. 10.
Fig. 10. Intensity modulation through ‘sinc’ RF function generates a square pulse in the temporal domain. Experimental cross-correlation of the generated square pulse is represented: (a) Cross-correlation of the output from only ‘sinc’ RF modulation; (b) same ‘sinc’ RF with a linear sweep and pre-compensated ‘sinc’ RF modulation with a linear sweep.
Fig. 11.
Fig. 11. (a) The plot of the imaginary and real part of the pre-compensated RF signal that is generated by adding a third-order phase with the linear frequency sweep. Experimental cross-correlation signal and expected square pulse (using αRF) are presented in (b) that are quite good agreement with each other.

Equations (9)

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

ξout(ωl)=ξin(ωl)M(ωl)
λ=d(sinϕsinθ)andx=f×tan(θλ±θλ0)
f(t+α)FourierTransforme2πiωαF(ω)
f1(tRF)=sin(2πνRF1tRF)+0.35sin(2πνRF2tRF)+0.25sin(2πνRF4tRF)+0.1sin(2πνRF6tRF)
f2(tRF)=sin(2πνRF2tRF)+0.2sin(2πνRF4tRF)+0.4sin(2πνRF6tRF)whereνRF1=1MHz, νRF2=2MHz,νRF4=4MHz,νRF6=6MHz
f1(pRF)=sin(2π.1.pRF)+0.35sin(2π.2.pRF)+0.25sin(2π.4.pRF)+0.1sin(2π.6.pRF)
f1(pω)=sin(2π.1.aNpω)+0.35sin(2π.2.aNpω)+0.25sin(2π.4.aNpω)+0.1sin(2π.6.aNpω)
fsinc(tRF)=sin(πνRF10tRF)πνRF10tRF.exp[i(2πνRF10tRF20πνRF10tRF3)]whereνRF10=10MHz
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.