# An all-Optical Photoacoustic Sensor for the Detection of Trace Gas

Thomas Lauwers, Alain Glière, Skandar Basrour

► To cite this version:

Thomas Lauwers, Alain Glière, Skandar Basrour. An all-Optical Photoacoustic Sensor for the Detection of Trace Gas. Sensors, 2020, 20 (14), 10.3390/s20143967. [https://doi.org/10.33383/s20143967](https://doi.org/10.33383/s20143967)

HAL Id: hal-03338389

https://hal.science/hal-03338389v1

Submitted on 26 Sep 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# An all-Optical Photoacoustic Sensor for the Detection of Trace Gas

Thomas Lauwers <fnref n="1" />*, Alain Glière <fnref n="1" /> and Skandar Basrour <fnref n="2" />

Received: 23 June 2020; Accepted: 13 July 2020; Published: 16 July 2020

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** A highly sensitive Fabry-Perot based transduction method is proposed as an all-optical alternative for the detection of trace gas by the photoacoustic spectroscopy technique. A lumped element model is firstly devised to help design the whole system and is successfully compared to finite element method simulations. The fabricated Fabry-Perot microphone consists in a hinged cantilever based diaphragm, processed by laser cutting, and directly assembled at the tip of an optical fiber. We find a high acoustic sensitivity of 630 mV/Pa and a state-of-the-art noise equivalent pressure, as low as ~ 2 µPa/√Hz at resonance. For photoacoustic trace gas detection, the Fabry-Perot microphone is further embedded in a cylindrical multipass cell and shows an ultimate detection limit of 15 ppb of NO in nitrogen. The proposed optical trace gas sensor offers the advantages of high sensitivity and easy assembling, as well as the possibility of remote detection.

**Keywords:** Extrinsic Fabry-Perot interferometer; photoacoustic sensor; trace gas sensor; lumped model; all-optical sensor; all-optical PA sensor; multipass PA sensor

## 1. Introduction

The photoacoustic (PA) effect consists in the generation of an acoustic wave originating from the specific absorption of laser radiation by the medium of interest. Based on this phenomenon, the photoacoustic spectroscopy technique enables high sensitivity probing [1,2] and sensing [3,4]. The targeted molecules absorb a part of the optical energy emitted from the laser and release it by thermal relaxation. The modulation of the optical power or wavelength creates a variation of the heat produced, acting as a harmonic source for the photoacoustic pressure wave, whose amplitude indicates the targeted molecules concentration. For gas sensing, the mid-infrared range, often called the molecule fingerprint region, is a very interesting zone for photoacoustic spectroscopy, where most molecules of interest (polluting or toxic gases, biomarkers, etc.) have a specific spectral signature.

A generic trace gas photoacoustic sensor can be composed of a laser, a closed volume, called the cell, where the acoustic wave is generated, and a transducer for the detection of the pressure variation. Most systems use electromechanical transductions. Among them, the standard condenser microphone or the piezoelectric quartz tuning fork are known for their robustness, high stability, low-cost and miniaturization capabilities. Due to the low pressure level generated by PA effect, the signal is often amplified by means of acoustic [5], mechanical, [6,7], or optical [8,9] resonance.

Optical readout of the acoustic wave can be an interesting alternative to electromechanical transductions for PA detection in harsh environment, where high temperatures or electromagnetic noise can severely damage the transducer or decrease its sensitivity. A wide variety of optical transductions have already been proposed in the literature, ranging from the simple position sensor [10] to interferometric readout systems, like Mach Zehnder [11] or free space or fibered Michelson

<footnote n="1">¹ Université Grenoble Alpes, CEA, LETI, F38000 Grenoble, France; alain.gliere@cea.fr</footnote>

<footnote n="2">² Université Grenoble Alpes, CNRS, Grenoble INP, TIMA, 38000 Grenoble, France; skandar.basrour@univ-grenoble-alpes.fr</footnote>

<footnote>* Correspondence: thomas.lauwers@cea.fr</footnote>
configurations [12,13]. However, those systems are bulky and require either a precise alignment of the optical elements or long optical fibers, causing a loss of robustness. Among the interferometric transducers, a variant of the traditional Fabry-Perot (FP) cavity, called fiber-optic extrinsic Fabry-Perot interferometer (EFPI) [14,15], can settle the size problem by using the air/fiber interface as first mirror and a reflective diaphragm as second mirror. This low finesse optical cavity, requiring no specific alignment, enables remote detection with optical fiber [16]. Different materials and geometries have been proposed, involving mainly silver [17], graphene [18] or parylene [19] circular diaphragms, or more recently stainless steel cantilevers [20]. The latter are more sensitive, due to the enhanced displacement obtained at the free end. The photoacoustic signal in such systems is mainly amplified by means of mechanical resonances, sometimes coupled with optical amplification, as has been shown in the recent work by Zhang et al. [21], with the use of Heriot cells.

In this paper, we develop a lumped element model to assist in the design of an EFPI based photoacoustic sensor. Compared to the finite element methods (FEM), this analytical model provides much faster problem solving with good predictive capability. The second part of the article is focused on the experimental implementation of such a sensor. We propose a highly sensitive all-optical PA sensor based on an EFPI transducer. The latter builds on a previous work on the resonant hinged cantilever [22], with enhanced sensitivity through an improved laser cutting process and the use of a thicker structure. To create the PA sensor, the transducer is embedded in a cylindrical multipass cell, thus increasing the interaction length between the light and the gas of interest. We finally present a complete characterization of our all-optical PA sensor prototype and assess the sensor performance on a weak NO absorption line.

## 2. Design of the Photoacoustic Sensor

### 2.1. *Description of the System*

The all-optical photoacoustic sensor proposed here (Figure 1) is composed of a cylindrical cell filled with gas, illuminated (red light beam) by a mid-infrared quantum cascade laser (QCL) to generate the acoustic wave (lower, front volume). The pressure wave generated is measured by the fiber optic Fabry-Perot acoustic transducer (upper, back volume), with a second laser diode at 1.55 µm (green light beam) used to probe the position of a reflective diaphragm sensitive to acoustics.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Diagram (not to scale) of the all-optical multipass PA trace gas sensor with the PA cell (lower volume), the transducer part with the optical cavity (upper volume), the probe and excitation laser (respectively green and red light beam).</figcaption></figure>

The transducer (Figure 2a) is composed of two mirrors. The first one ($M_1$), with low reflectivity, is constituted by the air/fiber interface and the second one ($M_2$) is a reflective mirror, designed to be highly sensitive to acoustic perturbations. The light emitted by a 1.55 µm laser diode is injected in the FP cavity through an optical fiber and the reflected light is measured by a photodetector.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** (a) Diagram of the fiber-optic Fabry–Perot based transducer; (b) interferometric pattern of the reflected optical signal with quadrature point Q.</figcaption></figure>

The reflectivity of the FP cavity $R_{FP}$ can be analytically calculated considering a two light—beam interferometer:

$$R_{FP}(\lambda, L_{cav}) = R_1 + R_2 + 2\sqrt{R_1 R_2} \cos\left(\frac{4\pi L_{cav} n}{\lambda}\right). \quad (1)$$

Here $R_1$ and $R_2$ correspond respectively to the reflectivity of the air/fiber interface and the reflective diaphragm, $L_{cav}$ is the length of the FP cavity, $n$ is the refractive index of the medium and $\lambda$ is the wavelength of the probe laser. With the assumption that the laser wavelength $\lambda$ is stable and that the refractive index variations are slow with respect to the acoustic frequency, the optical power reflected back from the Fabry–Perot cavity oscillates according to the variation $\delta L_{cav}(p)$ of the length of the optical cavity:

$$\delta I_r(p) = I_i \frac{dR_{FP}}{dL} \delta L_{cav}(p) \sin(\omega_{ac} t), \quad (2)$$

where $I_i$ is the incident power of the laser, $p$ the amplitude of the photoacoustic pressure wave and $\omega_{ac}$ the acoustic angular frequency. To obtain a linear relation between the displacement and the reflected power with high sensitivity, the probe laser wavelength is set to the quadrature point (point Q on Figure 2b). The displacement $\delta L_{cav}$ induced by the pressure wave depends on the mechanical structure and the photoacoustic cell geometry. Due to the cylindrical symmetry of the fiber, the diaphragm is a circular plate clamped on its circumference. To benefit from a large displacement of the structure, it has been decided to embed a cantilever structure in the plate (Figure 3), so that the free end of the cantilever is situated at the center, where the probe laser thus measures the maximum deflection (red spot on Figure 3).

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Diagram of the cantilever-based diaphragm (**left**) and the hinged cantilever-based diaphragm (**right**). Red spots indicate the probe laser location.</figcaption></figure>

The problem arising for this particular structure is the coupling between the unwanted damped plate resonance and the sensitive cantilever resonance. In order to decouple these mechanical modes we propose a hinge geometry (Figure 3, right) detailed in the next section. Moreover, this structure allows us to obtain a resonant frequency in the kHz range while using a relatively thick structure (100 µm) which is easy to fabricate and assemble.
## 2.2. Mechanical Model

The mechanical structure is designed to resonate at a few kHz, well below the cut-off frequency induced by the relaxation time of the molecules of interest (few tens of kHz for simple molecules, such as NO, in air [23]). In order to simplify the study of the diaphragm's dynamics, it is possible to isolate the clamped-free cantilever and the clamped circular plate fundamental resonance modes, respectively calculated with the Euler beam and Kirchhoff-Love plate theories:

$$ \omega_{m,c} = \alpha_c \frac{t_m}{l_1^2} \sqrt{\frac{E}{12\rho}} \quad \text{and} \quad \omega_{m,p} = \alpha_p \frac{t_m}{R^2} \sqrt{\frac{E}{12\rho(1-v^2)}} \quad (3) $$

Here $\omega_{m,c}$ (respectively $\omega_{m,p}$) refers to the first mechanical resonance angular frequency of a clamped-free cantilever (respectively a clamped circular plate), $\alpha_c = 3.52$ and $\alpha_p = 10.22$ are the coefficients corresponding to the fundamental modes of resonance, $t_m$ is the thickness of the structure, $l_1$ the length of the cantilever, $R$ the radius of the plate, $E_m$ the Young modulus of the material, $\rho_m$ its mass density and $v_m$ its Poisson's coefficient.

The hinged cantilever resonant angular-frequency $\omega_{m,hc}$ can also be calculated, by the Rayleigh's quotient method [24]:

$$ \omega_{m,hc} = \sqrt{\frac{w_2[l_2(3l_1^2 + 3l_1l_2 + l_2^2)w_1 + l_1^3w_2]}{B}} t_m \sqrt{\frac{E_m}{\rho_m}}, \quad (4) $$

where $w_1, l_1$ refer to the width and length of the cantilever, $w_2, l_2$ to the width and length of the hinge, and $B$ is a polynomial function of the cantilever geometry (for more details, see Appendix A).

One way to ensure that the first bending mode of the cantilever of the coupled structure will not be disturbed by the fundamental plate mode is to ensure that $\omega_{m,p} \gg \omega_{m,c}$. In this way, only the higher cantilever mode orders will be coupled with the plate. With this condition fulfilled, we also expect a correct prediction of the isolated case on the final system dynamics. To identify a suitable cantilever length $l_1$ we calculate the resonance frequency of the three structures and we introduce the reduced frequency $x$ defined as the ratio $x = \omega_{m,c}/\omega_{m,p}$.

For the calculations, a plate thickness $t_m = 100$ µm with a radius $R = 4$ mm are fixed to obtain a high plate resonance frequency and to meet with the integration constraints of the diaphragm inside the photoacoustic cell. Then we chose a cantilever width $w_1 = 1$ mm and the hinge parameters $w_2 = 100$ µm and $l_2 = 300$ µm. This set of parameters allows to obtain a flexible hinge and gives a large flexibility for the choice of $l_1$ while remaining feasible by the diaphragm fabrication process. The material parameters are $E_m = 185$ GPa and $\rho_m = 7850$ kg/m³. Figure 4 compares the reduced frequency for the case of the normal clamped-free (CF) cantilever (blue curve) and the hinged CF cantilever (red curve). The hinged cantilever resonance frequency calculated with the Rayleigh's quotient method is very close to the finite element method (FEM) results (yellow curve) performed with a commercial software (COMSOL Multiphysics, COMSOL AB, Sweden). We clearly see on the figure that adding a hinge helps to lower the cantilever resonance frequency and thus to lie in the zone where $x < 1$. For a length of $l_1 = 2$ mm, the analytical calculation gives $x_{\text{cantilever}} \sim 1$ while $x_{\text{hinged-cantilever}} \sim 0.3 < 1$, enabling a weak coupling between the hinged cantilever and the plate mode.

This result is confirmed by taking in account the complete mechanical structure with a FEM eigenmode study, presented on the right side of Figure 4. As expected, the eigenmode resolution for the simple cantilever based diaphragm (first column) shows a first cantilever mode (19.4 kHz) that occurs after the fundamental plate mode ($x > 1$), while for the hinged cantilever based diaphragm (second column) the first cantilever mode (5.1 kHz) is well below the fundamental plate mode ($x \sim 0.3$).
<figure><img src="image_7.png" /><figcaption>**Figure 4.** On the left, reduced frequency of the isolated cantilever beam (blue) and hinged cantilever beam (analytical model in red and FEM result with yellow squares) as a function of the length $l_1$. On the right, first two eigenmodes of the cantilever-based diaphragm (first column) hinged cantilever-based diaphragm (second column).</figcaption></figure>

However, the cantilever mode of the complete structure differs notably from the isolated cantilever mode: this is because the cantilever is not perfectly clamped at the hinge end, which is situated at around 2.5 mm from the clamped circumference. The analytical calculation of the isolated hinged cantilever resonance frequency gives $f_{m,hc} = 6.7$ kHz which is above the frequency obtained for the complete structure. To match the isolated case resonance frequency and the one obtained with the FEM resolution on the complete structure, an effective hinge length $l_{2,eff}$ is introduced in the mechanical resonance formula from Equation (4), which is found to be 480 µm instead of 300 µm.

## 2.3. Photoacoustic Model

To provide a complete description of the PA sensor, the next section focuses on the photoacoustic generation and the acoustic coupling between the considered medium (ambient air or nitrogen in our case) and the mechanical structure. In the following sections, the fluid is considered as an ideal gas with equilibrium pressure, temperature and velocity $P_0$, $T_0$, $v_0$ and described by small variations of the pressure field $p$, temperature $\tau$ and velocity $v$, which are a solution of the fully linearized Navier-Stokes equations (FLNS).

### 2.3.1. Heat Source

The generation of the photoacoustic wave originates from a periodic heat source due to the specific absorption of laser radiation. Using the energy equation of the FLNS model, the temperature fluctuation inside the PA cell at the position $r = (r, \theta, z)$ and angular frequency $\omega$ can be expressed as:

$$ \tau(r, \omega) = (1 - F(r, \omega)) \frac{1}{\rho_0 c_p} \left[ p + \frac{Q}{\dot{j}\omega} \right], \quad (5) $$

where $Q$ refers to the deposited power density in W·m⁻³, $\rho_0$ the air density and $c_p$ the heat capacity at constant pressure. $F$ is a function taking into account the thermal loss originating from conduction [25] and is a solution of the Helmholtz equation. The temperature expression obtained in Equation (5) can be introduced in the FLNS continuity equation in the integral form:

$$ \int_V \left( \dot{j} \omega \frac{p}{P_0} \right) dV = \int_{\partial V} \mathbf{v} \cdot (-\mathbf{n}) dS + \int_V \left( \dot{j} \omega \frac{\tau}{T_0} \right) dV, \quad (6) $$
where V indicates the volume of the cell and n is the outward pointing unit vector, normal to the volume boundary ∂V. If we consider a closed cell with volume V<sub>fv</sub> that is small enough to have a uniform pressure distribution, the Equation (6) leads after integration to:

$$p = i_{PA} Z_{fv}(\omega) \text{ with } Z_{fv}(\omega) = \frac{1}{\int_{\omega}^{\omega+\langle F \rangle} \int_{\langle F \rangle}^{\langle F \rangle+\langle F \rangle} \int_{\langle F \rangle+\langle F \rangle}^{\langle F \rangle+\langle F \rangle+\langle F \rangle} \int_{\langle F \rangle+\langle F \rangle+\langle F \rangle}^{\langle F \rangle+\langle F \rangle}} \text{ and } i_{PA} = (\gamma - 1)C_{fv}Q, \quad (7)$$

where γ is the heat capacity ratio, C<sub>fv</sub> = V<sub>fv</sub>/(γP<sub>0</sub>) is called acoustic compliance of the cell, ⟨F⟩ depends on ω and refers to the volume averaged value of F(r, ω). This equation can be represented by an electric equivalent model, where the heat source Q is responsible for gas expansion and acts as a current source i<sub>PA</sub> (m<sup>3</sup>·s<sup>-1</sup>), while the pressure increase p depends on the impedance Z<sub>fv</sub>(ω). The analytical expression of Z<sub>fv</sub>(ω) is not straightforward because it is expressed with the ⟨F⟩ function which cannot be simply solved in the case of a cylindrical volume. The approach used here consists in considering two asymptotic cases of the cylindrical cell of radius R<sub>cell</sub> and height h<sub>cell</sub>. The thermal loss occurring on the cylinder sidewalls (denoted R<sub>r</sub>) is obtained by considering an infinitely long cylinder (R<sub>cell</sub> ≪ h<sub>cell</sub>) and by evaluating the real part of lim<sub>ω→0</sub> Z<sub>fv</sub>(ω). In a similar way, the loss occurring on the cylinder ends (denoted R<sub>z</sub>) is obtained by considering an infinitely flat cylinder (R<sub>cell</sub> ≫ h<sub>cell</sub>). The analytical expressions of R<sub>r</sub> and R<sub>z</sub> are given in Appendix B. Combining the two asymptotic cases gives a final thermal resistance R<sub>th</sub>, corresponding to R<sub>r</sub> in parallel with R<sub>z</sub>. The cell impedance can then be written in term of the thermal resistance in parallel with the acoustic compliance C<sub>fv</sub>:

$$R_{th} = \frac{\rho_0 c_p P_0}{8 \kappa \pi} \frac{h_{cell}}{h_{cell}^2 + 3/2 R_{cell}^2} \text{ and } C_{fv} = \frac{V_{fv}}{\gamma P_0}. \quad (8)$$

The photoacoustic pressure generation thus behaves like a low-pass filter, with a cut-off angular frequency ω<sub>h</sub> = 1/R<sub>th</sub>C<sub>fv</sub>. This effect is due to the coupled effect of the volume and the thermal resistance taking into account the losses due to thermal conduction to the walls. The Figure 5 shows the equivalent circuit with the current source and a comparison between the lumped model and a FEM resolution of the FLNS equation in the case of a cylinder with a radius R<sub>cell</sub> = 10 mm and a height h<sub>cell</sub> = 10 mm. These dimensions correspond to the size of the already available PA cell used for the experimental implementation.

<figure><img src="image_8.png" /><figcaption>**Figure 5.** On the left, equivalent electrical circuit responsible for the pressure generation in the closed front volume p<sub>fv</sub>. On the right, transfer function of the circuit (generated pressure per deposited power) and comparison with the FEM simulation using the thermoacoustic model.</figcaption></figure>

With the lumped model we can calculate the transfer function corresponding to the increase of pressure for a given deposited power in Pa/W. We observe a capacitive behavior of the pressure inside the cell volume p<sub>fv</sub> at the working frequencies from 1 kHz to 10 kHz, while a plateau is noticeable
at very low frequencies below 1 Hz, due to the thermal losses $R_{th}$. Both asymptotic behaviors are correctly described by the analytical lumped model.

### 2.3.2. Acoustic Coupling

Once the photoacoustic transfer function is defined the behavior of the fluid surrounding the diaphragm must be described. The problem is divided in three parts: the external environment, the ventilation slit, which corresponds to the gap between the cantilever frame and the rest of the plate, and the back volume constituting the FP cavity (Figure 6a). The diaphragm is supposed to be rigid (no displacement). By reducing the three FLNS equations to a one-dimensional slit problem we can find an analytical solution of the pressure propagation as a solution of the Helmholtz equation [26]:

$$ \frac{\partial^2 p}{\partial z^2} + k_l^2 p = 0 \qquad (9) $$

where $k_l$ is a thermoviscous wave vector that depends on the air parameters and slit geometry. This harmonic behavior can be formulated with a lumped element model, where the differential pressure across the slit is described with a thermoviscous impedance $Z_v$. In a similar way to the previous section, $Z_v$ is defined as the ratio between the differential pressure between both parts of the slit $\Delta p = p_{bv} - p_{ext}$ and the fluid velocity flow through the slit, normal to the diaphragm surface, denoted $i = v S$ (see Figure 6a), where S is the ventilation slit area. In the slit configuration, $Z_v$ is mainly dominated by its real part, which represents the viscous resistance $R_v$, that takes in account the viscous losses experienced by the fluid through the thin slit. On the other hand, the volume of the back chamber introduces an acoustic compliance $C_{bv}$ (see Section 2.3.1):

$$ R_v = 12\mu_0 \frac{t}{w_v^3 l_c} \quad \text{and} \quad C_{bv} = \frac{V_{bv}}{\gamma P_0}. \qquad (10) $$

<figure><img src="image_9.png" /><figcaption>Figure 6. (a) Sectional schematic of the slit problem; (b) Pressure behavior inside the slit and the front volume for a slit length of 100 µm with $p_{ext} = 1$ Pa. (the mechanical structure and the back volume are not shown for the sake of clarity).</figcaption></figure>

Here $\mu_0$ and $P_0$ correspond to the viscosity and atmospheric pressure, and $V_{bv}$ is the volume of the cavity. The cantilever slit is assimilated to a rectangular slit with a width $w_v$ and a length $l_c$ equal to the perimeter of the hinged cantilever: $l_c = 2(l_1 + l_2) + 2w_1 - w_2$.

A FEM resolution of the FLNS equation is performed to confirm the results obtained with the lumped model for the acoustic response of the transducer. Here, the external environment consists in a harmonic plane wave boundary condition, the mechanical structure is assimilated to a rigid wall and the back volume has isothermal walls boundary conditions. The slit has the shape of the hinged cantilever, a height of 100 µm and a width of 35 µm. The volume of the slit surrounding the hinged
cantilever is reported on Figure 6b. We can note the pressure variation along the height of the slit from 0.1 Pa in the back volume up to 1.1 Pa in the front volume.

Coupled to the back volume, this impedance induces a high-pass filter, meaning that the fluid leaks through the slit at low frequencies and thus slowly equilibrates the differential pressure between both sides of the diaphragm.

As can be seen from the Figure 7, a good agreement is obtained between the FEM simulations and the lumped element model. The acoustic response of the system shows that up to a certain cut-off frequency $f_v$, the gas leak through the slit diminishes, leading to an increase of the differential pressure and thus the measurable signal. Therefore, the slit should be chosen small enough to give a cut-off frequency $f_v$ situated below the mechanical resonance frequency of the system, denoted $f_m$.

<figure><img src="image_10.png" /><figcaption>**Figure 7.** On the left, equivalent electrical circuit responsible for the pressure differential across the slit, given an applied pressure $p_{ext}$. On the right, transfer function of the circuit and comparison with the FEM simulation using the FLNS model.</figcaption></figure>

## 2.4. Global Resolution of the System

The complete lumped elements model, represented on Figure 8, can be used to describe quantitatively the harmonic response of the photoacoustic sensor.

<figure><img src="image_11.png" /><figcaption>**Figure 8.** Complete equivalent electrical circuit of the PA sensor with (from left to right) the current source and the front volume, the acoustic slit with a velocity flow $i_v$ and the mechanical diaphragm with a velocity flow $i_m$ and, finally, the back volume.</figcaption></figure>

The mechanical parameters $C_m, L_m$ have been calculated with the expression of the stiffness and mass (see Appendix A) and the effective hinge length calculated in Section 2.2:

$$ \left\{ \begin{array}{l} C_m = S_{eff,m}^2 / k_{eff,m} \\ L_m = m_{eff,m} / S_{eff,m}^2 \end{array} \right. , \qquad (11) $$
where $S_{\text{eff,m}}$ is the effective area of the cantilever given in Appendix A. The thermoacoustic parameters $R_{\text{th}}$, $C_{\text{fv}}$, $R_{\text{v}}$, $C_{\text{bv}}$ have been calculated with the formula from Section 2.3. The fluid flow induced by the heat source $i_{\text{PA}}$ is calculated for a deposited power of 1 mW, and the cantilever displacement $\Delta z$ is given by the potential difference $\Delta p_C$ across the mechanical compliance $C_m$, as $\Delta p_C = i_m / (jC_m\omega) = (\Delta zS_{\text{eff,m}})/C_m$.

On the other hand, a FEM simulation taking into account the precise geometry of the system can be implemented. The front and back volumes are surrounded with isothermal walls and are separated from each other by the mechanical structure, modelled by a 3D plate. Finally, the laser excitation is modelled by a uniform source term distributed in the entire front volume. The two color maps on Figure 9 indicate the pressure field amplitude in the front and back volume, and the displacement amplitude of the cantilever. Let us note that the small volume approximation made in the previous section to establish the lumped element model is verified here, since the pressure fields in both volumes are uniform.

<figure><img src="image_12.png" /><figcaption>**Figure 9.** Results of the FEM simulations for the complete PA sensor near the resonance frequency for a deposited power of 1 mW. The left color bar indicates the pressure value in the front volume (dark red) and back volume (blue), while the right color bar indicates the displacement of the cantilever.</figcaption></figure>

Figure 10 shows the transfer function obtained by FEM resolution of the FLNS model for the different blocks studied: the blue curve corresponds to the photoacoustic generation (Section 2.3.1), the red curve to the acoustic transfer of the transducer (Section 2.3.2) and the yellow curve shows the mechanical transfer function of the diaphragm (Section 2.2). Finally, the FEM resolution taking into account the fluidic and mechanical coupling in purple is expressed in nm/W (cantilever displacement per deposited unit of power) and compared to the transfer obtained by analytical resolution of the lumped elements model (dotted black). A very good agreement between the analytical model and the FEM simulation is obtained.

Based on the modelling and simulation results, a set of rules can be established for further designs. Firstly, it is advantageous to consider a small ventilation slit width to obtain the condition $\omega_v < \omega_m$. Secondly, from Section 2.3.1, it appears that a reduced height $h_{\text{cell}}$ leads to a decrease of the energy lost by conduction mechanism through the walls ($R_{\text{th}}$), meaning a higher value of $\omega_h$, and thus an increase of the generated photoacoustic pressure at the working frequency $f_m$. Finally, a large cylinder radius is preferred, in order to increase the gas-laser light interaction and consequently the deposited power in the PA cell.
<figure><img src="image_13.png" /><figcaption>**Figure 10.** Calculated transfer function of the different blocks of the PA sensor and complete transfer in nm/W with FEM and comparison with the lumped element model (dotted black curve).</figcaption></figure>

# 3. Fabrication and Experimental Set-Up

## 3.1. Acoustic Transducer

The diaphragm is fabricated from a 100 µm thick stainless-steel foil with laser micromachining process. A water jet laser technology (Synova SA, Duillier, Switzerland) has been employed to obtain flat structures with reduced residual stress (Figure 11a). The minimal cutting size is around 35 µm. This sets the ventilation slit width around the cantilever and induces a cut-off frequency around $f_v = 1.3$ kHz.

<figure><img src="image_14.png" /><figcaption>**Figure 11.** (a) Microscope image of the manufactured hinged cantilever; (b) 3D view of the fabricated PA sensor with (from left to right) the optical fiber, a threaded adapter, the FP transducer support, the cantilever-based diaphragm, a threaded ring to hold the diaphragm and the cylindrical multipass PA cell.</figcaption></figure>

The acoustic transducer part is made with an aluminum support, which is simply screwed onto a fiber (Figure 11b). To enable fast prototyping and test of different mechanical resonator geometries the diaphragm is aligned by two posts and held on the support with a screwed ring to obtain a circularly clamped condition.

The last element of the sensor is the photoacoustic cell, which is fixed to the acoustic transducer. A cylindrical metallic cell available “off-the-shelf” was used, with a 10 mm radius and a height of 10 mm. Its walls have been polished in order to achieve maximum reflections and thus maximize the length of the laser interaction with the gas. Finally, a tapered aperture is drilled in the cylinder to let the excitation laser light enter without exit.
## 3.2. Experimental Bench Set-Up

The all-optical PA sensor has been used to detect traces of NO in N₂, with the following experimental set-up (Figure 12). The probe DFB laser (Eblana photonics, Ireland), emitting at 1.55 µm, is launched into the fiber optic FP cavity. The reflected light is retrieved with an optical circulator and measured by a photodiode with a high frequency response (500 kHz). The pump laser is a Quantum Cascade Laser (QD5250CM1, Thorlabs Inc, Newton, NJ, USA) with a power of 100 mW emitting in the mid-infrared at 5.2 µm, where a weak NO absorption peak lies (absorbance of 1.7 × 10⁻⁵ cm⁻¹ for 1 ppm of NO). The QCL is stabilized in temperature and the current is controlled for a slow wavelength scan across the NO absorption peak and a wavelength modulation scheme [27]. The NO is injected in the PA sensor through 1/16" tubes from a calibrated gas cylinder containing 25 ppm of NO in N₂. The injection and dilution are performed with a gas diluter (GasMix, Alytech, Juvisy-sur-Orge, France).

<figure><img src="image_15.png" /><figcaption>**Figure 12.** Diagram of the photoacoustic measurement setup. Optical fibers are drawn in green and electric connections in black.</figcaption></figure>

The wavelength modulation is performed with a lock-in amplifier (HF2LI, Zurich Instrument, Switzerland). The AC part of the signal retrieved from the photodetector is demodulated at the same frequency for a 1f detection or at twice the modulation frequency for a 2f detection. While the 1f signal has the greatest amplitude and is used to determine the location of an absorption peak, the 2f detection method is preferred for gas calibration as it allows eliminating the background signal.

Due to the high sensitivity of the FP cavity, the location of the quadrature working point can be perturbed by environmental parameters like temperature, static pressure or humidity. This leads to a drift of the working point and as a consequence a time variation of the sensitivity. In order to stabilize this drift, a proportional integral (PI) correction is implemented to control the emission wavelength of the probe laser and to keep the DC optical signal at its quadrature value. The tuning of the wavelength is directly obtained by controlling the drive current of the laser. This eliminates the need for an expensive and bulky piezo controlled external cavity diode laser, but results in unwanted variation of the emitted power, which is also a function of the drive current. In order to remain insensitive to the emitted power variations, the DC signal, on which the PI correction is made, is normalized with respect to the power of the probe laser. If the drift reaches the upper (or lower) limit of the laser wavelength tunability, a new quadrature point, lying in the laser tuning range, is recycled [28]. This closed loop operation helps to obtain a long-term stability of the system, as well as high sensitivity and immunity to low frequency external perturbations.
# 4. Experimental Results

The complete characterization of the sensor has been performed in two steps. Firstly, the acoustic sensitivity of the FP transducer alone has been measured by using it as a microphone, placed in front of a calibrated speaker. Secondly, the FP transducer has been mounted on a multipass photoacoustic cell and, with the help of a commercial gas diluter, further characterized (resonance frequency, NO limit of detection and Allan deviation).

## 4.1. Microphone Characterization

During this first stage of characterization, the acoustic transducer is placed in an insulation box, together with a calibrated speaker, to assess the microphone operation mode. Before each measurement, an initialization step is carried out as follows to determine the optimal working point. A sweep on the drive current of the probe laser is made to scan a few periods of the Fabry-Perot interferences. After a normalization by the emitted power, the signal is fitted to a sinusoid in order to retrieve the phase shift, the interfringe period and fringe amplitude and offset. The quadrature working point can then be calculated and is used to lock the DC signal, using a PI correction loop. Figure 13a shows the retrieved reflectivity of the cavity and the sinusoidal fit using the analytical FP reflectivity of Equation (1). In abscissa, the wavelength shift is relative to the working wavelength of 1.55 µm. The reflectivity of the first mirror R₁ corresponds to the air/fiber interface, while the second reflectivity R₂ takes into account not only the reflectivity of the cantilever but also the overall losses of this low-finesse cavity (divergence of the fiber output, misalignment of the diaphragm, etc.).

<figure><img src="image_16.png" /><figcaption>Figure 13. (a) Interference pattern of the transducer FP cavity as a function of the relative wavelength shift $\Delta\lambda$ around 1.55 µm; (b) Acoustic calibration of the FP transducer at the hinged cantilever resonance frequency (5.1 kHz).</figcaption></figure>

Figure 13b shows the acoustic sensitivity measured at the resonance frequency of the hinged cantilever (5.1 kHz). Each point is calculated as the mean value for a time duration of 10 s and is associated with its standard deviation (too small to be represented on the figure). We observe a linear acoustic response and a high sensitivity of 630 mV/Pa. Given the standard deviation obtained when no acoustic signal is applied and a measurement bandwidth of 1 Hz, a LOD of ~2µPa/√Hz is obtained. The LOD is probably overestimated as it corresponds in fact to the white noise present in the experiment room, which is not completely attenuated by the homemade acoustic insulation box used.

## 4.2. Photoacoustic Measurements

The FP transducer has then been mounted on a multipass cell for further characterization in photoacoustic detection mode. To identify the resonance frequency $f_{res}$ in the closed PA cell,
the response of the transducer is probed on a strong atmospheric water absorption line at 1906 cm⁻¹. To do this, the wavelength of the QCL is modulated around the absorption peak. The wavelength modulation is realized via the drive current of the QCL, with an amplitude of 10 mA corresponding to a modulation depth 0.16 cm⁻¹, and a frequency varying from 1 to 10 kHz. The normalized response (Figure 14) shows an excellent agreement with the FEM model with a resonant frequency of 5.16 kHz and a quality factor of 70 in the closed photoacoustic cell.

<figure><img src="image_17.png" /><figcaption>**Figure 14.** Normalized frequency response of the PA sensor in presence of atmospheric water at 1906 cm⁻¹.</figcaption></figure>

Once the exact resonance frequency is known, a second harmonic detection (2f) strategy is used to identify the NO absorption peak [27]. In the 2f wavelength scheme, the modulation frequency is at $f_{res}/2$ and the signal demodulated at $f_{res}$. As shown on the Figure 15a, the 2f signal is maximal when the QCL wavelength matches the maximum of the NO absorption line at 1906.15 cm⁻¹, which corresponds to a laser drive current of 531 mA. The previously determined experimental parameters (QCL drive current, wavelength modulation depth and frequency) are kept and the time constant of the lock-in is fixed at 1 s to perform a gas concentration calibration measurement. A gas mixture (NO in N₂) is injected through with a total flow of 20 mL/min and different dilution ratios to make concentrations steps. The measurement lasts around 2 h, with a stable signal of the acoustic transducer (Figure 15b).

<figure><img src="image_18.png" /><figcaption>**Figure 15.** (a) 2f scan of the NO absorption line (blue) compared with the Hitran database (orange); (b) Calibration of the PA sensor for a measurement of 6600 s.</figcaption></figure>
The step-by-step concentration measurement gives a linear sensor response and a sensitivity of 4.7 µV/ppm with a background noise of 11.7 µV. The noise measurement leads to a limit of detection (LOD) of 410 ppb/√Hz and a normalized noise equivalent absorption (NNEA) of 5.1 × 10⁻⁷ W·cm⁻¹·Hz⁻¹/². The Allan deviation measurement (Figure 16) confirms a stability up to 58 s, which corresponds to a noise equivalent bandwidth of 1.3 mHz and thus an ultimate LOD of 15 ppb.

<figure><img src="image_19.png" /><figcaption>**Figure 16.** Allan deviation of the PA sensor in presence of a 5 ppm of NO in N₂ (blue curve) and theoretical white noise curve in 1/√t (dotted black).</figcaption></figure>

**5. Conclusions**

An all-optical multipass photoacoustic trace gas sensor, based on a Fabry-Perot hinged cantilever transducer, allowing remote detection and offering high sensitivity, has been presented.

As a first step, we proposed a complete analytical lumped elements model to comprehensively describe the sensor's frequency behavior and performances. With respect to the complete and lengthy FEM resolution, the model is much faster and provides however correct predictions. In addition, it can be easily transposed to different types and geometries of PA transducers.

Then, we presented the fabrication process and the experimental results obtained on our Fabry-Perot acoustic transducer, characterized in a microphone configuration. Optical readout of the hinged cantilever displacement was performed with a 1.55 µm DFB laser. The mechanical structure was micromachined by laser cutting from a 100 µm stainless steel foil and assembled on a metallic support to create the FP cavity. The FP transducer has an optical cavity length of 1.6 mm and a mechanical resonance frequency at 5.1 kHz, for an acoustic LOD of ~ 2 µPa/√Hz, which is comparable to noise floors of the current state-of-the-art optical microphones.

Based on this transducer, a PA sensor prototype, assembled with an “off-the-shelf” PA cell, has been proposed. The PA cell consists in a cylindrical multipass cavity with reflective walls to enhance the PA signal generation. A complete characterization of the FP based PA sensor was carried out with a QCL on a weak NO absorption line at 1906.15 cm⁻¹. The measurements showed a good stability and a LOD of 410 ppb/√Hz, with a NNEA of 5.1 × 10⁻⁷ W·cm⁻¹·Hz⁻¹/². Finally, we observed a good agreement between the experimental and calculated responses of the PA sensor.

In the near future, further improvements are expected, based on new prototypes with optimized PA chamber geometries and functioning at lower mechanical resonance frequency.
**Author Contributions:** Conceptualization, T.L., A.G. and S.B.; methodology, T.L. and A.G.; software, T.L.; validation, T.L.; formal analysis, T.L.; investigation, T.L.; data curation, T.L.; writing—original draft preparation, T.L.; writing—review and editing, T.L., A.G. and S.B.; visualization, T.L.; supervision, A.G. and S.B. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Acknowledgments:** The authors would like to thank Thierry Verdot for the fruitful discussions and advices concerning the development of the lumped model, and Brice Ensenat for the machining of the mechanical parts of the sensor.

**Conflicts of Interest:** The authors declare no conflict of interest.

## Appendix A. Hinged Cantilever Effective Mass and Stiffness

The effective mass, stiffness and area of the first bending mode of the hinged cantilever are given by:

$$ egin{cases} k_{\text{eff},m} = E_m \left( \int_0^{l_1} I_1 \left[ \frac{d^2 g_1(x)}{dx^2} \right]^2 dx + \int_{l_1}^{l_1+l_2} I_2 \left[ \frac{d^2 g_2(x)}{dx^2} \right]^2 dx \right) \\ m_{\text{eff},m} = t_m \rho_m \left( \int_0^{l_1} w_1 [g_1(x)]^2 dx + \int_{l_1}^{l_1+l_2} w_2 [g_2(x)]^2 dx \right), \\ S_{\text{eff},m} = t_m \left( w_1 \int_0^{l_1} g_1(x) dx + w_2 \int_{l_1}^{l_1+l_2} g_2(x) dx \right) \end{cases} \quad (A1) $$

where $g_1$ (resp. $g_2$) is the distribution function corresponding to the cantilever segment (resp. the hinge segment) and $w_i, l_i, I_i = \frac{w_i t_i^3}{12}$ are the width, length and moment of area of the considered segment:

$$ g_i(x) = a_i + b_i x + c_i x^3 \quad \text{with} \quad \left\{ \begin{array}{l} a_1 = 1 \\ b_1 = -\frac{3}{2} \frac{w_2 l_1^2 + w_1 l_2 (2l_1 + l_2)}{c_1 = \frac{1}{2} \frac{w_2}{A}} \end{array} \right. \quad \text{and} \quad \left\{ \begin{array}{l} a_2 = \frac{w_1 (l_1 + l_2)^3}{A} \\ b_2 = -\frac{3}{2} \frac{w_1 (l_1 + l_2)^2}{c_2 = \frac{1}{2} \frac{w_1}{A}} \\ c_2 = \frac{1}{2} \frac{w_1}{A} \end{array} \right. \quad (A2) $$

where $A = w_2 l_1^3 + w_1 l_2 (3l_1^2 + 3l_1 l_2 + l_2^2)$. The resonance frequency of the first bending mode is given by:

$$ \omega_{m,\text{hc}} = \sqrt{\frac{k_{\text{eff},m}}{m_{\text{eff},m}}} = \sqrt{\frac{w_2[l_2(3l_1^2 + 3l_1l_2 + l_2^2)w_1 + l_1^3w_2]}{B}} t \sqrt{\frac{E}{\rho_m}}, \quad (A3) $$

with $B = 35l_1l_2^2(12l_1^4 + 30l_1^3l_2 + 33l_1^2l_2^2 + 18l_1l_2^3 + 4l_2^4)w_1^2 + l_2(231l_1^6 + 273l_1^5l_2 + 105l_1^4l_2^2 + 63l_1^2l_2^4 + 91l_1l_2^5 + 33l_2^6)w_1w_2 + 33l_1^7w_2^2$.

## Appendix B. Lumped Model for the Heat Source

The $F$ function is a solution of the Helmholtz equation:

$$ \Delta F + k^2 F = 0 \quad \text{with} \quad k^2 = (-i) \frac{\rho_0 \omega}{\kappa/c_p} \quad \text{and} \quad r \in \partial V : F = 1, \quad (A4) $$

where $\kappa$ is the thermal conductivity of air. The sidewall thermal losses are given by the $R_r$ impedance calculated with $F(r, \omega)$ considering an infinitely long cylinder and the upper and lower thermal losses are given by the $R_z$ impedance calculated with $F(z, \omega)$ considering an infinitely flat cylinder:

$$ \left\{ \begin{array}{l} F(r, \omega) = \frac{J_0(kr)}{J_0(kR_{\text{cell}})} \quad \text{with } r \in [0, R_{\text{cell}}] \\ F(z, \omega) = \frac{\cos(kz)}{\cos\left(\frac{kh_{\text{cell}}}{2}\right)} \quad \text{with } z \in [-h_{\text{cell}}/2, h_{\text{cell}}/2] \end{array} \right. \quad \Rightarrow \quad \left\{ \begin{array}{l} R_r = \frac{\rho_0 c_p P_0}{8 \kappa \pi} \frac{1}{h} \\ R_z = \frac{\rho_0 c_p P_0}{8 \kappa \pi} \frac{2h}{3R^2} \end{array} \right. , \quad (A5) $$
where $J_0$ denotes the Bessel function of first kind. The total thermal losses $R_{th}$ are given by the relation:

$$ \frac{1}{R_{th}} = \frac{1}{R_r} + \frac{1}{R_z}, \qquad (A6) $$

References

1. Zhang, P.; Li, L.; Lin, L.; Shi, J.; Wang, L.V. In vivo superresolution photoacoustic computed tomography by localization of single dyed droplets. *Light Sci. Appl.* **2019**, *8*, 36. [CrossRef]

2. Wissmeyer, G.; Pleitez, M.A.; Rosenthal, A.; Ntziachristos, V. Looking at sound: Optoacoustics with all-optical ultrasound detection. *Light Sci. Appl.* **2018**, *7*, 1–16. [CrossRef]

3. Zheng, H.; Dong, L.; Patimisco, P.; Wu, H.; Sampaolo, A.; Yin, X.; Li, S.; Ma, W.; Zhang, L.; Yin, W.; et al. Double antinode excited quartz-enhanced photoacoustic spectrophone. *Appl. Phys. Lett.* **2017**, *110*, 021110. [CrossRef]

4. Chen, K.; Yu, Q.; Gong, Z.; Guo, M.; Qu, C. Ultra-high sensitive fiber-optic Fabry-Perot cantilever enhanced resonant photoacoustic spectroscopy. *Sens. Actuators B Chem.* **2018**, *268*, 205–209. [CrossRef]

5. Rouxel, J.; Coutard, J.-G.; Gidon, S.; Lartigue, O.; Nicoletti, S.; Parvitte, B.; Vallon, R.; Zéninari, V.; Glière, A. Miniaturized differential Helmholtz resonators for photoacoustic trace gas detection. *Sens. Actuators B Chem.* **2016**, *236*, 1104–1110. [CrossRef]

6. Kosterev, A.A.; Tittel, F.K.; Serebryakov, D.V.; Malinovsky, A.L.; Morozov, I.V. Applications of quartz tuning forks in spectroscopic gas sensing. *Rev. Sci. Instrum.* **2005**, *76*, 043105. [CrossRef]

7. Patimisco, P.; Borri, S.; Sampaolo, A.; Beere, H.E.; Ritchie, D.A.; Vitiello, M.S.; Scamarcio, G.; Spagnolo, V. A quartz enhanced photo-acoustic gas sensor based on a custom tuning fork and a terahertz quantum cascade laser. *Analyst* **2014**, *139*, 2079–2087. [CrossRef] [PubMed]

8. Patimisco, P.; Borri, S.; Galli, I.; Mazzotti, D.; Giusfredi, G.; Akikusa, N.; Yamanishi, M.; Scamarcio, G.; De Natale, P.; Spagnolo, V. High finesse optical cavity coupled with a quartz-enhanced photoacoustic spectroscopic sensor. *Analyst* **2015**, *140*, 736–743. [CrossRef]

9. Wang, Z.; Wang, Q.; Zhang, W.; Wei, H.; Li, Y.; Ren, W. Ultrasensitive photoacoustic detection in a high-finesse cavity with Pound–Drever–Hall locking. *Opt. Lett.* **2019**, *44*, 1924. [CrossRef]

10. Dostál, M.; Suchánek, J.; Válek, V.; Blatoňová, Z.; Nevrlý, V.; Bitala, P.; Kubát, P.; Zelinger, Z. Cantilever-Enhanced Photoacoustic Detection and Infrared Spectroscopy of Trace Species Produced by Biomass Burning. *Energy Fuels* **2018**, *32*, 10163–10168. [CrossRef]

11. Gallego, D.; Lamela, H. High-sensitivity ultrasound interferometric single-mode polymer optical fiber sensors for biomedical applications. *Opt. Lett.* **2009**, *34*, 1807. [CrossRef] [PubMed]

12. Kuusela, T.; Kauppinen, J. Photoacoustic Gas Analysis Using Interferometric Cantilever Microphone. *Appl. Spectrosc. Rev.* **2007**, *42*, 443–474. [CrossRef]

13. Murray, M.J.; Davis, A.; Redding, B. Fiber-Wrapped Mandrel Microphone for Low-Noise Acoustic Measurements. *J. Lightwave Technol.* **2018**, *36*, 3205–3210. [CrossRef]

14. Yu, Q.; Zhou, X. Pressure sensor based on the fiber-optic extrinsic Fabry-Perot interferometer. *Photonic Sens.* **2011**, *1*, 72–83. [CrossRef]

15. Wang, Q.; Wang, J.; Li, L.; Yu, Q. An all-optical photoacoustic spectrometer for trace gas detection. *Sens. Actuators B Chem.* **2011**, *153*, 214–218. [CrossRef]

16. Chen, K.; Guo, M.; Liu, S.; Zhang, B.; Deng, H.; Zheng, Y.; Chen, Y.; Luo, C.; Tao, L.; Lou, M.; et al. Fiber-optic photoacoustic sensor for remote monitoring of gas micro-leakage. *Opt. Express* **2019**, *27*, 4648. [CrossRef]

17. Liu, B.; Lin, J.; Liu, H.; Jin, A.; Jin, P. Extrinsic Fabry-Perot fiber acoustic pressure sensor based on large-area silver diaphragm. *Microelectron. Eng.* **2016**, *166*, 50–54. [CrossRef]

18. Tan, Y.; Zhang, C.; Jin, W.; Yang, F.; Ho, H.L.; Ma, J. Optical Fiber Photoacoustic Gas Sensor With Graphene Nano-Mechanical Resonator as the Acoustic Detector. *IEEE J. Sel. Top. Quantum Electron* **2017**, *23*, 199–209. [CrossRef]

19. Gong, Z.; Zhou, X.; Chen, K.; Zhou, X.; Yu, Q. Parylene-C diaphragm-based fiber-optic Fabry-Perot acoustic sensor for trace gas detection. In Proceedings of the 2017 International Conference on Optical Instruments and Technology: Advanced Optical Sensors and Applications, Beijing, China, 28–30 October 2017; Dong, L., Zhang, X., Xiao, H., Arregui, F.J., Eds.; SPIE: Bellingham, WA, USA, 2018; p. 5.
20. Chen, K.; Gong, Z.; Guo, M.; Yu, S.; Qu, C.; Zhou, X.; Yu, Q. Fiber-optic Fabry-Perot interferometer based high sensitive cantilever microphone. *Sens. Actuators A Phys.* **2018**, *279*, 107–112. [CrossRef]

21. Zhang, B.; Chen, K.; Chen, Y.; Yang, B.; Guo, M.; Deng, H.; Ma, F.; Zhu, F.; Gong, Z.; Peng, W.; et al. High-sensitivity photoacoustic gas detector by employing multi-pass cell and fiber-optic microphone. *Opt. Express* **2020**, *28*, 6618. [CrossRef]

22. Lauwers, T.; Glière, A.; Basrour, S. Resonant optical transduction for photoacoustic detection. In Proceedings of the Photonic Instrumentation Engineering VII, San Francisco, CA, USA, 4–6 February 2020; Soskind, Y., Busse, L.E., Eds.; SPIE: Bellingham, WA, USA, 2020; p. 25.

23. Glière, A.; Rouxel, J.; Brun, M.; Parvitte, B.; Zéninari, V.; Nicoletti, S. Challenges in the Design and Fabrication of a Lab-on-a-Chip Photoacoustic Gas Sensor. *Sensors* **2014**, *14*, 957–974. [CrossRef] [PubMed]

24. Lobontiu, N. *Dynamics of Microelectromechanical Systems*; Microsystems: Santa Clara, CA, USA; Springer: New York, NY, USA, 2007; ISBN 978-0-387-36800-9.

25. Benade, A.H. On the propagation of sound waves in a cylindrical conduit. *J. Acoust. Soc. Am.* **1968**, *44*, 616–623. [CrossRef]

26. Kampinga, W.R. Viscothermal Acoustics Using Finite Elements—Analysis Tools for Engineers. Ph.D. Thesis, University of Twente, Enschede, The Netherlans, 2010. [CrossRef]

27. Schilt, S.; Thévenaz, L. Wavelength modulation photoacoustic spectroscopy: Theoretical description and experimental results. *Infrared Phys. Technol.* **2006**, *48*, 154–162. [CrossRef]

28. Mao, X.; Yuan, S.; Zheng, P.; Wang, X. Stabilized Fiber-Optic Fabry-Perot Acoustic Sensor Based on Improved Wavelength Tuning Technique. *J. Lightwave Technol.* **2017**, *35*, 2311–2314. [CrossRef]

<figure><img src="image_20.png" /></figure>

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).