Triangular cross-section grating couplers for integrated quantum nanophotonic hardware in silicon carbide
Abstract
Triangular cross-section nanodevices are among the leading approaches for integrating color centers with photonics for applications in quantum information processing. We design periodic and aperiodic fishbone triangular grating couplers in silicon carbide. We optimize the designs for achieving up to collection efficiency from color center integrated triangular devices to a microscopy system. Using an ion beam angle etching process, we demonstrate proof-of-principle fabrication of the designed devices for future implementation in wafer-scale quantum nanophotonics.
I Introduction
Quantum networks require a medium that can store quantum information for a long time and provide an interface to transfer the information reliably over long distances. Photons are excellent quantum information carriers as they experience low loss and less decoherence when transmitted through optical fiber networks, making it easier to integrate quantum networks into existing infrastructure Kimble (2008); Monroe (2002); Azuma et al. (2023). Color centers, luminous defects in wide band gap semiconductors Norman et al. (2021), have sparked interest within the quantum community due to their spin-photon interface where the electron or nuclear spin acts as a stationary qubit and spin-entangled photon acts as a flying qubit Hensen et al. (2015). For quantum networks and distributed quantum information processing (QIP), color centers in silicon carbide (SiC) are particularly interesting as they generate spin-entangled photons near telecom band wavelengths Majety et al. (2022) with seconds long spin coherence times Anderson et al. (2022) and ultra-narrow inhomogeneous broadening Cilibrizzi et al. (2023). Moreover, compatibility with silicon-based device fabrication processes and availability of high-purity wafer-scale substrates make SiC a potential platform for scalable QIP hardware.
Color centers are required to be integrated with nanophotonic devices to perform on-chip QIP Castelletto et al. (2022); Majety et al. (2022). In this regard, suspended photonics paves the way for higher optical confinement due to maximum refractive index contrast with the surrounding medium. Previous efforts in fabricating color center integrated SiC suspended devices either limited the performance of color center optical coherence Lukin, Guidry, and Vučković (2020) or encountered scaling issues Majety et al. (2022). As color centers exhibit excellent spectral stability in bulk substrates Udvarhelyi et al. (2019), the angle-etching method has emerged as a promising technique to fabricate bulk freestanding nanostructures for color center integration. Angle-etching can be implemented either by a Faraday cage or by ion beam etching (IBE). In this method, ions are directed at an angle toward the substrate to undercut the nanophotonic structures which results in a triangular cross-section profile. SiC color center photonics and spintronics in this non-conventional triangular geometry have demonstrated robust performance for QIP applications Majety et al. (2021); Babin et al. (2022); Majety et al. (2023); Saha, Majety, and Radulaski (2023); Majety et al. (2024a).
Recently, we have developed a reactive ion beam etching (RIBE) process for wafer-scale integration of triangular cross-section photonic devices with color centers in SiC Majety et al. (2024b). An on-chip wafer-scale testing mechanism is required for characterizing color center interaction with triangular shaped nanophotonic hardware. Although edge couplers offer highly efficient coupling over a large bandwidth, they do not support wafer-level testing and have complex post-fabrication steps with limited access to the circuits Cheng et al. (2020); Marchetti et al. (2019). Vertical grating coupler (GC) is a convenient choice for flexible integration with other photonic elements on the same chip, consequently, allowing for wafer-scale testing. Fiber-based GCs exhibit excellent coupling efficiency at the cost of large footprints and precise fiber alignment Marchetti et al. (2019). As cryogenic environment is required for color center characterization Norman et al. (2024), free-space GCs in combination with microscope objective have proven to be a very practical solution Faraon et al. (2008); Zhou et al. (2018). Typical microscopy systems have a higher numerical aperture (NA) compared to optical fibers and provide wider bandwidth Zhu et al. (2017). Moreover, due to the absence of mode mismatch issues between the fiber core and the GC, free-space GCs can be optimized for a compact layout.
In this work, we conduct a thorough study on free-space SiC triangular cross-section suspended GC in light of color center quantum nanophotonics. Fig. 1(a) shows a perspective illustration of color center coupling characterization with integrated GC. First, we describe the important parameters of the GC and evaluate their performance using Finite-Difference Time-Domain (FDTD) simulations. As we intend to utilize the wafer-scale RIBE process for device fabrication, we incorporate the fabrication constraints in our design. Next, we optimize the design to maximize the collection efficiency from color center emission through a microscope objective and discuss coupling dynamics in the triangular geometry. In the end, we demonstrate the fabrication of the triangular cross-section GC using the RIBE process.

II Triangular cross-section fishbone grating coupler
In this section, we analyze the design parameters for triangular cross-section suspended GC in SiC. Our triangular cross-section GC consists of (i) a triangular waveguide, (ii) an adiabatically tapered region, and (iii) a fishbone like grating structure. We choose the fishbone grating design in accordance with the fabrication compatibility of the RIBE process to produce wafer-scale triangular cross-section devices. Fig. 1(b) shows the top view schematic of the GC with waveguide width , etch angle (half angle at the apex) , taper angle , taper length , maximum taper width , grating period , and the central bone width . The refractive index of SiC is kept throughout the analysis. We parameterize the scales w.r.t. wavelength of interest to make the design scalable for the variety of color centers SiC hosts. We optimize each component individually to achieve the most efficient design.
II.1 Single-mode waveguide
Single-mode propagation is an important requirement in the waveguide design for performing on-chip QIP Caves and Drummond (1994); Humphreys et al. (2013). For each , there exists an optimal width for single-mode propagation in the triangular cross-section waveguides Babin et al. (2022). In this work, we choose as it offers higher confinement Saha, Majety, and Radulaski (2023) and coupling Babin et al. (2022); Majety et al. (2024a) compared to other etch angles that have been fabricated with the RIBE process in SiC Majety et al. (2024b). Based on the dipole orientation in the SiC crystal lattice, color center emission can couple to fundamental transverse electric (f-TE), fundamental transverse magnetic (f-TM), and other higher-order modes supported by the triangular waveguide. Hence, we optimize the triangular waveguide width for maintaining single-mode propagation of both horizontal and vertical electric dipole radiations with high coupling efficiency () using Lumerical FDTD package as shown in Fig. 1(c). We position the dipole sources with emission wavelength at the centroid of the triangle which is the optimum color center position in triangular geometry for the best coupling to the waveguide modesMajety et al. (2021); Babin et al. (2022); Majety et al. (2023). We pre-select the range using Lumerical MODE (Finite-Difference Eigensolver) package so that the horizontal (vertical) dipole couples only to the f-TE (f-TM) mode. Although the proposed design is a 1D-GC structure in which we can optimize only for a single state of polarization Marchetti et al. (2019), we choose to obtain high for both dipole orientations.

II.2 Adiabatic taper
An adiabatically tapered region is instrumental to coupling light in and out of the chip Son et al. (2018). It allows for a smooth transition of light from the waveguide mode to the Gaussian mode by gradually changing the size and shape of the optical mode along the taper. The slow transition of the effective refractive index () of the waveguide cross-section ensures minimal mode conversion to higher-order or radiation modes resulting in high coupling efficiency. Traditional fiber-coupled GCs require very long taper regions to match the mode field diameter of the single-mode optical fibers Cheng et al. (2020). On the other hand, light coupling through the microscope objective enables very compact designs for such tapers as the focused Gaussian beam spot size becomes orders of magnitude smaller compared to the optical fiber mode.
In the SOI platform, it has been observed that linear tapers with very small taper angle () result in higher transmission Fu et al. (2014). We apply a similar strategy in designing the triangular cross-section taper region by linearly increasing the taper width for , shown in Fig. 2(a). To ensure gradual change in of the waveguide mode, we calculate the of the f-TE and f-TM modes (mode #1 and mode #2, respectively) of the triangular waveguide along . Fig. 2(b) demonstrates that the adiabatic condition for triangular waveguides Burek et al. (2017): , where is the effective index of the waveguide cross-section with maximum taper width and is the effective index of the th waveguide mode along , is satisfied for all three . However, owing to the etch selectivity of the mask in the developed RIBE wafer-scale processMajety et al. (2024b), we restrict to be , larger than , for undercutting both and cross-section profiles in a single etch. Although all three designs exhibit unity transmission for the f-TE mode injection from the waveguide, we choose leading to for .

II.3 Fishbone grating
Refractive index variation with periodic configuration in the direction of light propagation is the most preferred approach for breaking the device symmetry to achieve out-of-plane coupling through diffractionCheng et al. (2020). Our fishbone grating is constructed by varying the width of the triangular cross-section periodically to implement refractive index contrast. The central bone ensures the mechanical stability of the suspended structure. The diffraction mechanism for a GC is governed by Bragg phase matching condition: sin = , where is the effective index of the guided mode in the grating medium, is the refractive index of air, is the diffraction angle, is the grating order, is the target wavelength, and is the grating period. For 1D-GC, the 1st order diffracted mode is vertically emitted from the grating Marchetti et al. (2019). Hence, with and , we can achieve vertical outcoupling and the phase matching condition reduces to . Due to the 1D periodicity, we can optimize for only one polarization state, in our case the f-TE mode. We estimate the using the analytical expression Zhang et al. (2015): , where is the duty cycle, or the fill factor, of the wider section, and and are the effective indices of the triangular cross-section with widths and , respectively. The grating period and duty cycle that satisfy the phase matching condition for different values are presented in Fig. 2(d). We observe that there exists no physical mode in the triangular cross-section for at which means the effective index of the narrow regions is . Therefore, we choose for obtaining maximum index contrast.
III Design Optimization
The GC design optimization is focused on enhancing the collection efficiency from color center integrated nanophotonic devices for characterizing light-matter interaction. Fig. 3(a) shows the side view of the 3D FDTD simulation setup for calculating the collection efficiency. We use a mode source to inject the f-TE mode of the waveguide (width ) at into the grating. A field and power monitor is placed above the triangular cross-section GC at a distance for measuring the fraction of power going in the upward direction. As we are interested in maximizing the light collection through a microscope objective, we add another field and power monitor very close to the top surface of the GC to calculate the far field radiation pattern and subsequently the fraction of power collected by an objective with a given NA (0.65 in this work). The collection efficiency is obtained by multiplying the power fractions calculated from these two monitors. We use a field and power monitor before the source to check the amount of light reflected back into the waveguide.

For a given design, the collection increases with the number of gratings and saturates for 5 unit cells. As a result, we use 5 grating periods for all the GCs in our analysis. Fig. 3(b) exhibits the collection efficiencies () from the injected f-TE waveguide mode with variation in and . We observe that reflection decreases with the increase in . The maximum is , acquired from a design with and . The obtained value is comparable to another compact GC design realized in silicon nitride platform for microscope objective collection where the theoretical coupling efficiency reaches up to ( with a reflective mirror coating on the substrate) Zhu et al. (2017). Notably, the most efficient design does not satisfy the phase matching condition which can stem from the broken symmetry of the triangular profile compared to rectangular structures. This phenomenon is observed in another theoretical work with triangular cross-section GCs in gallium nitride Hadden et al. (2022). Fig. 3(c)-(d) shows the far field radiation pattern and side view cross-section of electric field intensity () distribution, respectively, at for the most efficient design. We also simulate the collection efficiency of the optimized design as a function of normalized wavelength . The collection efficiency for the f-TM mode at is , lower than half of , which is expected. Conventional 1D grating couplers have to deal with bandwidth issues Andreani et al. (2016); Vitali et al. (2023) whereas the triangular cross-section GC broadband response appears flat, meaning the same design can be utilized for studying different zero-phonon line (ZPL) emissions from a color center in SiC. To characterize the incoupling efficiency of the optimized design, we employ a Gaussian source focused from above by an objective with 0.65 NA and a power monitor in the waveguide region as shown in Fig. 4(b). For higher incoupling, the source needs to be focused on the first grating period. Focusing on other unit cells results in lower incoupling as a large amount of light readily couples to the guided mode of the adiabatic taper as seen in Fig. 4(d). The s-polarized Gaussian beam mostly couples to the f-TE mode while the p-polarized beam couples mostly to the f-TM mode. Although we observe some initial coupling to high-order modes, they diminish as the incoupled light propagates through the waveguide. Fig. 4(c) illustrates the broadband response for the incoupling efficiency of the optimized design. The incoupling efficiency is much lower than outcoupling as our optimization technique prioritizes the collection efficiency.

Non-uniform GCs have proven to be very effective in enhancing the overall efficiency Marchetti et al. (2019); Ding, Ou, and Peucheret (2013); Lee et al. (2016). Non-uniformity is generally introduced by linearly apodizing the or along the direction of propagation. In search of the most optimized design, we opt for a similar technique. Keeping constant, we linearly change the of each unit cell with the following equation: , where is the apodization factor, , and . However, we observe the slowly decreasing for both cases with variation in from 0 to 0.02. On the other hand, we linearly change for each grating section, with constant , by the following equation: , where is the aperiodic factor, , and . Fig. 5(a) demonstrates that gradually increases (decreases) with increasing (decreasing) and starts to decrease (increase) after reaching an inversion point. We notice that even for the most optimized aperiodic design, with , only increases by . However, we observe a significant change in the broadband response, Fig. 5(b), and the far field radiation pattern, Fig. 5(d), of the aperiodic GC. The far field radiation appears less dispersive compared to the periodic design. While collection and incoupling efficiencies for the f-TM mode remain invariant for the aperiodic optimization, the broadband response for the f-TE mode becomes a normal distribution. Based on the specific application requirement, one can adopt either the periodic or the aperiodic design as both designs exhibit comparable collection efficiencies.
IV Proof-of-Principle Fabrication and Discussion

Designed grating couplers can readily be implemented with the wafer-scale quantum nanophotonic fabrication processes. As illustrated in Fig. 1(a), we envision fabricating a suspended triangular waveguide with grating couplers positioned on both ends. One coupler is used to direct the laser into the waveguide and stimulate the color center with off-resonant excitation. The other one couple both the excitation light and the color center emission out into a microscopy system where the excitation beam is selectively filtered out. With this mechanism, light-matter interaction in the triangular geometry can be experimentally characterized for all the suspended active and passive photonic devices.
We test the key step in the fabrication of the studied GC designs in SiC with the RIBE process. In this method, a collimated and uniform beam of high energy reactive and inert gas ions are directed toward a substrate holder that can be rotated and tilted at an angle. The tilt of the substrate holder defines the of the triangular cross-section and the rotation ensures a uniform device undercut. The presence of physical etch, due to the inert gas ions, reduces the mask etch selectivity as well as causes higher re-deposition from the sputtered material which makes the etching process slower. Hence, we have considered these constraints in our designs as explained in Section II.
We choose nm, one of the ZPL emissions from NV center in 4H-SiC, for both the periodic and aperiodic GC fabrication. Fig. 6(a) shows the nanofabrication process flow. We use electron beam lithography for defining the GC patterns on 4H-SiC substrate. Then, we deposit 120 nm thick nickel as a hard metal mask using electron beam evaporation and transfer the GC patterns on to the mask through a lift-off process. The triangular cross-section waveguides and GCs are then fabricated by the RIBE process using SF6, O2, and Ar gas ions. During the etching process, we deliberately set the substrate tilt angle to for obtaining Majety et al. (2024b). The SEM images of the fabricated GCs are presented in Fig. 6(b)-(e). We observe that the undercut in these devices are nm, after an hour of etching, which can be a bottleneck for experimental characterization as a significant amount of light can couple to the bulk. This issue can be resolved by performing a prolonged etch. However, a practical robust solution can be the inclusion of inductively coupled plasma-reactive ion etching (ICP-RIE) in the fabrication process flow. The etch rate for ICP-RIE, with the same SF6 and O2 chemistry, is much higher compared to the RIBE process. As a result, performing a vertical etch with ICP-RIE, before the angle-etching step, can result in a deeper undercut essential for experimental characterization of the fabricated devices.
Triangular cross-section photonic devices can be fabricated at a wafer-scale in SiC with the RIBE process. This process preserves color center optical properties which facilitates the implementation of on-chip QIP with integrated quantum photonic mesh architecture. With the addition of the proposed GCs on the same chip, wafer-level testing can be executed on the fabricated devices. Our simulation results show that the triangular cross-section fishbone grating design can yield high collection efficiency () for studying the color center coupling with the integrated nanophotonic devices and initial fabrication results demonstrate the feasibility of producing the optimized designs with the RIBE process.
Acknowledgements.
We acknowledge support from NSF CAREER (Award 2047564) and AFOSR Young Investigator Program (Award FA9550-23-1-0266). Work at the Molecular Foundry was supported by the Office of Science, Office of Basic Energy Sciences, of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. Part of this study was carried out at the UC Davis Center for NanoMicro Manufacturing (CNM2). The authors would like to acknowledge the valuable support provided by Elnaz Hamdarsi and CNM2 staff member Vishal Narang.Author declaration
Conflict of Interest
The authors have no conflict of interest.
Author Contributions
PS modeled, analyzed and fabricated devices, SM developed the angle etching process, SD performed electron beam lithography, MR supervised the project.
Data Availability Statement
The data that support the findings of this study are available upon reasonable request from the corresponding author Pranta Saha ([email protected]).
References
References
- Kimble (2008) H. J. Kimble, Nature 453, 1023 (2008).
- Monroe (2002) C. Monroe, Nature 416, 238 (2002).
- Azuma et al. (2023) K. Azuma, S. E. Economou, D. Elkouss, P. Hilaire, L. Jiang, H.-K. Lo, and I. Tzitrin, Reviews of Modern Physics 95, 045006 (2023).
- Norman et al. (2021) V. A. Norman, S. Majety, Z. Wang, W. H. Casey, N. Curro, and M. Radulaski, InfoMat 3, 869 (2021).
- Hensen et al. (2015) B. Hensen, H. Bernien, A. E. Dréau, A. Reiserer, N. Kalb, M. S. Blok, J. Ruitenberg, R. F. Vermeulen, R. N. Schouten, C. Abellán, et al., Nature 526, 682 (2015).
- Majety et al. (2022) S. Majety, P. Saha, V. A. Norman, and M. Radulaski, Journal of Applied Physics 131 (2022).
- Anderson et al. (2022) C. P. Anderson, E. O. Glen, C. Zeledon, A. Bourassa, Y. Jin, Y. Zhu, C. Vorwerk, A. L. Crook, H. Abe, J. Ul-Hassan, et al., Science advances 8, eabm5912 (2022).
- Cilibrizzi et al. (2023) P. Cilibrizzi, M. J. Arshad, B. Tissot, N. T. Son, I. G. Ivanov, T. Astner, P. Koller, M. Ghezellou, J. Ul-Hassan, D. White, et al., Nature Communications 14, 8448 (2023).
- Castelletto et al. (2022) S. Castelletto, A. Peruzzo, C. Bonato, B. C. Johnson, M. Radulaski, H. Ou, F. Kaiser, and J. Wrachtrup, ACS Photonics 9, 1434 (2022).
- Lukin, Guidry, and Vučković (2020) D. M. Lukin, M. A. Guidry, and J. Vučković, PRX quantum 1, 020102 (2020).
- Udvarhelyi et al. (2019) P. Udvarhelyi, R. Nagy, F. Kaiser, S.-Y. Lee, J. Wrachtrup, and A. Gali, Physical Review Applied 11, 044022 (2019).
- Majety et al. (2021) S. Majety, V. A. Norman, L. Li, M. Bell, P. Saha, and M. Radulaski, Journal of Physics: Photonics 3, 034008 (2021).
- Babin et al. (2022) C. Babin, R. Stöhr, N. Morioka, T. Linkewitz, T. Steidl, R. Wörnle, D. Liu, E. Hesselmeier, V. Vorobyov, A. Denisenko, et al., Nature materials 21, 67 (2022).
- Majety et al. (2023) S. Majety, S. Strohauer, P. Saha, F. Wietschorke, J. J. Finley, K. Müller, and M. Radulaski, Materials for Quantum Technology 3, 015004 (2023).
- Saha, Majety, and Radulaski (2023) P. Saha, S. Majety, and M. Radulaski, Scientific reports 13, 4112 (2023).
- Majety et al. (2024a) S. Majety, P. Saha, Z. Kekula, S. Dhuey, and M. Radulaski, MRS communications , 1 (2024a).
- Majety et al. (2024b) S. Majety, V. A. Norman, P. Saha, A. H. Rubin, S. Dhuey, and M. Radulaski, arXiv preprint arXiv:2405.07498 (2024b).
- Cheng et al. (2020) L. Cheng, S. Mao, Z. Li, Y. Han, and H. Fu, Micromachines 11, 666 (2020).
- Marchetti et al. (2019) R. Marchetti, C. Lacava, L. Carroll, K. Gradkowski, and P. Minzioni, Photonics Research 7, 201 (2019).
- Norman et al. (2024) V. A. Norman, S. Majety, A. H. Rubin, P. Saha, J. Simo, B. Palomarez, L. Li, P. B. Curro, S. Dhuey, S. Virasawmy, et al., arXiv preprint arXiv:2401.10509 (2024).
- Faraon et al. (2008) A. Faraon, I. Fushman, D. Englund, N. Stoltz, P. Petroff, and J. Vučković, Optics express 16, 12154 (2008).
- Zhou et al. (2018) X. Zhou, I. Kulkova, T. Lund-Hansen, S. L. Hansen, P. Lodahl, and L. Midolo, Applied Physics Letters 113 (2018).
- Zhu et al. (2017) Y. Zhu, J. Wang, W. Xie, B. Tian, Y. Li, E. Brainis, Y. Jiao, and D. Van Thourhout, Optics express 25, 33297 (2017).
- Caves and Drummond (1994) C. M. Caves and P. D. Drummond, Reviews of Modern Physics 66, 481 (1994).
- Humphreys et al. (2013) P. C. Humphreys, B. J. Metcalf, J. B. Spring, M. Moore, X.-M. Jin, M. Barbieri, W. S. Kolthammer, and I. A. Walmsley, Physical review letters 111, 150501 (2013).
- Son et al. (2018) G. Son, S. Han, J. Park, K. Kwon, and K. Yu, Nanophotonics 7, 1845 (2018).
- Fu et al. (2014) Y. Fu, T. Ye, W. Tang, and T. Chu, Photonics Research 2, A41 (2014).
- Burek et al. (2017) M. J. Burek, C. Meuwly, R. E. Evans, M. K. Bhaskar, A. Sipahigil, S. Meesala, B. Machielse, D. D. Sukachev, C. T. Nguyen, J. L. Pacheco, et al., Physical Review Applied 8, 024026 (2017).
- Zhang et al. (2015) J. Zhang, J. Yang, H. Lu, W. Wu, J. Huang, and S. Chang, Infrared Physics & Technology 71, 542 (2015).
- Hadden et al. (2022) J. Hadden, C. Maynard, D. M. Beggs, R. A. Taylor, and A. J. Bennett, Journal of Physics D: Applied Physics 55, 474002 (2022).
- Andreani et al. (2016) L. C. Andreani, D. Gerace, M. Passoni, A. Bozzola, and L. Carroll, in 2016 18th International Conference on Transparent Optical Networks (ICTON) (IEEE, 2016) pp. 1–4.
- Vitali et al. (2023) V. Vitali, T. Domínguez Bucio, C. Lacava, R. Marchetti, L. Mastronardi, T. Rutirawut, G. Churchill, J. Faneca, J. C. Gates, F. Gardes, et al., Photonics Research 11, 1275 (2023).
- Ding, Ou, and Peucheret (2013) Y. Ding, H. Ou, and C. Peucheret, Optics letters 38, 2732 (2013).
- Lee et al. (2016) M. H. Lee, J. Y. Jo, D. W. Kim, Y. Kim, and K. H. Kim, Journal of the Optical Society of Korea 20, 291 (2016).