Infrared spectrometer

<u>Field of the Invention</u>

0001The present invention relates to a infrared, IR, spectrometer and in particular to a tunable IR laser spectrometer.

0002The present invention also relates to a method of isotopic ratio determination.

<u>Background of the Invention</u>

0003Stable isotopomer concentration and flux measurements have become a vital underpinning technique in many areas of science. To date, mass spectrometry (MS) has been the preferred technique to perform isotope analysis. For any application requiring real time, high frequency data, and/or field measurements, MS suffers from drawbacks: tedious sample taking and preparation, difficulty of real time analysis. Furthermore, MS does not resolve small mass differences particularly well. The use of MS is confined to a dedicated laboratory and field applications are therefore difficult. There is a growing need to overcome these drawbacks and to seek alternative instrumentation for high precision (∼0.1‰) isotopic ratio determination both in the academic and industrial sectors. Applications encompass atmospheric studies, geology, ecology, medical research, planetary exploration, combustion science, fundamental analytical chemistry, and food industry quality check.

0004Therefore, there is requirement for a method and apparatus for providing high precision isotopomer concentration and ratio measurements without the above drawbacks.

0005A doctoral dissertation by Michael Haisch (Heinrich Heine University of Dusseldorf, 1995) describes a spectrometer including a light source for generating a first and second beam. The first beam is passed through a first sample cell and the second beam is passed through a second sample cell. The beams are then passed to a chopper such that the first and second beams are alternately received at a detector.

0006<patcit id="pcit0001" dnum="JP2003014631A"><text>JP 2003-14631</text></patcit> describes a spectrometer having a light source, wherein the beam from the light source is split. The first portion of the split beam is received by the sample cell, and the second portion of the split beam follows a pathway which avoids the sample cell. A chopper causes the first and second portions of the beams to be alternately received by a detector.

SUMMARY OF THE INVENTION

0007Generally, the present invention provides a method and apparatus for detecting, by absorption spectroscopy, an isotopic ratio of a sample, by passing first and second laser beams of different frequencies through the sample. In a particular example, this may be applied to the determination of ¹²CO₂ :¹³CO₂ or similar. In a first aspect of the invention, there is an infrared, IR, spectrometer according to claim 1.

0008Typically, two IR absorption cells are used, a first containing a reference gas of known isotopic ratio and the second containing a sample of unknown isotopic ratio.

0009An interlacer, switcher or reflective chopper is used so that as the laser frequencies are scanned the absorption of the sample cell and the reference cell are detected alternately. This ensures that the apparatus is continuously calibrated and rejects the baseline noise when phase sensitive detection is used.

0010Preferably, the rate at which the sample cell and reference cell are alternatively detected is faster than several hundred Hz and even more preferably, above 1 KHz. This improves the signal-to-noise of the system. The wavelength of the lasers may be scanned in synchronisation with this alternate detection frequency. For instance, the full range (or a range of wavelengths between particular absorption bands of the sample and reference) of laser wavelengths may be scanned between each detection cycle. Phase sensitive detection may be used to maintain this synchronisation.

0011According to an example there is provided a method of selecting a pair of isotopic spectroscopic lines for a sample suitable for use in isotopic ratio determination comprising the step of: (a) selecting a pair of isotopic spectroscopic lines of similar intensities such that a suitable pair of isotopic spectroscopic lines is selected.

0012Preferably, the method further comprises the step of: (b) comparing the thermal characteristics of said pair of isotopic spectroscopic lines and rejecting said pair if the thermal characteristics are significantly different.

0013Advantageously, the thermal characteristics are determined to be significantly different if the line intensity or wavelength varies significantly with temperature.

0014Preferably, the isotopic spectroscopic lines are spaced far apart.

0015Preferably, the sample is selected from the group consisting of ¹²CO₂/¹³CO₂; C¹⁶O₂/¹⁶OC¹⁸O, H₂¹⁶O/H₂¹⁸O, ¹²CH₄/¹³CH₄, ¹²CH₄/¹²CDH₃, ¹⁴N₂O/¹⁴N¹⁵NO and H₂O/HDO.

0016According to a further aspect of the present invention there is provided a method for detecting an isotopic ratio of a sample according to claim 12. In an illustrative example, not forming part of the claimed invention, there is a method comprising the steps of (a) passing a first laser beam and a second laser beam of different frequencies through a sample; and (b) detecting the optical absorption due to the sample such that a first absorption line and a second absorption line are measured.

0017According to a further illustrative example, there is provided a method for detecting an isotopic ratio of a sample by measuring the relative intensities of at least one of the following pairs of absorption lines: 3601.4210 cm^-1 / 2294.4811 cm^-1, 3599.7027 cm¹/ 2295.8456 cm^-1 and 3597.9626 cm^-1 / 2297.1862 cm^-1.

0018Further details of the present invention are described in the accompanying claims.

<u>Brief description of the Figures</u>

0019The present invention may be put into practice in a number of ways and an embodiment will now be described by way of example only and with reference to the accompanying drawings, in which: <ul id="ul0001" list-style="none"><li><figref idref="f0001">FIG. 1</figref> shows a schematic diagram of an IR laser spectrometer according to a first embodiment of the present invention, given by way of example only;</li><li><figref idref="f0002">FIG. 2</figref> shows a timeline of the data acquired from the IR laser spectrometer of <figref idref="f0001">FIG. 1</figref>;</li><li><figref idref="f0003">FIG. 3</figref> shows a flow-chart of a method of collecting data from the IR laser spectrometer of <figref idref="f0001">FIG. 1</figref>;</li><li><figref idref="f0004">FIG. 4</figref> shows an IR laser spectrometer according to an illustrative example, not forming part of the claimed invention and given by way of example only;</li><li><figref idref="f0005">FIG. 5</figref> shows a set of graphs depicting pairs of suitable spectral lines to be investigated by the IR laser spectrometer of <figref idref="f0001">FIG.1</figref> or <figref idref="f0004">FIG. 4</figref>;</li><li><figref idref="f0006">FIG. 6a</figref> shows a graph illustrating conventional phase sensitive detection in tunable spectroscopy using a chopper; and</li><li><figref idref="f0006">FIG. 6b</figref> shows a graph illustrating the signal obtained using the IR laser spectrometer of <figref idref="f0001">FIG. 1</figref>.</li></ul>

<u>Detailed description of an embodiment</u>

0020The inventor of the present invention has described carbon isotopomer ratio measurements in "<nplcit id="ncit0001" npl-type="s"><text>Carbon isotopomers measurement using mid-IR tunable laser sources" Damien Weidmann, et al, Isotopes in Environmental and Health Studies, Vol. 41, No. 4, December 2005, pp293-302</text></nplcit> (DOI: 10.1080/10256010500384325).

0021Absorption spectroscopy techniques, using infrared laser sources, offer an improvement over MS. A useful spectral region for absorption spectroscopy is known as the "fingerprint" region lying in the mid infrared range (2-25 µm), where most molecules exhibit intense vibrational transitions.

0022As an example, identification and measurement of gases released by volcanoes provide information on magmatic sources and so aid interpretation of volcanic behaviour and eruption forecasting. As well as information on real-time trace gas concentrations in volcanic emissions, isotopic ratio measurements of species including carbon, hydrogen, nitrogen, oxygen and sulphur provide additional information concerning the gas sources and their geochemical history. Stable carbon isotopomers, especially ¹²CO₂ and ¹³CO₂, are of particular interest. Carbon dioxide is typically the most abundant gas released by volcanoes after H₂O. ^13/12CO₂ monitoring can be used as a tracer to discriminate the sources contributing to observed CO₂ emissions.

0023<figref idref="f0001">FIG. 1</figref> shows a schematic diagram of an IR spectrometer 10 according a first embodiment of the present invention.

0024Diode laser DL is a room temperature tunable antimonide laser operating around 3601 cm^-1. Quantum cascade laser QCL is a room temperature tunable laser operating around 2294 cm^-1. Off-axis parabolic mirror OAPM1 collimates the beam from the QCL and this collimated beam is combined with the beam from DL using a dichroic mirror DM1.

0025The combined beam is then split into two paths by a beam splitter BS before each portion of the beam enters a dual channel absorption cell providing the sample-laser interaction volume. Mirror M1 directs one of the split beams into channel 20, whilst the portion passing through the beam splitter enters channel 30 directly. The beams may be split 50:50 or by any other ratio. Channel 20 contains a reference sample of known composition and channel 30 contains a sample under investigation. For instance channel 20 may contain a mixture of ¹²CO₂/¹³CO₂ of known proportion and channel 30 may contain a mixture of ¹²CO₂/¹³CO₂ of unknown proportion. The pressure, flow rate and temperature of the sample and reference may be controlled and monitored during absorption measurements.

0026As two wavelengths probe each cell channel, the system performs measurements over four spectroscopic channels.

0027Detector 1 and detector 2 detect transmitted laser light originating from the DL and QCL respectively. Detector 1 and detector 2 are each respectively sensitive to the wavelengths provided by the DL and QCL lasers. Dichroic mirror DM2 separates the combined beam ensuring that detector 1 receives light originating from the DL and detector 2 receives light originating from the QCL. Off-axis parabolic mirrors QAPM1 and QAPM2 focus the collimated beams onto detectors 1 and 2 respectively.

0028At the output of the dual channel absorption cell a 10KHz chopper 40 with reflective blades plays the role of an interlacer. Mirror M2 directs the output of channel 20 onto the chopper blades. With the chopper 40 presenting an opening 50 towards the output of the dual absorption cell the output of channel 20 (reference) passes straight through and is not directed towards the two detectors. However, the output from channel 30 (sample) also passes straight through the chopper and is directed onto the two detectors.

0029With the chopper 40 presenting a reflective blade towards the output of the dual absorption cell the output of channel 30 is blocked by the back of the mirror and the output from channel 20 (reference) is directed onto the two detectors. In this way detectors 1 and 2 sequentially receive IR radiation transmitted through each channel 20, 30 in turn. An acquisition unit collects and stores data from each of the two detectors. The acquisition unit may be in turn connected to a computer system for data manipulation and storage. Alternately, the computer system itself may act as the acquisition unit.

0030The DL and QCL lasers are each tunable across a range of wavelengths. DL control controls the output wavelength of diode laser DL. QCL control controls the output wavelength of the QCL laser. The wavelengths of each laser are scanned across a particular range in order to obtain an absorption spectrum from the reference and sample with intensity being measured by the two detectors. The tuning of each laser is synchronised such that a step change intensity occurs after each cycle (open, reflective blade) of the chopper. Alternatively, instead of a step change, a smooth continuous wavelength change may be made.

0031The alternate sampling of the reference cell 20 and sample cell 30 provides an effective onboard calibration.

0032<figref idref="f0006">FIG. 6a</figref> shows a schematic illustration of conventional phase sensitive detection in tunable spectroscopy using a chopper. The signal peak to peak amplitude is S₀ and only a small variation dS of this amplitude contains useful information.

0033<figref idref="f0006">FIG. 6b</figref> shows a schematic illustration of the signal obtained from the present invention using a reflective chopper. The signal peak to peak amplitude is dS. AC levels of amplification can be much higher as only the change is monitored. Other drifting parameters cancel out using such a technique with phase sensitive detection.

0034<figref idref="f0002">FIG. 2</figref> shows the timeline showing the synchronisation between the laser tuning for each laser (having wavelength λ₁ and λ₂, respectively) and the chopper state, i.e. which of cell 20 (reference) or 30 (sample) is being measured by detectors 1 and 2. Note that each laser wavelength changes at the same time once a reference/sample cycle has been completed. It is not necessary that the wavelength of each laser changes at the same rate. Alternatively, a smooth variation of wavelength may occur throughout the cycle.

0035Isotopic ratio is usually expressed in terms of delta value defined by equation 1: <maths id="math0001" num="Equation 1"><math display="block"><mrow><mi>δ</mi><mo>=</mo><mfenced separators=""><mfrac><mrow><msub><mi>R</mi><mi>x</mi></msub></mrow><mrow><msub><mi>R</mi><mi>s</mi></msub></mrow></mfrac><mo>−</mo><mn>1</mn></mfenced><mo>×</mo><mn>1000</mn></mrow></math><img file="EP2247929B1_D0001.tif" /></maths> where R_x denotes the ratio of the heavier isotope to the lighter one in the sample to be measured, and R_s refers to the ratio of a calibrated reference sample. Units are "per mil" or ‰.

0036<figref idref="f0003">FIG. 3</figref> shows a flow diagram of a method of obtaining a relative delta value, δ_rel of a CO₂ sample. In this case: <maths id="math0002" num="Equation 2"><math display="block"><mrow><msub><mi>δ</mi><mi mathvariant="italic">rel</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>P</mi><mn>13</mn></msub></mrow><mrow><msub><mi>P</mi><mrow><mi>S</mi><mn>12</mn></mrow></msub></mrow></mfrac><mo>−</mo><mn>1</mn></mrow></math><img file="EP2247929B1_D0002.tif" /></maths>

0037<figref idref="f0004">FIG. 4</figref> shows a QCL-based two channel spectrometer 200 which is provided as an illustrative example only, and is not an embodiment of the claimed invention. This two channel spectrometer allows isotopomer measures to be made for certain pairs of spectroscopic lines but has limited functionality when compared with spectrometer 10 of <figref idref="f0001">FIG. 1</figref>. The identified target lines for this sensor may be located at 2311.105566 cm^-1 (¹²CO₂, v3) and 2311.398738 cm^-1 (¹³CO₂, v3), respectively. With this selection, the δ¹³C measurements are insensitive to gas temperature variations. However, the line intensities differ by a factor of about 100, and therefore a dual path length absorption cell is required. The sensor has already been described in "<nplcit id="ncit0002" npl-type="s"><text>Development of a compact quantum cascade laser spectrometer for field measurements of CO2 isotopes", Appl. Phys. B, 80, pp255-260 (2005</text></nplcit>), incorporated by reference herein. The spectroscopic source is a pulsed, thermoelectrically cooled, DFB QC laser excited by 25 ns current pulses. A dual path length absorption cell 110 is an astigmatic Herriott cell modified to incorporate an additional short path. After exiting the cells, the two beams are directed to two thermoelectrically cooled HgCdZnTe detectors. The complete optical subsystem is housed inside a sealed Delrin case and purged with dry nitrogen to avoid contamination from atmospheric CO₂ or corrosive gases. The associated sensor electronics are located below the optical platform. A gas control subsystem alternating the flow of sample and calibration gases into the absorption cell is connected to the instrument. The operating wavelength of the QCL obtained for this system was found to be ∼8 cm^-1 distant from the selected frequencies. Also, an unusually large threshold current (9 A) prevented the use of a fast (∼400 Hz) subthreshold current ramp to rapidly tune the laser wavelength. A¹²CO₂ line at 2320.7501 cm^-1 and a ¹⁶O¹²C¹⁸O line at 2320.4599 cm^-1 may be targeted with this device. The intensity ratio of these lines is ∼20 The extrapolated precision of δ¹³C is ∼1 ‰. Because a single laser source is used in spectrometer 200 only close absorption lines may be investigated. Furthermore, the use of a single cell requires separate reference calibration once a sample has been purged.

0038Before spectroscopic measurements may be taken a suitable pair of spectroscopic absorption lines must be identified. In order to avoid the need for different path lengths within the absorption cells (such as used in spectrometer 200) lines of similar intensity may be chosen. As this technique relies on measuring absolute absorption line intensities the overriding requirement is to avoid lines which behave differently to temperature changes. Ideally, the relative change in intensity for each line in a pair when the temperature varies will be the same. With single laser spectrometers the choice of lines is limited to lines very close together. This leads to difficulty in resolving the lines and usually results in large differences in absolute line intensities. However, with the availability of the dual lasers system lines very far apart in wavelength terms may be used as a pair.

0039A line selection method, which consists of choosing optimum absorption lines for the two isotopomers under study, involves the consideration of several issues.

0040The line intensities have to be considered. Ideally, similar intensities would be sought so that the same cell can be used (i.e. same path length), and to limit potential non-linearities in the detection system. The line intensities should also be high enough to ensure the best sensor sensitivity and to allow the use of a compact absorption cell. The temperature stability requirement depends on the Boltzmann distribution of the tworespective isotopic transitions. A similar evolution of the two line intensities with temperature is important to limit sensitivity to temperature fluctuations. This condition is expressed equation 3: <maths id="math0003" num="Equation 3"><math display="block"><mrow><mi mathvariant="normal">Δ</mi><mi mathvariant="italic">T</mi><mo>•</mo><mi mathvariant="normal">Δ</mi><mi mathvariant="italic">E</mi><mo>=</mo><mi mathvariant="normal">Δ</mi><mi mathvariant="italic">δ</mi><mo>•</mo><mi>k</mi><mo>•</mo><msup><mi>T</mi><mn>2</mn></msup></mrow></math><img file="EP2247929B1_D0003.tif" /></maths> where Δδ is the target accuracy on δ, k the Boltzmann constant, T the temperature, ΔE the energy difference between the lower levels of the two transitions and ΔT the temperature variation. Hence, in terms of spectroscopic parameters to obtain temperature insensitivity, the lower energy levels must be identical.

0041The potential spectral and collisional interferences by other species have to be considered. In particular, it is important to avoid interference due to water vapour. Well isolated lines are preferred, making concentration retrievals easier. Other parameters may be dictated by the spectroscopic source to be used. For example, ensuring that the two lines are within the available tuning range and that the source linewidth matches the absorption widths of the selected isotopic line pair.

0042To implement line selection for any pair of isotopomers of any molecule, an algorithm based on the aforementioned method may be used. First, depending on the experimental conditions, the minimum line intensities to be considered may be defined (cutoff intensity, S_min). Then, precision parameters may be entered into the algorithm: the minimum acceptable temperature stability ΔT , the relative difference allowed between the two line intensities ΔS and the minimum frequency separation required between the two lines Δv. All the possible line pairs of the two selected isotopomers within the available database are then processed. The identified line pairs then have to be critically examined, taking into account the other aforementioned selection criteria.

0043The algorithm may be applied to potential ¹²CO₂, ¹³CO₂ line pairs or any other suitable molecule. The minimum intensity may be set to 10-21 cm^-1/molec cm^-2. This intensity criterion is low for atmospheric studies, but sufficient for volcanic environments where CO₂ concentrations can exceed several percent.

0044The parameters may be chosen as follows: a target delta value accuracy of 0.1‰, a temperature stability of better than 0.5 K, a relative difference in intensities of &lt;10% and line separation of at least 0.05 cm^-1. A stringent condition on line intensities is required to achieve the expected accuracy on the delta value.

0045<figref idref="f0005">FIG. 5</figref> shows the results of applying the above algorithm to the spectroscopic lines of ¹²CO₂ and ¹³CO₂. The lower graph provides details of the three ro-vibrational band intersections appearing in the main graph: (ν3, ν3), (ν1, ν3), (2ν1 + ν3, ν2), from left to right. One hundred and eighty-eight pairs of transitions satisfy these criteria. The points are located in three regions, detailed in the lower part of <figref idref="f0005">FIG. 5</figref>, corresponding to the intersection of the ¹²CO₂ and ¹³CO₂ bands indicated on the upper graph.

0046<figref idref="f0005">FIG. 5</figref> indicates that the optimum selection involves two lines from two different bands, thus requiring the two laser sources of the present invention. With the parameters used to establish this selection, only two pairs belong to the same band (ν3) and the smallest frequency separation is found to be ∼25 cm^-1 for 2336.5590 cm^-1 (¹³CO₂) and 2311.3990 cm^-1 (¹²CO₂). Such a tuning range is not accessible from a single semiconductor mid-IR sources.

0047By choosing more stringent parameters such as S_min = 10-20 cm^-1/molec cm^-2, ΔT = 10K, ΔS = 1%, only three line pairs are found. They are listed in table 1. <tables id="tabl0001" num="0001"><table frame="all"><title>Table 1. Optimum line selection for ^23/12CO₂ ratio measurements.</title><tgroup cols="6"><colspec colnum="1" colname="col1" colwidth="29mm" /><colspec colnum="2" colname="col2" colwidth="41mm" /><colspec colnum="3" colname="col3" colwidth="14mm" /><colspec colnum="4" colname="col4" colwidth="29mm" /><colspec colnum="5" colname="col5" colwidth="41mm" /><colspec colnum="6" colname="col6" colwidth="13mm" /><thead><row><entry namest="col1" nameend="col3" align="center" valign="middle">¹²CO₂</entry><entry namest="col4" nameend="col6" align="center" valign="middle">¹³ CO₂</entry></row><row><entry align="center" valign="middle">Frequency (cm^-1)</entry><entry align="center" valign="middle">Strength (cm^-1/molec cm^-1)</entry><entry align="center" valign="middle">Band</entry><entry align="center" valign="middle">Frequency (cm^-1)</entry><entry align="center" valign="middle">Strength (cm^-1/molec cm^-2)</entry><entry align="center" valign="middle">Band</entry></row></thead><tbody><row><entry align="center" valign="middle">3601.4210</entry><entry align="center" valign="middle">3.513e-020</entry><entry align="center" valign="middle">ν1 + ν3</entry><entry align="center" valign="middle">2294.4811</entry><entry align="center" valign="middle">3.547e-020</entry><entry align="center" valign="middle">V3</entry></row><row><entry align="center" valign="middle">3599.7027</entry><entry align="center" valign="middle">3.573e-020</entry><entry align="center" valign="middle">ν1 + ν3</entry><entry align="center" valign="middle">22955.8456</entry><entry align="center" valign="middle">3.574e-020</entry><entry align="center" valign="middle">ν3</entry></row><row><entry align="center" valign="middle">3597.9626</entry><entry align="center" valign="middle">3.525e-020</entry><entry align="center" valign="middle">ν1 + ν3</entry><entry align="center" valign="middle">2297.1862</entry><entry align="center" valign="middle">3.499e-020</entry><entry align="center" valign="middle">ν3</entry></row></tbody></tgroup></table></tables>

0048These three line pairs offer the best strategy to achieve a high accuracy for δ¹³C measurements.

0049The line pair selection may be made so that both the two intensities and the two absorption frequencies are close. Inherently, the lower state energies of the selected line pair are different, ΔE ∼ 1160 cm^-1. The expected required temperature stability, calculated with equation 3, is 6 mK. As a consequence, pressure change has a strong effect on the peak amplitude and overlapping wing contributions from one line to the other.

0050The present invention may be used in: <ul id="ul0002" list-style="none" compact="compact"><li>Geophysical studies: Real time continuous monitoring of isotopomers released by volcanic fumes. The knowledge of isotope concentrations allows the deduction of information about the different sources feeding a volcanic system. The present invention strengthens the reliability of volcanic eruption forecasting. Also oil prospecting can be facilitated by identifying carbon sources using this method.</li><li>Atmospheric &amp; environmental sensing: The present invention allows the provision of high temporal and high spatial resolution ¹³C/¹²C data, thus providing further constraint to carbon budget models. This allows a fuller understanding and account for carbon sources and sinks.</li><li>Planetary exploration: The present invention may be used in planetary exploration, especially Mars exploration.</li><li>Medical diagnostics: Real time, continuous monitoring of trace molecule enables non-invasive breath diagnostics. For instance, the activity of Pylobacter Pylori (bacteria responsible of stomach ulcer) may be identified by the measurement of ¹²CO₂/¹³CO₂ ratio in breath. ¹²C/¹³C ratio may also discriminate between a catabolic and anabolic state of living cells.</li><li><b>Agronomy and food industry:</b> Isotopomers concentration contained in a given product may indicate the product origin. Real time, continuous monitoring of isotopic ratio by the present invention may provide on-line implementation, and detection of artificial substitutes.</li></ul>

0051As will be appreciated by the skilled person, details of the above embodiments may be varied without departing from the scope of the present invention.

0052In a particular illustrative example, the chopper may be placed before the absorption cells 20, 30 to alternately block or switch the incident IR beams rather than blocking the transmitted beams. In this alternative example the beam splitter may not be required resulting in all of the available laser intensity being incident on one cell at a time. Obviously, the detector arrangement will require adjustment accordingly.

0053In an illustrative example only, a single detector may be used provided that it can detect both wavelengths.

0054Although a 10 kHz chopper frequency has been described, any suitable frequency may be used.

0055Preferred embodiments and methods of the present invention discussed in the foregoing are to be understood as descriptions for illustrative purposes only, and it will be appreciated that numerous changes, substitutions, omissions, and updates thereof are possible without departing from the scope of the claims.