Detecting parts-per-billion carbon monoxide with an ultra-enhanced near-infrared photoacoustic sensor
Yaopeng Cheng, Ting Chen, Ruili Zhang, Sailing He

TL;DR
A new near-infrared sensor detects carbon monoxide at very low concentrations, offering a cost-effective alternative to more expensive sensors.
Contribution
The integration of three enhancement techniques enables ultra-sensitive carbon monoxide detection at parts-per-billion levels.
Findings
The sensor achieved a minimum detection limit of 190 ppb at 10 seconds.
The detection limit can be improved to 11.4 ppb based on Allan analysis.
The sensor is cost-effective, with a price less than one-third of mid-infrared photoacoustic sensors.
Abstract
An ultra-enhanced near-infrared (NIR) photoacoustic gas sensor was developed by integrating three enhancing techniques: (a) boosting the excitation power up to 2 W via a custom-built large-mode erbium doped fiber amplifier (EDFA), (b) exploiting the acoustic resonance amplification of a hyperbolic nonlinear resonator (HNR), and (c) increasing the effective absorption path length by using a near-concentric multipass cavity (MPC) with 20 reflections. A weak CO absorption line at 1566.64 nm with the intensity of 2.074 × 10−23 cm/molecule was selected. The photoacoustic signal was enhanced 396 times. A minimum detection limit (MDL) of 190 ppb at 10 s was achieved and can be improved to be 11.4 ppb according to the Allan analysis, which was comparable to a mid-infrared (MIR) photoacoustic sensor. The ultra-enhanced NIR photoacoustic sensor is a cost-effective solution for the ppb-level trace…
Genes, proteins, chemicals, diseases, species, mutations and cell lines named across the full text — each resolved to its canonical identifier and authoritative record.
Click any figure to enlarge with its caption.
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Figure 7
Figure 8
Figure 9
Figure 10
Figure 11Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsSpectroscopy and Laser Applications · Photoacoustic and Ultrasonic Imaging · Combustion and flame dynamics
Introduction
1
Carbon monoxide (CO) is a by-product of incomplete combustion of carbon-containing fuels, such as gasoline, natural gas, and coal that are widely used for petrochemical refining, transportation, power generation and other applications [1]. CO is colorless and odorless but highly toxic to humans. Exposure to high concentrations of CO can result in a range of adverse effects, including cardiac and pulmonary disorders, central nervous system involvement, fatigue, coma, and even death [2]. Furthermore, in atmospheric chemistry, CO is recognized as a major air pollutant due to its active role in tropospheric ozone formation and global warming [3], [4]. Therefore, the sensitive detection of CO is highly desired in monitoring emissions from carbon-containing fuels.
Photoacoustic spectroscopy (PAS) has been widely used in trace gas detection due to its high sensitivity and high selectivity [5], [6], [7]. Traditionally, PAS relied on microphones for the detection of acoustic waves. However, the quartz tuning fork (QTF) has recently emerged as a highly sensitive alternative. Consequently, the field has evolved from the use of conventional microphones to commercial QTFs [8], and more recently, to novel custom-designed QTFs [9] coupled with various quartz-enhancement techniques [10], [11], [12]. These advancements have substantially elevated the detection performance of the PAS-based sensors. The PAS-based CO sensors demonstrated minimum detectable limits (MDL) down to the ppb level using the quantum cascade lasers (QCLs) excite the strong fundamental rotation-vibration band (4.6 μm) in the mid-infrared (MIR) of CO [13], [14], [15], [16], [17], [18]. However, the price of QCLs is several thousand dollars [19]. Also, QCL has a large threshold current, and thus a high power consumption in the order of few watt [20]. Therefore, the operation of QCL requires an additional active air or water cooling to manage the substantial heat generated by the power consumption [21], which increases the overall size, weight, and power (SWaP) of the sensor.
Near-infrared (NIR) diode lasers are attractive alternative to the QCLs in the gas sensing application. Typical emission wavelengths of commercial NIR diode lasers are in the spectral range of 0.5–2 μm which covers the gas vibrational overtones spectra [22]. NIR diode lasers have long service life (> 10 years) and low cost [23]. They can be easily operated at room temperature without any active cooling, simplifying the entire system and reducing the SWaP of the sensor. This makes them more suited for practical, field-deployable applications. Their narrow linewidths (<1 MHz), wide tuning ranges and great modulation capability are ideal for high selective and sensitive trace gas detection. Typically, they come with pigtailed single-mode fibers, which facilitate the beam delivery and enable efficient beam coupling into different gas sensor configurations [24]. The output power of the NIR diode lasers is on the order of tens of milliwatt, but can be readily boosted up to the watt level by using telecommunication optical fiber amplifiers. However, the gas absorption line intensity in the NIR region is typically 2–4 orders of magnitude lower than those in the MIR region. Fortunately, various techniques for enhancing photoacoustic signals have been extensively developed and reported [23]. By strategically exploiting and integrating these techniques, the limitation of weak gas absorption in the NIR region can be effectively compensated.
In this paper, a NIR diode laser was employed for the trace CO detection. Three techniques to enhance the photoacoustic (PA) signal were adopted: (a) boosting excitation power via a custom-built large-mode erbium doped fiber amplifier (EDFA), (b) exploiting acoustic resonance amplification by a hyperbolic nonlinear resonator (HNR) with a high setup constant, and (c) increasing the effective absorption path length by using a near-concentric multi-pass cavity (MPC). The cost of developing an ultra-enhanced NIR PA sensor was less than one-third that of a MIR PA sensor. In trace CO detection experiments, an MDL of 190 ppb at 10 s was achieved and can be improved to be 11.4 ppb according to the Allan analysis.
Materials and methods
2
Selection of the CO absorption line
2.1
As shown in Fig. 1(a), in the infrared region, the CO molecule has three main absorption bands: a. the second overtone band (1.56 μm), b. the first overtone band (2.33 μm), and c. the fundamental band (4.6 μm). In terms of the light source, NIR diode lasers are preferred over QCLs. Therefore, the two NIR absorption bands of CO (1.56 μm and 2.33 μm) are considered. The second overtone band (1.56 μm) of CO is two orders of magnitude weaker than the first overtone band (2.33 μm). However, it located at the wavelength window for optical fiber communication, so an optical fiber amplifier can be employed to boost the optical power. Hence, the second overtone band (1.56 μm) is a chosen to implement the NIR trace CO detection. The potential spectral interference from atmospheric water vapor (H_2_O) and carbon dioxide (CO_2_) was evaluated. As shown in Fig. 1(b), absorption coefficients of 500 ppm CO, 500 ppm CO_2_ and 1 % H_2_O were simulated based on the HITRAN 2020 database [25]. Accordingly, the CO absorption line at 1566.64 nm with the line intensity of 2.074 × 10^−23^ cm/molecule was selected in this work.Fig. 1. Selection of CO absorption line. (a) Three main infrared absorption bands of CO. (b) Absorption coefficients of 500 ppm CO, 500 ppm CO_2_ and 1 % H_2_O from 1562 to 1574 nm. All absorption coefficients are simulated at 1 atm and 296 K based on the HITRAN database.Fig. 1
Techniques to enhance the photoacoustic signal
2.2
For resonant operation, the laser modulation frequency is tuned to one of the resonance frequencies of the PA cell. The PA signal measured by the microphone is given by [23]:
where is the adiabatic coefficient of the gas, is the absorption path length, is the quality factor of the PA cell, is the resonance frequency of the PA cell, is the volume of the PA cell, is the laser power, is the gas absorption coefficient, and is the microphone sensitivity.
Eq. (1) indicates the three main techniques to enhance the PA signal: (a) boosting the laser power by using an optical amplifier; (b) exploiting the acoustic resonance amplification by employing a resonant PA cell with a high quality factor ; (c) increasing the effective absorption path length by introducing an MPC that allows the laser beam to reflect multiple times within the PA cell.
Design of the HNR and the near-concentric MPC
2.3
The nonlinear resonator, an axisymmetric closed tube with a varying cross section, exhibits better acoustic resonance effects than a traditional cylindrical resonator [26], [27], [28]. Therefore, a hyperbolic nonlinear resonator (HNR) was designed and used as a resonant PA cell to enhance the PA signal through acoustic resonance amplification. The sketch of the HNR is shown in Fig. 2(a), mainly including a hyperbolic resonator tube (colored in orange), two buffer chambers (colored in gray), gas inlet and outlet (colored in green), and four microphones. As shown in Fig. 2(b), finite element analysis (FEA) was performed based on the harmonic full linearized Navier–Stokes (FLNS) equations to simulate the pressure distribution in the HNR [29], [30]. A Gaussian distributed heat source with an amplitude of 0.35 W/m^3^ was introduced along the axis of symmetry of the HNR to simulate the gas absorption by a single-pass modulated laser beam. Through simulation, the Lres, Rres, δ, and Lbuf of the HNR were optimized to 40 mm, 3 mm, 120 m^−1^, and 12 mm, respectively. The resultant HNR obtained a strong of 1.83 × 10^−5^ Pa, leading to a high setup constant of 3420 Pa·cm/W.Fig. 2(a) A sketch of our hyperbolic nonlinear resonator. (b) Geometric features of the cavity and simulated distribution of the absolute pressure at resonance.Fig. 2
A near-concentric MPC was designed to further enhance the PA signal. The near-concentric MPC consisted of two silver-coated concave mirrors with a focal length of 25 mm and a diameter of 25.4 mm. The reflectance of the mirrors at a wavelength of 1566 nm was above 97.3 %. The two concave mirrors were coaxially positioned 100 mm apart to satisfy the near-concentric condition. The number of reflections was determined by the distance of incidence, the angle of incidence, and the tilt angle of the mirror. We simulated the optical path of the near-concentric MPC, with comprehensive details available in Appendix A. By adjusting these three parameters to 0.6 mm, 13.48°, and 0.19°, respectively, up to 20 reflections were achieved. The optical path and irradiance maps of the near-concentric MPC are shown in Fig. 3(a) and (b) respectively. The laser beam entered the MPC, reflected 10 times on each mirror, and then exited the MPC after 21 passes.Fig. 3. Simulation of the near-concentric multipass cavity. (a) Ray tracing. Irradiance maps of (b) the left concave mirror, (c) the center of the multipass cavity, and (d) the right concave mirror.Fig. 3
Design of the large-mode erbium-doped fiber amplifier
2.4
A custom-built large-mode erbium-doped fiber amplifier (EDFA) is configured as Fig. 4. The amplifier system consisted of a pre-amplifier followed by a large-mode fiber amplifier. The optical isolators (ISOs) were employed to suppress back-reflection and ensure unidirectional operation, while wavelength division multiplexers (WDMs) employed for pump-signal combination. In the pre-amplifier, a ∼ 60 cm segment of commercial single-cladding erbium-doped fiber (LIEKKI, Er-110–4–125) was used as the gain medium. The EDF was bidirectionally pumped by a commercial 980 nm laser with a maximum output power of 600 mW. A polarization controller (PC) was placed before the large-mode fiber amplifier to optimize the signal polarization state. In the large-mode fiber amplifier stage, a ∼ 3.5 m segment of commercially available large-mode-area Er/Yb co-doped double-cladding fiber (EYDCF, IXF-2CF-EY-O-25–250-COM) was used as the gain medium. This fiber features a core diameter of 25 μm, which significantly reduces the accumulated nonlinear effects and thereby facilitates high-power laser generation [31]. A commercial WDM was utilized to simultaneously couple a 915 nm multimode laser diode pump into the cladding and the energized the seed laser into the core of the EYDCF. The 915 nm laser provided an adjustable pump power ranging from 0 to 25 W.Fig. 4. Configuration of the large-mode erbium-doped fiber amplifier. ISO: isolator; WDM: wavelength division multiplexer; EDF: erbium-doped fiber. EYDCF: large-mode-area Er/Yb co-doped double-cladding fiber.Fig. 4
Experimental section
3
System configuration and sample gas preparation
3.1
A schematic configuration of the NIR CO sensor is depicted in Fig. 5. In our sensor, a 1566 nm distributed feedback laser (DFBL) with an optical power of 50 mW (OptoChip Optoelectronics, China, Model OC45650SMF) was chosen as the optical excitation source and was subsequently amplified by the custom-built large-mode EDFA to 2 W. The output laser beam was collimated by a fiber collimator (Thorlabs, USA, Model F230FC-1550) and coupled into the near-concentric MPC via a high-precision XYR translation stage. One mirror was mounted on a manual kinematic mount (Thorlabs, USA, Model KS1), while the second mirror was installed on a 2-axis piezo-driven mount (Thorlabs, USA, Model PIM1/M). The tilt angle of the mirror was finely controlled by a piezo inertia motor controller (Thorlabs, USA, Model KIM101), allowing for precise alignment to achieve 20 reflections. Wavelength modulation spectroscopy with second harmonic detection (WMS-2f) was adopted. A modulation signal composed of a sawtooth wave of 0.1 Hz and a sine wave of 3.232 kHz (half of the resonance frequency of the HNR) was sent to the laser controller to drive the DFBL. A temperature controller (Thorlabs, USA, Model TED200C) and a current driver (Thorlabs, USA, Model LDC205C) were employed to stabilize the laser temperature and convert the modulation signal into the laser driving current. To improve the acoustic detection performance, we used four microphones to detect the PA signal. The PA signal scales with the number of the microphones, while the uncorrelated noise increases only by a factor of [32]. Thus, the signal-to-noise ratio improves with . A homemade printed circuit board (PCB) amplified and summed the four microphone signals. The summed signal was connected to the lock-in amplifier (LIA, AMETEK, USA, Model 7270) to demodulate the second harmonic component (2 f). The equivalent noise bandwidth (ENBW) of the LIA was set to 0.33 Hz to reduce noise. A data acquisition card (DAQ, National Instruments, USA, Model USB-6229) along with a custom-written LabVIEW program were utilized to provide the modulation signal, generate the reference signal, collect the 2 f signal, and display the measured results. All the experiments were performed in the air-conditioned laboratory with the room temperature of 23℃. The gas preparation system consisted of two standard gas tanks of 60 ppm and 3000 ppm CO diluted in N_2_, a gas tank of pure N_2_, and two high-precision mass flow controllers (MFCs). Various CO/N_2_ gas mixtures of defined concentrations were prepared for detection experiments by means of MFCs. The standard gas sample of CO was flow-control diluted by the pure N_2_ via the two MFCs and the CO concentrations of gas mixtures were varied by changing the ratio of the flow rates of the two MFCs.Fig. 5. Simplified block diagram of the ultra-enhanced NIR PAS-based CO sensor. EC mic: electret condenser microphone; PCB: printed circuit board; 2 f: second harmonic signal; DAQ: data acquisition; EDFA: erbium doped fiber amplifier; MFC: mass flow meter.Fig. 5
Characterization of the large-mode erbium-doped fiber amplifier
3.2
The performance of the custom-built large-mode EDFA was systematically characterized. As illustrated in Fig. 6(a), the output power at 1566 nm as a function of the 915 nm pump power was investigated using an optical power and energy meter consoles (Thorlabs, USA, Model PM100D2). The amplifier exhibited a highly linear input-output relationship. The seed laser, with a power of 50 mW was first amplified to 162 mW by the pre-amplifier. Then the output power scaled efficiently from 162 mW to over 3 W by the large-mode fiber amplifier as the pump power increased from 1 W to 16 W. This linear trend indicated an excellent pump conversion efficiency and the absence of significant parasitic lasing or thermal roll-over within this operating range. An output of 2 W, generated at a pump power of 11 W, was employed for the subsequent gas sensing experiments.Fig. 6. Performance characterization of the custom-built large-mode EDFA. (a) Output power at 1566 nm as a function of pump power. (b) Output spectra measured at different amplification stages.Fig. 6
The spectral evolution through the amplifier stages was recorded by an optical spectrum analyzer (Thorlabs, USA, Model OSA203C). Fig. 6(b) presents the optical spectra of the seed laser, after the pre-amplifier, and after the large-mode fiber amplifier. The center wavelength remained stable throughout the amplification stages, with no observable shift or broadening within the instrument resolution of 0.06 nm. Furthermore, the side-mode suppression ratio (SMSR) was maintained at a high level throughout the system, showing a slight increase from 47.0 dB for the seed laser to 49.5 dB after the pre-amplifier, and reaching 50.3 dB after the large-mode fiber amplifier. No severe spectral distortion or emergence of new spectral components was observed. The results demonstrated the amplifier's superb spectral stability and its linear operating regime, effectively suppressing excessive amplified spontaneous emission (ASE) and nonlinear spectral broadening.
Modulation optimization for the HNR
3.3
The laser modulation depth has a significant impact on the 2 f signal. Hence, the dependence of the 2 f signal amplitude on the modulation depth was investigated. Experiments were conducted with a certified CO/N_2_ concentration of 3000 ppm. The modulation frequency was set to half of the resonance frequency of the HNR (3.232 kHz). The peak values of the 2 f signals were recorded as modulation depths varied from 0.02 to 0.09 nm with an increment of 0.01 nm. As shown in Fig. 7, with increasing modulation depths, the amplitude of the 2 f signal first increased, peaked at the modulation depth of 0.07 nm, and then decreased. Two representative PA 2 f curves at the modulation depths of 0.02 and 0.07 nm are illustrated in the inset graph of Fig. 7. Accordingly, the modulation depth of 0.07 nm was selected in this work to achieve the strongest 2 f signal amplitude.Fig. 7. Photoacoustic 2 f signal peak amplitude as a function of the modulation depth for 3000 ppm CO/N_2_ at 1 atm and 296 K. Inset: representative PA 2 f curves at the modulation depths of 0.02 and 0.07 nm. The amplitude was normalized to one.Fig. 7
Enhancement of the photoacoustic signal
3.4
We tested the effectiveness of the three adopted signal-enhancing techniques. Experiments were conducted under the same experimental conditions for fair comparisons. The 2 f signal curve was recorded in the standard 3000 ppm CO/N_2_ with three different system configurations: (a) the single-pass HNR; (b) the multi-pass HNR; (c) the multi-pass HNR with the power boosted by the EDFA. The results are shown in Fig. 8. The 2 f signal peak amplitude of the multi-pass HNR was 2.88 mV, which was approximately 11 times stronger than that of the single-pass HNR (0.248 mV). As for the enhancement by EDFA, the 2 f signal peak amplitude of the multi-pass HNR with the power boosted was up to 103 mV, resulting in a 36-fold signal enhancement than that of the multi-pass HNR (2.88 mV). By adopting these techniques, the PA signal was enhanced by a factor of 396 altogether; the weak CO gas absorption line strength in the NIR region was compensated. An increase in baseline level, from 0.042 mV to 0.12 mV, was observed. This rise is due to the use of the near-concentric MPC and the EDFA, which can contribute to a higher background PA signal. The enhancement by the near-concentric MPC can be further improved by increasing the number of reflections, the reflectance of the concave mirrors, and the transmission of the optical windows of the HNR, while the enhancement by the EDFA can be improved by simply increasing its output power.Fig. 8. Photoacoustic 2 f curves for 3000 ppm CO/N_2_ at 1 atm and 296 K measured by three different configurations: (a) bare hyperbolic nonlinear resonator (HNR), (b) HNR with the near-concentric multi-pass cavity (MPC), and (c) HNR with the near-concentric MPC and custom-built large-mode EDFA.Fig. 8
Photoacoustic detection of CO diluted in N2
3.5
The concentration response of the developed PA sensor was examined through the calibration. A series of CO/N_2_ gas mixtures of known concentrations from 10 to 3000 ppm were successively introduced into the PA cell for the sensor calibration. Fig. 9(a) depicts the measured PA 2 f signal curves at different concentrations of CO. The corresponding peak values of the 2 f signals are plotted in Fig. 9(b) as a function of CO concentration. According to the linear fit, the sensor showed an excellent linear response to CO concentrations in the range of 10–3000 ppm with an R-square value of 0.99962 and a calibration coefficient of 0.0338 mV/ppm.Fig. 9. The calibration of the ultra-enhanced NIR PA sensor for the detection of CO at 1 atm and 296 K. (a) Measured 2 f signal curves at different CO concentrations. (b) Photoacoustic 2 f signal peak amplitude as function of CO concentration.Fig. 9
As a common tool to demonstrate the stability performance of the sensor, Allan–Werle deviation analysis was performed to further investigate the long-term stability and optimal achievable MDL of the developed NIR PA sensor. The pure N_2_ was continuously flushed into the PA cell at a flow rate of 300 mL/min, and the 2 f noise was collected every 10 s without averaging for approximately two continuous hours. As shown in the inset of Fig. 10, with the integration time of 10 s, the MDL was tested to be 190 ppb. The Allan–Werle deviation analysis was then performed on the 2 f noise in N_2_ and the results are shown in Fig. 10. The deviations in millivolt units were converted into concentration units (ppb) using the calibration coefficient of the sensor. As the integration time was extended and the white noise was suppressed, the Allan–Werle deviations decreased as the function of 1/t^0.56^. The MDL can be improved to 11.4 ppb with an integration time of 1270 s. According to the Allan–Werle analysis, the developed PA sensor exhibited excellent MDL performance and long-term stability.Fig. 10. The Allan–Werle deviation analysis. The dashed line represents the 1/t^0.56^ slope. Inset: the raw data of the PA 2 f noise measured in pure nitrogen for 2 h.Fig. 10
Table 1 shows the comparison between the developed ultra-enhanced NIR PA CO sensor and previously reported PAS-based CO sensors. In the NIR region, our ultra-enhanced NIR PA sensor achieved the best MDL of 11.4 ppb, which was comparable to that obtained by MIR PA sensors. However, the cost of the ultra-enhanced NIR PA sensor was merely one-third that of the MIR PA sensors. In the NIR region, reference [33], [34] also utilized the EDFA to boost the laser power. However, they all were based on a traditional H-type PA cell and operated without a multi-pass cavity and they obtained less competitive MDLs than our ultra-enhanced NIR PA sensor.Table 1. Comparison of the developed ultra-enhanced NIR PA CO sensor and previously reported PAS-based CO sensors. MDL: minimum detectable limit. PAS: photoacoustic spectroscopy; HPAS: heterodyne photoacoustic spectroscopy; QEPAS: quartz enhanced photoacoustic spectroscopy; DFBL: distributed feedback laser; QCL: quantum caseate laser.Table 1. BandMethodsTarget linesLaserPowerMDLMIRHPAS4.61 μmQCL21 mW26 ppb [13]PAS4.61 μmQCL21 mW8 ppb [14]QEPAS4*.57 μmQCL32 mW12 ppb [16]QEPAS4.59 μmQCL40 mW22 ppb [18]QEPAS4.65 μmQCL71 mW2 ppb [17]PAS4.60 μmQCL175 mW0.8 ppb [15]NIRQEPAS2.33 μmDFBL3.3 mW1.3 ppm^a^[35]CEPAS2.33 μmDFBLNot given5.1 ppm [36]PAS1.568 μmDFBL1 W0.47 ppm^b^[34]PAS1.566 μmDFBL10 W110 ppb [33]PAS1.566 μmDFBL2 W11.4 ppb [This work]*: amplified by EDFA. a: H_2_O as the promotor; b: sulfur hexafluoride (SF_6_) as the promotor.
Conclusion
4
In conclusion, a NIR diode laser was employed for the trace CO detection. Three techniques to enhance the PA signal were adopted: (a) boosting excitation power up to 2 W via the custom-built large-mode EDFA, (b) exploiting acoustic resonance amplification by the HNR and (c) increasing the effective absorption path length by using the near-concentric MPC with 20 reflections. By adopting these techniques, the PA signal was enhanced by a factor of 396 altogether. In trace CO detection experiments, the ultra-enhanced NIR PA sensor compensated the weak CO absorption line strength in the NIR spectral region. An MDL of 190 ppb at 10 s was achieved and can be improved to be 11.4 ppb according to the Allan analysis, which was comparable to that obtained using an MIR PA sensor. The ultra-enhanced NIR PA sensor was a cost-effective solution for the ppb-level detection of trace gases, offering a price point that is less than one-third that of MIR PA sensors. The present PA gas sensing system can also be utilized for other applications such as fast and accurate measurement of exhaled breath gases for medical applications In medicine, the concentration of CO in exhaled human breath is analyzed as an indicators of health and for the diagnosis of certain ailments, such as asthma, anemia and inflammation [37], [38].
CRediT authorship contribution statement
Ting Chen: Visualization, Software, Investigation, Data curation. Yaopeng Cheng: Writing – original draft, Visualization, Validation, Software, Methodology, Data curation, Conceptualization. Sailing He: Writing – review & editing, Supervision, Resources, Project administration, Funding acquisition, Formal analysis. Ruili Zhang: Supervision, Resources, Investigation, Formal analysis.
Declaration of Competing Interest
We declare that we have no financial and personal relationships with other people or organizations that can inappropriately influence our work. We also declare no conflicts of interests with each other on this study.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Holloway T.Levy H.Kasibhatla I.I., P.Global distribution of carbon monoxide J. Geophys. Res. Atmospheres 105D 1020001212312147
- 2Mattiuzzi C.Lippi G.Worldwide epidemiology of carbon monoxide poisoning Hum. Exp. Toxicol.39420203873923178906210.1177/0960327119891214 · doi ↗ · pubmed ↗
- 3Logan J.A.Prather M.J.Wofsy S.C.Mc Elroy M.B.Tropospheric chemistry: a global perspective J. Geophys. Res. Oceans 86C 8198172107254
- 4Khalil M.A.K.Rasmussen R.A.Carbon monoxide in the Earth's atmosphere: increasing trend Science 2244644198454561778352110.1126/science.224.4644.54 · doi ↗ · pubmed ↗
- 5Wang F.Wu J.Cheng Y.Fu L.Zhang J.Wang Q.Simultaneous detection of greenhouse gases CH 4 and CO 2 based on a dual differential photoacoustic spectroscopy system Opt. Express 3121202333898339133785915910.1364/OE.503454 · doi ↗ · pubmed ↗
- 6Wang F.Fu L.Zhang J.Han Z.Pang S.Xue Q.Cao D.Li Q.Wang Q.TT-type resonator-based differential photoacoustic spectroscopy for trace gas detection Opt. Lett.4982024217321763862110410.1364/OL.520154 · doi ↗ · pubmed ↗
- 7Fu L.Wang F.Zhang J.Han Z.Pang S.Wang Q.Miniaturization study on helmholtz-based photoacoustic cell for PAS-based trace gas detection Measurement 2362024
- 8Ma Y.He Y.Tong Y.Yu X.Tittel F.K.Quartz-tuning-fork enhanced photothermal spectroscopy for ultra-high sensitive trace gas detection Opt. Express 2624201832103321103065067610.1364/OE.26.032103 · doi ↗ · pubmed ↗
