Doubly-resonant $\chi^{(2)}$ nonlinear photonic crystal cavity based on a bound state in the continuum

TL;DR

This paper introduces a novel doubly-resonant photonic crystal cavity utilizing a bound state in the continuum to significantly enhance second-order nonlinear effects, achieving high efficiency and collimated emission.

Contribution

It presents a new design approach for photonic crystal cavities using bound states in the continuum to improve nonlinear frequency conversion performance.

Findings

Enhanced nonlinear conversion efficiency compared to previous designs

Highly collimated far-field emission at both fundamental and second-harmonic frequencies

Potential for compact, efficient nonlinear optical devices

Abstract

Photonic nanostructures simultaneously maximizing spectral and spatial overlap between fundamental and second-harmonic confined modes are highly desirable for enhancing second-order nonlinear effects in nonlinear media. These conditions have thus far remained challenging to satisfy in photonic crystal cavities because of the difficulty in designing a band gap at the second-harmonic frequency. Here, we solve this issue by using instead a bound state in the continuum at that frequency, and we design a doubly-resonant photonic crystal slab cavity with strongly improved figures of merit for nonlinear frequency conversion when compared to previous photonic crystal designs. Furthermore, we show that the far-field emission at both frequencies is highly collimated around normal incidence, which allows for simultaneously efficient pump excitation and collection of the generated nonlinear signal. Our results could lead to unprecedented conversion efficiencies in both parametric down conversion and second harmonic generation in an extremely compact architecture.

Equations

$${\beta }_1=\frac{1}{4}\frac{\int \mathrm{d}\mathbf{r}{ϵ}_0{\sum }_{ijk}{\chi }_{ijk}^{(2)}(\mathbf{r})\left(E_{1i}^{\ast }{E}_{2j}{E}_{1k}^{\ast }+E_{1i}^{\ast }{E}_{1j}^{\ast }{E}_{2k}\right)}{\left(\int \mathrm{d}\mathbf{r}{ϵ}_0{ϵ}_1(\mathbf{r})|{\mathbf{E}}_1{|}^2\right){\left(\int \mathrm{d}\mathbf{r}{ϵ}_0{ϵ}_2(\mathbf{r})|{\mathbf{E}}_2{|}^2\right)}^{1/2}},$$

$${\mathbf{E}}_{n,p}(\mathbf{r})={\int }_{BZ}\mathrm{d}\mathbf{k}{f}_{n,p}(\mathbf{k})e^{i\mathbf{k}\cdot \mathbf{\rho }}{\mathbf{E}}_{\mathbf{k}n}(\mathbf{r}),$$

$$\bar{\beta }=\frac{{\lambda }_1^{3/2}\int \mathrm{d}\mathbf{r}\overline{ϵ}(\mathbf{r})\left(E_{1x}^2{E}_{2z}^{\ast }+E_{1y}^2{E}_{2z}^{\ast }\right)}{\left(\int \mathrm{d}\mathbf{r}{ϵ}_1(\mathbf{r})|{\mathbf{E}}_1{|}^2\right){\left(\int \mathrm{d}\mathbf{r}{ϵ}_2(\mathbf{r})|{\mathbf{E}}_2{|}^2\right)}^{1/2}},$$

$$\frac{P_o}{P_i^2}=\frac{8}{{\omega }_1}{\left(\frac{{\chi }_{\mathrm{eff}}^{(2)}}{\sqrt{{ϵ}_0{\lambda }_1}}\right)}^2|\bar{\beta }{|}^2\frac{Q_{1t}^4}{Q_{1r}^2}\frac{Q_{2t}^2}{Q_{2r}}{r}_{1c}^2{r}_{2c},$$

Full text

I Introduction

Highly efficient second-order nonlinear frequency conversion would heavily impact a number of applications in current nanophotonics research, such as biosensing Freund and Deutsch (1986); Campagnola and Loew (2003), nonlinear spectroscopy Heinz et al. (1982); Rakher et al. (2010), and efficient nonclassical radiation sources Caspani et al. (2017). For these purposes, developing compact devices with high conversion efficiency is of significant interest, and there has been substantial effort towards designing resonant cavity structures with the aim of enhancing the efficiency of nonlinear frequency conversion.

For second-order nonlinear materials, the cavities that offer the strongest enhancement of the frequency conversion are doubly-resonant, i.e. supporting resonances simultaneously at both the first harmonic (FH) and the second harmonic (SH) frequencies Berger (1997); Liscidini and Andreani (2006); Scaccabarozzi et al. (2006); Hayat and Orenstein (2007). Achieving such doubly-resonant cavities for nonlinear frequency conversion represents a long-standing challenge in optical device design. In addition to the stringent requirement of the double resonance, the nonlinear overlap factor, involving the fields of the first and the second harmonic modes and the spatial distribution of the $\chi^{(2)}$ tensor, must be sufficiently large to ensure an efficient frequency conversion process. Doubly-resonant cavities have been demonstrated in on-chip ring-resonator structures Pernice et al. (2012); Bruch et al. (2018), but with relatively large footprint. Recently, further progress was made by exploiting topology optimization to design small-footprint micropost and microring cavities with enhanced theoretical nonlinear efficiency Lin et al. (2016, 2017). Due to the topology optimization, however, these devices incorporate a number of very fine structural features which make them challenging to fabricate in practice.

Photonic crystal (PhC) defect cavities patterned in thin slabs support ultra-long lived resonant modes confined in volumes approaching the diffraction limit (Akahane et al., 2003; Song et al., 2005), which makes them attractive candidates for a number of applications, including nonlinear harmonic generation. In addition, they also provide a high degree of far-field radiation pattern engineering Tran et al. (2009); Portalupi et al. (2010), which could facilitate in-coupling from and out-coupling to free space. Singly-resonant photonic crystal slab cavities patterned in nonlinear materials have already shown very large nonlinear conversion efficiencies Galli et al. (2010); Buckley et al. (2014a); Yamada et al. (2014); Zeng et al. (2015); Mohamed et al. (2017), despite the lack of a confined resonance at the second harmonic. This promises that engineering doubly resonant modes in such structures could yield unprecedented enhancement. However, designing confined modes fulfilling the doubly resonant condition with either one-dimensional nanobeams or two-dimensional photonic crystal slab cavities has proven very challenging Rivoire et al. (2009, 2011); Buckley et al. (2014b); Andreani and Gerace (2006); Minkov and Savona (2014). Photonic crystal slab cavities are usually designed starting from a partial photonic band gap in the guided mode spectrum. This strategy can be successfully applied to design a high-$Q$ cavity at the first harmonic. However, the second harmonic frequency range generally lies entirely inside the light cone of the slab, such that not even a partial photonic band gap is possible. As a result, in spite of several decades of efforts in optimizing photonic crystal cavity designs, there has not been any report of a successful doubly-resonant photonic crystal cavity structure for second-order nonlinear frequency conversion.

In this paper, we introduce a design of a doubly resonant photonic crystal cavity with high modal overlap factor for $\chi^{(2)}$ nonlinear frequency conversion. This design is built upon the concept of bound states in the continuum (BIC), which has recently received renewed interest in photonics Hsu et al. (2016). BICs can be observed in various systems, but in PhC slabs in particular they correspond to modes that lie inside the light cone but are nonetheless non-radiative, fully-guided inside the slab. Thus, the key idea here lies in abandoning the commonly-held notion that having photonic band gaps at both the first and second harmonics is an essential prerequisite to design doubly-resonant PhC cavities. Indeed, it was recently shown that a heterostructure cavity can be realized in a two-dimensional PhC slab without a band gap Ge et al. (2018). The concept is illustrated in Fig. 1(a)-(b). The starting structure is a regular PhC made by a hexagonal lattice with lattice period $a$ of air holes of radius $r$ in a slab with thickness $d$ and refractive index $n$ (Fig. 1(a)). A heterostructure design can then be formed by introducing PhC regions with different hole radii, as shown in Fig. 1(b). Even with no band gap, modes localized in the core region of such a structure can still be observed close to a band-edge frequency of the photonic bands of the underlying PhC. When these band-edge modes have a diverging quality factor due to a BIC effect, a large $Q$ can also be expected for the heterostructure modes – in particular in the limit of increasing core size Ge et al. (2018). In our previous work Ge et al. (2018), we also found that having a transition region between the core and the outer regions (Fig. 1(b)) further improves the quality factor of the heterostructure cavity.

Below, we use this route to design a high-$Q$ localized mode at the SH frequency, and to precisely tailor the doubly resonant condition, starting from photonic bands that fulfill the right symmetry properties for efficient $\chi^{(2)}$ nonlinear frequency conversion. To show the results on a practical photonic platform, we explicitly choose the slab material to be Gallium Nitride (GaN), in which efficient continuous wave SH generation has been observed for singly resonant photonic crystal slab cavities with FH in the near infrared Zeng et al. (2015); Mohamed et al. (2017). The design strategy hereby introduced could be extended to different wavelength ranges as well as other material platforms. In addition to the applications in frequency conversion, these results could also prove important towards exploiting second-order nonlinearities to reach the regime of photon blockade Majumdar and Gerace (2013); Gerace and Savona (2014), thus opening the door to the engineering of single-photon Fock states in integrated quantum photonics.

II Recipe for band engineering

We begin our study by designing the photonic bands of the GaN photonic crystal slab as shown in Fig. 1(a) to satisfy the following conditions:

A band with a band edge mode outside the light cone at the FH frequency $\omega_{1}$. 2. 2.

A band with a singly-degenerate band edge mode at the $\Gamma$-point at a frequency $\omega_{2}\approx 2\omega_{1}$. 3. 3.

The band-edge modes in 1. and 2. are simultaneously either a maximum or a minimum of their corresponding bands. 4. 4.

The periodic parts of the electric field of these two band-edge modes have a nonzero nonlinear overlap factor.

Below we discuss the rationales behind these conditions. The nonlinear overlap factor between two resonant modes for second harmonic generation is defined as Rodriguez et al. (2007)

[TABLE]

where $\chi^{(2)}_{ijk}$ is the second-order nonlinear susceptibility tensor for the material, $\mathbf{E}_{1}$ and $\mathbf{E}_{2}$ are the mode profiles and $\epsilon_{1}(\mathbf{r})$ and $\epsilon_{2}(\mathbf{r})$ are the relative permittivity profiles at the FH and the SH frequencies, respectively. Here, we take into account that the permittivities in general are different for the first and second harmonic modes. Assuming the GaN slab is grown with its c-axis along the $z$-direction (orthogonal to the slab plane), the dominant $\chi^{(2)}$ components are $xxz$, $yyz$ and $zzz$ Sanford et al. (2005). In a PhC that is symmetric with respect to the $xy$-plane bisecting the slab, the positive- and negative-symmetry photonic modes are predominantly $E_{x,y}$ and $E_{z}$ polarized, respectively (quasi-TE and quasi-TM). Thus, the two modes in Conditions 1 and 2 can be either both quasi-TM (coupling through $\chi^{(2)}_{zzz}$), or quasi-TE at the FH and quasi-TM at the SH (coupling through $\chi^{(2)}_{xxz}$ and $\chi^{(2)}_{yyz}$). Since the hexagonal lattice PhC slab has a photonic band gap for the quasi-TE modes at low frequencies, which could help with the localization and the quality factor of the heterostructure modes, we will look for a quasi-TE FH and a quasi-TM SH mode. Of course, for other materials or other orientations of the GaN slab, the relevant $\chi^{(2)}$ components have to be considered when making this choice.

In the end of this section, we will also discuss how Conditions 1 and 2 ensure not only the resonance matching ($\omega_{2}\approx 2\omega_{1}$), but also that the heterostructure modes are derived from a photonic crystal mode with an infinite quality factor, and hence can have high $Q$-s at both the FH and the SH frequencies. Condition 3 is simply required to ensure that a heterostructure cavity will localize light in its central region for both frequencies, and needs no further discussion. Condition 4 imposes mode-matching for the periodic parts of the Bloch modes, and will be further discussed below.

We use the guided-mode expansion method Andreani and Gerace (2006) to scan the parameters $d$ and $r$ for the photonic crystal in order to find bands that match all four conditions above. The material dispersion is taken into account by setting two different values of the slab refractive index for the quasi-TE and the quasi-TM simulations. Namely, we set $n=2.32$ and $n=2.38$, as for GaN at wavelengths around $1.3$µm and $0.65$µm, respectively. In Fig. 1(c), we show the photonic band structure for the quasi-TE modes of the PhC with $d=0.28a$, $r=0.22a$ and $n=2.32$, while in panel (d) we show the quasi-TM modes of the same PhC but with $n=2.38$. As can be seen, there is a pair of bands that satisfy conditions 1, 2 and 3 from our list, highlighted in green in panel (c), and in red in panel (d).

Beyond the polarization matching, it is also crucial to consider the fact that $\beta_{1}$ is proportional to a spatial integral over the local field contributions that can add up constructively as well as destructively. This generally imposes phase-matching conditions for guided modes, while, for localized modes where the average $\mathbf{k}$-vector is zero, the mode profiles must still be engineered such that $\beta_{1}$ does not vanish. In other words, even for modes that are co-localized in the same region of space, we must also ensure that the integral in the numerator of eq. (1) is large. For a heterostructure cavity, a localized mode derived from a single band is expected to be well approximated by

[TABLE]

where $\bm{\rho}$ is the in-plane component of $\mathbf{r}$, $\mathbf{E}_{\mathbf{k}n}(\mathbf{r})e^{i\mathbf{k}\cdot\bm{\rho}}$ are the Bloch modes of the $n$-th band from which the heterostructure mode is derived, and $f_{n,p}(\mathbf{k})$ is some envelope function. The index $p$ reflects the fact that a number of envelope functions could be supported. Assuming that $f_{n,p}(\mathbf{k})$ is narrowly centered around the band-edge, as must be the case in the limit of a heterostructure size spanning multiple lattice periods, we can approximate $\mathbf{E}_{\mathbf{k}n}(\mathbf{r})\approx\mathbf{E}_{\mathbf{k}_{0}n}(\mathbf{r})$, where $\mathbf{k}_{0}$ is the band-edge Bloch vector. Note that, due to the $C_{6}$ and time-reversal symmetries, $\mathbf{E}_{\mathbf{k}n}(\mathbf{r})$ is the same for any of the possible $\mathbf{k}_{0}$ choices for a band edge at the M or the K points, and must also be real-valued.

Given these considerations, in Supplement 1 we study the overlap terms that need to be non-vanishing for FH and SH modes of the form of eq. (2). One of the terms is between the envelope functions $f(\mathbf{k})$, and we find it to be always non-zero when the SH envelope $f_{2,p^{\prime}}(\mathbf{k})$ is centered around $\Gamma$. Intuitively, this result stems from the fact that, in a homogeneous material, an SH guided mode at $\Gamma$ (i.e. constant $\mathbf{E}$-field in real space) is phase-matched with an FH standing-wave field of the form $\cos(\mathbf{k}\rho)$ for any $\mathbf{k}$. This is because the FH field would enter the overlap term in eq. (1) as $\cos^{2}(\mathbf{k}\rho)$. The second overlap integral that must be non-vanishing is between $\mathbf{E}_{\mathbf{k}_{0}1}(\mathbf{r})$ and $\mathbf{E}_{\mathbf{k}_{0}^{\prime}2}(\mathbf{r})$, i.e. the periodic parts of the PhC Bloch modes at the FH and SH frequency. Looking at eq. (1), we note that, regardless of the symmetries of the $\mathbf{E}_{\mathbf{k}_{0}1}$ profile, the spatial dependence of $E_{\mathbf{k}_{0}1x}^{2}$ and $E_{\mathbf{k}_{0}1y}^{2}$ is always even with respect to reflections both in the $xz$ and in the $yz$ planes, as can is also evident in Fig. 1(e). Thus, condition 4 translates to the requirement that the $E_{\mathbf{k}_{0}^{\prime}2z}$ profile must also be even with respect to both reflections. We took this into account in designing our PhC bands, as can be seen in Fig. 1(f).

Finally, since the mode in a heterostructure cavity is derived from the band-edge modes in a photonic crystal as shown in the discussions below Eq. (2), to form a high-$Q$ heterostructure cavity it is desirable that the corresponding band-edge mode has infinite $Q$. We note that there are two ways to achieve an infinite $Q$: first, band modes outside the light cone of a PhC slab (gray shaded region in Fig. 1(c)-(d)) have infinite intrinsic lifetime and are perfectly confined around the slab region Johnson et al. (1999). In contrast, inside the light cone there is a continuum of radiation modes, as well as a number of quasi-guided bands that are partially confined in the slab, but most generally have a finite radiative lifetime. Even in that region, however, it is possible to have non-radiative bound states in the continuum, either due to symmetry or due to an accidental destructive interference of the radiated light Hsu et al. (2016). In particular, any singly-degenerate mode at the $\Gamma$ point of the band structure must necessarily be non-radiative due to symmetry Ochiai and Sakoda (2001); Fan and Joannopoulos (2002). This is what motivates our particular requirement in Condition 2 for a singly-degenerate mode at $\Gamma$. To confirm our expectations, in Fig. 1(g), we plot the radiative quality factors corresponding to the two bands highlighted in Fig. 1(c)-(d). As can be seen, the $Q$ of the green band at the band minimum at the $M$-point goes to infinity since this mode lies outside the light cone, while the $Q$ of the red band at the band minimum at $\Gamma$ goes to infinity since the mode is a BIC.

III Cavity modes and figures of merit

On the basis of these regular PhC modes, we can now introduce the cavity modes through the heterostructure design schematically shown in Fig. 1(b). We define three hexagonal regions: core, transition, and outer, with radii $r_{c},r_{t},r_{o}$, respectively. The size of the regions is defined by the number of hexagonal ‘layers’ of holes in each of them $(N_{c},N_{t},N_{o})$, or, alternatively, through the side-length of the hexagons, as explained in Fig. 1(b). Because the two bands in Fig. 1(c)-(d) have a minimum at the band-edge frequency, the core region of the heterostructure is surrounded by holes with a larger hole radius. In the cavity design that we study, we set $(r_{c},r_{t},r_{0})=(0.22a,0.23a,0.25a)$ (see Supplement 1; also note that in the schematic of Fig. 1(b), the difference of the radii is larger for visualization purpose). The dotted lines in Fig. 1(c)-(d) thus show the bands for the regular PhC in the outer region of the cavity, which lie, as needed, at a higher frequency than the core-region (solid) bands.

To simulate the electromagnetic modes of the heterostructure, we use Lumerical FDTD solutions – a commercial finite-difference time-domain solver (see Supplement 1). These modes are shown in Fig. 2(a)-(c) for $(N_{c},N_{t},N_{o})=(6,4,14)$. In panel (a), we show the electric field of the SH mode at reduced frequency $\bar{\omega}_{2}=\omega_{2}a/(2\pi c)=0.889$, corresponding to wavelength $\lambda_{2}=652$nm for a lattice constant $a=580$nm. At the FH frequency, we find two degenerate modes at $\bar{\omega}_{1}=0.444$, shown in Fig. 2(b)-(c). These profiles can be understood as heterostructure modes described by eq. (2), with two envelope functions $f_{1,1}(\mathbf{k})$ and $f_{1,2}(\mathbf{k})$. In our FDTD simulation, we also find higher-order modes at the FH frequency, but there is more than $5$nm spectral difference between those and the ones shown in Fig. 2(c)-(d), even for the largest core size $N_{c}=14$ studied below. At the SH frequency, the only mode that we find to be resonant in the core region and with the correct mode-matching symmetries for SHG is the one shown in Fig. 2(a).

In Fig. 3(a), we show the radiative quality factors $Q_{r}$ of the three modes of the heterostructure as a function of the core region size $N_{c}$, while keeping $N_{t}=4$, $N_{o}=14$. As expected $Q_{r}$ increases with increasing $N_{c}$, and becomes particularly large for the FH modes – approaching one million. The SH mode has a quality factor ranging from $1,100$ to $5,900$, which is also remarkably high given the fact that the mode is deep inside the light cone of the PhC, and that there is no SH photonic band gap. In Fig. 3(b), we plot the two important figures of merit for nonlinear second-harmonic generation: the product $Q_{1r}^{2}Q_{2r}$, as well as a dimensionless version $\bar{\beta}$ of the nonlinear overlap factor of eq. (1) as defined in Ref. Lin et al. (2016):

[TABLE]

where $\lambda_{1}$ is the free-space wavelength at the FH frequency, and $\bar{\epsilon}(\mathbf{r})$ is defined to be 1 in the nonlinear dielectric, and 0 outside. This dimensionless overlap factor allows to compare this figure of merit independently of the magnitude of the $\chi^{(2)}$ values and of the FH wavelength. We note that both $Q_{r}$ and the $\bar{\beta}$ associated to the two modes $\omega_{1,1}$ and $\omega_{1,2}$ are expected to be the same due to symmetry, and the slight difference is likely due to the FDTD simulation mesh and boundary conditions, which break the hexagonal symmetry. Furthermore, note that with larger $N_{c}$, the cavity becomes multimode with a decreasing frequency separation between the modes. These additional modes correspond to higher-order envelope functions $f_{n,p}(\mathbf{k})$ in eq. (2).

Finally, we consider a realistic experimental realization of second harmonic generation in our cavity, and estimate the conversion efficiency in the undepleted pump limit. For this, we will need to estimate the coupling coefficients of the radiative channels that will be used for pump in-coupling and SH signal out-coupling. Thus, in Fig. 3(c)-(d) we show the far-field emission profiles of the FH mode Fig. 2(b) and the SH mode of Fig. 2(a), respectively. In Supplement 1, we also show these profiles for a range of core size values for all three modes. These plots illustrate one additional advantage of our design in terms of second harmonic generation: namely, all emission profiles have a strong lobe within a small angle around the vertical direction. This holds particularly true for the SH mode, because its envelope function $f_{2}(\mathbf{k})$ is centered around $\Gamma$. Note that the emission at exactly zero angle is vanishing because of the BIC at $\Gamma$. Importantly, however, all of the generated second harmonic signal can be collected with a small numerical aperture from above the cavity, which is a significant improvement to previous studies using a photonic crystal platform. Furthermore, the strong lobe of the FH mode at zero emission angle facilitates free-space in-coupling of the pump radiation. We note that the vertical emission at the FH frequency is generally not guaranteed, but it is nevertheless observed for our particular design. In general, this can also be induced (or strengthened) through cavity modifications that specifically introduce $\mathbf{k}=0$ components in the Fourier transform of the cavity modes Tran et al. (2009); Portalupi et al. (2010).

These considerations imply that efficient in- and out-coupling to the cavity can be achieved by exciting and collecting within a small angle around the vertical direction (schematically shown in Fig. 3(f)). To study the frequency conversion efficiency in such an experimental setup, we consider in-and out-coupling through part of the far field, assuming that a fraction of the radiative quality factor of the FH mode is due to coupling to the input channel. We define this fraction $r_{1c}=Q_{1r}/Q_{1c}$, where $Q_{1}^{c}$ is due to coupling to the pump mode and $Q_{1}^{r}$ is the radiative $Q_{r}$ of Fig. 3(a). Similarly, we write $r_{2c}=Q_{2r}/Q_{2c}$, where $Q_{2}^{c}$ is associated to the far-field modes from which the SH signal is collected. Furthermore, we include a phenomenological extrinsic quality factor $Q_{e}$, accounting for effects like linear loss or disorder-induced scattering, and write the total $Q_{t}$ of the modes as $Q_{it}^{-1}=Q_{ir}^{-1}+Q_{ie}^{-1}$ for $i=1,2$. In the limit of undepleted pump power $P_{i}$, the conversion efficiency $P_{o}/P_{i}^{2}$ for the SH output power $P_{o}$ can then be written as (see Refs. Rodriguez et al. (2007); Lin et al. (2016) and Supplement 1):

[TABLE]

where $\epsilon_{0}$ is the permittivity of free space and $\chi^{(2)}_{\mathrm{eff}}$ is determined by the magnitude of the relevant elements in the susceptibility tensor. The conversion efficiency is maximized when there are no extrinsic losses, $Q_{ie}\rightarrow\infty$, and when the two ratios $r_{ic}$ go to unity, which happens when the pump mode completely overlaps with the far-field radiation pattern at the FH, and all the radiated SH signal is collected.

In Fig. 3(g), we plot the conversion efficiency vs. the core size $N_{c}$, assuming no extrinsic losses, for three different values of $r_{ic}$, reflecting three different scenarios. In the conservative case (red), $r_{c}=[r_{1,c},r_{2,c}]=[0.2,0.5]$, we assume excitation and collection from only one side of the cavity, and only a partial overlap of the pump mode with the far field pattern. In the optimized case (blue) $r_{c}=[0.4,1.0]$, we still assume pumping from one side only, but with a better in-coupling of the pump mode, which can be achieved through far-field engineering. Furthermore, we assume that the SH signal is collected from both sides of the slab for optimal collection efficiency. Finally, the best-case scenario (green) assumes perfect in- and out-coupling at both frequencies, which would require either pumping the cavity from both sides or introducing a reflector on one side, and a perfect far-field overlap in either case. Notably, the values for the conversion efficiency exceed 100$\%/W$ even in the conservative case, and go beyond 10,000$\%/W$ in the best-case scenario and $N_{c}>10$.

In Fig. 3(h), we show the effect of extrinsic loss on the maximum achievable conversion efficiency for three different values of $Q_{e}$. This assumes $r_{c}=[1,1]$, as well as some far-field engineering such that $Q_{1r}$ matches $Q_{e}$ (for details see Supplement 1). One notable effect of the extrinsic losses is to set the optimal core size of the cavity. Specifically, while the nonlinear overlap factor decreases with $N_{c}$ because of the increasing cavity mode volume (Fig. 3(b)), in the absence of external losses this is more than compensated by the increasing $Q_{1r}^{2}Q_{2r}$ product, such that the efficiency in Fig. 3(g) is monotonically increasing with $N_{c}$. This is also the case for $Q_{e}=10^{6}$ in Fig. 3(h). However, this changes for lower extrinsic quality factors, which limit the maximum total $Q_{t}$. Thus, the efficiency peaks at $N_{c}=4$ for $Q_{e}=10^{5}$, and $N_{c}=2$ for $Q_{e}=10^{4}$. We note that a total (i.e. loaded) quality factor of $44,000$ has already been experimentally measured for a GaN PhC cavity Mohamed et al. (2017), suggesting that $Q_{e}$ approaching $10^{5}$ should be realistic in state-of-the-art systems.

IV Conclusion

We conclude by discussing the merits of the cavity presented here in comparison to previous works. The quantity $\bar{\beta}$ shown in Fig. 3(b) is defined in a way that is particularly useful for comparing to designs in other materials and wavelengths. A comparison to the devices summarized in Ref. Lin et al. (2016) shows that the cavity presented here is orders of magnitude better than previous PhC designs either in the $Q$-product, or in the nonlinear overlap $\bar{\beta}$, or in both. In fact, the only designs that outperform our cavity in terms of $\bar{\beta}$ are the topology-optimized micropillars Lin et al. (2016) and microrings Lin et al. (2017), which are, however, more challenging to fabricate in practice. Furthermore, the topology-optimized micropillars Lin et al. (2016) exhibit a significantly lower $Q_{1}^{2}Q_{2}$ product and therefore a much lower conversion efficiency as per eq. (4) and FOM1 defined in Lin et al. (2016), while the multi-ring designs of Ref. Lin et al. (2017) cannot be side-coupled to an input/output channel, and also offer little engineering possibility for out-of-plane emission. Alternatively, large microring cavities could support modes with a $Q_{1}^{2}Q_{2}$ product exceeding the one of our design, but at the expense of a much larger footprint Bruch et al. (2018). In short, our design combines a high nonlinear overlap with a high $Q$ product in a very compact cavity, with the additional advantage of strong vertical emission at both frequencies.

More broadly, we have outlined a general recipe for designing high-$Q$, doubly-resonant photonic crystal cavity modes with a high nonlinear overlap that can be applied to various materials and setups. This could be useful to further improve the design presented here, as the fact that a lot of the field of both modes is in air indicates that there is enormous room for improvement of the overlap factor. Thus, for example by finding another suitable pair of bands according to our recipe, by exploring PhCs formed by dielectric rods in air, or by infiltrating the PhC holes with a nonlinear material, the nonlinear overlap could be further increased to record-high values. Finally, we further notice that the doubly resonant cavity design proposed here naturally allows for strong field enhancement in the air holes, both at FH and SH. This unique characteristic might be exploited, e.g., for strongly enhanced nonlinear sensing of functionalized surfaces Estephan et al. (2008, 2010). All in all, we believe that the generality of our approach coupled with the well-established fabrication methods for photonic crystal slabs could enable extremely compact devices for second-order nonlinear conversion with unprecedented efficiency.

V Acknowledgements

This work was supported by a MURI project of the Air Force Office of Scientific Researches (FA9550-17-1-0002) and the Swiss National Science Foundation through Project No P300P2_177721. DG acknowledges enlightening discussions with M. Liscidini and M. Galli, and financial support from the EU H2020 QuantERA ERA-NET Cofund in Quantum Technologies project CUSPIDOR, cofunded by the Italian Ministry of Education and Research (MIUR).