Simulation of Light Coupling Reciprocity for a Photonic ...1.5 1.55 1.6 1.65 0 0.15 0.3 0.45 Uniform...

1
1.5 1.55 1.6 1.65 0 0.15 0.3 0.45 Uniform grating, coupled power Etch width distribution, coupled power Uniform grating, up scattered power Etch width distribution, up scattered power λ ( μm) P out / P in Simulation of Light Coupling Reciprocity for a Photonic Grating V. Kivijärvi 1 , M. Erdmanis 1 , I. Tittonen 1 1. Aalto University, Department of Micro- and Nanosciences, FI-00076 Aalto, Finland single mode fiber SOI waveguide port photonic grating absorbing layer 11.8 ° ρ -5.2 -2.6 0 2.6 5.2 0 0.25 0.5 0.75 1 position(μm) Uniform grating Etch width distribution Si = 3.45 SiO 2 = 1.44 air =1 Wave equation: Absorbing layer index: × × E - E =0 = + - / Gaussian beam electric field: Phase matching condition: E = z / ϕ = sin + /Λ Improvement in coupling E z / E zmax Introduction: SOI (Silicon on Insulator) technology utilizes silicon components on SiO 2 layer. Silicon waveguides carry signals on a SOI microchip. Propagating electric field distribution in a waveguide is called mode of the waveguide. Photonic gratings (Fig. 1) are formed by etching grooves on the top of a waveguide. Gratings can operate in two directions. They can guide incident beam into a waveguide or a waveguide mode out of the structure. Aim of the study: We study the grating operation by utilizing optical reciprocity (1), which manifests that coupling of light between a waveguide and an optical fiber should not depend on the propagation direction. This principle allows for study of operation in two directions and determination of the coupled and lost power. Acquired knowledge can be used to optimize the grating geometry. Computational Methods: We calculate the propagation of signal in two directions and compare the electric field profiles. The observed difference for the two cases allows for estimation of the mode mismatch between the coupled and scattered modes. Based on this comparison, we are able to calculate both the losses and coupled power. The losses originate from the mode mismatch between the signal scattered by grating and a fundamental fiber mode. We show that they can be minimized by implementing the proper profile into the grating. The simulations were performed by FD and BM analysis of COMSOL RF module. We used two excitation schemes: one with the power entering through the waveguide and another with an incident beam irradiating the grating (the beam was supposed to arrive from an optical fiber). In the fiber excitation scheme, the incident light was launched from the white rectangle, as shown in Fig. 4. The electric field profile was set to be Gaussian along the rectangle front side inducing a Gaussian beam (Fig. 4.) through the grating center. The coupled and scattered powers were calculated through the power flow integration over the model internal boundaries. The undesirable reflections from the model boundaries were suppressed by self constructed absorbing layer with the linearly increasing imaginary part of refractive index. The mesh density was adjusted to be dense near the grating and sparse in distant areas. In the absorbing layer, numerical reflections were minimized by using a mapped mesh. Matlab LiveLink TM was used to optimize the grating geometry in a semi-automated fashion. The code performed flipping of the geometry, stored electric field distributions, specified by unique ID string. Figure 2a. Uniform grang. Figure 2b. Opmized geometry. Figure 3. Signicant parameters. Figure 4. The model used in the simulaons. Results: By tracing the electric field profiles for the two directions of excitation and the corresponding power values, we identified the optimized grating geometry (Fig. 2b). In this case, the coupled power reached 42% at 1.55um, while before the optimization it was only 35% (Fig. 5). Similarly, the scattered field profiles (Fig. 7) verified the improvement in coupling efficiency was due to the reduced mismatch between the signal provided by the grating and the fundamental fiber mode. The scattered powers were almost equal for both geometries, as shown in Fig. 5. This supports the idea the improvements are caused by the mode matching enhancement but not increased overall scattered power. Governing equations: Figure 6. Coupling from the fiber (electric field z component). Figure 7. Scattered electric field profiles (normalized z component). Figure 8. Coupling from the waveguide (electric field z-component). Conclusions: The model represented in this study can be used to determine coupled power and coupling losses for various grating geometries. In addition, it is compatible with LiveLink TM based automatic calculation which allows for optimization of grating geometry. Following observations can be listed: The optimized geometry reaches higher coupling with smaller losses. The optimized geometry produces a focused beam of light in the fiber direction. The optimized geometry produces electric field distribution which resembles fundamental fiber mode. Significance of the work: Improved coupling can benefit optical technology because it allows for better communication between a SOI microchip and a network composed of optical fibers. Development of optical technology has became important because traditional electronics is fundamentally speed limited by narrow bandwidths and energy consumption. References: 1. R. J. Potton, Reciprocity in Optics, Reports on Progress in Physics, 67, pp. 719-749 (2004). 2. D. Taillaert, P. Bienstman et al., Compact efficient broadband grating coupler for silicon-on-insulator waveguides, Optics Letters, 29, pp. 75-86, (2007). 3. Tillack, B. (2013). Technology Solutions. Retrieved from IHP - Innovations for High Performance Microelectronics: http://www.ihp.microelectronics.com/en/solutions/technology-solutions.html. Acknowledgements: We acknowledge financial support from the Academy of Finland, project number 140009. +|E zmax | -|E | 0V/m zmax +|E zmax | -|E | 0V/m zmax Figure 1. SEM image of a photonic grating (3). Figure 5. Frequency responses. (similar curves for the uniform case are calculated in (2)).

Transcript of Simulation of Light Coupling Reciprocity for a Photonic ...1.5 1.55 1.6 1.65 0 0.15 0.3 0.45 Uniform...

  • 1.5 1.55 1.6 1.650

    0.15

    0.3

    0.45

    Uniform grating, coupled powerEtch width distribution, coupled powerUniform grating, up scattered powerEtch width distribution, up scattered power

    λ(μm)

    P out

    / P in

    E z /

    P in

    Simulation of Light Coupling Reciprocity for a Photonic GratingV. Kivijärvi1, M. Erdmanis1, I. Tittonen1

    1. Aalto University, Department of Micro- and Nanosciences, FI-00076 Aalto, Finland

    single mode �ber

    SOI waveguide port photonic grating

    absorbing layer

    11.8 °

    ρ

    -5.2 -2.6 0 2.6 5.20

    0.25

    0.5

    0.75

    1

    position(μm)

    Uniform gratingEtch width distribution

    Si = 3.45

    SiO 2 = 1.44

    air = 1

    Wave equation: Absorbing layer index:∇× ∇× E − E = 0 = + − /∆

    Gaussian beam electric �eld: Phase matching condition:

    E = z / ϕ = sin + /Λ

    Improvement in couplingE z

    / E z

    max

    Introduction: SOI (Silicon on Insulator) technology utilizes siliconcomponents on SiO2 layer. Silicon waveguides carry signals on a SOImicrochip. Propagating electric field distribution in a waveguide is calledmode of the waveguide. Photonic gratings (Fig. 1) are formed by etchinggrooves on the top of a waveguide. Gratings can operate in twodirections. They can guide incident beam into a waveguide or awaveguide mode out of the structure.

    Aim of the study: We study the grating operation by utilizing opticalreciprocity (1), which manifests that coupling of light between awaveguide and an optical fiber should not depend on the propagationdirection. This principle allows for study of operation in two directions anddetermination of the coupled and lost power. Acquired knowledge can beused to optimize the grating geometry.

    Computational Methods: We calculate the propagation of signal in two directions andcompare the electric field profiles. The observed difference for the two cases allows forestimation of the mode mismatch between the coupled and scattered modes. Based on thiscomparison, we are able to calculate both the losses and coupled power. The losses originatefrom the mode mismatch between the signal scattered by grating and a fundamental fiber mode.We show that they can be minimized by implementing the proper profile into the grating.

    The simulations were performed by FD and BM analysis of COMSOL RF module. We used twoexcitation schemes: one with the power entering through the waveguide and another with anincident beam irradiating the grating (the beam was supposed to arrive from an optical fiber). Inthe fiber excitation scheme, the incident light was launched from the white rectangle, as shownin Fig. 4. The electric field profile was set to be Gaussian along the rectangle front side inducinga Gaussian beam (Fig. 4.) through the grating center. The coupled and scattered powers werecalculated through the power flow integration over the model internal boundaries. Theundesirable reflections from the model boundaries were suppressed by self constructedabsorbing layer with the linearly increasing imaginary part of refractive index. The mesh densitywas adjusted to be dense near the grating and sparse in distant areas. In the absorbing layer,numerical reflections were minimized by using a mapped mesh.

    Matlab LiveLinkTM was used to optimize the grating geometry in a semi-automated fashion. Thecode performed flipping of the geometry, stored electric field distributions, specified by unique IDstring.

    Figure 2a. Uniform grating.

    Figure 2b. Optimized geometry. Figure 3. Significant parameters.

    Figure 4. The model used in the simulations.

    Results: By tracing the electric field profiles for the two directions of excitation and thecorresponding power values, we identified the optimized grating geometry (Fig. 2b). In thiscase, the coupled power reached 42% at 1.55um, while before the optimization it was only35% (Fig. 5). Similarly, the scattered field profiles (Fig. 7) verified the improvement in couplingefficiency was due to the reduced mismatch between the signal provided by the grating andthe fundamental fiber mode. The scattered powers were almost equal for both geometries, asshown in Fig. 5. This supports the idea the improvements are caused by the mode matchingenhancement but not increased overall scattered power.

    Governing equations:

    Figure 6. Coupling from the fiber(electric field z component).

    Figure 7. Scattered electric fieldprofiles (normalized z component).

    Figure 8. Coupling from the waveguide(electric field z-component).

    Conclusions: The model represented in this study can be used to determine coupled power and coupling losses for various grating geometries. In addition, it is compatible with LiveLinkTM based automatic calculation which allows for optimization of grating geometry. Following observations can be listed:• The optimized geometry reaches higher coupling with smaller losses.• The optimized geometry produces a focused beam of light in the fiber direction.• The optimized geometry produces electric field distribution which resembles fundamental

    fiber mode.

    Significance of the work: Improved coupling can benefit optical technology because it allows for better communication between a SOI microchip and a network composed of optical fibers. Development of optical technology has became important because traditional electronics is fundamentally speed limited by narrow bandwidths and energy consumption.

    References:1. R. J. Potton, Reciprocity in Optics, Reports on Progress in Physics, 67, pp. 719-749

    (2004).2. D. Taillaert, P. Bienstman et al., Compact efficient broadband grating coupler for

    silicon-on-insulator waveguides, Optics Letters, 29, pp. 75-86, (2007).3. Tillack, B. (2013). Technology Solutions. Retrieved from IHP - Innovations for High

    Performance Microelectronics: http://www.ihp.microelectronics.com/en/solutions/technology-solutions.html.

    Acknowledgements:We acknowledge financial support from the Academy of Finland, project number 140009.

    +|Ezmax|

    -|E |0V/m

    zmax

    +|Ezmax|

    -|E |0V/m

    zmax

    Figure 1.SEM image of aphotonic grating (3).

    Figure 5. Frequency responses.(similar curves for the uniformcase are calculated in (2)).