Two-dimensional plasmonic waveguides for nanolasing and four-wave mixing

TL;DR

This paper explores the design of two-dimensional plasmonic waveguides tailored for specific applications like nanolasers and four-wave mixing, emphasizing the importance of application-guided design for optimal performance.

Contribution

It introduces application-specific design principles for plasmonic waveguides, highlighting hybrid structures for FWM and wedge structures for nanolasers, with insights on buffer layer choices.

Findings

Hybrid plasmonic waveguides are optimal for four-wave mixing.

Metallic wedges are preferred for nanolasing.

High-index buffer layers improve nanolaser performance.

Abstract

Plasmonic waveguides are an essential element of nanoscale coherent sources, including nanolasers and four-wave mixing (FWM) devices. Here we report how the design of the plasmonic waveguide needs to be guided by the ultimate application. This contrasts with traditional approaches in which the waveguide is considered in isolation. We find that hybrid plasmonic waveguides, with a nonlinear material sandwiched between the metal substrate and a high-index layer, are best suited for FWM applications, whereas metallic wedges are preferred in nanolasers. We also find that in plasmonic nanolasers high-index buffer layers perform better than more traditional low-index buffers.

Figures

Equations

$$r_{\underline{\mathrm{D}}2M}\approx \left|\frac{{\chi }_{\underline{\mathrm{D}}}^{(3)}}{{\chi }_{\mathrm{M}}^{(3)}}\right|\frac{{\int }_{\underline{\mathrm{D}}}|\mathbf{e}{|}^4\mathrm{d}A}{{\int }_{\mathrm{M}}|\mathbf{e}{|}^4\mathrm{d}A},$$

$$A_{\mathrm{eff}}\equiv \frac{{\int }_{\infty }(\mathbf{e}\times {\mathbf{h}}^{\ast })\cdot \hat{z}\mathrm{d}A}{\max \{c{\epsilon }_0{n}_0|\mathbf{e}{|}^2{\}}_{\underline{\mathrm{D}}}},$$

$$k_0/g_{\mathrm{th}}=\frac{{\int }_{\underline{\mathrm{D}}}{n}_0^2|\mathbf{e}{|}^2\mathrm{d}A}{n_0{\int }_{\mathrm{M}}{\epsilon }_{\mathrm{r},\mathrm{M}}^{\prime \prime }|\mathbf{e}{|}^2\mathrm{d}A},$$

$$\mathcal{F}/\Delta n_{\max }=\frac{{\int }_{\underline{\mathrm{D}}}{n}_0^2|\mathbf{e}{|}^2\cdot U\mathrm{d}A}{n_0{\int }_{\mathrm{M}}{\epsilon }_{\mathrm{r},\mathrm{M}}^{\prime \prime }|\mathbf{e}{|}^2\mathrm{d}A},$$

$$k_0/g_{\mathrm{th}}$$

$$\mathcal{F}/\Delta n_{\max }$$

$${\Gamma }_{\mathrm{G}}=\frac{Z_0{\int }_{\underline{\mathrm{D}}}{n}_0^2|\mathbf{e}{|}^2\mathrm{d}A}{n_0{\int }_{\infty }(\mathbf{e}\times {\mathbf{h}}^{\ast })\cdot \hat{z}\mathrm{d}A};\frac{{\mathrm{EFF}}_{\mathrm{NL}}}{f_{\ell }}=\frac{Z_0{\int }_{\underline{\mathrm{D}}}{n}_0^2|\mathbf{e}{|}^2\cdot U\mathrm{d}A}{n_0{\int }_{\infty }(\mathbf{e}\times {\mathbf{h}}^{\ast })\cdot \hat{z}\mathrm{d}A},$$

$$k_0/g_{\mathrm{th}}$$

$${\mathrm{EFF}}_{\mathrm{NL}}/f_{\ell }$$

$${\overline{\rho }}_{\mathrm{G}}$$

$${\overline{\rho }}_{\mathrm{NL}}$$

Full text

1Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen 518055, Guangdong Province, China

2Institute of Photonics and Optical Science (IPOS), School of Physics, The University of Sydney, NSW 2006, Australia

3The University of Sydney Nano Institute, The University of Sydney, NSW 2006, Australia

*Corresponding author: gy.li@siat.ac.cn

1 Introduction

Waveguide-based nanoscale coherent light sources are essential to nanophotonic circuits in diverse applications including high-speed on-chip interconnects [1], quantum information processing [2], chemical and biomedical sensing [3]. Plasmonic waveguides are key elements of such devices since they allow subwavelength light confinement, and hence greatly enhanced light-matter interactions [4, 5]. We consider two different classes of plasmonic-waveguide based coherent nanoscale light sources: nanolasers [6, 7, 8], and nonlinear devices for frequency conversion, specifically degenerate four-wave mixing (DFWM) [9, 10, 11].

Although requirements on plasmonic waveguides for nanolasing and for DFWM applications are distinct, and are also different from those for pure waveguiding applications, plasmonic waveguides are traditionally characterized by their attenuation and by their mode confinement, independently of the ultimate application [12]. For example, after comparing four typical plasmonic waveguides, Oulton et al. found that hybrid plasmonic metal/low-index insulator/high-index semiconductor/air (known as MISA) configurations, exhibit low attenuation loss and high mode confinement [12]. Subsequently, MISA-based devices have been widely used in plasmonic nanolasers [13, 14, 15, 16, 17, 18]. Encouraged by the attractive waveguiding characteristics, as well as the successful applications in nanolasers, the use of hybrid plasmonic waveguides has been extended to nonlinear applications [19, 20, 9, 21], but without a rigorous comparison to other types of plasmonic waveguides.

We recently developed a comprehensive theory for the design and analysis of plasmonic waveguides for nanolasing, and applied this theory, after a slight modification, to the waveguide design for plasmonic FWM devices as well [22]. Our exploration based on one-dimensional (1D) plasmonic waveguide geometries not only provides deep understanding of currently used designs, but also led to potentially superior structures, with the potential to address longstanding challenges in plasmonic nanolasers. Although one- and two-dimensional plasmonic waveguides share the same physics, edge effects introduced by additional dimension may change the modal properties of plasmonic waveguides [23] and accordingly the nanolasing or DFWM performance.

Here we characterize four realistic 2D plasmonic waveguides for nanolasing and DFWM applications. The four structures, illustrated in Figure 1(a)–(d), are: a metallic wedge waveguide, a metal-slot waveguide, a dielectric loaded surface plasmon polariton (DLSPP) waveguide, and a hybrid plasmonic waveguide.

The outline of this paper is as follows. In Section 2 we compare the nonlinear contributions due to the metal with that due to a nonlinear dielectric, and show that the nonlinearity of the metal can be ignored in these DFWM devices. In Section 3 we then briefly introduce the theory we previously developed [22], highlighting the characteristic measures of plasmonic waveguides for nanolasing and DFWM applications. The characteristics of plasmonic waveguides for use in nanolasers and DFWM devices can then be compared and understood in Section 4, allowing us to identify the desirable plasmonic waveguides for these applications. Section 5 includes a discussion of the effects of the refractive index of buffer layers that are used in nanolasers to avoid quenching.

2 Nonlinear contributions from constituent materials

In nonlinear plasmonic waveguides all constituent materials contribute to some degree to the overall nonlinearity. However, since our theory [22] assumes that the contribution from a highly nonlinear dielectric dominates, we need to validate this assumption before we apply this theory to a specific plasmonic waveguide structure. This is especially important since recent studies showed that the third-order nonlinear susceptibility $\chi^{(3)}$ of metals can be large at optical wavelengths [24, 25, 26].

To carry out this comparison we follow the approach by Baron et al. for surface plasmon polaritons [27], which we recently extended to more complicated 1D geometries [28]. For simplicity, we take only one dielectric to be highly nonlinear, and neglect nonlinear contributions from other dielectrics in the plasmonic waveguide. The ratio of the nonlinear contribution from the highly nonlinear dielectric and that from metal can be written as [28]

[TABLE]

where we used $(2|\mathbf{e}|^{4}+|\mathbf{e}^{2}|^{2})/3\approx|\mathbf{e}|^{4}$, $\mathbf{e}$ is the modal electric field, $\chi^{(3)}$ is the third-order nonlinear susceptibility, and “D” and “M” refer to the highly nonlinear dielectric and the metal, respectively. The integrals are taken over the respective media.

The four 2D plasmonic waveguides we consider are shown in Figures 1(a)–(d), and are denoted as “MD” wedge, “MDM” slot, “MDA” DLSPP, and “MDHA” hybrid plasmonic waveguide, respectively. These denotations indicate the type of layers encountered in the direction perpendicular to the metal surface. In addition to “M” and “D” defined above, “A” means air, and “H” is a dielectric that has higher refractive index than “D”. For all these waveguides, we consider only the fundamental plasmonic mode. For the “MD” wedge, we take the radius of the metallic curvature to be $20\,$nm [29].

Throughout this work, we consider $\lambda=1.550\leavevmode\nobreak\ {\rm\mu m}$, and take metal to be silver with linear refractive index $n_{\rm M}=0.1450+11.4125i$ [30], and third-order nonlinear susceptibility $\left|\chi^{(3)}_{\rm M}\right|=2.8\times 10^{-19}\,{\rm m}^{2}/{\rm V}^{2}$ [31]. We consider the highly nonlinear dielectric to be DDMEBT with linear refractive index $n_{0}=1.8$, nonlinear refractive index $n_{2}=1.7\times 10^{-17}\,{\rm m}^{2}/{\rm W}$ [32], and thus $\chi^{(3)}_{\rm\underline{D}}=(4/3)c\varepsilon_{0}n_{0}^{2}n_{2}=1.95\times 10^{-19}\,{\rm m}^{2}/{\rm V}^{2}$. The field profiles of the fundamental modes are calculated using the finite element method as implemented in COMSOL, and are shown in Figure 1.

Figure 2 shows that for all the four plasmonic waveguides under study, $r_{\rm\underline{D}2M}\geq 2\times 10^{3}$. This means that the nonlinear contribution from metal is negligible compared to that from the highly nonlinear dielectric. This is because although $\chi^{(3)}_{\rm\underline{D}}<\left|\chi^{(3)}_{\rm Ag}\right|$, we find $\int_{\rm\underline{D}}|\mathbf{e}|^{4}\mathrm{d}A\gg\int_{\rm M}|\mathbf{e}|^{4}\mathrm{d}A$ since the electric field is mainly confined to the highly nonlinear dielectric region (see Figure 1(a)–(d)). This conclusion does not change even if silver is replaced by gold, which, according to a recent report [26] has a nonlinear susceptibility as large as $\left|\chi^{(3)}_{\rm Au}\right|=9.1\times 10^{-19}\,{\rm m}^{2}/{\rm V}^{2}$.

3 Characteristic measures for plasmonic nanolasers and DFWM devices

The theory we developed [22] shows that plasmonic waveguides used as nanolasers or as DFWM devices can be measured using a set of two characteristic dimensionless parameters. Surprisingly, these two sets of parameters are very similar, and each provides clear physical insights. We discuss these now.

The first characteristic measure for plasmonic waveguides that applies to both devices, is the effective area $A_{\rm eff}$ normalized by the diffraction-limited area $A_{0}\equiv(\lambda/(2n_{0}))^{2}$, i.e., $A_{\rm eff}/A_{0}$. Here $A_{\rm eff}$ is defined as [22]

[TABLE]

where $\{\mathbf{e},\mathbf{h}\}$ are the modal fields, ${\rm max}\{\cdot\}_{\rm\underline{D}}$ indicates the maximum over the active medium which has linear refractive index $n_{0}$. The active medium D denotes the gain medium in nanolasers, and the highly nonlinear layer in DFWM devices.

In nanolasers, $A_{\rm eff}$ measures the maximum Purcell factor of a waveguide mode [33], $F_{\rm m,max}$, through $F_{\rm m,max}=(3/\pi)A_{0}/A_{\rm eff}$, whereas in DFWM devices, $A_{\rm eff}$ measures the driving power $P_{0,{\rm max}}$ for which the nonlinear refractive index change of the highly nonlinear material saturates and reaches its maximum value $\Delta n_{\rm max}$, via $P_{0,{\rm max}}=A_{\rm eff}I_{\rm bulk,max}$. Here $I_{\rm bulk,max}$ is the maximum allowed intensity under plane wave illumination before optical damage occurs.

The second dimensionless measure for plasmonic waveguides in a nanolaser is

[TABLE]

where $g_{\rm th}$ is the threshold gain, i.e., the material gain needed to compensate the modal loss and to achieve lasing. Further, $k_{0}=2\pi/\lambda$, and $\varepsilon_{\rm r,M}=\varepsilon^{\prime}_{\rm r,M}+i\varepsilon^{\prime\prime}_{\rm r,M}$ is the relative permittivity of metal.

The second dimensionless measure for plasmonic waveguide in a DFWM device is

[TABLE]

where ${\cal F}$ is the Figure of Merit, and $U\equiv|\mathbf{e}|^{2}/{\rm max}\{|\mathbf{e}|^{2}\}$ expresses the field uniformity in the highly nonlinear medium.

A low-threshold waveguide-based nanolaser requires (i) enhanced spontaneous emission that is efficiently coupled into the guided mode, which is equivalent to requiring the mode to have a high Purcell factor in order to make full use of the gain medium; and (ii) a low threshold gain $g_{\rm th}$. Optimizing a nanoscale DFWM device requires (i) that the maximum nonlinearity of the material, i.e., the maximum achievable nonlinear index change is reached at a low driving power, $P_{0,{\rm max}}$; and (ii) a large Figure of Merit ${\cal F}$ in order to achieve a high conversion efficiency. We showed earlier that the maximum achievable conversion efficiency is $\eta_{\rm max}=4{\cal F}^{2}/27$ [34], which is defined as the ratio of the output idler power to the input signal power. Requirements (i) for both types of devices are equivalent to a small effective area $A_{\rm eff}/A_{0}$, whereas requirements (ii) for both devices correspond to large $k_{0}/g_{\rm th}$ and large ${\cal F}/{\Delta n_{\rm max}}$.

We can further write

[TABLE]

Here the gain confinement factor $\Gamma_{\rm G}$ [35] and the nonlinear effectiveness ${\rm EFF}_{\rm NL}$ [28] can be written as

[TABLE]

where $Z_{0}$ is the vacuum impedance, and $f_{\ell}=2/(3\ln{(3)})$ accounts for modal loss [34].

Comparing Eqs (3) and (4), or Eqs (7), we find that the only difference between $k_{0}/g_{\rm th}$ or $\Gamma_{\rm G}$ for plasmonic nanolasers and ${\cal F}/{\Delta n_{\rm max}}$ or ${\rm EFF}_{\rm NL}/f_{\ell}$ for plasmonic DFWM devices is that the latter has an additional factor $U$ in the integral kernel of the numerator. This factor highlights the importance of field uniformity in the active medium for nonlinear DFWM applications.

Compared with the parameters $A_{\rm eff}/A_{0}$ and $L_{\rm att}$ for characterizing plasmonic waveguides for pure waveguiding applications, it is clear that the measures for nanolasing and DFWM, $A_{\rm eff}/A_{0}$ and $k_{0}/g_{\rm th}$ or ${\cal F}/{\Delta n_{\rm max}}$, also include $\Gamma_{\rm G}$ or ${\rm EFF}_{\rm NL}/f_{\ell}$, corresponding to the overlap between the mode field and the gain or nonlinear medium. Therefore, plasmonic waveguides that are optimized for waveguiding applications may be suboptimal for nanolasing or DFWM.

4 Characteristics for nanolasing and DFWM

Figure 3 compares the nanolasing and DFWM characteristics of the four 2D plasmonic waveguides from Figure 1(a)–(d). We now discuss these in turn. At this stage we neglect the quenching effects in plasmonic nanolasers; we shall consider these in Section 5.

4.1 MD wedge waveguide

Figure 3(a) shows that, as the wedge angle $\theta$ for MD wedge waveguides increases from $5^{\circ}$ to $120^{\circ}$, $A_{\rm eff}/A_{0}$ increases from $0.035$ to $0.335$, thus remaining below the diffraction limit ($A_{\rm eff}/A_{0}<1$). Simultaneously, $k_{0}/g_{\rm th}$ increases from $150$ to $657$. In other words, the maximum modal Purcell factor $F_{\rm m,max}$ decreases from $27$ to $3$, while the gain threshold $g_{\rm th}$ decreases from $270\,{\rm cm}^{-1}$ to $62\,{\rm cm}^{-1}$. On the other hand, ${\cal F}/{\Delta n_{\rm max}}$ increases only slightly and remains below $30$. Thus the maximum achievable DFWM conversion efficiency $\eta_{\rm max}$ would be low: $\eta_{\rm max}<(400/3)(\Delta n_{\rm max})^{2}=0.33\%$ for $\Delta n_{\rm max}=0.5\%$, making MD wedge waveguides less suited for DFWM applications.

The large difference between the nanolasing characteristics ($k_{0}/g_{\rm th}$ versus $A_{\rm eff}/A_{0}$) and the DFWM characteristics (${\cal F}/{\Delta n_{\rm max}}$ versus $A_{\rm eff}/A_{0}$) can be explained by (3). Although the MD wedge waveguide combines long propagation length $L_{\rm att}$ and small effective area $A_{\rm eff}/A_{0}$, as shown by Figure 4(a), making it particularly attractive for waveguiding applications, the gain confinement factor $\Gamma_{\rm G}$ is large whereas the nonlinear effectiveness ${\rm EFF}_{\rm NL}$ is very low, as shown by Figure 4(b). This difference comes from the $U$ factor, which is important for MD wedge waveguides because of the extremely nonuniform modal electric field of the mode (see Figure 1(a)).

4.2 MDM slot waveguide

Figure 3(b) compares the nanolasing and DFWM characteristics of MDM slot waveguides. $A_{\rm eff}/A_{0}$ can reach the deep subdiffraction regime ($\sim 1\times 10^{-3}$) for narrow slot widths of $w=10$ nm, making the MDM slot waveguide particularly attractive for achieving ultra-high modal Purcell factor in nanolasing applications or ultra-low driving power in nonlinear DFWM applications. For narrow slots $k_{0}/g_{\rm th}$ and ${\cal F}/{\Delta n_{\rm max}}$ have close values. As the slot becomes wider, $k_{0}/g_{\rm th}$ increases whereas ${\cal F}/{\Delta n_{\rm max}}$ decreases in general.

These characteristics can also be understood using (3), where $L_{\rm att}$, $\Gamma_{\rm G}$ and ${\rm EFF}_{\rm NL}/f_{\ell}$ are shown in Figures 4(c)(d). As the slot width increases, more electric energy is confined to the slot, resulting in less attenuation loss or longer $L_{\rm att}$. However, the gain confinement factor $\Gamma_{\rm G}$ and the nonlinear effectiveness ${\rm EFF}_{\rm NL}$ both decrease unexpectedly; and more strikingly, they can be much larger than unity for narrow slot widths, indicating that the modal gain/nonlinearity is much larger than the gain/nonlinearity of the bulk material. This behaviour arises from slow light effects. We have shown [28] that both $\Gamma_{\rm G}$ and ${\rm EFF}_{\rm NL}/f_{\ell}$ can be factorized in an intuitive way as the product of factors that express the effects of losses, of slow light, and of the electric energy overlap with the active medium,

[TABLE]

Here $S\equiv v_{\rm pw}/v_{\rm e}$ is the slow-down factor, with $v_{\rm pw}\equiv c/n_{0}$ the velocity of a plane wave in bulk and $v_{\rm e}=P/W$ the energy velocity of the mode. Here the mode power $P=\frac{1}{2}\int_{\infty}(\mathbf{e}\times\mathbf{h}^{*})\cdot\hat{z}\mathrm{d}A$, and the mode energy $W=\frac{1}{4}\int_{\infty}\left[\mu_{0}|\mathbf{h}|^{2}+\varepsilon_{0}[\partial(\omega\varepsilon^{\prime}_{\rm r})/\partial\omega]_{\omega_{0}}|\mathbf{e}|^{2}\right]\mathrm{d}A$. $\overline{\rho}_{\rm G}$ ($0\leq\overline{\rho}_{\rm G}\leq 1$) quantifies the fraction of electric energy in the gain material, and $\overline{\rho}_{\rm NL}$ ($0\leq\overline{\rho}_{\rm NL}\leq 1$) quantifies fraction of the electric energy ratio in the nonlinear material, weighted by the field uniformity ($|\mathbf{e}|^{2}/{\rm max}\{|\mathbf{e}|^{2}\}$). Here

[TABLE]

where the approximations originate from $2W\approx\int_{\infty}[\partial(\omega\varepsilon^{\prime}_{\rm r})/\partial\omega]_{\omega_{0}}|\mathbf{e}|^{2}\mathrm{d}A$. For metal, which is dispersive, $[\partial(\omega\varepsilon^{\prime}_{\rm r,M})/\partial\omega]_{\omega_{0}}=\varepsilon^{\prime}_{\rm r,M}+\omega_{0}(\partial\varepsilon^{\prime}_{\rm r,M}/\partial\omega)_{\omega_{0}}$, whereas $[\partial(\omega\varepsilon^{\prime}_{\rm r})/\partial\omega]_{\omega_{0}}=n_{0}^{2}$ for the active material which is taken to be dispersionless. Equations (10) and (11) imply that a large $\overline{\rho}_{\rm G}$ requires strong electric energy confinement in the gain material, whereas a large $\overline{\rho}_{\rm NL}$ requires strong and uniform electric energy confinement in the nonlinear material. For MDM slot waveguides with narrow slots, the slow-down factor $S$ is very large, while the electric energy is strongly and uniformly confined to the slot [28], leading to larger-than-unity $\Gamma_{\rm G}$ and ${\rm EFF}_{\rm NL}/f_{\ell}$. As the slot width increases, $S$ strongly decreases, while $\overline{\rho}_{\rm G}$ and $\overline{\rho}_{\rm NL}$ changes slightly [28], resulting in decreasing $\Gamma_{\rm G}$ and ${\rm EFF}_{\rm NL}/f_{\ell}$.

As the MDM slot becomes deeper, i.e., as $h$ increases, $A_{\rm eff}/A_{0}$, $k_{0}/g_{\rm th}$ and ${\cal F}/{\Delta n_{\rm max}}$ all slightly increase. The argument is similar to that for varying $w$, and is thus not further discussed here.

4.3 MDA DLSPP waveguide

Figure 3(c) shows that for MDA DLSPP waveguides $A_{\rm eff}/A_{0}$ first decreases and then increases as the dielectric height increases. This can be understood following Holmgaard and Bozhevolnyi [36]: decreasing the height of the D layer squeezes the field inside the D region, but this behaviour does not continue indefinitely. Thus $A_{\rm eff}$ decreases until reaching a minimum; decreasing the D height further results in a rapid increase in the field outside the D region, resulting in increasing $A_{\rm eff}$. Similarly, $L_{\rm att}$ first decreases and then increases as the dielectric height increases, as shown by Figure 4(e). However, because of the complicated compromise between $L_{\rm att}$ and $\Gamma_{\rm G}$ (or ${\rm EFF}_{\rm NL}/f_{\ell}$) in Figure 4(f), $k_{0}/g_{\rm th}$ (or ${\cal F}/{\Delta n_{\rm max}}$) increases monotonously with height.

For MDA DLSPP waveguides, the distinct characteristics for waveguiding applications, i.e., $L_{\rm att}$ versus $A_{\rm eff}/A_{0}$, and for nanolasing or DFWM applications, i.e., $\Gamma_{\rm G}$ or ${\rm EFF}_{\rm NL}/f_{\ell}$ versus $A_{\rm eff}/A_{0}$, clearly illustrate that the plasmonic waveguide optimized for waveguiding applications may not be optimal for nanolasing or DFWM. For example, Figure 4(e) shows that MDA DLSPP waveguides with $h=200$ nm generally have better waveguiding performance than those with $h=400$ nm, i.e., have longer propagation lengths under the same confinement. However, because the latter have much larger $\Gamma_{\rm G}$ or ${\rm EFF}_{\rm NL}/f_{\ell}$, the net effect is that MDA waveguides with $h=400$ nm generally have better nanolasing or DFWM performance.

4.4 MDHA hybrid plasmonic waveguide

For MDHA hybrid plasmonic waveguides, Figure 3(d) shows that $k_{0}/g_{\rm th}$, ${\cal F}/{\Delta n_{\rm max}}$ and $A_{\rm eff}/A_{0}$ all increase, with increasing height of the active medium $h$. This is because $L_{\rm att}$ strongly increases, whereas $\Gamma_{\rm G}$ first increases and then saturates, and ${\rm EFF}_{\rm NL}/f_{\ell}$ first increases and then drops, as shown by Figures 4(g) and 4(h), respectively. The differences between $\Gamma_{\rm G}$ and ${\rm EFF}_{\rm NL}/f_{\ell}$, and thus between $k_{0}/g_{\rm th}$ and ${\cal F}/{\Delta n_{\rm max}}$ are relatively small. This is because even for large hight $h$, the electric field in the active medium is relatively uniform, as shown in Figure 1(d).

Similar to MDA DLSPP waveguides, the waveguiding characteristics and the nanolasing or DFWM characteristics for MDHA waveguides may be quite different. For example, Figure 4(g) shows that MDHA waveguides with $h=100$ nm have better waveguiding performance but worse nanolasing or DFWM performance than those with $h=200$ nm. The reason, again, is because the latter have much larger $\Gamma_{\rm G}$ or ${\rm EFF}_{\rm NL}/f_{\ell}$, as shown in Figure 4(h).

4.5 Comparison of four plasmonic waveguides

Having understood the nanolasing and DFWM characteristics of all four plasmonic waveguide structures, we now summarize their performance in Figure 5. Depending on the effective area, we define the deep-, moderate- and near-subdiffraction regions, as indicated by blue, green, and yellow in Figure 5.

Figure 5 shows that MDM slot waveguides perform best in the deep-subdiffraction region since it combines the highest modal Purcell factor in nanolasing, and the lowest driving power in DFWM applications. In the moderate-subdiffraction region, the MD wedge waveguide has the best performance for use as a nanolaser since it combines high modal Purcell factor and low gain threshold (Figure 5(a)), whereas the MDHA hybrid waveguide performs best for use as a DFWM device by combining low driving power and high conversion efficiency (Figure 5(b)). This is broadly consistent with the conclusion reached based on the waveguiding characteristics only [12]. As a example, metallic wedge waveguides with integrated reflectors and precisely placed colloidal quantum dots have been used to achieve high-quality quantum emitters [29]. In the near-subdiffraction region, or when we can tolerate a modest modal Purcell factor for nanolasers or high driving power for DFWM devices, the MDA DLSPP performs best.

In our previous work on 1D plasmonic waveguides [22], the 1D counterpart of the MD wedge waveguide is the 1D MD surface plasmon polariton (SPP) waveguide. Therefore, the pink curve in Figure 5 reduces to a point for its 1D counterpart [22]. Therefore, in the moderate-subdiffraction region the best performing structure for both nanolasing and DFWM applications is the MDHA hybrid plasmonic waveguide. In contrast, in the deep-, and the near-subdiffraction regions the wedge waveguide is the best performing 2D plasmonic waveguide configurations for both nanolasing and for DFWM, in line with our conclusions based on a 1D analysis [22].

5 Effect of buffer layers in plasmonic nanolasers

For plasmonic nanolasers, a buffer layer needs to be included between the metal “M” and the gain material “D” so as to avoid quenching. We recently showed that a high-index buffer outperforms a low-index buffer for 1D plasmonic waveguides, and illustrated the 2D extension of this with a specific MDA 2D structure [22]. We now show that this conclusion also applies to the four 2D plasmonic waveguides studied in this work.

Figure 6 shows that for all four 2D plasmonic waveguides the use of a high-index (indicated as “$\bar{\rm H}$” between “M” and “D”) buffer layer results in a smaller effective area than a low-index (“$\bar{\rm L}$”) buffer layer for the same lasing threshold. The thicknesses of all buffer layers, which can be deposited with atomic layer deposition or sputtering, are taken to be $g=5$ nm. In other words, in plasmonic nanolasers based on the four 2D plasmonic waveguides, a high-index buffer layer outperforms a low-index one. This conclusion is consistent with that for 1D configurations [22]. The reason is very similar: the use of a “$\bar{\rm L}$” buffer confines most of the electric field to the buffer layer (second row of Figure 1), thus greatly decreasing the gain confinement factor, although the loss is also reduced. In contrast, “$\bar{\rm H}$” buffers confine most of the electric field to the gain “D” medium (third row of Figure 1), thus greatly increasing the gain confinement factor, which outweighs the slightly larger loss [22].

6 Conclusions

In conclusions, we have investigated the nanolasing and DFWM characteristics of four 2D plasmonic waveguides. Our analysis shows that the MDM slot waveguide is the best performing structure for achieving the highest modal Purcell factor in nanolasing or the lowest driving power in DFWM applications. The MD wedge waveguide has the best performance for use as a nanolaser by combimning high model Purcell factor and low gain threshold, whereas the MDHA hybrid waveguide performs best for use as a DFWM device for combing low driving power and high conversion efficiency. We also show that for plasmonic nanolasers that require a buffer layer between the metal and the gain medium, high-index buffers outperform conventional low-index buffers. We expect this work will advance the design and understanding of plasmonic waveguides that are used in nanolasers and nanoscale DFWM devices.

Acknowledgments

This work was supported by the School of Physics, the University of Sydney, by the Shenzhen Research Foundation (Grant No. JCYJ20180507182444250), and by the State Key Laboratory of Advanced Optical Communication Systems and Networks, China (No. 2019GZKF2).