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

Stability of a monolithic integrated filtered-feedback laser

Open Access Open Access

Abstract

We experimentally and theoretically investigate the stability of a single-mode integrated filtered-feedback laser as a function of the electrically controlled feedback phase. We interpret the measurements in terms of feedback-induced dynamics, compare them with the results from a stability analysis model for conventional feedback, and find good qualitative agreement.

©2012 Optical Society of America

1. Introduction

Single-wavelength lasers with narrow linewidth are required in many photonic application fields, e.g. in optical communication and sensing. One way to realize such lasers is by utilizing feedback from an external mirror. Although the introduction of delayed optical feedback can be disadvantageous as it may cause instabilities to the laser like the undamping of relaxation oscillations and coherence collapse, under suitable conditions it can improve particular features of the laser, such as enhanced side mode suppression and linewidth narrowing [1].

Many semiconductor external-cavity lasers with different feedback configurations have been reported, where most of them have external cavities outside the chip and the external cavities are much longer than the solitary cavity. Unlike the free-space or fiber-coupled feedback laser systems, the external cavity of an integrated feedback laser is on-chip and with small footprint. Compared to an off-chip cavity, the external cavity is at least two orders of magnitude shorter, implying better stability, better reproducibility from batch to batch and operation on faster time scales.

Here we present the results of a theoretical and experimental investigation of the stability of a monolithically integrated filtered-feedback laser with electrically controlled feedback phase, which is part of a multi-wavelength laser device. A relevant earlier study involving a multi-section laser with phase control is published in [2], but it applies to a coupled-cavity laser in the strong coupling limit, whereas our laser operates in the weak-feedback regime, where external cavity modes are non-existent. In our device, each laser channel operates in single wavelength, which has been achieved by the presence of an arrayed waveguide grating (AWG) filter within the external feedback section. The experimental stability study is obtained by the analysis of the optical and electrical spectra, while the theoretical study is based on a rate equation approach for a single-mode laser with conventional external feedback and leads to results that are in good agreement with experimental findings.

2. Device description

The laser we investigate here is one channel of an integrated filtered-feedback 4-channel multi-wavelength laser (MWL) with channel spacing of 200 GHz and operating in the C-band, as reported by some of us earlier [3]. The operation principle of this MWL can be understood from its schematic representation in Fig. 1(a) . The device consists of an array of four Fabry-Pérot (FP) lasers, each of which is formed by a semiconductor optical amplifier (SOA) and two on-chip broadband multi-mode interference reflectors (MIRs) [4] as FP-mirrors. The MIRs provide reflection and transmission to and from the laser cavity. The physical length of each FP-cavity is 820 mm. One side of each FP-laser is connected to one common arrayed-waveguide grating (AWG) filter through a phase shifter, all on a single chip (Fig. 1(b)). Each FP-laser generates a broad spectrum of longitudinal modes, of which the spacing is determined by the FP-cavity length. The AWG multiplexes the light of the FP-laser channels and connects to another on-chip broadband reflector. This reflector will reflect the light and subsequently the AWG will demultiplex and route the reflected signals back to each original FP laser. Only the FP-mode with wavelength within the corresponding AWG pass band will have sufficiently low feedback-reduced loss to suppress all other longitudinal modes through gain competition. Hence, each laser will be operating in a single mode and the wavelength spacing between the lasers is determined by the spacing of the AWG pass bands, 200 GHz in this case. Thus, the AWG together with the external reflector effectively forms an external feedback section for each FP laser, as shown in Fig. 2 . The phase shifting section (PHS, see Fig. 1(a)) allows to individually control the feedback phase for each FP-cavity.

 figure: Fig. 1

Fig. 1 (a) Schematic of the device and (b) Microscopy photograph of the realized device. The chip dimensions are approximately 2 × 2.5 mm.

Download Full Size | PDF

 figure: Fig. 2

Fig. 2 Schematic of one channel: a FP laser and the external cavity with phase shifter inside; On the bottom, sketches of spectra feature are shown: FP modes (left) and the AWG passband (right).

Download Full Size | PDF

Each FP-laser has 49.5 GHz mode spacing, and the two-port MIRs forming the cavity have power reflectivity and transmittivity of R = T = ~0.32. The feedback section is formed by a common 1-port MIR with ~80% reflectivity [4], a single AWG channel and a separate phase shifter for each channel. In the following example we will show recorded spectra for a laser channel that has an external cavity of ~3.0 mm physical length. The AWG is estimated to have 5 dB loss per pass. The forward-biased phase shifter is 0.35 mm long with negligible loss and phase-shift efficiency of 15 degree per mA per pass, which was measured from a test structure with a Mach-Zehnder configuration. Including the waveguide and AWG losses back and forth, the effective strength of the feedback is estimated to be 2% of the light transmitted to the external cavity.

3. Theory and model

The lasers under study are subject to filtered feedback [5]. The filtering is achieved by the presence of the AWG. The main purpose of the AWG is to provide slightly reduced loss for the longitudinal mode that falls within a certain channel, so as to force that laser to operate in one particular mode. As the effective width (3-dB passband) of the AWG pass band filter is ~50 GHz and the laser operates in a single mode, while the feedback- induced frequency change is below ~2 GHz, the feedback can be modelled as conventional feedback [6], which simplifies the treatment and leads to good qualitative agreement with observations.

We will express the optical field inside the laser cavity asP(t)eiϕ(t)+iω0t+c.c, where the powerPand the phaseϕ depend slowly on time on the scale of the rapid oscillations implied by the optical frequency ω0 of the solitary laser. The power is in dimensionless units, such that Pis the number of photons in the cavity. The variable characterizing the gain is the slowly time-dependent number of electron-hole pairsN(t), taken with respect to its threshold value for the solitary laser (i.e. without feedback). In terms of the three variablesP, ϕand N, the equations describing the time evolution of the laser in the presence of feedback are

P˙=ξNP+2γP(t)P(tτ)cos[ϕ(t)ϕ(tτ)+ω0τ];
ϕ˙=12αξNγP(tτ)P(t)sin[ϕ(t)ϕ(tτ)+ω0τ],
N˙=(1T1+ξP)NΓ0(PP0).
whereξ is the linearized gain coefficient at threshold; γis the feedback rate, defined as [6]
γ=rext(1r2)τinr
with rextthe (amplitude) reflectivity of the external mirror, r the (amplitude) reflectivity of the laser mirror and τinthe internal roundtrip time of the laser cavity. The parameter τis the feedback delay time and αthe linewidth-enhancement factor; T1is the effective carrier lifetime for spontaneous recombination at threshold, Γ0the cavity decay rate (inverse photon lifetime) andP0 the solitary laser power, given by [6]
P0=JJthrΓ0
where Jthr is the solitary threshold current (in units of electrons per second).

This completes the model for the dynamical description of the laser with feedback section. In the next section we will study the steady-state solutions to the rate Eqs. (1)-(3), followed in sec.3.2 by a stability analysis for the corresponding stationary state.

3.1 Stationary state

Single-frequency laser operation corresponds to solutions of Eqs. (1), (2) and (3) satisfying

P=Ps;N=Ns;ϕ=Δωst
wherePs,Nsand Δωsare time independent. This yields
Δωs=γ1+α2sin(arctanα+ω0τ+Δωsτ),
Ns=2γξcos(ω0τ+Δωsτ),
Ps=P0NsΓ0T11+ξNsΓ0.
First, the frequency shift Δωsis numerically solved from Eq. (7). It was shown by Ref [7]. that, so long the feedback strength satisfies Cγτ1+α2<1, only one single solution to Eq. (7) exists. Substituting this solution in Eq. (8) yields a value forNs, where after the value for Ps can be evaluated with Eq. (9). Thus we obtain for each setting of the laser, determined by the value of the feedback phase delay ω0τand the solitary laser power P0, a set of values forΔωs, Nsand Ps. It is important to note that the applied phase ϕps in the phase shifting element will be accounted for in ω0τas well (see below, at the end of Sec.3.2).

3.2 Stability

Now the issue at stake is that the stationary state defined by Eqs. (7), (8) and (9) is not always stable. It is known that the transition from stable to unstable operation is associated with a so-called Hopf bifurcation, which indicates the birth of an undamped relaxation oscillation (RO) [6,8]. If this happens, there are, among others, the following possible scenarios:

  • 1) A self-sustained harmonic RO has appeared, notable as oscillations in P, N and ϕ, at a frequency approximately given by
    νR=12πξ(JJthr)

    with typical values ranging from 1 to 10 GHz.

  • 2) The oscillation period may double (frequency halves), often referred to as period-doubling bifurcation [9].
  • 3) As the amplitude of the RO increases, higher harmonics appear.
  • 4) The oscillations turn chaotic; in case the chaos is fully developed this state of operation is also referred to as coherence collapse [10].

The stability of the state characterized by Eqs. (7)-(9) is investigated by considering small deviations, i.e. putting

ϕ=Δωst+δϕ;N=Ns+δN;P=Ps+δP.
After substituting Eq. (11) in Eqs. (1)-(3) and Eqs. (7)-(9) we obtain, in lowest order, three coupled first-order differential delay equations for δϕ,δN and δP, which are analyzed for solutions with time dependence est for complex-valued s. This leads to a system determinant with zeroes in the complex s-plane. If at least one zero has a positive real part, the corresponding solution Eqs. (7)-(9) is unstable. This has been investigated using the principle of the argument. In this manner, we have obtained a stability diagram for a laser with delayed feedback indicating whether the laser is stable or unstable as functions of the applied feedback phase ϕps (between 0 and 2π) and the applied pump strength p, where the latter is defined asp(JJthr)/Jthr.

The stability diagram is displayed in Fig. 5(b) of Sec.5.2. Here, the total effective feedback phase has been taken as ω0τ=ϕps+ωipIthrτ, where ωi is the (empirical) injection-current-induced frequency shift per mA injection current and Ithris the threshold injection current (in mA), i.e. Ithr=103qJthr, with q the unit charge. We clearly observe tilted bands of stability and instability, respectively, which connect at2π. The negative slope of the bands is due to the negative current-induced frequency shift of the solitary laser (see the text pertaining to Fig. 5 in Sec.4).

4. Characterization of the device

The measurement setup is depicted in Fig. 3 . The device we investigated is one channel of the integrated filtered-feedback multi-wavelength laser with external cavity length ~3.0 mm. The laser chip is mounted on a copper chuck and its temperature is stabilized at 18 °C during the measurements. At this temperature the laser has a threshold current of 16 mA. The electrical currents applied to the SOA and to the phase shifter, respectively, are provided through two probe needles. A lensed fiber mounted on a nano-positioning stage is used to collect the light. An optical isolator is used to avoid the unwanted back reflection into the laser chip. The lightis analyzed using a 50 GHz bandwidth photodetector with an 8 GHz 30 dB gain low noise amplifier connected to a 50 GHz electrical spectrum analyzer [Agilent E4448A]. A fraction of the light is sent to an optical spectrum analyzer [APEX AP2041B] with a high resolution of 0.16 pm.

 figure: Fig. 3

Fig. 3 The measurement setup. Ipump is the laser injection current, and Iphase is the current on the phase shifter to tune the feedback phase; TEC denotes the temperature controller; ISO is an optical isolator. After the 50/50 coupler, one of the two branches of light is fed to a 50 GHz photodiode, amplified by a low noise amplifier (LNA) and measured with an electrical spectrum analyzer, while the other branch connects to an optical spectrum analyzer.

Download Full Size | PDF

4.1 Spectral lasing regimes

First, the phase shifter is unbiased, and the SOA injection current is varied until the laser shows a stable single wavelength lasing mode. Then keeping the injection current fixed to that value, the feedback phase shift is slowly swept from 0° to 210° by steps of 15°. The optical signal and the noise intensity signal are recorded only when there is a clear change in the signal appearance, and presented in Fig. 4 . Strong dynamics have been observed for certain feedback phase conditions, and the size of this region varies with the injection currents level applied to the SOA.

 figure: Fig. 4

Fig. 4 Intensity noise (dashed lines) and optical spectral (solid lines) from the ESA and OSA, respectively, for different feedback phase tuning 0°, 90°, 150°, 180°, 210°.

Download Full Size | PDF

The results presented in Fig. 4 are obtained by applying a fixed current of 72 mA to the SOA. When there is no extra phase shift applied to the feedback, the laser is very stable emitting single wavelength with narrow linewidth, where the latter could not be resolved by the 0.16 pm resolution of the OSA (<20 MHz). The electrical spectrum (Fig. 4(a)) presents a flat and relatively low response. The relaxation oscillation (RO) frequency is damped below the noise level of the spectrum analyzer. With the feedback phase tuned to 90° (Fig. 4(b)), the laser is not stable anymore. The RO becomes undamped and it is identified as a narrow peak in the electrical spectrum at 6.8 GHz, while in the optical spectrum two series of side peaks appear at distances equal to multiples of the RO-frequency on both sides of the central peak. When the phase is increased (Fig. 4(c)), the laser enters a period-doubling regime, identified by more side peaks in the optical spectrum and a new peak at half the RO frequency in the electrical spectrum. In the next three cases, Figs. 4(c)-4(e), the RO has developed into a more complicated and ultimately irregular oscillation indicating period doubling and higher-harmonic effects, but with the dominant frequency near the RO frequency. Gradually, after entering the instability region, the noise floor increases significantly, as seen in the noise intensity spectra and optical spectra. The last case presented in Fig. 4(e) (φ = 210°) shows coherence collapsed operation [10] dominated by relaxation oscillations. The different lasing regimes with different feedback phase are very reproducible.

To get a clearer picture of the overall stability of our filtered-feedback laser, the same phase tuning procedure is repeated for all injection currents above threshold, with a step of 2 mA. The situation as in Fig. 4(a), single wavelength lasing with no relaxation oscillation dynamics, is considered as stable. We have identified such stable regions and made a stability map plotted in Fig. 5(a) . The map displays the laser stability as a function of the feedback phase (horizontal) and injection current (vertical), where white indicates stability and black indicates instability, which includes all other unstable regimes. With injection current increase, the laser demonstrates different stable regions. The location of these regions changes with changing feedback phase. Unfortunately, the full 2π phase shift could not be achieved, because the phase shifter starts to saturate at currents beyond ~7 mA and is fully saturated at 10 mA (Fig. 5(a)), thus limiting the maximum feedback tuning to 210 °.

 figure: Fig. 5

Fig. 5 (a) Measured map of laser stability versus phase shifter current with different injection current. The dark region is the unstable region and the bright region is the stable region. (b) Theoretical stability diagram for a single-mode laser vs. pump strength p (vertical) and feedback phase (horizontal), with p defined as p=(JJthr)/Jthr. Color white denotes stable operation; black denotes unstable. Parameters are: feedback rate5×109s1; feedback delay time 7×1011s and α=2.6; the injection-current-induced frequency shift of the solitary laser is −0.6 GHz/mA. The most left column gives the relaxation oscillation frequency in GHz.

Download Full Size | PDF

4.2 Stability analysis

As was mentioned in Sec.3.2, we have applied a numerical stability analysis to the laser with feedback by investigating the dynamical system determinant. Using the value α = 2.6 for the linewidth enhancement factor, as obtained from injection locking measurements (see Appendix A), and the external cavity length and loss estimated as above, we find for the case of injection current of 72 mA (pump strength p~3) that over a phase range interval of about 100° the single frequency operation is stable, whereas it is unstable for other values of the applied feedback phase. Apparently, Fig. 4(a) is within this stability interval.

Figure 5(b) is the stability map generated from the theoretical model. The stability situation is recorded as a function of the feedback phase (0-2π; horizontal axis) and pump strength (injection current from threshold up to 80 mA, p = 4). The parameter values used for this diagram are given in the caption. It is to be noted that zero phase on the horizontal scale in Fig. 5(b) is arbitrary; it does not necessarily coincide with zero current in Fig. 5(a). Comparing this with the measured map Fig. 5(a), the structures of stable and unstable bands with negative slope show good qualitative agreement between experiment and theory, although theory predicts wider unstable bands. The slope is due to the injection current induced phase shift inside the FP laser cavity.

5. Discussion and conclusion

An integrated laser with a filtered-feedback configuration has been realized. The external cavity is integrated on the chip and the phase of the feedback can be controlled via an integrated phase shifter. We have measured the laser feedback dynamics as a function of the feedback phase and injection current. By changing the feedback phase, the laser can be set to operate in a stable single-frequency state as well as in various unstable states. Comparison of the measurements with the results from a stability model, confirms that the theory is applicable to a laser with short and integrated feedback. We are able to simulate the stability of the device with respect to pump current, feedback phase and delay time. This will allow further optimized designs for integrated filtered-feedback lasers that will show less or no feedback phase sensitivity, thus improving the laser stability and simplicity of the circuit.

Appendix: linewidth-enhancement factor

The linewidth-enhancement factor (α) is one of the key parameters for semiconductor lasers. This factor is also used in the model discussed in section 3, where it is taken as 2.6. This appendix will describe how the value 2.6 for α is measured. The determination of the α is based on measuring the injection-locking range of the laser [11]. The FP-laser in this measurement is a test laser with identical design of the solitary laser in the filtered-feedback multi-wavelength laser, and is fabricated under the same technology.

The measurement setup is shown in Fig. 6 . An optical spectrum analyzer (OSA, APEX) is used to measure the optical spectra of the test laser and the master laser. The tunable master laser (Agilent 86100B) and the MIR based FP laser were both temperature stabilized during the measurements. The tunable laser was set to a certain output power, and was attenuated by a variable optical attenuator. After the attenuator, the light was connected to a polarizing beam splitter (PBS) through a polarization controller (PC). The PC tunes the TE output from the PBS to the maximum, and the TE light is injected into the FP laser via a polarization-maintaining (PM) fiber. On the other side of the chip, a lensed fiber is aligned to the output facet of the FP laser to couple out the output light, and the light is sent to an APEX optical spectrum analyzer (OSA) for spectral analysis of the emission.

 figure: Fig. 6

Fig. 6 Injection locking measurement setup. DUT is the device under test. See the text for other details.

Download Full Size | PDF

The locking range of this FP laser at 35mA pump current for a longitudinal mode around 1545 nm is measured. The locked region is determined by both frequency detuning and the external injection light power at a fixed laser pump current. At each injecting power, the wavelength is swept from below and above the frequency of this mode, and the frequencies when the mode loses locking are recorded and plotted in Fig. 7 versus the square root of the injection power. The boundaries of the locking region are fitted by the two lines. Outside this area is the unlocked is unlocked region, where the injection level is too low or the detuning is too high to reach injection-locking condition. According to [11,12], the linewidth enhancement factor can be obtained as α=(Δfmax/Δfmin)21.The ratio of the two slopes of the boundaries of the locking region is independent of the optical powers. The linewidth enhancement factor of the active material can be estimated to be 2.6.

 figure: Fig. 7

Fig. 7 Frequency locking range dependence on injected power. FP laser is under injection current of 35 mA at 18 °C.

Download Full Size | PDF

References and links

1. C. H. Henry, “Phase noise in semiconductor laser,” J. Lightwave Technol. 4(3), 298–311 (1986). [CrossRef]  

2. S. Bauer, O. Brox, J. Kreissl, B. Sartorius, M. Radziunas, J. Sieber, H. J. Wünsche, and F. Henneberger, “Nonlinear dynamics of semiconductor lasers with active optical feedback,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 69(1), 016206 (2004). [CrossRef]   [PubMed]  

3. J. Zhao, P. J. Williams, M. K. Smit, and X. J. M. Leijtens, “ Monolithic integrated filtered-feedback multi-wavelength laser,” in Proc. Optical Fiber Communication Conf. (OFC 2012), Los Angeles, USA, Mar. 4–8, 2012, Paper OW1G.5.

4. E. Kleijn, M. K. Smit, M. J. Wale, and X. J. M. Leijtens, “New two-port multimode interference reflectors,” in Proc. 16th Eur. Conf. Int. Opt. (ECIO ’12), Sitges-Barcelona, Spain, Apr.18–20, (2012).

5. M. Yousefi and D. Lenstra, “Dynamical behavior of a semiconductor laser with filtered external optical feedback,” IEEE J. Quantum Electron. 35(6), 970–976 (1999). [CrossRef]  

6. G. H. M. Van Tartwijk, “Semiconductor lasers with optical injection and feedback,” Thesis Amsterdam, 76–80&88–90 (1994); see also J. Quantum Semiclass. Opt. 7, 87–143 (1995).

7. G. A. Acket, D. Lenstra, B. H. Verbeek, and A. J. Den Boef, “The influence of feedback intensity on longitudinal mode properties and optical noise in index-guided semiconductor lasers,” IEEE J. Quantum Electron. 20(10), 1163–1169 (1984). [CrossRef]  

8. J. Mork, B. Tromborg, and J. Mark, “Chaos in semiconductor lasers with optical feedback: Theory and Experiment,” IEEE J. Quantum Electron. 28(1), 93–108 (1992). [CrossRef]  

9. S. H. Strogatz, Nonlinear dynamics and chaos: with applications to physics, biology, chemistry, and engineering, (Perseus Books, Addison-Wesley Publishing, 1994).

10. D. Lenstra, B. H. Verbeek, and A. J. Den Boef, “Coherence collapse in single-mode semiconductor lasers due to optical feedback,” IEEE J. Quantum Electron. 21(6), 674–679 (1985). [CrossRef]  

11. G. Liu, X. Jin, and S. L. Chuang, “Measurement of inewidth enhancement factor of semiconductor lasers using an injection-locking technique,” IEEE Photon. Technol. Lett. 13(5), 430–432 (2001). [CrossRef]  

12. I. Petitbon, P. Gallion, G. Debarge, and C. Chabran, “Locking bandwidth and relaxation oscillations of an injection-locked semiconductor laser,” IEEE J. Quantum Electron. 24(2), 148–154 (1988). [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 (7)

Fig. 1
Fig. 1 (a) Schematic of the device and (b) Microscopy photograph of the realized device. The chip dimensions are approximately 2 × 2.5 mm.
Fig. 2
Fig. 2 Schematic of one channel: a FP laser and the external cavity with phase shifter inside; On the bottom, sketches of spectra feature are shown: FP modes (left) and the AWG passband (right).
Fig. 3
Fig. 3 The measurement setup. Ipump is the laser injection current, and Iphase is the current on the phase shifter to tune the feedback phase; TEC denotes the temperature controller; ISO is an optical isolator. After the 50/50 coupler, one of the two branches of light is fed to a 50 GHz photodiode, amplified by a low noise amplifier (LNA) and measured with an electrical spectrum analyzer, while the other branch connects to an optical spectrum analyzer.
Fig. 4
Fig. 4 Intensity noise (dashed lines) and optical spectral (solid lines) from the ESA and OSA, respectively, for different feedback phase tuning 0°, 90°, 150°, 180°, 210°.
Fig. 5
Fig. 5 (a) Measured map of laser stability versus phase shifter current with different injection current. The dark region is the unstable region and the bright region is the stable region. (b) Theoretical stability diagram for a single-mode laser vs. pump strength p (vertical) and feedback phase (horizontal), with p defined as p=(J J thr )/ J thr . Color white denotes stable operation; black denotes unstable. Parameters are: feedback rate 5× 10 9 s 1 ; feedback delay time 7× 10 11 s and α=2.6 ; the injection-current-induced frequency shift of the solitary laser is −0.6 GHz/mA. The most left column gives the relaxation oscillation frequency in GHz.
Fig. 6
Fig. 6 Injection locking measurement setup. DUT is the device under test. See the text for other details.
Fig. 7
Fig. 7 Frequency locking range dependence on injected power. FP laser is under injection current of 35 mA at 18 °C.

Equations (11)

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

P ˙ =ξNP+2γ P(t)P(tτ) cos[ϕ(t)ϕ(tτ)+ ω 0 τ];
ϕ ˙ = 1 2 αξNγ P(tτ) P(t) sin[ϕ(t)ϕ(tτ)+ ω 0 τ],
N ˙ =( 1 T 1 +ξP)N Γ 0 (P P 0 ).
γ= r ext (1 r 2 ) τ in r
P 0 = J J thr Γ 0
P= P s ;N= N s ;ϕ=Δ ω s t
Δ ω s =γ 1+ α 2 sin(arctanα+ ω 0 τ+Δ ω s τ),
N s = 2γ ξ cos( ω 0 τ+Δ ω s τ),
P s = P 0 N s Γ 0 T 1 1+ ξ N s Γ 0 .
ν R = 1 2π ξ(J J thr )
ϕ=Δ ω s t+δϕ;N= N s +δN;P= P s +δP.
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.