Laser absorption spectroscopy isotopic gas analyzer

Abstract

The present invention provides systems and methods for measuring the isotope ratios of one or more trace gases based on optical absorption spectroscopy methods. The system includes an optical cavity containing a gas. The system also includes a laser optically coupled with the optical cavity, and a detector system for measuring absorption of laser light by the gas in the cavity.

Claims

1. A gas analyzer system for measuring an isotopic ratio of a gaseous chemical species by an optical absorption spectroscopy method, the system comprises: an optical cavity containing a gas with the chemical species to be measured; a laser optically coupled to an optical cavity; a detector system for measuring absorption of laser light by the gas; and an intelligence module comprising a processor adapted to determine a concentration ratio of two different isotopologues, wherein rotational-vibrational spectra of the chemical species are measured at least within two non-overlapping spectral intervals separated by a minimum of 50 nm between any two spectral points from each of the spectral intervals and selected in such a way that a first rotational-vibrational line of a first rotational-vibrational band of a less abundant isotopologue is located in a first spectral interval, and a second rotational-vibrational line of a second rotational-vibrational band of a more abundant isotopologue is located in a second spectral interval, and a line intensity of the strongest line of the first rotational-vibrational band of the less abundant isotopologue is two or more times stronger than a line intensity of the strongest line of the second rotational-vibrational band of the more abundant isotopologue, and a ratio of an absorption spectrum of the less abundant isotopologue to the sum of absorption spectra of all other isotopologues of the chemical species weighted by mole-fraction abundance figures exceeds two somewhere in the first spectral interval, and a ratio of an absorption spectrum of the more abundant isotopologue to the sum of absorption spectra of all other isotopologues of the chemical species weighted by mole-fraction abundance figures exceeds two somewhere in the second spectral interval.

2. The system of claim 1, wherein the pressure broadening coefficients of the first line and the second line are different by no more than 50%.

3. The system of claim 1, wherein the temperature dependence coefficients of the first line and the second line are different by no more than 50%.

4. The system of claim 1, wherein the optical absorption spectroscopy method comprising of at least one of the following methods: the cavity ring down spectroscopy method, the cavity enhanced absorption spectroscopy method, the cavity phase shift spectroscopy method, the integrated cavity output spectroscopy method, or the cavity enhanced photo-acoustic spectroscopy method.

5. The system of claim 1, wherein the optical cavity is a resonant optical cavity.

6. The system of claim 1, wherein the detector system includes one of a photo-detector configured to measure an intensity of the intra-cavity light or both a photo-acoustic sensor configured to measure photo-acoustic waves generated in the cavity and a photo-detector configured to measure an intensity of the intra-cavity light.

7. The system of claim 1, wherein one of the measured isotopologues is .sup.16O.sup.13C.sup.16O and another measured isotopologue is .sup.16O.sup.12C.sup.16O.

8. The system of claim 1, further comprising a temperature sensor for measuring a temperature of the gas in the optical cavity; and a pressure sensor for measuring a pressure of the gas in the optical cavity.

9. The system of claim 8, further comprising a temperature control element configured to control a temperature of the gas in the optical cavity and a pressure control element configured to control a pressure of the gas in the optical cavity.

10. A gas analyzer system for measuring an isotopic ratio of a gaseous chemical species by an optical absorption spectroscopy method, the system comprises: an optical cell containing a gas with the chemical species to be measured; a laser configured to emit light into an optical cell; a detector system for measuring absorption of laser light by the gas; and an intelligence module comprising a processor adapted to determine a concentration ratio of two different isotopologues, wherein rotational-vibrational spectra of the chemical species are measured at least within two non-overlapping spectral intervals separated by a minimum of 50 nm between any two spectral points from each of the spectral intervals and selected in such a way that a first rotational-vibrational line of a first rotational-vibrational band of a less abundant isotopologue is located in a first spectral interval, and a second rotational-vibrational lines of a second rotational-vibrational band of a more abundant isotopologue is located in a second spectral interval, and a line intensity of the strongest line of the first rotational-vibrational band of the less abundant isotopologue is two or more times stronger than a line intensity of the strongest line of the second rotational-vibrational band of the more abundant isotopologue, and a ratio of an absorption spectrum of the less abundant isotopologue to the sum of absorption spectra of all other isotopologues of the chemical species weighted by mole-fraction abundance figures exceeds two somewhere in the first spectral interval, and a ratio of an absorption spectrum of the more abundant isotopologue to the sum of absorption spectra of all other isotopologues of the chemical species weighted by mole-fraction abundance figures exceeds two somewhere in the second spectral interval.

11. A method of measuring an isotopic ratio of a gaseous chemical species, the method comprising: coupling laser light to an optical cavity containing a gas with the chemical species to be measured; measuring an absorption of the laser light by the gas; determining the concentration ratio of two different isotopologues, wherein rotational-vibrational spectra of the chemical species are measured at least within two non-overlapping spectral intervals separated by a minimum of 50 nm between any two spectral points from each of the spectral intervals and selected in such way that a first rotational-vibrational line of a first rotational-vibrational band of a less abundant isotopologue is located in a first spectral interval, and a second rotational-vibrational line of a second rotational-vibrational band of a more abundant isotopologue is located in a second spectral interval, and a line intensity of the strongest line of the first rotational-vibrational band of the less abundant isotopologue is two or more times stronger than a line intensity of the strongest line of the second rotational-vibrational band of the more abundant isotopologue, and a ratio of an absorption spectrum of the less abundant isotopologue to the sum of absorption spectra of all other isotopologues of the chemical species weighted by mole-fraction abundance figures exceeds two somewhere in the first spectral interval, and a ratio of an absorption spectrum of the more abundant isotopologue to the sum of absorption spectra of all other isotopoloques of the chemical species weighted by mole-fraction abundance figures exceeds two somewhere in the second spectral interval.

12. The method of claim 11, wherein the pressure broadening coefficients of the first line and the second line are different by no more than 50%.

13. The method of claim 11, wherein the temperature dependence coefficients of the first line and the second line are different by no more than 50%.

14. The method of claim 11, wherein the optical absorption spectroscopy method comprising at least one of the following methods: the cavity ring down spectroscopy method, the cavity enhanced absorption spectroscopy method, the cavity phase shift spectroscopy method, the integrated cavity output spectroscopy method, or the cavity enhanced photo-acoustic spectroscopy method.

15. The method of claim 11, wherein the optical cavity is a resonant optical cavity.

16. The method of claim 11, wherein one of the measured isotopologues is .sup.13C.sup.16O.sub.2 and another measured isotopologue is .sup.12C.sup.16O.sub.2.

17. The method of claim 11, further comprising: measuring a temperature of the gas in the optical cavity by a temperature sensor and measuring a pressure of the gas in the optical cavity by a pressure sensor.

18. The method of claim 17, further comprising: controlling a temperature of the gas in the optical cavity by a temperature control element and controlling a pressure of the gas in the optical cavity by a pressure control element.

Description

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWING

(1) FIG. 1. Line intensities of different rotational-vibrational transitions from one rotational level in the ground vibrational state to one rotational level in the vibrationally excited state v.sub.3 of .sup.12C.sup.16O.sub.2 isotopologue at 296 K.

(2) FIG. 2. Line intensities of different rotational-vibrational transitions from one rotational level in the ground vibrational state to one rotational level in the vibrationally excited state v.sub.1+v.sub.3 of .sup.12C.sup.16O.sub.2 isotopologue at 296 K.

(3) FIG. 3. Line intensities of different rotational-vibrational transitions from one rotational level in the ground vibrational state to one rotational level in the vibrationally excited state v.sub.1+2v.sub.2+v.sub.3 of .sup.12C.sup.16O.sub.2 isotopologue at 296 K.

(4) FIG. 4. Line intensities of different rotational-vibrational transitions from one rotational level in the ground vibrational state to one rotational level in the vibrationally excited state v.sub.1+4v.sub.2+v.sub.3 of .sup.12C.sup.16O.sub.2 isotopologue at 296 K.

(5) FIG. 5. Line intensities of different rotational-vibrational transitions from one rotational level in the ground vibrational state to one rotational level in the vibrationally excited state 2v.sub.1+2v.sub.2+v.sub.3 of .sup.12C.sup.16O.sub.2 isotopologue at 296 K.

(6) FIG. 6. Line intensities of different rotational-vibrational transitions from one rotational level in the ground vibrational state to one rotational level in the vibrationally excited state v.sub.3 of .sup.13C.sup.16O.sub.2 isotopologue at 296 K.

(7) FIG. 7. Line intensities of different rotational-vibrational transitions from one rotational level in the ground vibrational state to one rotational level in the vibrationally excited state v.sub.1+v.sub.3 of .sup.13C.sup.16O.sub.2 isotopologue at 296 K.

(8) FIG. 8. Line intensities of different rotational-vibrational transitions from one rotational level in the ground vibrational state to one rotational level in the vibrationally excited state v.sub.1+2v.sub.2+v.sub.3 of .sup.13C.sup.16O.sub.2 isotopologue at 296 K.

(9) FIG. 9. Line intensities of different rotational-vibrational transitions from one rotational level in the ground vibrational state to one rotational level in the vibrationally excited state v.sub.1+4v.sub.2+v.sub.3 of .sup.13C.sup.16O.sub.2 isotopologue at 296 K.

(10) FIG. 10. Line intensities of different rotational-vibrational transitions from one rotational level in the ground vibrational state to one rotational level in the vibrationally excited state 2v.sub.1+2v.sub.2+v.sub.3 of .sup.13C.sup.16O.sub.2 isotopologue at 296 K.

(11) FIG. 11. The ratio of the spectrum of the .sup.13C.sup.16O.sub.2 isotopologue to the sum of spectra of all other isopologues of CO.sub.2 from 2040.5 nm to 2043 nm at a pressure of 100 mBar and a temperature of 296 K. R.sub.x symbols show positions where the corresponding rotational-vibrational transitions, from one rotational level in the ground vibrational state to another rotational level in the vibrationally excited state v.sub.1+2v.sub.2+v.sub.3 of .sup.13C.sup.16O.sub.2 isotopologue, are located.

(12) FIG. 12. The ratio of the spectrum of the .sup.13C.sup.16O.sub.2 isotopologue to the sum of spectra of all isopologues of water from 2040.5 nm to 2043 nm at a pressure of 100 mBar and a temperature of 296 K. R.sub.N symbols show positions where the corresponding rotational-vibrational transitions, from one rotational level in the ground vibrational state to another rotational level in the vibrationally excited state v.sub.1+2v.sub.2+v.sub.3 of .sup.13C.sup.16O.sub.2 isotopologue, are located.

(13) FIG. 13. The ratio of the spectrum of the .sup.13C.sup.16O.sub.2 isotopologue to the sum of spectra of all isopologues of methane from 2040.5 nm to 2043 nm at a pressure of 100 mBar and a temperature of 296 K. R.sub.N symbols show positions where the corresponding rotational-vibrational transitions, from one rotational level in the ground vibrational state to another rotational level in the vibrationally excited state v.sub.1+2v.sub.2+v.sub.3 of .sup.13C.sup.16O.sub.2 isotopologue, are located.

(14) FIG. 14. The ratio of the spectrum of the .sup.12C.sup.16O.sub.2 isotopologue to the sum of spectra of all other isopologues of CO.sub.2 from 1572.5 nm to 1575 nm at a pressure of 100 mBar and a temperature of 296 K. R.sub.x symbols show positions where the corresponding rotational-vibrational transitions, from one rotational level in the ground vibrational state to another rotational level in the vibrationally excited state 2v.sub.1+2v.sub.2+v.sub.3 of .sup.13C.sup.16O.sub.2 isotopologue, are located.

(15) FIG. 15. The ratio of the spectrum of the .sup.12C.sup.16O.sub.2 isotopologue to the sum of spectra of all isopologues of water from 1572.5 nm to 1575 nm at a pressure of 100 mBar and a temperature of 296 K. R.sub.x symbols show positions where the corresponding rotational-vibrational transitions, from one rotational level in the ground vibrational state to another rotational level in the vibrationally excited state 2v.sub.1+2v.sub.2+v.sub.3 of .sup.13C.sup.16O.sub.2 isotopologue, are located.

(16) FIG. 16. The ratio of the spectrum of the .sup.12C.sup.16O.sub.2 isotopologue to the sum of spectra of all isopologues of methane from 1572.5 nm to 1575 nm at a pressure of 100 mBar and a temperature of 296 K. R.sub.x symbols show positions where the corresponding rotational-vibrational transitions, from one rotational level in the ground vibrational state to another rotational level in the vibrationally excited state 2v.sub.1+2v.sub.2+v.sub.3 of .sup.13C.sup.16O.sub.2 isotopologue, are located.

(17) FIG. 17. Spectrum of the .sup.13C.sup.16O.sub.2 isotopologue, the sum of spectra of all other isopologues of CO.sub.2, and the sum of spectra of all isopologues of water from 2040.5 nm to 2043 nm at a pressure of 100 mBar and a temperature of 296 K. R.sub.x symbols show positions where the corresponding rotational-vibrational transitions, from one rotational level in the ground vibrational state to another rotational level in the vibrationally excited state v.sub.1+2v.sub.2+v.sub.3 of .sup.13C.sup.16O.sub.2 isotopologue, are located. The total concentration of carbon dioxide is 400 part per million. The total concentration of water is 1%.

(18) FIG. 18. Spectrum of the .sup.12C.sup.16O.sub.2 isotopologue, the sum of spectra of all other isopologues of CO.sub.2, the sum of spectra of all isopologues of water, and the sum of spectra of all isopologues of methane from 1572.5 nm to 1575 nm at a pressure of 100 mBar and a temperature of 296 K. R.sub.x symbols show positions where corresponding rotational-vibrational transitions from one rotational level in the ground vibrational state to another rotational level in the vibrationally excited state 2v.sub.1+2v.sub.2+v.sub.3 of .sup.13C.sup.16O.sub.2 isotopologue are located. The total concentration of carbon dioxide is 400 part per million. The total concentration of water is 1%.

DETAILED DESCRIPTION OF THE INVENTION

(19) Here and further we use the term the spectral line intensity in units of [cm.sup.1/(molecule cm.sup.2)] similar to what was given in FIG. 1 in the Appendix to the article on the 1996 Edition of HITRAN in the Journal of Quantitative Spectroscopy and Radiative Transfer vol. 60, pp. 665-710 (1998). However, in our case the spectral line intensities are specified per an isotopologue molecule, while the spectral line intensities that appear in HITRAN are weighted according to the natural terrestrial isotopic abundances.

(20) Systems and methods described herein may include or employ one or more spectrometers measuring rotational-vibrational spectra of different isotolologues in the gas phase. The rotational-vibrational spectra are often resolved into lines due to transitions from one rotational level in the ground vibrational state to one rotational level in the vibrationally excited state. The lines corresponding to a given vibrational transition form a band. The gas analyzer system measures rotational-vibrational spectra of isotopologues in the gas phase at least at two non-overlapping spectral regions: in a first spectral region the system measures at least one rotational-vibrational line of a less abundant isotopologue, and in a second spectral region the system measures at least one rotational-vibrational line of a more abundant isotopologue. Other lines of other isotopologues and other chemical species can also be measured in both spectral regions. To improve the system performance, these two lines belong to different vibrational modes chosen in such a way that the spectral line intensity of the strongest rotational-vibrational line of the vibrational mode of the less abundant isotopologue is at least two or more times stronger than the spectral line intensity of the strongest rotational-vibrational line of the vibrational mode of the more abundant isotopologue. This approach permits to improve both the precision and the accuracy of the isotope ratio measurements because the line intensities of two measured lines might be closer to each other, in comparison with the case when two lines of the same vibrational modes of two isotopologues are measured. Notice, in general, for spectral analysis one can choose rotational-vibrational lines with quite different spectral line intensities from the same vibrational modes of different isotopologues. However, in that case these lines usually have quite different pressure broadening and temperature dependences. So, the precise isotopic ratio analysis would require extremely accurate temperature and pressure stabilizations of the tested gas.

(21) FIGS. 1-10 show examples of the line intensities of different vibrational modes of the two most abundant isotopologues of carbon dioxide in Mid-IR and NIR-IR spectral regions. Circles represent different rotational-vibrational transitions from one rotational level in the ground vibrational state to another rotational level in the vibrationally excited state. The line intensities were obtained from the HITRAN database. The energy change of rotation can be either subtracted from or added to the energy change of vibration, giving the P- and R-branches of the spectrum, respectively. Both P- and R branches are shown. Taking to account that abundance of carbon dioxide isopologues, according to HITAN is 98.42% and 1.106% for .sup.12C.sup.16O.sub.2 and .sup.13C.sup.16O.sub.2, respectively, the following pairs of rotational-vibrational bands can be considered for spectral analysis according to one of the embodiments: 1) v.sub.1+v.sub.3 vibrational mode of .sup.12C.sup.16O.sub.2 and v.sub.3 vibrational mode of .sup.13C.sup.16O.sub.2; 2) v.sub.1+2v.sub.2+v.sub.3 vibrational mode of .sup.12C.sup.16O.sub.2 and v.sub.3 vibrational mode of .sup.13C.sup.16O.sub.2; 3) v.sub.1+4v.sub.2+v.sub.3 vibrational mode of .sup.12C.sup.16O.sub.2 and v.sub.3 vibrational mode of .sup.13C.sup.16O.sub.2; 4) 2v.sub.1+2v.sub.2+v.sub.3 vibrational mode of .sup.12C.sup.16O.sub.2 and v.sub.3 vibrational mode of .sup.13C.sup.16O.sub.2; 5) v.sub.1+2v.sub.2+v.sub.3 vibrational mode of .sup.12C.sup.16O.sub.2 and v.sub.1+v.sub.3 vibrational mode of .sup.13C.sup.16O.sub.2; 6) v.sub.1+4v.sub.2+v.sub.3 vibrational mode of .sup.12C.sup.16O.sub.2 and v.sub.1+v.sub.3 vibrational mode of .sup.13C.sup.16O.sub.2; 7) 2v.sub.1+2v.sub.2+v.sub.3 mode of for .sup.12C.sup.16O.sub.2 and v.sub.1+v.sub.3 mode of .sup.13C.sup.16O.sub.2; 8) v.sub.1+4v.sub.2+v.sub.3 vibrational mode of .sup.12C.sup.16O.sub.2 and v.sub.1+2v.sub.2+v.sub.3 vibrational mode of .sup.13C.sup.16O.sub.2; 9) 2v.sub.1+2v.sub.2+v.sub.3 vibrational mode of .sup.12C.sup.16O.sub.2 and v.sub.1+2v.sub.2+v.sub.3 vibrational mode of .sup.13C.sup.16O.sub.2;
where v.sub.1, v.sub.2, and v.sub.3 represent normal modes of the CO.sub.2 molecule: symmetric stretch, bend, and asymmetric stretch, respectively.

(22) By comparing FIGS. 1-10 one can easily see that the bands in these pairs are spectrally well separated and the line intensity of the strongest line in the band of the less abundant isotopologue is at least as twice as strong as the line intensity of the strongest line in the band of the more abandon isotopologue.

(23) Line intensity is the integrated absorption cross section across a line. Because the line intensities are temperature dependent, the line intensities of different lines should be compared at the temperature of the measured gas. The temperature dependence coefficient is defined as a derivative of the line intensity over temperature divided by the line intensity itself. The temperature dependence coefficients depend on temperature and according to one embodiment they are compared at the temperature of a measured gas mixture.

(24) After a pair of rotational-vibrational bands corresponding to two different isotopologues was selected for spectral analysis, spectral ranges for measuring the isotopologue concentration ratios can be chosen according to another embodiment: at least one of the rotational-vibrational lines of the band of the less abundant isotopologue is located and measured in the first spectral region, and at least one of the rotational-vibrational lines of the band of the more abundant isotopologue is located and measured in the second spectral region.

(25) According to another embodiment the ratio of an absorption spectrum of the less abundant isotopologue in the first selected spectral region to the sum of the absorption spectra of all other isotopologues of the same chemical substance exceeds two somewhere in the first selected spectral region, and the ratio of an absorption spectrum of the more abundant isotopologue in the second selected spectral region to the sum of the absorption spectra of all other isotopologues of the same chemical substance exceeds two somewhere in the second selected spectral region. These spectra are compared at the temperature and at the pressure of the measured gas or gas mixture.

(26) As an example, the rotational-vibrational band corresponding to the 2v.sub.1+2v.sub.2+v.sub.3 vibrational mode of .sup.12C.sup.16O.sub.2 and the rotational-vibrational band corresponding to the v.sub.1+2v.sub.2+v.sub.3 vibrational mode of .sup.13C.sup.16O.sub.2 are selected to illustrate the method. FIG. 11 shows the ratio of an absorption spectrum of the less abundant isotopologue, .sup.13C.sup.16O.sub.2, in the first selected spectral region to the sum of the absorption spectra of all other isotopologues of carbon dioxide at a pressure of 100 mBar and at a temperature of 296 K. FIG. 14 shows the ratio of an absorptions spectrum of the more abundant isotopologue, .sup.12C.sup.16O.sub.2, in the second selected spectral region to the sum of spectra of all other isotopologues of carbon dioxide at a pressure of 100 mBar and at a temperature of 296 K. P stands for gas pressure. All figures from 11 to 18 were created for gas pressure P=100 mBar, and gas temperature T=296K. S( ) is an absorption spectrum of a particular isotopologue. stands for sum, while indexes i, j, and k describe different isotopes of hydrogen, oxygen, or carbon atoms. One can see that at some spectral points corresponding to some rotational-vibrational lines the ratios exceed number two. The sums of spectra were calculated taking into account the abundances of different isotopologues.

(27) As another example, the R12 rotational-vibrational line of the v.sub.1+2v.sub.2+v.sub.3 vibrational mode of .sup.13C.sup.16O.sub.2 and the R10 rotational-vibrational line of the 2v.sub.1+2v.sub.2+v.sub.3 vibrational mode of .sup.12C.sup.16O.sub.2 can be selected for isotopic ratio analysis. The first spectral range to measure lines of .sup.13C.sup.16O.sub.2 is from 2040.5 nm to 2043 nm. The second spectral range to measure lines of .sup.12C.sup.16O.sub.2 is from 1572.5 nm to 1575 nm. Notice, these rotational-vibrational lines correspond to quite different transitions from one rotational level in the ground vibrational state to one rotational level in the vibrationally excited state. For example, the R12 rotational-vibrational line of the v.sub.1+2v.sub.2+v.sub.3 vibrational mode of .sup.13C.sup.16O.sub.2 corresponds to the transition (v.sub.1=0, v.sub.2=0, v.sub.3=0, J=12).fwdarw.(v.sub.1=1, v.sub.2=2, v.sub.3=1, J=13) from the rotation level J=12 of the ground state to the rotation level J=13 of the vibrationally excited v.sub.1+2v.sub.2+v.sub.3 vibrational mode, or to (v.sub.1=+1, v.sub.2=+2, v.sub.1=+1, J=+1) transition, while the R10 rotational-vibrational line of 2v.sub.1+2v.sub.2+v.sub.3 vibrational mode of .sup.12C.sup.16O.sub.2 corresponds to the transition (v.sub.1=0, v.sub.2=0, v.sub.3=0, J=10).fwdarw.(v.sub.1=2, v.sub.2=2, v.sub.3=1, J=11) from the rotation level J=10 of the ground state to the rotation level J=1 of the vibrationally excited 2v.sub.1+2v.sub.2+v.sub.3 vibrational mode, or to (v.sub.1=+2, v.sub.2=+2, v.sub.1=+1, J=+1) transition.

(28) According to a general convention rotational-vibrational bands can be divided on three branches: R-branch, P-branch, and Q-branch: R branch: when J=+1, i.e. the rotational quantum number in the ground state is one more than the rotational quantum number in the excited state; P branch: when J=1, i.e. the rotational quantum number in the ground state is one less than the rotational quantum number in the excited state; Q branch: when J=0, i.e. the rotational quantum number in the ground state is the same as the rotational quantum number in the excited state.

(29) The following tables show some parameters of the selected vibration modes of two isotpologues.

(30) TABLE-US-00001 TABLE 1 v.sub.1 + 2v.sub.2 + v.sub.3 band of .sup.13C.sup.16O.sub.2: branch symbol (P or R), J is the quantum number associated with the total angular momentum, wavelength [nm], line strength S [cm.sup.1/(molecule cm.sup.2)], air-broadened half-width .sub.air [cm.sup.1 atm.sup.1], self-broadened half-width .sub.self [cm.sup.1 atm.sup.1], and temperature dependence coefficient dS/SdT [K.sup.1]. J S .sub.air .sub.self dS/SdT P 38 2060.5951 2.14E22 0.0681 0.083 5.17E03 P 36 2059.7220 2.69E22 0.0683 0.085 4.12E03 P 34 2058.8614 3.32E22 0.0685 0.087 3.34E03 P 32 2058.0131 4.02E22 0.0687 0.089 2.30E03 P 30 2057.1772 4.78E22 0.0690 0.091 1.55E03 P 28 2056.3536 5.56E22 0.0693 0.093 8.33E04 P 26 2055.5421 6.36E22 0.0699 0.094 1.46E04 P 24 2054.7428 7.11E22 0.0706 0.096 5.22E04 P 22 2053.9556 7.79E22 0.0715 0.098 1.19E03 P 20 2053.1805 8.32E22 0.0727 0.099 1.68E03 P 18 2052.4173 8.67E22 0.0741 0.101 2.15E03 P 16 2051.6662 8.80E22 0.0758 0.103 2.64E03 P 14 2050.9269 8.66E22 0.0778 0.105 3.01E03 P 12 2050.1996 8.22E22 0.0798 0.107 3.40E03 P 10 2049.4841 7.48E22 0.0820 0.109 3.62E03 P 8 2048.7805 6.43E22 0.0842 0.112 3.91E03 P 6 2048.0886 5.11E22 0.0861 0.115 4.02E03 P 4 2047.4086 3.55E22 0.0877 0.117 4.20E03 P 2 2046.7403 1.83E22 0.0891 0.120 4.09E03 R 0 2045.7599 9.22E23 0.0925 0.125 4.34E03 R 2 2045.1210 2.75E22 0.0883 0.118 4.42E03 R 4 2044.4938 4.46E22 0.0870 0.116 4.18E03 R 6 2043.8783 6.01E22 0.0852 0.113 4.04E03 R 8 2043.2746 7.31E22 0.0831 0.110 3.96E03 R 10 2042.6826 8.32E22 0.0808 0.108 3.70E03 R 12 2042.1023 9.04E22 0.0788 0.106 3.40E03 R 14 2041.5338 9.41E22 0.0768 0.104 2.96E03 R 16 2040.9770 9.50E22 0.0749 0.102 2.93E03 R 18 2040.4319 9.32E22 0.0733 0.100 1.99E03 R 20 2039.8987 8.94E22 0.0721 0.099 1.66E03 R 22 2039.3773 8.35E22 0.0710 0.097 1.11E03 R 24 2038.8677 7.62E22 0.0702 0.095 4.88E04 R 26 2038.3699 6.80E22 0.0696 0.094 1.36E04 R 28 2037.8841 5.95E22 0.0691 0.092 7.79E04 R 30 2037.4101 5.11E22 0.0688 0.090 1.63E03 R 32 2036.9481 4.30E22 0.0686 0.088 2.37E03 R 34 2036.4981 3.55E22 0.0684 0.086 3.12E03 R 36 2036.0600 2.88E22 0.0682 0.084 4.17E03 R 38 2035.6341 2.30E22 0.0680 0.082 5.23E03

(31) TABLE-US-00002 TABLE 2 2v.sub.1 + 2v.sub.2 + v.sub.3 band of .sup.12C.sup.16O.sub.2: branch symbol (P or R), J is the quantum number associated with the total angular momentum, wavelength [nm], line strength S [cm.sup.1/(molecule cm.sup.2)], air-broadened half-width .sub.air [cm.sup.1 atm.sup.1], self-broadened half- width .sub.self [cm.sup.1 atm.sup.1], and temperature dependence coefficient dS/SdT [K.sup.1]. J S .sub.air .sub.self dS/SdT P 38 1584.0323 3.79E24 0.0682 0.082 5.20E03 P 36 1583.5076 4.80E24 0.0686 0.084 4.10E03 P 34 1582.9901 5.97E24 0.0689 0.086 3.30E03 P 32 1582.4799 7.28E24 0.0695 0.088 2.42E03 P 30 1581.9770 8.71E24 0.0698 0.090 1.67E03 P 28 1581.4815 1.02E23 0.0703 0.092 9.13E04 P 26 1580.9933 1.18E23 0.0709 0.094 0.00E+00 P 24 1580.5127 1.32E23 0.0716 0.096 7.84E04 P 22 1580.0395 1.45E23 0.0727 0.098 1.43E03 P 20 1579.5739 1.56E23 0.0741 0.100 1.33E03 P 18 1579.1158 1.64E23 0.0749 0.101 1.90E03 P 16 1578.6654 1.67E23 0.0766 0.104 2.49E03 P 14 1578.2224 1.65E23 0.0781 0.105 3.15E03 P 12 1577.7872 1.57E23 0.0797 0.108 3.31E03 P 10 1577.3596 1.43E23 0.0813 0.109 3.63E03 P 8 1576.9396 1.23E23 0.0830 0.112 4.21E03 P 6 1576.5273 9.80E24 0.0854 0.115 4.02E03 P 4 1576.1226 6.82E24 0.0873 0.118 4.10E03 P 2 1575.7256 3.51E24 0.0920 0.123 4.14E03 R 0 1575.1445 1.78E24 0.0953 0.129 4.08E03 R 2 1574.7667 5.28E24 0.0884 0.120 4.32E03 R 4 1574.3965 8.58E24 0.0858 0.116 4.11E03 R 6 1574.0340 1.15E23 0.0838 0.114 3.60E03 R 8 1573.6790 1.40E23 0.0816 0.111 3.70E03 R 10 1573.3317 1.59E23 0.0800 0.108 3.91E03 R 12 1572.9920 1.72E23 0.0781 0.106 3.01E03 R 14 1572.6598 1.80E23 0.0766 0.104 2.89E03 R 16 1572.3352 1.81E23 0.0747 0.102 2.29E03 R 18 1572.0180 1.77E23 0.0738 0.100 2.35E03 R 20 1571.7083 1.68E23 0.0727 0.098 1.85E03 R 22 1571.4060 1.56E23 0.0718 0.097 1.33E03 R 24 1571.1112 1.42E23 0.0710 0.095 7.30E04 R 26 1570.8236 1.26E23 0.0703 0.093 0.00E+00 R 28 1570.5434 1.10E23 0.0700 0.091 9.44E04 R 30 1570.2704 9.35E24 0.0694 0.089 1.66E03 R 32 1570.0046 7.82E24 0.0689 0.087 2.39E03 R 34 1569.7459 6.41E24 0.0685 0.085 3.23E03 R 36 1569.4943 5.16E24 0.0689 0.084 4.21E03 R 38 1569.2497 4.08E24 0.0679 0.081 5.08E03

(32) Table 1 and Table 2 show that both selected lines (R12 rotational-vibrational line of v.sub.1+2v.sub.2+v.sub.3 vibrational mode of .sup.13C.sup.16O.sub.2 and R10 rotational-vibrational line of 2v.sub.1+2v.sub.2+v.sub.3 vibrational mode) satisfy to another embodiment: the pressure broadening coefficients of these lines are different by no more than 50%.

(33) Table 1 and Table 2 also show that both selected lines (R12 rotational-vibrational line of v.sub.1+2v.sub.2+v.sub.3 vibrational mode of .sup.13C.sup.16O.sub.2 and R10 rotational-vibrational line of 2v.sub.1+2v.sub.2+v.sub.3 vibrational mode) satisfy yet to another embodiment: the temperature dependence coefficients of these lines are different by no more than 50%.

(34) If other inference species are present in the gas or in the mixture, similar spectral analysis may help better choose spectral regions. FIG. 12 and FIG. 13 show the ratios of the absorption spectrum of .sup.13C.sup.16O.sub.2 in the first selected spectral region to the sum of absorption spectra of all water isotopologues and to the sum of absorption spectra of all methane isotopologues at a pressure of 100 mBar and at a temperature of 296 K. FIG. 15 and FIG. 16 show the ratios of the absorption spectrum of .sup.12C.sup.16O.sub.2 in the second selected spectral region to the sum of absorption spectra of all H.sub.2O isotopologues and to the sum of absorption spectra of all CH.sub.4 isotopologues at a pressure of 100 mBar and at a temperature of 296 K. Concentrations of CO.sub.2, H.sub.2O, and CH.sub.4 in FIGS. 12-16 are the same. The sums were calculated taking into account abundances of different isotopologues.

(35) FIG. 17 and FIG. 18 show simulated absorptions spectra at two selected spectral regions. These figures show that by choosing rotational-vibrational lines which belong to different bands in different spectral regions, the instrument accuracy can be significantly improved: the analytical lines are not only relatively spectrally clean, but they also have close absorption coefficients and very close pressure broadening parameters and temperature dependencies. The first conditions permit to optimize the measurement method for both isotopologues. The last two conditions make these measurement less affected by both the pressure and temperature uncertainties.

(36) According to another embodiment, a gas analyzer system is provided for measuring an isotopic ratio gas species. The system typically includes a resonant optical cavity having two or more mirrors and containing a gas having a chemical species to be measured, a laser optically coupled to the resonant optical cavity, and a detector system for measuring absorption of laser light by the gas in the cavity.

(37) According to another embodiment, a gas analyzer system is provided for measuring an isotopic ratio gas species. The system typically includes a gas cell containing a gas to be measured, a laser optically coupled to the cell, and a detector system for measuring the laser light transmitted through the cell.

(38) According to another embodiment, a gas analyzer system is provided for measuring an isotopic ratio of gas species. The system typically includes an optical cavity containing a gas having a chemical species to be measured, a laser optically coupled to the optical cavity, and a detector system for measuring absorption of laser light by the gas in the cavity. In certain aspects, the gas analyzer system utilizes the cavity ring-down spectroscopy method to measure absorption of the laser light by the gas in the cavity.

(39) According to another embodiment, a gas analyzer system is provided for measuring an isotopic ratio of gas species. The system typically includes an optical cavity containing a gas having a chemical species to be measured, a laser optically coupled to the optical cavity, and a detector system for measuring absorption of laser light by the gas in the cavity. In certain aspects, the gas analyzer system utilizes the phase shift spectroscopy method to measure absorption of the laser light by the gas in the cavity.

(40) According to another embodiment, a gas analyzer system is provided for measuring an isotopic ratio of gas species. The system typically includes an optical cavity containing a gas having a chemical species to be measured, a laser optically coupled to the optical cavity, and a detector system for measuring absorption of laser light by the gas in the cavity. In certain aspects, the gas analyzer system utilizes the cavity enhanced absorption spectroscopy method to measure absorption of the laser light by the gas in the cavity.

(41) According to another embodiment, a gas analyzer system is provided for measuring an isotopic ratio of gas species. The system typically includes an optical cavity containing a gas having a chemical species to be measured, a laser optically coupled to the optical cavity, and a detector system for measuring absorption of laser light by the gas in the cavity. In certain aspects, the gas analyzer system utilizes the photoacoustic spectroscopy method to measure absorption of the laser light by the gas in the cavity.

(42) According to another embodiment, a gas analyzer system is provided for measuring an isotopic ratio of gas species. The system typically includes an optical cavity containing a gas having a chemical species to be measured, a laser optically coupled to the optical cavity, and a detector system for measuring absorption of laser light by the gas in the cavity. In certain aspects, the gas analyzer system utilizes the tunable diode lasers spectroscopy method to measure absorption of the laser light by the gas in the cavity.

(43) According to another embodiment, a gas analyzer system is provided for measuring an isotopic ratio of gas species. The system typically includes an optical cavity containing a gas having a chemical species to be measured, a laser optically coupled to the optical cavity, and a detector system for measuring absorption of laser light by the gas in the cavity. In certain aspects, the gas analyzer system also includes a temperature sensor for measuring the temperature of the gas in the cavity, and a pressure sensor for measuring the pressure of the gas in the cavity. In certain aspects, the detector system includes one of a photo-detector configured to measure an intensity of the intra-cavity light or both a photo-acoustic sensor configured to measure photo-acoustic waves generated in the cavity and a photo-detector configured to measure an intensity of the intra-cavity light.

(44) According to yet another embodiment, a gas analyzer system is provided for measuring an isotopic ratio gas species. The system typically includes an optical cavity containing a gas having a chemical species to be measured, a laser optically coupled to the optical cavity, and a detector system for measuring the absorption of laser light by the gas in the cavity. In certain aspects, the gas analyzer system also includes a control element configured to control temperature of the gas in the optical cavity and a pressure control element configured to control pressure of the gas in the optical cavity.

(45) All references, including publications, patent applications, and patents, cited herein are hereby incorporated by reference to the same extent as if each reference were individually and specifically indicated to be incorporated by reference and were set forth in its entirety herein.

(46) The use of the terms a and an and the and at least one and similar referents in the context of describing the disclosed subject matter (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. The use of the term at least one followed by a list of one or more items (for example, at least one of A and B) is to be construed to mean one item selected from the listed items (A or B) or any combination of two or more of the listed items (A and B), unless otherwise indicated herein or clearly contradicted by context. The terms comprising, having, including, and containing are to be construed as open-ended terms (i.e., meaning including, but not limited to,) unless otherwise noted. Recitation of ranges of values herein are merely intended to serve as a shorthand method of referring individually to each separate value falling within the range, unless otherwise indicated herein, and each separate value is incorporated into the specification as if it were individually recited herein. All methods described herein can be performed in any suitable order unless otherwise indicated herein or otherwise clearly contradicted by context. The use of any and all examples, or example language (e.g., such as) provided herein, is intended merely to better illuminate the disclosed subject matter and does not pose a limitation on the scope of the invention unless otherwise claimed. No language in the specification should be construed as indicating any non-claimed element as essential to the practice of the invention.

(47) Certain embodiments are described herein. Variations of those embodiments may become apparent to those of ordinary skill in the art upon reading the foregoing description. The inventors expect skilled artisans to employ such variations as appropriate, and the inventors intend for the embodiments to be practiced otherwise than as specifically described herein.