Optical frequency comb Faraday rotation spectroscopy
Alexandra C. Johansson Jonas Westberg Gerard Wysocki Aleksandra Foltynowicz
physics.optics



















A. C. Johansson et al. 

Aleksandra Foltynowicz

		aleksandra.foltynowicz@umu.se

		

		Gerard Wysocki

		gwysocki@princeton.edu

		

		†These authors contributed equally to this work.

		

		[1]Department of Physics

		Umeå  University

		901 87 Umeå , Sweden

		

		[2]Department of Electrical Engineering

		Princeton University

		Princeton, New Jersey 08544, USA
                         







Optical frequency comb Faraday rotation spectroscopy



    Alexandra C. Johansson[1],†Jonas Westberg[2],†Gerard Wysocki[2]
    Received: 25 Janauray 2018 
===================================================================










We demonstrate optical frequency comb Faraday rotation spectroscopy (OFC-FRS) for broadband interference-free detection of paramagnetic species. The system is based on a femtosecond doubly resonant optical parametric oscillator and a fast-scanning Fourier transform spectrometer (FTS). The sample is placed in a DC magnetic field parallel to the light propagation. Efficient background suppression is implemented via switching the direction of the field on consecutive FTS scans and subtracting the consecutive spectra, which enables long term averaging. In this first demonstration, we measure the entire Q- and R-branches of the fundamental band of nitric oxide in the 5.2-5.4 range and achieve good agreement with a theoretical model.











§ INTRODUCTION


Laser-based spectroscopy for quantitative gas detection offers high sensitivity, high species specificity and the possibility of performing real-time measurements in situ without sample consumption. Laser-based gas sensors have found applications in fields ranging from fundamental science <cit.> and environmental monitoring <cit.> to medical diagnostics <cit.> and security <cit.>. However, the inherently limited optical bandwidth of these sensors has mostly restricted them to single species detection in spectral regions that are free from spectral interferences. The development of optical frequency combs has largely addressed this issue by massively extending the attainable spectral bandwidth, while retaining the high resolution <cit.>. Contemporary, commercially available, frequency comb systems now span hundreds of wavenumbers with optical power distributed across tens of thousands of equidistant comb lines. These sources have opened up for broadband, sensitive, and high-resolution spectroscopy, which allows retrieval of immense amounts of spectroscopic information in short acquisition times <cit.> and makes multi-species detection with a single coherent light source feasible <cit.>. Computational methods involving baseline removal together with multi-line fitting can be straightforwardly implemented to retrieve quantitative information from these multi-species spectra if accurate spectroscopic models of each species are known <cit.>. However, quantifying the concentrations of each constituent in a congested absorption spectrum is a challenging task when sufficiently accurate spectral models for the background constituents are unknown and their spectra cannot be assigned. The difficulties may be further exacerbated by drifting optical fringe backgrounds that prohibit long-term averaging. Congested absorption spectra are commonly encountered in e.g. combustion diagnostics, and are especially troublesome in wavelength regions where sufficiently detailed spectral models for water vapor are not yet available <cit.>.

In case of paramagnetic species any potential spectral interference can be effectively suppressed by implementing the interference- and background-free Faraday rotation spectroscopy (FRS) technique <cit.>, which relies on probing the opto-magnetic properties of paramagnetic species subjected to an external magnetic field. In this work, we merge the spectral bandwidth of an optical frequency combs with FRS to obtain a broadband high-resolution spectroscopy system that can provide simultaneous access to the entire ro-vibrational band of the target paramagnetic species. By placing the sample inside a solenoid (providing the external magnetic field) located between two polarizers, the polarization rotation induced by the molecules is translated into an intensity change that is directly proportional to the molecular number density. The direction and magnitude of the magnetic field can be controlled through the current of the solenoid, which allows for direct modulation of the Faraday rotation signal. The insusceptibility of diamagnetic species, such as H2O and CO2, to the Faraday effect infers that interferences from these species are efficiently suppressed (with the caveat that a sufficient number of photons still reach the photodetector). The contributions from the paramagnetic species can therefore be disentangled from a congested absorption spectrum, which significantly simplifies the spectroscopic assignments. Magnetic field modulation also eliminates the transmission baseline, thereby allowing prolonged averaging times even outside a temperature controlled laboratory environment. 

In this first demonstration of optical frequency comb Faraday rotation spectroscopy (OFC-FRS) we detect interference-free spectra of nitric oxide (NO) using a system based on a femtosecond doubly resonant optical parametric oscillator (DROPO) and a fast-scanning Fourier transform spectrometer (FTS). We measure the spectrum of the entire Q- and R-branches of NO at 5.2-5.4, and the measurement shows good agreement with a theoretical model.





§ THEORY


A general overview of the light propagation in an FRS system with DC magnetic field is shown in Fig. <ref>. The polarization of the incoming light is oriented horizontally (x-axis) as indicated by the red arrow. The magnetically induced polarization rotation due to the paramagnetic absorber is converted to an intensity change using a polarizer (analyzer) after the sample cell. The choice of the analyzer angle depends on the noise properties of the system, where small angles are used for systems with large relative intensity noise and 45 is used for a detector noise limited system <cit.>.




Consider a horizontally polarized (x-direction) electric field impinging on the sample, where a small ellipticity, ε , has been added to account for the finite extinction ratio of the polarizer <cit.>. The electric field before the interaction with the sample can be written in Jones’ notation <cit.> as

    E_i = E_0 [  1; iε ]e^i(kz-ω t),

where E_0 is the electric field amplitude, k the wave vector, ω the angular frequency, and t the time. This expression can be rewritten as a sum of the left-hand (LCP) and right-hand (RCP) circular components expressed in terms of their unit vectors

    ê_L = 1/√(2)[ 1; i ],    ê_R = 1/√(2)[  1; -i ],

which gives,

    E_i  = E_0/√(2)e^ik_0L-ω t[ê_L(1-ε)e^-δ_L + iϕ_L + ê_R(1+ε)e^-δ_R + iϕ_R],

where L is the interaction length, k_0 is the wave vector in the absence of absorbers, and  δ_L,R and ϕ_L,R  denotes the attenuation and phase shift experienced by the circular components of the electric field upon interaction with the absorber. δ_L,R and ϕ_L,R are related to the extinction coefficient, κ, and the refractive index, n, through

    δ_L,R = κ_L,Rω L/c,    ϕ_L,R = (n_L,R-1)ω L/c.

After interaction with the paramagnetic sample subjected to a magnetic field oriented along the light propagation axis, the light is transmitted through a second polarizer (analyzer) oriented at an uncrossing angle, θ, with respect to the first polarizer. The Jones’ matrix for this element is given by <cit.>

    J_pol =
       [
      asin^2(θ) + bcos^2(θ)
    
     asin(θ)cos(θ)-bsin(θ)cos(θ)
    .
               . asin(θ)cos(θ)-bsin(θ)cos(θ) 
    acos^2(θ) + bsin^2(θ)
    ],

where a and b denote the fractions of the electric field transmitted along the main and orthogonal polarization axis of the analyzer. For a high extinction ratio analyzer, a≈ 1 and b ≪ 1. Finally, the transmitted intensity, I_t, is given by the product of the transmitted electric field and its complex conjugate, which can be calculated from Eqs. (<ref>) and (<ref>). This gives

    I_t    = I_0/2e^-2δ̅{(a^2+ε^2b^2)[cosh(Δδ) - cos(Δϕ + 2θ)]
       + 2(a^2 + b^2)εsinh(Δδ)
       + (ε^2 a^2 + b^2)[cosh(Δδ) + cos(Δϕ+2θ)]},

where δ̅ is the average of the attenuations experienced by the LCP and RCP light, Δδ is the differential attenuation due to magnetic circular dichroism (MCD), and Δϕ is the differential phase shift due to magnetic circular birefringence (MCB). These entities are in turn given by

    2δ̅   = δ_L + δ_R = α_0L/2(χ^abs_L + χ^abs_R) = α̅ L 
    Δδ   = δ_L - δ_R = α_0L/2(χ^abs_L - χ^abs_R) = 2Φ
    Δϕ   = ϕ_L - ϕ_R = α_0L/2(χ^disp_L - χ^disp_R) = 2Θ,

where α_0 is the on-resonance absorption coefficient given by α_0 = S'c_relp, where S' is the integrated gas linestrength (cm-2/atm)
, c_rel is the relative concentration, p the total gas pressure (atm), and α̅ the average absorption coefficient. χ^abs,disp_L,R are the absorption and dispersion lineshape functions for the LCP and RCP components given by <cit.>

    χ^abs,disp_L,R = 3∑_M”_JM'_J[    J'     1    J”; -M'_J Δ M_J  M”_J ]^2 χ_V^abs,disp(ν_M”_JM'_J),

where the ∑(...)^2 term is the Wigner 3-j symbol for a M”_J → M'_J transition, and Δ M_J =  M'_J - M”_J (equal to +1 or -1 for LCP and RCP, respectively). χ_V^abs,disp(ν_M”_JM'_J) denotes the Voigt absorption and dispersion lineshapes and ν_M”_JM'_J  is the detuning frequency shifted by the external magnetic field, given by

    ν_M”_JM'_J = [ν - ν_0 - (M'_Jg'_J - M”_Jg”_J)μ_B B]√(ln(2))/δν_D,

where ν is the frequency of the light, ν_0 is the unperturbed molecular transition frequency, g_J the Landé g-factor, μ_B the Bohr magneton, B the external magnetic field, and δν_D the Doppler width (FWHM) of the transition. 

The Faraday rotation angle, Θ, is defined as half of the induced differential phase shift, i.e. Θ = Δϕ/2 [see Eq. (<ref>)]. Since the sign of Θ depends on the direction of the magnetic field <cit.>, efficient background suppression can be achieved by alternating the direction of the magnetic field and subtracting the spectra measured with opposite fields. Using Eq. (<ref>), the transmitted intensities for the two field directions, denoted I_t^+ and I_t^-, where the + and - signs indicate the direction of the magnetic field, can be written as

    I_t^±   = I_0/2e^-α̅ L{(a^2+ε^2b^2)[cosh(± 2Φ) - cos(± 2Θ + 2θ)]
       + 2(a^2 + b^2)εsinh(± 2Φ)
       + (ε^2 a^2 + b^2)[cosh(± 2Φ) + cos(± 2Θ+2θ)]}.

It follows that their average, I_t^0 = (I_t^++I_t^-)/2, is given by

    I_t^0   = I_0/4e^-α̅ L{(a^2+ε^2b^2)[2cosh(2Φ) - 2cos(2θ) cos(2Θ)] 
       + (ε^2 a^2 + b^2)[2cosh(2Φ) + 2cos(2θ)cos(2Θ)]},

which for small Φ and high extinction ratio polarizers (ε^2b^2 ≪ε^2a^2 ≈ b^2 ≪ a^2 ≃ 1) reduces to

    I_t^0   = I_0/2e^-α̅ L[1 - cos(2θ)cos(2Θ) ].

When the polarizer uncrossing angle θ is set to 45, this simply gives the Beer-Lambert law with half of the power transmitted through the analyzer. The difference of the transmitted intensities measured with opposite magnetic field directions, Δ I_t^± = I_t^+ - I_t^-, is given by

    Δ I_t^±   = I_0 e^-α̅ L{sin(2θ)sin(2Θ) [(a^2+ε^2b^2) -  (ε^2 a^2 + b^2)] 
       +2(a^2 + b^2)εsinh(2Φ)},

which for high extinction ratio polarizers reduces to 

    Δ I_t^±   = I_0 e^-α̅ L[ sin(2θ)sin(2Θ) +2εsinh(2Φ)].

A normalized FRS signal can be obtained by taking the ratio of the differential intensities and their mean, i.e. 

    S_FRS =Δ I_t^±/I_t^0≈2sin(2θ)sin(2Θ) +4εsinh(2Φ)/1-cos(2θ)cos(2Θ),

where the last equality assumes high extinction ratio polarizers. The MCD contribution Φ [see Eq. (<ref>)] is coupled through the extinction ratio of the first polarizer, ε, which introduces an asymmetry in the lineshape <cit.>. Note that the dependence on I_0e^-α̅ L is removed by the normalization procedure of Eq. (<ref>), which infers that any variations in intensity that are common for the two signals are suppressed. This includes interferences from diamagnetic species, such as water vapor, and from slowly varying optical fringes affecting the transmission baseline.


For the special case of a system limited by detector noise, the highest signal-to-noise ratio is obtained when the uncrossing angle is set to 45  <cit.>. This further simplifies the expression above to

    S_FRS = 2sin(2Θ) + 4εsinh(2Φ).

For small rotation angles and high extinction ratio polarizers the normalized FRS signal is directly proportional to the Faraday rotation angle, which can be used to extract the sample concentration <cit.>. For the cases where the aforementioned approximations are not justified, the full expressions of Eqs. (<ref>) and (<ref>) must be used to model the FRS spectrum.







§ EXPERIMENTAL SETUP AND PROCEDURES


The experimental setup, shown in Fig. <ref>(a), is based on a femtosecond doubly resonant optical parametric oscillator (DROPO), two polarizers, a gas cell placed inside a DC-solenoid, and a fast-scanning Fourier transform spectrometer (FTS). A detailed description of the light source is provided in Ref. <cit.>. The DROPO is synchronously pumped with a mode-locked Tm:fiber laser with a repetition rate of 125. The idler of the DROPO operating in the degenerate configuration is tuned to cover 5.1-5.4 (1850-1960) with a total optical power of 5 by adjusting the length of the DROPO cavity. A long-pass filter after the DROPO is used to filter out the pump and signal light. The polarization of the idler light is cleaned by a Wollaston prism placed in front of a 17.5-long gas cell equipped with uncoated, tilted, and wedged CaF_2 windows. The cell is placed inside a 12.7-long DC-solenoid, which provides an axial magnetic field of 260±7G either parallel or antiparallel to the light propagation depending on the direction of the supplied current. The cell is operated with a continuous flow of 1% NO diluted in N_2 at a total pressure of 100Torr. The flow and pressure are actively controlled through a combination of flow and pressure controllers (providing 1.3Torr precision on the gas pressure). The magnetic field strength was determined from the intensity of the current supplied to the solenoid using a calibration curve obtained earlier by direct measurements with a Gaussmeter. A wire-grid polarizer, with the a and b parameters measured to be a = 0.67, and b = 0.40, is placed after the cell and rotated by 45 with respect to the Wollaston prism to convert the polarization rotation induced by the interaction with the sample to an intensity change. The transmitted light is measured with a fast-scanning FTS, 



illustrated in Fig. 2(b), equipped with two HgCdTe detectors in balanced configuration <cit.>. The two out-of-phase outputs from the FTS are subtracted, which allows for common-mode noise suppression (and higher signal-to-noise ratios). An interferogram with a spectral resolution of 250 is measured in 3.5 and a stable continuous-wave laser at 1.56 is used for optical path difference calibration. Consecutive interferograms are collected with opposite magnetic field by synchronizing the direction of the supplied solenoid current with the movement of the FTS mirrors. The corresponding spectra, I_t^±, are retrieved from fast Fourier transform of the interferograms. The normalized FRS spectrum is obtained by subtracting spectra recorded with opposite magnetic field, Δ I_t^± = I_t^+ - I_t^-, and normalizing to their mean value, I_t^0 = (I_t^+ + I_t^-)/2, which makes the FRS spectrum baseline- and calibration-free. 







§ RESULTS


Figure <ref>(a) shows an average of 500 transmission spectra measured with opposite magnetic field, I_t^+ in red and I_t^- in green, and their mean, I_t^0, in blue, zoomed in on the ^2Π_3/2(ν=1→ 0) Q-branch of NO at 1875.8. The rotation of the polarization of the light due to the magnetic field is manifested as either an increase or a decrease of intensity depending on the direction of the field intensity (compared to the case with no magnetic field, which corresponds to the mean value). Figure 3(b) shows the resulting normalized FRS spectrum. Subtraction of the consecutive transmission spectra removes any structure in the baseline, which can clearly be seen by comparing the upper and lower panels of Fig. <ref>. 



Figure <ref> further demonstrates the water and optical fringe suppression capabilities, where (a) displays a part of the transmission spectrum (500 averages) and (b) shows the resulting normalized FRS spectrum. This region is visibly affected by optical fringes, originating from the analyzer, and by water absorption, caused by the presence of atmospheric water in the unpurged beam path between the DROPO and the FTS. The water absorption reaches almost 20% for the strongest transition. However, the diamagnetic water molecules do not induce any polarization rotation and hence two consecutive transmission spectra (with opposing magnetic fields) exhibit nearly perfect overlap as shown in Fig. <ref>(a). Figure <ref>(b) shows the corresponding normalized FRS spectrum, in which the water absorption and optical fringes are efficiently cancelled. The noise on the baseline over this range is ∼ 5×10^-4, and is primarily limited by the detector noise. The normalized FRS spectrum does not show any significant structure remaining from the water absorption, however there is a small residual baseline caused by optical fringes that do not cancel out completely due to baseline drift that occurs on a time scale faster than that required to perform two FTS scans. The fluctuation of the baseline caused by these residual fringes is ∼ 2× 10^-3, which corresponds to a reduction of optical fringes observed in the transmission spectrum by a factor of ∼ 200. It should also be noted that the measurement is possible only if sufficiently many photons reach the photodetectors. Therefore, in spectral regions where water absorption approaches 100% the system will exhibit a severe degradation of the signal-to-noise ratio. In this case, water vapor removal by optical path purging is still recommended although the demand on its efficiency is relaxed. 


 
Figure <ref>(a) shows the normalized FRS spectrum of 1% NO in N_2 at 100Torr, covering 1850-1930 (500 averages). The residual baseline caused by uncorrelated structure between two transmission spectra has been corrected in the post-processing, and water lines with 100% absorption have been masked out. Figure <ref>(b) shows a zoomed view of the normalized FRS spectrum for the Q-branch (black) together with a fit (red). The strongest normalized FRS signal comes from the Q(3/2) transition of the ^2Π_3/2 subsystem, and higher J-quantum numbers yield diminishing signal strengths. Figure <ref>(c) displays the normalized FRS spectrum (black) for two lines in the R-branch, the R(13/2) and the R(15/2) transitions of the ^2Π_1/2 and ^2Π_1/2 subsystems, respectively, together with a fit (red).
 
The fitted model is based on the full expressions of Eqs. (<ref>) and (<ref>), to account for imperfect polarizers <cit.>. The spectral line parameters are taken from the HITRAN database <cit.> and the Landé g-factors are calculated from Refs. <cit.>. For the Q-branch in Fig. <ref>(b) the fitting parameters are the NO concentration, c_NO, the magnetic field strength, B, and the unbalancing term between the right-hand and the left-hand circularly polarized components for the first polarizer, ε. The values returned from the fit are c_NO = 1.055(2 )%, B = 228.8(5)G, and ε = 0.0507(3). The discrepancy from the measured peak field intensity of 260G is likely due to a non-uniform field distribution over the interaction length, and the frequency dependence of the polarizer parameters a and b that is not taken into account. For the fit to the R-branch, B and ε were fixed to the values returned from the fit to the Q-branch and only the NO concentration was used as a fitting parameter. Here, the fit returned c_NO = 0.944(4)%. The residuals in Figs <ref>(b) and (c) indicate that the model agrees well with the measured data, although some discrepancies can be observed. These are likely caused by experimental limitations (poor extinction ratio of the wire-grid analyzer and non-linearities in the FTS detectors), as well as imperfections in the theoretical model (i.e. the use of the Voigt lineshape function and uncertainties in the values of the Landé g-factors <cit.>).



 
It should be noted that the optimum magnetic field strength (i.e. the one that maximizes the signal-to-noise ratio) varies for the different branches due to their different Landé g-factors. For the NO spectrum, the field strength needed to maximize the signal from the R-branch is higher than that needed to maximize the signal from the Q-branch. The magnitude of the magnetic field used in the experiment was set to the technical maximum that could be sustained for long measurement durations (limited by the heat dissipation of the coil). At this field, the R-branch transitions were strongly undermodulated, while the Q-branch transitions were slightly overmodulated. In general, the magnetic field and sample gas pressure should be optimized for each FRS measurement by taking into account the difference in susceptibility to the magnetic field of the branches, magnetic field limitations, and the noise sources in the system. 

The stability of the system is characterized by an Allan-Werle plot, where the standard deviation of the fitted concentration as a function of averaging time is calculated. The concentrations are extracted from multiline fits to consecutive individual FRS spectra of the Q-branch, with a model based on Eqs. (<ref>) and (<ref>). The resulting Allan-Werle plot is shown in Fig. <ref>, and indicates a white-noise limited behavior for the full durations of the measurement (>1000 seconds). As mentioned previously, this is not unexpected given the fact that the normalized FRS spectrum is based on the difference between two consecutive transmission spectra, which is equivalent to time-multiplexed differential detection, albeit with long acquisition times. Based on the Allan deviation presented in Fig. <ref>, the 1σ minimum concentration detection limit of the system was estimated to 1.4ppm·m after 1000 seconds of integration time.









§ CONCLUSIONS


In conclusion, we presented for the first time optical frequency comb Faraday rotation spectroscopy (OFC-FRS) by measuring the spectrum of the entire Q- and R-branches of the fundamental vibrational band of NO at ∼5.3 using a femtosecond DROPO and a fast-scanning FTS. Switching the direction of the magnetic field for consecutive FTS scans and subtracting the resulting transmission spectra enables efficient background suppression and allows measurements over long averaging times. Moreover, the normalization by the mean of consecutive spectra provides a calibration-free signal. A theoretical model of the normalized FRS spectrum is presented, which shows good agreement with the measurements. A minimum detection limit for NO of 1.4ppm·m at 1000 is achieved, and the system remains stable for more than 1000 seconds. Further improvements in sensitivity can be obtained by improving the extinction ratio of the polarizers or increasing the interaction length through multiple passes <cit.> or by cavity enhancement <cit.>. While in this first demonstration we used an FTS to acquire broadband OFC-FRS spectra, the technique is compatible with other detection schemes of comb spectroscopy, for example with continuous Vernier filtering <cit.>, which should allow faster acquisition and thus faster modulation of the magnetic field. 

The broadband OFC-FRS technique may be useful for fundamental science applications, such as assessing Landé g-factors, where the calibrated frequency axis and immunity to instrumental lineshapes <cit.> are advantages compared to FRS based on tunable cw lasers in the mid-infrared <cit.>. The OFC-FRS technique is also applicable when targeting trace amounts of gaseous paramagnetic species masked by excessive quantities of interfering diamagnetic compounds (H_2O, CO_2, etc.). Such conditions are common in combustion diagnostics, where the large optical bandwidth provided by the comb in combination with the Faraday rotation technique will provide a new tool for probing high temperature chemical reactions involving magnetically sensitive species, e.g. free radicals such as NO, NO_2, OH and HO_2. Also in such applications any reduction of the spectral interference clatter from diamagnetic species coupled with the broadband coverage of the OFC-FRS technique opens up possibilities for measurements requiring access to multiple transitions simultaneously (e.g. precise multi-line thermometry). 


The authors thank Amir Khodabakhsh for help with operating the DROPO. The work at UmU was supported by the Knut and Alice Wallenberg Foundation (KAW 2015.0159) and the Swedish Foundation for Strategic Research (ICA12-0031). The work at Princeton was supported by the DARPA SCOUT program (W31P4Q161001).   






applphysb





