We demonstrate by numerical simulations of coupled generalized nonlinear Schrödinger equations that the input pulse duration as well as the retarded material response have a crucial impact on the properties of an ultrafast nonlinear optofluidic fiber coupler. This device is composed of two waveguides in close vicinity embedded in a photonic crystal fiber which are filled with a highly nonlinear liquid. We show that in particular the critical peak power above which the coupling between the waveguides is suppressed increases dramatically for short input pulses and long characteristic response times of the liquid. We establish a simple model which describes these effects with high accuracy.
© 2011 Optical Society of America
With the recent advances in preparing selectively liquid-filled photonic crystal fibers (PCFs) [1,2] it has become possible to fabricate new coupling devices  based on two or more nearby waveguides filled with highly nonlinear liquids embedded in the photonic structure of a PCF  (see Fig. 1a). These devices are of high interest for the studies of the nonlinear coupling behavior since they offer the opportunity to tune the coupling strength, the dispersion, as well as the nonlinear parameter by adjusting the fiber geometry and the used materials accordingly. In addition to these advantages the PCF forms a convenient scaffold which ensures a constant core-to-core distance over several centimeters of propagation length, hardly achievable by usual liquid-filled capillaries put in close vicinity of each other.
The Kerr index of the highly nonlinear liquids, among them CS2, CCl4, nitrobenzene, or toluene, is up to 100 times larger than in the case of fused silica, which reduces the required nonlinear switching power so that low-power nonlinear coupling becomes possible.
Earlier experimental or theoretical investigations of nonlinear fiber couplers have mainly focused on a continuous-wave (cw) input [3, 5] or on nonlinear soliton switching between two waveguides made of fused silica [6, 7]. In the work presented in  and , the retarded response of fused silica due to the Raman effect and its influence on the coupling was investigated. However, when using highly nonlinear liquids as the carrier medium, their retarded responses differ considerably from each other and that of fused silica in strength as well as in their characteristic time scales . Due to this reason a general analysis of the effects induced by the retarded response without focussing on a specific medium is desirable for the design of an optofluidic coupler.
In this paper we investigate the influence of the input pulse duration and the retarded material response on the properties of an ultrafast nonlinear optofluidic fiber coupler. We find that in particular the critical input peak power PC, above which the coupling is suppressed, may increase dramatically, depending strongly on both parameters. We present a simple model capable to explain the findings.
By introducing a hypothetic medium with fixed dispersion and nonlinearity but with a variable retarded response we are able to separate the effects induced by the retarded response from other influences . The nonlinear propagation in two identical nearby waveguides filled with this hypothetic medium is known to be well described by two coupled generalized nonlinear Schrödinger equations (NLSEs) :15].
The response of the medium is governed via R(t), defined by
This phenomenological design of R(t) has been shown to represent a good approximation  to measured retarded responses of commonly used media (Fig. 1b), allowing us to model all regarded liquids with only one free parameter apart from the fraction fR to keep the parameter space small. This free parameter we call tR, namely the characteristic time, and use it to scale the time axis.
The PCF in which we embed the hypothetic medium has a hole diameter of 2.7 μm and a hole-to-hole distance of Λ = 5.6 μm. This leads to a waveguide separation of . For the calculation of the dispersion coefficients  we use the Kerr-corrected refractive indices of CCl4 , a central wavelength of λ0 = 1040 nm and include the fiber geometry. This results in β2 = 5.5787 × 10−27 s2/m, β3 = 9.6736 × 10−42 s3/m, β4 = −3.8758 ×10−57 s4/m, and β5 = 1.6628 × 10−72 s5/m for the linear coefficients. Since β2 is positive, we operate in the normal dispersion regime, and soliton formation is not expected. With the Kerr index of CCl4 being n2 = 15 × 10−20 m2/W , we obtain for the nonlinear parameter γ = 0.0265 (Wm)−1. The coupling strength κ is set to 174.53 m−1, leading to a linear coupling length LC = π/2κ = 9 mm . The propagation length z0 is set to 4LC = 36 mm. The cross-phase modulation coefficient becomes fj(3–j) = 8 ×10−4 .
In principle, the waveguides can also be operated in the anomalous dispersion regime. However, in that case due to the high nonlinearity of the regarded liquids, soliton fission and supercontinuum generation is likely. A full analysis of the propagation and switching behavior in the anomalous dispersion regime including the supercontinuum generation requires a more detailed discussion which we will present in a future publication.
To solve Eq. 1 we implement the symmetrized split-step Fourier method , which is well suited for the pulse durations and spectral widths we expect for the results. As the initial pulse we choose a sech-pulse given by
Fig. 2 shows a selection of evolutions of the pulse energy contained in waveguide 1, normalized to the input pulse energy, for different input pulse durations and retarded response parameters. Common for all plots is the following behavior: Starting from low input powers, the coupling length increases for increasing input peak power. This corresponds to a decreasing coupling strength, which can be explained by a higher confinement and hence lower overlap of the fields. This higher confinement is caused by the increasing index contrast between the fiber core and cladding due to the Kerr effect. Once a certain input peak power is reached, namely the critical power PC, the coupling is suppressed and the pulse energy remains in waveguide 1. In case of cw input, the critical power can be calculated by P̃C = 4κ/γ ≈ 26.4 kW , which in contrast to the results of Fig. 2 is independent of T0 and the retarded response. We observe a shift to higher values of PC if the pulse duration T0 is decreased. Furthermore, the shift of the critical power to higher values becomes stronger for a more distinctive retarded response, i.e., for increasing fR and tR.
Performing more simulations with varying input pulse durations and retarded response parameters, we show in Fig. 3a the critical powers PC extracted manually from the simulation results and the corresponding evolution plots. Clearly a strong deviation from the cw-case is observed. PC may exceed the cw-case P̃C by a factor of more than 3 for short pulses and high retardation.
To quantify our findings we expand the cw-case critical power by taking the retarded response into account:
Inserting Eq. 2, we obtain for the correction divisorFig. 3b shows PC obtained by this simple model, which exhibits an excellent agreement with the critical powers extracted by the simulations for a small retarded contribution, and gives still a good qualitative agreement for larger retarded contributions. We attribute the occurring deviations to the fact that our model does not include dispersion terms and self-steepening, and that it takes into account only the input field envelope at z = 0.The correction divisor DC is plotted in dependence of tR/T0 and fR in Fig. 4.
To summarize, we have shown by numerical simulations that the input pulse duration as well as the retarded material response have a large influence on the properties of an ultrafast nonlinear optofluidic fiber coupler. In particular, the critical power may have to be corrected by a factor larger than 3 for commonly used laser sources and liquids. These findings are supported by a simple model and might aid further investigations on design and experimental realization of an ultrafast nonlinear optofluidic fiber coupler.
The authors would like to thank the Landesgraduiertenförderung of Baden-Württemberg and DFG (Open Access Publishing) for support of this work.
References and links
1. M. Vieweg, T. Gissibl, S. Pricking, B. J. Eggleton, D. C. Wu, B. T. Kuhlmey, and H. Giessen, “Ultrafast nonlinear optofluidics in selectively liquid-filled photonic crystal fibers,” Opt. Express 18, 25232–25240 (2010). [CrossRef] [PubMed]
2. B. Kuhlmey, B. J. Eggleton, and D. K. Wu, “Fluid-filled solid-core photonic bandgap fibers,” J. Lightwave Technol. 27, 1617–1630 (2009). [CrossRef]
3. St. M. Jensen, “The nonlinear coherent coupler,” IEEE Trans. Microwave Theory Tech.30, 1568–1571(1982). [CrossRef]
4. P. St. J. Russell, “Photonic crystal fibers,” Science 17, 358–362 (2003). [CrossRef]
8. B. A. Malomed, I. M. Skinner, and R. S. Tasgal, “Solitons in a nonlinear optical coupler in the presence of the Raman effect,” Opt. Commun. 139, 247–251 (1997). [CrossRef]
9. Y. Wang and W. Wang, “Study of ultrafast pulse coupling dynamics considering retarded nonlinear response and self-steepening effects,” J. Lightwave Technol. 24, 1041–1047 (2006). [CrossRef]
12. P. Wiewior and C. Radzewicz, “Dynamics of molecular liquids studied by femtosecond optical Kerr effect,” Opt. Appl. 30, 103–120 (2000).
13. K. Itoh, Y. Toda, R. Morita, and M. Yamashita, “Coherent optical control of molecular motion uUsing polarized sequential pPulses,” Jpn. J. Appl. Phys. 436448–6451 (2004). [CrossRef]
14. T. F. Laurent, H. Hennig, N. P. Ernsting, and S. A. Kovalenko, “The ultrafast optical Kerr effect in liquid fluoroform: an estimate of the collision-induced contribution,” Phys. Chem. Chem. Phys. 2, 2691–2697 (2000). [CrossRef]
15. G. P. Agrawal, Nonlinear Fiber Optics (Academic Press, 1995).
16. H. H. Marvin, “The selective transmission and the dispersion of the liquid chlorides,” Phys. Rev. 34, 161–186 (1912).
17. P. P. Ho and R. R. Alfano, “Optical Kerr effect in liquids,” Phys. Rev. A 20, 2170–2187 (1979). [CrossRef]
18. K. Okamoto, Fundamentals of Optical Waveguides (Academic Press, 2005).