Note: the cost of running the entire notebook is around 20 FlexCredits.

Photonic integrated circuits that operate at visible wavelengths (400-700nm) are gaining interest for emerging applications beyond traditional telecommunications. Compared to PICs designed for infrared wavelengths, visible light PICs can potentially enable new applications in spectroscopy, bio-medical imaging, quantum computing, augmented reality, and more. However, efficiently interfacing visible light between optical fibers and nanophotonic components on a chip remains challenging. Due to the shorter wavelengths, photonic waveguides must be smaller to confine light, leading to significant mode mismatch with optical fibers. Efficient coupling techniques like inverse tapers reach the limits of fabrication capabilities at visible wavelengths. Novel solutions are needed to address these challenges in fiber-to-chip coupling for visible light PICs.

In this notebook, we demonstrate a SiN bi-layer edge coupler design for efficient and broadband fiber-to-chip coupling in visible frequency. The design uses a thin top SiN layer for mode field diameter expansion at the facet to improve fiber coupling and a thicker bottom SiN layer for tight mode confinement and routing. The design achieves a coupling loss ≤4 dB over a broad wavelength range, which is 3-5 dB lower compared to conventional single-layer edge couplers. The example is adapted from Yiding Lin, Jason C. C. Mak, Hong Chen, Xin Mu, Andrei Stalmashonak, Youngho Jung, Xianshu Luo, Patrick G.-Q. Lo, Wesley D. Sacher, and Joyce K. S. Poon, "Low-loss broadband bi-layer edge couplers for visible light," Opt. Express 29, 34565-34576 (2021) DOI:10.1364/OE.435669.

In another case study, we investigate inverse taper edge couplers at the telecom frequency. In addition, different grating coupler designs including the uniform grating coupler, the focusing apodized grating coupler, and the inverse designed compact grating coupler are also explored in various case studies.

Simulation Setup¶

The wavelength range of interest for this edge coupler is 400 nm to 640 nm.

The bi-layer structure is surrounded by a 2 $\mu m$ top SiO$_2$ cladding and a 3$\mu m$ bottom SiO$_2$ layer. We do not model the air region above the top cladding and the silicon region below the SiO$_2$ layer since the mode diameter of the single mode fiber is only about 3.7 $\mu m$. In this case, we first define two materials: SiN and SiO$_2$, both of which are described by a constant refractive index given in the reference. Using non-dispersive media also helps to reduce the computational cost.

Based on the reference, the optical fiber used is a Nufern S405-XP single-mode fiber, which has a core diameter of 3 $\mu m$ and a numerical aperture (NA) of 0.12. Here we assume The cladding material is the same SiO$_2$ and calculate the core refractive index based on the NA and $n_{SiO_2}$.

Next we define the geometric parameters of the edge coupler. The physical meanings of the parameters are schematically shown in Fig.1 of the publication. We assume a 500 nm air gap between the optical fiber facet to the beginning of the edge coupler.

Given the shapes of the SiN layers, the most convenient way to define them is by using PolySlab. The vertices can be easily calculated from the geometric parameters defined above.

Since we already know the core radius and the refractive indices of the fiber core and cladding materials, we model the fiber explicitly and then use a ModeSource to launch the fiber mode. Alternatively, we can use a GaussianBeam with the same beam size to mimic the fiber mode, as demonstrated in the previous inverse taper edge couplers example. We will run two separate simulations to investigate the coupling efficiencies of TE and TM modes. The two polarizations correspond to two different symmetry conditions.

To compute the coupling efficiency, we add a ModeMonitor at the straight SiN waveguide. Although not necessary, it's still quite helpful to add a FieldMonitor to visualize the field distribution. Since the simulation domain is quite large, we don't need to store the field at each grid point to have a good field distribution plot. Therefore, we set interval_space=(3,1,1) to downsample the stored field by a factor of 3, which helps reduce the size of the downloaded data significantly.

First we define the simulation with TE polarization. This polarization corresponds to symmetry=(0,-1,0). To learn how to properly utilize symmetry to reduce computational cost, please refer to our symmetry tutorial. The simulation domain sizes in the $y$ and $z$ directions need to be sufficiently large to ensure the optical fiber mode as well as the waveguide mode field at the tip of the edge coupler decays to zero at the boundaries.

In addition, we zoom in to the optical fiber region to ensure it is set up correctly.

The TM polarization simulation can be defined by copying the TE simulation and updating the symmetry.

Each simulation contains over 1 billion grid points since the length of the coupler is over 500 free space wavelengths. For Tidy3D, such large simulations can be handled easily.

Batch Sweeping Different Polarizations¶

We can run the two simulations sequentially. However, for maximum efficiency, we will run them concurrently by making a simulation Batch. To do so, we first put simulations in a Python dictionary.

Then we submit the simulation as a Batch to run on the server.

Result Visualization¶

After the simulations are complete, we can calculate and plot the coupling efficiencies. For both polarizations, the coupling losses are below 4 dB in a broad wavelength range, which is a great improvement compared to the simple single-layer edge coupler.

Lastly, we can plot the field norm distribution to visualize the coupling process.