1. Introduction
Photoacoustic spectroscopy (PAS) is based on the absorption of modulated laser radiation and the generation of a sound wave. In gaseous samples, light energy is absorbed by transforming it into molecular vibrations, rotations or electronic excitations, depending on the wavelength of the light (i.e., its energy). The modulated excitation causes local pressure changes that generate a sound wave with the same frequency as the light modulation.
A variety of acoustic detectors are available to record photoacoustic signals. High-precision measuring microphones, MEMS microphones or ultrasonic detectors, and piezoelectric sensors are used in photoacoustic sensors. Piezoelectric sensors can be made of different materials like polyvinylidene fluoride, lead zirconate titanate or polypropylene films and quartz crystals [
1]. Quartz tuning forks (QTFs) are micromechanical tuning forks made of silicon dioxide. They vibrate at their maximum when exited with their resonance frequency, generating a voltage corresponding to the charge displacement caused by the vibration. Standard tuning forks for quartz watches exhibit a resonance frequency of 32,768 Hz. QTFs have the following advantages [
1,
2,
3,
4]:
They separate noise from the rest of the signal and, therefore, can reach a higher specific sensitivity than microphones.
They are temperature-stable over a wide frequency range and insensitive to magnetic fields.
They are inexpensive and miniaturised, making them suitable for mass production.
They are more robust than other piezoelectric sensors.
They typically have high resonance quality factors.
Quartz-Enhanced PAS (QEPAS) combines the advantages of QTFs with those of photoacoustic spectroscopy. So far, this technique has been exclusively used for gaseous samples. A variety of configurations have been thoroughly tested, with the most common being the on-beam configuration. This involves the light passing through the gap between the prongs of the QTF. To amplify the signal, the QTF is often positioned between two identical cylindrical resonance tubes.
In the off-beam configuration, the sound pressure wave is detected at a certain distance from the laser beam. The QTF is here typically positioned outside a resonance tube at a narrow opening. This slit or hole has a width on the order of the spacing of the prongs. There are other off-beam approaches with two tubes arranged in a T-shape [
5]. Implementations using two parallel tubes with facing slits and a single QTF [
6] aim for increased amplification. Alternatively, two T-shaped resonators with opposite openings were tested with the QTF in between. Since each prong was facing an opening, the sound hit the QTF from both sides [
7]. Other tested resonator geometries include H-shaped cells [
8], Helmholtz resonators [
9], elliptical tubes [
10] or semicylindrical resonators [
11].
Apart from theoretical considerations of a tuning fork immersed in a liquid medium [
12] and using a QTF to optically analyse solid samples with so-called light-induced thermoelastic spectroscopy (LITES) [
13], there are currently no studies of QEPAS on solid or liquid samples. To our knowledge, we performed the first QEPAS measurements on a solid sample and designed a special cell for this purpose. The requirements for such a cell differ from those for gaseous samples.
The following section briefly describes the photoacoustic effect of solids. After that, we present the new cell design, which is supported by a three-dimensional Finite Element simulation.
Section 4 shows the experimental methods and results, including the characterisation of the cell resonance using an ultrasound (US) microphone and the QEPAS measurements. A discussion and an outlook are presented in the concluding
Section 5.
2. Photoacoustic Effect of Solids
In solid samples, the absorption of periodically modulated optical energy leads to lattice vibrations or phonons. The resulting heat modulation has the same frequency as the light modulation and propagates to the sample boundary. A thin layer of the coupling medium in front of the sample acts like a piston on the rest of the volume, causing periodic pressure fluctuations in the gas. Thus, a photoacoustic (PA) signal is generated by a periodic temperature flow from the sample surface to an adjacent gaseous medium, which should exhibit low absorption in the spectral excitation range [
14,
15]. In contrast to traditional absorption spectroscopy, PAS is not affected by scattered light, which can cause inaccurate results when investigating solid samples. Additionally, PAS can also be used to analyse completely non-transparent samples [
16].
In the case of an optically opaque and thermally thick sample, the acoustic sound pressure fluctuation
as a result of the movement of the gas piston depends on the time
and the circular modulation frequency
[
14].
The factor is as a function of the thermal diffusivity of the gas in front of the sample (air) and the following properties of the sample:
optical absorption coefficient ,
thermal diffusivity of material ,
thermal diffusion length in the sample ( = ),
thermal conductivity of the material
[
17].
The factor
depends on the incident light power
, the length of the gas piston
being equal to the length of the cylindric resonator adjacent to the sample, the ratio of the specific heats in the adiabatic gas law
and the overall temperature
(ambient temperature and the DC-component at the solid surface due to temperature rise) and the ambient pressure
[
17].
This shows the photoacoustic signal is dependent on the absorption and thermal properties of the sample as well as on ambient parameters such as air temperature and thermal air diffusivity. The signal increases with increasing absorption in the sample and increasing light intensity.
The depth of penetration depends on the absorption coefficient of the sample and the incident light energy, which decreases proportional to
. The thermal energy generated by absorption contributes only up to the simple thermal diffusion length
to the photoacoustic signal [
14].
3. QEPAS Cell for Solid Samples
A semi-open cylindrical cell for the QEPAS analysis of solid samples has been developed. The manufacturing process for the cell is not only simpler compared to cuboid or T-/H-shaped cells [
18,
19,
20,
21]; the PA signal generation is also more efficient because the aforementioned geometries exhibit pipe diameter changes, which results in additional damping due to impedance discontinuities, leading to significant energy losses [
22].
3.1. Cell Design
The new cylindrical cell is tightly sealed at one end by the solid sample and thus represents a half-open cylinder. The open end enables free access to the sample for the laser beam. Closed ends result in a 180-degree phase shift of the reflected sound wave, while open ends do not cause a phase shift. The resonance frequencies
of the half-open cylinders are described by the equation [
23].
where
is the sound velocity,
is the geometric length of the cylinder and
is the end correction which mainly depends on the radius of the cylinder
and can be calculated from [
24].
For a maximum photoacoustic signal, modulating the radiation with a frequency that corresponds to the most intense acoustic resonance is recommended. Here, we use the first overtone, which is also called the second harmonic (
), as the fundamental leads to a maximum directly at the sample and higher modes have an increasingly lower intensity [
25].
Figure 1a shows the corresponding sound pressure distribution.
In the experiments described in
Section 4, a commercial QTF will serve as a sound detector for the PA signal. Its nominal resonance frequency is
kHz. Since the sample is excited by a laser beam, the diameter of the cylinder must be larger than the beam width. The laser beam should not hit the wall of the tube. Otherwise, a wall signal will be generated that will interfere with the desired sample signal. In order to fulfil this requirement, a tube radius of
was selected, which is significantly larger than the laser spot size of ca. 0.1 mm.
Using Equations (1) and (2) and taking into account the second harmonic , which corresponds to the QTF resonance frequency , a geometric length of is obtained. The value was used for the speed of sound.
3.2. Finite Element Simulation
We calculated the length of a cylindrical cell that can serve as a half-open resonator for solid samples. Previously, a 2D model was developed employing COMSOL Multiphysics
® to determine the length of the cylinder and the location of the opening [
20]. Subsequently, we performed a 3D finite element (FE) simulation. The homogeneous Helmholtz equation was solved according to the acoustic sound pressure field. One condition was that the sound pressure at the walls must be maximum, and therefore, its normal derivative must be zero. It was assumed that the sample represented a soundproof seal of the tube. The particular difficulty of the open end was solved with perfectly matched layers (PML). Using an iterative approach, we determined the length of the cylinder at which the first overtone matched the resonance frequency of the QTF.
The additional dimension required more computational effort but is considered to be more accurate.
Figure 1b shows the 3D sound pressure distribution of the second harmonic longitudinal resonance of the half-open cylinder.
The hemisphere attached to the top of the resonator represents the external space that integrates the PML. The sound pressure of the acoustic mode spreads into the external space, which makes the end correction necessary for the analytical calculation. The sound pressure distribution shows that the second harmonic results in two maximum amplitudes. One is located 1.37 mm from the open end and is determined as the location of the signal detection.
The 3D simulation confirms the previous 2D results [
23]. There was no significant difference in the sound pressure distribution of the first overtone. However, the accuracy has improved due to fewer approximations and the consideration of azimuthal and radial modes at higher frequencies.
Table 1 shows the results of the analyses regarding the dimensions of the cylinder geometry. Compared to the analytical calculation, the required cell is, according to the FE simulation,
(0.037 mm) larger.
4. Experimental Investigations
4.1. Experimental Setup
The cylindrical PA cell was made of stainless steel with a wall thickness of 0.100 mm. The accuracy of the cell length equals 0.005 mm. The optimal location for the sound detector is where the maximum sound pressure occurs. In order to achieve this as best as possible with minimum distortion of the sound field, a hole of 0.250 mm diameter is laser-drilled into the cell wall, 1.370 mm from the open end. The sound detector is then placed outside the cell, closest to the hole.
For the verification of the cell resonance presented in
Section 4.2, an ultrasound (US) micro-electromechanical system (MEMS) microphone is utilised [
26]. For the QEPAS measurements presented in
Section 4.3, the microphone is replaced by a QTF. Custom-made electronics based on high-precision differential amplifiers boost the respective signals.
In both investigations, a graphite disk serves as a test sample. This carbon material represents a spectrally broadband absorber. At the wavelength of 3380 nm, the material exhibits a relatively high absorption value for a solid, which is advantageous for the initial functionality test of our sensor. The absorption of the material for the given frequency was measured using FTIR-ATR (ThermoFisher Scientific Nicolet iS5, Waltham, MA, USA), resulting in an absorbance of 0.265, equivalent to 54.33% transmission. Its surface was repeatedly sanded with fine grit to achieve a soundproof seal with the cell.
A mounting frame for the resonance cell and the sample was designed and 3D-printed from clear resin material. The version for the QTF is schematically depicted in
Figure 2. The laser-drilled hole is marked with an arrow.
Figure 3 shows the experimental setup. A continuous wave (cw) distribute feedback (DFB) interband cascade laser (ICL) is used to excite the sample. Its emission is focused onto the sample surface.
The laser is regulated to a constant temperature of 20 °C by a thermoelectric cooler (TEC). A laser diode driver, in combination with a signal generator, modulates the laser between threshold and maximum current/maximum optical output power. A customised Matlab script automatically tunes the modulation frequency. A lock-in amplifier detects the PA signal. The complete setup is controlled by a PC that records the PA signal as a function of frequency.
Table 2 provides specifications of the experimental setup.
4.2. Cell Resonance Verification
Prior to the QEPAS measurements, the resonance frequency of the new PA cell, which was determined by FE simulation (
Section 2), should be verified. For this, we measured the PA signal generated by the sample with the US MEMS microphone as a function of the modulation frequency of the laser. The cell was filled with air at 22 °C and 48% relative humidity and ambient air pressure of 101.325 kPa. None of the air components show significant absorption in the laser’s spectral range. Subsequently, we normalised the measurement according to the microphone’s transfer function [
26].
Figure 4 displays the results. The normalised PA signal is represented by crosses. A Lorentz profile was fitted to the result, and the QTF resonance frequency is indicated by the dashed vertical line. Due to the challenges of the open resonator, a relatively low-quality factor of q = 4 is still achieved (3 dB-bandwidth of 4 kHz).
The Lorentz profile’s maximum corresponds almost precisely with the QTF resonance. The slight shift of only 0.1% to higher frequencies is likely due to the increased speed of sound at higher temperatures. The measurement, therefore, strongly supports the simulation results, demonstrating the accuracy of our methods.
4.3. QEPAS on Solid Samples
A standard QTF with a resonance frequency of 32,768 Hz was used for the QEPAS investigation. Its electronic original q-factor listed by the manufacturer (cf.
Table 2) is 60,000 in a vacuum. This is lowered to around 10,000 due to the damping of air [
13]. The precise positioning of the QTF is particularly important in this special off-beam configuration. The sound pressure wave should hit the upper third of the prongs [
27]. Furthermore, both prongs should have the exact same distance to the sound source. Otherwise, one prong is more deflected than the other, which leads to asymmetries in the resonance function [
28]. Additionally, the QTF should be positioned close to the hole in the cylinder to minimise the loss of sound pressure energy into the external space. However, if the QTF is placed too close, the signal could be reduced due to squeeze damping [
29]. In order to achieve the best possible positioning, the QTF was mounted to a three-axis displacement unit and then placed approximately 0.5 mm from the laser-drilled hole parallel to the resonator wall. The centre of the hole was 2.6 mm below the QTF tip. To prevent direct excitation of the QTF by the laser (light-induced thermoelastic spectroscopy—LITES), it was shielded with an iris (Thorlabs ID15/M).
Figure 5 displays the QEPAS signal of the solid sample as a function of the modulation frequency of the laser. The resonance curve is consistently smooth and exhibits only low noise.
Figure 5 also shows a Lorentz fit, which is in good agreement with the measured data. It results in a resonance frequency of 32,754.2 Hz and a resonance width of 4.72 Hz (FWHM). This corresponds to a quality factor of
6.939. The QTF’s nominal resonance frequency is 32,768 Hz, with the manufacturer stating deviations of up to 5 ppm.
To investigate the dependency of the sensor signal on the sample absorbance, we measured the QEPAS signal of the following three different solid synthetics: phenol formaldehyde resin, polyamide and polypropylene. They cover a relatively wide range of absorbance [
30].
Figure 6 shows the results together with a linear fit. The signal increases continuously with increasing absorbance.
5. Conclusions
We presented the first QEPAS measurement on a solid sample. A cw DFB ICL with a wavelength of around 3380 nm served as a light source. A disk made of graphite represented the sample.
A special cell has been developed based on the second harmonic of the fundamental longitudinal resonance of a half-open cylinder. The cell exhibits an acoustical resonance at the resonance frequency of the QTF in order to additionally amplify the PA signal and, thus, enhance the detection sensitivity. The required cell length was determined analytically and verified numerically using an FE simulation. A measurement with an ultrasound MEMS microphone confirmed the resonance frequency of the cell.
A tiny hole in the cell wall at the location of the pressure anti-node allowed the sound detector to access the sound pressure maximum outside the cell with only minimal distortion of the sound field. By tuning the modulation frequency of the laser, we were able to measure the ultrasound resonance function of the QTF around 32,754 Hz. The resonance quality was approximately 7000. A subsequent investigation of the QEPAS signals of three different solid synthetics shows a linear relationship between absorbance and signal amplitude.
The results open up completely new areas of application for the powerful QEPAS technology, as the new sensor can be adapted to a variety of scientific or industrial applications that require the investigation of solid samples.
However, there are still some challenges that simultaneously represent opportunities for future research. Precise sensor positioning under different environmental conditions is, for instance, crucial to achieving optimal performance and application-specific calibration. In order to adapt the QTF to a given cell, it could also be advantageous to use a custom-made tuning fork [
31,
32,
33].
Author Contributions
Conceptualization, M.W., J.F. and B.B.; methodology, J.F., M.W., B.B. and M.-S.B.; software, J.F. and B.B.; investigation, J.F. and M.-S.B.; writing—original draft preparation, J.F.; writing—review and editing, J.F., M.W. and M.-S.B.; visualisation, J.F.; supervision, M.W. All authors have read and agreed to the published version of the manuscript.
Funding
This research received no external funding.
Institutional Review Board Statement
Not applicable.
Informed Consent Statement
Not applicable.
Data Availability Statement
The data presented in this study are available on request from the corresponding author (accurately indicate status).
Conflicts of Interest
The authors declare no conflict of interest.
References
- Lv, H.; Luo, H.; Zheng, H.; Zhu, W.; Fang, J.; Yu, J.; Tittel, F.; Chen, Z. Application of piezoelectric transducers in photoacoustic spectroscopy for trace gas analysis. Microw. Opt. Technol. Lett. 2021, 63, 2040–2051. [Google Scholar] [CrossRef]
- Zifarelli, A.; Menduni, G.; Giglio, M.; Sampaolo, A.; Patimisco, P.; Wu, H.; Dong, L.; Spagnolo, V. Latest Advances in Quartz Enhanced Photoacoustics Spectroscopy for Environmental and Industrial Applications. In Proceedings of the 23rd International Conference on Transparent Optical Networks (ICTON), Bucharest, Romania, 2–6 July 2023; IEEE: Piscataway, NJ, USA, 2023; pp. 1–3. [Google Scholar]
- Ayache, D.; Rousseau, R.; Kniazeva, E.; Charensol, J.; Seoudi, T.; Bahriz, M.; Gouzi, F.; Spagnolo, V.; Vicet, A. Commercial and Custom Quartz Tuning Forks for Quartz Enhanced Photoacoustic Spectroscopy: Stability under Humidity Variation. Sensors 2023, 23, 3135. [Google Scholar] [CrossRef] [PubMed]
- Sampaolo, A.; Patimisco, P.; Giglio, M.; Zifarelli, A.; Wu, H.; Dong, L.; Spagnolo, V. Quartz-enhanced photoacoustic spectroscopy for multi-gas detection: A review. Anal. Chim. Acta 2022, 1202, 338894. [Google Scholar] [CrossRef] [PubMed]
- Rousseau, R.; Loghmari, Z.; Bahriz, M.; Chamassi, K.; Teissier, R.; Baranov, A.N.; Vicet, A. Off-beam QEPAS sensor using an 11-μm DFB-QCL with an optimized acoustic resonator. Opt. Express 2019, 27, 7435–7446. [Google Scholar] [CrossRef] [PubMed]
- Lin, C.; Xueyang, Y.; Yongchao, H. An all-optical off-beam quartz-enhanced photoacoustic spectroscopy employing double-pass acoustic microresonators. Opt. Commun. 2022, 503, 127447. [Google Scholar] [CrossRef]
- Rousseau, R.; Ayache, D.; Maurin, N.; Trzpil, W.; Bahriz, M.; Vicet, A. Monolithic double resonator for quartz enhanced photoacoustic spectroscopy. Appl. Sci. 2021, 11, 2094. [Google Scholar] [CrossRef]
- Ma, Y.; Hong, Y.; Qiao, S.; Lang, Z.; Liu, X. H-shaped acoustic micro-resonator-based quartz-enhanced photoacoustic spectroscopy. Opt. Lett. 2022, 47, 601–604. [Google Scholar] [CrossRef] [PubMed]
- Luo, H.; Wang, C.; Lin, H.; Wu, Q.; Yang, Z.; Zhu, W.; Zhong, Y.; Kan, R.; Yu, J.; Zheng, H. Helmholtz-resonator quartz-enhanced photoacoustic spectroscopy. Opt. Lett. 2023, 48, 1678–1681. [Google Scholar] [CrossRef] [PubMed]
- Shang, Z.; Wu, H.; Li, S.; Tittel, F.K.; Dong, L. Elliptical-tube off-beam quartz-enhanced photoacoustic spectroscopy. Appl. Phys. Lett. 2022, 120, 171101. [Google Scholar] [CrossRef]
- De Carlo, M.; Menduni, G.; Sampaolo, A.; De Leonardis, F.; Spagnolo, V.; Passaro, V.M. Modeling and design of a semi-integrated QEPAS sensor. J. Light. Technol. 2021, 39, 646–653. [Google Scholar] [CrossRef]
- Campanale, A.; Putignano, C.; De Carolis, S.; Patimisco, P.; Giglio, M.; Soria, L. A theoretical-experimental framework for the analysis of the dynamic response of a QEPAS tuning fork device immersed in a fluid medium. Mech. Syst. Signal Process. 2021, 149, 107298. [Google Scholar] [CrossRef]
- Van Neste, C.W.; Senesac, L.R.; Thundat, T. Standoff photoacoustic spectroscopy. Appl. Phys. Lett. 2008, 92, 234102. [Google Scholar] [CrossRef]
- Thakur, S.N.; Virendra, N.R.; Singh, J.P. (Eds.) Photoacoustic and Photothermal Spectroscopy: Principles and Applications; Elsevier: Amsterdam, The Netherlands, 2022. [Google Scholar]
- Sarode, A.P.; Mahajan, O.H. Theoretical aspects of Photoacoustic Signal Generation with Solid Crystals. Bulletin of the Karaganda University. Phys. Ser. 2021, 4, 104. [Google Scholar]
- Raghu, O.; Philip, J. A dual channel photoacoustic cell for imaging experiments on solid samples. J. Instrum. Soc. India 2003, 33, 155–158. [Google Scholar]
- Rosencwaig, A.; Gersho, A. Theory of the photoacoustic effect with solids. J. Appl. Phys. 1976, 47, 64–69. [Google Scholar] [CrossRef]
- Rabasović, M.D.; Nikolić, M.G.; Dramićanin, M.D.; Franko, M.; Markushev, D.D. Low-cost, portable photoacoustic setup for solid samples. Meas. Sci. Technol. 2009, 20, 095902. [Google Scholar] [CrossRef]
- El-Busaidy, S.; Baumann, B.; Wolff, M.; Duggen, L.; Bruhns, H. Experimental and numerical investigation of a photoacoustic resonator for solid samples: Towards a non-invasive glucose sensor. Sensors 2019, 19, 2889. [Google Scholar] [CrossRef] [PubMed]
- El-Busaidy, S.; Baumann, B.; Wolff, M.; Duggen, L. Modelling of open photoacoustic resonators. Photoacoustics 2020, 18, 100161. [Google Scholar] [CrossRef] [PubMed]
- El-Busaidy, S.; Baumann, B.; Wolff, M.; Duggen, L. Shape optimization of an open photoacoustic resonator. Appl. Sci. 2021, 11, 2571. [Google Scholar] [CrossRef]
- Peat, K.S. The acoustical impedance at discontinuities of ducts in the presence of a mean flow. J. Sound. Vib. 1988, 127, 123–132. [Google Scholar] [CrossRef]
- Falkhofen, J.; Bernd, B.; Wolff, M. New Photoacoustic Cell Design for Solid Samples. Eng. Proc. 2023, 58, 86. [Google Scholar]
- Levine, H.; Schwinger, J. On the radiation of sound from an unflanged circular pipe. Phys. Rev. 1948, 73, 383. [Google Scholar] [CrossRef]
- Wu, D.; Gao, S.; Li, J.; Yao, L.; Yu, X.; Zhang, Z.; Wang, X. Amplitude and phase relation of harmonics in nonlinear focused ultrasound. AIP Adv. 2022, 12, 065317. [Google Scholar] [CrossRef]
- Falkhofen, J.; Wolff, M. Near-Ultrasonic Transfer Function and SNR of Differential MEMS Microphones Suitable for Photoacoustics. Sensors 2023, 23, 2774. [Google Scholar] [CrossRef] [PubMed]
- Patimisco, P.; Sampaolo, A.; Zheng, H.; Dong, L.; Tittel, F.K.; Spagnolo, V. Quartz–enhanced photoacoustic spectrophones exploiting custom tuning forks: A review. Adv. Phys. X 2017, 2, 169–187. [Google Scholar] [CrossRef]
- Lee, M.; Kim, B.; An, S.; Jhe, W. Dynamic responses of electrically driven quartz tuning fork and qPlus sensor: A comprehensive electromechanical model for quartz tuning fork. Sensors 2019, 19, 2686. [Google Scholar] [CrossRef] [PubMed]
- Olivieri, M.; Zifarelli, A.; Menduni, G.; Di Gioia, M.; Marzocca, C.; Passaro, V.M.; Patimisco, P. Influence of air pressure on the resonance properties of a t-shaped quartz tuning fork coupled with resonator tubes. Appl. Sci. 2021, 11, 7974. [Google Scholar] [CrossRef]
- Sadlowsky, B.; Hamburg University of Applied Sciences, Hamburg, Germany. Personal communication, 2024.
- 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. [Google Scholar] [CrossRef] [PubMed]
- Duquesnoy, M.; Aoust, G.; Melkonian, J.-M.; Lévy, R.; Raybaut, M.; Godard, A. Quartz enhanced photoacoustic spectroscopy based on a custom quartz tuning fork. Sensors 2019, 19, 1362. [Google Scholar] [CrossRef] [PubMed]
- Patimisco, P.; Sampaolo, A.; Dong, L.; Tittel, F.K.; Spagnolo, V. Recent advances in quartz enhanced photoacoustic sensing. Appl. Phys. Rev. 2018, 5, 011106. [Google Scholar] [CrossRef]
| Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 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 (https://creativecommons.org/licenses/by/4.0/).