Review Papers

Breath analysis using external cavity diode lasers: a review

[+] Author Affiliations
Ismail Bayrakli

Suleyman Demirel University, Biomedical Engineering, Isparta, Turkey

J. Biomed. Opt. 22(4), 040901 (Apr 17, 2017). doi:10.1117/1.JBO.22.4.040901
History: Received December 17, 2016; Accepted March 28, 2017
Text Size: A A A

Open Access Open Access

Abstract.  Most techniques that are used for diagnosis and therapy of diseases are invasive. Reliable noninvasive methods are always needed for the comfort of patients. Owing to its noninvasiveness, ease of use, and easy repeatability, exhaled breath analysis is a very good candidate for this purpose. Breath analysis can be performed using different techniques, such as gas chromatography mass spectrometry (MS), proton transfer reaction-MS, and selected ion flow tube-MS. However, these devices are bulky and require complicated procedures for sample collection and preconcentration. Therefore, these are not practical for routine applications in hospitals. Laser-based techniques with small size, robustness, low cost, low response time, accuracy, precision, high sensitivity, selectivity, low detection limit, real-time, and point-of-care detection have a great potential for routine use in hospitals. In this review paper, the recent advances in the fields of external cavity lasers and breath analysis for detection of diseases are presented.

Figures in this Article

Technologies for clinical diagnoses and monitoring are always needed for discriminating patients from healthy subjects. To identify diseases and to monitor therapy and/or the progression of diseases, most used methods—such as biopsy, computer tomography, magnetic resonance imaging, endoscopy, and blood analysis—are time-consuming, expensive, invasive, painful, stressful, and difficult to perform. Therefore, noninvasive tools are needed to monitor the diseases in their early stages and to holistically measure the physiological status of patients. There are different potential noninvasive techniques, such as breath, nail, hair, mucus, saliva, sweat, and/or tear analysis. Among them, the detection of volatile organic compounds (VOCs) present in breath at trace concentrations with an acceptable accuracy and precision is most promising. This method may provide noninvasive, point-of-care (PoC), easily repeatable, painless, and real-time diagnosis of a number of disease states and therefore is preferable. Breath biomarkers can be analyzed either after administration of a drug or substrate or without any prior administration of a chemical.

A number of breath molecules have been detected and quantified so far using the different methods, such as proton transfer reaction mass spectrometry (PTR-MS), selected ion flow tube mass spectrometry (SIFT-MS), gas chromatography-mass spectrometry (GC-MS), e-noses, and laser-based sensors. Laser spectroscopic detection techniques—such as tunable diode laser absorption spectroscopy (TDLAS), cavity-enhanced absorption spectroscopy (CEAS), integrated cavity output spectroscopy (ICOS), cavity ringdown spectroscopy (CRDS), photoacoustic spectroscopy (PAS), quartz-enhanced photoacoustic spectroscopy (QEPAS), cavity leak-out spectroscopy, and optical frequency comb cavity-enhanced absorption spectroscopy—provide the features of high-sensitivity, high-selectivity, low costs, real-time measurements, and PoC function. Therefore, these methods are candidates for performing accurate breath diagnostics in hospitals.

To identify and quantify the molecules in breath, narrow linewidth (<1  MHz), mode-hop-free (1 to 2  cm1) wavelength-tunable, single-mode light sources are needed. These features can be obtained using external cavity lasers (ECLs). ECL-based spectroscopic techniques may become very suitable for monitoring diseases by analyzing exhaled breath air samples. Many review papers about breath analysis using laser-based techniques have been published. Ciaffoni et al.1 focused on the current status of volatile sulfur compounds sensing via laser absorption spectroscopy. The review article of Wang and Sahay2 is related to breath analysis of almost all of the biomarkers that have been analyzed in actual human breath using high-sensitivity laser spectroscopic techniques. Wojtas et al.3 presented recent advances in human breath analysis using optical methods. Risby and Tittel4 reviewed the current status of mid-infrared quantum and interband cascade lasers for clinical breath analysis. The review article of Bielecki et al.5 presented an overview on high-sensitivity laser absorption spectroscopic techniques used for trace gas detection. Mid-infrared spectroscopic sensing techniques were reviewed for potential application in breath diagnostics.6 In the review paper,7 techniques such as mass spectrometry-based (GC-MS, PTR-MS, SIFT-MS), laser absorption spectroscopy-based (CRDS and TDLAS), and other spectroscopic techniques for breath analysis applications were compared in terms of advantages/disadvantages, versatilities and plausibility to be transformed in clinical applications. McCurdy et al.8 summarized some of the recent advances in laser absorption spectroscopy based on semiconductor lasers and optical detection techniques for clinically relevant exhaled gas analysis in breath, specifically molecular biomarkers such as nitric oxide, ammonia, carbon monoxide, ethane, carbonyl sulfide, formaldehyde, and acetone.

Regarding the issues above mentioned, the present review article focuses on the latest status of breath analysis using the ECL-based methods. The article is structured as follows: first, we will explain the respiratory system and breath sampling. The absorption of electromagnetic waves by breath molecules is discussed next. The latest developments and functional principle of ECLs are described in detail as well. Finally, the laser-based techniques used for breath analysis and limitations of breath analysis are presented in detail. Regarding this content, the present review article differentiates itself from previous reviews related to the laser-based breath analysis.

The lung is primarily responsible for gas exchange between the atmosphere and the blood. The molecules present in blood can diffuse into the alveoli and subsequently appear in the exhaled breath. Transport of metabolites from the blood to the exhaled breath air can make it possible to utilize breath analysis as a noninvasive diagnostic tool. Low and high blood-soluble gases present in the breath air:9 Gases with low blood-solubility like oxygen (λb:a0.7) and carbon dioxide (λb:a3) exchange in the alveoli. For these gases, there is a simple relationship between the breath and blood, which depends predominately on two factors: the blood solubility of the gas (λb:a) and the distribution of ventilation-to-perfusion in all gas exchange units. Gas exchange of high blood-soluble gases like ethanol (λb:a1756) takes place in the airways, and there is a more complex relationship between end-exhaled breath measurement and blood. Airway gas exchange is related to multiple factors, such as airway temperature, bronchial blood flow, and blood–air partition coefficient. Soluble gases with a blood:air partition coefficient less than 10 (λb:a<10) exchange almost solely with the alveoli, whereas gases with a blood:air partition coefficient greater than 100 (λb:a>100) exchange almost exclusively with the airways. Gases in between these two extremes (10<λb:a<100) exchange partially with the airways and partially with the alveoli.9,10

Figure 1 shows expirograms in which gases that exchange in the alveoli and those that exchange in the airways can be seen: phase 1 shows that the anatomic dead space empties. When fresh air is inhaled, low blood-soluble gases in the blood are transfered into the inspired air entering the alveoli. This action continues throughout the entire breathing cycle (inspiration and expiration). When the air is exhaled, molecules present in the conducting airways exit the lung first (phase 1). The airways take part in the gas exchange. Therefore, phase 1 does not exist in the expirogram for gases that exchange completely in the airways (e.g., ethanol). Therefore, an anatomical dead space cannot be given for such gases. For these gases, the end-exhaled breath concentration is always less than that in the alveoli. Phase 2 represents the transition from air residing in the dead space to alveolar air containing high partial pressures of gas. Phase 3 corresponds to air from the alveoli for low soluble gases and airways for high soluble gases. Alveolar air is expired via the conducting airways (phase 3). Molecules with high blood solubility are absorbed from the tissue and mucus lining the airway wall. In contrast to low blood-soluble gases, therefore, these gases present in large concentrations in the airway tissue and mucus for a given partial pressure. The partial pressure of soluble gas in the end-exhaled breath is always less than that in the alveoli due to the absorption–desorption process, which is the major mechanism of pulmonary gas exchange for any gas with a blood–air partition coefficient greater than 100.10

Graphic Jump Location
Fig. 1
F1 :

Phases of an expirogram from a low blood-soluble gas that exchanges in the alveoli (red line) and from a high blood-soluble gas that exchanges in the airways (blue line).9

Breath sampling procedures are threefold: alveolar, dead-space, and mixed sampling. Dead space air contains the first 150 mL of expiration air, and here no gas exchange takes place. Alveolar breath is the part of exhaled air coming after the dead-space air, whereas end-tidal air is the last fraction of expired air. Alveolar gas exchange is dependent on ventilation, pulmonary perfusion, and the blood:air partition.10 The term “end-exhaled breath” should be utilized for highly blood soluble VOCs, whereas the term “alveolar breath” may be used for low blood-soluble compounds.

It is often tacitly assumed that end-tidal air will reflect the alveolar concentration CA, which is related to the concentration of the VOC in mixed venous blood Cv, the substance-specific blood:gas partition coefficient λb:air, alveolar ventilation V˙A, and cardiac output Q˙c:11,12Display Formula

where λb:air states the diffusion equilibrium between alveoli and capillaries, V˙A describes the volume of VOC per unit time that arrives to the alveoli, and Q˙c describes the blood volume, which is pumped by the heart.

In mixed sampling, total breath (alveolar + dead space) is collected. Concentrations of endogenous VOCs in alveolar air are two to three times higher than those found in mixed expiratory samples because there is no dilution by dead space gas.13

There may be discrepancies between the true alveolar breath and the measured levels of water-soluble VOCs because their concentrations may be diluted on their way (a dilution effect). A single breath or multiple breath cycles can be applied for breath sampling.14 Breaths may considerably change from each other. The alveolar gas sample may not be represented by the composition of a single breath. Therefore, multiple breaths may be preferable.

Analysis of the breath molecules using laser-based methods is based on the absorption of electromagnetic waves by molecules under investigation. Therefore, we will explain the absorption phenomena in this section. The overall excited energy of a molecule is a sum of the electronic (Eel), vibrational (Evib), and rotational transitions (Erot) energy: Display Formula


Rotational spectra originate from the transitions between rotation levels of a given vibration level in an electronic state [Fig. 2(a)]. In this process, only the rotational quantum number (j) changes. These spectra lie in the spectral range from the microwave to the far infrared. The rotational–vibrational spectra consist of the transitions from rotation levels of a vibration level to rotation levels of another vibration level in the same electronic state [Fig. 2(b)]. The rotation and vibration (ν) quantum number change in these transitions. These spectra lie in the infrared spectral region. Electronic spectra consist of the transitions between the rotation levels of different vibration levels of an electronic state and rotation or vibration levels of another electronic level [Fig. 2(c)]. There are changes in all three quantum numbers (J, ν, and electronic quantum number). These spectra lie in the near-infrared or visible range or in the UV range.

Graphic Jump Location
Fig. 2
F2 :

Electronic, vibrational, and rotational quantum states and corresponding transitions: (a) rotational transition, (b) rotational–vibrational transition, and (c) electronic transition.

The number f of normal vibrations of a molecule is given as Display Formula

f=3N5(for linear molecules),(3)
Display Formula
f=3N6(for nonlinear molecules),(4)
where N is the number of internal degrees of freedom. For water molecules, N is equal to 3: there are three modes of normal vibrations (Fig. 3). Transitions from the ground state to the first excited state are known as fundamental transitions (e.g., 0,0,0 → 1,0,0). There are also other transition types: harmonic transitions (overtones) from the ground state (e.g., 0,0,0 → 2,0,0), combination transitions from the ground state (e.g., 0,0,0 → 1,2,1), harmonic transitions between excited states only (e.g., 1,0,0 → 3,0,0), and combination transitions between excited states only (e.g., 1,2,0 → 2,3,1).

Graphic Jump Location
Fig. 3
F3 :

Normal vibration modes of a water molecule and its characteristic ground and first excited vibrational energy levels.

The ratio of the transmitted intensity I and initial intensity I0 of laser beam transmitting through an absorbing medium at a particular frequency is given using the Beer–Lambert equation: Display Formula

or Display Formula
where α is absorption coefficient, ϕ is lineshape function, P is total pressure of the medium, x is mole fraction of the absorbing species, L is pathlength, n is the number-density of the absorbing species, and Si (S) is linestrength, which is defined as the total absorption per molecule. The linestrength of an absorption transition depends on the population in the lower quantum state, which is a function of the Boltzmann fraction, and the probability of the transition.15 The linestrength Slu for a transition lu normalized to the total number density N is expressed (in cgs units) by16Display Formula
where Nl and Nu are the number density of molecules in the states l and u, respectively, with their corresponding statistical weights gl and gu, h denotes the Planck constant, |μlu| is the transition dipole moment, and νlu (in cm1) is the frequency of the transition. In the case of thermal equilibrium of the system at temperature T, the equation above can be rearranged using the Maxwell–Boltzmann distribution law (SI units): Display Formula
with the energy of the lower state El and the total internal partition function Q(T). The equation for nu/gu is analogue. This leads to a form of S (in cm/molecule), which accords with the notation of the HITRAN database: Display Formula
where Ia is the abundance. Since the linestrength in HITRAN is defined for the reference temperature Tref=296  K, S(T) must be corrected using the tabulated values S(Tref), El, νlu, and Q(T):16Display Formula
where S is in units of [cm1/(mol·cm2)]. For units of [cm2atm1], the following temperature scaling can be used:15Display Formula

The conversion between S(T) and Si(T) is given as follows:15Display Formula


A typical lineshape of an isolated absorption line is centered at ν0 and has a linewidth Δν. The absorption line can be broadened due to natural broadening, Doppler broadening, and collision broadening. The homogeneous natural and collision broadening result in Lorentzian line shapes, whereas the Doppler broadening gives rise to a Gaussian line profile. The natural linewidth is determined by the lifetimes of the energy states (τ1 and τ2): Display Formula


The linewidth can be broadened due to the Doppler effect. The distribution of random velocities of molecules in gas is described by the Maxwellian velocity distribution function. Each velocity group gives rise to its own Doppler shift. The Doppler broadening is given by15Display Formula

with Doppler halfwidth [full-width at half-maximum (FWHM)]: Display Formula
ΔνD=ν0(8  kTln2mc2)1/2.(15)

The perturbations that shorten energy level lifetimes cause collisional broadening. The linewidth broadening occurs due to the Heisenberg uncertainty principle; shortening the energy level lifetime results in greater uncertainty of the energy of molecule, and thus the absorption lineshape becomes broader. The lineshape for collisional broadening takes the form of a Lorentzian function:15Display Formula

with collision halfwidth (FWHM): Display Formula
where m is the mass of molecule, d is the diameter of molecule, and p is the pressure of gas.

At low pressure, the near-Gaussian line shape is obtained because thermal motion is the dominant broadening mechanism. The line shape becomes near-Lorentzian with the increasing pressure. When the Doppler contribution can be neglected (ΔνcpΔνD), the line shape can be considered pure Lorentzian. In the case of ΔνcpΔνD, a pure Gaussian line shape is received. Molecular ro-vibrational transitions exhibit pressure broadened FWHM of 3  MHz (0.0001  cm1) at low pressures <5  mbar and 1.5  GHz (0.05  cm1) at atmospheric pressure, where Doppler widths of 30 … 90 MHz (0.0010.003  cm1) at room temperature are typically observed for molecular transitions in the infrared.16 When the condition ΔνcpΔνD is fulfilled (between 5 and 100 Torr for most lightweight gases), the line shape will be a combination of the Gaussian and Lorentzian, which is known as Voigt profile.

Different methods can be utilized to realize a wavelength-tunable single-mode operation: tunable CO2 lasers, tunable dye lasers, tunable solid-state lasers, tunable excimer lasers, nonlinear optical frequency conversion lasers, tunable free-electron lasers, and tunable external cavity semiconductor lasers. Within the framework of this work, the focus will only be on the tunable external cavity semiconductor lasers (ECLs). Further information on the other tunable lasers can be found in, e.g., Ref. 17.

Grating-coupled ECLs can be realized in Littrow or Littman configurations.18 In the Litrow configuration [Fig. 4(a)], there are a gain chip, a lens (or a parabolic mirror), and a diffraction grating. The gain chip serves as the gain medium for ECL setup. The lens (or parabolic mirror) is used to collimate the beam emitted from the gain chip. The collimated beam is diffracted by a grating. The first order of the diffracted beam is fed back into the gain medium. The wavelength is tuned by rotating the grating. Because two optical components (mirror or lens and grating) are used, the alignment is easy and the feedback is strong.

Graphic Jump Location
Fig. 4
F4 :

(a) Littrow and (b) Littman configuration.

In the Littman configuration [Fig. 4(b)], a mirror is added to the setup. The diffracted beam is sent back to the grating by a mirror and diffracted again. Finally, the beam reaches to the gain chip and it is amplified. The tuning is usually achieved by rotating the mirror. Because of the double pass of the beam on the grating, the optical feedback strength is reduced. However, the Littman setup has an advantage of the increased wavelength selectivity.

Generally, the Littrow setup is preferable for maximizing the tuning range. To optimize the optical feedback, the incident beam and first order of the diffracted beam have to be colinear. The directions of the various orders of the diffracted beam with respect to the normal are determined by Display Formula

where φi is the angle of the incident beam, φr,m is the angle of the diffracted beam of the m, order, d is the period of the grating, and λ is the wavelength. The wave whose wavelength satisfies the condition: Display Formula
is selected by the grating and is amplified. For the other wavelengths, it can be written as follows: Display Formula
where Δφ is the angle between incident and diffracted beam [see Fig. 4(a)] and is defined as Display Formula

Finally, the wave, which makes an angle Δφ, is focalized on the facet at a distance: Display Formula

where f is the focal length of the lens.

The reflectivity of the facet inside the external cavity setup is very important for a properly operating ECL. This reflectivity should be as low as possible to effectively suppress the Fabry–Perot modes of the gain chip, so a pure single mode can be selected and amplified. The reflectivity of the outcoupling facet of the gain chip can be reduced by various ways, such as the tilted gain stripes, the buried facets, and the antireflection (AR) coating. The tilting of the waveguide decreases in the internal Fresnel reflection from the facet and therefore increases the losses.19,20 In the case of buried facets, a semi-insulating material between the end of the guide and the facet, which reduces the coupling of the beam back into the active region, is utilized.21 When a facet is AR-coated, a phase shift of π and, therefore, an optical path difference of λ/2 occur after the reflection in the optical thick medium. The reduction of the reflectivity of the AR-coated facet is achieved by the destructive interference between the waves reflected from the front and rear sides of the layer. The reflectivity reduction causes a decrease in the competition between the modes of the small cavity (gain chip) and of the large cavity (ECL) for oscillation. This situation leads to enhancement of the coupling efficiency to the EC, to increased output power of the ECL, and to enlarged coarse/fine tuning range and the range of the stable operating region.

Wavelength Selection

Tuning mechanism of an ECL can be explained by considering a sketch for the gain and losses in the ECL setups (Fig. 5).17 As one can see in Fig. 5, there are three loss mechanisms: (1) Filter loss due to the wavelength selectivity property of the wavelength selective element (e.g., grating), (2) coupling and output loss between the guided waves in the gain chip and vacuum in ECL, and (3) facet etalon loss, which is modulated at the period of the internal Fabry Perot mode spacing Δνint. The facet etalon loss is caused by the interference between the reflections from facets.

Graphic Jump Location
Fig. 5
F5 :

Gain and losses for ECLs.

The injection current is controlled to tune the gain profile of the gain chip and the internal cavity mode frequency. To vary the external mode frequency, the grating can be rotated and/or the EC length can be changed by applying the voltage to a piezoelectric transducer controlling the grating position. For a perfect operation of the ECLs, the facet etalon loss ripple should be weak and the bandwidth of the wavelength selective element loss should be narrow compared with the period of the ripple. The bandwidth of the grating can be estimated by the following equation: Display Formula

where λ is the wavelength, a is the diameter of the collimated laser beam, and θ is the incident angle of the beam. Only the EC mode, which is the closest to the loss minimum of the wavelength selective element, can oscillate. The other modes see higher losses and cannot be amplified. Tuning to the next EC mode can be achieved by varying the loss minimum of the wavelength selective element. In the case of the grating-coupled ECLs in Littrow configuration, this variation occurs by rotating the grating [Fig. 6(a)] and/or changing the cavity length [Fig. 6(b)]. A strong facet etalon loss and a wider bandwidth of the wavelength-selective element loss lead to tuning with EC mode jumps approximately equal to Δνint. This discontinuity is called “tuning gap” [Figs. 6(a) and 6(b)]. The shift in the external mode frequency exceeding half the free spectral range of the external cavity gives rise to mode hops. At this point, the adjacent external mode has moved closer to the peak of the internal mode, and the laser oscillation hops to this new mode.22 To achieve the tuning without any gap (tuning on the EC mode, mode-hop-free tuning) [Fig. 6(c)], there are different methods: (1) the injection current, the grating angle, and/or the external cavity length should be simultaneously controlled;23 (2) the grating has to be mounted on a pivot arm, which is rotated about an axis lying near the intersection of the grating and laser output facet plane;24 and (3) an electro-optically tunable single-mode ECL has to be built.25

Graphic Jump Location
Fig. 6
F6 :

Possibility of changing the grating position and corresponding wavelength tuning. (a) The grating is rotated around its the center, (b) the grating is moved back and forth, and (c) the grating is rotated around an axis lying near the intersection of the grating and laser output facet plane.

Spectral Tuning Characteristics

To consider the gain chip and the ECL as a combined laser system with a compound cavity, the intracavity facet should have a low reflectivity, and the coherence length should be longer than the distance from the gain chip to the grating. The tuning range is limited by the gain bandwidth of the used gain chip that is an inherent characteristic of the gain medium. The separation between the FP and EC modes must be larger than the bandwidth of the grating for a single-mode operation. This can be achieved, for instance, using an appropriate AR coating.

The theoretically reachable tuning range of an ECL for wavelength is given by26Display Formula

and, equivalently, the frequency tuning range: Display Formula
where λ is the emission wavelength, 2γgain is the FWHM of the modal gain spectrum, E0 is the photon energy at the gain maximum, g0 is the apex of the modal gain spectrum, h is the Planck constant, L is the chip length, RG is the first-order efficiency of the grating, and TL is the lens transmission. The relative tuning is given by Δλ/λc, where λc=(λmin+λmax)/2 is the center wavelength and λmin and λmax are the shortest and the largest achievable wavelengths, respectively, in the ECLs. From these equations, one can conclude that the following should be considered to increase the tuning range: a gain chip with a broad gain bandwidth should be used in an ECL setup, and reflectivity from the facet inside the ECL should be decreased by one of the methods mentioned above. In Table 1, some ECLs reported in the spectral range between 1.0 and 11.5  μm are listed.

Table Grahic Jump Location
Table 1Some ECLs in the spectral range between 1 and 11  μm.

The laser linewidth (FWHM of the lineshape) describes the spectral purity of a laser (Fig. 7). The linewidth of a semiconductor laser is given by the Schawlow–Townes formula modified by Henry (Ref. 37 and references therein): Display Formula

where ν is the lasing frequency, α is the total losses of the cavity, αm is the mirror losses of the output facet, P is the output power, Δνc=αc/(2πn) is the FWHM of the cavity resonance, n1t and n2t are the populations of the upper and lower states of the laser transition at threshold, respectively, and αe=(dnr/dN)/(dni/dN) is the linewidth enhancement factor introduced by Henry. The imaginary (nr) and real parts (ni) of the refractive index change during the relaxation oscillations of the field intensity. Varying of the imaginary part of the index causes the gain to change, whereas the change in the real part of the refractive index gives rise to an additional phase shift. Therefore, the laser linewidth is broadened by a factor (1+αe2) compared with the Schawlow–Townes linewidth. The lower energy level of semiconductor lasers is not empty, which is different from most other lasers. As laser photons in the valence band may reradiate by spontaneous emission, phase fluctuations of the laser field increase. This is another phenomenon that leads to the linewidth broadening of semiconductor lasers.

Loh et al.38 investigated experimentally the influence of the grating reflectivity, grating resolution, and AR coating on the intrinsic linewidth of a grating-coupled ECL. They found that the linewidths are little affected by the presence of an AR coating on the diode facet and that the grating resolution has a larger effect on the linewidth than the grating reflectivity.38

The laser linewidth can be determined using several methods, such as optical cavity measurements,29 heterodyne measurements,37 homodyne/heterodyne interferometry using a long interferometer delay,39 measurements from the frequency noise power spectral density,40 measurements from root-mean-square frequency deviation,41 and saturation spectroscopy.42

Tkach and Chraplyvy43 reported that there are five distinct regimes of feedback effects, with well-defined transitions between them as a function of feedback power ratio. At the lowest levels of feedback down to 90  dB (regime I), linewidth narrowing or broadening is observed depending on the phase of the feedback. Regime II is seen at a feedback level (up to around 45  dB), which depends on the distance to the external reflector, and frequency splitting of the mode and mode hopping are observed. Regime III is found at a level (about 39 to 45  dB), which does not depend on the distance to the reflection. In this regime, the mode hopping is suppressed, and the laser is observed to operate on a single narrow line. Regime IV starts at 39  dB and is independent from the distance to the reflection. In this regime, line broadening is observed due to the coherence collapse effect. Relaxation oscillation sidebands that present in each mode in the optical spectrum increase in strength giving rise to the line broadening. Consequently, coherence time decreases. Regime V is independent of distance too and observed at the highest levels of feedback (usually greater than 10  dB). AR coating of output facet is needed to achieve this regime. In this regime, the operation of the laser occurs on a single longitudinal mode with narrow linewidth for all phases of the feedback and the laser is relatively insensitive to additional external optical perturbations.43

Frequency Stabilization

There are different noise sources, which can affect on the laser’s linewidth: (1) fundamental noises (white, shot, and thermal noises): the intrinsic Lorentzian linewidth can be reduced using optical feedback44 and (2) flicker (1/f) noises and technical noises (noises from the laser driver, piezo drift, mechanical vibrations, temperature fluctuations, acoustic noises, air pressure fluctuations). These noise sources broaden the laser linewidth. This Gaussian linewidth can be measured using the methods given above. For reduction of these noises, with other word, stabilization of a laser’s frequency over time and consequently narrowing the linewidth of laser, there are different stabilization techniques, such as side-of-fringe stabilization,45 on-peak-locking, Paul–Drever–Hall (PDH) method.46 Using these techniques, the sub-kHz-linewidth47 and sub-Hz-linewidth48 ECLs have been reported.

Rate Equations

A set of rate equations can describe the time evolution of the photon flux S and the electronic populations of the upper and lower states of the radiative transition (n3,n2). Yang et al.26 conducted rate equations analysis of EC quantum cascade lasers. They considered a four-level quantum cascade laser structure. The EC system was simplified to an effective system with both the internal cavity and the EC having the same cavity length Leff. To consider the longer round-trip time of the EC, they used the effective speed of light introduced by Maulini.49 The EC can be simplified to have an effective cavity length of Leff but with a reduced speed of light c/ρcav, where ρcav=(τint+τext)/τint, τint and τext are the round-trip times of the internal cavity and the EC, respectively. The rate equations then are given by26,49,50Display Formula

Display Formula
Display Formula
Display Formula
where J is the injection current density, e is the unit charge, n3 and n2 are the sheet carrier density of upper laser level 3 and lower laser level 2, τij represents the nonradiative scattering lifetime from level i to level j (i, j=1, 2, 3), 1/τ3=1/τ32+1/τ31, τ2=τ21, τsp is the spontaneous emission time (calculated to be 1.4×107  s), β0.001 is the fraction of the spontaneous light emitted in the lasing mode, SFP is the photon flux density (per unit length per second) generated by the laser chip facets, SEC is the photon flux density for a certain mode selected by the external cavity, αFP and αEC represent the total losses corresponding to the frequencies of FP mode and the EC mode, gFP and gEC denote the values of the gain cross-section corresponding to the frequencies of FP mode and the EC mode and are determined by26Display Formula
Display Formula
where z322 is the matrix dipole element of the optical transition, Γ is the optical confinement factor of the laser, ε0 is the free-space permittivity, Lp is the period length of the QCL, λ0 is the emission wavelength, 2γ32 is the FWHM of the gain spectrum, ν is the lasing frequency, and h is the Planck constant.

The ECL-based techniques can be effectively used for the detection and quantification of breath VOCs at concentrations ranging from parts per million by volume (ppmv) to the parts per trillion (pptv) level. The laser sources should have high sensitivity and selectivity for breath molecule sensing. For selectivity, an absorption line that is free of interference from other species in gas mixture should be used. For the best sensitivity, a strong molecular absorption line should be chosen. To achieve a higher signal-to-noise ratio (SNR), either the noise levels should be decreased and/or the signal should be enhanced. For decreasing the noise levels, there are different methods, such as wavelength/frequency modulation spectroscopy. The signal can be further increased by obtaining an effective long pathlength using multipass cells or cavity enhancement methods. ECL-based spectroscopic detection techniques, such as multipass cell absorption method (MCAB), CRDS, CEAS, ICOS, PAS, and QEPAS, provide the features of high-sensitivity, high-selectivity, low costs, real-time measurements, and PoC function. Furthermore, laser-based techniques can be implemented without the need for radioactive labels. Therefore, these methods are candidates for performing accurate breath diagnostics in hospitals. In Table 2, some important breath molecules, detection wavelength, detection limit, and methods are summarized.

Table Grahic Jump Location
Table 2Some important breath molecules, detection wavelengths, detection limits, and methods.
Multipass Absorption Spectroscopy

In the MCAB,75 the absorption signal can be enhanced via longer optical pathlengths using optical multipass cells. There are four types of multipass cells most commonly applied: White cells,76,77 Herriott cells,78 Chernin cells,79 and astigmatic mirror multipass cells.80 Beer–Lambert law of linear absorption is used to determine the absorption coefficient α(ν): Display Formula

where I(ν) is transmitted intensity of the radiation after a distance Leff through the medium, I0 is initial intensity of the radiation, Leff is effective path-length, σ(ν) is absorption cross-section of the absorbing species, N is molecular concentration (the number density of the absorbing species), S is linestrength (i.e., the total absorption per molecule) of the absorbing species at temperature T, ϕ(νν0) is line shape function for the particular absorption line.

Using Eq. (33), the absolute number density N can be calculated by integrating the absorption coefficient: Display Formula


The Herriott cell consists of two mirrors (Fig. 8). The input mirror has a small hole. A laser beam is injected into a Herriott cell through the hole and it is reflected multiple times between the two mirrors. The light path might be lengthened 10 or 100 times inside the cell. The laser beam leaving the cell through the same or the other hole reaches the detector. The sensors employing a multipass gas cell have been reported for many trace gases, such as acetaldehyde,51 ethane,61 formaldehyde,67 and methane.71,72

Graphic Jump Location
Fig. 8
F8 :

Herriott multipass cell.

Cavity Ringdown Spectroscopy

CRDS81 is based on the observation of the decay rate of a laser beam stored in an optical cavity formed from ultra-high reflective mirrors (Fig. 9) and is first demonstrated by O’Keefe and Deacon.82 In this case, the absorption coefficient is determined using the following equation: Display Formula

where τ0 and τ are the ringdown time without and with the absorbing gas present, c is the light velocity.

Graphic Jump Location
Fig. 9
F9 :

Cavity ringdown absorption cell.

A laser beam, which is coupled into the optical cavity, is reflected back and forth inside the cavity. When beam is reflected every time, a small fraction of this light leaks out of the cavity. Many trace gases, such as acetone,52 carbon dioxide,57 ethanol,63 isoprene,70 and nitric oxide,73,74 have been analyzed using the CRDS technique.

Integrated Cavity Absorption Spectroscopy and Cavity-Enhanced Absorption Spectroscopy

ICOS83 and CEAS84 are the same method. ICOS is used in the context of pulsed excitation, and CEAS is used in the context of continuous-wave cavity excitation. This method is based on analyzing the intensity of light transmitted through a high finesse optical cavity as a function of wavelength (Fig. 10). To reduce the noise levels and therefore to improve the detection limit, the laser beam can be coupled into the optical cavity at a small angle with respect to the cavity axis or away from the central axis of the optical cavity. This is known as “off-axis alignment,” which has advantage over on-axis alignment: it is almost insensitive to vibrations. In the case of ICOS and CEAS, Eq. (36) is applied to extract the absorption coefficient: Display Formula

where R is mirror reflectivity and d is the distance between mirrors.

Graphic Jump Location
Fig. 10
F10 :

High finesse optical cavity for ICOS or CEAS.

The gases—such as acetone,53 ammonia,54 carbon monoxide,54 carbonyl sulfide,60 and ethane60—have been analyzed using the ICOS or CEAS techniques.

Photoacoustic Spectroscopy and Quartz-Enhanced Photoacoustic Spectroscopy

In the PAS,85 the acoustic waves resulting from the absorption of radiation (photoacoustic effect) are detected and analyzed using a sensitive microphone (Fig. 11). The PAS response S to optical absorption is given by Eq. (37): Display Formula

where α is the absorption coefficient per unit concentration of the target species, C is the concentration of the target species, P is the optical power, f is the PA sound frequency, A is the resonator cross-section area, and k is a constant describing the microphone transfer function and other system parameters. Usually, PA resonators are designed for f values in the 500 to 4000 Hz range and have Q factors of 20 to 200.86

Another approach to photoacoustic detection of trace gases is QEPAS,86,87 which utilizes a quartz tuning fork (QTF) as a sharply resonant acoustic transducer. The QTF is a piezoelectric element. To further increase the QEPAS signal, a microresonator can be added to the QEPAS sensor system.88 Also, to achieve the best QEPAS detection sensitivity and selectivity, it is important to choose the optimum operating pressure for a QEPAS sensor system. The PAS- or QEPAS-based sensors have been used for the detection of many gases, such as ammonia,55,56 carbonyl disulfide,58 carbon monoxide,59 ethanol,62 ethylene,6466 and hydrogen sulfide.68,69

Owing to its noninvasiveness, ease of use, easy repeatability, and painlessness, breath analysis has attracted considerable attention for scientific and clinical studies. Until now, many studies have demonstrated that the human breath air contains more than 1000 different molecules.2,89 The matrix of exhaled breath air is much less complicated compared to serum or urine. A few atmospheric molecules in % concentration (N2, O2, H2O, CO2), inorganic gases (e.g., NO, NH3, and CO), VOCs at the levels of parts per million (ppm), parts per billion (ppb), or parts per trillion (ppt) (e.g., CH4, HCO2, N2O), and nonvolatile substances measured in breath condensate (e.g., isoprostanes, cytokines, leukotrienes) have been found to be present in exhaled human breath (Table 3). Spectral fingerprints of these molecules span from the UV to the IR spectral regions. The mixture is exhaled at temperatures between 34°C and 37°C while relative humidity may range from 91% to 96% in oral exhalations and from 82% to 85% in nasal exhalations.90 The mixture of exhaled breath air varies from person to person both quantitatively and qualitatively.91 Analysis of some VOCs, which present in everyone, can provide a safe method for studying human metabolism, detecting diseases, clinical diagnosis, therapeutic monitoring, and investigation of potential adverse effects of new drug therapy. In Table 4, some important breath molecules, concentration, physiological basis, and related diseases are summarized. The 872 breath VOCs from apparently healthy humans have been reported in the review article of de Lacy Costello et al.89

Table Grahic Jump Location
Table 3Main breath molecules.3
Table Grahic Jump Location
Table 4Some important breath molecules, concentration, physiological basis, and related diseases.

It has been known for a long time that the smell of exhaled breath is an indicator for diseases.122 The work of Pauling et al., in which hundreds of VOCs were determined in normal human breath at the levels of parts per billion (ppb) or lower, has been accepted as the beginning of the modern breath analysis. In spite of the long investigation time, this method has not yet been introduced as a standard tool into clinical diagnosis. Very few diseases are actually diagnosed using exhaled breath air analysis; for example, breath urea analysis. There are several potential confounders and issues blocking the introduction of breath analysis into clinical applications. To widely use the breath test for more specific diseases in hospitals, the following issues should be addressed:123125

  • The ease of use, low response time, low detection limit, high sensitivity, high selectivity, portability, compactness, small weight, reliability for long time periods, low energy consumption, and low construction and maintenance cost devices with good precision and accuracy being able to perform the online measurements, point-of-care detection, operating at room temperature are very important for widespread clinical use of the breath tests. Such devices are currently being developed. The availability of such devices can represent a breakthrough for clinical use of breath analysis.
  • The correlation between diseases and breath biomarkers should be understood in depth. The issues of biointerferences should be addressed: a molecule present in breath can be established as a biomarker for more than one disease; one specific disease can be detected by more than one breath molecule.
  • Exhalation physiology of VOCs should be understood exactly. The sources of the breath biomarkers inside the human body should be known very well. Also, we should be able to explain whether the breath molecules are exogenous (not biochemically produced in the body and therefore not a part of metabolic processes; e.g., taken up via skin absorption, inhalation from the ambient air, ingested foods, and beverages) or endogenous (generated by metabolic processes at various organs within the human body, excreted into blood, exchanged with the inhaled air in the lung via the alveolar interface, and ultimately exhaled) in origin. The other alternative sources of VOCs are bacteria in the gut, viruses, fungi, saliva, and mucus created in the respiratory tract. VOCs may also be generated in the airways and the oral cavity by bacterial infections. Cells or tissues in the mouth, nose, sinuses, airway, and the gastrointestinal tract can also generate molecules in exhaled breath that are not found in circulating blood (such as nitric oxide).123 Therefore, the concentrations of breath biomarkers can measure higher when there are also exogenous sources, which contribute molecules to exhaled breath. The levels of some molecules are dynamic, e.g., the levels of ammonia, which is the by-product of protein metabolism. Because protein metabolism can change based on many factors, such as pH, food/beverage consumption, medications, and exercise, the concentrations of ammonia are changeable depending on the protein metabolism.
  • The acceptable standardization and normalization of procedures for sampling, analysis, and background correction should also be determined. The correction for background concentrations of VOCs is necessary for distinguishing endogenous substances from exogenous contaminants, which have no diagnostic value. There are different approaches for background correction. One approach is to calculate “alveolar gradients” (expiratory concentrations—ambient or inspiratory concentrations).126 Calculating alveolar gradients is easy to realize, but this method does not consider the complexity of the inhalation/exhalation function, as given above and below. Another approach is to try washing out the lung from ambient concentrations by subjects breathing pure air for a certain time before measurement.127 However, the washout of the entire body may take days or weeks, depending on the identity of the molecule.4 There is no consensus for a standard method to enable the background levels to be subtracted.4 Currently, accepted standardized techniques for breath air sampling and analysis do not exist. Therefore, some results obtained by different groups are not always consistent. The American Thoracic Society published the first guidelines in 1999128 and updated guidelines in 2005129 for the online and offline measurement of exhaled lower respiratory nitric oxide and nasal nitric oxide in adults and children.
  • An acceptable unit should be used. It is common to express the concentration of breath VOCs as ppm by volume (ppmv), ppbv, or pptv. An alternative is to express exhaled concentrations in mol or weight per volume.
  • Increased or decreased levels of biomarkers in breath can be an indicator for a large number of diseases. Therefore, the range between disease or wellness states should be defined exactly. For this purpose, a larger number (several hundreds) of human subjects should be investigated.
  • How the concentrations of breath molecules under investigation depend on the parameters, such as ethnicity, gender, pulmonary function, age, body mass index, smoking and alcohol consumption, etc, should also be defined. Dependence of diffusion of VOCs from blood to alveolar air across the alveolar–capillary membrane on their physicochemical properties, such as polarity, solubility in fat, Henry partition constant, and volatility, should also be considered.130 The low blood-soluble gases exchange in the alveoli, while the gas exchange of highly blood-soluble VOCs takes place in the airways rather than alveoli.9 A lipid-soluble molecule could be stored in tissues, not well perfused by blood, and be released more slowly than a similar molecule with hydrophilic properties that is not stored.123 The physiological parameters, such as blood pressure, heartbeat rate, blood:gas partition coefficient, cardiac output, ventilation/perfusion ratio of the lung, and alveolar ventilation, may also have an effect on the exhalation of breath biomarkers.131 These effects may be negligible when the subject is healthy. However, in patients with severe hemodynamic or pulmonary dysfunction, the normalization or correction of data can be needed by means of end tidal and/or arterial pCO2 measurement.13
  • Data interpretation is a crucial step in breath analysis.125 Confounding variables, which have a real statistical correlation with the disease and a breath marker, are environmental compounds, medical history, gender, age, weight, physiological parameters like cardiac output and respiratory flow, effects of breath sampling, sample storage, and the analytical technique. These effects have to be considered and addressed during data interpretation.125

To address all these issues and to move breath analysis forward, collaboration among optical engineers, experts in breath analysis, clinicians, and spectroscopists is required.

Several issues should be addressed for advances in the field of breath monitoring and for widespread usage of this method as a routine clinical tool. Breath analysis has not yet been used in hospitals as a clinical tool due to lack of generally accepted evaluation criteria for data. Multiple potential confounders, such as interference from other molecules, smoking and alcohol consumption, dermal absorption, inhalation from the ambient air, and so on, should be carefully considered when performing measurements of breath biomarkers.

Breath research may be performed as (1) exploratory research (the discovery of new breath molecules, e.g., breath markers for chronic kidney disease), (2) focused research (the in-depth evaluation of a molecule whose biochemistry is already well-understood, e.g., breath acetone), and (3) breath condensate (the analysis of nonvolatile molecules).123

Breath analyses conducted by the techniques, such as GC-MS, SIFT-MS, and PTR-MS, are still in laboratory research. The development of inexpensive, real-time, portable breath monitors is a key factor in the advancement of breath tests. Such devices should be moved from laboratory research to commercial reality. Laser-based spectroscopic instruments are a good candidate for this field. The current laser-based methods are not satisfactory for medical applications as routine diagnostics tool. Technologies in laser physics and optical engineering should be further improved. The development of absorption spectroscopy instruments with lasers may provide certain advantages for sensitive examination of not only small molecules but also heavy molecules with relatively broad spectral structures. The issue of spectral interferences can be eliminated using single-mode tunable laser. Another type of interference, which should be considered, is biointerferences; for instance, ammonia is a biomarker for both asthma and chronic kidney diseases. This challenging field requires the extensive interdisciplinary collaborations.

The availability of a cheap breath analyzer can even enable patients to measure their own biomarker levels at home. This situation can result in performing more large cross-sectional and longitudinal in-depth studies with many data points, instead of making just small pilot studies limiting the gain of widespread usage of breath analysis. One may expect that the issues mentioned in the present manuscript may be soon addressed and breath analysis may be introduced into clinical practice.

The author has no relevant financial interests and no other potential conflicts of interest to disclose.

The author thanks Yasar Kemal Erdogan for helping to compile Table 2.

Ciaffoni  L., , Peverall  R., and Ritchie  G. A. D., “Laser spectroscopy on volatile sulfur compounds: possibilities for breath analysis,” J. Breath Res.. 5, , 024002  (2011).CrossRef
Wang  C., and Sahay  P., “Breath analysis using laser spectroscopic techniques: breath biomarkers, spectral fingerprints, and detection limits,” Sensors. 9, , 8230 –8262 (2009). 0746-9462 CrossRef
Wojtas  J.  et al., “Ultrasensitive laser spectroscopy for breath analysis,” Opto-Electron. Rev.. 20, , 26 –39 (2012). 1230-3402 CrossRef
Risby  T. H., and Tittel  F. K., “Current status of mid-infrared quantum and interband cascade lasers for clinical breath analysis,” Opt. Eng.. 49, , 111123  (2010).CrossRef
Bielecki  Z.  et al., “Application of quantum cascade lasers to trace gas detection,” Bull. Pol. Acad. Sci. Tech. Sci.. 63, , 515 –525 (2015).CrossRef
Kim  S. S.  et al., “Potential and challenges for mid-infrared sensors in breath diagnostics,” IEEE Sens. J.. 10, , 145 –158 (2010). 1530-437X CrossRef
Chow  K. K., , Short  M., and Zeng  H., “A comparison of spectroscopic techniques for human breath analysis,” Biomed. Spectrosc. Imaging. 1, , 339 –353 (2012).CrossRef
McCurdy  M. R.  et al., “Recent advances of laser-spectroscopy-based techniques for applications in breath analysis,” J. Breath Res.. 1, , 014001  (2007).CrossRef
Anderson  J. C., and Hlastala  M. P., “Breath tests and airway gas exchange,” Pulm. Pharmacol. Ther.. 20, , 112 –117 (2007).CrossRef
Anderson  J. C., , Babb  A. L., and Hlastala  M. P., “Modeling soluble gas exchange in the airways and alveoli,” Ann. Biomed. Eng.. 31, , 1402 –1422 (2003). 0090-6964 CrossRef
Farhi  L. E., “Elimination of inert gas by the lung,” Respir. Physiol.. 3, , 1 –11 (1967).CrossRef
King  J.  et al., “A mathematical model for breath gas analysis of volatile organic compounds with special emphasis on acetone,” J. Math. Biol.. 63, , 959 –999 (2011). 0303-6812 CrossRef
Miekisch  W., , Schubert  J. K., and Noeldge-Schomburg  G. F. E., “Diagnostic potential of breath analysis—focus on volatile organic compounds,” Clin. Chim. Acta. 347, , 25 –39 (2004). 0009-8981 CrossRef
O’Hara  M. E.  et al., “Development of a protocol to measure volatile organic compounds in human breath: a comparison of rebreathing and on-line single exhalations using proton transfer reaction mass spectrometry,” Physiol. Meas.. 29, , 309 –330 (2008). 0967-3334 CrossRef
Webber  M. E., “Diode laser measurements of NH3 and CO2 for combustion and bioreactor applications,” Dissertation, Stanford University (2001).
Welzel  S., “New enhanced sensitivity infrared laser spectroscopy techniques applied to reactive plasmas and trace gas detection,” Dissertation, Ernst-Moritz-Arndt Universität Greifswald (2009).
Duarte  F. J., Eds., Tunable Lasers Handbook (Optics and Photonics). ,  Eastman Kodak Company ,  Rochester, New York  (1995).
Littman  M. G., and Metcalf  H. J., “Spectrally narrow pulsed dye laser without beam expander,” Appl. Opt.. 17, , 2224 –2227 (1978). 0003-6935 CrossRef
Scifres  D. R., , Streifer  W., and Burnham  R. D., “GaAs/GaAlAs diode lasers with angled pumping stripes,” IEEE J. Quantum Electron.. 14, , 223 –227 (1978).CrossRef
Robbins  M. C., , Buus  J., and Robbins  D. J., “Analysis of antireflection coatings on angled facet semiconductor laser amplifiers,” Electron. Lett.. 26, , 381  (1990). 0013-5194 CrossRef
Dutta  N. K.  et al., “Fabrication and performance characteristics of buried-facet optical amplifiers,” J. Appl. Phys.. 67, , 3943 –3947 (1990). 0021-8979 CrossRef
Dutta  S., , Elliott  D. S., and Chen  Y. P., “Mode-hop-free tuning over 135 GHz of external cavity diode lasers without antireflection coating,” Appl. Phys. B. 106, , 629 –633 (2012).CrossRef
Schremer  A. T., and Tang  C. L., “External-cavity semiconductor laser with 1000 GHz continuous piezoelectric tuning range,” IEEE Photonics Technol. Lett.. 2, , 3 –5 (1990). 1041-1135 CrossRef
Favre  F., and Le Guen  D., “82 nm of continuous tunability for an external cavity semiconductor laser,” Electron. Lett.. 27, , 183 –184 (1991).CrossRef
Levin  L., “Mode-hop-free electro-optically tuned diode laser,” Opt. Lett.. 27, , 237  (2002). 0146-9592 CrossRef
Yang  Q. K.  et al., “Rate equations analysis of external-cavity quantum cascade lasers,” J. Appl. Phys.. 107, , 043109  (2010).CrossRef
Bayrakli  I., “Frequency stabilization at the sub-kilohertz level of an external cavity diode laser,” Appl. Opt.. 55, , 2463 –2466 (2016).
Bayrakli  I., “Actively frequency-stabilized external cavity diode laser with a linewidth of 2.9 kHz,” Opt. Quantum Electron.. 48, , 1 –8 (2016). 0306-8919 CrossRef
Bayrakli  I., and Akman  H., “Ultrasensitive, real-time analysis of biomarkers in breath using tunable external cavity laser and off-axis cavity enhanced absorption spectroscopy,” J. Biomed. Opt.. 20, (3 ), 037001  (2015). 1083-3668 CrossRef
Bennetts  S.  et al., “External cavity diode lasers with 5 kHz linewidth and 200 nm tuning range at 1.55  μm and methods for linewidth measurement,” Opt. Express. 22, , 10642 –10654 (2014). 1094-4087 CrossRef
Tan  S.  et al., “Low-threshold, high SMSR tunable external cavity quantum cascade laser around 4.71  μm,” Opt. Quantum Electron.. 45, , 1147 –1155 (2013). 0306-8919 CrossRef
Wysocki  G.  et al., “Widely tunable mode-hop free external cavity quantum cascade laser for high resolution spectroscopic applications,” Appl. Phys. B. 81, , 769 –777 (2005).CrossRef
Dougakiuchi  T.  et al., “Broadband tuning of external cavity dual-upper-state quantum-cascade lasers in continuous wave operation,” Appl. Phys Express. 4, , 102101  (2011). 1882-0778 CrossRef
Maulini  R.  et al., “External cavity quantum-cascade laser tunable from 8.2 to 10.4  μm using a gain element with a heterogeneous cascade,” Appl. Phys. Lett.. 88, , 201113  (2006). 0003-6951 CrossRef
Hugi  A.  et al., “External cavity quantum cascade laser tunable from 7.6 to 11.4  μm,” Appl. Phys. Lett.. 95, , 061103  (2009). 0003-6951 CrossRef
Wittmann  A.  et al., “Heterogeneous high-performance quantum-cascade laser sources for broad-band tuning,” IEEE J. Quantum Electron.. 44, , 11  (2008). 0018-9197 CrossRef
Maulini  R.  et al., “Continuous-wave operation of a broadly tunable thermoelectrically cooled external cavity quantum-cascade laser,” Opt. Lett.. 30, , 2584  (2005). 0146-9592 CrossRef
Loh  H.  et al., “Influence of grating parameters on the linewidths of external-cavity diode lasers,” Appl. Opt.. 45, , 9191 –9197 (2006). 0003-6935 CrossRef
Okoshi  T., , Kikuchi  K., and Nakayama  A., “Novel method for high resolution measurement of laser output spectrum,” Electron. Lett.. 16, , 630 –631 (1980). 0013-5194 CrossRef
Di Domenico  G., , Schilt  S., and Thomann  P., “Simple approach to the relation between laser frequency noise and laser line shape,” Appl. Opt.. 49, , 4801 –4807 (2010). 0003-6935 CrossRef
Elliott  D. S., , Rajarshi  R., and Smith  S. J., “Extracavity laser band-shape and bandwidth modification,” Phys. Rev. A. 26, , 12 –18 (1982).CrossRef
Mukherjee  N.  et al., “Linewidth measurement of external grating cavity quantum cascade laser using saturation spectroscopy,” Appl. Phys. Lett.. 92, , 111116  (2008). 0003-6951 CrossRef
Tkach  R. W., and Chraplyvy  A. R., “Regimes of feedback effects in 1.5-μm distributed feedback laser,” J. Lightwave Technol.. 4, , 1655 –1661 (1986). 0733-8724 CrossRef
Samutpraphoot  P.  et al., “Passive intrinsic-linewidth narrowing of ultraviolet extended-cavity diode laser by weak optical feedback,” Opt. Express. 22, , 11592 –11599 (2014). 1094-4087 CrossRef
Kourogi  M., and Ohtsu  M., “Novel optical frequency discriminator for FM noise reduction of semiconductor lasers,” Opt. Commun.. 81, , 204 –208 (1991). 0030-4018 CrossRef
Drever  R. W. P.  et al., “Laser phase and frequency stabilization using an optical resonator,” Appl. Phys. B. 31, , 97 –105 (1983).CrossRef
Zhao  Y.  et al., “100-Hz linewidth diode laser with external optical feedback,” IEEE Photonics Technol. Lett.. 24, , 1795 –1798 (2012). 1041-1135 CrossRef
Hirata  S.  et al., “Sub-hertz-linewidth diode laser stabilized to an ultralow-drift high-finesse optical cavity,” Appl. Phys. Express. 7, , 022705  (2014). 1882-0778 CrossRef
Maulini  R., “Broadly tunable mid-infrared quantum cascade lasers for spectroscopic applications,” PhD Thesis, Université de Neuchâtel (2006).
Faist  J.  et al., Semiconductors and Semimetals Quantum Cascade Laser. , , Liu  H. C., and Capasso  F., Eds., Vol. 66, ,  Academic Press ,  San Diego  (2000).
Kamat  P. C.  et al., “Measurement of acetaldehyde in exhaled breath using a laser absorption spectrometer,” Appl. Opt.. 46, , 3969 –3975 (2007). 0003-6935 CrossRef
Sun  M.  et al., “A fully integrated standalone portable cavity ringdown breath acetone analyzer,” Rev. Sci. Instrum.. 86, , 095003  (2015). 0034-6748 CrossRef
Ciaffoni  L.  et al., “Demonstration of a mid-infrared cavity enhanced absorption spectrometer for breath acetone detection,” Anal. Chem.. 85, , 846 –850 (2013). 0003-2700 CrossRef
Thorpe  M. J.  et al., “Cavity-enhanced optical frequency comb spectroscopy: application to human breath analysis,” Opt. Express. 16, , 2387 –2397 (2008). 1094-4087 CrossRef
Kosterev  A. A., and Tittel  F. K., “Ammonia detection by use of quartz-enhanced photoacoustic spectroscopy with a near-IR telecommunication diode laser,” Appl. Opt.. 43, , 6213 –6217 (2004). 0003-6935 CrossRef
Lewicki  R.  et al., “Real time ammonia detection in exhaled human breath using a distributed feedback quantum cascade laser based sensor,” Proc. SPIE. 7945, , 79450K  (2011). 0277-786X CrossRef
Crosson  E. R.  et al., “Stable isotope ratios using cavity ring-down spectroscopy: determination of 13C/12C for carbon dioxide in human breath,” Anal. Chem.. 74, , 2003 –2007 (2002). 0003-2700 CrossRef
Waclawek  J. P., , Moser  H., and Lendl  B., “Compact quantum cascade laser based quartz-enhanced photoacoustic spectroscopy sensor system for detection of carbon disulfide,” Opt. Express. 24, , 6559 –6571 (2016). 1094-4087 CrossRef
Yufei  M.  et al., “Sensitive detection of carbon monoxide based on a QEPAS sensor with a 2.3  μm fiber-coupled antimonide diode laser,” J. Opt.. 17, , 055401  (2015). 0150-536X CrossRef
Wojtas  J., “Application of cavity enhanced absorption spectroscopy to the detection of nitric oxide, carbonyl sulphide, and ethane-breath biomarkers of serious diseases,” Sensors. 15, , 14356 –14369 (2015). 0746-9462 CrossRef
Li  C.  et al., “Compact TDLAS based optical sensor for ppb-level ethane detection by use of a 3.34  μm room-temperature CW interband cascade laser,” Sens. Actuators B. 232, , 188 –194 (2016). 0925-4005 CrossRef
Ivascu  I. R.  et al., “CO2 laser photoacoustic measurements of ethanol absorption coefficients within infrared region of 9.2-10.8  μm,” Spectrochim. Acta Part A. 163, , 115 –119 (2016).CrossRef
Zhou  S., , Han  Y., and Li  B., “Simultaneous detection of ethanol, ether and acetone by mid-infrared cavity ring-down spectroscopy at 3.8  μm,” Appl. Phys. B. 122, , 187  (2016).CrossRef
Navas  M., , Jiménez  A., and Asuero  A., “Human biomarkers in breath by photoacoustic spectroscopy,” Clin. Chim. Acta. 413, , 1171 –1178 (2012). 0009-8981 CrossRef
Cernat  R.  et al., “Laser photoacoustic spectroscopy method for measurements of trace gas concentration from human breath,” Rom. Rep. Phys.. 62, , 610 –616 (2010). 1221-1451 
Wang  Z., , Geng  J., and Ren  W., “Quartz-enhanced photoacoustic spectroscopy (QEPAS) detection of the ν 7 band of ethylene at low pressure with CO2 interference analysis,” Appl. Spectrosc.. 1 –8 (2017). 0003-7028 CrossRef
Ren  W., , Luo  L., and Tittel  F. K., “Sensitive detection of formaldehyde using an interband cascade laser near 3.6  μm,” Sens. Actuators, B. 221, , 1062 –1068 (2015).CrossRef
Szabo  A.  et al., “In situ and wide range quantification of hydrogen sulfide in industrial gases by means of photoacoustic spectroscopy,” Meas. Sci. Technol.. 24, , 065501  (2013).CrossRef
Viciani  S.  et al., “A quartz-enhanced photoacoustic sensor for H2S trace-gas detection at 2.6  μm,” Appl. Phys. B. 119, , 21 –27 (2015).CrossRef</