Phase-synchronous undersampling in nonlinear spectroscopy Lukas Bruder,∗ Marcel Binz, and Frank Stienkemeier Institute of Physics, University of Freiburg, 79104 Freiburg, Germany (Dated: November 9, 2018) We introduce the concept of phase-synchronous undersampling in nonlinear spectroscopy. The respective theory is presented and validated experimentally in a phase-modulated quantum beat experiment by sampling high phase modulation frequencies with low laser repetition rates. The advantage of undersampling in terms of signal quality and reduced acquisition time is demonstrated and breakdown conditions are identified. The presented method is particularly beneficial for exper- imental setups with limited signal/detection rates. Ultrafast nonlinear spectroscopy in the time domain conditions and show explicitly how PSU improves the has emerged as a powerful tool to study photophysi- signal quality in nonlinear spectroscopy. In particular, cal and photochemical processes in real-time with high we demonstrate undersampling by more than two orders spectro-temporal resolution [1, 2]. To meet the demand of magnitude without loss of performance in signal re- on interferometric phase stability, necessary in this type covery, hence, providing the opportunity to use the PM of spectroscopy, numerous stabilization methods have technique with sampling rates down to the 10 Hz regime. been developed in the past years [3]. Among these, one Our work makes the PM approach thus accessible to a particularly efficient concept is the phase modulation much wider range of experimental setups. (PM) technique, as demonstrated in pump-probe and To demonstrate and characterize the PSU concept, we multidimensional spectroscopy [4{11]. This method com- performed quantum beat measurements utilizing the PM bines acousto-optical PM with lock-in detection which technique in pump-probe configuration (Fig. 1). An am- drastically reduces the demands on phase stability and plified femtosecond (fs) laser system of variable repeti- greatly improves the overall sensitivity. As an additional tion rate (!rep ≤ 5 kHz) was used, tuned to λ = 795 nm advantage, it can be combined with 'action' signals for to drive the D1 transition in a low density rubidium va- detection, e.g. fluorescence [6{8], photocurrent [9{11] and por. The vapor was contained in a spectroscopy cell of mass-resolved photoion detection [5], which improves se- which the fluorescence was detected with a photo mul- lectivity and facilitates studies at low particle densities. tiplier tube (PMT). The PM technique and the experi- An integral part of the PM technique is the lock-in mental setup is described in detail elsewhere [4, 6]. detection which relays on modulating a signal with a Briefly, a collinear pump-probe pulse sequence is gener- well-defined frequency, followed by demodulation with a ated in a Mach-Zehnder interferometer. Thereby, a signal lock-in amplifier (LIA). This scheme works most efficient modulation (Ω) is introduced with two acousto-optical with high modulation frequencies. However, due to the modulators (AOMs) driven phase-locked at the frequency Nyquist theorem, the modulation frequency is generally difference Ω. The demodulation of the signal S is done limited by the sampling rate of the experimental appa- with a digital LIA, Stanford Research Systems, model SR ratus. In spectroscopy, limiting factors typically are the 810. To this end, a reference waveform R is constructed signal count rate, the speed of detectors or the laser repe- from a continuous wave (cw) laser (λ = 780 nm) superim- tition rate. This constrains the combination of lock-in de- tection with CCD camera based detection schemes (e.g. dispersed fluorescence [12], velocity map imaging detec- diode laser tion [13]) and, more severely, limits the technique to high BS repetition rate laser systems (! > 100 kHz [4, 7]). Yet, Ω1 DL rep fs laser PMT many applications still require operation at low sampling AOM1 S PH rates (!sa < 5 kHz). In particular, in the XUV spectral BS SR sample Ω2 LIA arXiv:1712.09109v2 [physics.chem-ph] 4 Jan 2018 range, where the PM technique shows great promise in facilitating nonlinear spectroscopy [14], laser sources op- AOM2 PH R erate mostly at 0.01-1 kHz. PD A common strategy to overcome the Nyquist limit FIG. 1. Optical PM setup. The pump-probe signal S is mod- is undersampling, as proposed in some lock-in detec- ulated using two AOMs driven at the frequency difference tion schemes, e.g. in lock-in thermography [15], electrical Ω = Ω2 − Ω1. A cw reference signal R is generated with a impedance spectroscopy [16], in combination with ran- diode laser and detected with a photo diode (PD). Pinholes dom sampling [17] and in case of very low signal count (PHs) block the diode/fs laser at either of the beam splitter's rates [18, 19]. In this letter, we add to this work by (BS) exit port, respectively. Signal demodulation is done with introducing phase-synchronous undersampling (PSU) to a LIA, referenced to the cw waveform R. The pump-probe the PM technique. We identify optimum undersampling delay is controlled with a motorized delay line (DL). 2 a) LIA dominant Lock-In sampling S(τ,Ωt) noise sources bandpass filter Nyquist rate LP SR(τ,t) SR (τ) R(τ,Ωt) 푆(휏) 푆(휏, Ω푡) 1/f noise Tavg noise floor b) 휔 SR up shift Ω 휔Nyquist 휔sa a FIG. 2. Signal-to-noise advantage by up-shifting S from low SR frequencies to Ω where the noise floor is much smaller and S SR is not dominated by lab noise. SRus posed with the optical path with a slight vertical offset t (Fig. 1), as also done in Ref. [9]. Throughout the data SR 1/2Ω 1/Ωa 1/ωsa run, the pump-probe delay τ is incremented in discrete c) steps ∆τ and the demodulated signal is recorded as a function of τ. This yields the temporal evolution of elec- tronic coherences induced in the system and a Fourier transform provides the system's absorption spectrum. t d) As an important feature, the phase-synchronous detec- tion of S and R leads to cancellation of phase jitter δφS. Thus, yielding passive phase stabilization while the de- sired information is deduced from the relative phase shift t of S and R, denoted φ (τ). As a second signal-to-noise SR FIG. 3. a) Schematic description of the lock-in amplifica- advantage, the imprinted modulation Ω shifts S from the tion scheme. LP denotes low-pass filter. b) Representa- low frequency spectrum, which is often dominated by lab tion of the product signal SR, as constructed in the LIA. noise, to a spectral region with less noise, where filtering Green: fully sampled signal (SR), black: undersampled sig- is much more efficient (Fig. 2). For this purpose, suf- nal (SRus), dashed black: resulting aliased waveform (SRa), ficiently large modulation frequencies (Ω & 1 kHz) are and red: mean value of SRus scaled according to the sampling highly desirable. However, if the sampling rate is small, duty cycle (hSRi). Undersampling of SR yields the correct average value/output signal, if the averaging interval T is i.e. ! ≤ 2Ω, the modulation is effectively shifted to its avg sa sufficiently large. c) Optimum undersampling case. d) Break- aliased frequency, that is down case: the output value hSRi strongly depends on the relative phase between the signal modulation and the sam- Ωa = minjn!sa − Ωj ≤ 0:5!sa; n 2 Z ; (1) pling points. and the signal-to-noise advantage of using large Ω would vanish. In our PSU scheme, this case is circumvented due frequency component: to the cw reference waveform used for the demodulation. SR(t; τ) / cos [2Ωt + φS(τ) + φR(τ)] + cos [φSR(τ)] ; (4) For a more detailed discussion of the PSU concept, see also Fig. 3b. While the modulation vanishes in the it is sufficient to consider the excitation of a model two- difference-frequency component, the sum-frequency term level system and exemplary examine the demodulation of exhibits a 2Ω-modulation. The subsequent RC-type low the in-phase signal component. In this case, the relevant pass filter can be described by a moving average over a signals S and R are given by [4] given time interval, denoted Tavg [17]. This removes the sum frequency component and thus the residual demod- S(t; τ) = A(τ) cos [Ωt + φ (τ)] (2) S ulated signal is R(t; τ) = cos [Ωt + φR(τ)] ; (3) hSRi(τ) = 0:5A(τ) cos [φSR(τ)] : (5) where A(τ) and φS(τ) reflect the system properties and φR(τ) denotes the pump-probe dependent phase func- Note, that hSRi contains the desired information, that is tion of the reference. In the standard lock-in algorithm, A(τ) and φSR(τ), while correlated phase jitter is removed S and R are multiplied and subsequently low-pass fil- (δφS(t) − δφR(t) ≈ 0). tered (Fig. 3a). The multiplication of the two waveforms Now, if operating in the undersampling regime, S yields a signal SR that consists of a sum- and difference- would appear at the respective aliased frequency (ΩS = 3 a) Ωa) and since R exhibits a different frequency (ΩR = Ω), 250 one would intuitively expect that the LIA blocks S with 225 high extinction ratio. However, in the lock-in amplifi- 200 cation process, undersampling of S simply leads to an 175 undersampled product waveform of SR (Fig. 3b). The 150 resulting aliased waveform SR exhibits the correct av- 125 a SNR erage value of SR, scaled according to the sampling duty 100 cycle. Hence, as long as the mean value is calculated over 75 SR 50 a sufficient number of aliased periods, i.e. Tavg 1=Ωa , the LIA will return a correct output signal hSRi.