We characterized the sputtered Sb2S3 thin films to obtain the refractive indices and verify the crystallization capability (annealed under 325°C in Nitrogen environment for 10 mins) of the as-deposited a-Sb2S3 (Supplementary Section S1). The silicon photonic devices (Fig. 1-Fig. 3) were designed to operate at the telecommunication O-band (1260–1360 nm) and fabricated on a standard silicon-on-insulator (SOI) wafer with 220 nm silicon and 2 µm buried oxide. The 500-nm-wide waveguides are fabricated by partially etching 120-nm silicon. We then deposited 450-nm-wide Sb2S3 onto the SOI chip via sputtering. The slightly smaller width than the waveguide is to compensate for electron beam lithography (E-beam) overlay tolerance. Our simulation results show a change in effective index \({{\Delta }n}_{eff}\approx 0.018\) between a- and c-Sb2S3 (Supplementary Section S2). The Sb2S3 films are electrical controlled via on-chip silicon PIN micro-heaters20,22. The p++ and n++ doping regions were designed 200 nm away from the waveguide to avoid free-carrier absorption loss22, indicated in the scanning electron microscope (SEM) images with false colors in Fig. 1c, Fig. 2c, and Fig. 3c. The Sb2S3 stripes are encapsulated with 40 nm of Al2O3 grown by atomic layer deposition (ALD) under 150°C. This conformal encapsulation is critical to prevent Sb2S3 from oxidation and thermal reflowing, and thus is essential to attain high endurance. To show that our Sb2S3-clad silicon photonic platform is versatile and compatible with most PIC components, we demonstrate three widely used PIC components: (1) a microring resonator to show low-loss tuning of cavities, (2) a balanced MZI to demonstrate a full π phase shift and a broadband operation, and (3) an asymmetric directional coupler to create a compact programmable unit (see simulation results in Supplementary Section S2).
Nonvolatile microring switch integrated with Sb2S3 phase shifter
We deposited 10-µm-long, 20-nm-thick Sb2S3 on a micro-ring resonator with 30 µm radius (Figs. 1a-c). The free spectral range (FSR) is ~ 2.42 nm (Fig. 1d), and the bus-ring gap is 280 nm to achieve a near-critically coupled device. We switched the as-deposited a-Sb2S3 to c-Sb2S3 on the microring resonator by applying three 1.6 V, 200-ms-long pulses with ~ 3.4 mJ energy (or SET pulses) separated by 1 second and then re-amorphized the material via three 7.5 V, 150-ns-long pulses with ~ 56 nJ energy also separated by 1 second (RESET pulses). The unit energy density (energy/total Sb2S3 volume) for switching Sb2S3 is then estimated as 0.6 \(fJ/n{m}^{3}\) (38 \(pJ/n{m}^{3}\)) for amorphization (crystallization). We note that the 200-ms-long SET pulse is indeed significantly longer than other reported PCMs, such as GST (50 µs22 or 100 µs21) and Sb2Se3 (5 µs16 or 100 µs21), and causes a large crystallization energy and energy density. But we found that the SET pulse duration can be reduced to around 100 µs after the first few cycles (see Method and Supplementary Section S12). Therefore, the crystallization power is reduced to ~ 1.7 \(\mu\)J (19 \(fJ/n{m}^{3}\)) after the initial conditioning, comparable to other PCMs16,22. With the \(100 \mu s\) crystallization pulse condition, this device could be operated at ~ kHz speed for a complete SET/RESET cycle. The slow crystallization, however, can allow amorphizing a large volume of Sb2S3 (Supplementary Section S13). Figure 1d shows a resonance shift of ~ 0.394 nm upon switching the 10 µm Sb2S3 from the a- (blue) to the c-phase (orange), corresponding to a π-phase shift length Lπ of ~ 30.7 µm, significantly shorter than the 1-millimeter Lπ of ferroelectric non-volatile phase shifter32. The SET and RESET processes were repeated for 10 cycles. The shift in resonance is highly repeatable, as suggested by the slight standard deviation (Fig. 1d). The excess loss from a-Sb2S3 is negligible24, and in c-Sb2S3 hybrid waveguides the loss is estimated to be 0.024 dB/µm (0.72 dB/π), which is three times larger than our simulation result (0.26 dB/\(\pi\)) (Supplementary Section S2). We attribute this excess loss to the scattering from the Sb2S3 thin film due to non-uniform deposition/liftoff and local crystal grains in c-Sb2S3 (Supplementary Section S3). We verified that the loss due to mode mismatch at the transition between the bare silicon waveguide and the Sb2S3-loaded waveguide is small (~ 0.013 dB/facet), consistent with the fact that the thin Sb2S3 film should not significantly change the mode shape (Supplementary Section S2).
Nonvolatile Mach-Zehnder switch integrated with Sb2S3 phase shifter
Figures 2a-c show a balanced MZI operating at wavelengths between 1,320 nm and 1,360 nm with both arms covered with 30-µm-long 20-nm-thick Sb2S3. A multimode interferometer with a 50:50 splitting ratio was designed and fabricated, as shown in Fig. 2d (simulation results in Supplementary Section S2). Initially, the light comes out mainly from the bar port with an extinction ratio of ~ 13 dB (Fig. 2e). The light coming out from the bar port instead of the cross port in this balanced MZI can be explained by the random phase errors in two arms due to fabrication imperfection, especially that in the S-bend. One can overcome such imperfection by exploiting a wider waveguide to improve fabrication robustness33. Alternatively, this random phase error can be corrected using Sb2S3 for post-fabrication trimming, as we show later in this paper. The Sb2S3 on one arm was switched by two 1.7 V, 200-ms SET electric pulses with an energy of 11.6 mJ (density: 42 \(pJ/n{m}^{3}\)) to provide a full π phase shift. An 8.1 V, 150-ns short RESET electric pulse with an energy of 197 nJ (density: 0.73 \(fJ/n{m}^{3}\)) switched the device back to the initial state. Figures 2e and 2f show the transmission spectra normalized to a reference waveguide when the Sb2S3 film is in the amorphous and crystalline phases, respectively. The c-Sb2S3 displays a complete spectrum flip, showing a bar state with an extinction ratio of 15 dB. We then recorded the bar port transmission at 1,330 nm for 100 switching events without device degradation (Supplementary Section S4).
Compact asymmetric directional coupler switch
We also designed and fabricated a compact asymmetric directional coupler (coupling length Lc ≈ 79 µm) (simulation in Supplementary Section S2), as shown in Figs. 3a-c. The coupler consists of two waveguides with different widths. The narrower 409-nm-wide waveguide (hybrid waveguide) was capped with 20-nm-thick Sb2S3 and designed to allow phase match with the wider 450-nm-wide waveguide (bare waveguide) for c-Sb2S334. As such, the input light could completely couple to the cross port in one coupling length. Once Sb2S3 is switched to the amorphous state, the effective index of the hybrid waveguide changes while the bare waveguide remains the same. The resulting phase mismatch changes the coupling strength and coupling length. Then, co-optimizing the gap and waveguide length permits a complete bar transmission. The unique c-Sb2S3 phase matching approach, instead of a-Sb2S318,20,34,35, allows a more symmetric performance regardless of the input port (Supplementary Section S2), crucial for a 2 × 2 device. If phase mismatch happens in the c-Sb2S3 state, the slight loss of c-Sb2S3 on one of the waveguides will result in different bar state insertion loss when the light goes from different input ports. We note that the 79-µm coupling length can be potentially reduced (~ 34 µm) by depositing a thicker (50 nm) Sb2S3 to provide stronger refractive index modulation.
Figures 3d and 3e show the transmission spectra for a- and c-Sb2S3, switched with three 9.6 V, 500 ns, 922 nJ RESET, and 2.7 V, 200 ms, 29.2 mJ SET pulses, respectively. The energy density for amorphization (crystallization) is 1.28 \(fJ/n{m}^{3}\) (40 \(pJ/n{m}^{3}\)). The insertion losses are 2 dB (0.5 dB), and the extinction ratios are around 10 dB (11 dB) for a (c)-Sb2S3). The unexpected high insertion loss when the Sb2S3 is in the amorphous state can be attributed to several factors, including gap discrepancy, Sb2S3 overlay deviation, or the cross-port grating coupler fabrication imperfection. To estimate the actual loss of the device, we apply the c-Sb2S3 loss extracted from the ring resonator to the simulation and calculate this device’s insertion loss to be ~ 0.1 dB (0.9 dB) for a(c)-Sb2S3 (Supplementary Section S2).
Figure 4 shows 1,600 switching events for the asymmetric directional coupler. Limited by our measurement setup, we separately measured the cross (Fig. 4a) and bar ports (Fig. 4b). The higher insertion loss (~ 1 dB) at around event 500 was due to optical fiber misalignment. We note that almost no performance degradation occurred at the end (Supplementary Section S5); hence, 1,600 switching events are not the limit of this device. The cross-port transmission shows a relatively large variation for a-Sb2S3 (Fig. 4a blue scatterers, from ~ -15 dB to ~ -35 dB), which was caused by incomplete amorphization or thermal reflow of the Sb2S3 film. Since the plot is on a logarithmic scale, such a large variation in a-Sb2S3 (due to higher transmission) is not visible in Fig. 4b.
Multilevel 5-bit operation with dynamic electrical control
Our Sb2S3-Si integrated structures further show a stepwise multilevel operation up to 32 levels with dynamic pulse control. In Fig. 5a, we show the multilevel transmission at both cross and bar ports of an asymmetric directional coupler while sending in RESET 10 V, 550 ns, 1.1 \(\mu\)J pulse every other second to amorphized the Sb2S3. We started the experiment with “coarse” tuning, where unoptimized, identical pulses were sent to partially amorphize the c-Sb2S3 device to demonstrate multiple levels. The asymmetric directional coupler was originally in the “cross state” (red region). After one partial amorphization pulse, it was reconfigured into an intermediate state (orange region), where light comes out from both cross and bar ports. After six pulses, a complete “bar state” (green region) was achieved. We repeated this experiment five times for each port and plotted the average transmission levels and the standard deviation. The variation is attributed to the stochastic phase change process using electrical controls36. We also experimentally tested the partial crystallization of the device and multilevel operation (Supplementary Section S5), which is based on the growth-dominant nature of Sb2S3 crystallization process34. In the following experiments, we mainly focused on partial amorphization because of lower energy consumption and more operation levels with finer resolution.
Such stepwise multilevel operation by applying identical pulses is distinctly different from previously reported multilevel operations in GST20,22 and Sb2Se316,21, where different voltage amplitudes or pulse duration were used to access multiple levels during amorphization. The pulse-number-dependent behavior is quite counterintuitive: one expects that after the first amorphization pulse, the thin Sb2S3 film would have reached its new equilibrium phase. Moreover, since the thermal processes relax within 10 µs (Supplement Section S6), each pulse is independent due to the relatively long one-second interval. As a result, the subsequent identical and separated pulses should not further change the material phase. To understand the origin of the multi-level operation, we closely inspected four partially amorphized Sb2S3 devices under the microscope. We observed a few separate patches (Supplement Section S7) and a region that grew with more voltage pulses. As reported in some literature, one possibility could have been that a- and c-Sb2S3 have significantly different thermal conductivities and specific heat capacities. But we measured the thermal conductivities to be similar (a-Sb2S3: 0.2 W/m/K; c-Sb2S3: 0.4 W/m/K, see Methods), and hence, we ruled out this as a possible explanation. We hypothesize that this unique behavior comes from Sb2S3’s multiple crystalline phases. Sb2S3 has at least two distinct crystalline phases37, which may differ in the amorphization conditions. The partial amorphization pulse can cause amorphization in the hottest region, but at the lower temperature region, it may cause phase transition to the other crystalline phase. These regions get amorphized in subsequent pulses, resulting in multilevel operation.
An even finer multilevel operation was realized by monitoring the transmission level and dynamically changing the pulse conditions slightly. Here, we demonstrate on-demand 5-bit operation in a quasi-continuously tunable directional coupler, as shown in Fig. 5b. We dynamically controlled the partial amorphization pulses to have slightly lower, near-identical voltages (ranging from 9.65 V to 9.85 V) and obtained up to 32 levels. Figure 5b demonstrates 5-bit operation (32 distinct levels) with a target resolution of 0.5 dB per level step at 1,340 nm (see the detailed pulse conditions in Supplementary Section S8). We emphasize that dynamic control is necessary to mitigate the stochastic nature of electrically controlled PCMs36, hence essential for a reliable many-level operation. In Fig. 5b, a linear fit shows a slope of -0.50 dB per step and a standard deviation of 0.16 dB among five experiments, indicating a repeatable operation. While the 5-bit operation of GST was shown using laser pulses38, our demonstration is thus far the highest number of operating levels reported using electrical control in PCMs-based photonics. Moreover, our multilevel operation does not require sophisticated heater geometry engineering, such as the segmented doped silicon heater design16, and solely relies on the unique phase-change dynamics of Sb2S3.
Random phase error correction in balanced MZIs exploiting multilevel operation
Finally, exploiting the multilevel operation, we corrected random phase errors in a balanced MZI. A perfectly balanced MZI should initially be in an all-cross state. However, random phase errors due to fabrication imperfections can easily build up to a phase error of π, making the initial state unpredictable. Therefore, balanced MZIs usually require extra calibration33. For example, Fig. 6a shows the transmission spectra of a phase error corrupted balanced MZI, indicating a high bar transmission. Both arms of the MZI are loaded with 40-µm-long Sb2S3 film to guarantee a phase tuning range of more than ± π. The corrected MZI spectra by multilevel tuning of Sb2S3 are demonstrated in Fig. 6b, showing a pure “cross state” with a high extinction ratio of 24 dB. The trimming process is shown in Fig. 6c. We sent in a partial amorphization pulse (8.8 V, 150 ns, 232 nJ) every other second, which gradually increased the portion of amorphous Sb2S3, resulting in quasi-continuous changes of the bar (blue) and cross (orange) transmission (at 1,340 nm). The correction finishes once a bar transmission minimum is reached, indicated by the red arrows in Fig. 6c, and the spectra reported in Fig. 6b were then measured. Further pulses increase bar transmission because of the over-compensated phase. We performed the same experiment three times, indicated by different colored regions in Fig. 6c. Complete phase error correction was observed in all three instances, exhibiting excellent repeatability of our trimming process. Note that binary tuning cannot accomplish this task due to the random initial phase error. Even multilevel operations with limited discrete-level resolution can cause over- or under-corrected phase error, ultimately determining the trimming resolution. We highlight that this method requires zero static energy supply once the phase error is corrected, as the phase transition in Sb2S3 is non-volatile (> 77 days, Supplementary Section S9). Thanks to the relatively fine operation levels, slight over-tuning does not significantly affect the performance. The trimming resolution can be further improved using dynamically changed, near-identical electrical pulses, as shown in the previous Section. Moreover, if the phase error is over compensated, the device could be tuned back with partial crystallization pulses (Supplementary Section S5). In the future, our trimming process can potentially be fully automated by real-time adjusting the pulse numbers according to the measured transmission.