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

Verilog-A behavioral model for resonance-modulated silicon micro-ring modulator

Open Access Open Access

Abstract

We present an accurate behavior model for Si micro-ring modulators (MRM) based on Verilog-A, a standard simulation tool for electronic system design. Our model describes the electrical characteristics of the Si MRM using an equivalent circuit and the optical characteristics based on the couple-mode theory. The accuracy of our model is confirmed by comparing simulation results of our behavior model with the measurement results of a fabricated Si MRM. With this behavior model, co-simulation of Si MRM and electronic driving circuits in the standard electronic design environment can be easily performed.

© 2015 Optical Society of America

1. Introduction

There are growing interests in optical interconnects as they can provide higher bandwidth and more energy-efficient solutions to many interconnect applications than is possible with electrical interconnects [1,2]. In particular, Si photonics is very attractive for optical interconnect applications as it is capable of providing high-performance integrated photonic devices in a cost-effective manner [3,4]. Furthermore, with Si photonics, it is possible to integrate photonic devices and Si electronics on the same Si platform, with which new and more powerful integrated solutions can be realized than is possible with photonics or electronics alone. Several such electronic-photonic integrated circuits (EPICs) have been reported [5–7].

For successful realization of EPICs, it is very important to be able to co-simulate photonic devices and electronic circuits on a single simulation platform in the initial design stage. Such co-simulation can be challenging as simulation techniques for photonics and electronics are usually based on completely different platforms and, consequently, not easily compatible with each other. For this reason, there are growing research interests in developing simulation techniques that can handle both photonic devices and electronic circuits [8–11].

We present a simple and easy-to-use co-simulation method that can be used for EPIC design. In particular, we take an approach in which behavior models for photonic devices are implemented in Verilog-A. Verilog is a popular hardware description language used for high-level behavioral description of electronic systems, and Verilog-A is an extended version of Verilog that supports analog signals. The characteristics of any photonic component can be modeled in the behavior level as long as the values of model parameters can be accurately determined. Although there are many different methods for implementing such behavior models, our approach based on Verilog-A is particularly useful for EPIC applications since Verilog-A behavior models can be easily co-simulated with any electronic circuits in the standard electronic design environment. Establishing such compatibility is very important in order to successfully introduce Si photonic devices into the Si electronics world.

The photonic device of our interest is depletion-type Si micro-ring modulators (MRMs), which are expected to play an important role for high-speed small-footprint optical interconnect systems [12]. Si MRMs can provide very large modulation bandwidth with much smaller footprint and are capable of wavelength division multiplexing (WDM) [13].

This paper is organized as follows. In chapter 2, we give a brief introduction of the Si MRM modulator and describe its electrical and optical characterizations, from which our Verilog-A behavior model is developed. In chapter 3, we verify the accuracy of our behavioral model by comparing model simulation results with measurement results of a fabricated Si MRM. Finally, we conclude this paper in chapter 4. A very brief version of this paper has been reported in [14].

2. Modeling a Si micro-ring modulator

Figure 1 shows the structure of a Si MRM. In this work, we focus on a resonance-modulated depletion-type Si MRM. It consists of ring and bus waveguides, which are located close to each other forming a directional coupler. About three quarters of the ring waveguide are doped to form a lateral PN junction so that its effective group index can be modulated by the reverse-bias voltage applied to the PN junction.

 figure: Fig. 1

Fig. 1 Structure of Si MRM.

Download Full Size | PDF

Si MRMs have both electrical and optical characteristics, which should be modeled separately. For the electrical characteristics, an equivalent circuit model based on passive electrical components can be used [3]. Figure 2 shows the equivalent circuit model where Cox, Cj, Rs, and Rsi represent capacitance through the buried oxide layer, of the depletion region in PN junction, diode series resistance, and resistance through the Si substrate, respectively. CPAD represents the parasitic capacitance of the signal pad used for measurement.

 figure: Fig. 2

Fig. 2 Equivalent circuit model for Si MRM.

Download Full Size | PDF

For optical characteristics of resonance-modulated Si MRM, a model based on the coupled-mode theory can be used, which describes the rate of change in time for the energy amplitude, a(t), inside the ring in the following manner [15,16],

ddta(t)=(jωr1τ)a(t)jμEi(t).
Here, Ei(t) is the incident optical field in the bus waveguide, ωr is the ring resonance angular frequency given as ωr=2πmcηeffL, where m is an integer representing the mode number, ηeff is the effective index for the ring waveguide, and L is the circumference of the ring. τ is decay time constant given as 1τ=1τl+1τe, where τl and τe are the decay time constants due to the round-trip loss and ring-bus coupling, respectively. µ is the mutual coupling coefficient between the ring and the bus and related to τe by μ2=2τe. The output optical field, Et(t), is given as

Et(t)=Ei(t)jμa(t).

When the MRM is electrically modulated by v(t), the values of such parameters as ηeff, τl, and τe change with the applied voltage, and a(t) as well as Et(t) change dynamically according to Eq. (1) and Eq. (2). In order to model such dynamics, it is necessary to solve the differential equation given in Eq. (1) numerically in an efficient manner. The solution for Eq. (1) with the incident optical field Ei(t)=E0exp(jωt) is the sum of the homogeneous and the steady-state solutions, or

a(t)=Pexp(jωrt1τt)+Qexp(jωt),
where P and Q are, respectively, the homogeneous and steady-state coefficients that need to be determined. For numerical solutions, v(t) can be approximated with the series of constant values shown in Fig. 3 assuming the time step is sufficiently small, or
v(t){v0=v(t0),t0t<t1v1=v(t1),t1t<t2vk=v(tk),tkt<tk+1.
Then, corresponding a(t) can be represented as
a(t)={a0(t)=P0exp(jωr,0t1τ0t)+Q0exp(jωt),t0t<t1a1(t)=P1exp(jωr,1t1τ1t)+Q1exp(jωt),t1t<t2ak(t)=Pkexp(jωr,kt1τkt)+Qkexp(jωt),tkt<tk+1,
where Pk and Qk are the homogeneous and steady-state coefficients, and ωr,k and τk are the optical parameter, all of which are influenced by vk when tkt<tk+1.

 figure: Fig. 3

Fig. 3 Step change in the modulation voltage and numerical calculation of a(t).

Download Full Size | PDF

Since the energy in the ring should be continuous at each boundary of time steps, ak-1(tk) should be equal to ak(tk). Assuming a(t0) has reached the steady-state, we can set P0 = 0 and a0(t)=Q0exp(jωt). Q0 can be determined by substituting a0(t) into Eq. (1) as

Q0=j2τe,0j(ωωr,0)+1τ0E0.

a1(t)can be determined if P1 and Q1 are known. Since Q1 is the steady-state coefficient, it should have the same expression as in Eq. (6) with corresponding values for τe, τ, and ωr. P1 can be evaluated at t=t1 from Eq. (3) as

P1=[a0(t1)Q1exp(jωt1)]exp(jωr,1t1+1τ1t1),
where a1(t1)=a0(t1)is used. Above process can be repeated for higher values of index k, and we have following iterative equations:
Pk=[ak1(tk)Qkexp(jωtk)]exp(jωr,ktk+1τktk).
Qk=j2τe,kj(ωωr,k)+1τkE0.
ak(t)=Pkexp(jωr,kt1τkt)+Qkexp(jωt)=[ak1(tk)Qkexp(jωtk)]exp[(jωr,k1τk)(ttk)]+Qkexp(jωt).
When the time step, Δt, is very small, ak(t) for tk1t<tk can be approximated by a constant value, Ak, or

Ak=ak(tk)=ak1(tk)=[ak2(tk1)Qk1exp(jωtk1)]exp[(jωr,k11τk1)Δt]+Qk1exp(jωtk)=[Ak1Qk1exp(jωtk1)]exp[(jωr,k11τk1)Δt]+Qk1exp(jωtk).

These iterative equations can be implemented in Verilog-A as shown in Fig. 4. The input parameters for the model are the voltage applied to the PN junction (vin), wavelength of the incident optical field (vlambda), and clock signal that provides time-step for the calculation (vclk). At every rising edge of the clock signal, the transient response given as Eq. (11) is calculated with the corresponding optical parameters, n_m, tl and te, where each represents the ratio between ηeff and m, τl and τe, respectively. The parameter r is the radius of the ring. The period of the clock signal is the same as the time-step for the calculation, dt. The dependence of such optical parameters as n_m, tl and te on the input voltage should be experimentally determined. In our model, the dependence is modeled with a linear equation for n_m and a quadratic equation for tl and te, as will be described in in Section III.

 figure: Fig. 4

Fig. 4 Behavioral model in Verilog-A.

Download Full Size | PDF

Since Verilog-A does not support complex numbers, variables for the real and the imaginary parts are separately defined and processed. Ra and Ia, Ra_temp and Ia_temp, and Rq and Iq are the real and imaginary parts for ak, ak-1, and Qk, respectively. Exponential functions having imaginary numbers are converted into the sinusoidal functions.

3. Model verification

The Si MRM device used for verifying our behavior model is fabricated through OpSIS-IME multi-project wafer foundry service. Figures 5(a) and 5(b) show the microphotograph and the cross-section of the fabricated device. The radius of the ring waveguide is 8-μm and the nominal gap between ring and bus waveguides is 300-nm. It is implemented with 220-nm thick Si layer on 2-μm thick buried oxide layer. 500-nm wide rib waveguide is used for ring and bus waveguides. A lateral PN junction is formed in the ring waveguide except around the directional coupler. The nominal doping concentrations for the junction are 2 x 1018 /cm3 for P region and 3 x 1018 /cm3 for N region.

 figure: Fig. 5

Fig. 5 (a) Microphotograph and (b) cross-section of fabricated Si MRM.

Download Full Size | PDF

As the first step, numerical values of parameters used in our model have to be extracted from the fabricated Si MRM. Electrical model parameters are extracted by fitting simulated S11 coefficients of the equivalent circuit shown in Fig. 2 to measured S11. Figures 6(a) and 6(b) show the magnitude and phase responses of measured (circles) and simulated by ADS (solid lines) S11. The extracted numerical values of the electrical parameters are given in Table 1.

 figure: Fig. 6

Fig. 6 Curve-fitting of the simulated S11 to the measurements: (a) Magnitude response (b) Phase response.

Download Full Size | PDF

Tables Icon

Table 1. Extracted parameter values of the Si MRM

Optical model parameter values are extracted by fitting the measured optical power transmission characteristics of the Si MRM to the steady-state responses of our model obtained from Eq. (1) and Eq. (2) by letting ddta(t)=jωa(t), or

|EtEi|2=|j(ωωr)+1τl1τej(ωωr)+1τl+1τe|2.
Specifically, ηeffm is determined from the measured resonance wavelength using the resonance condition mλres=ηeffL, where λres represents the resonance wavelength. The values for τl and τe are determined with the minimum mean square error curve-fitting technique [17]. Parameter extraction is done for three different reverse bias voltages of 0, 1, 2V. The extracted values at each bias voltage are 0.0308674, 0.0308679, 0.0308682 for ηeffm, 18.7081 ps, 19.2456 ps, 19.5853 ps for τl, and 21.8929 ps, 21.8932 ps, 21.8934 ps for τe. The voltage dependence of ηeffm, τl and τe are approximated by a linear equation for ηeffm, and quadratic equations for τl and τe as shown in Figs. 7(b) and 7(c).

 figure: Fig. 7

Fig. 7 (a) Measured and fitted optical transmission characteristics. (b) Extracted ηeff / m. (c) Extracted τl and τe.

Download Full Size | PDF

Using our model with the extracted parameters, the Si MRM dynamic characteristics can be easily simulated with Spectre in Cadence, one of the most widely-used CAD tools for electronic circuit design. Figure 8(a) shows the schematic used for such simulation. As mentioned before, a clock signal is used for in providing the time-step for model simulation.

 figure: Fig. 8

Fig. 8 (a) Simulation setup. (b) Measurement setup.

Download Full Size | PDF

In order to confirm the accuracy of our model, the simulation results for the Si MRM transient responses are compared with the measurement results. The measurement setup is shown in Fig. 8(b). The tunable laser source provides the light into the Si MRM through a grating coupler. The PPG provides the electrical signals to the modulator with the desired DC offset voltage. In order to avoid any self-heating of the Si MRM [18], the optical power injected into the Si MRM is minimized during measurement. With −5-dBm injection into the grating coupler and about 5-dB coupling loss, the optical power level inside the waveguide is about −10-dBm. This requires amplification of output optical signal with an erbium-doped fiber amplifier. The optical receiver converts the optical output into the electrical signals, which are observed with an oscilloscope.

For accurate comparison between simulation and measurement results, the measurement condition should be accurately accounted for. For this, rising and falling times of electric signals provided by the PPG as well as the frequency response of the optical receiver used in measurement are measured and their equivalent circuits are used in simulation. The influence of the optical amplifier is not considered in our simulation since it does not influence the dynamics of our measurement.

Figures 9(a), 9(b) and 9(c) show the measured and simulated step responses of the Si MRM at three different input wavelengths of 1551.45 nm, 1551.50 nm, and 1551.55 nm, respectively. The resonance wavelength of Si MRM at 0V bias is 1551.565 nm as can be seen in Fig. 7(a). For this measurement, a voltage step from 0V to 2V is applied to Si MRM by the PPG. A time step of 200-fs is used for simulation. Different amounts of peaking are observed for Si MRMs at different wavelengths, which is due to the wavelength-dependent frequency response of the Si MRM [19–21]. Figure 10 shows the normalized modulation frequency responses measured at three different wavelengths as well as the simulated results with our model. For this measurement, an RF source is used instead of PPG as an input, and the output signals are observed with a spectrum analyzer. When the separation between input and resonance wavelengths is larger, the MRM has a peaking at higher frequency, which results in the larger overshoot for the step response. As can be seen in the figures, our behavioral model accurately describes such wavelength-dependent step and frequency responses.

 figure: Fig. 9

Fig. 9 Measured and simulated step response at the wavelength of (a) 1551.45 nm (b) 1551.50 nm (c) 1551.55 nm.

Download Full Size | PDF

 figure: Fig. 10

Fig. 10 Measured and simulated frequency response at the wavelength of 1551.45 nm, 1551.50 nm, and 1551.55 nm.

Download Full Size | PDF

For further verification, simulated and measured eye diagrams are compared for 28-Gb/s 231-1 PRBS data. For obtaining simulated eye diagrams, first 40 bits before the Si MRM reaches the steady-state are discarded. Similar to the step responses, we observe the wavelength dependent eye diagrams as shown in Figs. 11(a), 11(b) and 11(c). The maximum on-off ratio is observed for at λ = 1551.55 nm, but the maximum eye opening is observed at λ = 1551.5 nm, where the SI MRM has larger modulation bandwidth. With our model, which can accurately describe the wavelength-dependent Si MRM dynamics, we can take the Si MRM as one of the easily accessible function blocks in the design environment and use it along with other electronic blocks for realizing optimally-performing optical transmitters for various optical interconnect applications.

 figure: Fig. 11

Fig. 11 Measured and simulated eye-diagrams for 28-Gb/s PRBS 231-1 pattern at the wavelength of (a) 1551.45 nm (b) 1551.50 nm (c) 1551.55 nm.

Download Full Size | PDF

4. Conclusion

We developed an accurate behavioral model for the Si MRM in Verilog-A. The accuracy of our model was experimentally verified with a fabricated Si MRM device. In addition, we provided the parameter extraction technique, whose accuracy is extremely important for any device model, and details of the numerical calculation process. Our approach of implementing the Si MRM model in Verilog-A is highly valuable as it allows accurate simulation of the Si MRM in the well-established electronic design environment. We believe such compatibility is very important for successful realization of Si photonic devices integrated with Si electronic circuits.

Acknowledgments

This work was supported by the National Research Foundation of Korea funded by the Korean Government (MEST) under Grant 2012R1A2A1A01009233.

References and links

1. I. A. Young, E. Mohammed, J. T. S. Liao, A. M. Kern, S. Palermo, B. A. Block, M. R. Reshotko, and P. L. D. Chang, “Optical I/O technology for tera-scale computing,” IEEE J. Solid-State Circuits 45(1), 235–248 (2010). [CrossRef]  

2. M. A. Taubenblatt, “Optical interconnects for high-performance computing,” IEEE/OSA J. Lightwave Technol. 30(4), 448–457 (2012). [CrossRef]  

3. X. Zheng, D. Patil, J. Lexau, F. Liu, G. Li, H. Thacker, Y. Luo, I. Shubin, J. Li, J. Yao, P. Dong, D. Feng, M. Asghari, T. Pinguet, A. Mekis, P. Amberg, M. Dayringer, J. Gainsley, H. F. Moghadam, E. Alon, K. Raj, R. Ho, J. E. Cunningham, and A. V. Krishnamoorthy, “Ultra-efficient 10 Gb/s hybrid integrated silicon photonic transmitter and receiver,” Opt. Express 19(6), 5172–5186 (2011). [CrossRef]   [PubMed]  

4. G. Kim, J. W. Park, I. G. Kim, S. Kim, S. Kim, J. M. Lee, G. S. Park, J. Joo, K.-S. Jang, J. H. Oh, S. A. Kim, J. H. Kim, J. Y. Lee, J. M. Park, D.-W. Kim, D.-K. Jeong, M.-S. Hwang, J.-K. Kim, K.-S. Park, H.-K. Chi, H.-C. Kim, D.-W. Kim, and M. H. Cho, “Low-voltage high-performance silicon photonic devices and photonic integrated circuits operating up to 30 Gb/s,” Opt. Express 19(27), 26936–26947 (2011). [CrossRef]   [PubMed]  

5. J. F. Buckwalter, X. Zheng, G. Li, K. Raj, and A. V. Krishnamoorthy, “A monolithic 25-Gb/s transceiver with photonic ring modulators and Ge detectors in a 130-nm CMOS SOI process,” IEEE J. Solid-State Circuits 47(6), 1309–1322 (2012). [CrossRef]  

6. L. Zimmermann, D. J. Thomson, B. Goll, D. Knoll, S. Lischke, F. Y. Gardes, Y. Hu, G. T. Reed, H. Zimmermann, and H. Porte, “Monolithically integrated 10Gbit/sec Silicon modulator with driver in 0.25μm SiGe:C BiCMOS,” in Proceedings of European Conference on Optical Communication (ECOC), pp. 1–3 (2013).

7. J.-S. Youn, M.-J. Lee, K.-Y. Park, H. Rücker, and W.-Y. Choi, “An integrated 12.5-Gb/s optoelectronic receiver with a silicon avalanche photodetector in standard SiGe BiCMOS technology,” Opt. Express 20(27), 28153–28162 (2012). [CrossRef]   [PubMed]  

8. T. Smy, P. Gunupudi, S. Mcgarry, and W. N. Ye, “Circuit-level transient simulation of configurable ring resonators using physical models,” J. Opt. Soc. Am. B 28(6), 1534–1543 (2011). [CrossRef]  

9. W. Bogaerts, M. Fiers, and P. Dumon, “Design challenges in silicon photonics,” IEEE J. Sel. Top. Quantum Electron. 20(4), 8202008 (2014). [CrossRef]  

10. L. Chrostowski, J. Flueckiger, C. Lin, M. Hochberg, J. Pond, J. Klein, J. Ferguson, and C. Cone, “Design methodologies for silicon photonic integrated circuits,” Proc. SPIE 8989, 89890G (2014).

11. K. Zhu, V. Saxena, and W. Kuang, “Compact Verilog-A modeling of silicon traveling-wave modulator for hybrid CMOS photonic circuit design,” in Proceedings of IEEE International Midwest Symposium on Circuits and Systems, pp. 615–618 (2014). [CrossRef]  

12. D. J. Thomson, F. Y. Gardes, J.-M. Fedeli, S. Zlatanovic, Y. Hu, B. P. P. Kuo, E. Myslivets, N. Alic, S. Radic, G. Z. Mashanovich, and G. T. Reed, “50-Gb/s silicon optical modulator,” IEEE Photon. Technol. Lett. 24(4), 234–236 (2012). [CrossRef]  

13. X. Zheng, F. Y. Liu, J. Lexau, D. Patil, G. Li, Y. Luo, H. D. Thacker, I. Shubin, J. Yao, K. Raj, R. Ho, J. E. Cunningham, and A. V. Krishnamoorthy, “Ultralow power 80 Gb/s arrayed CMOS silicon photonic transceivers for WDM optical links,” IEEE/OSA J. Lightwave Technol. 30(4), 641–650 (2012). [CrossRef]  

14. J. Rhim, Y. Ban, J.-M. Lee, S.-H. Cho, and W.-Y. Choi, “A behavior model for silicon micro-ring modulators and transmitter circuit-level simulation using it,” in Proceedings of IEEE International Conference on Group IV Photonics, pp. 7–8 (2014). [CrossRef]  

15. B. E. Little, S. T. Chu, H. A. Haus, J. Foresi, and J.-P. Laine, “Microring resonator channel dropping filters,” J. Lightwave Technol. 15(6), 998–1005 (1997). [CrossRef]  

16. L. Zhang, Y. Li, J.-Y. Yang, M. Song, R. G. Beausoleil, and A. E. Willner, “Silicon-based microring resonator modulators for intensity modulation,” IEEE J. Sel. Top. Quantum Electron. 16(1), 149–158 (2010). [CrossRef]  

17. B.-M. Yu, J.-M. Lee, Y. Ban, S.-H. Cho, and W.-Y. Choi, “Model parameter extraction for Si micro-ring modulators,” presented at OptoElectronics and Communication Conference and the Australian Conference on Optical Fibre Technology (OECC/ACOFT), Melbourne, Australia, 6–10 July, 2014.

18. X. Zheng, Y. Luo, G. Li, I. Shubin, H. Thacker, J. Yao, K. Raj, J. E. Cunningham, and A. V. Krishnamoorthy, “Enhanced optical bistability from self-heating due to free carrier absorption in substrate removed silicon ring modulators,” Opt. Express 20(10), 11478–11486 (2012). [CrossRef]   [PubMed]  

19. W. D. Sacher, W. M. J. Green, S. Assefa, T. Barwicz, H. Pan, S. M. Shank, Y. A. Vlasov, and J. K. S. Poon, “Coupling modulation of microrings at rates beyond the linewidth limit,” Opt. Express 21(8), 9722–9733 (2013). [CrossRef]   [PubMed]  

20. J. Müller, F. Merget, S. S. Azadeh, J. Hauck, S. R. García, B. Shen, and J. Witzens, “Optical peaking enhancement in high-speed ring modulators,” Sci. Rep. 4, 6310 (2014). [CrossRef]   [PubMed]  

21. Y. Ban, J.-M. Lee, B.-M. Yu, S.-H. Cho, and W.-Y. Choi, “Small-signal frequency responses for Si micro-ring modulators,” in Proceedings of IEEE Optical Interconnects Conference, pp. 47–48 (2014). [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 (11)

Fig. 1
Fig. 1 Structure of Si MRM.
Fig. 2
Fig. 2 Equivalent circuit model for Si MRM.
Fig. 3
Fig. 3 Step change in the modulation voltage and numerical calculation of a(t).
Fig. 4
Fig. 4 Behavioral model in Verilog-A.
Fig. 5
Fig. 5 (a) Microphotograph and (b) cross-section of fabricated Si MRM.
Fig. 6
Fig. 6 Curve-fitting of the simulated S11 to the measurements: (a) Magnitude response (b) Phase response.
Fig. 7
Fig. 7 (a) Measured and fitted optical transmission characteristics. (b) Extracted ηeff / m. (c) Extracted τl and τe.
Fig. 8
Fig. 8 (a) Simulation setup. (b) Measurement setup.
Fig. 9
Fig. 9 Measured and simulated step response at the wavelength of (a) 1551.45 nm (b) 1551.50 nm (c) 1551.55 nm.
Fig. 10
Fig. 10 Measured and simulated frequency response at the wavelength of 1551.45 nm, 1551.50 nm, and 1551.55 nm.
Fig. 11
Fig. 11 Measured and simulated eye-diagrams for 28-Gb/s PRBS 231-1 pattern at the wavelength of (a) 1551.45 nm (b) 1551.50 nm (c) 1551.55 nm.

Tables (1)

Tables Icon

Table 1 Extracted parameter values of the Si MRM

Equations (12)

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

d dt a(t)=( j ω r 1 τ )a(t)jμ E i (t).
E t (t)= E i (t)jμa(t).
a(t)=Pexp( j ω r t 1 τ t )+Qexp( jωt ),
v(t){ v 0 =v( t 0 ), t 0 t< t 1 v 1 =v( t 1 ), t 1 t< t 2 v k =v( t k ), t k t< t k+1 .
a(t)={ a 0 (t)= P 0 exp( j ω r,0 t 1 τ 0 t )+ Q 0 exp( jωt ), t 0 t< t 1 a 1 (t)= P 1 exp( j ω r,1 t 1 τ 1 t )+ Q 1 exp( jωt ), t 1 t< t 2 a k (t)= P k exp( j ω r,k t 1 τ k t )+ Q k exp( jωt ), t k t< t k+1 ,
Q 0 = j 2 τ e,0 j( ω ω r,0 )+ 1 τ 0 E 0 .
P 1 =[ a 0 ( t 1 ) Q 1 exp( jω t 1 ) ]exp( j ω r,1 t 1 + 1 τ 1 t 1 ),
P k =[ a k1 ( t k ) Q k exp( jω t k ) ]exp( j ω r,k t k + 1 τ k t k ).
Q k = j 2 τ e,k j( ω ω r,k )+ 1 τ k E 0 .
a k (t)= P k exp( j ω r,k t 1 τ k t )+ Q k exp( jωt ) =[ a k1 ( t k ) Q k exp( jω t k ) ]exp[ ( j ω r,k 1 τ k )( t t k ) ]+ Q k exp( jωt ).
A k = a k ( t k )= a k1 ( t k ) =[ a k2 ( t k1 ) Q k1 exp( jω t k1 ) ]exp[ ( j ω r,k1 1 τ k1 )Δt ]+ Q k1 exp( jω t k ) =[ A k1 Q k1 exp( jω t k1 ) ]exp[ ( j ω r,k1 1 τ k1 )Δt ]+ Q k1 exp( jω t k ).
| E t E i | 2 = | j( ω ω r )+ 1 τ l 1 τ e j( ω ω r )+ 1 τ l + 1 τ e | 2 .
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.