Abstract
Two-dimensional electronic spectroscopy based on passive phase stabilization methods is now well known and widely employed worldwide. In the most recent fully non-collinear implementations, a great phase stability is often achieved at the expense of the independent control over the pulse timings, limiting the full potential of the technique. Here we propose several modifications in the experiment geometry, calibration procedures, and data acquisition and processing routines. The setup is easily tunable to record different phase-matching directions, such as rephasing, non-rephasing and double quantum signals, still maintaining high levels of phase control and phase stability. The performances of the proposed setup are exemplified by measures on the standard dye zinc phtalocyanine and porphyrin J-aggregates.
Published by The Optical Society under the terms of the Creative Commons Attribution 4.0 License. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI.
1. INTRODUCTION
Two-dimensional electronic spectroscopy (2DES) is a third-order non-linear technique which is successfully spreading in the latest years thanks to its capability of unraveling fine details of the ultrafast dynamics of the systems under investigation [1–3]. 2DES is based on the interaction of three ultrafast laser pulses with an optical sample. The experiment consists of probing the resulting third-order optical polarization as a function of the time delays between the three different interactions. A typical pulse sequence with the definition of the time variables is reported in Fig. 1(a). The capability of 2DES to scan independently the three time delays between pairs of pulses allows disentangling overlapped signals and unraveling correlations between different transitions. Suitable pulse sequence definition allows indeed the access to rephasing (R), non-rephasing (NR), and double quantum (2Q) signals, each one providing different and complementary information [4,5]. This flexibility enables investigation of molecular systems with previously unavailable detail; for example, molecular couplings and energy relaxation pathways can be directly followed and the temporal evolution of both coherent and non-coherent processes can be characterized.
The first 2DES spectrometer [6] was based on a fully non-collinear BOXCARS geometry. Since then, other geometries have been proposed [7–11] and recently reviewed by Fuller and Ogilvie [12]. The fully non-collinear geometry remains, however, the most adopted configuration and most efforts are still being paid in the development of new technological solutions to overcome the well-known issues and to improve the quality of the data.
One of the major sources of signal degradation in fully non-collinear setups is the lack of phase stability, i.e., the capability to keep constant the phase difference between pairs of optical pulses for a given time delay [13].
Many different approaches have been used to increase the phase stability in fully non-collinear 2DES experiments, including active [14] and passive phase stabilization methods, through diffractive optics (DO) [15,16], pulse shaping [17] and, more in general, particular geometries that allow the phase-locking of specific pulse pairs [18,19]. These last methods, in particular, achieve a great phase stability working in a full optical rotating frame, but they require coupling the motion of pairs of pulses, partially losing the independent control of timings. Therefore, a better phase stability is obtained at the cost of a less versatile setup, unable to exploit the full potential of 2DES.
Here we report an experimental scheme for a fully non-collinear 2DES whose optical layout has been specifically designed to provide full independent control over all the time intervals allowing an easy retrieval of all the different phase-matching signals (R, NR, and 2Q), still maintaining excellent reliability and particularly high short- and long-term phase stability.
Experimental data recorded for a tetrahydrofuran (THF) solution of zinc phthalocyanine (ZnPc) are used to elucidate the crucial features of the proposed setup and of the corresponding data processing routines. The full potential of the described approach is exemplified through the characterization of the full 2D response of solutions of J-aggregates of the diacid form of the water-soluble tetra-( sulphonato)-phenyl-porphyrin (TPPS).
2. OPTICAL LAYOUT AND PULSE SEQUENCE
The experimental setup is illustrated in Fig. 1. The scheme is inspired to the BOXCARS geometry setup proposed by Nemeth et al. [16]. A Coherent Libra laser system generates a continuous pulse train centered at 800 nm with a repetition rate of 3 KHz. The central wavelength of laser pulses is converted in the visible range using a commercial NOPA (Light Conversion TOPAS White).
Before entering the optical setup, the pulses from the NOPA pass through a stage of pulse shaping [with an acousto-optic programmable dispersive filter (AOPDF), Fastlite Dazzler] and prisms compression to maximize the time resolution of the experiment and to guarantee full control over the phase of the different frequency components in the broadband pulse. A detailed description of the pulse shaping stage is reported in Section 3.
Four identical beams in BOXCARS geometry are obtained focusing the incoming beam into a 2D-DO element with spherical mirror SM3. The DO is a 2 mm thick phase mask that is etched in a 2D cross-hatched pattern where the width and the depth of the features are optimized for first-order diffraction of the incident wavelength. Since the diffraction properties of the DO are wavelength dependent, different suitably optimized DOs could be necessary to achieve maximum performances in different spectral ranges. The four first-order diffraction beams, each about 10% of the input power, are arranged in the BOXCARS configuration, as depicted in Fig. 1(b). The outgoing divergent beams are collimated and parallelized by a donut-shaped spherical mirror (DSM1), with the same focal length of the focusing mirror SM3 (500 mm).
The main advantage of using a 2D-DO is the more easy retrieval of the BOXCARS phase-matching conditions, naturally emerging from the 2D diffraction pattern, whereas the use of a 1D-DO requires an additional preliminary step of beam splitting through conventional optics and the subsequent alignment of pairs of beams to form a square. A drawback of this configuration is that the time delays between pulses can be introduced conveniently only through antiparallel wedge pairs (WP). Typically, a 1° fused silica wedge pair is employed, which is capable of providing a time delay of the order of hundreds of femtoseconds. This is enough for but it can be limiting for , where one might be interested in observing phenomena with longer characteristic time scales (of the order of picoseconds). In the proposed setup this complication is overcome through the use of 4° 25 mm long wedges, providing time delays up to 2 ps. Instead of fused silica, the choice of is more suitable since it induces less phase distortion in broadband pulses. Such wedges require anyway a more careful calibration procedure, as described in Section 4. The insertion of wedges in the beam path could produce small deviations in the beam propagation. This is detrimental for the fulfillment of the phase-matching conditions, because any deviation, even of a small angle, on the direction of one of the beams can produce artifacts in the collected signal [20]. To correct these possible deviations, the geometry and spatial propagation of the beams is controlled with a charge-coupled device (CCD) camera, allowing the optimization of the wedges orientation.
One wedge from each pair is mounted on a computer-controlled linear delay stage (Aerotech Ant 95), whose position has to be accurately controlled to guarantee repeatability and precision of the time delays between pulses (in our setup the repeatability is of about 0.07 fs and the phase stability is higher than ). The fourth pulse is used as a local oscillator (LO) and its power is conveniently attenuated with a graduated neutral filter in order to maximize the heterodyne signal. Its temporal position is kept fixed during the 2DES experiment and thus it can be conveniently set as reference time for the definition of the time delays.
This configuration enables us to scan independently all the different time delays of the three exciting beams, easily obtaining the R, NR, and 2Q signals simply designing the correct pulse sequence. We found that the most suitable configuration consists of fixing the time position of the pulse generating the third interaction with the system, i.e., for the R and NR signals and for the 2Q signal. The other two pulses are scanned accordingly to construct the full 2DES spectra. This choice has the major benefit of guaranteeing a phase-locking relation between () and the LO, leading to a better phase stability. Moreover, in this way, the signal is always generated at about the same time delay from the LO and its spectral interference does not change much during the experiments, thus avoiding all the artifacts arising when the signal is moving in time. For example, when the recorded spectral interference is pushed near the limit resolution of the detection system, false decay of the signals can be possibly generated.
It must be recognized, however, that this choice can generate different issues. It indeed implies scanning the first two pulses in different ranges of their associated wedge positions, possibly resulting in different chirp in pulse 1 and 2 for each population time . This complication can be overcome through a careful use of the pulse shaper that can be suitably tuned so to recompress the pulses at every delay (Section 3). Furthermore, if the phase difference between the third interaction and the LO is kept fixed, then the phase stability depends critically only on the phase difference between the first two interactions. It is thus necessary to achieve a strict control over the delay (and also for 2Q signals) and therefore particularly reliable wedge calibration procedures have been set up (Section 4).
Once the wanted time ordering is set, the four beams are focused onto the sample by DSM2 with a focal length of 200 mm. The third-order signal propagates in the LO direction and their ensuing interference is delivered to the spectrograph (Andor Shamrock 303i) equipped with a CCD camera (Andor Newton 920) for detection. The CCD is calibrated exploiting the capability of the AOPDF to produce holes in the pulse spectrum. The AOPDF is calibrated using certified absorptive filters (holmium oxide solution and didymium glass from Hellma).
During the 2DES experiment, the CCD camera collects other spurious contributions besides the signals. A suitable acquisition method and a post-processing procedure are thus necessary to remove the unwanted contributions and gain a high signal-to-noise ratio. The acquisition methodology is based on the double modulation lock-in method proposed by Augulis and Zigmantas [21] and will be discussed in detail in Section 5.
3. COMPRESSION AND SHAPING STAGE
The pulse properties are determined by two factors: the shape of the laser spectrum and the phases of its spectral components. The shape of the laser spectrum affects the shape of the final 2DES maps [22,23] while a distortion of the phase (like in chirped pulses) leads to longer pulse duration, reducing the resolution of the experiment [12], and to altered 2DES maps [24,25]. Therefore, it is necessary to work with transform-limited (TL) pulses, defined as the shortest possible pulses for a given optical laser spectrum, to obtain high-quality 2DES data. This condition is achieved when the phase is a linear function of the frequency [26]. Possible non-linear deviations of the phase, mainly caused by propagation in transmissive optics, have to be compensated to retrieve the TL (and thus temporally compressed) pulse at the sample position. In the implemented setup, this is performed by a pulse-shaping stage consisting of an AOPDF and a prism compressor [27,28], used to relieve the self-induced chirp of the Dazzler crystal.
The AOPDF employs a Taylor expansion of the phase of the electric field around the central frequency of the pulse spectrum, . The terms of the expansion at different orders describe different properties of the temporal profile of the pulse. The linear term is related to the arrival time of the pulse, while the second order describes the linear chirp. Higher even (odd) terms correspond to phase distortions with point (axial) symmetry with respect to [29].
The measurement of the pulse profile is performed with the frequency-resolved optical gating (FROG) technique at the sample position [30] that reveals the frequency-resolved auto-correlation of the pulse. In our setup, we can take advantage of the BOXCARS geometry to perform FROG by replacing the sample with a non-zero third-order material, typically a solvent. The use of the same cuvette in the same position for sample and solvent guarantees identical operating conditions in FROG and 2DES measurements, keeping the same pulse properties in the two experiments. The procedure to recover a TL pulse consists of manually modifying the coefficients of the expansion via the AOPDF and checking the FROG result until the best conditions are found.
In the experiment the delay times are set modulating the thickness of the wedge medium traversed by the beams. This implies that, at different time delays, the pulses can be affected by different chirp and thus the optimal compression conditions are position dependent. Therefore we must evaluate if and when the chirp is negligible or how we can correct its impact on the final signal. This issue is particularly relevant in the scan during R and NR experiments, where the full thickness of the wedges is exploited. The chirp has two main effects: a decrease of the maximum signal intensity and the presence of distortions because the medium dispersion leads to frequency-dependent time delays.
Simulation of FROG time traces [29] can be used to predict the loss of heterodyne third-order signal as a consequence of the chirping of the pulses. This procedure is applied to monitor in particular how the maximum intensity of the simulated FROG decays as a function of the chirp induced during the scan. The pulse duration, specified as full width at half-maximum (FWHM), must be accounted for in the discussion; indeed the shorter the pulse duration, the larger the bandwidth of the laser spectrum and the bigger the effect of the dispersion. The contour plots in Fig. 2 report the ratio between the maximum of the FROG traces simulated with chirped pulses and FROG traces simulated with optimal TL pulses, as a function of the pulse duration and the imposed delay. Deviations from the ideal fully compressed condition are studied as function of wedge material [panel (a)], different pulse shaping operating parameters [panel (b)], and different pulse spectrum central frequency [panel (c)]. It is clear that the greatest aberrations are expected for short pulses at long time delays. Therefore, the capability to manage large bandwidth, and thus short pulses, sets the limit of the maximum delay that can be investigated with a negligible loss of signal. In this context the use of wedges allows to reliably scan a time interval about 30% longer than fused silica wedges [Fig. 2(a)]. Moreover the capability of managing a specific bandwidth depends sensibly upon the central wavelength of the pulse spectrum [Fig. 2(c)]. The same simulations can also be used to find the best operating conditions of the pulse shaper at each investigated value. In our setup the AOPDF is indeed programmable to minimize the chirp effect at each . Since the AOPDF is positioned before the DO, it affects in the same way all the pulses. Thus it cannot operate independently on the different pulses of the experimental sequence and the optimal TL conditions for pump pulses ( and ) do not necessarily correspond to optimal conditions for the probe () pulse. The simulations in Fig. 2(b) confirm that the best choice is to set the pulse shaper apparatus so to maximize the pulse compression of the pump pulses at each . The Sellmeier expression of refractive index [31] is used to obtain the expansion coefficient to program the AOPDF. The choice of programming the AOPDF to keep compressed the pump pulses with respect to the probe pulse allows for another gain of about 30% in the exploitable time interval, see Fig. 2(b). The effect of the medium-induced chirp and how it can be effectively corrected through a smart programming of the AOPDF is illustrated in Fig. 3, reporting an example of FROG measurement recorded at different wedge positions.
The procedure described above is able to correct and minimize the loss of signal due to the chirp but it leaves unmodified the frequency-dependent timing distortions due to the dispersion induced by the wedges. It follows that the higher-frequency components are subjected to a longer delay than the lower-frequency components, resulting in a frequency shift after Fourier transforming from the time to the frequency domain. This leads to possible distortions of the signal along the excitation axis in 2DES maps and frequency shifts of the coherence beating detected along the axis. In order to estimate and possibly correct these aberrations, the time delay imposed through wedge motion can be expressed as a function of the frequency:
where and are the refractive indices of and air, respectively and is the time delay determined for the central frequency of the laser spectrum. The pulse chirp can be estimated using the group velocity dispersion of the medium, the thickness of the interposed material and the bandwidth taken as FWHM. Equation (1) can be used to evaluate the timing distortions due to the dispersion and to correct them, if needed, during the signal processing operations. For example, for a TL laser pulse centered at 600 nm and with 15 fs pulse duration, a chirp of 10 fs is estimated when a 2000 fs time delay is imposed by wedges, corresponding to a timing error of 0.5% at the wedges extreme positions.4. TIME DELAYS CALIBRATION
A crucial step in setting a 2DES experiment is the calibration of translation stages that regulate the wedge positions and the amount of that the pulse passes through. A sub-femtosecond precision of the time delay within the entire scanned range is necessary to achieve a full control over the phases of the pulses and therefore obtain a meaningful 2DES map. In a non-collinear 2DES measurement the phase of the signal depends indeed on the phase of all pulses. When the time delay is controlled by a translation stage, its accuracy should be adequate to conserve the phase difference between pulses. The methodology for time delay calibration here proposed is based on the determination of the sub-cycle time delays from the local oscillator [17].
The starting point is the definition of a linear relation that associates the position [in millimeters (mm)] of the th linear stage, and thus the position of the optical wedge, to a time delay [in femtoseconds (fs)] of the th pulse: , where is the linear calibration coefficient (fs/mm) that has to be determined and (mm) are the zero time positions at which all the pulses have the same time delay from the LO. The LO is taken naturally as a reference because it is the only pulse whose timing () is fixed since it is not controlled by a WP.
The procedure of time delay calibration consists of three main tasks: (i) determining time zero of the experiment and therefore the positions of the three linear stages that satisfy the equation ; (ii) finding the linear coefficients to convert a displacement of the stage in a time delay from the zero of the experiment; and (iii) finding the precise value of .
The typically employed calibration procedure consists of recording the spectral interference (SI) produced by two pulses [32,33]. To this purpose, a 25 μm pinhole is placed where the four beams are focused at the sample position and the SI is recorded by the spectrometer. The Fourier transform of the SI is the simplest method to correlate the position of the stages with the time delay between the LO and the selected pulse: when a spectral interference between two pulses is Fourier transformed, a peak in the time domain is obtained, which is centered at the delay time between the interfering pulses. As this delay increases, the number of fringes in the resulting SI increases shifting the peak in the time domain at longer times. This method works only until the CCD camera and the spectrograph can spectrally separate the fringes of the interference; therefore, it could fail at long values of delay times such as 2 ps.
Moreover, the adopted pulse sequence with the third interaction fixed and phased-locked with the LO, makes particularly crucial a full control over the relative time and phase of the first two pulses. Therefore, a new procedure has been proposed to fulfill these requirements. With respect to previously proposed approaches [33,34] where each wedge is calibrated independently, our method allows a more robust determination of time-zero positions thanks to the delay correlations plots , that correlate the positions of couples of stages (Section 4.A). Moreover, a further step of refinement of the linear calibration coefficients has been added, ensuring a reliable phase-locking condition between the first two pulses, effectively scanned in our experimental configuration (Section 4.B).
The first step consists of estimating a set of positions of the three stages that roughly induce the same delay from the LO. This evaluation can be performed by adjusting the stages in order to produce similar interference patterns. We subsequently scan each stage in a short range, about 20% of the length of the wedges, around and obtain three matrices of spectral interferograms, recorded for different positions and resolved in frequency by the spectrograph [see Fig. 4(a)]. Every column of the matrix represents an interference spectrum, which evolves as the delay between the interfering beams changes due to the motion of the linear stage. At this point we have three matrices, , , and , that contain all the information about the evolution of the pulse delays with respect to the fixed LO as the stages move.
A. Evaluation of Zero Positions
The procedure is based on differences of interference spectra. One of the exciting beams, for example , is selected as the reference pulse. One at a time, the remaining two beams are compared with the reference. Here the procedure is illustrated for beam , but the same will be applied for beam . The main objective of this processing is to produce, from the interference patterns and , a delay correlation plot called , which correlates the positions of the two stages, showing a clear signature when the two pulses are synchronous. The elements of matrix are defined as
where the index runs over the frequency axis and the indices and run over the positions of the stages and , respectively. For a specific value, the data matrix with elements is retrieved. This matrix presents a characteristic checkerboard pattern, as shown in Fig. 4(b). The absolute value of this matrix is spectrally integrated, obtaining an oscillating function of [solid black line in Fig. 4(c)]. This oscillatory trace has its biggest amplitude in the correspondence of the position of the stage that assures an equal delay for the compared pulses. The envelope of this oscillating curve is computed subtracting the constant part and taking the absolute value of the analytic signal [35], obtaining the dashed black line in Fig. 4(c). The matrix is then built collecting all the envelopes derived from all the possible indices [Fig. 4(d)]. A linear function can be used to fit the trend observed in matrix . The function is defined as where and are the slope and the intercept of the linear function. Each point on the line can be used as actual zero position of the calibration. In other words, with this linear function we can choose an arbitrary position as zero for the reference stage and automatically get the zero position for the other stages. In order to practically and reliably obtain the parameters and of the linear function, we can fit matrix with a Gaussian surface represented in the matrix . The center of the Gaussian function is modeled with the linear function in Eq. (3) obtaining a Gaussian tunnel shape. The matrix elements of are defined as where is a scaling factor, is the width of the Gaussian function, and the indices and are associated with the -th and -th positions and , respectively. Then the following unconstrained least-square minimization problem is solved: given a good set of starting parameters the Gaussian surface is rapidly fitted. The same procedure is repeated for the third pulse, so that for any arbitrary choice of , the zero stage positions and , in which the three exciting pulses are contemporary, are determined. This condition is used as the time-zero of the experiment.B. Evaluation of the Calibration Coefficients
The second part of the calibration procedure consists of determining the linear coefficients , , and that allow accurately moving the pulses from the zero positions. The optimal procedure should provide conversion coefficients able to guarantee, at the maximum excursion of the linear stages, a time error much smaller than the optical cycle of the electromagnetic radiation. This accuracy assures that the final signal, recorded during experiments, is correctly phase-locked. The idea at the basis of this procedure is to mathematically compensate for the evolution of the spectral interference coupled to the movement of the wedge. In other words, we want to determine which is the exact delay to be mathematically subtracted in order to achieve a perfectly constant interference pattern. We will describe the procedure for the pulse ; the same applies for the other pulses.
First, matrix is Hilbert transformed along the frequency axis in order to operate on the complex analytic signal. This allows subtracting the delay from the spectral interference multiplying by an oscillating exponential function as
where for the correct value the matrix is a collection of constant interference spectra. This operation cancels out the time delay induced by moving the WP from the zero position . In order to retrieve the exact parameter , we can solve a maximization problem for an objective scalar function which exploits the condition of the constant interference pattern that we want to achieve. The problem is defined as where the first summation over the index , which runs over the stage positions, is going to output non-null values only if there are no oscillations along that direction. The summation over serves to obtain the final scalar objective function of the maximization. All the parameters are then recovered repeating the procedure for the three stages.The interference spectra in matrices used for the previous calculation are sampled only within a limited portion of the stage movement range. This is because at large delay values the spectrometer is not able to efficiently resolve all the fringes of the SI. This is not detrimental for the accurate evaluation of the zero positions of the stages but it may generate inaccuracy at the extreme positions of the stages. This is particularly critical for pulses and , whose position in time is effectively scanned during the experiment, exploiting the full length of the optical wedges. The third pulse and LO are instead fixed during R and NR measurements, and are not delayed relevantly in 2Q experiments. Therefore, we verify that the relative phase between and is conserved also at the extreme positions and possible small errors, even of fractions of an optical cycle, are corrected.
In practice, using the calibration coefficients previously determined, and are moved simultaneously so to keep constant. The interference pattern generated by interaction of and as a function of the position (or ) and is recorded and the matrix is retrieved. If the spectral interference is constant for the full scanned range, then the calibration is accurate and ready for the experiments. When a small linear drift of the phase of the interference is observed, it means that the values of the linear coefficients are not fully rigorous in the whole range and they must be corrected to avoid phase drift of the third-order signal.
To this purpose, it is possible to define a corrected interference pattern as
where is the small delay correction to be applied, and is an adimensional scaling parameter. Analogously to problem (7) we can maximize a scalar function solving the following problem: The phase drift is caused by an error in the relative time positions and not in the absolute time positions of the two pulses. Assuming that the correction is small, it is convenient to include the whole corrective factor only in the calibration coefficient of one of the two pulses, e.g., we can take as exact and impose the correction factor only on . In these conditions, it can be demonstrated that the corrected coefficient can be expressed by This correction is usually smaller than the 0.05%, but it is fundamental for removing any phase trend of the third-order signal induced by calibration uncertainty. Figure 5 compares 2DES maps at three different population times , obtained with [Figs. 5(d)–5(f)] and without [Figs. 5(a)–5(c)] the linear correction of the calibration coefficients. A clear spurious phase trend is introduced for long values if the correction in Eq. (10) is not accounted for. We verified that in our setup the phase drift has a predominant linear trend as a function of the positions. After the correction described in Eq. (10) we retrieved the trend shown in Fig. 5(g) that demonstrates that the linear trend is indeed well suited. However, manufacturers produce high but not absolute precision wedges. Indeed, Fig. 5(h) shows that, after the correction of the linear coefficient in Eq. (10), small phase deviations are still present and have to be corrected during data processing.C. Evaluation of the Delay from the LO
The last step consists of retrieving the precise value of . This task cannot be accomplished by simply looking at the SI generated by the 25 μm pinhole used for the calibration, because this configuration does not reproduce conveniently the exact working conditions of 2DES experiments. We therefore suggest to exploit the heterodyne non-resonant FROG of a solvent in a cell at the sample position. The procedure consists of collecting the interference patterns produced by the LO and the third-order non-resonant signal of the solvent generated when all the pulses are superimposed in time. During this experiment, all the stages positions are scanned simultaneously around their values, so to keep . Therefore, we obtain a matrix of interference patterns which evolves as a function of , similar to matrices . This matrix is Fourier transformed along the frequency axis obtaining a new matrix, which has a tunnel-like shape similar to matrix . A minimization procedure analogous to the one reported in Eq. (5) can be applied and a more accurate value of at can be extrapolated.
The determination of the exact value of is a crucial parameter in the processing of the 2DES spectrum (see the next section). The time distance of the LO from the zero of the experiment () controls indeed how dense the fringes of the interference pattern in the collected spectra are, usually is set at about 400 fs. As an example, Fig. 6 shows the effects of an error in the estimation of in the final 2DES rephasing map.
5. DATA ACQUISITION AND PROCESSING
A double lock-in modulation acquisition [21] is implemented in order to suppress the spurious contributions in the recorded spectra. In this acquisition procedure, and exciting pulses are modulated by two phase-locked optical choppers at two different frequencies, and , in order to isolate the signal from almost all the other contributions after a Fourier transform operation. The laser pulse train, the CCD camera, and the optical choppers are reciprocally triggered and phase-locked. To fully exploit the potential of this method, each modulation frequency must be carefully chosen. First, the signal should be extracted at a frequency lower than the Nyquist limit, so that , with being the sampling frequency of the detector. Second, to avoid phase instability it is convenient to set the higher modulating frequency as a multiple of the lower one. However, since the choppers modulate the beams with square waves, the scattering contributions modulated at the lower frequency are distributed in all their odd harmonic frequencies. Therefore, taking the higher frequency as the odd multiple of the lower one is a convenient approach, since it guarantees that the modulation frequencies of the signal are always even multiples of the slower modulation. The signal can be eventually retrieved at both sum and difference frequencies, and , respectively. The signal is recovered as the real part of the Fourier transform of the series of spectra acquired by the detector. There is a particular combination of relative phases between the two modulations and the CCD sampling for which the real part of the signal is maximized at both sum and difference frequencies and the maximum signal-to-noise ratio can be achieved. In our setup, , , and ; thus the signal can be extracted at or or both. With this method, most of the scattering contributions are removed, since they reach the CCD camera with the frequency of modulation of a single chopper or they are not modulated at all.
However, the acquired signal still contains the unwanted homodyne third-order signal, , the spurious contribution arising from the interference between the signal and the scattering of the non-modulated beam, , and the interference between the scattering of the modulated beams, , so that the raw signal acquired by the CCD can be expressed as . An example of rephasing data collected at a fixed population time is reported in Fig. 7(a).
First, to isolate the pure signal from the contribution , it is necessary to enforce the time ordering between the signal and the local oscillator. The recorded interference can be expressed as , where and are the square root of the spectra of the corresponding fields. The local oscillator arrives at the CCD after the signal, thus we need to impose that is negative. This operation is performed (i) by inverse-Fourier transforming along , (ii) by retaining only the negative part of the resulting time axis, and (iii) by Fourier transforming back. After this set of operations, the signal becomes a complex valued signal : the real part is equal to the data recorded by the CCD and the imaginary part is a copy of the original data with a phase shift of . The delay can now be added and at the same time the data are divided by the square root of the LO spectrum giving
as reported in Fig. 7(b). Equation (11) highlights how an accurate estimation of is crucial to retrieve artifacts-free 2D maps, as already discussed in Section 4.C.At this point a rotating frame (RF) can be applied, as previously suggested by Schlau-Cohen et al. [36]. For R and NR signals, this is done detuning the optical frequency during by subtracting a reference frequency , which is usually taken as the central frequency of the laser bandwidth:
After this operation, the resulting R and NR signals, [Fig. 7(c)] and , respectively, retain only slow modulations along the coherence time . It follows that their sampling frequency in can be greatly reduced. The maximum frequency observed in the RF signal corresponds to half of the bandwidth of the laser spectrum, which is of the order of few thousands of . Note that the RF is introduced here a posteriori after the data acquisition and only at the level of data processing calculation. The procedure does not require specific pulse scan schemes, leading to maximum freedom and flexibility in the control of pulse timings. This post-processing approach allows acquiring and managing undersampled data, greatly reducing the acquisition time. It does not improve the quality of the data but speeds up the acquisition of a single map that can be anyway averaged in case of a low signal-to-noise ratio. Note that the undersampling of 2Q signals is not necessary because the final response is a single map and therefore the acquisition time is already very fast [37].In order to remove the remaining spurious contributions from the signal, a time-filtering in the dimension is applied using a super-Gaussian function. The filter is tuned to keep only the portion of the axis between zero and about 100–150 fs, after which the signal is completely decayed. This operation removes most of the remaining spurious contributions in the recorded signal, including the homodyne third-order signal, because the operation in Eq. (11) shifts all these contributions by and thus outside the filter window. The initial transient of the signal can be slightly affected by and when the delay between signal and and the delay between and are inside the filter window, respectively. This condition is verified only at early values of population time and does not arise any concern for the rest of the data.
The small dispersion of the pulse spectral components affecting and delays is corrected applying Eq. (1) and interpolating the new data. This step of data processing is necessarily performed after the RF because the interpolation of the axis is not feasible when dealing with undersampled data.
The next step is the Fourier transform along (or in the case of 2Q maps), to obtain the final frequency-frequency maps. The resulting excitation frequency axis must then be shifted by to account for the rotating frame previously applied. Finally, the arbitrary phase of the signal with respect to the LO must be removed to isolate the real from the imaginary part. As suggested in [2], this last step is performed exploiting the projection-slice theorem. The same phase shift is employed for the phasing of R, NR, and 2Q signals, as long as these measurements are performed with the same calibration conditions.
At this point, the processed signal and its time evolution along is ready to be examined with the various data analysis procedures proposed to extract information about populations and coherent dynamics [38–48]. The versatility of this setup is exemplified in Fig. 8 that reports the pump-probe spectrum, the R and NR maps recorded at , and the 2Q map for a solution of J-aggregates of TPPS. This particular system represents an ideal example of how the synergic use of high-quality R, NR, and 2Q response can provide a clear picture of very subtle details of the excitonic structure and dephasing mechanisms in complex systems, not accessible by single measurements [5].
6. CONCLUSIONS
We reported a detailed description of a versatile fully non-collinear setup for 2DES spectroscopy. Novel implementations in the optical layout and in calibration and data processing procedures have been proposed to efficiently obtain high-quality and artifacts-free 2DES spectra along different phase-matching directions. With respect to previously proposed fully non-collinear wedge-based setups, we introduced the use of 4° wedge pairs, allowing longer scanning times (up to 2 ps). The particular optical layout of the setup allows an independent control the three pulse delays, giving access to rephasing, non-rephasing, and double quantum signals simply imposing the correct pulse sequence. Possible aberrations introduced especially at long delay times by the dispersion of the wedge medium have been estimated and discussed. Correction methodologies based on a fine control of the phases of the pulses with an AODPF pulse shaper and a robust wedge calibration procedure have been proposed. Reduced acquisition times were achieved thanks to the application of an a posteriori rotating frame approach that allows a drastic reduction of the number of data points along to acquire a single map.
Funding
H2020 European Research Council (ERC) (QUENTRHEL (278560)); Seventh Framework Programme (FP7) (MULTI (317707), PAPETS (323901)).
REFERENCES
1. E. Collini, “Spectroscopic signatures of quantum-coherent energy transfer,” Chem. Soc. Rev. 42, 4932–4947 (2013). [CrossRef]
2. D. M. Jonas, “Two-dimensional femtosecond spectroscopy,” Annu. Rev. Phys. Chem. 54, 425–463 (2003). [CrossRef]
3. A. M. Brańczyk, D. B. Turner, and G. D. Scholes, “Crossing disciplines—a view on two-dimensional optical spectroscopy,” Ann. Phys. 526, 31–49 (2014). [CrossRef]
4. A. Nemeth, F. Milota, T. Mančal, T. Pullerits, J. Sperling, J. Hauer, H. F. Kauffmann, and N. Christensson, “Double-quantum two-dimensional electronic spectroscopy of a three-level system: experiments and simulations,” J. Chem. Phys. 133, 094505 (2010). [CrossRef]
5. L. Bolzonello, F. Fassioli, and E. Collini, “Correlated fluctuations and intraband dynamics of J-aggregates revealed by combination of 2DES schemes,” J. Phys. Chem. Lett. 7, 4996–5001 (2016). [CrossRef]
6. J. D. Hybl, A. W. Albrecht, S. M. Gallagher Faeder, and D. M. Jonas, “Two-dimensional electronic spectroscopy,” Chem. Phys. Lett. 297, 307–313 (1998). [CrossRef]
7. 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, 748–761 (2009). [CrossRef]
8. J. A. Myers, K. L. M. Lewis, P. F. Tekavec, and J. P. Ogilvie, “Two-color two-dimensional Fourier transform electronic spectroscopy with a pulse-shaper,” Opt. Express 16, 17420–17428 (2008). [CrossRef]
9. E. Harel, A. F. Fidler, G. S. Engel, and B. Alexander Pines, “Real-time mapping of electronic structure with single-shot two-dimensional electronic spectroscopy,” Proc. Natl. Acad. Sci. USA. 107, 16444–16447 (2010). [CrossRef]
10. C. Manzoni, D. Brida, and G. Cerullo, “Phase-locked pulses for two-dimensional spectroscopy by a birefringent delay line,” Opt. Lett. 37, 3027–3029 (2012). [CrossRef]
11. F. D. Fuller, D. E. Wilcox, and J. P. Ogilvie, “Pulse shaping based two-dimensional electronic spectroscopy in a background free geometry,” Opt. Express 22, 1018–1027 (2014). [CrossRef]
12. F. D. Fuller and J. P. Ogilvie, “Experimental implementations of two-dimensional Fourier transform electronic spectroscopy,” Annu. Rev. Phys. Chem. 66, 667–690 (2015). [CrossRef]
13. V. I. Prokhorenko, A. Halpin, and R. J. D. Miller, “Coherently-controlled two-dimensional photon echo electronic spectroscopy,” Opt. Express 17, 9764–9779 (2009). [CrossRef]
14. A. D. Bristow, D. Karaiskaj, X. Dai, T. Zhang, C. Carlsson, K. R. Hagen, R. Jimenez, and S. T. Cundiff, “A versatile ultrastable platform for optical multidimensional Fourier-transform spectroscopy,” Rev. Sci. Instrum. 80, 073108 (2009). [CrossRef]
15. M. L. Cowan, J. P. Ogilvie, and R. J. D. Miller, “Two-dimensional spectroscopy using diffractive optics based phased-locked photon echoes,” Chem. Phys. Lett. 386, 184–189 (2004). [CrossRef]
16. A. Nemeth, J. Sperling, J. Hauer, H. F. Kauffmann, and F. Milota, “Compact phase-stable design for single-and double-quantum two-dimensional electronic spectroscopy,” Opt. Lett. 34, 3301–3303 (2009). [CrossRef]
17. D. B. Turner, K. E. Wilk, P. M. G. Curmi, and G. D. Scholes, “Comparison of electronic and vibrational coherence measured by two-dimensional electronic spectroscopy,” J. Phys. Chem. Lett. 2, 1904–1911 (2011). [CrossRef]
18. U. Selig, F. Langhojer, F. Dimler, T. Löhrig, C. Schwarz, B. Gieseking, and T. Brixner, “Inherently phase-stable coherent two-dimensional spectroscopy using only conventional optics,” Opt. Lett. 33, 2851–2853 (2008). [CrossRef]
19. I. A. Heisler, R. Moca, F. V. A. Camargo, and S. R. Meech, “Two-dimensional electronic spectroscopy based on conventional optics and fast dual chopper data acquisition,” Rev. Sci. Instrum. 85, 063103 (2014). [CrossRef]
20. B. M. Cho, M. K. Yetzbacher, K. Kitney, E. R. Smith, and D. M. Jonas, “Propagation and beam geometry effects on 2D Fourier transform spectra of multi-level systems,” J. Phys. Chem. A 113, 13287–13299 (2009). [CrossRef]
21. R. Augulis and D. Zigmantas, “Two-dimensional electronic spectroscopy with double modulation lock-in detection: enhancement of sensitivity and noise resistance,” Opt. Express 19, 13126–13133 (2011). [CrossRef]
22. F. V. de A. Camargo, L. Grimmelsmann, L. Anderson, S. R. Meech, and I. A. Heisler, “Resolving vibrational from electronic coherences in two-dimensional electronic spectroscopy: the role of the laser spectrum,” Phys. Rev. Lett. 118, 033001 (2017). [CrossRef]
23. V. Perlík, J. Hauer, and F. Šanda, “Finite pulse effects in single and double quantum spectroscopies,” J. Opt. Soc. Am. B 34, 430–439 (2017). [CrossRef]
24. P. F. Tekavec, J. A. Myers, K. L. M. Lewis, F. D. Fuller, and J. P. Ogilvie, “Effects of chirp on two-dimensional Fourier transform electronic spectra,” Opt. Express 18, 11015–11024 (2010). [CrossRef]
25. N. Christensson, Y. Avlasevich, A. Yartsev, K. Müllen, T. Pascher, and T. Pullerits, “Weakly chirped pulses in frequency resolved coherent spectroscopy,” J. Chem. Phys. 132, 174508 (2010). [CrossRef]
26. R. Paschotta, Encyclopedia of Laser Physics and Technology (Wiley-VCH, 2008).
27. W. Dietel, J. J. Fontaine, and J.-C. Diels, “Intracavity pulse compression with glass: a new method of generating pulses shorter than 60 fsec,” Opt. Lett. 8, 4–6 (1983). [CrossRef]
28. R. L. Fork, O. E. Martinez, and J. P. Gordon, “Negative dispersion using pairs of prisms,” Opt. Lett. 9, 150–152 (1984). [CrossRef]
29. K. W. DeLong, R. Trebino, and D. J. Kane, “Comparison of ultrashort-pulse frequency-resolved-optical-gating traces for three common beam geometries,” J. Opt. Soc. Am. B 11, 1595–1608 (1994). [CrossRef]
30. D. J. Kane and R. Trebino, “Characterization of arbitrary femtosecond pulses using frequency-resolved optical gating,” IEEE J. Quantum Electron. 29, 571–579 (1993). [CrossRef]
31. M. Daimon and A. Masumura, “High-accuracy measurements of the refractive index and its temperature coefficient of calcium fluoride in a wide wavelength range from 138 to 2326 nm,” Appl. Opt. 41, 5275–5281 (2002). [CrossRef]
32. T. Brixner, T. Mancal, I. V. Stiopkin, and G. R. Fleming, “Phase-stabilized two-dimensional electronic spectroscopy,” J. Chem. Phys. 121, 4221–4236 (2004). [CrossRef]
33. R. Augulis and D. Zigmantas, “Detector and dispersive delay calibration issues in broadband 2D electronic spectroscopy,” J. Opt. Soc. Am. B 30, 1770–1774 (2013). [CrossRef]
34. D. B. Turner, R. Dinshaw, K.-K. Lee, M. S. Belsley, K. E. Wilk, P. M. G. Curmi, and G. D. Scholes, “Quantitative investigations of quantum coherence for a light-harvesting protein at conditions simulating photosynthesis,” Phys. Chem. Chem. Phys. 14, 4857–4874 (2012). [CrossRef]
35. S. L. Marple, “Computing the discrete-time “analytic” signal via FFT,” IEEE Trans. Signal Process. 47, 2600–2603 (1999). [CrossRef]
36. G. S. Schlau-Cohen, A. Ishizaki, and G. R. Fleming, “Two-dimensional electronic spectroscopy and photosynthesis: fundamentals and applications to photosynthetic light-harvesting,” Chem. Phys. 386, 1–22 (2011). [CrossRef]
37. J. Kim, S. Mukamel, and G. D. Scholes, “Two-dimensional electronic double-quantum coherence spectroscopy,” Acc. Chem. Res. 42, 1375–1384 (2009). [CrossRef]
38. I. H. M. Van Stokkum, D. S. Larsen, and R. Van Grondelle, “Global and target analysis of time-resolved spectra,” Biochim. Biophys. Acta 1657, 82–104 (2004). [CrossRef]
39. E. E. Ostroumov, R. M. Mulvaney, J. M. Anna, R. J. Cogdell, and G. D. Scholes, “Energy transfer pathways in light-harvesting complexes of purple bacteria as revealed by global kinetic analysis of two-dimensional transient spectra,” J. Phys. Chem. B 117, 11349–11362 (2013). [CrossRef]
40. F. V. d. A. Camargo, H. L. Anderson, S. R. Meech, and I. A. Heisler, “Time-resolved twisting dynamics in a porphyrin dimer characterized by two-dimensional electronic spectroscopy,” J. Phys. Chem. B 119, 14660–14667 (2015). [CrossRef]
41. A. Gelzinis, V. Butkus, E. Songaila, R. Augulis, A. Gall, C. Büchel, B. Robert, D. Abramavicius, D. Zigmantas, and L. Valkunas, “Mapping energy transfer channels in fucoxanthin-chlorophyll protein complex,” Biochim. Biophys. Acta 1847, 241–247 (2015). [CrossRef]
42. J. A. Myers, K. L. M. Lewis, F. D. Fuller, P. F. Tekavec, C. F. Yocum, and J. P. Ogilvie, “Two-dimensional electronic spectroscopy of the D1-D2-cyt b559 photosystem II reaction center complex,” J. Phys. Chem. Lett. 1, 2774–2780 (2010). [CrossRef]
43. D. B. Turner, K. W. Stone, K. Gundogdu, and K. A. Nelson, “Three-dimensional electronic spectroscopy of excitons in GaAs quantum wells,” J. Chem. Phys. 131, 144510 (2009). [CrossRef]
44. H. Li, A. D. Bristow, M. E. Siemens, G. Moody, and S. T. Cundiff, “Unraveling quantum pathways using optical 3D Fourier-transform spectroscopy,” Nat. Commun. 4, 1390 (2013). [CrossRef]
45. J. R. Caram and G. S. Engel, “Extracting dynamics of excitonic coherences in congested spectra of photosynthetic light harvesting antenna complexes,” Faraday Discuss. 153, 93–104 (2011). [CrossRef]
46. J. Prior, E. Castro, A. W. Chin, J. Almeida, S. F. Huelga, and M. B. Plenio, “Wavelet analysis of molecular dynamics: efficient extraction of time-frequency information in ultrafast optical processes,” J. Chem. Phys. 139, 224103 (2013). [CrossRef]
47. A. Volpato and E. Collini, “Time-frequency methods for coherent spectroscopy,” Opt. Express 23, 20040–20050 (2015). [CrossRef]
48. A. Volpato, L. Bolzonello, E. Meneghin, and E. Collini, “Global analysis of coherence and population dynamics in 2D electronic spectroscopy,” Opt. Express 24, 24773–24785 (2016). [CrossRef]