• Photonics Research
  • Vol. 8, Issue 3, 381 (2020)
Roy Avrahamy1, Moshe Zohar2, Mark Auslender1、*, Benny Milgrom3, Shlomo Hava1、2, and Rafi Shikler1、4
Author Affiliations
  • 1School of Electrical and Computer Engineering, Ben-Gurion University of the Negev, Beer-Sheva 84105, Israel
  • 2Electrical and Electronics Engineering Department, Shamoon College of Engineering, Beer-Sheva 84100, Israel
  • 3School of Electrical Engineering, The Jerusalem College of Technology, Jerusalem 91160, Israel
  • 4Ilse Katz Institute for Nanoscale Science & Technology, Ben-Gurion University of the Negev, Beer-Sheva 84105, Israel
  • show less
    DOI: 10.1364/PRJ.375802 Cite this Article Set citation alerts
    Roy Avrahamy, Moshe Zohar, Mark Auslender, Benny Milgrom, Shlomo Hava, Rafi Shikler, "In-depth investigation and applications of novel silicon photonics microstructures supporting optical vorticity and waveguiding for ultra-narrowband near-infrared perfect absorption," Photonics Res. 8, 381 (2020) Copy Citation Text show less

    Abstract

    We propose a novel concept of designing silicon photonics metamaterials for perfect near-infrared light absorption. The study’s emphasis is an in-depth investigation of various physical mechanisms behind the 100% ultra-narrowband record peak absorptance of the designed structures, comprising an ultrathin silicon absorber. The electromagnetic power transport, described by the Poynting vector, is innovatively explored, which shows combined vortex and crossed-junction two-dimensional waveguide-like flows as outcomes of optical field singularities. These flows, though peculiar for each of the designed structures, turn out to be key factors of the perfect resonant optical absorption. The electromagnetic fields show tight two-dimensional confinement: a sharp vertical confinement of the resonant-cavity type combined with a lateral metasurface supported confinement. The silicon-absorbing layer and its oxide environment are confined between two subwavelength metasurfaces such that the entire design is well compatible with silicon-on-insulator microelectronics. The design concept and its outcomes meet the extensive challenges of ultrathin absorbers for minimum noise and an ultra-narrowband absorptance spectrum, while maintaining an overall very thin structure for planar integration. With these materials and such objectives, the proposed designs seem essential, as standard approaches fail, mainly due to a very low silicon absorption coefficient over the near-infrared range. Tolerance tests for fabrication errors show fair tolerability while maintaining a high absorptance peak, along with a controllable deviation off the central-design wavelength. Various applications are suggested and analyzed, which include but are not limited to: efficient photodetectors for focal plane array and on-chip integrated silicon photonics, high-precision spectroscopic chemical and angular-position sensing, and wavelength-division multiplexing.

    1. INTRODUCTION

    Enhancement of optical absorption is sorely sought in various light-management devices, such as solar photovoltaic and heat-storage cells and photodetectors (PDs), and it may also contribute to cutting-edge optical camouflage implementations. Silicon (Si) has such physical properties, including optical and electronic ones, that made it dominant in microelectronics and also a favorite for many photonic and optoelectronic applications, in particular those mentioned above. Broadly, Si photonics terms the field of prediction, design, fabrication, and implementation of photonic devices made of Si, its compounds with O and N, and with other IV-group elements. Often, this term is used in a narrower sense to refer to optical solutions to problems of connectivity and communication within a Si chip package by integrating Si-based optics and optoelectronics within the chip. Research on Si photonics began in the 1990s [15], proceeded through the 2000s [69], and has continued ever since [1014]. On-chip optical interconnects require high-power lasers to convert the electrical data signals into light pulses, optical waveguides/fibers to propagate the light, and fast PDs to turn the optical signals back into electrical signals. Many breakthroughs have already been achieved [1014] in the realm of diverse Si-based photonic and optoelectronic components.

    From now on, the term PD is used for brevity to present any photo-absorbing device. Assume for simplicity that each photon absorbed by the photoactive layer of the PD generates an electron–hole pair. Then η=A, where η and A are the quantum efficiency and optical absorptance, respectively, of the PD. In a monolithic PD, the detected light traverses the absorber once, so η may be increased only by increasing its thickness ta, but the thicker the absorber, the stronger the thermal noise and the longer the transit time (in a diode PD). For some applications, e.g., PDs’ focal plane arrays (FPAs) [15], high η and low noise are of paramount importance in order to show high responsivity while avoiding faulty high-density packaging. For others, e.g., the on-chip integrated Si-PDs (which should also be thin enough to be compatible with current planar-chip technology), the short transit time is also crucial.

    A drastic decrease of ta to decrease the noise and transit time while maintaining high η requires the confinement of incoming light around the absorber layer. An approach adopted to this end is based on the solid-state resonant Fabry–Perot cavity, in which this layer is sandwiched between two multilayered sections. In each of these, the distributed Bragg reflector (DBR) acts as a mirror and two layers adjacent to the absorber layer form an optical cavity. There, even very thin active layer may absorb most of the incident light power. Such a resonant-cavity enhanced PD (RCE-PD) [16] enables much more efficient (and faster, when relevant) detection, as well as operation at higher temperatures with a much thinner absorber than in its monolithic counterpart, although these benefits are traded off against a smaller spectral bandwidth.

    The field of RCE-PDs has matured over the past decades regarding different spectral ranges, absorber materials, and physical devices; see, e.g., Refs. [36,17,18] for Si-PDs. However, the restriction on the DBRs’ reflectances [16] for attaining maximal peak η may require deposition of too many layers, such that even with quite a thin absorber the entire device may prove too thick for on-chip integration. Such a high number of layers, as a rule, exceeds a technology-defined critical value, beyond which growth defects inevitably emerge and impair the device’s performance even with lossless DBRs; see, e.g., Ref. [19] for the SiO2/Si3N4 ones. The reported on-chip Si-PDs, fabricated with thin but absorbing SiO2/Si DBRs and rather large ta [36], showed very low NIR efficiency ηmax<50%. Thus, reducing the Si RCE-PD thickness below some threshold clearly conflicts with the perfect-η objective.

    Due to the reasons outlined above, to concurrently meet the challenges of a high η, thin Si layer, and small overall thickness of a NIR Si-PD, a new design solution, quite different from the RCE approach, is needed. We suggest using metasurfaces [20], which enable us to redesign the known photonic devices to make them thinner and planar while concurrently being of a highly reduced complexity and a high performance. For our purposes, we use the simplest kind of metasurface, a subwavelength (SW) grating-on-layer (GL) metasurface. So far, an upgrade of the RCE-PD with the replacement of the front DBR by a one-dimensional (1D) GL (not necessarily an SW one) was proposed, and proof-of-concept examples of the thus modified HgCdTe-PDs were reported [21,22]. The round-trip phase control [16] extended to account for the presence of the grating and an optimization allowed for [21,22] fabrication-tolerant designs with peak η100% and a back DBR twice as thin as that of the best performing RCE-PD composed of the same materials. However, due to the remaining DBR, such a design solution proves yet insufficient for downscaling the Si-PD at goal.

    In Si photonics, a 0,±1 diffraction order Si-SiO2 grating was theorized [1] to bend normally incident light into the depletion region of the on-chip Si-PD. Another efficient light coupling into a CMOS-compatible Si-PD was demonstrated; see, e.g., Ref. [7], with a Si–air grating under oblique incidence, at which only one nonspecular diffraction order is supported. A similar strategy is followed [13,14,23] in contemporary Si photonics. In several important applications, the high-refractive-index (RI) contrast SW-GL, despite its structural simplicity, may function like a much more complex metasurface [20]. Due to this fact, such SW-GL nowadays attracts much interest [24] as a new platform for integrated optoelectronics, in particular, Si photonics [13,14].

    In this study, novel ultrathin Si-photonics structures are proposed and designed to show strikingly efficient light absorption within a narrow portion of the NIR range, around the center-design wavelength (CDW) λ0=0.8  μm. In theoretical studies [1], this wavelength was allocated as one from operation bandwidths of on-chip NIR Si photodiode PDs [35,7,9], implying that a short-haul light guiding with Si [1013] is still feasible at λ=λ0. Remarkably, these structures exhibit strong vortex and waveguide flows of the internal electromagnetic (EM) power as outcomes of the EM fields’ singularities [25].

    The rest of the paper is organized into six sections. Section 2 states the problem, outlines the suggested structures, and discusses their modeling. In Section 3, an optimization-based design procedure is described, and the considered structural parameters and computed absorptance spectra are presented. In Section 4, the simulations of the EM-fields-related physical quantities at the CDW are reported, and their results are discussed. Section 5 reports the designs’ tolerance testing results. Section 6 considers possible applications of the designed structures, and Section 7 concludes the study.

    2. STATEMENT OF PROBLEM, PROPOSED ENHANCED ABSORPTION STRUCTURES, AND PRINCIPLES OF MODELING

    In the quest to thin the Si-photonics structure as much as possible, it will perfectly absorb light at the designated NIR CDW while surpassing the bounds on the mirrors’ and cavity’s thicknesses noted in the introduction. We redesign the RCE structure by replacing both DBR mirrors by 1D GL with the same grating period Λ. Since SiO2 and Si3N4 are well technologically compatible with each other and Si for composing the DBR and the cavity [19], respectively, we employ these materials in the redesign, while benefiting from the previous design technique [21,22] to attain nearly perfect peak absorptance.

    Figure 1 shows a schematic of the dual-GL structures under study. There, the cavity and substrate material is SiO2, and the front air–Si3N4/Si3N4 with the back SiO2-Si3N4/Si3N4 GL encloses the cavity. As shown in Fig. 1, we consider the through-substrate backside illumination (BSI) at normal incidence. Benefiting from the absence of the backside electrical contacts, BSI has the essential advantage of reducing the PD pixel size without decreasing the amount of input light power [26], which is desirable for on-chip applications.

    Sketch of the enhanced light-absorption structure that comprises a cavity-embedding Si layer, two GLs enclosing the cavity, and a substrate. The coordinate system is shown, where the grating periodicity and grooves/lines are along the x and y axes, respectively, and the layer stacking and light impinging directions are along the z axis.

    Figure 1.Sketch of the enhanced light-absorption structure that comprises a cavity-embedding Si layer, two GLs enclosing the cavity, and a substrate. The coordinate system is shown, where the grating periodicity and grooves/lines are along the x and y axes, respectively, and the layer stacking and light impinging directions are along the z axis.

    In this paper, we adopt the η=A modeling, which simplifies the design to one in which only the optics need to be considered. Herewith, rather than being restricted to a specific PD technology, we analyze the structures optically, characterizing the constitutive layers by complex RIs n+ik, where n are the real RIs and kn are the extinction coefficients. This modeling still has widespread use for initial designs of the electronic PDs [16,18,21,22,26,27]. In such cases, optical modeling can be complemented and refined by an electron-device simulation [18,28,29]; however, it is not expected to disprove the optical modeling outcomes. The n and k spectra for the involved materials are found in Refs. [3032]; their values at the CDW are given in Table 1.

    λ0 [μm]Si3N4SiO2Si
    0.81.9962+i  01.4533+i  03.6925+i  0.0065

    Table 1. RIs at the CDW of the Materials that Are Set in the Text

    Here, the design emphasis is put on obtaining the maximal CDW efficiency ηmax=maxη(λ0), while maintaining ta constant per design, and keeping the back-GL cladding’s thickness (similarly to the DBR layers’ thicknesses) fixed at the quarter-wave value tqw100.2  nm. Such dimensions as Λ and the groove width W (constrained for simplicity to be the same for both GLs), etch depths tfg and tbg, as well as the thicknesses of the front-grating cladding tSi3N4 and front and back cavities, df and db, respectively, are the design variables. We observe the SW regime, limiting ab initio the operation wavelengths from below by the Rayleigh wavelength λR=ΛnSiO2, which suppresses the non-specular orders of the grating diffraction into the cavity and so prevents losses due to light escaping from the structure’s sides. Note that this is not compulsory since the non-specular orders can be suppressed with extra constraints; see, e.g., Refs. [21,33].

    3. DESIGN PROCEDURE, OPTIMAL STRUCTURES, AND THEIR EFFICIENCY SPECTRA

    We perform the computer-aided designs with an in-house software run in MATLAB environment. The software consists of two modules: “Simulation,” which contains the codes of an in-house recast rigorous coupled-wave analysis (RCWA) for simulations of the structures merging smooth and grating-patterned layers as described in detail elsewhere (see Ref. [21] and references therein); and “Optimization,” which calls a multi-start optimization algorithm from MATLAB Optimization Toolbox [3436] and a graphical user interface (GUI), which forces the modules to interact and thus drives the simulations and designs as shown in Fig. 2.

    Sketchy flowchart of the optimization process. The GUI Module inputs the trial parameters, assessed as described in the text, fixed parameters, and constraints to the Optimization Module. There, the trial-and-error multi-start algorithm generates the next start points and inputs them into the Simulation Module, which feeds the algorithm back and loops until attaining an optimum.

    Figure 2.Sketchy flowchart of the optimization process. The GUI Module inputs the trial parameters, assessed as described in the text, fixed parameters, and constraints to the Optimization Module. There, the trial-and-error multi-start algorithm generates the next start points and inputs them into the Simulation Module, which feeds the algorithm back and loops until attaining an optimum.

    In order that the optimization would not be blind, we are guided by an approach to obtaining initial (start, trial) design parameters that, in modern terms, are geometric-phase control [20]. Homogenizing the gratings as described in Appendix A, we first design the front and back GLs to be standalone and to be highly and lowly reflective (subsections A and B), respectively, at a normal CDW incidence from a semi-infinite SiO2. In this way, we assess trial values Λ*, W*, tfg*, tbg*, tSi3N4* for the GLs’ parameters shown in Fig. 1. With these, our RCWA code returns the GLs’ reflection-amplitudes’ phases φf and φb, the sum of which Φ gives a rough idea on the geometric phase of an unloaded dual-GL cavity. Φ leads us to estimate tc* for such cavity’s thickness with the round-trip phase condition 4πnSiO2(λ0)tc*/λ0+Φ=2πm; see, e.g., Refs. [16,27]. Though the Si layer is not involved in the start design, the cavity in the final design does adopt it (see Section 2), so given a ta, we hold the bound tc*>ta by tuning m. Finally, the front and back cavity thicknesses, but not tc (see Fig. 1), are optimized with initial values df*, db*, chosen compatible with the relation df*+db*=tc*, where we take db*=df* for simplicity. The start equality proves a posteriori far from the relation between db and df for the optimized dual-GL structures but much closer to that for the optimized dual-DBR structures.

    Due to the very small extinction of Si in the operation range [see kSi(λ0) in Table 1], it might be expected that the Si layer would not be too thin; see, e.g., Refs. [35]. In fact, due to Ref. [27], we figured out two initial values of ta for optimizing the dual-DBR structure λ0/(2nSi)0.108  μm and λ0/(4nSi)0.054  μm. To ultimately maximize ηmax, these values were finally refined by the optimization process to ta1=0.150  μm and ta2=0.025  μm. These were used among the fixed input parameters in optimizing the dual-GL structures (see Fig. 2) to attain ηmax100%. The results of the computer-aided optimal design, with two Si layer thicknesses indicated above are shown in Table 2; the η(λ) spectra are shown in Figs. 3 and 4.

    StructuretbgdbdftSi3N4tfgΛWηmaxQλRNt
    ta1=0.150
    2DBR10.0690.0300.808 (0.995)30333 (63)3.779 (7.346)
    2GLTE10.1200.1301.7100.7920.3050.4390.2520.999197.530.64273.307
    2GLTM10.1730.0321.7641.3800.2150.4680.1850.9991012.60.68473.814
    ta2=0.025
    2DBR20.0940.0530.353 (0.985)23029 (71)3.226 (8.220)
    2GLTE20.1090.2030.9680.2710.3220.4430.2510.999249.220.64471.998
    2GLTM20.5860.0031.2620.0470.4370.4790.1750.9955925.90.69672.460

    Table 2. Parameters of the Structures Designed Shown Below with Two Si Layer Thicknesses Outlined in the Texta

    Efficiency spectra of the designed structures with a 0.15 μm thick Si layer. Broken lines: red, 2GL−TE1; green, 2GL−TM1; full blue line, 2DBR1. The structural parameters are found in Table 2.

    Figure 3.Efficiency spectra of the designed structures with a 0.15 μm thick Si layer. Broken lines: red, 2GLTE1; green, 2GLTM1; full blue line, 2DBR1. The structural parameters are found in Table 2.

    Efficiency spectra of the designed structures with a 0.025 μm thick Si layer. Broken lines: red, 2GL−TE2; green, 2GL−TM2; full blue line, 2DBR2. The structural parameters are found in Table 2.

    Figure 4.Efficiency spectra of the designed structures with a 0.025 μm thick Si layer. Broken lines: red, 2GLTE2; green, 2GLTM2; full blue line, 2DBR2. The structural parameters are found in Table 2.

    While the planar-interface layers are insensitive to the polarization of normally incident light [37], the 1D gratings (see Appendix A) are strongly sensitive to it. Hence, given a ta, the design yields at least two different dual-GL structures, i.e., those for the TE- or TM-polarized irradiation, when the incident electric or magnetic field is overall along the grating grooves, respectively. To have an idea up to what extent they improve the efficiency over the RCE structure, we optimized the latter with the SiO2/Si3N4 DBRs at both sides, using the same technique, while keeping its overall thickness comparable with those of the dual-GL structures. Releasing this constraint, the optimal dual RCE structures reach 100% CDW absorptance at the cost of a huge number of layers (parenthesized in Table 2), which makes them technologically unfeasible. Further, we refer to the thus optimized structures as 2GL−TE, 2GL−TM, 2DBR, and 2DBRtu.

    Table 2 and Figs. 3 and 4 show that the thinner the absorber, the higher the relative gain in ηmax of the designed dual-GL structures compared to 2DBR1 and 2DBR2. It can be inferred from Figs. 3 and 4 that, as regards the full width at half-maximum (FWHM) Δλ, the η(λ) spectra of the 2GL−TE and 2GL−TM are the widest and narrowest, respectively, whereas for the RCE structure, Δλ is in between them, with the same absorber thickness ta. It seems that the η(λ) spectra of the optimized dual-GL structures get sharper and narrower as ta decreases; the extremely small Δλ spectrum of the 2GLTM2 in Fig. 4 is notable. Due to the definition of the damped-resonator quality (Q) factor Q=λ0/Δλ, the Q factors for the designed structures are easily assessed from the data of Figs. 3 and 4. The estimated values are displayed in Table 2 and prove rather high, especially for the TM structures. Note that on-chip Si RCE-PDs were fabricated [5,6] with ta that are quite larger than those of the 2DBR1,2; see Table 2. An RCE-PD with the MgF2/ZnSe front and SiO2/polySi back DBRs [5] showed ηmax34% at λ0, and for one with only SiO2/Si back DBR [6] it was 44% at 0.86 μm, which is notably smaller than the ηmax of 2DBR1; see Table 2.

    4. ELECTROMAGNETIC FIELDS AND POWER FLOW

    It is worth noting that all the studied structures, laterally unlimited (actually, huge times λ0 wide) and several times λ0 thick, are three-dimensional (3D) structures. In this regard, our simulations are 3D, though they are simplified by the independence of the EM fields on y for the dual-GL structures (as with, e.g., all RCWA-based simulations), or on x, y for the dual-DBR ones, which does not imply any two-dimensional (2D) or 1D spatial confinement. The light confinement in the cavity–absorber region, mentioned in the introduction for the RCE structures, and that for the dual-GL ones—both discussed and quantified in this section—is entirely related to the EM fields.

    In the RCE structures, the EM confinement is built up due to the optical resonance in a cavity enclosed by plane-parallel mirrors [16,27] and so are 1D in essence. Since the dual-GL structures under study are assemblies of identical vertical (z) top-to-substrate microscale cuts of width Λ from the structure shown in Fig. 1, another physics, which is an EM field resonant recirculation in micro/nanoscale areas or volumes on the sub-wavelength scales [38], might be more suitable for them. We view the above cut as a “microcavity”; the quotes indicate that it has no real vertical boundaries. Inside the microcavity, the lamella-groove boundary line marks an actual materials interface in the grating regions; continued outside them is a dummy, reminiscent of the above real boundary. The nano/microcavity reasoning has been adopted previously for a simpler grating-based optical resonator [39] and 2D photonic crystal slabs [40]. While the physics in the RCE devices, behind the EM confinement, is well understood, e.g., Refs. [16,27], its current understanding for the micro/nanocavity devices is much less mature and is more particular per structure. In these cases, visualizing the electric (E) and magnetic (H) fields and Poynting vector (S), computed at the CDW, may shed light on the near-field EM recirculation and thus portray the confinement and power flow ad hoc. Doing that for the dual-DBR structures aims primarily at confirming the adopted vision [16,27]. For any layer, a rigorous relation between A and S states that A is the difference between surface-averaged normal components of S (power fluxes) at the layer’s exit and entry.

    For the 2DBR1,2, in which the only nonzero components Ey, Hx, and Sz depend only on z, we plotted |Ey|2, |Hx|2, and |Sz| versus z in Fig. 5. As is seen, while |Sz| is slightly sub-linear in one half and super-linear in the other half of the absorber of 2DBR1, it is visibly linear throughout the absorber of 2DBR2. Overall, sup|Ey|2<sup|Hx|2, which is due to high RI contrast between SiO2 and Si; see Table 1. In the absorber layers, the maxima of |Ey|2 and minima of |Hx|2, and vice versa, are complementary.

    Electric and magnetic fields’ amplitudes squared (left ordinate) and Poynting vector modulus (right ordinate), normalized to the input |Sz|, across the cavity–absorber range of the dual-DBR structures: (a) 2DBR1; (b) 2DBR2. The intersections of the dashed lines with the abscissa mark the Si layer ends. Simulated at the CDW.

    Figure 5.Electric and magnetic fields’ amplitudes squared (left ordinate) and Poynting vector modulus (right ordinate), normalized to the input |Sz|, across the cavity–absorber range of the dual-DBR structures: (a) 2DBR1; (b) 2DBR2. The intersections of the dashed lines with the abscissa mark the Si layer ends. Simulated at the CDW.

    In the dual-GL TE or TM structures, Ey,Hx,Hz0 or Hy,Ex,Ez0, respectively, and S0. Hence, we plotted the surface graphs of the appropriate squared modules of these fields’ components versus (x,z) in Figs. 6 and 7. These show the continuity of H [Figs. 6(c), 6(d), 6(e), 6(f), 7(c)], Ey [Figs. 6(a) and 6(b)], and Ex [Fig. 7(a)], and the discontinuity of Ez [Fig. 7(b)] when crossing the cavity–absorber interfaces, which is in fair agreement with the EM boundary conditions. Also, the quiver maps of S are plotted in Figs. 8 and 9. Due to the a priorix periodicity of |Ei|, |Hi|, and |Si|, restricting the x to within any segment of length Λ, not necessarily [0,Λ], is sufficient to pick out the microcavity. For reasons of scale and convenience, z is also restricted to a lateral slice of the microcavity. Due to a perfect-mirror action of the front GL, the EM fields strongly decay there, and hence it is overall excluded from the scope. In Figs. 6 and 7, the z range covers such parts of the absorber–cavity region, in which the front cavities are top cut to 10% and 25% of their thickness for ta1 and ta2, respectively. In Figs. 8 and 9(a), the z range is extended to also include the back GL. The z range in Fig. 9(b) includes a part of the absorber–cavity region, where the front cavity is top cut to retain 5% of its thickness and top 5% of the back grating z range.

    Normalized to the input values, squared amplitudes of the EM fields versus 0≤x/Λ≤1 and z as defined in the text. (a), (c), (e) For 2GL−TE1: 0.22 μm≤z≤0.671 μm; x/Λ=0.426 is the internal grating lamella-groove interface derived line; (0,1)×(0.35 μm,0.5 μm) is the absorber cut within the microcavity. (b), (d), (f) For 2GL−TE2: 0.209 μm≤z≤0.679 μm; x/Λ=0.433 is the internal grating lamella-groove interface derived line; (0,1)×(0.412 μm,0.437 μm) is the absorber cut within the microcavity. Simulated at the CDW.

    Figure 6.Normalized to the input values, squared amplitudes of the EM fields versus 0x/Λ1 and z as defined in the text. (a), (c), (e) For 2GLTE1: 0.22  μmz0.671  μm; x/Λ=0.426 is the internal grating lamella-groove interface derived line; (0,1)×(0.35  μm,0.5  μm) is the absorber cut within the microcavity. (b), (d), (f) For 2GLTE2: 0.209  μmz0.679  μm; x/Λ=0.433 is the internal grating lamella-groove interface derived line; (0,1)×(0.412  μm,0.437  μm) is the absorber cut within the microcavity. Simulated at the CDW.

    Normalized squared amplitudes of the (a), (b) electric and (c) magnetic fields’ components for the 2GL−TM1 structure shown versus 0≤x/Λ≤1 and 0.273 μm≤z≤0.631 μm as defined in the text. For a simulation convenience, the grating cell was inverted so that x/Λ=0.396 is now the grating lamella-groove interface derived line; (0,1)×(0.305 μm,0.455 μm) is the absorber cut within the microcavity. The normalization and simulation wavelengths are the same as in Fig. 6.

    Figure 7.Normalized squared amplitudes of the (a), (b) electric and (c) magnetic fields’ components for the 2GLTM1 structure shown versus 0x/Λ1 and 0.273  μmz0.631  μm as defined in the text. For a simulation convenience, the grating cell was inverted so that x/Λ=0.396 is now the grating lamella-groove interface derived line; (0,1)×(0.305  μm,0.455  μm) is the absorber cut within the microcavity. The normalization and simulation wavelengths are the same as in Fig. 6.

    Maps of the EM power flow at the CDW in the structures (a) 2GL−TE1 and (b) 2GL−TE2 (Table 2). The quiver ranges are defined in Section 4, and the horizontal dashed lines show the material interfaces. The structure cut is shown in the inset to (a).

    Figure 8.Maps of the EM power flow at the CDW in the structures (a) 2GLTE1 and (b) 2GLTE2 (Table 2). The quiver ranges are defined in Section 4, and the horizontal dashed lines show the material interfaces. The structure cut is shown in the inset to (a).

    Maps of the EM power flow at the CDW in the structures (a) 2GL−TM1 and (b) 2GL−TM2 (Table 2). The quiver ranges are defined in Section 4, and the horizontal dashed lines show the material interfaces.

    Figure 9.Maps of the EM power flow at the CDW in the structures (a) 2GLTM1 and (b) 2GLTM2 (Table 2). The quiver ranges are defined in Section 4, and the horizontal dashed lines show the material interfaces.

    It appears that, when simulating E, H, and S for the microcavity frame [0,Λ], the |Ei|2, |Hi|2 surfaces and S flow maps prove asymmetric as regards their boundaries-adjoined critical features (see further in the present section). An analysis shows that the symmetry is restored while preserving the graphs’ topology, with applying a rigid lateral shift of the microcavity frame by the distance Δ=Λ/4W/2; note that Δ=0 for symmetric grating structures with W=Λ/2. Table 2 affords the following values of Δ/Λ: 0.03702 (2GLTE1), 0.03352 (2GLTE2) and 0.06891 (2GLTM1), 0.06769 (2GLTM2), so the shift is to the right for 2GLTE1,2, and the shift is to the left for 2GLTM1,2. Thereby, the EM fields and Poynting vector were recomputed for the shifted frame, and Figs. 69 show the resulting |Ei|2, |Hi|2 surface graphs and S quivers. In all the following figures and discussions, x still denotes the lateral coordinate, and bare values of the microcavity lateral boundaries, i.e., 0, Λ, are assigned to the respective shifted boundaries.

    A. Dual Grating-On-Layer TE Structures

    1. Electric and Magnetic Field Amplitudes

    For 2GLTE1 and 2GLTE2, the plots of |Ey|2, |Hx|2, and |Hz|2 shown in Figs. 6(a), 6(c), 6(e), and Figs. 6(b), 6(d), 6(f), respectively, display highly peaked surfaces. It is found that the minimum x-direction distance between each two peaks is 0.5, in units of Λ, as may be expected from the x periodicity. The bell-shaped |Ey|2 and |Hx|2 surfaces are strictly confined inside the microcavity. They have four and two peaks located at the cavity–absorber interface lines and inside the absorber, respectively, for 2GLTE1, and vice versa for the 2GLTE2.

    As seen in Figs. 6(e) and 6(f), the |Hz|2 surfaces show three pairs of peaks for 2GLTE1 and three peaks for 2GLTE2, respectively. One of the peaks (pairwise or singular) is internal, whereas the two others adjoin the vertical boundaries of the microcavity. The surface(s) adjoined to one boundary is (are), as expected, strictly symmetric to the one(s) adjoined to the other boundary, attaining there their maximums. In 2GLTE1, the peaks of |Ey|2 and |Hz|2 lie on the cavity–absorber interface lines [see Figs. 6(a) and 6(e)]. The peaks of |Hx|2 lie on the midline of the absorber (z=0.425  μm); see Fig. 6(c). In 2GLTE2, the situation is just the opposite; here the absorber is much thinner than in the 2GLTE1, so that in order to provide 100% resonant absorptance, |E| tends to be as large as possible, which pushes H out of the absorber. The maximums of |Ey|2, |Hx|2 at the vertical microcavity boundaries, over the corresponding z ranges in Fig. 6, are about 1.24, 1.29 and 1.58, 1.50 for 2GLTE1 and 2GLTE2, respectively. In 2GLTE1, |Ey|2 has the absolute maxima at the crossing points of two vertical midlines between the lamella-groove interface-derived line and microcavity boundaries with the horizontal absorber–cavity interface lines, while |Hx|2 has the absolute maxima at the crossing points of the same vertical lines with the absorber’s midline. In 2GLTE2, the maxima points of these functions are swapped.

    In view of the location and sharpness of the |Ey|2 and |Hx|2 peaks, the tight vertical-by-lateral 2D EM fields’ confinement occurs in each microcavity of both structures. Moreover, because the above boundary values of |Ey|2, |Hx|2, relative to their global maximums, are about 1.5%, 0.3% for 2GLTE1, and 0.5%, 0.9% for 2GLTE2, the microcavities in both structures come close to electric and lateral-magnetic isolation. In Figs. 6(e) and 6(f), we observe that the internally peaked |Hz|2 surface’s parts are confined similarly to the |Ey|2 and |Hx|2 surfaces. Because of the boundary-peaked parts of the |Hz|2 surface, Hz penetrates from a given microcavity to any one adjacent to it, so there is no vertical magnetic isolation of the microcavities.

    2. The EM Power Flow

    In Figs. 8(a) and 8(b), we visualize the EM power flow in 2GLTE1 and 2GLTE2, respectively, using the quiver plots of S. The figures display the Poynting vector distributions of those structures, demonstrating vortex-like flows. Both flows consist of two vortices, a central one with counterclockwise S circulation and the other, which is split into two symmetric parts adjacent to the microcavity boundaries with clockwise S circulation, so that the total circulation appears to be zero. The distance between the vorticity centers is 0.5Λ, as expected from the Λ periodicity of |Si|. When moving along the x axis through the microcavity, Sx consistently switches its sign at the vertical lines of x=0.25Λ=x1 and x=0.75Λ=x2. The continuity of Sx dictates that these lines are also the lines of the |Sx| minima. The inner vortices’ centers lie on the midline of the microcavity x=Λ/2, and the boundary-adjoined vortices are centered at the vertical boundaries of the microcavity. For both structures, the vortices are strictly confined between the lateral boundaries of the chosen microcavity slices. Contrary to the inner vortices, the boundary-adjoined ones are unconfined laterally and percolate continuously from the given microcavity to the neighboring ones.

    Both vortices in 2GLTE1 [Fig. 8(a)] consist of two linked sub-vortices sustaining the above noted respective S circulations, which correlates clearly with the pairwise peaks of |Ey|2 in Fig. 6(a) and of |Hz|2 in Fig. 6(e). Figuratively speaking, the vortices start the circulation at the back cavity–absorber interface, stretch vertically to the absorber layer, and make the counter-circulation turn at the absorber–front cavity interface. Between the vortices and apart from their vicinity, the EM power flow inside the absorber layer is dominantly vertical. It crosses the absorber back and forth as a result of the reflection from the GLs. The combination of the vertical flow alongside the lateral one at the vortices’ closure seems to be the key mechanism of the 100% peak absorptance of this structure. Four additional singularities appear at the back grating region, two of which are vortices, and the other two at x=x1 and x=x2 are seemingly flow junctions. The left junction seams to channel a lateral power that flows into it from opposite directions (left and right) to a vertical one that flows away from it, also in opposite directions (up and down), whereas the right junction does exactly the opposite.

    The vortices in 2GLTE2 [Fig. 8(b)] do not cluster into sub-vortices, and the S flow they carry circulates entirely inside the back cavity. Between the vortices and apart from their close vicinity, the S flow is dominantly vertical and is pumped into and out from the absorber layer via the vertical junctions, the midlines of which are again x=x1 and x=x2, respectively. Such an EM power transfer produces lateral waveguide-like right and left S flows inside the absorber layer slice within the microcavity. As seen in Fig. 8(b), the right flow occurs between the junctions’ entrances and is restricted within the microcavity, whereas the left flow is piecewise inside [x2,Λ] and [0,x1] and is responsible for the EM power transfer between the periodically arranged microcavities. This lateral back and forth S flow is the main absorption enhancement mechanism for 2GLTE2. The reason why the lateral S flow occurs in the absorber layer of this structure and does not in that of 2GLTE1 is the very small thickness of the former layer, 6 times smaller than that in 2GLTE1; it appears that with such a thin Si layer, a planar SiO2/Si/SiO2 trilayer will guide only a single light mode. Although the light-absorption enhancement mechanism in the 2GLTE2 structure differs notably from that in the previously discussed 2GLTE1 one, it also provides 100% peak absorptance. Note that the S flow at the back grating region is very similar to that of 2GLTE1, which seems to be due to the closeness of their dimensions. Originating at the back grating, this mixed lateral-vertical flow supports the EM energy confinement and channels the EM power in the microcavity of both TE structures.

    A spectral deviation from the optimal CDW or an alternation of the structure’s optimal dimensions, such as the grating parameters and the absorber thickness, which result in a change of the optical properties, affects the overall and absorber flows (direction and/or magnitude) and the singular points (location and/or existence). In the 2GLTE2 structure, for example, the decrease in absorptance due to these variations is noted as breaking the absorber’s balanced lateral optimal flow and the vortices’ existence and location. Note that the vortices do not appear only in the optimized structures and not discretely at the CDW. Still, moving further away from peak resonances and confinement wavelengths, by changing the incident wavelength or optimal dimensions, conventional flows of no vorticity are likely to be observed.

    B. Dual Grating-On-Layer TM Structures

    1. Electric and Magnetic Field Amplitudes

    The plots of |Ex|2, |Ez|2, and |Hy|2 for 2GLTM1 shown in Figs. 7(a)7(c), respectively, also display the surfaces that are highly peaked but topologically very different from those for the dual-GL TE structures considered above. It is seen that the |Hy|2 surface (disregarding the back cavity layer), similarly to the |Ey|2 one in Fig. 6(a), consists of four highly peaked regions. However, here two of the |Hy|2 surface peaks are relatively weak and sharp, with maximums located at the absorber–front cavity interface, and two others are strong and less sharply peaked inside the absorber layer. In contrast with the |Ey|2 plot, the peaks of the |Hy|2 plot in Fig. 7(c) are asymmetric relative to the absorber layer. Nevertheless, with an appropriate interchange EH, several properties of the present plots in Fig. 7 seem to comply with those in Figs. 6(a),6(c), and 6(e). Namely, the shortest x-direction distance between each two maximum pairs of both |Hy|2 and |Ex|2, and minima of |Ez|2, is 0.5Λ. Also, a strong resemblance between the |Ex|2 and |Hx|2 plots can be seen when comparing Fig. 6(c) to Fig. 7(a). Similarly to the 2GL−TE cases, the two |Hy|2, |Ex|2 maximums lines are at x=x1 and x=x2, and the same holds for the |Ez|2 minima lines.

    The above asymmetry along the z axis may be attributed to the very small thickness (32 nm) of the back cavity layer (see Table 2), in which the following dissimilarity between the TE and TM structures is observed. While all the above discussed EM field amplitudes in the dual-GL cavity TE structures (Fig. 6) increase when progressing up along the z axis, i.e., moving right to left in the figure inside the back cavity layer towards the absorber interface, the EM fields for 2GLTM1 decrease along the similar path. An additional dissimilarity is found in the above-noted discontinuity of |Ez| as seen in Fig. 7(b). Inside the absorber |Ez|2 drops 45 times at most, but it still peaks there up to 20 relative units, which concurrently with high |Ex|2 provides 100% peak absorptance. Due to a close similarity to the 2GLTM1 fields and to avoid repetitions, the 2GLTM2 fields are not shown. However, the EM power flows and the absorption mechanisms of both structures, which appear to be quite different, are discussed in the next subsection.

    2. The EM Power Flow

    The EM power flow quiver of the 2GLTM1 is shown in Fig. 9(a); for convenience, the flow in the absorber and front cavity regions is scaled and displayed in the inset. This time the dominant vortices are vertically centered at the back GL region, one at the x-axis center and the other split between the boundaries, similar to the cavity–absorber region vortices in the TE case (Fig. 8). Inside the chosen microcavity, no strictly confined vortices are observed. In the present structure, the back cavity layer is the thinnest one, and it is not surprising that, similar to the case of the absorber layer in 2GLTE2, a mostly lateral S flow occurs within it. A similar flow proceeds at the lower [upper in Fig. 9(a)] part of the absorber layer, though it is about 6 times less intense due to the discontinuity of Sx in z, while crossing the back cavity–absorber interface z=0.305  μm, quantitatively resulting from the discontinuity |Ez|; see Fig. 7(b). For the same reason, a similar “jump” of Sx occurs when crossing the absorber–front cavity interface z=0.455  μm, which is seen in Fig. 9(a). This behavior contrasts with the overall continuity of Sx in 2GLTE1 and 2GLTE2.

    In this case, the strong lateral S flow in the back cavity layer originates from the EM power pumping out and into the vortices that persist in the back GL region (not shown here). As seen in Fig. 9(a), the inter-layer transfer of the EM power goes down and up via the above junctions along the x=x1 and x=x2 midlines, respectively. Although the S flow over the upper [lower in Fig. 9(a)] part of the absorber is negligibly small as compared to that in the rest, the predominantly lateral S flow over the latter proves sufficiently strong to provide the perfect optical absorption.

    Focusing mainly at the absorber and its closer regions, the EM power flow of the 2GLTM2 structure is shown in Fig. 9(b). As for the previous structures, a laterally centered and a boundary split vortices are displayed, only this time they are both vertically centered inside the absorber layer. Both vortices are stretched vertically throughout the entire absorber layer and close at the cavity layers. The vortices’ closure echoes vertically deep and laterally wide into the back GL and front cavity regions; actually, it appears that the stretched singularity inside the absorber is replaced by a wide lateral flow outside the absorber layer. Between the vortices and apart from their closest vicinity, the EM power flows almost completely vertically up and down inside the absorber. So although it is ultrathin, the high electrical fields and the back and forth flow inside the absorber layer fuel the 100% peak absorptance.

    To conclude the section, maximizing the absorptance peak at the CDW dictates the vortices centers’ locations and the waveguide-like flows within the structures in general, and inside the active layer in particular, which allows attaining the A100%. Let us stress a generic difference between the EM confinement and power flow in the dual-DBR (RCE) and dual-GL structures. Namely, in the former, due to the lateral interfaces, light paths are 1D vertical round-trips along the z axis; see Fig. 4. The lateral and vertical interfaces in the latter concurrently force light to flow along complex 2D paths in the xz plane; see Figs. 8 and 9. For the dual-DBR structures, the zero EM field singularities have little effect on Sz, which is a rather plain function of z, as seen in Fig. 5. However, the EM field zeros in the dual-GL structures, seen in Figs. 6 and 7, do complicate the S flow, which manifests itself in diverse topological structures, e.g. vortices, shown in Figs. 6 and 7.

    5. DESIGN TOLERANCES

    The issue of fabrication tolerances is of primary practical importance in the design of the proposed structures. Since in the manufacturing of these devices only few-layer deposition steps are required, the grating fabrication accuracy is critical. For the dual-GL cavity-enhanced photo-absorbers under study, we carried out the performance tolerability tests, varying tfg, tbg, and W around their optimal values in Table 2. We adopted ±10% relative error in all the etch parameters for the 2GLTE1,2 and 2GLTM1 structures. However, with these, and even twice as small etching errors, the 2GLSTM2 structure looses its absorptive performance. Due to that, we inquired about the performance of this structure against ±2% relative etching error.

    Tables 3 and 4 display the test results for the designed TE- and TM-polarization targeted structures, respectively. Table 3 shows that the structures 2GLTE1,2 sustain high performance, i.e., exhibit a small peak-efficiency deterioration Δηmax concurrently with a tiny peak-wavelength shift Δλ0 on a sub-nanometer scale. Particularly, apart from one exceptional case for each structure, such as with Δηmax=3.4% and Δηmax=1.1% for 2GLTE1 and 2GLTE2, respectively, all other combinations of the errors result in a less than 1% deterioration.

    Etching Errors2GLTE12GLTE2
    ΔWWΔtbgtbgΔtfgtfgΔλ0, nmΔηmax, %Δλ0, nmΔηmax, %
    +0.3+0.0+0.36+0.1
    ++0.3+0.0+0.26+0.1
    +−0.3−3.4−0.23−0.7
    ++−0.3−0.6−0.11−1.1
    ++0.4−0.3+0.32−0.6
    +++0.4−0.3+0.32+0.0
    ++−0.2−0.3−0.20−0.3
    +++−0.2−0.2−0.20−0.3

    Table 3. Performance Changes of the 2GL−TE Structures with Varying Grating Etch Parametersa

    Etching Errors2GLTM12GLTM2
    ΔWWΔtbgtbgΔtfgtfgΔλ0, nmΔηmax, %Δλ0, nmΔηmax, %
    +1.8−16.5−0.8−35.3
    ++1.8−8.3−0.7−53.3
    ++9.5−22.9+3.1−61.6
    +++9.5−18.9+3.6−7.5
    +−9.7−42.4−3.6−23.1
    ++−9.7−51.4−3.4−36.3
    ++−3.8−18.2+0.4−86.3
    +++−3.9−31.2+0.7−59.7

    Table 4. Performance Changes of the 2GL−TM Structures with the Varying Grating Etch Parametersa

    Similar testing of the structures 2GLTM1,2 shows that they still support resonant absorption all over the tests, though with essentially larger |Δλ0|, as seen in Table 4. It is also seen that with the same relative errors, 2GLTM1 deteriorates in the performance one and a half orders of magnitude stronger than 2GLTE1. For 2GLTM1, except for two severe deterioration cases with Δηmax=51.4% and 42.4% (at both Δλ0=9.7  nm), at the other errors’ combinations, |Δηmax|31.2%, which means a still reasonable tolerance. The high performance sensitivity of the 2GL−TM structures, sharper for the 2GLTM2, also correlates well with their superior Q factors (see Table 2) and may be attributed to an alternating charging of the line-groove interfaces by the light electric field component Ex that emerges only in the TM polarization. To shed more light on the above noted deterioration cases for 2GLTM1, we inquired how a +10% error in W only, with the other dimensions kept at their optimal values, affects the peak absorptance and S flow. We found that Δηmax  37%, so just such an error holds a dominant contribution to those cases, while the errors in tbg,fg when switched on concurrently can either enhance or reduce the deterioration as seen in Table 4. Also, at only ΔW/W=+10%, two vortices, managing the EM power flow back to the absorber at the optimal dimensions, disappear (not shown). Instead, the power continues to flow up, and the optimal confinement is disturbed to a large extent. Table 4 shows that for 2GLTM2 the performance tolerance issue is problematic; in a half of the etch errors’ combinations |Δηmax| it exceeds 50%, rising up to 86% in the worst case. This is not surprising due to its highest Q factor and very specific EM confinement mechanism, which are extra sensitive to the structural dimensions.

    6. APPLICATIONS

    Due to Table 2, the designed Si-photonics structures appear to be loaded in high-Q-factor and 100%-peak-absorptance resonators for the chosen λ0, and quite similar design can be carried out for any higher NIR wavelength below 1 μm. Hence, due to their high efficiency and ultra-small thickness, these proposed structures can be used as efficient and thin PDs in FPAs, on-chip integrated Si photonics, and precise spectroscopic chemical sensing. When unloaded, the designed double-GL cavities exhibit a 100% narrowband transmission notch at λ0, so they may be used as ultra-narrowband optical filters for a wide range of applications.

    For the PD applications, electrical contacts need to be connected to the Si layer. Close methods to this end were developed for RCE-PDs, both photodiode and photoconductor ones, e.g., Refs. [3] and [28,41], respectively, and waveguide-coupled PDs, e.g., Refs. [7,23]. Briefly, the technique amounts to ad hoc making the part below the front cavity wider than that above the Si layer, on both sides in the x direction, and then depositing the poly-Si/metal contacts atop the protruding Si strips as shown in Fig. 10. Since the structures are a huge times λ0 (in our case, also Λ) wide, the results obtained without the width offset prove valid with it as experimentally verified [26,28,41].

    Possible placement of electrical contacts to the photoactive Si layer in the PD applications of the dual-GL structure, sketched in Fig. 1. The placement is appropriate for the BSI considered in this study as described in the text.

    Figure 10.Possible placement of electrical contacts to the photoactive Si layer in the PD applications of the dual-GL structure, sketched in Fig. 1. The placement is appropriate for the BSI considered in this study as described in the text.

    Upon an offset of the light incidence from the normal to an oblique small angle θ, the absorption peaks split into two peaks at some λ0<λ0 and λ0+>λ0 for the TE structures and shift to positions at some λ0>λ0 for the TM structures; see the example in Fig. 11 for 2GLTM2. The splitting δTE=λ0+λ0 and shift δTM=λ0λ0 fit accurately with linear and quadratic regressions, respectively, viz. δTE=aθ and δTM=bθ2. Here a and b, each with less than 2% maximum standard relative error, equal 9.818 nm/deg for 2GLTE1, 16.19 nm/deg for 2GLTE2, 2.366  nm/deg2 for 2GLTM1, and 2.395  nm/deg2 for 2GLTM2. Both δTE and δTM performances can be easily calibrated, provided that the quantities are measurable. Thus, the considered dual-GL cavity structures may be used for high-precision angular-position sensing.

    Efficiency spectra θ dependence of the 2GL−TM2 structure. The structural parameters are found in Table 2.

    Figure 11.Efficiency spectra θ dependence of the 2GLTM2 structure. The structural parameters are found in Table 2.

    The above noted η(λ) sensitivity to a varying incidence angle θ may also be beneficial for the WDM, highly important in contemporary optical communication. The exponential growth in the need for ultra-fast and ultra-wideband communication stimulated the development of the WDM. It provides an increase of the capacity by multiplexing a number of optical carriers using multiple wavelengths onto a single optical guide. With the existing electrical and optical infrastructure for long-haul communication, the WDM approach has proven to be highly cost effective and to benefit, for an efficient operation, from such optical components as ultra-narrowband transmitters, receivers, and filters.

    The evolving DWDM requires channel separation that can reach higher than 50 GHz in frequency or lower than 0.8 nm in wavelength. This demands ultra-narrow and stable components acting as demultiplexers to separate each wavelength with minimum crosstalk. Structures such as those proposed and designed in this work can help address this challenge due to their high Q, ηmax, and θ sensitivity. A possible approach for separating the wavelength channels is directly by the dependence of the absorptance spectra on θ. For example, as can be inferred from Fig. 11, 2GLTM2 provides 0.1  nm separation around 0.8 μm. Though our structures can function as DWDM devices only for integrated-circuit short- to medium-haul communication in the range 0.8  μmλ<1  μm, principles quite similar to those developed in this work can be applied for a higher efficiency redesign of the detectors currently in use for the conventional communication band’s DWDM, such as the monolithic InGaAs.

    The emerging field of optical camouflage would also benefit from designs inspired by light-absorbing structures such as those considered in this study, which are ultrathin, subwavelength, and perfectly absorptive. For this purpose, managing the light power to flow laterally in a thin buried absorbing channel, as Figs. 8(b) and 9(a) demonstrate, seems indispensable.

    7. CONCLUSIONS

    We proposed and theoretically investigated a new type of enhanced light-absorption structures for NIR Si-photonics (see Table 2), comprising only Si-microelectronics-compatible materials (SiO2 and Si3N4) and Si (see Table 1) as the NIR light absorber, in contrast with, e.g., Refs. [5,6]. Inspired by the RCE-PD, this structure contains the cavity–absorber (Si) wrap enclosed with the subwavelength GLs. We argued that the assembly of the RCE-PD with the same materials and topmost absorptance is beyond the scope of feasibility. On the contrary, the subwavelength dual-GL cavity enhanced ultrathin structures, attaining 100% ultra-narrowband η(λ) peak at the near IR, were shown to be feasible with the in-house software design.

    We extensively shed light on the multiple physical mechanisms behind the record ηmax of the designed structures, visualizing the resonant EM field amplitudes and power flow over the important space regions inside the structures. We concluded that the EM fields in the cavity–absorber region are tightly 2D confined, which combines the sharpened vertical confinement of the resonant-cavity type with lateral confinement set by the GLs’ derived microcavity. Also, we gave novel comprehensive and structure-peculiar insights into the EM power transfer, which combines vortex and junction crossed waveguide-like flows. The flows point at optical singularities, at which E=0, such as the vorticity centers and junctions, which connects this study to singular optics [25]. These combined singular optical phenomena proved to be the keys to the perfect optical absorption of each of the designed and simulated structures. Furthermore, we tested the structures’ performance tolerances to the grating fabrication errors, and found them to be superlative for 2GLTE1,2; see Table 3. Among the TM-polarization targeted dual-GL enhanced structures, 2GLTM1 tolerates most of the errors an order of magnitude weaker than its TE-polarization counterpart, still exhibiting a practical performance.

    Finally, in Section 6, we discussed at length possible applications of the proposed structures. To our knowledge, using the metasurfaces proposed and scrupulously theoretically investigated in this study conceptually advances design of topmost optical-absorption structures. Last but not least, this study brings a new focus on striking features of the EM field distribution and Poynting vector flow over the metasurfaces-confined micro-structures. The design principles and subsequent analyses developed in this study may be implemented in other fields of the micro- and nano-optics research and applications.

    APPENDIX A: GRATING-ON-LAYER DESIGN

    Starting values of the GLs’ parameters, needed in the computer-aided design input as shown in Fig. 2, are obtained using a trial semianalytical design on the basis of effective-medium and grating-waveguide resonance theories. In these designs, we approximate the true standalone GLs shown in Figs. 12(a) and 13(a) by effective bilayers as shown in Figs. 12(b) and 13(b), respectively, in which the binary-grating layers are described by effective RIs. These RIs nν depend on one of the principal orthogonal light polarizations, viz. ν=s (TE), p (TM).

    Front standalone GL: (a) the original; (b) an effective-medium bilayer substitute.

    Figure 12.Front standalone GL: (a) the original; (b) an effective-medium bilayer substitute.

    Back standalone GL: (a) the original; (b) an effective-medium bilayer substitute.

    Figure 13.Back standalone GL: (a) the original; (b) an effective-medium bilayer substitute.

    Let permittivity of the grating line be εl and the groove be εg; the line duty cycle is f. Then one can use the effective-medium theory (EMT) to obtain nν=εν via a series in powers of (Λ/λ0)2 [42]: εν=εν,0+Δεν,2+o((Λ/λ0)2),in which the zeroth-order (achromatic) approximations read εs,0=fεl+(1 f)εg,εp,0=εgεlfεg+(1 f)εl,and the subsequent second-order corrections are given by Δεs,2=π33f2(1 f)2(εl εg)2(Λλ0)2,Δεp,2Δεs,2=εp,03εs,0εg2εl2.

    High-Reflection GL Design

    In this design, Λ/λ0 is given, and EMT εν, t1, and t2 are estimated with a grating-waveguide resonance theory [43], in which the grating is not completely smoothed out. Rather, it adds to the persistent non-diffracted zeroth-order wave two virtual counter-propagating ±1st-order ones with the lateral wavenumbers kx=±2π/Λ. These waves may be guided in the layers shown in Fig. 12(b), since nSi3N4>nSiO2, and for this grating nSi3N4>nν>1. Interference between the non-diffracted guided waves is the mechanism of the high reflectance [43].

    For example, at ν=s, complying the diffracted waves with the guided modes yields the following equation: iβ2(β12+β0β3)tan k0β1t1+iβ1(β22+β0β3)tan k0β2t2+(β12β3+β22β0)tan k0β1t1 tan k0β2t2=β1β2(β3+β0),where k0=2π/λ0, and the propagation constants are given by β0=1 (λ0/Λ)2,β3=εSiO2 (λ0/Λ)2,β1=εs (λ0/Λ)2,β2=εSi3N4 (λ0/Λ)2.Here εs relates to the front grating, β0 is imaginary, and, concurrently at nSiO2<λ0/Λ<ns, β3 is imaginary, and β1,2 are both real. Then, Eq. (A4) is satisfied by infinitely many pairs of real t1,t2. Solving Eq. (A4) supplies us the trial parameters tfg*,tSi3N4* to be used in the optimization.

    Low-Reflection GL Design

    Due to the theory of optical multilayers [37], for the antireflection of the bilayer with design t2=tqw=λ0/(4nSi3N4), it is uniquely required that t1=tqw and nν=nSi3N4. Yet, such RI matching seems physically dubious, whichever EMT approximation for εν is adopted. However, given a Λ* compatible with the solubility of Eq. (A4), one can increase εν,2 towards εSi3N4 using Eqs. (A1)–(A3) by tuning f. This may prove possible since, due to Eq. (A3), the second-order correction increases εν beyond the achromatic EMT approximation given by Eq. (2). The assessed Λ* and W*=(1 f*)Λ*, where f* is an adopted value of f and tbg*=tqw, all obtained with the data of Table 1, are used among the starting values in the optimization procedure.

    References

    [1] B. D. Clymer, D. Gillfillan. Corrugation gratings for fast integrated complementary metal-oxide semiconductor photodetectors: implementation and diffraction analyses. Appl. Opt., 30, 4390-4395(1991).

    [2] U. Hilleringmann, K. Goser. Optoelectronic system integration on silicon: waveguides, photodetectors, and VLSI CMOS circuits on one chip. IEEE Trans. Electron Dev., 42, 841-846(1995).

    [3] S. S. Murtaza, H. Nie, J. C. Campbell, J. C. Bean, L. J. Peticolas. Short-wavelength, high-speed, Si-based resonant-cavity photodetector. IEEE Photon. Technol. Lett., 8, 927-929(1996).

    [4] J. C. Bean, C. L. Schow, R. Li, H. Nie, J. Schaub, J. C. Campbell. High-speed polysilicon resonant-cavity photodiode with SiO2−Si Bragg reflectors. IEEE Photon. Technol. Lett., 9, 806-808(1997).

    [5] J. D. Schaub, R. Li, C. L. Schow, J. C. Campbell, G. W. Neudeck, J. Denton. Resonant-cavity-enhanced high-speed Si photodiode grown by epitaxial lateral overgrowth. IEEE Photon. Technol. Lett., 11, 1647-1649(1999).

    [6] M. K. Emsley, O. Dosunmu, M. S. Unlu. High-speed resonant-cavity-enhanced silicon photodetectors on reflecting silicon-on-insulator substrates. IEEE Photon. Technol. Lett., 14, 519-521(2002).

    [7] S. M. Csutak, S. Dakshina-Murthy, J. C. Campbell. CMOS-compatible planar silicon waveguide-grating-coupler photodetectors fabricated on silicon-on-insulator (SOI) substrates. IEEE J. Quantum Electron., 38, 477-480(2002).

    [8] L. Pavesi. Will silicon be the photonic material of the third millenium?. J. Phys. Condens. Matter., 15, R1169-R1196(2003).

    [9] B. Jalali, S. Fathpour. Silicon photonics. J. Lightwave Technol., 24, 4600-4615(2006).

    [10] W. N. Ye, Y. Xiong. Review of silicon photonics: history and recent advances. J. Modern Opt., 60, 1299-1320(2013).

    [11] K. Yamada, T. Tsuchizawa, H. Nishi, R. Kou, T. Hiraki, K. Takeda, H. Fukuda, Y. Ishikawa, K. Wada, T. Yamamoto. High-performance silicon photonics technology for telecommunications applications. Sci. Technol. Adv. Mater., 15, 024603(2014).

    [12] D. Thomson, A. Zilkie, J. E. Bowers, T. Komljenovic, G. T. Reed, L. Vivien, D. Marris-Morini, E. Cassan, L. Virot, J.-M. Fédéli, J.-M. Hartmann, J. H. Schmid, D.-X. Xu, F. Boeuf, P. O’Brien, G. Z. Mashanovich, M. Nedeljkovic. Roadmap on silicon photonics. J. Opt., 18, 073003(2016).

    [13] Y. Zou, S. Chakravarty, C.-J. Chung, X. Xu, R. T. Chen. Mid-infrared silicon photonic waveguides and devices. Photon. Res., 6, 254-276(2018).

    [14] R. Marchetti, C. Lacava, L. Carroll, K. Gradkowski, P. Minzioni. Coupling strategies for silicon photonics integrated chips. Photon. Res., 7, 201-239(2019).

    [15] A. Rogalski. Recent progress in infrared detector technologies. Infrared Phys. Technol., 54, 136-154(2011).

    [16] M. S. Ünlü, S. Strite. Resonant cavity enhanced photonic devices. J. Appl. Phys., 78, 607-639(1995).

    [17] Y. El-Batawy, J. Deen. Resonant cavity enhanced photodetectors (RCE-PDs): structure, material, analysis and optimization. Proc. SPIE, 4999, 363-378(2003).

    [18] Y. El-Batawy, J. Deen. Modeling and optimization of resonant cavity enhanced-separated absorption graded charge multiplication-avalanche photodetector (RCE-SAGCM-APD). IEEE Trans. Electron Dev., 50, 790-801(2003).

    [19] K. Hejduk, K. Pierściński, W. Rzodkiewicz, J. Muszalski, J. Kaniewski. Dielectric coatings for infrared detectors. Optica Applicata, 35, 437-442(2005).

    [20] N. Yu, F. Capasso. Flat optics with designer metasurfaces. Nat. Mater., 13, 139-150(2014).

    [21] M. Zohar, M. Auslender, L. Faraone, S. Hava. Novel resonant cavity-enhanced absorber structures for high-efficiency midinfrared photodetector application. J. Nanophoton., 5, 0518248(2011).

    [22] M. Zohar, M. Auslender, L. Faraone, S. Hava. New resonant cavity-enhanced absorber structures for mid-infrared detector applications. Opt. Quantum Electron., 44, 95-102(2012).

    [23] M. M. P. Fard, C. Williams, G. Cowan, O. Liboiron-Ladouceur. High-speed grating-assisted all-silicon photodetectors for 850 nm applications. Opt. Express, 25, 5107-5118(2017).

    [24] C. J. Chang-Hasnain, W. Yang. High-contrast gratings for integrated optoelectronics. Adv. Opt. Photon., 4, 379-440(2012).

    [25] G. J. Gbur. Singular Optics(2016).

    [26] J. G. A. Wehner, C. A. Musca, R. H. Sewell, J. M. Dell, L. Faraone. Mercury cadmium telluride resonant-cavity-enhanced photoconductive infrared detectors. Appl. Phys. Lett., 87, 211104(2005).

    [27] M. S. Ünlü, K. Kishino, H. J. Liaw, H. Morkoç. A theoretical study of resonant cavity-enhanced photodectectors with Ge and Si active regions. J. Appl. Phys., 71, 4049-4058(1992).

    [28] J. Wehner, C. Musca, R. Sewell, J. Dell, L. Faraone. Responsivity and lifetime of resonant-cavity-enhanced HgCdTe detectors. Solid-State Electron., 50, 1640-1648(2006).

    [29] Y. El-Batawy, M. Medhat. Biasing dependent circuit modeling and optimization of resonant cavity enhanced PIN photodetectors (RCE-PIN-PDs). Proc. SPIE, 9656, 96560P(2015).

    [30] E. D. Palik, H. Philipp. Silicon dioxide (SiO2) (glass). Handbook of Optical Constants of Solids, 749-763(1997).

    [31] E. D. Palik, H. Philipp. Silicon nitride (Si3N4) (noncrystalline). Handbook of Optical Constants of Solids, 771-774(1997).

    [32] D. F. Edwards, E. D. Palik. Silicon (Si). Handbook of Optical Constants of Solids, 547-569(1997).

    [33] C. Alonso-Ramos, X. Le Roux, J. Zhang, D. Benedikovic, V. Vakarin, E. Durán-Valdeiglesias, D. Oser, D. Pérez-Galacho, F. Mazeas, L. Labonté, S. Tanzilli, E. Cassan, D. Marris-Morini, P. Cheben, L. Vivien. Diffraction-less propagation beyond the sub-wavelength regime: a new type of nanophotonic waveguide. Sci. Rep., 9, 5347(2019).

    [34] F. Glover, R. Martí, G. A. Kochenberger. Multi-start methods. Handbook of Metaheuristics, 355-368(1997).

    [35] Z. Ugray, L. Lasdon, J. Plummer, F. Glover, J. Kelly, R. Martí. Scatter search and local NLP solvers: a multistart framework for global optimization. INFORMS J. Comput., 19, 328-340(2007).

    [36] MATLAB R2019a and optimization toolbox 8.3(2019).

    [37] O. S. Heavens. Optical Properties of Thin Solid Films(1991).

    [38] K. Vahala. Optical microcavities. Nature, 424, 839-846(2003).

    [39] W. Nakagawa, Y. Fainman. Tunable optical nanocavity based on modulation of near-field coupling between subwavelength periodic nanostructures. IEEE J. Sel. Top. Quantum Electron., 10, 478-483(2004).

    [40] P. Lalanne, C. Sauvan, J. Hugonin. Photon confinement in photonic crystal nanocavities. Laser Photon. Rev., 2, 514-526(2008).

    [41] J. Wang, J. Hu, P. Becla, A. M. Agarwal, L. C. Kimerling. Resonant-cavity-enhanced mid-infrared photodetector on a silicon platform. Opt. Express, 18, 1640-1648(2010).

    [42] S. Rytov. Electromagnetic properties of a finely stratified medium. Sov. Phys. JETP, 2, 466-475(1956).

    [43] S. S. Wang, R. Magnusson. Theory and applications of guided-mode resonance filters. Appl. Opt., 32, 2606-2613(1993).

    Roy Avrahamy, Moshe Zohar, Mark Auslender, Benny Milgrom, Shlomo Hava, Rafi Shikler, "In-depth investigation and applications of novel silicon photonics microstructures supporting optical vorticity and waveguiding for ultra-narrowband near-infrared perfect absorption," Photonics Res. 8, 381 (2020)
    Download Citation