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

Critical conditions of control light for a silicon-based photonic crystal bistable switching

Open Access Open Access

Abstract

In this paper, we investigate the critical conditions of control light for a silicon-based photonic crystal bistable switching. By establishing a time-dependent evolution equation, the critical pump power and pump time of the control light are derived, respectively. It is found that with the increase of the frequency detuning of the incident light, the critical power of the control light will rise, while the corresponding critical time will be shortened. It is also revealed that under the same conditions, the critical total power of the multiple-beam control light is less than the one of the single-beam control light. The theoretical predictions show perfect agreement with the simulation results.

©2010 Optical Society of America

1. Introduction

In the last decade, photonic crystal (PC) bistable switching has attracted many attentions for its great potential in all-optical signal processing [13]. Commonly, to observe the bistable switching action, one should perform up- and downsweeping power cycles by coupling a continuous wave (CW) signal light into the PC waveguide. Alternatively, this bistable action can also realized under pulsed injection of the control light when the power of the CW signal light is fixed at a certain value within the bistable region [35]. Obviously, from a practical point of view, the latter method is more feasible. However, what is the critical condition of the control light (i.e., the minimum values of the pump power and pump time of the control light to trigger a bistable switching to the “on” state)? So far, this question is rarely studied. Of course, one might obtain these critical values by experiments or simulations after many times of attempts. However, it is time consuming and the achieved results cannot be generalized. Therefore, a theoretical investigation on the critical pump conditions becomes very important in PC switching design work.

To date, almost all of the proposed control-light sources are high-peak-power Gaussian pulses [36]. Unfortunately, it is very difficult to obtain an analytic solution by this way. In this paper, we present a new trigger method: using a continuous wave (CW) source with suitable launching time to act as control light. We will show this new method is more efficient for the theoretical study on the critical pump conditions, and the achieved results can be readily generalized to the case of pulsed control light. To verify the theoretical predictions, nonlinear finite-difference time-domain (FDTD) method is employed to make simulations. The paper is organized as follows. In section 2, we establish a time-dependent evolution equation to describe the PC bistable switching under the pump of a CW signal light superposed with a CW control light. Based on this model, the critical conditions of the pump power and pump time of the control light are derived in sections 3 and 4, respectively, and FDTD simulations are performed to make verifications. Subsequently, in section 5, these critical conditions for single-beam-CW control light are generalized to the cases of multiple-beam-CW control light and pulsed control light. Finally, a brief summary is given in section 6.

2. Time-dependent evolution equation for a PC bistable switching: using a single-beam-CW control light

As proposed in Refs [7]. and [8], two-photon absorption (TPA) and Kerr effect are two dominant nonlinear factors in a Si system. In our previous work [9], we have established a simple model to study the influence of TPA on the characteristics of bistable switching, and the threshold power (up- and down-jump) of the signal light can be precisely calculated based on this model. Now, we will fix the power of the signal light, and further investigate the critical conditions of the control light.

To describe conveniently, we depict a typical structure of PC bistable switching in Fig. 1 , which is consist of one nonlinear microcavity with two waveguides (WG) connected to the two ends. Suppose the decay rates of the input WG (WG1) and output WG (WG2) are γ 1 and γ 2, respectively, and the intrinsic decay rate of the nonlinear cavity is γ C. Phenomenologically, γ C can be written as γ C = γ 0 + γ TPA, where γ 0 is the linear intrinsic decay rate, and γ TPA is the decay rate correlated with TPA. Therefore, the total decay rate is γ = γ 1 + γ 2 + γ 0 + γ TPA.

 figure: Fig. 1

Fig. 1 Sketch of a typical bistable switching composed of a resonator cavity coupled to both input and output WGs.

Download Full Size | PDF

For simplicity, we suppose the CW control light has the same frequency ω and initial phase (set to be zero in our case) with the CW signal light. When the signal light superposed with the control light is launched into the nonlinear cavity along WG1, according to the coupled mode theory (CMT) [10], the mode amplitude A can be written as

{dAdt=[j(ω0γ|sT|2p0)γ]A+2γ1(ss+sx)γ=γ0+γ1+γ2+γTPAsT=2γ2Ass+sx+sR=2γ1A,
where s s, s x, s T and s R are the field amplitudes of the incident signal light, control light, transmitted light and reflective light, respectively, ω 0 is the resonant frequency of the PC cavity, and p 0 is the characteristic power [4]. Noticing that γ TPA is proportional to the intensity of energy stored in the PC cavity, γ TPA can be expressed as [9]
γTPA=kγ1|sT|2,
where k is defined as TPA coefficient which denotes the strength of the TPA process, and can be calculated by theory or simulation [9].

Obviously, it is very difficult to solve Eq. (1) directly. However, if the power of the incident light is not very high, we can utilize a perturbation technology and obtain an approximate solution (suppose there is not any photon energy stored in the PC cavity at the initial moment)

A=2γ1(|ss|+|sx|)i(ωω0+γ|sT|2/p0)+γ[eiωteγteiω0t]
and
ST=4γ1γ2(|ss|+|sx|)i(ωω0+γ|sT|2/p0)+γ(eiωteγteiω0t),
where ω is the frequency of the incident light. Thus, the instant values of the energy stored inside the cavity and the transmission rate can be easily calculated by using Eqs. (3) and (4)
|A|2(t)=(ps+px)2/2γ2(1+kη|sT|2/2)2η(δ|sT|2/p0)2+1[e2γt2eγtcos(ω0ω)t+1],
T(t)=|sT|2|ss|2=(ps+px)2/ps(1+kη|sT|2/2)2η(δ|sT|2/p0)2+1[e2γt2eγtcos(ω0ω)t+1],
whereps=|ss|2andpx=|sx|2are the incident power of the signal light and control light, respectively; η=4γ1γ2/(γ1+γ2+γ0)2 is the linear peak transmission rate of the bistable switching; δ=(ω0ω)/(γ1+γ2+γ0) is the frequency detuning of the input CW with respect to the cavity mode.

Equations (5) and (6) are the time-dependent evolution equations for PC bistable switching, under the pump of a CW signal light superposed with a CW control light. From these equations, one can see that both the photon energy inside the cavity and the transmission rate will decay with time, and access a steady-state value ultimately.

If the control light is shut off at t 0, then the instant photon energy stored in the cavity |A|2(t0) can be calculated by Eq. (5). Simultaneously, |A|2(t0) can be treated as the initial condition for the following evolution of the energy in cavity, and its value will determine whether the final state of the bistable switching is “on” or “off” state.

3. Critical pump-power condition of the control light

3.1 Theory

Now, let us discuss the critical pump-power condition of the control light. To describe clearly, we depict a working principle diagram of a bistable switching in Fig. 2 , where ABCDA is the bistable region and the horizontal axis indicates the power of the signal light. The curve of AB denotes the “off” state, while DC denotes the “on” state. It is well known to all, when there is only signal light being incident, the “on” and “off” states can be achieved by upward and downward sweeps of the signal power level. We define A as the down-jump point (D→A), and the corresponding incident power of the signal light is p 1; B is the up-jump point (B→C), and the corresponding incident power of the signal light is p 2. Alternatively, we can also fix the power of the signal light at a certain value within the bistable region, and trigger the switching to the “on” state by a control light. Here we should keep in mind that the jump from the “off” state to “on” state is indirect, e.g., the jump from A to D must pass through B and C points, although this process may be very quick.

 figure: Fig. 2

Fig. 2 Working principle diagram of a bistable switching

Download Full Size | PDF

Noticing that the control light has been assumed to have the same frequency and initial phase with the signal light, we can treat them as a whole, i.e., an “equivalent signal light” with incident power of(p1+px)2. Therefore, to make the switching jump indirectly from A to D, the power of the “equivalent signal light” should be no less than p 2. Thus we have

(p1+px)2p2,
i.e.,

px(p2p1)2.

Equation (8) is the critical pump-power condition of the control light to trigger a bistable switching to “on” state when the signal power is p 1. That is to say, if p x is less than the critical value, one will not be able to trigger the bistable switching successfully even if the pump time of the control light is long enough.

3.2 FDTD simulation: verify the theoretical predictions

To verify the above mentioned theoretical predictions, we choose to study a typical structure of 2D PC bistable switching. The structure and the corresponding parameters are the same as the ones presented in our previous work [11], except that in this case the central defect cavity is nonlinear, and the real part of the nonlinear coefficient is 1 x 10−5 μm2/W, while the imaginary part is set to be 2.5 x 10−6 μm2/W, denoting TPA effect. To measure the transmission rate, one monitor is set at the output port of WG2. The grid sizes in the horizontal and vertical directions are chosen to be a/20 (where a = 1μm is the lattice constant), and a perfectly matched layer (PML) of 1μm is employed as boundary. Thus, we can begin numerical simulation and focus on the TM modes. By employing nonlinear FDTD technique [12], it is not difficult to obtain the values of p 1 and p 2 (these values can also be obtained by theoretical calculations, as presented in Ref [9].).

As an example, when the frequency detuning of the incident signal light is δ = 4.3478, p 1 and p 2 are calculated to be 0.120 kW/μm and 0.215 kW/μm, respectively. Thus, the critical pump power of the control light at the down-jump point (p s = p 1) can be obtained immediately by using Eq. (8): p x = 0.0136 kW/μm. To examine the precision of the critical value predicted by theory, we begin simulation by using p x = 0.0136 kW/μm and p x = 0.0134 kW/μm as the power of the CW control light, respectively, for comparison. When a CW signal light (δ = 4.3478, p s = p 1 = 0.120 kW/μm) superposed with the control light (which has the same frequency and initial phase with the signal light) is launched into the PC switching along WG1, the time-dependent evolution process of the transmission rate recorded by a monitor is shown in Fig. 3 (since the critical pump time keeps unknown at present,to ensure a long enough pump time in this simulation, the pump time of the control light is selected to be 4000(a/c), about 13.33ps). One can see when p x = 0.0136 kW/μm, the bistable switching is successfully triggered to the “on” state with T = 0.6; while when p x is slightly weakened to 0.0134 kW/μm, the switching can only works on the “off” state with T = 0.054.

 figure: Fig. 3

Fig. 3 Evolvement curves of bistable switching for different pump power of the control light. The upper curve: p x = 0.0136 kW/μm; The lower curve: p x = 0.0132 kW/μm.

Download Full Size | PDF

In Table 1 , we make a further comparison between the theoretical predictions and FDTD simulation results, with respect to different frequency detuning, where pxTheory is calculated by Eq. (8), and pxFDTDis obtained by FDTD simulation. One can see clearly that in a great range of frequency detuning, the theoretical predictions agree with the FDTD simulation results very well.

Tables Icon

Table 1. Critical pump power of the control light with respect to different frequency detuning.

4. Critical pump-time condition of the control light

In section 3, we have obtained the critical pump-power condition of the control light. However, this sole condition cannot ensure the bistable switching to be successfully triggered to an “on” state, because the decisive factor is the photon energy stored inside the cavity, which is surely related with the pump time. Now, let us make a discussion on the critical pump-time condition of the control light.

To make the switching jump indirectly from A to D, the photon energy inside the cavity before the control light is shut off should be accumulated to exceed the energy value at D point, say, |A|D2=p1T1/(2γ2), where T 1 is the steady-state transmission rate at D point after the control light is shut off. Thus, we obtain the critical condition of the energy stored inside the cavity for the incident signal power p 1

|A|2p1T1/(2γ2).
Obviously, Eqs. (8) and (9) can be readily generalized to the arbitrary signal power in the bistable region by substituting p 1 with p s.

Since the dynamic evolution process of the energy in cavity is determined by Eq. (5), thus, the critical pump-time condition of the control light can be directly derived by using Eqs. (5) and (9):

(ps+px)2(1+kη|sT|2/2)2η(δ|sT|2/p0)2+1[e2γt2eγtcos(ω0ω)t+1]Tsps,
where ps is the incident signal power within the bistable region and T s is the corresponding steady transmission rate of “on” state.

Specially, let us make a discussion on the critical pump time of the control light with a critical power at the down-jump point (p s = p 1, T s = T 1). In this case, we have

(p1+px)2=p2.
Substituting formula (11) into Eq. (10), we obtain
p2/p1(1+kη|sT|2/2)2η(δ|sT|2/p0)2+1[e2γt2eγtcos(ω0ω)t+1]T1.
From Eq. (12) the desired critical pump time of the control light can be calculated. Noticing that the left part of Eq. (12) is actually the time-dependent transmission rate under the pump of the signal light and control light, Eq. (12) can be rewritten as a more brief form

T(t)T1.

Obviously, it is rather difficult to obtain an analytic solution of the critical pump time of the control light directly from Eq. (12), since both|sT|2and γ are time dependent. Therefore, we have to fall back on a numerical method: plot the evolution curve of T(t) and a horizon line of T = T 1, and the first intersection point of them reads the critical value, as shown in Fig. 4(a) (δ = 4.3478). One can see clearly that the critical pump time of the control light is t s = 3522(a/c) [before we know this critical value, the pump time of the control light should be long enough, e.g., 4000(a/c) in our case]. To verify it, in Fig. 4(b) we make a FDTD simulation on the bistable switching when the control light is shut off at the critical value of t s = 3522(a/c) (about 11.74ps). One can see that the bistable switching is successfully triggered to the “on” state with transmission rate of 0.6.

 figure: Fig. 4

Fig. 4 Critical pump time of the control light for δ = 4.3478. (a) Illustration for how to obtain the critical pump time of the control light. (b) Evolvement curve of transmission when the control light is shut off at t = 3522(a/c).

Download Full Size | PDF

To further examine the precision of the critical pump time calculated by Eq. (12), we make a comparison between the theoretically-calculated value tsTheoryand FDTD simulation resulttsFDTDin a very wide range of the frequency detuning, as shown in Table 2 , and a perfect agreement is achieved between them.

Tables Icon

Table 2. Critical pump time of the control light with respect to different frequency detuning.

Form Tables 1 and 2, it is also found that with the increase of the frequency detuning, the critical pump power of the control light will rise while the corresponding critical pump time is shortened.

5. Critical condition of multiple-beam-CW control light and pulsed control light

So far, we have discussed the critical conditions for single-beam-CW control light. Before the further investigation on the case of pulsed control light, we would like to continue our discussion on the multiple-beam-CW-control-light case, since it is the basis for the former, as will be seen below.

When a bistable switching is simultaneously triggered by N beams of CW control light, for convenience, we suppose these beams of control light have the same frequency and initial phase with the signal light. Therefore, the total power of the incident light is

pin=(ps+i=1Npxi)2,
where pxidenotes the power of the ith beam of control light.

The following discussion is exactly the same with the one performed for single-beam-CW control light provided that px is substituted byi=1Npxi. Thus, Eq. (8) can be directly transformed as

(i=1Npxi)2(p2p1)2.
Equation (15) is the critical pump-power condition of the multiple-beam-CW control light, and the corresponding critical pump power is determined by
f(t)exp(iω0t)
Since(i=1Npxi)2i=1Npxi, from Eq. (16) we obtain
i=1Npxi(p2p1)2.
Comparing Eq. (17) with (11), it is not difficult to deduce that under the same conditions, the critical total power of the multiple-beam-CW control light is less than the one of the single-beam-CW control light, especially for a greater N. Enlightened by this find, if the power of the control-light source in lab is unfortunate to be lower than the critical value, we still have opportunity to trigger the bistable switching by dividing the control light into several beams. Obviously, this method is promising in actual applications.

Additionally, the critical pump time of multiple-beam-CW control light can also be calculated by Eq. (12). As for the other incident signal power in bistable region, the corresponding critical conditions can also be readily obtained only by substituting p 1with p s in Eqs. (12) and (15).

Finally, let us make a brief discussion on the case of pulsed control light.

We suppose the pulse can be expressed as “f(t)exp(iω0t)”, where f(t) is the envelopment function in time domain, and ω 0 is the central frequency. By using Fourier transformation, the pulse can be approximately regarded as a superposition of a series of CWs with frequencies near ω 0 and appropriate weights:

f(t)exp(iω0t)=ig(ωi)exp(iωit)Δω,
where ω i and g(ω i) are the frequency and field amplitude of the ith CW component, respectively, Δω is the step width in frequency domain. Thus, we can use “|g(ω i)|2Δω” as the corresponding weight factor, and a smaller Δω will lead to a more precise result.

Specially, if the control pulse is relatively wide (e.g., tens of ps), the frequency of each CW component is approximately equal to ω 0. In this case, the problem can be converted to the one in the multiple-beam-CW-control-light case, which has been above discussed. While for an ultrashort control pulse (e.g, several fs), the method presented above is not a good approximation anymore (due to the further broadening of the frequency spectrum), and a more precise model including multiple-beam-CW control light with deferent frequencies is required. However, we can still predict that the critical pump power of CWs with deferent frequencies must be greater than the one with same frequency, as has been testified by our FDTD simulations.

6. Conclusion

In summary, by establishing a time-dependent evolution equation for PC bistable switching, we have derived the critical conditions of the pump power and pump time of the control light, respectively, and the theoretical predictions are precisely verified by FDTD simulations. It is found that with the increase of the frequency detuning of the incident light, the critical power of the control light will rise, while the corresponding critical time will decrease. It is also found that under the same conditions, the critical total power of the multiple-beam-CW control light is less than the one of the single-beam-CW control light. The theory and technology presented in this paper will be helpful for actual designs of PC switching.

Acknowledgments

This paper is supported by Key Technologies R&D Program of Guangdong Province, Major Research Plan (No. 2009A080301013), and Key Technologies R&D Program of Guangzhou City, Major Plan Program (No. 2009A1-D081), and National Natural Science Foundation of China (No. 60835001). Chao Li and Hong Wang contributed equally to this work.

References and links

1. K. Ogusu and K. Takayama, “Optical bistability in photonic crystal microrings with nonlinear dielectric materials,” Opt. Express 16(10), 7525–7539 (2008). [CrossRef]   [PubMed]  

2. M. F. Yanik, S. H. Fan, and M. Soljacic, “High-contrast all-optical bistable switching in photonic crystal microcavities,” Appl. Phys. Lett. 83(14), 2739–2741 (2003). [CrossRef]  

3. M. K. Kim, I. K. Hwang, S. H. Kim, H. J. Chang, and Y. H. Lee, “All-optical bistable switching in curved microfiber-coupled photonic crystal resonators,” Appl. Phys. Lett. 90(16), 161118 (2007). [CrossRef]  

4. M. Soljačić, M. Ibanescu, S. G. Johnson, Y. Fink, and J. D. Joannopoulos, “Optimal bistable switching in nonlinear photonic crystals,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 66(5), 055601 (2002). [CrossRef]  

5. M. Belotti, J. F. Galisteo Lòpez, S. De Angelis, M. Galli, I. Maksymov, L. C. Andreani, D. Peyrade, and Y. Chen, “All-optical switching in 2D silicon photonic crystals with low loss waveguides and optical cavities,” Opt. Express 16(15), 11624–11636 (2008). [PubMed]  

6. T. Tanabe, M. Notomi, S. Mitsugi, A. Shinya, and E. Kuramochi, “All-optical switches on a silicon chip realized using photonic crystal nanocavities,” Appl. Phys. Lett. 87(15), 151112 (2005). [CrossRef]  

7. G. Priem, P. Bienstman, G. Morthier, and R. Baets, “Impact of absorption mechanisms on Kerr-nonlinear resonator behavior,” J. Appl. Phys. 99(6), 063103 (2006). [CrossRef]  

8. D. Brissinger, B. Cluzel, A. Coillet, C. Dumas, P. Grelu, and F. de Fornel, “Near-field control of optical bistability in a nanocavity,” Phys. Rev. B 80(3), 033103 (2009). [CrossRef]  

9. C. Li, J. F. Wu, and W. C. Xu, “Influence of two-photon absorption on bistable switching in a silicon photonic crystal microcavity,” Opt. Commun. 283(14), 2957–2960 (2010). [CrossRef]  

10. H. A. Haus, Waves and Fields in Optoelectronics (Prentice-Hall, New Jersey, 1984).

11. C. Li, J. F. Wu, and W. C. Xu, “Precise control of light frequency via a linear photonic crystal microcavity,” Appl. Opt. 49(14), 2597–2600 (2010). [CrossRef]  

12. A. Taflove, and S. C. Hagness, Computational Electrodynamics (Artech House, Norwood, MA, 2000)

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 (4)

Fig. 1
Fig. 1 Sketch of a typical bistable switching composed of a resonator cavity coupled to both input and output WGs.
Fig. 2
Fig. 2 Working principle diagram of a bistable switching
Fig. 3
Fig. 3 Evolvement curves of bistable switching for different pump power of the control light. The upper curve: p x = 0.0136 kW/μm; The lower curve: p x = 0.0132 kW/μm.
Fig. 4
Fig. 4 Critical pump time of the control light for δ = 4.3478. (a) Illustration for how to obtain the critical pump time of the control light. (b) Evolvement curve of transmission when the control light is shut off at t = 3522(a/c).

Tables (2)

Tables Icon

Table 1 Critical pump power of the control light with respect to different frequency detuning.

Tables Icon

Table 2 Critical pump time of the control light with respect to different frequency detuning.

Equations (18)

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

{ d A d t = [ j ( ω 0 γ | s T | 2 p 0 ) γ ] A + 2 γ 1 ( s s + s x ) γ = γ 0 + γ 1 + γ 2 + γ T P A s T = 2 γ 2 A s s + s x + s R = 2 γ 1 A ,
γ T P A = k γ 1 | s T | 2 ,
A = 2 γ 1 ( | s s | + | s x | ) i ( ω ω 0 + γ | s T | 2 / p 0 ) + γ [ e i ω t e γ t e i ω 0 t ]
S T = 4 γ 1 γ 2 ( | s s | + | s x | ) i ( ω ω 0 + γ | s T | 2 / p 0 ) + γ ( e i ω t e γ t e i ω 0 t ) ,
| A | 2 ( t ) = ( p s + p x ) 2 / 2 γ 2 ( 1 + k η | s T | 2 / 2 ) 2 η ( δ | s T | 2 / p 0 ) 2 + 1 [ e 2 γ t 2 e γ t cos ( ω 0 ω ) t + 1 ] ,
T ( t ) = | s T | 2 | s s | 2 = ( p s + p x ) 2 / p s ( 1 + k η | s T | 2 / 2 ) 2 η ( δ | s T | 2 / p 0 ) 2 + 1 [ e 2 γ t 2 e γ t cos ( ω 0 ω ) t + 1 ] ,
( p 1 + p x ) 2 p 2 ,
p x ( p 2 p 1 ) 2 .
| A | 2 p 1 T 1 / ( 2 γ 2 ) .
( p s + p x ) 2 ( 1 + k η | s T | 2 / 2 ) 2 η ( δ | s T | 2 / p 0 ) 2 + 1 [ e 2 γ t 2 e γ t cos ( ω 0 ω ) t + 1 ] T s p s ,
( p 1 + p x ) 2 = p 2 .
p 2 / p 1 ( 1 + k η | s T | 2 / 2 ) 2 η ( δ | s T | 2 / p 0 ) 2 + 1 [ e 2 γ t 2 e γ t cos ( ω 0 ω ) t + 1 ] T 1 .
T ( t ) T 1 .
p i n = ( p s + i = 1 N p x i ) 2 ,
( i = 1 N p x i ) 2 ( p 2 p 1 ) 2 .
f ( t ) exp ( i ω 0 t )
i = 1 N p x i ( p 2 p 1 ) 2 .
f ( t ) exp ( i ω 0 t ) = i g ( ω i ) exp ( i ω i 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.