Breaking the reciprocity of light propagation in a nanoscale photonic integrated circuit (PIC) is a topic of intense research, fostered by the promises of this technology in areas ranging from experimental research in classical and quantum optics to high-rate telecommunications and data interconnects. In particular, silicon PICs fabricated in processes compatible with the existing complementary metal-oxide-semiconductor (CMOS) infrastructure have attracted remarkable attention. However, a practical solution for integrating optical isolators and circulators within the current CMOS technology remains elusive. Here, we introduce a new non-reciprocal photonic circuit operating with standard single-mode waveguides or optical fibers. Our design exploits a time-dependent index modulation obtained with conventional phase modulators such as the one widely available in silicon photonics platforms. Because it is based on fully balanced interferometers and does not involve resonant structures, our scheme is also intrinsically broadband. Using realistic parameters we calculate an extinction ratio superior to 20dB and insertion loss below 3dB.
© 2013 OSA
Materials typically used to guide and manipulate light in optical fibers and photonic integrated circuits (PIC) have symmetric permittivity and permeability tensors and in the linear regime light propagation satisfies the Lorentz reciprocity theorem . One way to break reciprocity is to operate in the non-linear regime and harvest bi-stability . A more versatile approach is to guide light through a material exhibiting strong magneto-optical effect, which causes a Faraday rotation of the polarization dependent on the propagation direction [3–5]. In practice, this effect is used in commercial optical isolators and circulators, which are essential components in today’s fiber-optics systems. For example, protecting lasers from back-scattered light is critical to ensure their integrity and stability. Yet, whereas fiber-optics isolators and circulators relying on the magneto-optic effect have become widely available, a practical solution for PICs has yet to be demonstrated. Obtaining non-reciprocity with the help of magneto-optic materials by constructing a hybrid chip [6, 7] suffers from much increased fabrication complexity. This approach suppresses most benefits of silicon PICs, which rely on their entire fabrication in a CMOS-compatible process that ensures scalability, yield and low cost.
Recently, the first electrically driven non-reciprocal device in a silicon PIC was reported , based on the concept of interband photonic transitions developed earlier by Yu and Fan [9, 10] (an idea similar to the one already used for non-reciprocal mode-conversion in optical fibers ). In this device, a traveling-wave radio-frequency (RF) signal induced a time-varying and spatially non-homogeneous modulation of the refractive index in a silicon waveguide specially engineered to support two TE modes with opposite symmetries and different propagation constants and wavevectors. While the concept is elegant and no hybrid technology is needed, the implementation is so far prohibitively complicated, and the fabricated device exhibited >70dB insertion loss .
It is therefore intriguing to ask whether a simpler scheme, relying only on existing components that can be readily fabricated in a CMOS process, can be exploited to obtain non-reciprocal light propagation in a monolithic silicon PIC. Others have shown interesting progress in this direction by employing two “tandem” phase modulators to imprint a non-reciprocal frequency shift . Two limitations of this scheme are its intrinsic narrow-band operation and its modest extinction ratio of 10.8dB. Passive resonant structures can be employed to enhance intrinsic silicon non-linearities , but this approach also suffers from narrow optical bandwidth and the performance intrinsically depends on the input light power. Finally, non-reciprocal light modulation can be achieved in traveling-wave modulators, at the cost of very long devices operating at very high-speed, and with only modest extinction ratio .
In this Letter, we propose and analyze a novel non-reciprocal device utilizing balanced interferometers and multi-stage phase modulation at a frequency of a few GHz or less. Our scheme is broadband by design and can achieve an extinction ratio of better than 20dB with insertion loss smaller than 3dB, while being readily manufacturable in a standard CMOS photonic process. Our proposal relies on phase modulators separated by passive waveguide sections and driven by delayed versions of the same RF signal. The engineered RF and optical delays between modulator sections ensures non-reciprocity , and thus our scheme does not require very long traveling-wave modulators nor very high operation frequency, as was the case in previous publications [8, 13, 14], reducing simultaneously the technological complexity, the insertion loss and the achievable footprint . In a single stage, our scheme acts as a non-reciprocal modulator that can be used to encode return-to-zero (RZ) symbols. Total extinction of counter-propagating light is achieved independently of the driving signal shape and amplitude. Cascading two stages, we achieve isolation with predicted performances comparable to what is obtainable with fiber-based components, while the transmitted signal’s phase and amplitude remain non-modulated. We emphasize that none of our arguments relies on mode conversion  and that the device operates on single-mode waveguides, as the ones used in the dominant silicon-on-insulator PIC technology . To our knowledge, these are the first practical designs proposed for on-chip broadband optical non-reciprocity and they should open new avenues in the development and applications of integrated photonics.
2. Single-stage non-reciprocal modulator
2.1 System design
The design schematized in Fig. 1(a) consists of two Mach-Zehnder modulators (MZMs) driven by delayed versions of the same RF signal and separated by an optical delay line introducing a quarter-period retardation in the signal driving MZM a (left-hand side) with respect to MZM b (right-hand side). (An alternative and equivalent design based on a four-stage modulator is presented in section 5.3.) We note that tunable RF delay lines capable of 100ps delay with little distortion on 10Gb/s data stream can readily be implemented in a standard CMOS circuit  that could be wire- or flip-chip-bonded to the photonic chip. We consider a periodic function F(t): F(t) = F(t + T); with period T = 1/f; also satisfying F(t ± T/2) = -F(t), and normalized to have peak-to-peak amplitude ± 1. Examples of such functions are sine and cosine waves, as well as square waves with 50% duty cycle. An MZM is implemented by modulating the optical phase in each arm of a Mach-Zehnder interferometer. For chirp-free operation, the phase modulators are driven in push-pull mode: the upper arm experiences a phase shift φ(t) while the lower arm is driven symmetrically with a phase shift -φ(t) with respect to a constant offset. In the laboratory time-frame of reference, the optical phase modulations in MZM a and b can be written φa(t) = γ(1 ± F(t-T/4)) and φb(t) = γ(1 ± F(t)), respectively, with γ the effective phase modulation amplitude (in radians) and +/− for the upper/lower arm. By choosing the waveguide length between the two MZMs to be Lopt = T/4 c/ng (ng is the group index and c the speed of light in vacuum) light incoming from the right (ports b1, b2) travels in phase with the RF signal and therefore experiences twice the same modulation in MZMs b and a. On the contrary, for light incoming from the left (portsa1, a2), the modulation functions in MZMs a and b exhibit a π phase shift (i.e. have opposite signs). In this configuration, non-reciprocity is ensured by the presence of two phase-shifted periodic signals. To obtain high-visibility interferences the optical path difference between the two arms of the delay line must be adjustable (see below) and thus resistive heaters shall be used to finely tune their relative index shift . In the following simulations we will take f = 4 GHz, which for a group index of ng = 4.2  leads to Lopt = 4.46 mm.
2.2 Transfer matrix calculations
To give an intuitive understanding of the functioning of our device, we present first a simplified calculation of its idealized transfer matrix. Fully realistic simulations are used next to compute the accurate behavior (see section 5.1 for details). Within the coupled-mode formalism, the relationship between complex field amplitudes a1, a2and b1, b2 at the input and output of a four-port device, written in vector notation , is expressed through a 2x2 matrix. The transfer matrix of an ideal 50/50 beam splitter (directional coupler or multimode interferometer) writes:Fig. 1(a) for light propagating from a to b (left to right):Figure 1(b) presents how our system can be conceptualized as an optical circulator by assigning ports 1, 2 and 3 to a1,b2and a2, respectively (port 0 = b1 can be used for monitoring or can be terminated by a taper to avoid reflection). From the above expression for it can be seen that port 1 is perfectly isolated from light incident in port 2, independently of the precise shape and amplitude of the RF signal F(t), as long as the condition Lopt = T/4 c/ng is fulfilled.
2.3 Numerical simulations
We show the simulated behavior of a realistic device in Fig. 1(c), for a phase modulation amplitude γ = π/4 (π/2 peak-to-peak) and two different waveforms (a cosine and a bandwidth-limited square wave). In the calculations we included physically relevant effects such as waveguide loss (0.3dB/cm)  and beam splitter insertion loss (0.1dB) [21, 22], and the dynamic loss intrinsically linked to phase modulation using free-carrier dispersion effect  (see section in section 5.2). Even in the presence of losses, and independent of the signal shape, transmission from port 2 to 1 is exactly zero, demonstrating robust isolation. In contrast, transmission from port 1 to 2 is non-zero, with a maximal value of ~0.74 (or −1.3dB) and a time-averaged value of –3.27dB for a cosine signal. Improved averaged transmission is achieved by driving the modulators with a square-wave. For example, assuming a modulation bandwidth of 5f (Fig. 1(c), right panel), time-averaged insertion loss decreases to 2.17dB.
This first result is by itself remarkable and potentially useful in real systems. For example, the output at port 2 can be directly used as a clock signal, or as an information carrier in return-to-zero encoding schemes. The device can even be used to simultaneously perform the encoding by modulating the amplitude γ. Our system therefore integrates a return-to-zero modulator and a high extinction isolator into a single compact device, making use of only conventional components. We emphasize that the modulation frequency can be chosen arbitrarily, as low as permitted by waveguide propagation loss in the delay line and desired footprint, and as fast as permitted by the modulators and drivers bandwidths.
3. Double-stage isolator
3.1 System design
We now consider the device used in reversed direction (or equivalently with an opposite sign of the RF delay). As seen in Fig. 1(c), we have non-modulated transmission from port 2 to 3, independently of the driving signal parameters, yet there are spikes of transmission from port 3 to 2, preventing complete isolation. We therefore propose to cascade two identical devices and drive them with the same RF source, with a quarter-period RF delay between them as shown in Fig. 2(a). The two devices should be positioned immediately next to each other so that negligible optical delay is introduced between them. In this scheme, light passing through the transmission spikes of the first device is always rejected by the second one.
3.2 Numerical simulations
The simulations in Figs. 2(b)-2(c) demonstrate that this configuration indeed achieves non-modulated transmission of right-to-left propagating light, both in amplitude (Fig. 2(b)) and phase (Fig. 2(c)), with insertion loss of 2.9dB. Extinction is reduced compared to Fig. 1(c) but still better 20dB in the case of square-wave modulation and perfect arm balancing (14dB for cosine modulation). This figure can be further improved by increasing the modulator bandwidth-to-frequency ratio. When the optical paths in the arms of each delay line are not perfectly balanced (modulo 2π), the performances are degraded, as shown by the dashed lines in Figs. 2(b)-2(c) for an optical phase mismatch of π/10 between upper and lower arm in each stage’s delay section.
To estimate the impact of relevant parameters and imperfections on the system, we study two figures of merit: the insertion loss (IL), defined as the time-averaged transmission in the passing direction, and the extinction ratio (ER), equal to the peak transmission value in “blocking” direction divided by the IL. In Fig. 3(a), we show the dependence of IL and ER on the modulation amplitude for the single-pass and the cascaded configurations; the optimal modulation amplitude in both cases is close to γ = π/4, which yields the lowest IL for single-pass configuration and the highest ER for cascaded configuration.
Figure 3 also reports the sensitivity of the performances on the relative phase difference accumulated in the two arms of the delay line (modulo 2π). Given the thermo-optic coefficient of silicon dnSi/dT = 1.9 x 10−4 K−1 around 1.55µm, we calculate the temperature-dependent phase shift (in rad) per unit length of waveguide to be less than π/4 mm−1K−1. Controlling the phase difference within π/10 can thus be achieved by tuning the temperature over a 0.4 mm section with 1 K accuracy, well within reach of existing technology . Yet some feedback circuit may be needed to ensure stable operation.
To conclude, we estimate the expected optical bandwidth. As the system relies on a series of fully balanced interferometers, it is by design wavelength insensitive. Yet, several second-order effects may limit the actual operating wavelength range. Directional couplers usually have limited bandwidth, so we opt for MMIs [21, 26] (see Fig. 1(a)), for which uniform splitting ratio over 94 nm have been demonstrated . Second, due to group velocity dispersion the optical delay between MZM a and b is actually wavelength dependent. Using the dispersion measured in  and a waveguide length of 4.5 mm, we estimate a delay variation of less than 1 ps over more than 100 nm bandwidth around 1550 nm. This is much smaller than the modulation period (250 ps here) and has therefore negligible impact on performance (the ratio ~1/250 is independent of the particular modulation frequency). This variation could be further reduced by tailoring the dispersion . The limiting factor may eventually come from the wavelength dependence of the plasma dispersion effect , but this too could be easily compensated for by tuning the modulation amplitude γ.
5.1 Transfer matrix calculations
In this section we detail the calculations of the transfer matrices used to produce the simulations shown in the paper. Within the coupled-mode formalism the complex field amplitudes a1, a2and b1, b2 at the input and output of a four-port device are written as vectors ; and the input-output relationship takes the form of a complex 2x2 matrix. We note that all following calculations are up to a phase factor common to both arms of the circuit. We list below the elementary components into which our non-reciprocal devices (Fig. 1(a) in main text) can be decomposed, and derive their respective transfer matrices.
- • In a Mach-Zehnder modulator (MZM) driven in push-pull mode, the upper arm experiences a phase shift φ(t) while the lower arm is driven symmetrically with a phase shift -φ(t). This is obtained by applying an offset bias and modulating each arm with opposite voltage signs around this offset. We also introduce an arm imbalance quantified by a phase delay ∂φin the upper arm relative to the lower arm. We furthermore account for dynamic losses caused by free carrier absorption, which always accompanies phase modulation by plasma dispersion effect (see next section). They are characterized by an excess loss β (in dB) per π phase shift. Under increasing free carrier concentration, the refractive index (and therefore the optical phase delay) decreases, so that, with reference to the intrinsic waveguide delay, the phase modulation term writes . The dynamic loss increases under carrier injection and is thus expressed as ; where we have defined . In the opposite arm the factor 1 + F(t) is replaced by 1-F(t) to simulate push-pull operation. The transfer matrix for the modulator section in an MZM is then:
Our model can be used to study all relevant optical effects in the linear regime, as well as arbitrary RF signals (including bandwidth limitation, by appropriate choice of the function F(t)), and of course the impact of non-ideal physical implementation such as MZ arm imbalance, asymmetric splitting ratios, etc.
5.2 Dynamic loss and free carrier modulation
Here we detail how we obtain the value for the dynamic modulation loss that we use in the simulations. High-speed (>20 GHz) phase modulation in Si is achieved through the “plasma dispersion” effect, i.e. the change in the complex dielectric function due to a change in free carrier concentration. Soref and Bennett derived the magnitude of this effect in 1987, using a combination of previous experimental data:
Using these relations we can derive the excess loss [in dB] caused by free carrier absorption for a change of index over a propagation length L, when only free electrons are considered:
5.3 Alternative proposal
We propose a second scheme pictured in Fig 4. Here, a single MZI is built with four phase modulators in each arm, pair-wise driven in push-pull mode by the signal F(t). Again, the electric signal is delayed from right to left by a quarter-period between each successive pair of modulators, and an optical path length Lopt = T/4c/ng is introduced to ensure that light launched from the right keeps a fixed relationship with the phase of the electric drive, thus experiencing a total optical phase modulation . Light propagating from left to right sees a π retardation in the electric signal between each successive modulator, thus accumulating a zero net phase shift relative to the other arm:. The transfer matrices for this configuration are therefore:
We now recall the transfer matrices obtained for the design of Fig. 1(a):
In order to give an intuitive understanding of the second device operation, we consider an ideal square-wave signal (no bandwidth limitation) and split each period into 4 equal time intervals Δti, i = 1..4 (see Fig. 5) with the convention Δti+4 = Δti, i any integer. We denote the phase shift applied on the upper (resp. lower) arm by modulator number j at time i by the matrix element φij. (resp. -φij). With this notation, light traveling from left to right accumulates a phase shift in the upper/lower arm of (with the convention φi,j+4 = φi,j). This corresponds to summing over the diagonals of the square matrix φij (i,j=1..4). For light traveling in “backward” direction (from right to left here), the accumulated phase shift writes: , corresponding to a sum over the anti-diagonals. With this insight, it is readily seen that the following matrix
Using this matrix approach, it is easy to see that 4 modulation sections is the minimum number to achieve optical isolation with this scheme.
The authors would like to thank Gernot Pomrenke, of the Air Force Office of Scientific Research, for his support under the OPSIS (FA9550-10-1-0439) PECASE (FA9550-10-1-0053) and STTR (FA9550-12-C-0079) programs, and would like to thank Mario Paniccia and Justin Rattner, of Intel, for their support of the OpSIS program. Professor Hochberg would like to acknowledge support from the Singapore Ministry of Education ACRF grant R-263-000-A09-133.
References and links
1. S. Fan, R. Baets, A. Petrov, Z. Yu, J. D. Joannopoulos, W. Freude, A. Melloni, M. Popović, M. Vanwolleghem, D. Jalas, M. Eich, M. Krause, H. Renner, E. Brinkmeyer, and C. R. Doerr, “Comment on “Nonreciprocal Light Propagation in a Silicon Photonic Circuit”,” Science 335(6064), 38, author reply 38 (2012). [CrossRef] [PubMed]
6. R. L. Espinola, T. Izuhara, M.-C. Tsai, R. M. Osgood Jr, and H. Dötsch, “Magneto-optical nonreciprocal phase shift in garnet/silicon-on-insulator waveguides,” Opt. Lett. 29(9), 941–943 (2004). [CrossRef] [PubMed]
7. L. Bi, J. Hu, P. Jiang, D. H. Kim, G. F. Dionne, L. C. Kimerling, and C. A. Ross, “On-chip optical isolation in monolithically integrated non-reciprocal optical resonators,” Nat. Photonics 5(12), 758–762 (2011). [CrossRef]
9. Z. Yu and S. Fan, “Complete optical isolation created by indirect interband photonic transitions,” Nat. Photonics 3(2), 91–94 (2009). [CrossRef]
13. L. Xu and H. K. Tsang, “Nonreciprocal Optical Modulation for Colorless Integrated Optical Transceivers in Passive Optical Networks,” Opt. Commun. Netw. 2(3), 131–136 (2010). [CrossRef]
14. S. Bhandare, S. K. Ibrahim, D. Sandel, Z. Hongbin, F. Wust, and R. Noe, “Novel nonmagnetic 30-dB traveling-wave single-sideband optical isolator integrated in III/V material,” IEEE J. Sel. Top. Quantum Electron. 11(2), 417–421 (2005). [CrossRef]
16. V. Liu, D. A. B. Miller, and S. Fan, “Ultra-compact photonic crystal waveguide spatial mode converter and its connection to the optical diode effect,” Opt. Express 20(27), 28388–28397 (2012). [CrossRef] [PubMed]
17. R. Soref, “The past, present, and future of silicon photonics,” IEEE J. Sel. Top. Quantum Electron. 12(6), 1678–1687 (2006). [CrossRef]
18. J. Buckwalter and A. Hajimiri, “An Active Analog Delay and the Delay Reference Loop,” Proceedings of the IEEE Radio Frequency Integrated Circuits (RFIC) Symposium, 17–20 (2004). [CrossRef]
19. R. L. Espinola, M. C. Tsai, J. T. Yardley, and R. M. Osgood Jr., “Fast and low-power thermooptic switch on thin silicon-on-insulator,” IEEE Photon. Technol. Lett. 15(10), 1366–1368 (2003). [CrossRef]
20. E. Dulkeith, F. Xia, L. Schares, W. M. J. Green, and Y. A. Vlasov, “Group index and group velocity dispersion in silicon-on-insulator photonic wires,” Opt. Express 14(9), 3853–3863 (2006). [CrossRef] [PubMed]
21. Z. Sheng, Z. Wang, C. Qiu, L. Li, A. Pang, A. Wu, X. Wang, S. Zou, and F. Gan, “A Compact and Low-Loss MMI Coupler Fabricated With CMOS Technology,” IEEE Photon J. 4(6), 2272–2277 (2012). [CrossRef]
22. R. Halir, I. Molina-Fernandez, A. Ortega-Monux, J. G. Wanguemert-Perez, X. Dan-Xia, P. Cheben, and S. Janz, “A Design Procedure for High-Performance, Rib-Waveguide-Based Multimode Interference Couplers in Silicon-on-Insulator,” Lightwave Technology, Journalism 26, 2928–2936 (2008).
23. P. Dong, W. Qian, S. Liao, H. Liang, C.-C. Kung, N.-N. Feng, R. Shafiiha, J. Fong, D. Feng, A. V. Krishnamoorthy, and M. Asghari, “Low loss shallow-ridge silicon waveguides,” Opt. Express 18(14), 14474–14479 (2010). [CrossRef] [PubMed]
24. R. A. Soref and B. R. Bennett, “Electrooptical effects in silicon,” IEEE J. Quantum Electron. 23(1), 123–129 (1987). [CrossRef]
26. L. B. Soldano and E. C. M. Pennings, “Optical multi-mode interference devices based on self-imaging: principles and applications,” J. Lightwave Technol. 13(4), 615–627 (1995). [CrossRef]
28. A. C. Turner, C. Manolatou, B. S. Schmidt, M. Lipson, M. A. Foster, J. E. Sharping, and A. L. Gaeta, “Tailored anomalous group-velocity dispersion in silicon channel waveguides,” Opt. Express 14(10), 4357–4362 (2006). [CrossRef] [PubMed]