Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Experimental demonstration of multi-parameter sensing based on polarized interference of polarization-maintaining few-mode fibers

Open Access Open Access

Abstract

Polarization-maintaining few-mode fibers (PM-FMFs) have found applications in sensors with attractive features that traditional single-mode fibers do not possess. We propose a measurement method for multi-parameter sensing based on PM-FMFs. This method is mainly based on the polarized interference effect of PM-FMF. We experimentally demonstrate the principle of the simultaneous sensing of temperature, strain and transverse force. The sensor has been verified by three groups of multi-parameter measurement experiments showing that the changes of multi-parameter derived from the wavelength shift of the polarized interference spectrum are consistent with the multi-parameter changes set in the experiments. The proposed PM-FMF based sensor has a temperature sensitivity of 0.3 nm/°C, a strain sensitivity of 0.01 nm/µε and a transverse force sensitivity of 0.0065 nm/(N/m).

© 2020 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Recently, few-mode fiber (FMF) has attracted significant attention due to its potential to solve the capacity crunch problem through mode-division multiplexing (MDM) [16]. It is important to mention that FMF also becomes a research focus in the field of fiber sensing [710]. Employing the spatial dimension, sensors based on FMF has the ability to measure a variety of physical disturbances (temperature, stain and transverse force, etc.) with high sensitivity and insulation. The existing FMF sensing systems based on Brillouin frequency shift (BFS) in each mode have achieved distributed temperature and strain discrimination [11,12]. It has solved the problem in single-mode fiber (SMF) sensors that the sensor performance can be jointly affected by parameters. However, the FMF sensing systems are quite complicated and the DSP algorithms are indispensable. Besides, the FMF sensors can not be applied for simultaneous measurement of multiple parameters. In this case, researchers turn their attention to special FMFs [1317]. The application of PM-FMF in sensing fields is an interesting area, which is only briefly investigated at best. It’s both the polarization-maintaining and few-mode property of PM-FMF that make it more suitable for multi-parameter sensing. PM-FMF sensor is simpler than FMF-based Brillouin sensor and it can achieve simultaneous multi-parameter sensing.

In this work, we first propose a multi-parameter sensor based on the polarized interference of the PM-FMF, which can achieve the high-precision sensing of temperature, strain and transverse force simultaneously. Because of the linearly independence, the temperature, strain and transverse force sensing coefficients of the LP01, LP11a and LP11b modes can be measured respectively. We experimentally demonstrate the effectiveness with a high temperature sensitivity of 0.3 nm/°C, strain sensitivity of 0.01 nm/µε and transverse force sensitivity of 0.0065 nm/(N/m), respectively.

2. Principle and experimental setup

The basic principle supporting the multi-parameter sensing measurement method is that the LP modes in PM-FMF have different beat lengths [18], due to the refractive index difference of two orthogonal polarization modes is different. The changes of external parameters cause the wavelength shift of polarized interference spectrum. The sensing coefficients of the single parameters in each mode can be measured in the single parameter sensing experiment, and the variation of the parameters in the multi-parameter sensing experiment can be calculated by the wavelength shift of the spectrum of each mode. The effect of the external parameter fields on the polarized interference spectrum can be expressed by the following equation:

$$A = k \cdot B$$
where $A = {[{\Delta {\lambda_{\textrm{01}}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \Delta {\lambda_{\textrm{11a}}}{\kern 1pt} {\kern 1pt} {\kern 1pt} \Delta {\lambda_{11b}}} ]^T}$ , $B = {[{\Delta T{\kern 1pt} {\kern 1pt} {\kern 1pt} \Delta S{\kern 1pt} {\kern 1pt} {\kern 1pt} \Delta F} ]^T}$ and $k = \left[ \begin{array}{l} {k_{T01}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {k_{S01}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {k_{F01}}\\ {k_{T11a{\kern 1pt} {\kern 1pt} {\kern 1pt} }}{k_{S11a}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {k_{F11a}}\\ {k_{T11b}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {k_{S11b}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {k_{F11b}} \end{array} \right]$. A is the wavelength shifts in the polarized interference spectrum of each mode. B is the changes of the external parameter fields. ${k_{T01}}$, ${k_{T11a}}$ and ${k_{T11b}}$ are the temperature sensing coefficients, ${k_{S01}}$, ${k_{S11a}}$ and ${k_{S11b}}$ are the strain sensing coefficients, ${k_{F01}}$, ${k_{F11a}}$ and ${k_{F11b}}$ are the transverse force sensing coefficients of the LP01, LP11a and LP11b mode respectively.

For PM-FMF, the beat length of each LP mode is different, and the effects of the external parameters are independent of each other, so the row vector or column vector of the matrix k acting as multi-parameter coefficients obtained by the polarized interference spectrum is linearly independent. In other words, the multi-parameter coefficients matrix k is invertible. By measuring the wavelength shift of the polarized interference spectrum of each mode, the variation of the multi-parameters can be calculated according to the following equation:

$$B = {k^{\textrm{ - 1}}} \cdot A$$

Figure 1 shows the experimental setup of multi-parameter sensing system. In this PM-FMF sensing system, an amplified spontaneous emission (ASE) broadband light source is coupled into one end of the PM-FMF which is called free-space fiber coupling. The phase plate converts the LP01 mode into a higher-order mode, and the linear polarizer 1 is inserted to generate linearly polarized light. After coupling into the PM-FMF, linear polarizer 2 is inserted to generate polarized interference effect between two orthogonal polarizing modes in the same LP mode. The same phase plate can be used to convert back the excited higher-order mode into the LP01 mode, then the free-space light is coupled into the single-mode fiber through a collimator, and the polarized interference spectrum is obtained by an optical spectrum analyzer (OSA). The heater in Fig. 1 plays the role of heating the PM-FMF. The fiber rotators are used for holding the PM-FMF, and at the same time, they can rotate the angle to change the force direction on the PM-FMF. Each fiber rotator is placed on a 3-axis flexure stage, we can tune the fine-threaded drives on the stage to stretch the PM-FMF.

 figure: Fig. 1.

Fig. 1. Experimental setup for multi-parameter sensing measurement. a: amplified spontaneous emission (ASE) light source; b: collimator; c: polarizer 1; d: phase plate; f: collimator; g: fiber rotator; h: heating panel; n: optical spectrum analyzer (OSA).

Download Full Size | PDF

The polarized interference intensity of each mode i can be written as

$${I_i} = A_i^2({R^2} + {S^2} + 2RS\cos {\phi _i})$$
where ${A_i}$ is the optical field amplitude, $R = \cos (\alpha )\cos (\gamma )$, $S = \sin (\alpha )\sin (\gamma )$, $\alpha$ is the angle between the principal axis of the polarizer 1 and the principal axis of the PM-FMF, $\gamma$ is the angle between the principal axis of the polarizer 2 and the principal axis of the PM-FMF, the phase difference ${\phi _i}$ accumulated by sensing of two polarization modes of mode i is given as
$${\phi _i} = \frac{{\textrm{2}\pi }}{\lambda }{B_{ai}}L + \frac{{\textrm{2}\pi }}{\lambda }{B_{bi}}{L_E}$$
where ${B_{ai}} = {n_{e - i}} - {n_{o - i}}$, ${n_{e - i}}$ and ${n_{o - i}}$ are the initial refractive index of mode i in the slow and fast axis of PM-FMF respectively, $L$ is the fiber length, $\lambda$ is the optical wavelength, ${B_{bi}}$ is the parameter-induced excess birefringence, ${L_E}$ is the length of the PM-FMF affected by the external parameter field.

The temperature T and phase ${\phi _i}$ satisfies the following relationship

$$\frac{{d{\phi _i}}}{{dT}} = \frac{{2\pi }}{\lambda }\left( {{L_H} \cdot \frac{{d{B_{bi}}}}{{dT}} + {B_{bi}} \cdot \frac{{d{L_H}}}{{dT}}} \right)$$
where ${L_H}$ is the length of heated fiber. The birefringence of polarization mode ${B_{bi}}$ is temperature dependent, and this phenomenon is called thermo-optic effect. Figure 2 shows the wavelength shift as the temperature increases and the reference transmission spectrum around the transmission minimum we selected. The effect of temperature on fiber length is negligible, that is $d{L_H}/dT = 0$. According to the equation in [19], when the temperature changes, the wavelength shift can be expressed as
$$\frac{{d\lambda }}{{dT}} = {L_H} \cdot {\frac{{d{B_{bi}}}}{{dT}}} \left/ {\left( {L \cdot \frac{{d{B_{ai}}}}{{d\lambda }} - \frac{{{B_{ai}}L + {B_{bi}}{L_H}}}{\lambda }} \right)}\right.$$

 figure: Fig. 2.

Fig. 2. The polarized interference spectrum of the PM-FMF and the wavelength shift as the temperature increases.

Download Full Size | PDF

The strain S and the phase ${\phi _i}$ satisfies the following relationship

$$\frac{{d{\phi _i}}}{{dS}} = \frac{{2\pi }}{\lambda }\left( {{L_S} \cdot \frac{{d{B_{bi}}}}{{dS}} + {B_{bi}} \cdot \frac{{d{L_S}}}{{dS}}} \right)$$
where $S = d{L_S}/{L_S}$ is the applied strain, $d{L_S}$ is the stretched length of PM-FMF and ${L_S}$ is the length of the PM-FMF between fiber holders. The birefringence of polarization mode ${B_{ei}}$ is strain dependent, and this phenomenon is called elasto-optic effect. The effect of strain on fiber length is negligible, that is $d{L_S}/dS = 0$, therefore when the strain changes, the wavelength shift can be expressed as
$$\frac{{d\lambda }}{{dS}} = {L_S} \cdot {\frac{{d{B_{bi}}}}{{dS}}} \left/ {\left( {L \cdot \frac{{d{B_{ai}}}}{{d\lambda }} - \frac{{{B_{ai}}L + {B_{bi}}{L_S}}}{\lambda }} \right)}\right.$$

The transverse force W and the phase ${\phi _i}$ satisfies the following relationship

$$\frac{{d{\phi _i}}}{{dW}} = \frac{{2\pi }}{\lambda }\left( {{L_W} \cdot \frac{{d{B_{bi}}}}{{dW}} + {B_{bi}} \cdot \frac{{d{L_W}}}{{dW}}} \right)$$
where ${L_W}$ is the length of the fiber under transverse force. The effect of transverse force on fiber length is negligible, that is $d{L_W}/dW = 0$, therefore when the transverse force changes, the wavelength shift can be expressed as
$$\frac{{d\lambda }}{{dW}} = {L_W} \cdot {\frac{{d{B_{bi}}}}{{dW}}} \left/ {\left( {L \cdot \frac{{d{B_{ai}}}}{{d\lambda }} - \frac{{{B_{ai}}L + {B_{bi}}{L_W}}}{\lambda }} \right)}\right.$$
then we define $F = W/{L_W}$ is the transverse force per unit length.

3. Experimental results and discussions

The fiber used in this experiment is a 3-mode (LP01, LP11a and LP11b) PANDA PM-FMF fabricated by Fiberhome. Figures 3(a), 4(a) and 5(a) show the shift of the inference spectrum with the parameter changes. We first conducted the sensing experiments with temperature as a single variable. Figure 3(b) shows the wavelength shift fitting curves of LP01 mode, LP11a mode and LP11b mode versus temperature, the spheres, stars and up triangles are the scatter points of the experimental results. These scatter points are linear fitted and three wavelength-temperature relationship curves of LP01 mode, LP11a mode and LP11b mode are formed respectively. We can see that the wavelength shift of each mode changes linearly with temperature, and the corresponding temperature sensing coefficients of each mode are ${k_{T01}} ={-} 0.\textrm{28693}$, ${k_{T11\textrm{a}}} ={-} 0.\textrm{28534}$, ${k_{T11\textrm{b}}} ={-} 0.\textrm{30364}$.

 figure: Fig. 3.

Fig. 3. (a) The transmission spectrum around the same transmission minimum of the LP01 mode at different temperatures; (b) Wavelength versus temperature.

Download Full Size | PDF

 figure: Fig. 4.

Fig. 4. (a) The transmission spectrum around the same transmission minimum of the LP01 mode at different strains; (b) Wavelength versus strain.

Download Full Size | PDF

We then conducted the sensing experiments with strain as a single variable, the corresponding strain sensing coefficients of each mode showed in Fig. 4(b) are ${k_{S01}} = 0.01001$, ${k_{S11a}} = 0.00830$ and ${k_{S11b}} = 0.01057$. Then we conducted the sensing experiments with transverse force as a single variable, and the slow and fast axis of the PM-FMF were aligned to the direction of force respectively. After testing we found that the fast axis is more sensitive than the slow axis, so in the multi-parameter sensing experiment the fast axis were chosen to be the direction of the force. Figure 5(b) shows the wavelength shift fitting curves of LP01 mode, LP11a mode and LP11b mode versus transverse force per unit length. The wavelength shift of each mode also changes linearly with transverse pressure, and the corresponding force sensing coefficients of each mode are ${k_{F01}} = 0.00\textrm{658}$, ${k_{F11a}} = 0.00629$ and ${k_{F11b}} = 0.00\textrm{603}$.

 figure: Fig. 5.

Fig. 5. (a) The transmission spectrum around the same transmission minimum of the LP01 mode at different transverse forces per unit length; (b) Wavelength versus transverse force per unit length.

Download Full Size | PDF

After measuring the sensing coefficients of three single parameter fields of each mode, we can write out the multi-parameter coefficients matrix $k = \left[ \begin{array}{l} - 0.2896\textrm{3}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} 0.01001{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} 0.00\textrm{658}{\kern 1pt} {\kern 1pt} \\ - 0.2853\textrm{4}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} 0.00830{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} 0.00\textrm{629}{\kern 1pt} {\kern 1pt} \\ - 0.3036\textrm{4}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} 0.01057{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} 0.00\textrm{603} \end{array} \right]$, and the changes of the three parameter ΔT, ΔS and ΔF can be calculated according to Eq. (2).

We conducted the sensing experiment of simultaneous changes in temperature, strain, and force to verify the accuracy of the multi-parameter sensor. Figure 6 shows the wavelength shift under different temperature, strain and transverse force combinations in the multi-parameter sensing experiment. The wavelength shifts of the three LP modes are $\Delta {\lambda _{\textrm{01}}} = \textrm{0}\textrm{.52}nm$, $\Delta {\lambda _{\textrm{11}a}} = 0.32nm$, $\Delta {\lambda _{\textrm{11}b}} = \textrm{0}\textrm{.24}nm$. According to the equation we can calculate that $\mathrm{\Delta }T = \textrm{4}.75{\kern 1pt} {\kern 1pt} ^\circ C$, $\mathrm{\Delta }S = \textrm{107}\textrm{.88}{\kern 1pt} {\kern 1pt} \mu \varepsilon$, $\mathrm{\Delta }F = \textrm{123}\textrm{.95}{\kern 1pt} {\kern 1pt} N/m$ respectively. The data obtained in the experiment and the results obtained according to the equation are roughly consistent, the deviation is caused by the resolution and accuracy of the experimental instrument. Table 1 shows the other two groups of multi-parameters sensing data, where ΔTS, ΔSS and ΔFS represent the parameter changes we set in the experiment. It can be seen from the table that the calculated changes of the temperature differ from the changes set in the experiment by -0.36 °C to -0.06 °C, the calculated changes of the strain differ from the changes set in the experiment by -0.11 µε to 2.24 µε and the calculated changes of the transverse force differ from the changes set in the experiment by -1.67 N/m to -1.05 N/m.

 figure: Fig. 6.

Fig. 6. The wavelength shift of LP01, LP11a and LP11b mode under different temperature, strain and transverse force.

Download Full Size | PDF

Tables Icon

Table 1. Experiment data of multi-parameter sensing

4. Conclusion

We propose a multi-parameter sensor using PM-FMF based on polarized interference and deduce the equation to inversely solve the change of external parameters according to the wavelength shift of the interference spectrum by using a 3-mode PANDA PM-FMF to measure the simultaneous changes of temperature, strain and transverse, we experimentally demonstrated the effectiveness of the measurement method. The multi-parameter sensing system we used in this paper is economical and practical, which requires less system operation, and the sensor has a temperature sensitivity of 0.3 nm/°C, a strain sensitivity of 0.01 nm/µε and a transverse force sensitivity of 0.0065 nm/(N/m).

Disclosures

The authors declare that there are no conflicts of interest related to this article.

References

1. A. D. Ellis, J. Zhao, and D. Cotter, “Approaching the non-linear Shannon limit,” J. Lightwave Technol. 28(4), 423–433 (2010). [CrossRef]  

2. Z. Yang, J. Zhao, N. Bai, E. lp, T. Wang, Z. Li, and G. Li, “Experimental demonstration of adaptive VFF-RLS-FDE for long-distance mode-division multiplexed transmission,” Opt. Express 26(14), 18362–18367 (2018). [CrossRef]  

3. N. Zhao, X. Li, G. Li, and J. M. Kahn, “Capacity limits of spatially multiplexed free-space communication,” Nat. Photonics 9(12), 822–826 (2015). [CrossRef]  

4. SÖ Arik, D. Askarov, and J. M. Kahn, “Adaptive frequency-domain equalization in mode-division multiplexing systems,” J. Lightwave Technol. 32(10), 1841–1852 (2014). [CrossRef]  

5. E. M. Ip and J. M. Kahn, “Fiber impairment compensation using coherent detection and digital signal processing,” J. Lightwave Technol. 28(4), 502–519 (2010). [CrossRef]  

6. P. Sillard, M. Bigot-Astruc, D. Boivin, H. Maerten, and L. Provost, “Few-mode fiber for uncoupled mode-division multiplexing transmissions,” 37th European Conference and Exhibition on Optical Communication, Geneva, pp. 1–3 (2011).

7. K. A. Murphy, M. S. Miller, A. M. Vengsarkar, and R. O. Claus, “Elliptical-core two mode optical-fiber sensor implementation methods,” J. Lightwave Technol. 8(11), 1688–1696 (1990). [CrossRef]  

8. A. Kumar, N. K. Goel, and R. K. Varshney, “Studies on a few-mode fiber-optic strain sensor based on LP01-LP02 mode interference,” J. Lightwave Technol. 19(3), 358–362 (2001). [CrossRef]  

9. X. L. Li, W. Zhang, J. Ruan, and S. S. Zhang, “Temperature- and strain-insensitive torsion sensor based on phase-shifted ultra-long-period grating,” Electron. Lett. 48(4), 235–236 (2012). [CrossRef]  

10. J. Li, W. Zhang, S. Gao, Z. Bai, L. Wang, H. Liang, and T. Yan, “Simultaneous force and temperature measurement using S fiber taper in fiber Bragg grating,” IEEE Photonics Technol. Lett. 26(3), 309–312 (2014). [CrossRef]  

11. A. Li, Y. Wang, J. Fang, M. Li, B. Y. Kim, and W. Shieh, “Few-mode fiber multi-parameter sensor with distributed temperature and strain discrimination,” Opt. Lett. 40(7), 1488–1491 (2015). [CrossRef]  

12. A. Li, Y. Wang, Q. Hu, and W. Shieh, “Few-mode fiber based optical sensors,” Opt. Exp 23(2), 1139 (2015). [CrossRef]  

13. X. Dong, H. Y. Tam, and P. Shum, “Temperature-insensitive strain sensor with polarization-maintaining photonic crystal fiber based Sagnac interferometer,” Appl. Phys. Lett. 90(15), 151113 (2007). [CrossRef]  

14. Y. Liu, B. Liu, X. Feng, W. Zhang, G. Zhou, S. Yuan, G. Kai, and X. Dong, “High-birefringence fiber loop mirrors and their applications as sensors,” Appl. Opt. 44(12), 2382–2390 (2005). [CrossRef]  

15. Y. Yang, J. Li, W. Duan, X. Zhang, W. Jin, and M. Yang, “An embedded pressure sensor based on polarization maintaining photonic crystal fiber,” Meas. Sci. Technol. 24(9), 094004 (2013). [CrossRef]  

16. G. Hu and D. Chen, “Side-hole dual-core photonics crystal fiber for hydrostatic pressure sensing,” J. Lightwave Technol. 30(14), 2382–2387 (2012). [CrossRef]  

17. L. Yu, J. Zhao, Q. Mo, and G. Li, “Simultaneous measurement of temperature and strain based on a polarization-maintaining few-mode fiber,” Optical Fiber Communications Conference and Exhibition,” Los Angeles, 1–3 (2017).

18. Y. Liu, Z. Han, J. Zhao, Z. Yang, L. Yu, and W. Wang, “Simultaneous beat-length measurement of a polarization-maintaining few-mode fiber,” IEEE Photonics J. 11(5), 1–6 (2019). [CrossRef]  

19. J. Li, L. Sun, S. Gao, Z. Quan, Y. Chang, Y. Ran, L. Jin, and B. Guan, “Ultrasensitive refractive-index sensors based on rectangular silica microfibers,” Opt. Lett. 36(18), 3593–3595 (2011). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (6)

Fig. 1.
Fig. 1. Experimental setup for multi-parameter sensing measurement. a: amplified spontaneous emission (ASE) light source; b: collimator; c: polarizer 1; d: phase plate; f: collimator; g: fiber rotator; h: heating panel; n: optical spectrum analyzer (OSA).
Fig. 2.
Fig. 2. The polarized interference spectrum of the PM-FMF and the wavelength shift as the temperature increases.
Fig. 3.
Fig. 3. (a) The transmission spectrum around the same transmission minimum of the LP01 mode at different temperatures; (b) Wavelength versus temperature.
Fig. 4.
Fig. 4. (a) The transmission spectrum around the same transmission minimum of the LP01 mode at different strains; (b) Wavelength versus strain.
Fig. 5.
Fig. 5. (a) The transmission spectrum around the same transmission minimum of the LP01 mode at different transverse forces per unit length; (b) Wavelength versus transverse force per unit length.
Fig. 6.
Fig. 6. The wavelength shift of LP01, LP11a and LP11b mode under different temperature, strain and transverse force.

Tables (1)

Tables Icon

Table 1. Experiment data of multi-parameter sensing

Equations (10)

Equations on this page are rendered with MathJax. Learn more.

A = k B
B = k  - 1 A
I i = A i 2 ( R 2 + S 2 + 2 R S cos ϕ i )
ϕ i = 2 π λ B a i L + 2 π λ B b i L E
d ϕ i d T = 2 π λ ( L H d B b i d T + B b i d L H d T )
d λ d T = L H d B b i d T / ( L d B a i d λ B a i L + B b i L H λ )
d ϕ i d S = 2 π λ ( L S d B b i d S + B b i d L S d S )
d λ d S = L S d B b i d S / ( L d B a i d λ B a i L + B b i L S λ )
d ϕ i d W = 2 π λ ( L W d B b i d W + B b i d L W d W )
d λ d W = L W d B b i d W / ( L d B a i d λ B a i L + B b i L W λ )
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.