A quantum cascade laser-based Mach–Zehnder interferometer for chemical sensing employing molecular absorption and dispersion

We employ a novel spectroscopic setup based on an external cavity quantum cascade laser and a Mach–Zehnder interferometer to simultaneously record spectra of absorption and dispersion of liquid samples in the mid-infrared. We describe the theory underlying the interferometric measurement and discuss its implications for the experiment. The capability of simultaneously recording a refractive index and absorption spectrum is demonstrated for a sample of acetone in cyclohexane. The recording of absorption spectra is experimentally investigated in more detail to illustrate the method’s capabilities as compared to direct absorption spectroscopy. We find that absorption signals are recorded with strongly suppressed background, but with smaller absolute sensitivity. A possibility of optimizing the setup’s performance by unbalancing the interferometer is presented.


Introduction
Mid-infrared (Mid-IR) absorption spectroscopy of liquids finds a wide range of applications in chemical sensing due to its simplicity, speed and outstanding selectivity. More recently, quantum cascade lasers (QCL) were employed successfully to replace notoriously weak thermal emitters as Mid-IR sources. In particular, external cavity (EC) QCLs proved to be an ideal match for many liquid sensing applications owing to their large spectral tuning range of typically 200 cm −1 [1][2][3][4][5][6][7][8][9]. To date, the vast majority of measurements are performed using direct absorption spectroscopy, in which the laser beam's intensity is recorded after passing through the sample or being totally reflected at its surface. Since this technique relies on recording small changes of the laser intensity, its sensitivity is ultimately limited by fluctuations of the laser power and noise from the detection system.
As of today, documented advantages of EC-QCL based systems over Fourier transform infrared spectrometers with thermal light sources are mainly increased ruggedness made possible by an enlarged optical path [10]. For example, recently an EC-QCL was successfully employed for the determination of protein secondary structures and protein folding in water using an optical path of 38 µm instead of 8 µm as required in FTIR spectrometry [10]. The so far achieved sensitivities of fully room temperature operated EC-QCL based systems compare well with those documented with high-end FTIR spectrometers which, however, employ liquid nitrogen cooled MCT detectors [11,12].
For gas phase measurements, monolithic distributed feedback (DFB) QCLs have already been established as reliable light sources void of mechanical imperfections limiting the performance of EC-QCLs, albeit of a far reduced tuning range. Nevertheless, considering the narrow ro-vibrational transitions in small gas molecules, a range of new measurement approaches expanding on direct absorption spectroscopy have been demonstrated to feature important advantages in terms of ruggedness as well as sensitivity [13,14]. Some of these novel approaches exploit molecular dispersion, rather than absorption, to decouple the recorded quantity (phase or frequency) from the laser intensity [15,16]. For liquids, these developments are yet pending. Linewidths of tens of cm −1 , the accordingly required lasers with tuning ranges beyond those achievable with DFB-QCLs, highly absorbing matrices and generally less well defined samples preclude a direct technology transfer from gas-to liquid spectroscopic sensing.
We present a novel, alternative approach to direct absorption spectroscopy for liquid samples that simultaneously and independently senses a sample's absorption and dispersion (refractive index). This is made possible by employing an interferometric principle of operation in a Mach-Zehnder type optical setup. Although single wavelength (integrated) Mach-Zehnder interferometers (MZI) are widely employed for (bio-)chemical sensing, especially in the visible spectral range [17,18], a report on directly recording a liquid sample's absorption and dispersion spectrum in the Mid-IR using a MZI is, to the best of our knowledge, yet pending. The reader is, however, pointed to infrared spectroscopic ellipsometry [19] and dispersive Fourier transform spectroscopy [20,21,22] that were used in (far-) infrared spectroscopic measurements of the complex dielectric function, mostly in reflection measurements.
We derive simple expressions relating the signals recorded in our setup to refractive index and absorption. Based on theory and experimental results we discuss distinct differences between direct absorption spectroscopy and our setup and how they can be exploited to extend the scope of QCL-based Mid-IR spectroscopy.

Setup
The setup is depicted schematically in Fig. 1. It is based on a free space MZI illuminated by an EC-QCL. A sampleand a reference liquid cell are introduced into the arms of the MZI, making the setup analogous to dispersive Fourier transform spectroscopy. The MZI is adjusted to equal optical path lengths in the two arms. To minimize fringing, all transparent optics, except for the liquid cells, are placed in the beam path in Brewster's angle and the light is p-polarized, resulting in the rhombic shape of the interferometer.
The two interfering partial beams leaving the interferometer to the right in Fig. 1 (one from the reference arm and one from the sample arm of the interferometer) are recorded at a thermoelectrically cooled MCT detector, while a BaF 2 -window is tilted to adjust the phase φ rel between them to π, i.e. to best possible destructive interference. The window is twice as thick as the beam splitters which compensates for dispersion picked up by the reference (upper) beam in Fig. 1 within the beam splitters' substrate. With no or equal samples in the liquid cells, the partial beams' amplitudes are equal, resulting in full destructive interference and no signal from the MCT-detector. Equal amplitudes are ensured by both beams being once transmitted and reflected at beam splitters of equal beam splitting ratios.
The basic principle of operation is as follows: When the sample liquid cell is filled with an absorbing species, the amplitudes of the interfering fields differ. Hence, destructive interference will be incomplete and a detectable signal is recorded that is a measure of the sample's absorption. At the same time, dispersion in the sample will cause a change in the relative phase of the two fields which needs to be compensated by the BaF 2 window's tilt. By reading the tilt, the sample's refractive index is recorded. The tilt of the window is measured via the position of a reflected HeNe-laser beam on a two zone photodiode. Optionally, the constructively interfering output of the MZI can be recorded by a second MCT detector and used to reference the destructively interfering signal to the laser power.  Fig. 1 Experimental setup. Light from the EC-QCL enters the Mach-Zehnder-type interferometer through a beam splitter and passes the reference-or the sample liquid cell. The intensity I des of the destructively interfering fields leaving the MZI is measured by two mercury-cadmium-telluride detectors and is a measure of the sample's absorption. Destructive Interference is ensured by regulating I des to its minimum by tilting a BaF 2 window in the sample beam path. The position sensitive photo diode facilitates the measurement of tilt, yielding the sample's refractive index During the experiment, the relative phase is set to π by tilting the BaF 2 window. In a symmetric interferometer with equal samples in both arms, rel does not change while scanning the wavelength and the window's tilt is constant. When the sample cell is loaded with a sample other than the reference, dispersion in the sample causes a spectrally varying additional phase φ that needs to be actively compensated by tilting the window to set rel = . This phase φ, recorded via the window's tilt, is a measure of the sample's refractive index n s as compared to the reference n ref and is given as Herein, the film thickness d was assumed equal for both cells. As can be seen from (2), the recorded phase is directly proportional to the sample's refractive index and hence, for diluted samples, its concentration c.

Basic equations
In the configuration described above, the intensities I s and I ref are equal, except for absorption in the sample (decadic absorption is used throughout).
Herein, I 0 is the laser intensity, T ref the transmission of the reference and R = 1 − T the reflectivity of the beam splitters.
Combining (1) and (3) and setting rel = yields For calibration purposes using Beer's law it is useful to write out A explicitly from (4).

Discussion of interferometric signals
When recording I des to measure concentrations or absorptions, some pronounced differences to direct absorption spectroscopy need to be taken into account. Most of them stem from the interferometric setup being sensitive to the amplitude of the transmitted E-field ('amplitude spectroscopy'), while direct absorption measures its power ('power spectroscopy') [22]. This behaviour is indicated by the square root of the fraction I des /I con in the logarithm in (5) rather than the fraction of two intensities (as is the case in direct absorption spectroscopy). The differences between amplitude-and power spectroscopy can be seen in the graphs of I s (direct absorption signal, compare (3)) and I des (Fig. 2) as a function of A as well as from Taylor series expansion of I des in A around A 0 = 0.
One important feature of the interferometric signal is that it is background free, i.e. I des (0) = 0 (no constant (3) (6)). However, as indicated by the vanishing linear term of (6), its sensitivity dI des ∕dA is zero too. The sensitivity increases with A until a maximum is reached at A = 2 ln (2)∕ ln (10) ≈ 0.6 , with a shallow role off towards larger A. This behaviour is in strong contrast to the direct absorption signal I s , which has its largest magnitude (background) and largest sensitivity at A = 0 , decreasing with A.
Due to the vanishing sensitivity at small A the MZI should in general not be operated with equally intense reference and sample beams. By increasing I ref relative to I s , the sensitivity can be strongly enhanced. This can be implemented by using beam splitters of different splitting ratios or-when dealing with samples featuring broadband background absorptions-by simply using a sample liquid cell of larger film thickness as compared to the reference. The latter approach potentially increases the method's sensitivity due to increased film thickness, but imbalances the optical path lengths in the interferometer, which requires additional care when interpreting the recorded spectra of I des and φ. A generalized theory for unbalanced interferometers is derived in the Appendix.
For large absorptions A > 1 , the sensitivity of the interferometric signal is much larger than the quickly decreasing sensitivity dI s ∕dA of the direct absorption signal. This makes the technique an interesting candidate when high concentrations need to be measured at good precision or for strongly absorbing solvents such as water. Also, the good sensitivity for large A yields a significantly larger dynamic range of I des as compared to direct absorption spectroscopy.
Just like in conventional FTIR spectroscopy, I des needs to be referenced to the laser power and the absorption of the Signals and Sensitivities / a.u.
Absorption / dec.AU In comparison to direct absorption spectroscopy, however, referencing is much less critical because the interferometer inherently references the sample beam to the reference beam. Drifting laser power solely scales I des (and retrieved A and c), but does not add constant errors to it, as is the case for power spectroscopy. A detailed comparative discussion of LODs in the presence of detector noise, shot noise and laser intensity noise is given in Sect. 6. The simultaneously recorded refractive index is independent from laser power and provides additional robustness towards laser power fluctuations. For diluted samples, refractive index spectra scale linearly with concentration and show highly characteristic features in the Mid-IR, making them a useful tool for chemical analysis.

Spectra acquisition and instrumentation
The optical setup was described above. For all measurements presented herein, two close to identical beam splitters with a 50/50 splitting ratio were employed. The liquid cells' film thickness was 1.9 mm. Uncoated BaF 2 windows were employed as windows for the liquid cells. With the liquid cells filled with the same sample, the optical path length in the two arms was found to be equal on the order of the wavelength used. This was confirmed by scanning the wavelength over approximately 200 cm −1 without changing the BaF 2 -window's tilt and recording the destructive signal. The absence of spectral oscillations in the recorded interferometric signal indicates equal optical path lengths. The modulation efficiency M of the interferometer, defined as the contrast between constructively and destructively interfering signals, M = I con (A = 0)∕I des (A = 0) , exceeded 23 dB throughout the experiments (fringes from the liquid cells not taken into account). It was limited by slightly different intensities of the reference-and sample beam as well as imperfect overlapping of the beams at the detector.
To record spectra of I des and φ, a simple program on a personal computer was used to tune the EC-QCL's wavelength across the desired range stepwise. At each step, the BaF 2 -window, mounted in a piezo-motor driven kinematic mount (Newport Agilis), was tilted in a closed loop to minimize the signal I des at the MCT detector. The closed loop differential controller was implemented in software to step the window's angle by a given value for a predefined number of times. Each step's direction was determined by the change in destructive signal upon the previous step to move towards smaller I des . The smallest recorded value of I des as well as I con and its corresponding tilt-signal was stored.
The EC-QCL, developed and manufactured at Fraunhofer IAF Freiburg, provided a tuning range from 1080 to 1370 cm −1 in pulsed mode with an average power of 250 mW within a 120-ns pulse at 1200 cm −1 [23]. Spectra were recorded by tuning the EC-QCL in steps of 0.4 cm −1 at a rate of 0.2 steps/s. The acquisition speed was limited by the slow closed loop for tilting the window, especially the actuator, which was by no means optimized towards speed. I des as well as the attenuated second output of the interferometer, I con , was measured by thermoelectrically cooled MCTdetectors (PVI-10.6, Vigo Systems S.A.) which were read out by a boxcar integrator at the laser repetition frequency of 9 kHz (limited by the maximum frequency of the gated integrator) and averaged over 300 pulses. I con was used to reference I des to the laser power and the reference absorption spectrum by calculating I des /I con . This is only valid for small absorptions and sufficient to demonstrate the characteristics of the interferometric setup within the context of this work. For a rigorous calibration, a separately recorded spectrum of I ref must be used to calculate A from I des using (5).
Reference absorption spectra were recorded on a Bruker Vertex 80v-FTIR spectrometer.

Chemicals and samples
Samples of acetone in cyclohexane at varying concentrations were prepared from reagent grade chemicals. Methanol added to some of the samples was of the same grade.

Simultaneous measurement of absorptionand refractive index spectra
The simultaneously measured spectra of I des /I con and the BaF 2 window's tilt of 2.4 mmol/L acetone in cyclohexane are depicted in Fig. 3. The spectra were smoothed using a Savitzky-Golay filter to remove fringes that could be assigned to the liquid cells. Inverse fast Fourier transformation of the raw spectra revealed multiple etalons of round trip optical path lengths between 4.5 and 7.5 mm, well separated from the signal spanning out from DC to 0.2 mm (FWHM of studied absorption profile is 8 cm −1 ). The BaF 2 window's tilt resembles the sigmoid anomalous dispersion profile of the refractive index. Kramers Kronig transformation yields an expected change in refractive index (peak to trough) of Δn max,min ≈ 7 × 10 −5 , corresponding to a phase difference of 6°, around the molecular resonance.

Dilution series of acetone in cyclohexane
Absorption spectra of different concentrations of acetone in cyclohexane were recorded. To demonstrate that no second detector is required to derive a reference for normalization, the interferometric spectra presented in Fig. 4 were normalized as I des − I des (c = 0) ∕I des (c = 0) . I des (c = 0) was recorded for pure cyclohexane in the sample-and reference position. Although this signal should be zero according to the theory derived above, due to limited modulation efficiency, I des (c = 0) represents a background spectrum that was found to be a good measure for I ref .
The normalized spectra of the recorded intensity I des in Fig. 4 qualitatively exhibit the expected behaviour of the interferometric signal discussed above with the sensitivity clearly increasing for larger concentrations, i.e. larger absorption A. After converting concentration to absorption via A = c × 0.0359 L/mmol , the parameters a 1 and a 2 of the function I des = a 1 × 1 + 10 −(A+a 2) − 2 × 10 −(A+a 2 )∕2 were fitted to the signal at the absorption maximum (Fig. 4, bottom). a 2 represents an offset absorption that resembles the small imbalance in the intensities of the reference-and sample beam.

Increased reference intensity
As discussed above, the absorptive signal exhibits its smallest sensitivity dI des /dA when the intensities of the reference and the sample beam are equal (A = 0) and increases as the two are unbalanced (increasing A, compare Fig. 2). To demonstrate the effect of unbalancing the signal and reference intensity, varying amounts of a spectrally flat absorber (methanol) were added to a sample of 3 mmol/L acetone in cyclohexane and their interferometric spectra were recorded with cyclohexane as a reference (Fig. 5, top left). The absorption spectrum of 3 mmol/L acetone in cyclohexane, recorded on an FTIR spectrometer, was used to replace the wavenumber axis of the interferometric spectra by an absorption axis (compare Fig. 5). The so derived curves of I des /I con versus A were analysed by fitting second-order polynomials to them. The fitted coefficients are given in Table 1. The coefficient K 1 , representing the sensitivity of the absorptive signals for small A, is increased by a factor of seven by reducing the intensity of the sample beam by 11 percent. At the same time, the background-signal (K 0 ) is increased by a factor of five, giving rise to a trade-off between background intensity and sensitivity. As discussed in detail in Sect. 6, this trade-off must be made based on the sources of noise in the system.
As discussed above, the noticeable offset of I des /I con for I s = I ref and A = 0 in Fig. 5 arises from the limited modulation efficiency M and the slightly differing intensities in the two interferometric arms (Fig. 7).

Modelling limits of detection
Experimental shortcomings of the current setup (fringes, poor phase adjustment, limited ruggedness) prohibit a direct experimental assessment of figures of merit as compared to established direct absorption spectroscopy (e.g. FTIR). Here, we give a comparative discussion based on the equations presented above and a noise model including shot noise, detector noise and power noise of the light source. For the Mach-Zehnder interferometer, noise associated with the interferometric phase is considered negligible, as can be expected for a rugged design and well-optimized feedback loop that keeps the phase at π. For both, direct absorption and interferometric measurements, fringes are indistinguishable from absorption (and dispersion) in the equations and effect the measurements equally. The noise equivalent power (NEP), given in WHz −1/2 , overlying the recorded signals, is modelled as Herein, P shot , NEP det and P PN are shot noise, the noise equivalent power of the detector and power noise from the laser (or other light source), respectively. Shot noise is given as  Laser power noise P / P Hz with h representing Planck's constant and ν the optical frequency. The power of the detected light P sig is given by (compare (4)) for the interferometric signal and by for direct absorption, with P 0 representing the laser power.
A off represents an intentionally introduced imbalance between the intensities of the destructively interfering fields (compare Sect. 5.3). P PN is derived by considering the effect of normally distributed laser power noise of power P N on the recorded signals P sig for A = 0.
NEP det , often specified inversely as detectivity D * , usually scales with the power of light to be detected. Typical NEP det for Mid-IR detectors are on the order of 10 −6 Hz −1/2 of the saturation power. Assuming an amplification chosen to match the highest expected signals, NEP det is given as From here on, we assume A max = 0.2 (corresponding to P sig,int = 10 −2 P 0 ) for the strongest interferometric signal. a is assumed 10 −6 Hz −1/2 .
To calculate limits of detection (LOD), we employ The sensitivities dP sig /dA read The so estimated LOD of the interferometric absorption measurement for three different offset absorptions A off as well as the LOD of the direct absorption measurement are plotted in Fig. 6 against the laser intensity noise P N . For small P N , signals are detector noise limited and independent from P N . With an offset absorption A off = 0.1, the interferometric measurement yields slightly smaller LODs due to the reduced detector noise associated with the reduced dynamic range of the detector. As P N increases, the interferometric  measurement yields smaller LODs. This is due to the interferometric referencing of the light from the sample arm of the interferometer to the reference arm, which cancels common mode noise. As A off increases, the increased sensitivity dP sig,int ∕dA yields lower LODs in the detector noise limited regime, but higher LODs for increased laser power noise. This is explained by incomplete destructive interference in the unbalanced interferometer and hence only partial removal of power noise.
Assuming no additional noise sources for the interferometric phase, LODs for the refractive index can be estimated based on the noise floor of the signal since this signal's minimum with respect to φ defines the refractive index (see (2)). An estimate of LODs equivalent to (15) using the sensitivity dP sig,int /dφ is prohibited by the vanishing sensitivity at φ = 0. Instead, we propagate the error NEP of the P sig,int measurement to the measurement of phase by solving NEP = P sig,int (Δ ) − P sig,int (0) for Δφ setting A = 0.
The measurement uncertainty Δφ and resulting Δn are plotted in Fig. 7 against laser power noise P N using the same parameters for NEP as in Fig. 6, a bandwidth BW of 1 Hz, a liquid film thickness of 1.9 mm and a wavelength of 8.3 µm. For a comparison, refractive index and phase can be referenced to absorption using Kramers-Kronig transformation and assuming an isolated Gaussian absorption profile of peak absorbance A. For this profile, Kramers-Kronig transformation yields an anomalous dispersion profile spanning, from top to bottom, n max − n min = A∕d × 1.86 × 10 −3 mm −1 , or max − min = A × 80 • ∕AU . For the Mach-Zehnder interferometer, this comparison yields a much better performance of absorptive measurements, which is a result of the vanishing sensitivity dP sig,int /dφ at φ = 0.

Conclusions
In this contribution, we introduced an interferometric setup that facilitates the simultaneous measurement of absorption and refractive index spectra of liquid samples. We discussed the theory relating the recorded signals to absorption and dispersion and presented proof of principle measurements supporting theoretical predictions.
For the absorptive signal, the theoretical considerations and experimental data demonstrate important differences to (18) P sig,int ( ) = P 0 RT 1 + 10 −(A+A off ) − 2 × 10 − A+A off 2 cos ( ) , � direct absorption spectroscopy. Spectra are recorded close to background free and are inherently referenced to the sample (e.g. solvent or matrix) introduced into the reference position. Although the sensitivity, i.e. change of recorded signal with absorption or concentration, is reduced as compared to direct absorption spectroscopy, this does not imply higher limits of detection (LOD). The latter depend on the noise floor, which, assuming the same sources of noise, is very different for the two techniques. While the reduced sensitivity (smaller signals) makes interferometric absorption measurements more prone to detector noise, they are much less susceptible to intensity noise due to the reduced levels of background intensity at the detector. This is of particular interest for external cavity laser based spectroscopy, since these lasers are known for intensity noise associated with mode hops [10,24]. A trade-off between sensitivity and immunity to intensity noise can be made by unbalancing the intensities of the reference and sample beam.
So far, an experimental investigation of these interesting properties was hampered by fringes in the sample cells as well as the slow actuator and feedback loop for the phase adjustment. Using wedged and/or anti reflection-coated windows could reduce the fringes. Adjustment of the phase can be greatly speeded up by using one of the interferometer's mirrors on a piezo element as the phase actuator [25].