Phase Matching via Plasmonic Modal Dispersion for Third Harmonic Generation.

The overall effectiveness of nonlinear optical processes along extended nonlinear media highly depends on the fulfillment of the phase-matching condition for pump and generated fields. This is traditionally accomplished by exploiting the birefringence of nonlinear crystals requiring long interaction lengths (cm-scale). For nonbirefringent media and integrated photonic devices, modal phase matching can compensate the index mismatch. Here, the various interacting waves propagate in transverse modes with appropriate phase velocities, but they suffer from a low refractive index contrast and cm-scale interaction lengths. This work harnesses modal phase matching for third-harmonic generation (THG) in plasmonic waveguides using an organic polymer (poly[3-hexylthiophene-2,5-diyl]) as the nonlinear medium. One demonstrates experimentally an effective interaction area as small as ≈ 0.11 µm2 and the phase-matched modal dispersion results in THG efficiency as high as ≈ 10-3 W-2 within an effective length scale of ≈ 4.3 µm. THG also shows a strong correlation with the polarization of the incident laser beam, corresponding to the excitation of the antisymmetric plasmonic modes, corroborating that plasmonic modal phase matching is achieved. This large reduction in device area of orders of magnitude is interesting for various applications where space is critical (e.g., device integration or on-chip applications).

Introduction

Nonlinear optical processes have garnered a renewed interest in a multitude of applications ranging from all‐optical signal processing,[ , ] sensing,[ ] multiphoton microscopy,[ ] photocatalysis,[ , ] to quantum optics.[ ] Recently, optical frequency conversion, third‐harmonic generation (THG) in particular, has gained extensive research significance in spectroscopy[ ] and biological and medical applications,[ , , ] along with the demand for telecom to visible wavelength conversion capabilities.[ , ] Compared to conventional optical methods, an integrated micrometer‐scale device for nonlinear conversion, however, demands a major improvement in the design due to the inevitable long interaction length requirement between the pump and generated fields. In harmonic generation processes, these fields are typically separated by hundreds of nanometers, and an in‐phase build‐up of the field amplitudes is usually achieved by birefringent nonlinear crystals at a centimeter scale.[ ] Therefore, the potential for high‐density integration for on‐chip nanoscale optical conversions is compromised in these methods. Here we demonstrate phase‐matched modal dispersion in plasmonic waveguides as an alternative solution to improve the coherence length l coh = π/Δk, by minimizing the momentum mismatch Δk between the interacting fields, and consequently the interaction length, to attain good nonlinear conversion efficiencies for THG. We show that the material's dispersion is compensated by a judiciously chosen geometry to match the phase velocity of the pump and generated fields, allowing us to demonstrate an effective interaction area as small as ≈ 0.11 µm2 and THG efficiency as high as ≈ 10–3 W‐2 which is comparable to other devices based on nonlinear optical polymers inside short waveguides (≈ 4.3 µm long). This is an order‐of‐magnitude reduction in device footprint compared to photonic devices[ , ] and is potentially useful for applications with stringent space requirements and device integration.

Plasmonics offers a versatile route to downscale the nonlinear components for the realization of active photonic circuitry with device designs in subwavelength dimensions.[ , ] The field confinements and the associated large field amplitudes allow for coherent enhancement of inherently weak nonlinear light–matter interaction in nanoscale systems.[ , , ] While plasmonic metals, in particular gold, exhibit very high third‐order susceptibilities,[ ] their nonlinear response is generally limited by the intrinsic electron damping and electromagnetic screening, especially in the visible and near‐infrared energy ranges.[ ] Plasmon‐mediated nonlinear conversion has been achieved in the past in 2D material systems[ ] and colloidal nanorods,[ , , , ] but at the cost of active interaction volume or limited flexibility in the device design. In comparison, nonlinear optical polymers (NOPs) offer solution‐based processing and large nonlinear susceptibilities via chemical modification or changing the structure of the polymer films. Their fabrication protocols are further suited for integrated nonlinear devices,[ , , , ] such as silicon photonics structures for applications in nonlinear optics, plasmonic phase modulation,[ ] and all‐optical signal processing.[ , ] Recently, four‐wave mixing was demonstrated using NOP as the nonlinear medium in a plasmonic slot waveguide.[ ] In this case, however, phase matching was not a concern as l coh, for small Δk, was much larger than the waveguide length.

In this study, we demonstrate THG in plasmonic waveguides in a two‐wire transmission line (TWTL) configuration,[ ] where two Ag slabs are placed in nanoscale proximity with a NOP cladding as the active medium, as illustrated in Figure . This type of platform supports symmetric (S) and antisymmetric (AS) plasmonic modes[ ] and is in favor of modal phase matching (MPM),[ , ] with operation wavelengths extending from the optical to the near‐infrared spectral range. We used poly(3‐hexylthiophene‐2,5‐diyl) (P3HT) as NOP for the conversion of the telecommunication wavelength (1550 nm) to its third harmonic (TH) (517 nm). This NOP has a large third‐order nonlinear susceptibility χ (3) ≈ 10−19 m2V−2.[ ] We fabricated a series of devices with different waveguide lengths (L) and gap widths (W), carefully designed to achieve MPM with a good overlap between the two AS eigenmodes: one at the fundamental wavelength (FW) 1550 nm and the other at 517 nm. By focusing a 1550 nm wavelength pump laser on one of the antennas (Antenna‐in, Figure 1), we found a significant THG output at 517 nm at the other end of the waveguide (Antenna‐out, Figure 1). The polarization angle (θ in) of the incident laser beam shows a strong correlation with the THG signal intensity, an observation consistent with the excitation of the AS plasmonic modes in nonlinear energy conversion and further verified with the nonlinear finite difference time domain (FDTD) simulations. The results presented here pave the way to a nanophotonic device element for telecom‐to‐visible nonlinear optical conversion, but are also useful in other areas where space is limited.

Figure 1

Schematic illustration of the device consisting of a TWTL plasmonic waveguide and two plasmonic antennas. The incident laser polarization angle θ in is defined with respect to the orientation of the waveguide. The inset shows the cross‐section of the waveguide with P3HT layer as the nonlinear medium.

Results and Discussion

Waveguide Design and Modal Phase Matching

The plasmonic waveguides investigated here consist of two Ag slabs in TWTL configuration (Figure 1 inset) on SiO2 substrate (285 nm thick on a Si wafer), with the width and height of the slabs fixed as 400 and 80 nm respectively. The gap width W is varied from 100 to 200 nm (see Section S1, Supporting Information). For the slabs, Ag is preferred as the plasmonic material, which has a lower interband absorption loss around the THG wavelength (517 nm in this work) than Au.[ ] The NOP (here P3HT) was spin‐coated over the structure to a thickness of around 200 nm, filling the gap between the two slabs (see Figure S1, Supporting Information for AFM images). In the TWTL configuration, the electric field inside the slot is dominated by the AS modes, whereas the S modes (not being investigated in this work) are delocalized outside the Ag slabs.[ ] Therefore, we consider only the first‐order AS mode for FW and the second‐order AS mode for THG, which are depicted in Figure as AS1, FW (top panel) and AS2, TH (middle panel), respectively. The eigenmode analysis and mode dispersion calculations were performed using MODE‐solutions of the commercial FDTD package Lumerical.[ ] The mode orders are designated according to their charge distribution in the Ag slabs (see Figure S3, Supporting Information for the charge distribution plot). We note that the field intensity inside the Ag slabs is suppressed by the electromagnetic screening, leading to a negligible contribution to the THG signal. For the THG process, the effective interaction area between the AS1, FW mode and the P3HT medium is crucial and is quantified as the ratio of the integrated electromagnetic power in the transverse plane (xy‐plane) to the effective intensity in the nonlinear medium (P3HT):[ ] where Z 0 is the free space impedance, is the unit vector orthogonal to the xy‐plane and n P3HT is the linear refractive index of the NOP. Plasmonic mode confinement warrants A eff as small as ≈ 0.11 µm2 for W = 100 nm, as estimated from the FDTD calculations (A eff for all values of W are listed in Table S1, Supporting Information). Good overlap between the two AS modes inside the waveguide gap (Figure 2a, bottom panel) enables efficient coupling to the AS2, TH mode of THG photons radiated by the nonlinear polarization created in the P3HT by the FW field propagating in the AS1, FW mode.

Figure 2

Eigenmode analysis and the mode excitation for the waveguide segment. a) Modal field profile (E x) of the eigenmodes supported by the waveguide for AS1, FW (top) and AS2, TH (mid). The bottom panel shows the E x obtained from the top and middle panels at y = 0 (middle of the waveguide) plotted as a function x. b) Real part of the effective refractive index (n eff) for the eigenmodes AS1, FW (bottom x‐axis) and AS2, TH (top x‐axis) as a function of wavelength. c) Simulated electric field (E x) in the waveguide excited by a Gaussian beam with θ in = 90°. d) Polar plot for the simulated in‐coupled power (normalized to the maxima) for λ = 1550 nm as a function of the laser polarization angle θ in, which follows a sin2 θ in with a maximal mode power coupled at θ in = 90° and 270°.

Along with this mode overlap, the efficiency of the THG process depends critically on phase matching between FW and TH eigenmodes and we turn now to the discussion of the MPM mechanism. Important results obtained from the FDTD calculation are the complex effective indices for the two AS modes with α standing for FW and TH. Figure 2b shows as a function of wavelength for the AS1, FW mode (bottom x‐axis) and the AS2, TH mode (top x‐axis) for the case of W = 100 nm. The phase mismatch between the FW and the TH eigenmodes is defined as Δk = 3k FW − k TH, where is the wavenumber of the α‐eigenmode, with λ the corresponding wavelength. As for THG λ FW = 3λ TH, then , showing that the MPM condition ΔkL ≪ 1 is attained when the two modes have the same effective index, where L represents the length of the waveguide. Figure 2b shows that the MPM is achieved for the AS1, FW and the AS2, TH at the crossing point of the dispersion curves at around 1600 nm. In this study, the AS1, FW and AS2, TH modes wavelengths are 1550 and 517 nm and the difference between the effective indices is about ≈ 0.015, which corresponds to l coh ≈ 17 µm, for W = 100 nm.

To couple the pump beam into the AS modes from free space and to couple out the THG signal from the waveguide end, we use dipole‐like plasmonic antennas (Antenna‐in and Antenna‐out respectively, Figure 1) oriented perpendicularly to the waveguide. Such an antenna structure offers two advantages: i) it has a polarization angle θ in (defined with respect to the waveguide orientation in Figure 1) dependent mode selectivity for the excitation of AS or S modes,[ ] and ii) it allows for nearly vertical outcoupling for the guided THG modes, which can be captured by an objective lens in the reflection mode (see Figure S4d, Supporting Information). A linearly polarized pump beam is introduced through Antenna‐in at normal incidence. For θ in = 90°, the coupling to the AS modes is maximal and a propagating plasmonic field is established inside the waveguide as shown in Figure 2c, where E is the x component of the complex electric field amplitude evaluated at 1550 nm. The polar plot in Figure 2d shows the electromagnetic power coupled into the waveguide (normalized to the maximum at θ in = 90°) for λ = 1550 nm calculated as a function of θ in; it follows a sin2 θ in dependence with two maxima at 90○ and 270○, which corresponds to the direction colinear with the antenna dipole.

Nonlinear Optical Measurements

Figure shows the scanning electron microscope (SEM) image of a representative device with L = 4 µm and W = 100 nm, and the Antenna‐in is located at z = 0. The devices were fabricated on a SiO2‐Si substrate with standard electron‐beam lithography and evaporation techniques and the P3HT polymer film was spin‐coated on the waveguide (See Experimental Section for details). The schematic illustration of the experimental setup for the nonlinear optical measurements is shown in Figure 3b. To obtain the spatially resolved nonconfocal mapping for the THG signals, we focus the laser beam (1550 nm, FW) on the Antenna‐in through an objective lens while scanning the THG signal collection spot along the waveguide (i.e., z‐axis). The acquisition of signal collection location is realized using a Y‐fiber equipped with a laser diode as an indicator. The input 1550 nm pump laser intensity (P in) is monitored and controlled by the combination of a half‐wave plate, a polarizing beam splitter, and a power meter. The angle of excitation polarization (θ in) is controlled by another half‐wave plate before the objective (see Experimental Section for details).

Figure 3

Spectroscopic THG measurements. a) SEM image for the representative device with a waveguide length L = 4.0 µm and a gap width W = 100 nm. b) Schematic illustration of the optical setup for recording the THG spectra and the images. The nonconfocal THG spectroscopy is realized by moving the fiber coupler on the image plane using a piezo‐controlled Y‐fiber detection system. The injection of a 532 nm CW laser in one of the ends of the Y‐fiber allows precise determination of the location where the signal is acquired through direct visualization of the green spot in the sCMOS camera. Part of the THG signal is routed to a spectrometer by the other Y‐end of the fiber. c) Intensity map of the line‐trace spectra (top). Integrated THG signal as a function of detection position along the waveguide (bottom). d) THG spectra from the Antenna‐out location, corresponding to the detection position at z = 4.5 µm.

The THG intensity (I TH) for θ in = 90°, collected along the length of the waveguide (corresponding to the dashed‐line box shown in Figure 3a) as a function of the wavelength is plotted in Figure 3c (top). In the bottom panel of Figure 3c the spectrally integrated I TH along the waveguide length is plotted to show a maximum intensity at the Antenna‐out location. Figure 3d shows the wavelength resolved spectra of the I TH, collected at z = 4.5 µm (Antenna‐out) and centered at 517 nm. Under the laser excitation of 1550 nm with θ in = 90°, the AS1, FW mode of the waveguide is excited and the third‐order P3HT isotropic susceptibility induces a nonlinear polarization[ ] in the waveguide transverse direction (where ε 0 is the permittivity of free space and E(ω) is the x component of the complex amplitude of the electric field of the AS1, FW mode at frequency ω), which generates the AS2, TH mode, phase‐matched with the AS1, FW mode, thus resulting in the I TH at the Antenna‐out centered at the wavelength 517 nm.

The third‐order character of the nonlinear process is further confirmed by measuring the pump intensity (P in) dependence of the integrated output THG power (P out) as shown in Figure : the power‐law fit to the experimental data yields β ≈ 3 as expected.[ ] For P in in the range of 0.22–1.1 mW, we observe P out ranging between 11 and 1300 fW from the Antenna‐out location. This gives an average THG power conversion efficiency η (P out/(P in)3) ≈ 10–3 W‐2 (see Experimental Section for details). We restrict the unit of η to W–2 instead of W–2 m–1 to avoid inconsistent normalization of η with L, as the plasmonic losses make the FW power along the waveguide an attenuating function of L. Figure 4b shows the angular dependence of P out as a function of θ in. The minimum in the polar plot corresponds to P out for polarization along the z‐axis ( θ in = 0° or 180°), where only the S modes of the plasmonic waveguide are excited.[ ] For θ in = 90° or 270°, the AS1, FW mode (Figure 2a, top) is efficiently excited, with major field intensity localized within the waveguide gap, leading to the conversion to the AS2, TH mode (Figure 2a, middle) from the nonlinear interaction with the P3HT medium. Therefore, we ascribe this strong dependence of the laser polarization as a direct evidence of the plasmonic MPM in our waveguide. The third order character of P out with the variation in θ in is further evidenced from the fitting (solid line), in contrast with the sin2 θ in∝|E x|2 dependence (dashed line) which would be expected for a linear response. Here E x represents the x‐component of the incident laser field.

Figure 4

a) Output THG power (P out) from the Antenna‐out as a function of average P in. The black solid line is a power‐law fit . b) P out as a function of θ in for P in = 0.6 mW. The black solid line is (sin3 θ in)2 function fitted to the experimental data and the gray dashed curve shows sin2 θ in function. The error bars show the standard deviation of five repeated measurements on the same device. c) Nonlinear FDTD simulation results for the power (P FDTD) at the Antenna‐out location, evaluated as a function of the θ in and λ. P FDTD is normalized to the maximum value at 1550 nm. d) Comparison of the output power from the experiment (P out, top) and the nonlinear FDTD simulations (P FDTD, bottom) as a function of θ in. Black solid and dashed lines are the fits to (sin3 θ in)2 and sin2 θ in functions respectively.

The plasmonic MPM efficacy was also checked by nonlinear FDTD simulations. A Gaussian beam was used as the input source for the FW excitation (1550 nm) with a pulse duration of 100 fs and a field peak amplitude of | | ≈ 107 V m‐1, such that the third‐order term (χ (3)| |2 ≈ 10–5) is a perturbation to the linear response (see Experimental and Simulation Section for more details). Figure 4c shows the power outcoupled (P FDTD) from Antenna‐out as a function of the wavelength (x‐axis) and θ in (y‐axis) for the waveguide with L ≈ 4 µm and W = 100 nm. P FDTD is normalized to the maximum value at 1550 nm. The maxima observed at θ in = 90° and 270° for wavelengths 1550 and 517 nm respectively correspond to the outcoupling of the FW and TH. To compare with the experimental data, P FDTD values along the dotted lines in Figure 4c (at 1550 and 517 nm and normalized to unity individually) are plotted as a function of θ in in the bottom panel of Figure 4d. The third order nature of the THG process due to the antisymmetric mode AS1, FW is evident from the (sin3 θ in)2 fitting (solid line). The peak positions (at θ in = 90° and 270°) correspond to the excitation of the AS modes in the plasmonic waveguide, consistent with the experimental result (blue spheres, top panel). The variation of P FDTD at FW (red circles, bottom panel) follows the projection of the laser polarization on the antenna direction in a sin2 θ in form (dashed line, bottom panel), consistent with the FW mode coupling (Figure 2d).

Nonlinear Wave Equation Model

The effect of the plasmonic MPM in the growth of the THG wave amplitude in its propagation along the waveguide is theoretically modeled by solving the nonlinear Helmholtz equation assuming the validity of the slowly varying amplitude approximation.[ ] In this case, the second‐order derivative with respect to z, the coordinate in the propagation direction, of the THG electric field complex amplitude A TH(z) is neglected, and the equation is simplified to: where A FW(z) is the FW field amplitude, κ is a proportionality constant, is the TH propagation length, and Δk is the phase mismatch between FW and THG. Neglecting the depletion of the FW pump due to the parametric process but keeping linear absorption, the spatial variation of A FW(z) is simply with as the FW propagation length and as the pump amplitude at z = 0 (Antenna‐in, Figure 2a). Solving Equation 2 with boundary condition A TH (z = 0) = 0 gives the detected signal, proportional to S (L, Δk) = |A TH(z = L)|2, as: with L as the waveguide length. Within the limit of lossless propagation, Equation 4 reduces to the familiar phase mismatch expression L 2sinc2[Re(Δk) L/2].[ ]

To compare the variation of S(L, Δk) with L and Δk, we measured the THG conversion efficiency η for devices with L ranging from 3 to 5.5 µm, with increments of 0.5 µm, with three representative gap widths W = 100, 150, and 200 nm. In general, as W is increased from 100 to 200 nm, the phase mismatch Δk between AS1, FW and AS2, TH becomes larger, compromising the spatial overlap of the modes and diminishing the THG efficiency (see Section S4, Supporting Information for the effective index and phase matching plots). Figure shows η (spheres) plotted along with the theoretical results (solid lines) obtained by evaluating Equation 4. Individual theoretical curves are rescaled to match the maxima of the experimental data set corresponding to each W. Good qualitative agreement is observed between the theory and the experiment. Slight deviations can be attributed to fabrication imperfections. For W = 100 nm case, η increases with L, until a maximum is reached at L = L eff (upward arrows in Figure 5a), the most well‐suited waveguide length for the coherent accumulation of the THG amplitude to compensate for the propagation losses of the combined TH and FW modes. After this point, the observed signal power starts to decrease as the parametric process ceases to dominate over the losses. For the W = 150 nm and 200 nm cases, we also observe similar trends, but the THG conversion peaks are shifted to lower values of L eff with reduced power, because of the higher phase mismatch between the TH and FW modes, resulting in faster decoherence. This behavior with varying W is obvious in the contour map of S(L, Δk) plotted in Figure 5b. The values of Δn ( = Δk × λ TH/2π) reflecting the phase mismatch (in Figure S4a, Supporting Information), obtained from the FDTD mode calculations for 100 < W < 200 nm are used for computing the contour plot. The white dashed line highlights the trend in the THG maximum and indicates the variation of the coherent signal accumulation length when MPM is fulfilled with a decrease in W. By including the propagation and effective index information obtained from the FDTD simulations, the general trend observed in the THG intensity from the theory shows a consistent agreement with the experiment. We note that by decreasing W to below 100 nm the propagation losses start dominating the THG process, which will effectively reduce L eff to smaller values, with lower η. It is envisioned that other structural parameters such as the antenna coupling factors can be optimized to further enhance THG efficiencies. We note that the efficiency of the Antenna‐out used for the present study is limited to ≈ 0.17% for the THG outcoupling at 517 nm (see Section S4, Supporting Information). Nevertheless, for milliwatt input powers (P in), we observe η ≈ 10–3 W‐2 over a length scale (L eff) ≈ 4.3 µm (see Table S1, Supporting Information). Taking into account the antenna efficiencies, it indicates a high internal (within the TWTL waveguide) TH conversion in the plasmonic waveguide before the Antenna‐out, as a result of the MPM achieved here.

Figure 5

a) THG conversion efficiency (η) experimentally evaluated for the waveguides of W = 100 nm, 150 nm, and 200 nm, plotted as a function of L. Upward arrows represent L eff corresponding to the local maxima of η. Error bars represent the standard deviation of measurements of five devices considering fabrication deviations. Solid lines show the theoretical evaluation S(L, Δk) from Equation 4. b) Contour map of S(L, Δk) from Equation 4, with values of Δn and l p from Figure S4a,b (Supporting Information), for 100 < W < 200 nm and 0 < L < 10 µm. The white dashed line indicates the maxima of S(L, Δk) values.

Conclusions

In summary, we have demonstrated MPM for THG using plasmonic waveguides in the TWTL configuration, with P3HT as the nonlinear medium for the wavelength conversion from 1550 to 517 nm. The gap width between the slabs is chosen to maximize the MPM between the AS modes (AS1, FW and AS2, TH) while keeping the incident FW laser beam effectively coupled to the plasmonic modes, resulting in a coherent nonlinear interaction in the waveguide to accumulate the THG signal. As the MPM is strongly correlated with the symmetry of the modes, THG can be directly modulated with the incident laser polarization angle θ in. We observe maxima in the THG signal for θ in= 90○ and 270○ when the AS modes are efficiently excited in the waveguide. The emitted power follows the input power in the third order, consistent with the typical TH response. Numerical calculations performed using nonlinear FDTD and wave‐equation studies corroborate the experimental observations. The case of W = 100 nm with A eff ≈ 0.11 µm2 and L eff ≈ 4.3 µm ensuing a local maximum in η ≈ 10–3 W‐2 (this efficiency is similar to that of other devices based on NOPs[ , ]) demonstrates that the tradeoff between the THG and the plasmonic propagation losses within µm scale interaction length can be mitigated. In other words, based on our analysis we anticipate that the device performance can be further improved for even smaller devices provided Δk is further minimized (see Figure 2b). In addition, the overall efficiency of the devices can be improved by optimizing the impedance matching between guided THG modes and radiation at the Antenna‐out. Our present study provides a solution for nonlinear optical conversion based on compact plasmonic waveguides that enable the tunability of phase‐matching via geometry and polarization control.

Experimental and Simulation Section

Fabrication

Plasmonic waveguides were fabricated on a 285 nm SiO2–Si substrate, using electron beam lithography (EBL) and e‐beam evaporation. The waveguides were defined by EBL (JEOL, JBX‐6300FS), with bilayer PMMA resist (PMMA 495 A3 and PMMA 950 A5, prebaked at 180 °C for 5 min and 2 min respectively), followed by resist development in MIBK:IPA (1:3) for 30 s and IPA for 30 s. After development, the Ag layer (80 nm thick) with a Ti adhesion layer (3 nm thick) was deposited using e‐beam evaporation (AJA, Orion 8E) followed by lift‐off in acetone. The P3HT layer was first created by dissolving 20 mg of regioregular P3HT (average molecular weight MW = 70000–100000 g mole‐1, Sigma‐Aldrich) in 1 mL of chlorobenzene. The solution was then spin‐casted on the samples at 500 rpm for 10 s, then at 1000 rpm for 60 s to create a 200 nm film. The devices were thermally annealed in air at 100 °C for 10 min to stabilize the optical performance.

Numerical Simulations

The eigenmode profiles for the plasmonic waveguides were simulated using the FDTD method (Lumerical[ ] FDTD Solutions, Ansys Canada Ltd.). The optical properties of Ag (80 nm), SiO2 (285 nm), and Si were taken from Palik.[ ] The optical properties of the P3HT cladding (200 nm) were obtained from the ellipsometry measurements as detailed in Section S5 (Supporting Information). For the effective index calculation, a 2D cross‐section of the waveguide in the xy‐plane was used (see Figure 1 inset), with a simulation domain size of 1.5 µm × 0.75 µm (x × y). A minimum mesh size of dx = dy = 2 nm was applied across the simulation area. To solve for the AS modes, antisymmetric boundary conditions were used in the ±x direction, whereas perfectly matched layers were used as the ±y boundary conditions with 12 Standard absorbing layers. The THG response of the device was calculated using the nonlinear FDTD method for a 3D domain (domain size 7 µm × 5 µm × 1.5 µm (x × y × z)). Nonlinear material implementation[ ] of the Lumerical was used to represent the P3HT material. For THG, the polarization P in the small field amplitude limit can be written as . For the convergence of the FDTD model, nonlinear terms have to be small. In other words, should be much smaller than . A value of ≈ 10–19 m2 V‐2 was used for the P3HT film and a field peak amplitude of ≈ 107 V m‐1 for the incident Gaussian beam, which gives ≈ 10–5. This value is much smaller than (typically ≈ 1) and guarantees the convergence of the numerical simulations. Perfectly matched layers were used as the boundary conditions in all three directions and non‐uniform meshing was applied across the simulation domain for the nonlinear FDTD simulations.

Optical Measurements

The laser beam at 1550 nm (FW) with an fwhm spot size of ≈ 1.2 µm and average power of ≈ 10 mW produced from a femtosecond fiber laser system (Toptica, FemtoFiber pro, 1550 nm, 90 fs, 80 MHz) was focused on the Antenna‐in (Figure 1) through an objective (Leica, 100×, NA 0.85). To spatially monitor the THG along the waveguide (i.e., z‐axis) one used a Y‐fiber (Thorlabs, TW560R1F1) mounted on a 3D piezo stage (Piezojena Tritor 100) tandem on a 3D manual mechanical stage. The input port of the fiber was mounted on the microscope to collect the signal focused by a lens (Lens A, Nikon 20×, NA 0.5). The signal port (99%) was connected to the spectrometer, while the tap port (1%) was coupled to a 532 nm laser. The location and area (with a detection spot diameter of ≈ 0.7 µm) of the signal being collected and measured are known by observing the 532 nm laser spot on the sample with the imaging branch of the microscope. By moving the input port of the fiber controllably with the 3D piezo stage, the line‐trace THG spectra were recorded. The spectrometer (Andor, Kymera 328i) equipped with a Vis‐detector (Andor, iDus 401A) was used to measure the spectra. With the visible sCMOS camera (PCO, Edge 4.2) on the other branch, the distribution of light emission intensity was captured. The collection efficiency of the setup to the sCMOS (QE = 0.75, gain = 0.46 e‐ per count[ ]) for 517 nm is estimated to be 16.9%, hence the detection sensitivity (S d) is 1.4×10–3 fW s count–1. THG power (P out) is converted from the time‐averaged intensity (I TH) collected by the sCMOS through P out = S d × I TH.

Conflict of Interest

The authors declare no conflict of interest.

Supporting Information

Click here for additional data file.