Abstract

Laser-spectrometric methods to derive absolute and traceable carbon monoxide (CO) amount fractions in exhaled human breath could be of advantage for early disease detection as well as for treatment monitoring. As proof-of-principle laboratory experiment, we employed intra-pulse and continuous wave (cw) quantum cascade laser spectroscopy (QCLAS), both at 4.6 μm. Additional experiments were carried out applying cw cavity ring-down spectroscopy (CRDS) with a CO sideband laser and a QCL. We emphasize metrological data quality objectives, thatis, traceability and uncertainty, which could serve as essential benefits to exhaled breath measurements. The results were evaluated and compared on a 100 μmol/mol CO level using the two QCLAS spectrometers, and the cw CO sideband laser CRDS setup. The relative standard uncertainties of the pulsed and the cw QCLAS CO amount fraction results were ±4.8 and ±2.8%, respectively, that from the CO sideband laser CRDS was ±2.7%. Sensitivities down to a 3 nmol/mol CO level were finally demonstrated and quantified by means of cw CRDS equipped with a QCL yielding standard uncertainties of about ±2.5 that are exclusively limited by the available line strength figure quality. With this study we demonstrate the achieved comparability of CO quantifications, adhering metrological principles.

1. Introduction

The last years exhibited the development of new laser sources such as quantum cascade lasers [1] operating in the mid infrared where the fundamental bands of most infrared active molecules are located. In breath analysis for instance, mid infrared light sources have been used to measure the amount fraction of biomarkers such as acetone or carbon monoxide (CO) found in exhaled human breath using different laser spectroscopic techniques with detection limits down to the pmol/mol level [2].

Carbon monoxide, which we focus on in this work, is in air a pollutant resulting from the incomplete burning of carbon-containing fuels. As a process product it can reach quite large concentration levels of several tens of μmol/mol, expressed as amount fractions. CO in breath of smokers could be found in a range of up to a few μmol/mol, which is well above atmospheric levels (typically 100 nmol/mol). Patients suffering from anaemias, oxidative stress, and respiratory infections have been found with abnormal levels of carbon monoxide. Therefore, CO is discussed as being a biomarker for these diseases. However, because of the low μmol/mol level of CO reported to be present in exhaled breath and because of the presence of several molecules at the same or even larger amount of substance fraction levels, a very selective and sensitive technique is required to measure the amount fraction of breath CO [2, 3]. Laser spectrometric techniques are seemly unique to perform such measurements because of their robustness from sample preparation to measurement, selectivity, and sensitivity [2].

Laser spectroscopic techniques such as tunable diode laser absorption spectroscopy (TDLAS), quantum cascade laser absorption spectroscopy (QCLAS), and cavity ring-down spectroscopy (CRDS) applied to gas detection had been reported to be selective, sensitive, and have near real-time and point of care capabilities [2].

Employing QCLAS and CRDS, which is the focus of this manuscript, to perform absolute amount fraction measurements of different molecular species had been reported by several groups [49]. Some scientists have reported CO amount fraction results based on QCLAS and CRDS [2, 10, 11]. What we found missing in most of the reports is information on metrological data qualifiers, that is, traceability statements or uncertainty budgets and a methodological standardization attempt where performed measurements could be referred to.

In this work, we report laboratory-based amount fraction measurements of CO in gravimetrically prepared gas mixtures, as artificial breath samples, for example, comprised of CO2, O2, CH4, and N2 as matrix gas. We applied QCLAS operated with both, intrapulse and contiunous wave (cw) QCLs, and cw-CRDS operated by a QCL and a CO sideband laser setup, respectively. The CO sideband laser setup has proved to be a suitable tool for CO breath analysis as shown in [12]. The measurements were performed in the framework set by an iMERA-Plus joint research project on “breath analysis” [13]. The basic idea of the measurement method used was promoted in a previous EURAMET project and referred to as “traceable infrared laser spectrometric amount fraction measurement” (TILSAM) method [14]. Its application to spectroscopy techniques with potentials for breath analysis (BA) was projected in an extended BA-TILSAM descriptive document [13].

According to the TILSAM method, the aim of this paper was to rely on the quality of the input quantities such as the gas pressure (p), gas temperature (T), and the optical path length through the absorbing medium (L) to directly derive absolute CO amount fractions rather than to calibrate with reference gas mixtures. The uncertainties were estimated referring to the “guide to the expression of uncertainty in measurement” (GUM) [15].

This manuscript is structured as follows. We briefly outline the background of QCLAS similar to what we reported in one of our earlier publications [9]. We then present the experimental setups, measurements, and the results of pulsed-QCLAS, cw-QCLAS, and CRDS in Sections 3, 4, and 5. The results from the three spectrometric techniques are compared and discussed in Section 6. Some conclusions are presented in the last section. In contrast and in addition to studies aiming towards improved sensitivities and assessments of noise levels [16, 17], this paper aims mainly to compare CO quantification results on a 100 μmol/mol level, first, and second, to demonstrate some few nmol/mol CO measurement capabilities. All the quantifications refer to the same absolute method.

2. Conceptual Background

For absorption spectroscopy, the interaction of the gas molecules and the sensing radiation at wave number is modeled by the Beer-Lambert law with and being the incident and transmitted radiant powers, respectively, of which the SI unit is W. The absorption is governed by the molecular transition line strength at gas temperature T, the respective normalized absorption profile centered at , and the absorption path length . The quantity is the isotopic composition factor, given for example, as , for a probed 12C16O excitation, where and are the abundances of 12C16O in the sample and the conventional abundance value given by HITRAN [18], respectively. The line strength is specific for the probed molecular transition. Relying on the ideal gas law, the molecular density of the absorbing species can be expressed in terms of the partial pressure of the absorbing molecular species and the gas temperature. The partial pressure can be related to the total pressure using the amount of substance fraction of the absorbing species, . Thus, (1) becomes While probing a certain molecular transition measuring and for a known path length, as well as and , leads to the amount of substance fraction of the species Introducing the spectral absorbance , which in some cases (see, e.g., http://goldbook.iupac.org/A00028.html) can also be called extinction, and making use of the normalization of , (2) can also be written in its integral form where is the line area obtained by integration of the measured absorbance data over wave number. Direct absorption spectroscopy is described by the Beer-Lambert law, as expressed in (2). All input quantities are SI traceable. In turn, using (4) to determine the amount of substance fraction, , in a given gas mixture would deliver results that are SI traceable. An uncertainty budget for the measurement according to the ISO Guide to the expression of uncertainty in measurement (GUM) [15] can be evaluated with (4) as model function.

3. Intrapulse QCLAS

Pulsed QCLs can be operated in inter- and intrapulse mode [11, 1921]. A comparison of the two modes of operation is found in [22, 23]. Here, we focus on the intrapulse mode operation of a QCL to perform absolute amount fraction measurements of CO. Measurements and results in this section are reported in a similar manner as in [9]. A single pass gas cell was used with an optical path length of 0.82 m [9].

Briefly, the intrapulse chirp mode of operation utilizes the intrinsic temperature rise inside the active area of the QCL chip within each single pulse caused by the current driven through the chip [24]. Right from the beginning of each pulse the emission wavelength increases within several nanoseconds. Thus, the laser emission spans a certain pulse-length-dependent spectrum within a single pulse. The laser wavelength sweeps from the blue towards the red side of the spectrum. By setting the laser temperature and the laser voltage, chirp-onset wavelength and laser power can be tuned. The sweep width is then set by the pulse length. In principle the intrapulse mode spectroscopy allows us to study the molecular spectrum of interest within one single laser pulse of a few nanoseconds, what could be important, for example, in fast process studies.

Figure 1 depicts a schematic of the experimental QCLAS setup. The QCL (Fraunhofer-IAF) emitted at 4.6 μm and was operated with pulse length of 255 ns (2 kHz repetition rate, substrate temperature 292 K). After collimation, the QCL beam was split into a reference beam and a second sample beam was directed through a 0.82 m stainless steel absorption cell. The reference beam was used for intensity normalization and for fixing the wavelength axis by means of a removable etalon. Both channels were terminated into TE-cooled HgCdZnTe detectors (Vigo PDI-2TE-10.6). The electrical signals were high-speed digitized at 5 GS/s by an 8-bit A/D-converter (Compuscope 85 G) to resolve the temporal evolution within single pulses.

For gas detection, the P(1) line of CO at 2139.43 cm−1 [18] was selected. The choice of the P(1) line was based on practical criteria, since this line was easily accessible with the used QCL, first, and second, because this line has some intermediate line strength of 9.072·10−20 cm−1/(molecule·cm−2) [18] providing neither the ultimate limit in terms of sensitivity nor the smallest absorption feature. A reference gas mixture containing a nominal level of 100 μmol/mol CO, 5% CO2, and 15% O2 in nitrogen, gravimetrically prepared by the Instituto Portugues da Qualidade (IPQ) [25], was filled into the gas cell for analysis. The spectrometric measurements were performed at a total gas pressure of 407.1 hPa referring to [14]. Prior to the measurements, effects such as those of the “rapid passage” [26], a non constant chirp rate along the QCL pulse, the dependence of the instantaneous line width of the QCL on the chirp rate [27] and the nonlinearity of the TE-cooled HgCdZnTe detectors were checked and, if necessary, accounted for as pointed out in [9].

Figure 2(a) shows typical signals from the intrapulse mode QCLAS setup. The mid-IR QCL radiation was transmitted through a removable etalon with a free spectral range (FSR) of about 0.05 cm−1 placed on the beam path to detector 1. As evidenced by 37 measured etalon fringes depicted in Figure 2(b), the intrapulse chirp for this condition was shown to be 1.85 cm−1. Figure 2(a) represents the signal of the sample channel detector with CO absorption (solid line) and the reference channel detector (dashed line) with the etalon removed, respectively. One hundred QCL pulses were averaged and measured in the time domain in units of sample points (SP).

The data in Figure 2(a) were converted to absorbances as depicted in Figure 2(d) using the detector signals of the two channels, representing a simultaneous two channel measurement scheme. The wave number axis in Figure 2(d) was accomplished by means of the chirp rate [9]. is the known free spectral range of the etalon in wave numbers and FSRSP is the experimental fringe separation measured in the time domain (SP) which is visible in panel (c) of Figure 2. was determined by calculation, using the refractive index of the etalon material and its length; FSRSP by fitting the measured etalon transmission spectrum with a multipeak function.

Fitting a Voigt profile to the measured absorbance data in Figure 2(d) by means of a nonlinear least square fit [28] (Levenberg-Marquardt [29]) delivered the absorption peak area. The uncertainties in the absorption peak area from fitting ranged from 0.5 to 1% relative.

Figure 3 depicts spectrometric CO amount fractions, derived directly using (4), as a function of their sequence number. From them, the mean CO amount fraction was computed to be 93.7 μmol/mol. The expanded uncertainty of each measurement evaluated using the GUM workbench [30] is in the ±9 μmol/mol range. As presented in Figure 3, the directly retrieved CO amount fractions agree with the gravimetric reference value of () μmol/mol certified by [25].

4. Cw-QCLAS

The intrapulse QCLAS measurements were performed with a path length of 82 cm. Because of the moderate line strength value of the probed P(1) line at 2139.43 cm−1 [18], we could only derive measurable absorbance values (Figure 2) at 407.1 hPa for the 100 μmol/mol CO gas mixture. However, for these conditions the signal-to-noise ratio of 3 is poor as seen in Figure 2(d). Therefore, to enhance the sensitivity, we turned to a laser spectrometer comprising a continuous wave (cw) mode QCL (Hamamatsu). By means of the cw-QCL which was mounted in a high-heat load package, the R(12) line of CO at 2190.02 cm−1 (see Figure 4) was probed with a line strength about three-times stronger than that of the P(1) [18]. To slightly compensate for the larger line strength we used a single pass gas cell with reduced optical pass length of L = 46.1 cm to perform the spectrometric measurements.

The cw-QCL package was mounted on a water-cooled heat sink and an AR-coated aspheric lens (Lightpath) was used to collimate the highly diverging beam. The experimental setup is shown in Figure 5. The temperature of the cw-QCL was kept constant using a Lightwave LDT-5545B temperature controller (long-term stability of 0.01°C) and the laser current was regulated by a Lightwave LDX-3232 high compliance laser driver.

Tuning of the laser wavelength was accomplished by applying a voltage ramp to the external modulation input of the laser driver. In order to make use of the low-noise AC PbSe detector the laser beam was chopped (amplitude modulated) and the detector signal subsequently demodulated by means of a lock-in amplifier before feeding into the A/D conversion.

For the CO quantification, a similar gravimetrically prepared gas mixture prepared by IPQ [25] of nominal 100 μmol/mol CO, 150 mmol/mol O2, 50 mmol/mol CO2, and 3 μmol/mol CH4 in N2 balance was measured. The measurements were performed referring to the TILSAM method [14]. While probing the R(12) line, interferences from nearby lines of H2O, CO2, and CH4, which were at least 104-times weaker, were negligible compared to the CO line strength of = 2.876·10−19 cm−1/(molecule·cm−2) [18].

To derive the absorbance from the measured detector signal, first a signal was measured with a nonabsorbing gas (nitrogen) in the cell, after which, the gas mixture, subsequently filled into the gas cell, was measured. The absorbance was then evaluated using the ratio of these two signals.

Spectral data measured at atmospheric pressure are shown in Figure 6(a). By comparing the measured absorbance to a simulation based on HITRAN [18], displayed in Figure 6(b), a fairly good agreement is observed. A Lorentz function is fitted to the experimental data to derive the line area. The respective CO amount fraction result of = 96.84 μmol/mol calculated by means of (4), using the derived line area and the HITRAN line strength figure (stated uncertainty of 2–5%) [18] is in good agreement with the nominal reference of 100 μmol/mol. The uncertainty of a single measurement is in the ±5 μmol/mol range. The latter is pretty much limited by the availability of the line strength uncertainty.

Besides this quantitative agreement, the most striking feature of the experiments described in the last two sections is the much better raw data quality of the cw-QCLAS instrumentation. This is to some smaller extent due to the somewhat larger “absorption signal strength,” that is, peak absorbance of 0.18 and 0.07 for cw- and intrapulse QCLAS, respectively, and to a much larger extent due to the much smaller noise level of our cw-QCL experiment.

5. Cavity Ring-Down Spectroscopy (CRDS)

So far, the measurements above demonstrate the feasibility of pulsed and cw-QCLAS for absolute CO amount fraction quantification based on the TILSAM method on a 100 μmol/mol level of CO. However, to reach a sensitivity of 0.8 μmol·mol−1, equivalent to the CO amount fraction in exhaled human breath, one has to enhance the sensitivity of the QCLAS instruments. This could have been done by increasing the absorption path length. However, we decided to change to another spectroscopy technique. Cavity ring-down spectroscopy has been proven to provide a very versatile tool in trace gas detection as well as in applications like breath analysis or process controlling [2, 31, 32]. First of all, we have compared QCLAS as demonstrated above with cw-CRDS on the 100 μmol/mol CO level. Then, we were applying QCL-based cw-CRDS to much lower amount fraction levels to demonstrate sensitivity and accuracy at relevant breath levels.

The theory of CRDS employed to the measurements shown in this section has been published in some very prominent studies [11, 3336]. Briefly, there is not that much of a difference to the final equations presented in Section 2 of this work. However, compared to (4) (TDLAS and QCLAS), in CRDS is given in units of cm−2 [4] and the path length has to be omitted in (3) and (4) for CRDS.

Figure 7 depicts a schematic of a typical CRDS spectrometer. In contrast to QCLAS, CRDS is based on the measurement of light absorption inside a high finesse cavity. Once some reasonable light intensity is builtup inside the cavity the incoupling is interrupted, for example, by means of an acousto-optical modulator (AOM). The decay of intracavity light intensity is measured, and referred to as ring-down signal, by a detector at the outcoupling side of the cavity. The out-coupling of light is due to the finite reflectivity (less than 1) of the cavity mirrors. The temporal decay constant of the ring-down signal is derived and referred to as ring-down time . Ring-down times without absorbing gas inside the cavity, and with gas inside the cavity , are then compared to derive the absorption spectrum. In case absorption of light occurred, will be measurably larger than [37]. Due to the trapping of the light inside the cavity the effective interaction length with gas molecules will extensively exceed the physical length of the cavity. Consequently, compared to single pass QCLAS at least, CRDS yields a higher sensitivity for given physical footprints of gas cells, depending of course on the mirror reflectivity.

For the present study we have employed two different setups based on cavities, both with mirror reflectivities of 0.9998 in the 4.6 μm range. With this, effective path lengths of few kilometres were realized. The first setup, referred to as the CALOS spectrometer, is based on a continuous wave CO sideband laser as was presented in [7, 36]. Its performance in CO detection in breath analysis application was demonstrated in [12]. The high finesse cavity is used to confine the absorbing gas sample while continuously flowing through the cavity ring-down cell at constant and controlled flow rates of about 1000 cm3/s. With typical ring-down times of 30 μs its path length equivalent results roughly to about 10 km.

For quantitative CRDS-based detection of CO, the CO sideband laser was tuned to probe the 13C16O P(25) line at 1994.7 cm−1 [18]. The measurements were also performed in accordance with [14]. Figure 8(a) shows cw-CRDS data obtained by the CALOS setup from a single ring-down trace (a), where the ring-down time is fitted on, via processed absorption coefficient data of the probed line for a certain value of (b), to the resultant regression-based retrieval [14] on a set of values (c). Here, summarizes the experimental parameters, . The slope of the linear dependence versus , derived by means of a generalized regression analysis is yielding the amount fraction result.

The ring-down times, measured before and after the CO gas mixture was injected into the cavity, were used to derive the absorption coefficients at respective wave numbers. The measurements were performed at total gas pressures between 30 and 70 hPa and a gas temperature of 292 K. By fitting the absorption coefficient with a Voigt profile, the area of the absorption line was determined. As evidenced by typical fingerprints in the residuals, Figure 8(c), the Voigt profile does not completely describe the actual experimental line shape. This might be caused by some additional effects like the Dicke-narrowing which is not taken into account applying a Voigt profile. However, by means of a regression-based retrieval [9], depicted in Figure 8(d), the determined CO amount fraction evaluated from these data was () μmol·mol−1. The reproducibility was estimated to be ±3% based on the standard deviation of individual measurements. The direct retrieval yields a mean amount fraction of () μmol·mol−1 as averaged on the different experimental runs. The 5.2 μmol/mol uncertainty figure derived for this direct retrieval is combined from the uncertainties of each individual result. Given the comparably larger uncertainty figure of the regression-based retrieval both results are agreeing with each other. Just to note, of course, the two retrievals are not independent, since they have been performed on the same identical data set. However, the regression-based result is based on a linear model whereas the direct retrieval is the mean of the different experimental runs. For the linear model approach a generalized linear regression has been performed by means of the BLeast software [38].

In order to go down to a few ppb of CO, a second, cw-QCL-based CRDS setup was used to perform low-concentration CO measurements. Details on the principle design and operation of this type of CRDS setup are given in [32]. Low CO amount fraction gas mixture standards were prepared dynamically by VSL by means of dilution of reference gas mixtures [39].

For the cw-QCL CRDS-based CO quantification we chose again the R(12) line at 2190 cm−1. Measurements on mixtures of CO in nitrogen were complemented by data on mixtures of CO in hydrogen. The aim of these efforts was to demonstrate the quality and sensitivity of the CO quantification applied in this work down to breath-relevant amount fraction levels. Figure 9(a) depicts the two extremes of the CO-in-nitrogen-mixture data, spanning the range from 100 nmol/mol CO down to 10 nmol/mol. The lines were measured at () K and () hPa, and subsequently, fitted by Lorentz profiles. The linearity of the applied quantification can be derived from Figure 9(b) where results on CO/N2 and CO/H2 mixtures are displayed, spanning the range from 200 nmol/mol down to 2.4 nmol/mol of CO.

The reproducibility of the CO quantification achieved with cw-QCL CRDS was 0.6% based on the derived standard deviation for the range of 200 nmol/mol to 10 nmol/mol CO.

6. Discussions

The spectrometric CO amount fraction results from three different spectrometers agree with the respective reference values at the 100 μmol/mol CO level. The relative standard uncertainties of these spectrometric CO amount fraction results are in the few percent range. Comparability was achieved among the direct absorption and the cavity ring-down results as summarized in Table 1 for direct retrievals.

For each of the three results the value component of the degree of equivalence calculated as the difference of the respective result and the gravimetric reference value is smaller than its uncertainty component U(D).

Although direct absorption spectroscopy, as applied in the present study by means of intrapulse and cw-QCLAS, could be improved in terms of sensitivity by increasing the optical path length, a much higher potential in breath analysis has to be attributed to the cavity ring-down spectroscopy for our purposes. This is demonstrated by the CALOS and the cw-QCL CRDS setups used in this work. However, to keep on going with direct absorption, it is clearly shown that based on the technology employed for this work, cw-QCLs are superior to pulsed systems in terms of noise and in turn on sensitivity for a given path length. With respect to the spectrometer’s path length, one could argue that an increase would be possible, thus increasing our sensitivity towards that of our CRDS. In fact, we might have gained up to two orders of magnitude, since multi pass Herriott-type cells with path lengths of 100 or even 200 meters are available. Instead, the path length equivalent of CRDS systems is easily on the kilometer scale.

Following another spectrometer comparing study [17], Table 2 compares the instrument’s sensitivity of a previously reported cw-QCL cavity ring-down spectrometer with those of our two CRDS setups to provide some performance characteristics. Numbers compared in Table 2 were computed referring to [17, 40], with quantities and being the average ring-down time and its standard deviation, respectively, for an empty cavity or a nonabsorbing gas. By means of these quantities the absorption coefficient can be expressed as With this absorption coefficient an instrument’s base line sensitivity can be defined by where is the acquisition rate of succeeding ring-down events. Although not squeezed to its outmost performance, our cw-QCL-based CRDS setup has some comparable performance ( sensitivity level) to that of Brumfield et al. which was reported in [17], being higher than those of previous studies [40, 41]. However, the CALOS system used in this work has some similar instrument sensitivity to that of [40] being both in the range.

Nevertheless, sensitivities of both instruments used in this work were found to be very sufficient for the CO amount fraction measurements as presented above. Although not driving at impressive sensitivity figure of merits in this work, it is maybe worthwhile to have such information which details somehow the quality of used instrumentation.

The ideas of the TILSAM method [14] were applied in this work to two different spectroscopic techniques. Its potential relevance to breath analysis is evidenced by the present study with respect to the achieved sensitivity down to a few ppb, shown for the cw-QCL CRDS setup. With this sensitivity, CO levels typically to be measured in applied breath analysis could be matched, while quality flags of the results as delivered by GUM [15] compliant uncertainty statements are still provided. The uncertainty for the cw-QCL CRDS is mostly limited by the availability of line strength data with uncertainties better than 5%, as, for example, reported by HITRAN 2008 [18] on the probed R(12) line of CO at 2190 cm−1. At the 100 μmol/mol level a reproducibility of 0.5%, relative, was achieved. The relative deviation of the spectrometric results from the gravimetric reference was determined to be in the range of 1% to 0.3% on the 200 to 10 nmol/mol level, what would totally be covered by a line strength uncertainty figure of 2%.

With a documented method, for example, that of TILSAM, spectrometric amount fraction results could be made comparable. Ensuring the comparability of measurement results is crucial, for example, in medical treatment or environmental measurements for decision-making procedures. The comparability of measurement results is ensured if the results are traceable. Traceability is also in line with the stated aims of the international quality assurance standard ISO EN 17025 [42].

7. Conclusions

The goals of the presented study were threefold, namely, first, to demonstrate comparability of CO quantification by means of absolute laser spectroscopy referring to the TILSAM method. This goal was achieved at the 100 μmol/mol CO level for three different experimental techniques and setups used. The overall achieved level of comparability can be expressed by means of a maximum (see Table 1) as to be within ±6%. Second, to demonstrate CO quantification concepts based on metrological aspects such as uncertainty assessments, as it may provide impact to practical breath analysis. This was achieved best by means of the estimated uncertainty levels of the cw-QCL results at some few percent, mostly limited by the available line strength data. Finally, the third goal was to demonstrate CO quantification capabilities down to a few nmol/mol level as being relevant for the breath matrix.

Acknowledgments

Parts of this paper were performed in an EURAMET joint research project and received funding from the European Union Seventh Framework Programme, ERA-NET Plus, under the iMERA-Plus Project—Grant Agreement no. 217257. J. A. Nwaboh acknowledges also the support from the Braunschweig International Graduate School of Metrology (http://www.igsm.tu-bs.de).