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

Study of differential cross-polarization modulation in a semiconductor optical amplifier

Open Access Open Access

Abstract

Differential cross-polarization modulation (DXPoM) wavelength conversion using a semiconductor optical amplifier has been proposed to outperform substantially traditional cross-polarization modulation (XPoM). This work presents analytical small-signal models of XPoM and DXPoM. The transfer function of DXPoM is used to elucidate the significant improvement in bandwidth and the relationships among the modulation bandwidth, the delay, the operating points and other device parameters of a semiconductor optical amplifier (SOA).

©2005 Optical Society of America

1. Introduction

All optical wavelength conversion (AOWC) is considered to be an indispensable function in next-generation wavelength routing network architecture [1]. AOWCs based on semiconductor optical amplifiers (SOAs) are particularly attractive because they provide the advantage of monolithic integration and the potential of low-cost production. Various SOA-based AOWCs have been proposed, such as cross-gain modulation (XGM) [2], cross-phase modulation (XPM) [2], cross-polarization modulation (XPoM) [3], four-wave mixing (FWM) [4] and delay-interference signal-wavelength converter (DISC) [5,6]. Based on an interferometric principle that is similar to that exploited by the Mach-Zehnder interferometer (MZI), the XPoM employs the phase difference between the TE and the TM modes as an optical beam passes through an SOA. Accordingly, an additional delay line in the MZI-based XPM [7,8] also works in the XPoM scheme. Restated, simply adding an extra birefringence delay line to the standard XPoM considerably improves the conversion performance. This novel scheme is called differential cross-polarization modulation (DXPoM) [9]. Although the architecture of DXPoM is similar to that of DISC, DXPoM utilizes the interference between two distinct and orthogonal modes achieving wavelength conversion even without an extra delay line. However, DISC employs self-interference, and works only for the RZ format.

Although large-signal simulations have already established the difference between the performance of XPoM and that of DXPoM, and yielded results that agree with the experimental results [9], this study presents the analytical small-signal model to explain the significant performance improvement offered by DXPoM, and to clarify the relationships among modulation bandwidth, the delay, the operating points and other device parameters of an SOA. Furthermore, the optimum delay which corresponds to not only a largest conversion bandwidth but also a least phase distortion is derived analytically from the small-signal model.

2. Implementing DXPoM

Figure 1 presents the configuration of DXPoM. A signal pump laser beam at wavelength λ1 and a CW probe laser beam at wavelength λ2 are fed into an SOA, as in a typical XPoM. Properly controlling the polarization states of λ1 and λ2 allows the injected pump light to introduce additional birefringence in the SOA, resulting in a change in the difference between the refractive indices of the TE and TM modes of the probe beam. At the polarizer, these two orthogonal modes are partially combined coherently. Restated, the XPoM exploits the phase difference between the TE and TM modes, when the probe beam passes through an SOA, to rotate the polarization state. This phase difference, controlled by the signal power, determines the output power of the CW beam after it has passed through a polarizer. However, an extra birefringence delay line is added in front of the polarizer in the DXPoM to improve the conversion speed of the traditional XPoM.

 figure: Fig. 1.

Fig. 1. Configuration of a DXPoM wavelength converter

Download Full Size | PDF

 figure: Fig. 2.

Fig. 2. Illustration of a differential sinusoidal wave

Download Full Size | PDF

A simplified and intuitive example helps to explain why this extra delay line improves the conversion speed. Figure 2 shows two sinusoidal waves with the same angular frequency, Ω, but different amplitudes, φa and φb. The difference between these two waves is also a sinusoidal wave with frequency, Ω, and amplitude, φc. If the time offset, ∆t, as presented in Fig. 2 is applied, then φc2 =φa2 +φb2 - 2φaφb cos(Ω∆t) . For Ω ≤ π/∆t, the differential amplitude increases with frequency. Namely, properly selecting ∆t increases the differential amplitudes at some frequencies. Accordingly, adding an extra delay between the TE and TM modes may amplify some high-frequency components of the phase difference in DXPoM, to compensate for the insufficient frequency response associated with the long lifetime of the carriers in an SOA. The polarization state of the output CW beam after the delay line is passed is rotated more rapidly as the signal power varies. Consequently, DXPoM has a higher conversion speed and a better performance than XPoM.

3. Small-signal models of XPoM

Based on the equations derived in references [10,11], the propagating equations of the optical fields in an SOA are,

dikdz=[Γkaik(ηkNNi)αlossk]ik,k=TE,TM
dNdt=Nτsi,k[aik(ηkNNi)ikħωiAeff]+IeV

where i =1, 2 represent the signal and the CW beam, respectively. ωi and Ni are the optical angular frequency and the carrier density at transparency of the i beam; ik and aik are the optical intensity and the differential gain coefficient of the i beam polarized in k mode; Γk and αlossk are the confinement factor and the waveguide loss of an SOA in k mode; N is the carrier density; Aeff is the effective area of the waveguide; ħ is the reduced Planck constant; I is the injected current; e is the elementary charge; V is the active volume, and τs is the lifetime of the carriers governed by spontaneous emission and non-radiative recombination. ηk is the modified imbalance factor used to describe the asymmetry of optical transitions between the TE and TM modes when the tensile strain is built into the active layer of an SOA. If the tensile strain is not applied, then ηk =1. Using Eq. (1), which neglects waveguide loss, and integrating both sides of Eq. (2) with respect to z yield,

dσdt=σσ0τsi,k0,ikΓkħωiAeff·{exp[Γkaik(ηkσNi)]1}

where σ 0 = τsI/(eAeff ) and z,ik = ik (z) · σ = 0 N(t,z)dz and ℓ are the integrated carrier density and the length of an SOA, respectively. The phase of the output converted signal is,

ψ2k=αH2Γka2k[ηkσN2]

where αH is the linewidth enhancement factor.

The optical power, the integrated carrier density and the output phase associated with harmonically modulated input optical fields, can be expressed as z,ik = Pz,ik + ∆Pz,ik e jΩt + c.c, σ = ζ + ∆ζ e jΩt + c.c and ψ2k = ϕ2k+∆ϕ2k e jΩt + c.c, respectively. The CW beam exhibits no input modulation and ∆Pk 0,2 =0. Considering only the first-order terms in Eq. (3) yields the small-signal response of the integrated carrier density,

Δζ=k(G1k1)Γka1kηkPsat,1kΔP0,1kjΩτs+1+i,kGikP0,ikPsat,ik

where Gik = exp[Γk aik (ηk ζ-Ni)] and Psat,ik = ħωiAeff /(τsaikηk ) . Similarly, the small-signal response of the intensity and phase of the converted CW beam are,

ΔP,2k=G2kP0,2kΓka2kηkΔζ
Δϕ2k=αH2Γka2kηkΔζ

Notably, Eq. (6) is the small-signal response of XGM. In XPoM, the amplitude conversion is associated with interference between the TE and TM modes. Therefore, the output power of the CW beam behind a polarization controller (PC) and a polarizer, ℘Pol2, is given by,

pol,2=12cosθ·,2TEsinθ·,2TM·ej(Φ+ΔΦ)2

where θ and Φ are the parameters that specify the polarization state, which could be adjusted by a PC. ∆Φ = ∆ϕ2TM -∆ϕ2TE is the relative phase difference between the TE and TM modes of the output CW beam. To maximize the extinction ratio (ER) after the polarizer, θ must be chosen to ensure that the TE and TM modes beyond the polarizer have equal DC intensity, such that cos2 θ · P,2TE = sin2 θ · P,2TM = Pθ and Φ is selected to ensure that ℘pol,2 has a minimum zero level. For example, Φ + ∆Φ(0) = 0 should be applied for the non-inverting operation, where ∆Φ(0) is the phase difference when the input signal is logical zero. Thus, under these conditions, the small-signal response of XPoM, ∆P Pol,2, given by Eq. (8), could be written as,

ΔPpol,2=PθsinΦ·ΔΦ+12(1cosΦ)(cos2θ·ΔP,2TE+sin2θ·ΔP,2TM)

The first and second terms on the right-hand side of Eq. (10) are associated with XPM and XGM, respectively. Equations (5)–(7) and (9) yield,

ΔPpol,2=[χ(1+rme)+(1rme)]sinΦjΩ+τs1+τstim,11+τstim,21×αH2τsk[ΓTEa2TEηTEPθΓka1kηkPsat,1k(G1k1)ΔP0,1k]

where τstim,i1kGikP0,ik(τsPsat,ik) is the lifetime of the carriers governed by the stimulated emission of the i beam; χ ≡ (1-cosΦ)/(αH sinΦ) = αH1tan(Φ/2) is a parameter that is related to the operating point, Φ; rme = (ΓTM a2TM ηTM )/(ΓTE a2TE ηTE ) represents the strength of the dependence of an SOA on the polarization, and is smaller than unity in general. In Eq. (10), the factor defined as ξ ≡ [χ(1+rme )+(1-rme )]sinΦ influences conversion logic and the efficiency of XPoM. Therefore, ξ > 0 and ξ < 0 represent the implementation of XPoM according to an inverted and a non-inverted conversion scheme, respectively. However, Eqs. (6) and (10) reveal that the small-signal bandwidths of XGM and XPoM are the same, and are constrained by the lifetime of the carriers.

4. Small-signal models of DXPoM

The small-signal response of DXPoM can be obtained by simply modifying Eq. (10). If the delay, ∆t, is added to the TM mode, then Eq. (10) becomes,

ΔPpol,2={[χ(1+rmeejΩΔt)+(1rmeejΩΔt)]jΩ+τtot1}
×{αHsinΦ2τsk[ΓTEa2TEηTEPθΓka1kηkPsat,1k(G1k1)ΔP0,1k]}

Where τtot1τs1+τstim,11+τstim,21 . Clearly, the first bracket in Eq. (11) determines the modulation bandwidth of DXPoM. Therefore, this term divided by (χ+1)τtot is defined as a transfer function T(Ω):

Γ(Ω)1γejΩΔt1+jΩτtot

where γ=rme (1-χ)/(1+χ) . Equation (12) can be rewritten as a normalized transfer function to take into account explicitly the improvement in the bandwidth of DXPoM:

TN(Ω')1γe'Δt'(1+jΩ')(1γ)

where Ω' = Ωτtot and ∆t' = ∆t/τtot are the normalized angular frequency and the normalized delay, respectively.

Although Eq. (13) cannot be solved analytically to yield the modulation bandwidth with arbitrary delay, the minimum and maximum bandwidths that correspond to ∆t' = 0 and Ω'|∆t'| = π can be derived as,

Ω'n,min=n21
Ω'n,max=(n1+γ1γ)21n1+γ1γ

where n is the parameter used to describe the 1/n bandwidth, Ω'n. Equation (15) shows the improvement in bandwidth depends on γ, which is related to the operating point and the properties of an SOA. However, the minimum bandwidth represents the traditional XPoM, in which the corresponding bandwidth of Eq. (14) is independent of γ.

The γ value of an SOA determines the performance of the transfer function, so a numerical simulation is carried out to illustrate the concept of γ value. The parameters presented in Table 1 and the basic propagation equations, Eqs. (1) and (2), indicate the polarization-dependent gain (PDG) is 1.3 dB; τtot is around 30 ps and the values of γ are 0.882 and 0.852, for non-inverted and inverted conversions, respectively. Additionally, no tensile strain is applied on the SOA, so ηk are unity. Then, PDG is 8 dB, and γ are 0.724 and 0.665 for non-inverted and inverted conversions, respectively.

Tables Icon

Table 1. Parameters Definitions and Their Values

Figure 3 depicts the normalized 1 dB bandwidth, n =1.26, with various values of γ. As γ approaches unity, Ω'1-dB,max becomes greater and the corresponding normalized delay, |&t'|, decreases. Moreover, the square marks in Fig. 3 clearly show that the approximation in Eq. (15) is reasonable. The bandwidth is maximized by choosing γ close to unity, indicating that rme and χ must be close to one and zero, respectively. Therefore, the SOA exhibits minimum PDG. However, according to Eq. (11), this will reduce the conversion efficiency. Therefore, a trade-off exists between maximizing the bandwidth and maximizing the conversion efficiency. Furthermore, compared with a non-inverted conversion scheme, an inverted conversion scheme always suffers less bandwidth improvement, but has larger conversion efficiency due to the different values of γ governed by different operating points

 figure: Fig. 3.

Fig. 3. Normalized 1 dB angular bandwidth of DXPoM with different values of γ.

Download Full Size | PDF

 figure: Fig. 4.

Fig. 4. Amplitude and delay responses of XPoM, DXPoM with TM delay and DXPoM with TE delay and corresponding simulated eye-diagrams

Download Full Size | PDF

When ∆t is positive (negative), DXPoM operates with a TM (TE) delay. Figure 4 compares T(Ω) = |T(Ω)|exp(jΘ(Ω)) of the TE delay with that of the TM delay, with γ = 0.85, τtot = 30 ps and |∆t| = 6.1 ps, corresponding to a maximum 1 dB bandwidth in Eq. (15). Figure 4 demonstrates that although DXPoMs with TM and TE delay exhibit identical amplitude responses, they exhibit different delay responses, -dΘ/dΩ. The TE delay has the worst phase response which drastically distorts the converted signal. This fact explains why DXPoM with TE delay underperforms [9]. DXPoM with TM delay not only has a larger bandwidth but also a flatter delay response than XPoM. The 10 Gbps PRBS (pseudo random binary sequence) input signal spectrum is directly multiplied by the transfer functions plotted in Fig. 4 to elucidate the effects of the phase response on the signal distortion. The results after inverse Fourier transform back to the time domain are illustrated in the insets in Fig. 4. The eyes closure is clearly observed in the TE delay because of the poor phase response.

5. Optimum delay of DXPoM

Figure 4 reveals that the conversion performance is determined by not only the bandwidth of the amplitude response but also the linearity of the phase response, which corresponds to the flatness of the delay response. A natural question is how to obtain the optimum delay, ∆t'opt, at which DXPoM exhibits the least amplitude and the least phase distortion. This question is answered firstly by considering the relationship between ln|TN (Ω')| and ΘN(Ω'), which are the real and imaginary parts of ln(TN (Ω')), respectively. First, s = jΩ′ is set and all of the zeros of TN (s) in the complex s plane are found as zT = (ln γ)/∆t′ + j · (2m π)/∆t′, where m are all integrals. When ∆t' > 0, all zeros are in the left half of the s plane and this system is a minimum phase system (MPS) [12]. An interesting characteristic of an MPS is that the amplitude and phase responses are Hilbert transforms of each other. It has been demonstrated that a constant logarithmic amplitude response implies a linear phase response, as these responses are a Hilbert transform pair [12]. Hence, only the flatness of the amplitude response of TN (Ω') must be considered in optimizing the extra delay. Approximating exp (-jΩ' ∆t') in Eq. (13) as 1-jΩ'∆t'-(Ω'∆t')2/2 and forcing |TN (Ω')|=1, yield the optimum delay,

Δt'opt=1γγ

which corresponds to simultaneous a flat amplitude response and a linear phase response. According to Eq. (15), this delay corresponds to a maximum bandwidth of n = πγ/(1+γ). Figure 5 plots the amplitude and delay responses at delays of ∆t'opt, 0.95∆t'opt and 1.05 ∆t'opt. The dotted lines in Fig. 5 represent the second-order approximations of the amplitude responses. The good match at low frequency indicates that this approximation can be applied to determine the optimum delay. Furthermore, a nearly constant delay response and the flattest amplitude response are achieved simultaneously by applying ∆t'opt.

6. Time domain performance

Much information concerning the time domain can be obtained from the impulse response, which is the inverse Fourier transform of the frequency response, as well as from the frequency response, itself, of the small-signal model. The inverse Fourier transform and Eq. (13) yield the following impulse response.

𝔽1{T(Ω)2π}=et'u(t')γ·e(t'Δt')u(t'Δt')

where u(t′) is a unit step function and t' is the normalized time. Figure 6 plots several impulse responses with different values of γ and ∆t'. Figure 6(a) presents the impulse response of XPoM without an extra delay, and Fig. 6(b) presents that of DXPoM with the optimum TM delay, as determined by Eq. (16). Figure 6(b) clearly presents a narrower response, because the extra delay cancels the relaxation tail that is associated with a finite lifetime of the carriers. However, increasing the delay beyond the optimum value may increase the response time, as shown in Fig. 6(c). In Fig. 6(d), γ is closer to unity than in Fig. 6(b), so the response is narrower and the delay required optimizing the response time is shorter. These phenomena are consistent with the tendency in Fig. 3. Furthermore, unlike Fig. 6(b), Fig. 6(e) illustrates the case of TE delay. Although the TE delay and the TM delay are associated with the same bandwidth in Fig. 4, the impulse responses of Figs. 6(b) and 6(e) are such different due to their different phase performance. Restated, the impulse response of Fig. 6(e) is so poor that the TE delay cannot improve the conversion performance.

 figure: Fig. 5.

Fig. 5. Amplitude and phase responses with delay of ∆t'opt , 0.95 ∆t'opt and 1.05 ∆t'opt.

Download Full Size | PDF

7. Conclusion

This study presents small-signal models of XPoM and DXPoM for the first time. In these models, the relationships between the bandwidth improvement provided by DXPoM and several parameters are readily observed. The substantial difference between the conversion performance of TM delay and that of TE delay is explained by the phase responses of the small-signal model. The optimum delay that corresponds to both the flattest amplitude and the optimal delay responses is analytically determined. In addition to those in the frequency domain, the impulse responses of the DXPoM in the time domain are derived to provide further insight into the performance improvement.

 figure: Fig. 6.

Fig. 6. Impulse responses of several conditions.

Download Full Size | PDF

Acknowledgment

The authors would like to thank the National Science Council, Republic of China, Taiwan for financially supporting this research under Contract No. NSC 93-2215-E-009-027.

References and links

1 . S. J. B. Yoo , “ Wavelength conversion technologies for WDM network applications ,” J. Lightwave Technol. 14 , 955 – 966 ( 1996 ). [CrossRef]  

2 . T. Durhuus, B. Mikkelsen, C. Joergensen, S. L. Danielsen, and K. E. Stubkjaer , “ All-optical wavelength conversion by semiconductor optical amplifiers ,” J. Lightwave Technol. 14 , 942 – 954 ( 1996 ). [CrossRef]  

3 . Y. Liu, M. T. Hill, E. Tangdiongga, H. de Waardt, N. Galabretta, G. D. Khoe, and H. J. S. Dorren , “ Wavelength conversion using nonlinear polarization rotation in a single semiconductor optical amplifier ,” IEEE Photonics Technol. Lett. 15 , 90 – 92 ( 2003 ). [CrossRef]  

4 . S. Diez, C. Schmidt, R. Ludwig, H. G. Weber, K. Obermann, S. Kindt, I. Koltchanov, and K. Petermann , “ Four-wave mixing in semiconductor optical amplifiers for frequency conversion and fast optical switching ,” IEEE J. Sel. Top. Quantum Electron. 3 , 1131 – 1145 ( 1997 ). [CrossRef]  

5 . Y. Ueno, S. Nakamura, K. Tajima, and S. Kitamura , “ 3.8-THz wavelength conversion of picosecond pulses using a semiconductor delay-interference signal-wavelength converter (DISC) ,” IEEE Photonics Technol. Lett. 10 , 346 – 348 ( 1998 ). [CrossRef]  

6 . J. Leuthold, B. Mikkelsen, R. E. Behringer, G. Raybon, C. H. Joyner, and P. A. Besse , “ Novel 3R regenerator based on semiconductor optical amplifier delay-interference configuration ,“ IEEE Photonics Technol. Lett. 13 , 860 – 862 ( 2001 ). [CrossRef]  

7 . B. Mikkelsen, K. S. Jepsen, M. Vaa, H. N. Poulsen, K. E. Stubkjaer, R. Hess, M. Duelk, W. Vogt, E. Gamper, E. Gini, P. A. Besse, H. Melchior, S. Bouchoule, and F. Devaux , “ All-optical wavelength converter scheme for high speed RZ signal formats ,” Electron. Lett. 33 , 2137 – 2139 ( 1997 ). [CrossRef]  

8 . S. Nakamura and K. Tajima , “ Bit-rate-transparent non-return-to-zero all-optical wavelength conversion at up to 42 Gb/s by operating Symmetric-Mach-Zehnder switch with new scheme ,” OFC’04 , paper FD3, 2004 .

9 . C. C. Wei, M. F. Huang, and J. Chen , “ Enhancing the frequency response of cross-polarization wavelength conversion ,” IEEE Photonics Tenhnol. Lett. 17 , 1683 – 1685 ( 2005 ). [CrossRef]  

10 . D. Marcenac and A. Mecozzi , “ Switches and frequency converters based on cross-gain modulation in semiconductor optical amplifiers ,” IEEE Photonics Tech. Lett. 9 , 749 – 751 ( 1997 ). [CrossRef]  

11 . H. J. S. Dorren, D. Lenstra, Y. Liu, M. T. Hill, and G. D. Khoe , “ Nonlinear polarization rotation in semiconductor optical amplifiers: theory and application to all-optical flip-flop memories ,” IEEE J. Quantum Electron. 39 , 141 – 148 ( 2003 ). [CrossRef]  

12 . G. Lenz, B. J. Eggleton, C. R. Giles, C. K. Madsen, and R. E. Slusher , “ Dispersive properties of optical filters for WDM systems ,” IEEE J. Quantum Electron. 34 , 1390 – 1402 ( 1998 ). [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. Configuration of a DXPoM wavelength converter
Fig. 2.
Fig. 2. Illustration of a differential sinusoidal wave
Fig. 3.
Fig. 3. Normalized 1 dB angular bandwidth of DXPoM with different values of γ.
Fig. 4.
Fig. 4. Amplitude and delay responses of XPoM, DXPoM with TM delay and DXPoM with TE delay and corresponding simulated eye-diagrams
Fig. 5.
Fig. 5. Amplitude and phase responses with delay of ∆t' opt , 0.95 ∆t' opt and 1.05 ∆t' opt .
Fig. 6.
Fig. 6. Impulse responses of several conditions.

Tables (1)

Tables Icon

Table 1. Parameters Definitions and Their Values

Equations (18)

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

d i k d z = [ Γ k a i k ( η k N N i ) α loss k ] i k , k = TE , TM
d N d t = N τ s i , k [ a i k ( η k N N i ) i k ħ ω i A eff ] + I e V
d σ dt = σ σ 0 τ s i , k 0 , i k Γ k ħ ω i A eff · { exp [ Γ k a i k ( η k σ N i ) ] 1 }
ψ 2 k = α H 2 Γ k a 2 k [ η k σ N 2 ]
Δ ζ = k ( G 1 k 1 ) Γ k a 1 k η k P sat , 1 k Δ P 0,1 k j Ω τ s + 1 + i , k G i k P 0 , i k P sat , i k
Δ P , 2 k = G 2 k P 0,2 k Γ k a 2 k η k Δ ζ
Δ ϕ 2 k = α H 2 Γ k a 2 k η k Δ ζ
pol , 2 = 1 2 cos θ · , 2 TE sin θ · , 2 TM · e j ( Φ + Δ Φ ) 2
Δ P pol , 2 = P θ sin Φ · Δ Φ + 1 2 ( 1 cos Φ ) ( cos 2 θ · Δ P , 2 TE + sin 2 θ · Δ P , 2 TM )
Δ P pol , 2 = [ χ ( 1 + r me ) + ( 1 r me ) ] sin Φ j Ω + τ s 1 + τ stim , 1 1 + τ stim , 2 1 × α H 2 τ s k [ Γ TE a 2 TE η TE P θ Γ k a 1 k η k P sat , 1 k ( G 1 k 1 ) Δ P 0,1 k ]
Δ P pol , 2 = { [ χ ( 1 + r me e j Ω Δ t ) + ( 1 r me e j Ω Δ t ) ] j Ω + τ tot 1 }
× { α H sin Φ 2 τ s k [ Γ TE a 2 TE η TE P θ Γ k a 1 k η k P sat , 1 k ( G 1 k 1 ) Δ P 0,1 k ] }
Γ ( Ω ) 1 γ e j Ω Δ t 1 + j Ω τ tot
T N ( Ω ' ) 1 γ e ' Δt ' ( 1 + j Ω ' ) ( 1 γ )
Ω ' n , min = n 2 1
Ω ' n , max = ( n 1 + γ 1 γ ) 2 1 n 1 + γ 1 γ
Δ t ' opt = 1 γ γ
𝔽 1 { T ( Ω ) 2 π } = e t ' u ( t ' ) γ · e ( t ' Δ t ' ) u ( t ' Δ t ' )
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.