Disclosed is a linear absorption spectrometer to optically determine an absolute mole fraction of radiocarbon in a sample, the linear absorption spectrometer comprising: a laser light source that provides mid-infrared laser light for linear absorption by the radiocarbon in the sample; a high finesse optical resonator that is actively stabilized in a resonance frequency and comprising: a first supermirror comprising a first radius of curvature that provides cavity ring down reflection and that receives the mid-infrared laser light; a second supermirror comprising a second radius of curvature that provides cavity ring down reflection, the second supermirror in combination with the first supermirror comprises a relative difference of refractive index Δn/n from 1×10-8 to 6×10-6, such that the second supermirror transmits cavity ring down light from communicating the mid-infrared laser light through the sample in a sample cell; the sample cell interposed between the first supermirror and the second supermirror to contain the sample, the sample cell operating at a temperature from 180 K to 300 K during linear absorption of the mid-infrared laser light by the radiocarbon and comprising: a linear absorption optical path length that is greater than a kilometer (km); a first zero-pressure difference mirror mount on which the first supermirror is disposed and mechanically coupled to the sample cell; a second zero-pressure difference mirror mount on which the second supermirror is disposed and mechanically coupled to the sample cell; an optical switch interposed between the laser light source and the high finesse optical resonator such that the optical switch receives the mid-infrared laser light from the laser light source, modulates the mid-infrared laser light, and communicates modulated mid-infrared laser light to the first supermirror of the high finesse optical resonator; a photoreceiver in optical communication with the high finesse optical resonator and that receives the cavity ring down light from the second supermirror, the photoreceiver comprising a noise equivalent power that is less than a shot noise limit of the cavity ring down light, the linear absorption spectrometer providing the absolute mole fraction of the radiocarbon in the sample for the absolute mole fraction being from 1 part-per-quadrillion to 2.5 parts-per-trillion of radiocarbon in the sample.
The following descriptions should not be considered limiting in any way. With reference to the accompanying drawings, like elements are numbered alike.
A detailed description of one or more embodiments is presented herein by way of exemplification and not limitation.
It has been discovered that a linear absorption spectrometer herein provides optical detection of radiocarbon (14C) in a sample (e.g., gaseous carbon dioxide (CO2) that includes a combination of carbon isotopes) with a sensitivity below modern picomole per mole levels (e.g., 14C/12C<10−12). The linear absorption spectrometer includes a high-finesse optical resonator with a sample cell that includes cryogenic cooling, refrigerant cooling, and is vacuum compatible. As a result, the linear absorption spectrometer provides determination of the amount of the radiocarbon in a stable, low temperature gas environment during optical analysis and reduces spectral interferences. Optical detection of the radiocarbon in the sample is via linear absorption, wherein nonlinear molecular absorption does not compromise a measured cavity decay time of signal from the sample cell. Also, it is contemplated that a laser diode in an absence of high-bandwidth phase locking or laser stabilization is used in some embodiments of the linear absorption spectrometer.
Conventional radiocarbon analysis is performed at accelerator mass spectrometry (AMS) facilities that are expensive to construct and whose operation involves a team of technicians and scientists. The linear absorption spectrometer significantly reduces cost and provides automation of spectral analysis to reduce the time for radiocarbon analysis. The elimination of current financial and analysis time costs provide availability for radiocarbon analysis across geoscience, biology, atmospheric chemistry, archeology, and the like.
In an embodiment, with reference to
In an embodiment, with reference to
In an embodiment, with reference to
According to an embodiment, with reference to
In an embodiment, with reference to
In the linear absorption spectrometer 100, laser light source 2 provides mid-infrared laser light 4. Mid-infrared light source 2 can be a continuous wave laser, wherein a wavelength of mid-infrared laser light 4 is selected for absorption of radiocarbon in the sample. In a particular embodiment, the wavelength of mid-infrared laser light 4 is selected such that radiocarbon absorbs mid-infrared laser light 4 while other chemical species that may be present in simple cell 24 do not absorb at that wavelength or weakly absorb at that wavelength as compared to the absorption coefficient for radiocarbon. In a certain embodiment, the wavelength of mid-infrared laser light 4 is selected for absorption by radiocarbon in a P-branch in a rovibrational absorption spectrum for radiocarbon, based on a quantum mechanical selection rule for a change in rotational quantum number J, wherein ΔJ=−1.
It is contemplated that the wavelength of mid-infrared laser light 4 is from 1 micrometers (μm) to 50 μm, specifically from 1 μm to 30 μm, and more specifically from 2 μm 20 μm. An average power of mid-infrared laser light 4 can be from 1 microwatts (μW) to 10 watts (W), specifically from 100 μW to 100 milliwatts (mW), and more specifically from 1 mW to 30 mW. Moreover, mid-infrared laser light 4 in the high finesse optical resonator comprises a fractional frequency stability that is greater than 1 in 108.
Laser light source 2 can include a semiconductor laser. Exemplary laser light sources 2 include a quantum cascade laser, difference frequency generation laser, fiber laser, optical parametric oscillator, vertical-cavity surface-emitting laser, and the like.
In an embodiment, reference laser light 42 is combined with mid-infrared laser light 2. Reference laser light 42 can include a mid-infrared optical comb. The mid-infrared optical comb can include a fractional frequency instability that is less than 1 in 1011. A wavelength of the mid-infrared optical comb can be from 1 micrometers (μm) to 50 μm, specifically from 1 μm to 30 μm, and more specifically from 2 μm 20 μm. An average power of the mid-infrared optical comb can be from 1 μW to 10 W, specifically from 100 μW to 100 mW, more specifically from 1 mW to 30 mW.
It is contemplated that mid-infrared optical comb 41 comprises optical cavities, fiber laser oscillators, solid-state laser oscillators, microresonators, saturable absorbers, laser pump diodes, semiconductor lasers phase modulators, intensity modulators, optical amplifiers, optical compressors, nonlinear frequency conversion crystals, nonlinear photonic devices, waveguides, waveplates, highly nonlinear fiber, and the like.
In some embodiments, reference laser 40 is a stabilized frequency reference such as a helium neon laser (HeNe laser). The HeNe laser can include a fractional frequency instability that is less than 1 in 1010. A wavelength of the HeNe laser can be from 600 nanometers (nm) to 4000 nm, specifically 633 nm. An average power of the HeNe laser can be from 1 μW to 1 W, specifically from 100 μW to 1 mW, more specifically from 100 μW to 2 mW.
Mid-infrared laser light 4 from laser light source 2 is received by optical switch 6. Optical switch 6 modulates mid-infrared laser light 4. Here, modulation of mid-infrared laser light 4 includes intensity modulation, frequency modulation, phase modulation, or a combination thereof. In an embodiment, optical switch 6 modulates a frequency of mid-infrared laser light 4 to produce modulated mid-infrared laser light 4 that are communicated to high finesse optical resonator 10. A repetition rate of modulated mid-infrared laser light 4 produced by optical switch 6 can be from 1 hertz (Hz) to 50 kHz, specifically from 5 hertz (Hz) to 10 kHz, and more specifically from 25 hertz (Hz) to 2 kHz.
Exemplary optical switches include acousto-optic modulators, electro-optic modulators, mechanical switches, fiber optic switches, and the like. In an embodiment, optical switch 6 includes a germanium acousto-optic modulator to modulate mid-infrared laser light 4 at a frequency of 60 MHz.
High finesse optical resonator 10 receives mid-infrared laser light 8 from optical switch 6, reference laser light 42 from reference laser 40, the sample that includes radiocarbon, a cooling fluid, or a combination thereof. With reference to
A combination of first supermirror 12, first zero-pressure difference mirror mount 26, second supermirror 16, and second zero-pressure difference mirror mount are selected to provide a high finesse for high finesse optical resonator 10, wherein a reflectivity of first supermirror 12 and second supermirror 16 independently is from 99.9% to 99.99999% at a wavelength of mid-infrared laser light 20. In an embodiment, reference laser light 42 is present with mid-infrared laser light 20 in high finesse optical resonator 10, wherein a reflectivity of first supermirror 12 and second supermirror 16 independently is from 10% to 99.7%, specifically from 50% to 99.7% at a wavelength of reference laser light 40.
First supermirror 12 and second supermirror 16 are selected to transmit and reflect mid-infrared laser light (8, 20, 22). Moreover, first supermirror 12 and second supermirror 16 independently include a radius of curvature that provides cavity ring down reflection in high finesse optical resonator 10. The radius of curvature can be from micrometers (μm) to hundreds of meters (m), specifically from 5 mm to 5 m, and more specifically from 35 cm to 1.5 m. Exemplary supermirrors (12, 16) include alternating dielectric stacks, metallic mirrors, silicon mirrors, calcium fluoride mirrors, zinc selenide mirrors, and the like. Supermirrors (12, 16) can include a plurality of coating layers that provide a selected reflectivity, birefringence, and the like. Second supermirror 16 in combination with first supermirror 12 includes a relative difference of refractive index Δn/n from 3×10−11 to 3×10−2, specifically from 1×10−8 to 3×10−4, and more specifically from 1×10−8 to 6×10−6.
With reference to
In an embodiment, secondary fluid conduit 73 is isolated from fluid communication with primary fluid conduit 34, wherein the sample is present in an absence of contact with the cooling fluid. In some embodiments, secondary fluid conduit 73 is in fluid communication with primary fluid conduit 34, wherein the sample can be combined with the cooling fluid.
Supermirror (12, 16) is disposed in zero-pressure difference mirror mount (26, 28). In this manner, supermirrors (12, 16) are not subjected to strain such that the respective birefringence of supermirrors (12, 16) are maintained without variation when disposed in zero-pressure difference mirror mounts (26, 28) or subjected to a change in pressure of sample cell 24. With reference to
With reference to
It is contemplated that linear absorption spectrometer 100 provides frequency stabilization of mid infrared laser light (8, 20, 22). The frequency stabilization of mid infrared laser light (8, 20, 22) can be accomplished mechanically, electronically, optically (e.g., see
A number (e.g., 1, 2, 3, or more) of spacer members 97 can be selected to stabilize high finesse optical resonator 10 and can be used for alignment thereof. A cross-sectional shape of spacer member 97 can be any shape effective for mechanical stabilization and can be, e.g., circular, polygonal, square, rectangular, ellipsoidal, and the like. Exemplary materials for spacer member 97 include low coefficient of thermal expansion plastic, glass, composite, ceramic, metal, or a combination thereof, and the like. In an embodiment, spacer member 97 includes a low coefficient of thermal expansion alloy such as iron-nickel (available under the tradename INVAR). A length of spacer member 97 can be from 10 μm to 1 km, specifically from 1 cm to 10 m, and more specifically from 75 cm to 3 m.
In linear absorption spectrometer 100, the sample is disposed in sample cell 24 of high finesse optical resonator 10 and subjected to radiation by midinfrared laser light 20. In an embodiment, the sample includes a gas. The gas can be from a petrogenic source, biogenic source, medical source, radioactive source, and the like, or a combination thereof. The sample can include radiocarbon, and the radiocarbon absorbs midinfrared laser light 20. As used herein, “radiocarbon” refers to carbon-14 (14C). Exemplary species that include radiocarbon in the sample include an inorganic oxide (e.g., CO2, CO and the like), a hydrocarbon (e.g., an alkane such as methane, an alkene such as ethylene, an alkyne such as acetylene, and the like), a hydrocarbon (e.g., an organic acid such as acetic acid, formic acid, and the like), an ether, an ester, a ketone, and the like. It should be appreciated that carbon in the sample can include radiocarbon alone or in combination with carbon-12 (12C) or carbon-13 (13C).
Other species present in the sample can include an atom or molecule that include a nitrogen, oxygen, phosphorus, sulfur, and the like. Exemplary contaminant species in the sample include a phosphorus oxide (e.g., PxOy, wherein x and y are integers and independently can be 0, 1, 2, and the like), nitrogen oxides (NxOy, wherein x and y are integers and independently can be 0, 1, 2, and the like), sulfur oxides (SxOy, wherein x and y are integers and independently can be 0, 1, 2, and the like), and the like.
Radiocarbon can be present in the sample from 100 attomol per mol (amol/mol) to 10 micromole (μmol/mol), and more specifically from 100 amol/mol to 10 pmol/mol, based on a total number of moles in the sample. It is contemplated that the absolute mole fraction of the radiocarbon in the sample is from 1 parts-per-quadrillion (ppq) to 2.5 parts-per-trillion (ppt).
The sample can be cooled via thermal transfer with the sample when cooling member 73 of sample cell 24. The cooling fluid disposed in coolant reservoir 75 of cooling member 73 can have a temperature from 2 K to 330 K, specifically from 150 K to 300 K, and more specifically from 180 K to 220 K. Exemplary cooling fluids include helium, nitrogen, ethylene glycol, nanofluids, polydimethylsiloxaneor a combination thereof, and the like.
Cavity ring down light 22 is produced in response to absorption of mid infrared laser light 20 by the sample. Photo receiver 30 receives cavity ring down light 22 from high finesse optical resonator 10 and produces absorption signal 32. Photo receiver 30 can be followed a photodiode, photomultiplier, photodiode array, thermopile, and the like. In an embodiment, photo receiver 30 includes the photodiode and electronic amplifiers and filters. It is contemplated that photo receiver 30 is cooled to a temperature at which photo receiver 30 noise equivalent power is less than five times a shot noise limit of cavity ring down laser light 22. In this manner, photo receiver 30 that is sensitively requires cavity ring down light 22. It is contemplated that photoreceiver 30 is arranged to provide an etalon-immune distance. As used herein, “etalon-immune distance” refers to a distance from the second supermirror reflective surface 18 to photo receiver 30 that is an integer multiple of distance from the second supermirror reflective surface 18 to first supermirror reflective surface 14.
In an embodiment, a process for making linear absorption spectrometer 100 includes providing laser light source 2; disposing optical switch 6 in optical communication with laser light source 2; providing high finesse optical resonator 10 in optical communication with optical switch 6 to receive mid infrared laser light 8 modulated by optical switch 6; and disposing photo receiver 30 in optical communication with high finesse optical resonator 10 to receive cavity ring down light 22 therefrom. Here, the process also can include disposing first super mirror 12 in first zero-pressure difference mirror mount 26; disposing second super mirror 16 in second zero-pressure difference mirror mount 28; disposing first zero-pressure difference mirror mount 26 on sample cell 24; and disposing second zero-pressure difference mirror mount 28 on sample cell 24.
Linear absorption spectrometer 100 has numerous beneficial uses, including performing linear absorption spectrometry on a sample, determining an absolute mole fraction and number density of radiocarbon in the sample, and the like. In an embodiment, a process for determining an absolute mole fraction and number density of radiocarbon in the sample includes adjusting the wavelength of mid-infrared laser light 4 to be transmitted by linear absorption spectrometer 100 and probing radiocarbon absorption; measuring midinfrared laser frequency by referencing to mid-infrared optical frequency comb 41 or the like; measuring the decay signal of midinfrared laser light 8 exiting linear absorption spectrometer 100 using photoreceiver 30 and digitizer 108; determining the decay rate by fitting exponential decay plus baseline offset model to measured decay signal; repeating decay rate measurement a plurality of times to improve measurement precision; repeating absolute laser frequency and cavity ring-down decay rate measurements over multiple wavelengths yielding spectrum of base cavity losses plus absorption losses; fitting parameterized, mathematical model of the absorption spectrum to the measured spectrum to yield peak area associated with midinfrared light absorption by radiocarbon, treating the spectrum as a linear superposition of base spectrometer losses having a constant, linear or quadratic wavelength dependence with residual, sinusoidal features caused by parasitic coupled-cavity effects (etalons), an absorption line profile for radiocarbon centered on radiocarbon transition frequency, additional line profiles of interfering species; calculating number density of the radiocarbon as ratio of fitted peak area for radiocarbon to line intensity of radiocarbon transition evaluated at sample temperature; calculating molar fraction of radiocarbon by dividing radiocarbon number density by sample number density based on measured sample pressure and temperature.
Linear absorption spectrometer 100 has numerous advantageous and beneficial properties. In an aspect, linear absorption spectrometer 100 provides ultrasensitive detection of modern levels (e.g., less than 2 ppt) of radiocarbon in the sample. Moreover, linear absorption spectrometer 100 overcomes presence of parasitic etalons because etalons have been minimized by of optical element materials, coatings, and relative position of normal-incidence optics, photoreceiver and the like. Further, averaging of spectra provides reduction in effect due to randomly varying etalons. Additionally, linear absorption spectrometer 100 includes simple cell 24 with cooling member 73 to provide an ultralow temperature to the sample to reduce spectral congestion due to contaminants in the sample.
The articles and processes herein are illustrated further by the following Examples, which are non-limiting.
Optical detection of radiocarbon via linear absorption below contemporary levels was performed by cavity ring-down spectroscopy in the linear absorption spectrometer. Petrogenic and biogenic samples of CO2 were distinguished by optical measurements.
Sub-picomole per mole (pmol/mol) detection of radiocarbon in samples of CO2 originating from the combustion of either petrogenic (i.e. fossil fuel) or biogenic samples was performed with a low-temperature (T=220K) cavity ring-down spectrometer pumped by a quantum cascade laser. Repeated spectroscopic measurements of a representative 86% biogenic sample yielded a standard error of 130 fmol/mol (1 part-per-quadrillion, or ppq, is defined as 1 fmol/mol) in two hours of laboratory time. Multi-hour measurements repeated over several days showed a standard deviation of 60 fmol/mol, thus validating the potential for high-fidelity optical detection using linear absorption spectroscopy.
The mid-infrared cavity ring-down spectrometer included a continuous-wave (cw) quantum cascade laser (QCL), an acousto-optic modulator (AOM) to act as a fast optical switch, a high-finesse (F=50000) optical resonator whose axis resides at the cross-hairs of a cylindrical cold finger, and a liquid-nitrogen-cooled InSb photoreceiver. Cavity time constants τ were digitized and fit at rates≈100 Hz with an electronic bandwidth of 500 kHz using a low-bandwidth transmission lock of the laser to a given cavity mode. Molecular spectra were acquired by tuning the QCL current and thus discretely jumping the laser frequency from one mode of the optical resonator to another mode. With a cavity free spectral range of FSR≈100 MHz, 1.1 GHz-wide spectra of CO2 were recorded in approximately 70 s. The spectral sampling grid defined by the optical resonator length was shifted in sub-FSR increments by tuning a piezo-electric transducer mounted to one of the cavity mirrors.
A representative compilation of interleaved, baseline-corrected spectra is shown in panel A of
Panel B of
An ultra-sensitive cavity ring-down spectrometer which operates in the mid-infrared spectral region near 4.5 μm was constructed in used to acquire a noise-equivalent absorption coefficient of 2.6×10−11 cm−1 Hz−1/2 with less than 150 nW of optical power incident on a photodetector. Quantum noise was observed in the individual ring-down decay events, leading to quantum-noise-limited short-time performance. A combination of high sensitivity and robustness make it well suited for measurements of ultra-trace gas species as well as applications in optics and fundamental physics.
Cavity ringdown spectroscopy (CRDS) is a cavity-enhanced technique in which a high finesse Fabry-Pérot cavity is used to achieve a long path length (10's of km to 100's of km). During a CRDS measurement the cavity is first optically pumped by the laser source. This light is then extinguished and the corresponding passive decay of the intracavity optical power is monitored. The exponential decay time, the ring-down time constant (τ), is directly related to the absorption, α, within the cavity:
wherein c is the speed of light and τ0 is the empty-cavity ring-down time constant which accounts for mirror losses due to transmission, absorption, and scattering. CRDS is insensitive to laser amplitude fluctuations.
While CRDS is routinely performed in the near-infrared, working in the mid-infrared spectral region can be far more challenging due to generally lower component performance and higher component costs. Importantly, supermirror coatings are more difficult to produce in the mid-infrared, leading to much higher absorption and scattering losses. These significant losses lead to far lower build-up of optical power within the cavity and correspondingly lower optical transmission. As a result, very low noise photodetectors are used to achieve high signal-to-noise ratios on the ring-down decay events. In addition, high efficiency electro-optic modulators are not available in the mid-infrared
A schematic of the cavity ring-down spectrometer is shown in
The optical resonator had a nominal length of 1.5 m and a corresponding free spectral range of 100 MHz. The ring-down mirrors had a radius of curvature of 1 m and a power reflectivity of 99.99% leading to a finesse of 31 000 and an effective path length of 15 km. We estimate that the cavity mode matching efficiency was >90%. The ring-down decays were recorded on a 0.1 mm diameter liquid-nitrogen-cooled InSb photodetector with a field of view of 60°, a responsivity of 3.7 A/W, a trans impedance gain of 106 V/A, and a measured noise-equivalent power of 70 fW Hz−1/2 and digitized at 1 MSamples s−1 by a 22-bit acquisition board. The 3-dB bandwidth of the acquisition board was measured to be B=480 kHz. Given its high bit depth, the digitizer board did not contribute significantly to the overall technical noise. The laser was stabilized to the optical resonator via a low bandwidth (4 Hz) transmission lock which actuated the laser current. Spectral scanning was performed by stepping the laser temperature in increments of the cavity's free spectral range.
The expected root-mean-square (RMS) noise on a given ring-down decay event, σV(t) is a combination of two-sided quantum and technical sources as follows
wherein e is the electron charge, G is the trans impedance gain, V(t) is the observed ring-down signal amplitude, and Δt is the sampling interval.
The technical noise is given by σtech=(NEP/√{square root over (2)})RG√{square root over (B)}, where NEP is the one-sided detector noise-equivalent power, R is the detector responsivity, and B is the electronic bandwidth. The quantum noise component scales as the square-root of the signal amplitude while the technical noise component is constant.
The noise within each decay event can be measured by applying the Wiener-Khintchine theorem, in which we equate the noise power spectral density σV,d2(f) (V2/Hz) to the Fourier transform of the autocorrelation function of the fit residual r(t)=V(t)−Vfit(t). Here f is the Fourier frequency and Vfit(t) is an exponential fit to the decay signal. Symbolically,
from which the integrated noise power (variance) can be found as:
σV2=∫0∞σV,d2(f)df (4)
Equations (3) and (4) were evaluated using discrete summations for the autocorrelation and Fourier transform. To measure the time dependence of the noise power spectral density, σV,2(f) was computed for averaging bins of width T<<τ and at times ti relative to the beginning of the fitting window. In the absence of correlation, Eqs. (3) and (4) yield the variance of the residual signal within an averaging bin.
The signal noise was also measured by analyzing the distribution of fit residuals for an ensemble of decays. For each time tj relative to the beginning of the fit window, the noise power was given by the variance in the ensemble of residuals. The time-dependent noise calculated in this manner was in good agreement with the single-shot analysis described above, indicating that the noise behaves ergodically. Both approaches yield quantum-noise-limited results which are consistent with Eq. (2), and which decay exponentially with a time constant of τ/2 until the technical noise floor is reached. The agreement between the autocorrelation of the fit residuals and Eq. (2) can be observed in
In the quantum noise limit (QNL) the relative standard uncertainty in the determination of τ from an individual fit to a decay event is given by the inverse of the square root of the number of photoelectrons:
While in the technical noise limit (TNL) the relative standard uncertainty in τ is given by:
The quadrature sum of these two expressions can then be used to approximate the relative uncertainty in the regime where both quantum noise and technical noise are considerable. From Eqs. (5) and (6), the expected quantum-noise-limited and technical-noise-limited (i.e. detector limited) single-decay-event fit uncertainties were 0.008% and 0.015%, respectively, where V(0)=0.50(06) V (corresponding to 136(16) nW). Adding these uncertainties in quadrature yields an approximate total expected fit uncertainty of 0.017%. Experimentally, average fit uncertainties were 0.019%. A typical 1 s ensemble of ring-down time constants (i.e. the first 30 time constants from
The presence of quantum noise during a significant fraction of the cavity decays necessitates the use of a weighted least-squares fit with weighting factors w(t)=1/σV(t)2, where σV(t) is calculated from Eq. (2). An estimate of the proper weighting factors for any given ensemble requires V(t), which is approximated by the average fit of the first 10 cavity decay events using equal weights. When fitting the portion of each cavity decay that exhibits significant quantum noise (t<100 μs), equal weighting of the data yields an average fit uncertainty that is 20% higher than that observed with proper weighting. Each cavity decay event has a relatively constant offset arising from dark current in the InSb photodetector. For the purpose of calculating w(t), this dark-current offset has been subtracted from V(t) to ensure that the quantum noise contribution to Eq. (2) is not overestimated.
An Allan variance plot, which is a measure of the system stability, can be found in
With reference to
The relatively low cost instrument described herein is well suited for precise measurements of ultra-trace gas species including radiocarbon and atmospherically relevant free radicals. This instrument has been used to measure birefringence in supermirror coatings and can be applied in addressing other challenges in fundamental physics including the search for symmetrization postulate violations in molecular physics and measurements of absolute number densities for magnetically-trapped ultra-cold molecules using absorption-based techniques.
A linear absorption spectrometer that included a cavity ring-down spectrometer with a quantum cascade laser was constructed for optical detection of radiocarbon in the mid-infrared spectral region. The linear absorption spectrometer had a distributed-feedback quantum cascade laser to provide single frequency mid-infrared radiation to probe the 14CO2 v3 fundamental transition near 4.5 μm. We achieved a noise-equivalent absorption coefficient of 2×10−11 cm−1 Hz−1/2.
The linear absorption spectrometer included a low temperature absorption cell with a closed cycle refrigeration system to cool the cell down to 208 K, reducing the interferences from other isotopologues of CO2. The instrument was fully automated and ran unattended for up to 8 h. We made a series of measurements of anthropogenic and biogenic carbon dioxide samples and can quantify the relative abundance of 14C within them. Data acquired with the linear absorption spectrometer shown in
In this Example, equations are numbered beginning with Equation 1.
High-finesse optical resonators found in ultrasensitive laser spectrometers utilize supermirrors ideally consisting of isotropic high-reflectivity coatings. Strictly speaking, however, the optical coatings are often nonuniformly stressed during the deposition process and therefore possess birefringence. When physically mounted, the cavity mirrors can be additionally stressed in such a way that large optical birefringence is induced. Here we report a direct measurement of optical birefringence in a two-mirror Fabry-Perot cavity with R=99.99% by observing TEM00 mode beating during cavity decays. Experiments were performed at a wavelength of 4.53 with precision limited by both quantum and technical noise sources. We report a splitting of δv=618(1) Hz, significantly less than the intrinsic cavity line width of δcav≈3 kHz. With a cavity free spectral range of 96.9 MHz, the equivalent fractional change in mirror refractive index due to birefringence is therefore Δn/n=6.38(1)×10−6.
By performing CRDS with a high signal-to-noise ratio in the regime where cavity losses were greater than the round-trip phase retardance, we directly measured mirror birefringence via determination of the lowest-order transverse mode (TEM00) splitting. A high-finesse optical resonator (F≈31000) constructed using mid-infrared mirrors with observable birefringence resulted in a unique set of net slow and fast axes which supported orthogonal linear polarization states with different resonant optical frequencies. We observed this small difference in optical frequencies as a beating during the cavity decays for specific input photon polarization states in the presence of appropriate polarization analysis. This method of determining the supermirror birefringence by a measurement of the beat frequency from TEM00 mode splitting was, in a sense, a form of birefringence interferometry [4]. Our high-precision CRDS approach had high sensitivity for small amounts of birefringence (10−8) as well as more than three orders of magnitude of dynamic range. We also present a complementary generalized model for identifying an ideal input polarization state to optimize high precision retrieval of the cavity decay time constant (and therefore cavity losses) even in the presence of significant mirror birefringence.
Our model of supermirror birefringence in a high-finesse optical resonator begins with the projection of an electric field of arbitrary polarization onto an arbitrary set of spatially orthogonal axes 1 and 2:
E(t)=E1(t)exp[i(ω1t−φ1)]+E2(t)exp[i(ω2tφ2)]. (1)
Here, E1(t) and E2(t) are time-dependent electric field magnitudes, ω1 and ω2 are the optical angular frequencies of each field (ω1=2πv1 and ω2=2πv2, respectively), and φ1 and φ2 are the phases of each optical field. The natural choice for a basis set is the net slow and fast axes of our optical resonator under study, a linear combination of the slow and fast optical axes of each individual mirror with observable birefringence. We aim to measure the cavity decay after projection onto a linear polarization analyzer (PA) orientated at an angle γ relative to the cavity slow axis. The intensity of the decaying optical field incident on a photodetector after the PA is
I(γ,t)=E(γ,t)E*(γ,t)=a1f1I1(t)cos2(γ)+a2f2I2(t)sin2(γ)+√{square root over (a1f1I1(t)a2f2I2(t))}sin(2γ)cos(ωbt−δϕ), (2)
where a1 and a2 are the fraction of light intensity propagating along the slow or fast axes (where a1+a2=1), f1 and f2 are frequency-dependent mode coupling factors (where f1+f2=1), Im(t)=I0 exp(−t/τm) for m=1 and 2, and ωb and δφ are the difference in angular frequency and phase between the slow and fast electric fields, respectively. In the limit where the laser line width is much greater than both the cavity line width and the birefringence splitting (i.e., Δvlaser δcav and >>ωb/(2π)) on the timescale of cavity buildup, the frequency-dependent mode coupling factors become identical and constant.
If in addition we consider the mirror reflectivity to be the same for both axes (i.e., τ1=τ2≡τ), Eq. (2) simplifies to Eq. (3):
The factors A(γ) and B(γ) depend upon the polarization state by which the optical resonator is excited. Here we investigate in detail two specific cases: excitation with circularly polarized light and excitation with linearly polarized light. For circular excitation, a1=a2=1/2 and the absolute difference in phase is |δφ|=π/2. Thus, A(γ)=1, B(γ)=sin(2γ), and δφ=±π/2 for left-hand (LH) and right-hand (RH) excitation, respectively. Equation (3) can then be further simplified using the small-angle approximation sin(ωbt)≈ωbt:
In the limit of weak birefringence and using circular excitation, the cavity decay filtered by a linear PA will be nearly exponential with an effective time constant of τeff(γ)=τ/[1∓B(γ)ωbT] for LH and RH excitation. The maximum fractional deviation from τ is then defined as eff,max
τ≈ωbτ. (5)
For excitation with linearly polarized light, the following conditions are met: η=tan−1(√a2/a1) where η is the angle the incident linearly polarized field makes with the net slow axis, and δφ=0. In each of the experiments described here using linearly polarized light, η=γ. Making the substitution, A(γ)=cos4(γ)+sin4(γ) and B(γ)=sin2(2γ)/2, a similar expression to Eq. (4) can then be derived using the small-angle approximation cos(ωbt)≈1−(ωbt)2/2:
In SI units to I0=pc∈0Tc(<I1,0>+<I2,0>)/2, where the coefficient p=1/4 or p=1/2 for circular or linear excitation, respectively, is the permittivity of free space, Tc is the overall efficiency with which the incident spectrum is transmitted by the optical resonator, and denotes an averaging over optical cycles. Since our cavity design does not allow for one mirror to be rotated in a systematic fashion relative to the other, we must describe the two-mirror optical resonator and thus only model net cavity birefringence.
The above expressions provide an intuitive, parameterized description of a single-mode cavity decay with TEM00 mode beating that arises from the presence of mirror birefringence. As evidenced by the experiments reported here, the general time-domain approach quantitatively described cavity excitation by both circularly and linearly polarized light without ignoring polarization-dependent loss. We also introduce frequency-dependent mode coupling factors that, for the case of a swept laser-cavity locking scheme, renders these analytical expressions useful for cases of both weak and strong birefringence (relative to the round-trip cavity losses).
A cavity round-trip Jones matrix M contains the following parameters: the round-trip phase retardance ∈=ωbtr/2 and the angle α the slow axis makes with the x axis in the laboratory frame. We can now make a connection between the perturbed exponential expressions derived above from an intuitive description of transverse mode beating and the Jones matrix approach. Through Eq. (5) Δτ/τ=∈/(1−R and the above definition of round-trip phase retardance, we find for excitation of the cavity by circular polarization that) in the limit of weak birefringence. In its general form, M also treats polarization-dependent loss in each cavity mirror, a small perturbation on supermirror birefringence that has not been included in Eqs. (3)-(6).
A schematic of the linear absorption spectrometer that performed ultrasensitive cavity ring-down spectroscopy (CRDS) in the mid-infrared spectral region is shown in
The high-finesse optical resonator under interrogation included two supermirrors from the same batch coated for maximum power reflectivity R at a wavelength of 4.6 μm (CRD Optics 901-0008-4600). The supermirrors were created using ZnSe substrates of the following dimensions: a thickness of 5 mm, a diameter of 2 cm, and a radius of curvature of 1 m. The mirrors were glued to 2.54 cm outer diameter (OD) adapter rings using standard two-part epoxy and then mounted inside a threaded stainless steel knife-edge flange using a single aluminum retaining ring. Each of the two mirror mounts was adjusted using high-thread-count set screws to aid in initial alignment of the optical resonator. The optical resonator length was L=1.55 m, and the round-trip absorption and scattering losses were measured to be 1×10−4.
The individual zero-pressure difference mirror mounts were supported by stainless steel brackets, each of which was connected to the other by four spacer members (2.54 cm OD Invar-36 rods). The supermirror zero-pressure difference mirror mounts were connected by a stainless steel tube of 1.27 cm OD and a stainless steel bellows to reduce strain in the event of slow drifts in cavity length due to room temperature fluctuations. This tube housed the optical mode as well as any sample gas under study. Each mirror mount was capped with a CaF2 window with antireflective coatings which sealed the entire optical resonator from the laboratory environment. The mirrors in this design experience a zero-pressure difference between their antireflective and reflective faces, significantly reducing the potential for stress-induced birefringence. Experiments were performed under the vacuum of a turbomolecular pump at pressures <1.3 Pa.
Results of exciting the two-mirror optical resonator with circularly polarized light are presented. A quarter-wave plate (Altechna 2-IRPW-L/4-4500-C) and a high-contrast PA (ISPOpticsPOL-1.5-5-SI,extinction 10000:1) were arranged at positions WP and PA as shown in
At each angle θ the average empty-cavity decay constant τeff for a 3 s acquisition ensemble is shown in
For each respective handedness, the LH and RH data sets were fit to the model τeff(θ)=τ[(1∓ωbτsin [2(θ−α)])]−1, where ωb=2πδv. The fit parameters from
Of the fit τeff in
The observed cavity decays when στ is at its extrema is presented here.
At each value of θ, the NEM and JMM provided fitted values for the amplitude of the decay, a single global time constant τ and the beat frequency E. Panel E of
Excitation of the cavity by linear polarization can also be treated by both the NEM and JMM using the appropriate coefficients A(γ) and B(γ) in Eq. (3). When η=γ, rotating η and γ together is equivalent to rotating the cavity itself relative to a fixed linear polarization. We note that, at a constant trigger threshold, this approach does not suffer from any potential bias due to changes in the intracavity power as a function of the angle of the linear polarization analyzer. Exciting the cavity at θ=33° (near where τ/τ for circular polarization in
The Allan deviation for the linear excitation measurement shown in panel F of
In CRDS, we desire to simply avoid mirror birefringence altogether to achieve optimized fit statistics for the unperturbed cavity time constant τ. With knowledge of the spatial location of the slow and fast cavity axes in the laboratory frame we can selectively excite one of the modes and measure the Allan deviation to determine an NEA normalized to a 1 s acquisition. When θ=168° (η=γ≈3°, I0=203 nW) we report στ/τ=0.024% for the SEM, within a factor of two of the theoretical value of στ/τ=0.012% calculated from the quadrature sum of quantum and technical noise when the fitting is performed with appropriate weighting. This results in an NEA=5.5×10−11 cm−1 Hz−1/2 at an acquisition rate of only 9 Hz. If in the laboratory the exact condition θ=α (or θ=α±90°) is met with excitation by linear polarization, we can indeed achieve the theoretical limit of στ/τ set by both quantum and detector noise.
When (and therefore δv<δcav) both the NEM and JMM perform equally well in the limit where PDL is negligible (and at the experimental signal-to-noise ratio on the cavity decays of SNR 3000:1). To measure the PDL we have fit the average cavity decays at each value of θ for the LH and RH data sets to a global set of parameters that includes PDL in the JMM. Each averaged cavity decay was normalized by its previously fit amplitude and then collectively fit using the JMM+PDL with the following floated terms: τ, α, ∈, b, and β, where b=(rmax−rmin)/(rmax+rmin) (r is the net field reflectivity of the mirrors) and β is the angle rmax makes with the x axis. We report b/∈=0.057(16) with β=17(12)° and b/∈=0.09(7) with β=60(20)° for the LH and RH data sets, respectively. Clearly PDL is a small perturbation on the mirror birefringence, and we can safely place an upper bound of b/∈<0.1 on the net cavity PDL for this pair of mirrors.
The simultaneous fit over all θ within a data set returns large uncertainties in β, again suggesting that PDL is small. A TEM00 mode splitting approaches the cavity line width δv≈δcav), involving frequency−∈≈1−R dependent mode coupling factors f1 and f2 in the NEM when a ∈≈1−R narrow line width continuous-wave laser is coupled to the optical resonator. In this strong birefringence regime (δv≥δcav), f1 and f2 will change for every cavity decay recorded by our spectrometer due to the swept laser locking and thresholding used to trigger the optical switch. The small-angle approximations used to derive Eqs. (4)-(6) are not justified when, and the SEM fails even more dramatically than shown in
By replacing M1 of the optical resonator with a presumably identical supermirror we fortuitously observed strong net birefringence in the regime. For LH, RH, and linear excitation we identified several PA orientations in the laboratory frame where the SEM produced large fitted residuals. For 133 total cavity decays fitted by Eq. (2) (still under the assumption of no PDL, i.e., τ1=τ2≡τ) at various θ we report δv=3.1 kHz with a standard deviation of 0.4 kHz. We observed a maximum deviation from the SEM for LH excitation at θ=138° (i.e., α=93°) and proceeded to record 170 successive cavity decays at an acquisition rate of 17 Hz with the PA in that orientation. The average beat frequency for the θ=138° data set was δv=2.789 kHz with 9 Hz precision at 5 s of averaging. An Allan plot revealed that the √ deviation of δv remained inversely proportional to the N, where N is the number of cavity decays, for a minimum of 5 s. Therefore, with longer averaging times it may be possible to again achieve 1 Hz precision even in the presence of strong mirror birefringence.
The origin of the observed strong mirror birefringence that coincided with replacing a single mirror of the optical resonator could be reasoned in one of two ways. Since we measured net birefringence, the new M1 could simply have significantly larger birefringence than the first M1. Alternatively, the relative orientations of αM1 and αM2 could be very different in the two optical resonators, meaning that in our first optical resonator with δv=618 Hz the slow and fast axes of the individual supermirrors were oriented close to 90° from one another in the laboratory frame.
Mid-infrared supermirror birefringence using cavity ring-down spectroscopy was determined. Analysis of the beating observed during cavity decays yielded a Δn/n sensitivity on the order of 10−8 at a wavelength of 4.5 μm.
In an optical cavity, after every round-trip, photons leak from the cavity by a fractional amount equal to the sum of the mirror transmission T and losses L (the total of mirror absorption and scattering as well as intracavity losses). Remaining photons are reflected at the mirror surface. Therefore, the fate of any photon following one resonator round trip can be described by the expression R+T+L=1, where R is the power reflection coefficient. The rate at which photons leak from an optical resonator is τ=l/[c(1−R)], where l is the resonator length. By precisely measuring τ, detailed information regarding the total intracavity losses T+L is learned. In cavity ring-down spectroscopy (CRDS), the precision measurement of τ as a function frequency allows for the molecular spectroscopy of trace gases to be performed with ultrahigh sensitivity.
Determined here were perturbations to empty cavity τ, τ0, that arise from ultralow supermirror birefringence. We excited the mid-infrared optical resonator (to be used for CRDS) with either a circularly or a linearly polarized continuous-wave (CW) laser and thus identified the spatial location of each birefringence cavity eigenmode. This resulted in a measurement of the relative difference in their refractive index, Δn/n, with 10−8 precision. Identifying the polarization conditions that resulted in the highest precision retrievals of τ provided optical detection of radiocarbon dioxide, 14CO2, below its nature abundance of 1.2 parts-per-trillion (ppt).
Panel A of
Panel C of
By including mode beating in or model of the CRDS signals, we extracted a high-precision measurement of δv. Panel A of
The data represent direct measurement of supermirror birefringence by birefringence interferometry, i.e., the time-domain measurement of optical mode beating in the CRDS signals.
Cavity-enhanced mid-infrared spectroscopy provided by the linear absorption spectrometer provided trace detection of rare CO2 isotopologues, specifically, radiocarbon dioxide (14CO2) at a wavelength of 4.5 μm. Radiocarbon dioxide can be used to clearly differentiation) between CO2 that originated from either biogenic or anthropogenic sources. Here, we report the all optical, quantitative detection of 14CO2 using a QCL and CRDS as shown in
For the experiments in this Example, a linearly polarized CW quantum cascade laser (QCL) operated at a wavelength of 4.5 μm to excite one longitudinal TEM00 mode of a high finesse optical resonator including a two-mirror Fabry-Perot optical resonator. Wave plates, either λ/4 or λ/2, were used to control the laser polarization prior to the optical resonator, and a linear polarization analyzer was used after the cavity but before the photodetector to identify the spatial location of each linear birefringence eigenmode. An acousto-optic modulator was used as a fast optical switch to initiate CRDS signals. Laser-cavity coupling was maintained by a loose, low-bandwidth (4 Hz) transmission lock.
While one or more embodiments have been shown and described, modifications and substitutions may be made thereto without departing from the spirit and scope of the invention. Accordingly, it is to be understood that the present invention has been described by way of illustrations and not limitation. Embodiments herein can be used independently or can be combined.
Reference throughout this specification to “one embodiment,” “particular embodiment,” “certain embodiment,” “an embodiment,” or the like means that a particular feature, structure, or characteristic described in connection with the embodiment is included in at least one embodiment. Thus, appearances of these phrases (e.g., “in one embodiment” or “in an embodiment”) throughout this specification are not necessarily all referring to the same embodiment, but may. Furthermore, particular features, structures, or characteristics may be combined in any suitable manner, as would be apparent to one of ordinary skill in the art from this disclosure, in one or more embodiments.
All ranges disclosed herein are inclusive of the endpoints, and the endpoints are independently combinable with each other. The ranges are continuous and thus contain every value and subset thereof in the range. Unless otherwise stated or contextually inapplicable, all percentages, when expressing a quantity, are weight percentages. The suffix “(s)” as used herein is intended to include both the singular and the plural of the term that it modifies, thereby including at least one of that term (e.g., the colorant(s) includes at least one colorants). “Optional” or “optionally” means that the subsequently described event or circumstance can or cannot occur, and that the description includes instances where the event occurs and instances where it does not. As used herein, “combination” is inclusive of blends, mixtures, alloys, reaction products, and the like.
As used herein, “a combination thereof” refers to a combination comprising at least one of the named constituents, components, compounds, or elements, optionally together with one or more of the same class of constituents, components, compounds, or elements.
All references are incorporated herein by reference.
The use of the terms “a” and “an” and “the” and similar referents in the context of describing the invention (especially in the context of the following claims) are to be construed to cover both the singular and the plural, unless otherwise indicated herein or clearly contradicted by context. “Or” means “and/or.” Further, the conjunction “or” is used to link objects of a list or alternatives and is not disjunctive; rather the elements can be used separately or can be combined together under appropriate circumstances. It should further be noted that the terms “first,” “second,” “primary,” “secondary,” and the like herein do not denote any order, quantity, or importance, but rather are used to distinguish one element from another. The modifier “about” used in connection with a quantity is inclusive of the stated value and has the meaning dictated by the context (e.g., it includes the degree of error associated with measurement of the particular quantity).
This application claims the benefit of U.S. Provisional Patent Application Ser. Nos. 62/430,129, filed Dec. 5, 2016, the disclosure of which is incorporated herein by reference in its entirety.
This invention was made with United States Government support from the National Institute of Standards and Technology, an agency of the United States Department of Commerce. The Government has certain rights in the invention.
Number | Date | Country | |
---|---|---|---|
62430129 | Dec 2016 | US |