Efficient ultrathin all-dielectric nonlinear optical device

Abstract

An apparatus and method are provided for generating harmonic light from a pump beam that is impinged on a metasurface comprising a plurality of all-dielectric resonator bodies. A multiple quantum well structure formed in each resonator body includes asymmetric coupled quantum wells having intersubband transition frequencies that couple to Mie resonances of the resonator bodies.

Claims

1. Apparatus comprising a metasurface in which a plurality of elements is disposed on a substrate surface, wherein: each of the elements is an all-dielectric resonator body conformed to resonate at a first resonant frequency; each of the elements has at least some resonant response at a second resonant frequency that is a multiple of the first resonant frequency; each of the elements includes at least one multiple quantum well (MQW) structure constituted by a sequence of layers that defines multiple quantum wells; each MQW structure comprises at least one coupled quantum well (coupled QW) in which two or more dissimilar quantum wells are disposed sequentially; each of the coupled QWs is dimensioned to have at least a first, a second, and a third electronic eigenstate having respective energy levels; the eigenstate energy levels are separated from each other by transition energies, each corresponding to a respective optical transition frequency; the first resonant frequency is substantially equal to one of the optical transition frequencies; and the second resonant frequency is substantially equal to another of the optical transition frequencies.

2. The apparatus of claim 1, wherein the transition energies are energies of intersubband transitions of the MQW structures.

3. The apparatus of claim 1, wherein the elements are disposed in a periodic pattern on the substrate surface.

4. The apparatus of claim 1, wherein each of the resonator bodies is conformed to resonate at both the first and second resonant frequencies.

5. The apparatus of claim 1, wherein the first, second, and third electronic eigenstates have respective energy levels separated from each other by substantially equal separations in energy.

6. The apparatus of claim 1, wherein each of the resonator bodies comprises III-V semiconductor material.

7. The apparatus of claim 1, wherein each of the MQW structures is an InGaAs/AlInAs heterostructure.

8. The apparatus of claim 1, wherein the first resonant frequency belongs to a fundamental Mie resonance of the all-dielectric resonator bodies.

9. The apparatus of claim 1, wherein the first resonant frequency is a frequency of a magnetic dipole resonance of the all-dielectric resonator bodies.

10. The apparatus of claim 1, wherein the second resonant frequency lies within a portion of the electromagnetic spectrum extending from a vacuum wavelength of 1.5 μm to a vacuum wavelength of 4 μm.

11. The apparatus of claim 1, wherein each of the all-dielectric resonator bodies comprises Group-III antimonide material.

12. The apparatus of claim 1, wherein each of the MQW structures is an AlSb/InAs/AlSb heterostructure.

13. The apparatus of claim 1, wherein each of the coupled QWs has exactly two dissimilar quantum wells.

14. The apparatus of claim 1, wherein each of the coupled QWs has exactly three dissimilar quantum wells.

15. A method for generating harmonic light from a pump beam of light having a pump frequency, comprising: impinging the pump beam on a metasurface comprising a plurality of all-dielectric resonator bodies, such that a Mie resonance is excited in each of the said bodies; and collecting harmonic light from the metasurface, the harmonic light having a frequency that is a multiple of the pump frequency, wherein: the impinging of the pump beam excites a fundamental intersubband transition (IST) of MQW heterostructures contained within the said bodies, wherein the fundamental IST is coupled to a Mie resonance supported by the said bodies; and the impinging of the pump beam further excites at least one higher-order IST of the said MQW heterostructures, wherein the higher-order IST has a transition frequency that overlaps at least one Mie-resonant response of the said bodies having a frequency that is a harmonic of the pump frequency.

16. The method of claim 15, wherein the harmonic of the pump frequency is the second harmonic.

17. The method of claim 15, wherein the harmonic of the pump frequency is the third harmonic.

18. The method of claim 15, wherein the Mie resonance excited by impinging the pump beam is a fundamental Mie resonance of the all-dielectric resonator bodies.

19. The method of claim 15, wherein the Mie resonance excited by impinging the pump beam is a magnetic dipole resonance of the all-dielectric resonator bodies.

20. The method of claim 15, wherein the harmonic light lies within a portion of the electromagnetic spectrum extending from a vacuum wavelength of 1.5 μm to a vacuum wavelength of 4 μm.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a notional drawing, in perspective view, showing an example metasurface in operation for generating second-harmonic light.

(2) FIG. 2 is a perspective view of a single element of a metasurface in an example embodiment.

(3) FIG. 3 is an energy-band diagram for a quantum heterostructure used in the elements of an example metasurface. FIG. 3 also shows three wavefunctions superimposed on the energy-band diagram.

(4) FIG. 4 presents a graph of the numerically calculated value of the intersubband nonlinear susceptibility of a multi-QW structure as a function of the pump wavelength. Also presented on the same graph is an experimentally measured intersubband absorption spectrum of multi-QW heterostructures discussed here. An inset shows the waveguide configuration used for the absorption measurements with the input direction k indicated.

(5) FIG. 5 is an intensity map, in the xz-plane, of numerically calculated values of |E.sub.z/E.sub.inc| of the fundamental magnetic dipole resonance within a resonator body in an element of an example metasurface.

(6) FIG. 6 shows three plots of an experimentally measured linear reflectance spectrum. From left to right, these spectra were measured on resonators with different aspect ratios due to different radii of the cylinders. The reflectance peak at the MD resonance is labeled on each of the three spectral curves and scales to longer wavelength with increasing radius.

(7) FIG. 7 is a flowchart in which cross-sectional cartoon views illustrate several sequential stages in a fabrication process for metasurfaces of the kind described here.

(8) FIG. 8 is a graph of experimentally measured optical power versus pump wavelength, for the second-harmonic signal reflected from three respective Mie metasurfaces.

(9) FIG. 9 is an experimentally measured linear reflectance spectrum of a metasurface with cylinders of R=1.4 μm. Rabi splitting of the MD resonance can be seen in the figure.

(10) FIG. 10 is a numerically calculated graph of susceptibility plotted versus wavelength for an example metasurface of the kind described here.

(11) FIG. 11 is a numerically calculated graph of field enhancement plotted versus wavelength for an example metasurface of the kind described here.

(12) FIG. 12 is a graph of experimentally measured reflected SH power conversion efficiency as a function of pump wavelength. The measurements were for normal pump incidence on a metasurface with cylinders of R=1.4 μm.

(13) FIG. 13 is a graph illustrating the performance of an example metasurface. Reflected optical power in the second harmonic is plotted versus the square of the pump optical power. The data in the figure were from experimental measurements.

(14) FIG. 14 is another graph illustrating the performance of an example metasurface. Experimentally measured conversion efficiency is plotted versus pump power.

(15) FIG. 15 is a graph of conversion efficiency, in normalized units, numerically calculated for normal incidence on a metasurface as a function of pump wavelength. Several finite sizes of arrays of resonator elements are represented as separate plots in the figure.

(16) FIG. 16 is an intensity map in which the quantities plotted on the respective axes are pump angle of incidence and pump wavelength, and the intensity represents the SH power conversion efficiency.

(17) FIG. 17 is an energy-level diagram of an InGaAs/AlInAs coupled-well heterostructure with three equally spaced transitions, which can simultaneously generate both second-harmonic and third harmonic signals. Wavefunctions for the four electronic states are notionally represented in the figure, with vertical separations corresponding to their separations in energy.

(18) FIG. 18 is an energy-level diagram of an AlSb/InAs coupled well heterostructure. Wavefunctions for the lowest three electronic states are notionally represented in the figure, with vertical separations corresponding to their separations in energy. This structure can generate second harmonic light at 1.5 μm.

DETAILED DESCRIPTION

(19) In embodiments of the invention, multi-quantum-well (MQW) semiconductor heterostructures are embedded within all-dielectric Mie resonators in a metasurface. The embedded MQW heterostructures are designed to support intersubband transitions (ISTs). Accordingly, the physical response of each of the resulting resonator structures has two coupled components: A spectral component associated with the ISTs, and an optical, or photonic, component associated with the Mie resonances.

(20) A “Mie resonator”, in this regard, is an optical resonator having at least one spatial dimension that is comparable to or smaller than the vacuum wavelength of light at some fundamental resonant optical frequency. “Mie resonance” is a general term for an optical resonance that is supported by the resonator at such a fundamental frequency or at one of its higher harmonics that correspond to shorter wavelengths. These resonances are what we refer to here as Mie resonances or “Mie modes”, where the lowest order (smallest frequency) is the magnetic dipole Mie mode, and second-order modes are the quadrupole modes.

(21) In operation, polaritonic coupling of the ISTs to the Mie modes can result in high conversion factors for second-harmonic generation, and in some cases it can also produce third-harmonic generation and other higher-order nonlinear effects. In examples, conversion factors of 0.5 mW/W.sup.2 or even more may be achievable using the lowest-order magnetic dipole Mie mode of the resonator, or using the second-order magnetic quadrupole Mie mode.

(22) The metasurface is operated by impinging it with a pump beam. No special polarization of the pump beam is required. In particular, any polarization of the electric field along the x- or y-direction of the pump beam can be effective. Moreover, the pump beam can be directed either normally or off-normally to the metasurface. It should be noted, however, that normal emission of second-harmonic light is suppressed when the pump beam is at normal incidence. That phenomenon is discussed below.

(23) The design of the MQW heterostructure offers some control over the magnitude and sign of the second-order nonlinear susceptibility elements χ.sup.(2)(2ω). This degree of freedom, when combined with the ability to control the polarization and spatial mode profiles of all-dielectric Mie resonators, can be leveraged to further enhance the SHG efficiency using techniques similar to quasi-phase-matching.

(24) Our experimental demonstration was at mid-infrared wavelengths, for which the example embodiments described below were designed. However, the same approach for SHG can be scaled from near-IR to longwave-IR wavelengths, using suitable combinations of material systems and quantum well designs.

(25) In embodiments, the resonator bodies are fabricated from III-V semiconductors, and they include n-doped quantum wells (QWs). Each resonator body is conformed as a Mie resonator.

(26) For example, we fabricated and tested metasurfaces in which each resonator body comprises an InGaAs/AlInAs quantum heterostructure.

(27) FIG. 1 is a notional drawing in perspective view, showing an example metasurface in operation. As seen, the metasurface comprises a repeating pattern of elements 100, which in this example are arranged in a periodic array. Each of the elements is conformed as a cylindrical resonator body capped with a low-refractive-index layer 105 of silicon dioxide. A pump beam 110 is normally incident from the top, and generated second-harmonic light 115 is seen radiating in both reflection and transmission.

(28) FIG. 2 is a perspective view of a single resonator element 200 on a sapphire substrate 205, similar to an example embodiment that we fabricated and tested. As indicated in the figure, the element includes a resonator body comprising an InGaAs/AlInAs quantum coupled-well heterostructure 210, with a low-refractive-index cap 215 of silicon dioxide. As also indicated in the figure, the resonator body has a height h and a radius R. The aspect ratio of the resonator body is defined as R/h.

(29) Quantum coupled-well heterostructures are described, for example, in F. Capasso et al., “Coupled Quantum Well Semiconductors with Giant Electric Field Tunable Nonlinear Optical Properties in the Infrared,” IEEE J. Quantum Electr. 30 (1994) 1313-1325, cited below as “Capasso 1994” and hereby incorporated herein by reference in its entirety.

(30) FIG. 3 is an energy-band diagram for a quantum heterostructure used in our experiments. The respective layers of this heterostructure had the composition In.sub.0.53Ga.sub.0.47As/Al.sub.0.52In.sub.0.48As. The thicknesses of the layers in the sequence: (Al.sub.0.52In.sub.0.48As/In.sub.0.53Ga.sub.0.47As/Al.sub.0.52In.sub.0.48As/In.sub.0.53Ga.sub.0.47As/Al.sub.0.52In.sub.0.48As)
are 10/5/1.2/2.2/10 (all in nm). The height of the resonator element is 1.5 μm. An example aspect ratio R/h is 0.93, where R is cylinder radius and h is height.

(31) The heterostructure was designed to have electronic intersubband transitions (ISTs) at a fundamental Mie resonance of the resonator body with a wavelength of about 7.8 μm, and also at a second harmonic of the resonator body with a wavelength of about 3.9 μm.

(32) FIG. 3 also shows three wavefunctions 301, 302, 303 superimposed on the energy-band diagram. These are the wavefunctions for the three subbands that give rise to the intersubband transitions. The respective wavefunctions are drawn with vertical separations between them that correspond to the differences in their energy levels.

(33) For purposes of discussion, we will now refer to the three electronic states, in order of increasing energy, as <1>, <2>, and <3>. Accordingly, the three transitions that are involved in the process of SHG are 1.fwdarw.2, 2.fwdarw.3, and 1.fwdarw.3, with respective dipole matrix elements z.sub.12, z.sub.23, and z.sub.13. The energy difference corresponding to the 1.fwdarw.3 transition is hω.sub.3, where h is the reduced Planck constant.

(34) As explained in Capasso 1994, cited above, the second-order nonlinear susceptibility is theoretically predicted by the equation

(35) ${\chi }^{\left(2\right)}\left(2\omega \right)=\frac{e^3}{{\mathrm{.Math.}}_0}N.Math.\frac{.Math.z_{12}.Math..Math.z_{23}.Math..Math.z_{31}.Math.}{\left(2\mathrm{\hslash \omega }-{\mathrm{\hslash \omega }}_{13}-i{\Gamma }_{13}\right)},$
where e is the electronic charge, ε.sub.0 is the permittivity of free space, Nis the electron density in the quantum wells, and Γ.sub.13 is the half width at half-maximum of the 1.fwdarw.3 transition.

(36) One feature of Mie resonances that is favorable to SHG and other nonlinear phenomena is that at resonance, certain components of the electromagnetic field internal to the resonator are enhanced in value, relative to the incident field. Accordingly, afield enhancement f.sub.p is defined by the expression f.sub.p=|E.sub.z/E.sub.incident|, where E.sub.inc is the amplitude of the electric field at the pump wavelength in the incident wave, and E.sub.z is the electric field component on the axis normal to the metasurface (i.e., normal to the growth direction), which by convention is taken to be the z-axis. The z-directed field enhancement is important because quantum selection rules dictate that only photonic modes with E.sub.z components can couple to the ISTs.

(37) The SHG efficiency is proportional to the quantity |χ.sub.zzz.sup.(2)f.sub.p.sup.2|.sup.2, where f.sub.p is defined above, and χ.sub.zzz.sup.(2) is the (z,z,z) tensor element of the susceptibility χ.sup.(2)(2ω). For that reason, it is advantageous to engineer the Mie modes of the resonators to satisfy two conditions: (1) an efficient coupling of the incident pump light to the ISTs; and (2) a maximized volume-integrated product of xml and the field enhancement.

(38) To simultaneously satisfy both conditions, the resonators in the present example were designed to support a particular mode, namely, the lowest-order magnetic dipole (MD) mode (or, alternatively, the second-order magnetic quadrupole mode) at a wavelength slightly detuned from the fundamental IST wavelength of 7.8 μm, but still within the linewidth of the IST resonance. The IST resonance linewidth can be inferred from the experimentally measured normalized transmission spectrum of the multi-QW stack.

(39) For example, FIG. 4 presents a graph obtained from numerical modeling of our example MQW structure. Plotted on the left-hand scale is the intersubband nonlinear susceptibility 401 of the MQW structure as a function of the pump wavelength. In addition, we plotted on the right-hand scale the intersubband absorption spectrum 402 that we experimentally measured on the same structure. An inset 405 in the figure shows the waveguide configuration used for the absorption measurements. The input direction k is indicated. The absorption spectrum was produced by taking the ratio between two intensity values: The intensity of the transmitted transverse magnetic (TM) polarized input light, and the intensity of the transmitted transverse electric (TE) polarized input light.

(40) Because the MD mode has strong z-directed electric field components (due to its circulating out-of-plane electric fields), choosing the MD mode allowed us to strongly couple normally incident pump light to the ISTs. In addition, slight detuning of the MD mode with respect to the IST resonance helped to maximize |χ.sub.zzz.sup.(2f.sub.p.sup.2|.sup.2 by reducing absorption losses. In fact, as shown in FIG. 4, the maximum of |χ.sup.(2)| spectrally overlaps the maximum of the IST-induced absorption in the present example.

(41) FIG. 5 provides a typical example of field enhancement. It is an intensity map, in the xz-plane, of numerically calculated values of |E.sub.z/E.sub.inc| of the MD resonance of an illustrative resonator body similar to the resonator body of FIG. 2.

(42) Thus, in example embodiments, we generated a nonlinear response by introducing polaritonic coupling at the fundamental wavelength between the MD mode and the IST. It should be noted in this regard that the MD mode has the additional advantage that it is the mode for which the resonator body can be made as small as possible. We can also achieve similar coupling for the magnetic quadrupole mode and IST, which we have also experimentally demonstrated.

(43) It is also noteworthy that in embodiments, the height of the resonator body can be less than the vacuum wavelength of the pump light at resonance, and it can also be less than the in-material wavelength. Consequently, effective metamaterials can be made that are very thin, for example with total thicknesses that are less than twice the vacuum wavelength of the pump light.

(44) Turning now to FIG. 6, three plots are shown of experimentally measured linear reflectance spectra. From left to right, these spectra 601, 602, 603 were measured on resonators with different aspect ratios of 0.93, 1, and 1.06 respectively. The reflectance peak at the MD resonance is labeled on each of the three spectral curves.

(45) It can be seen in the figure that at an aspect ratio of 1, the MD peak is split. This is an indication that the MD resonance is aligned with the IST. This is an example of the vacuum-Rabi splitting that is the characteristic signature of strong coupling and formation of an intersubband polariton.

(46) It will be evident from this example that using InGaAs/AlInAs heterostructures gives us access to a spectral range at least from 7 μm to 8 μm, and beyond, for potential second harmonic conversion. This is significant, not least because it provides a path to generating light in the mid-wave infrared region, i.e., the portion of the electromagnetic spectrum from 3 μm to 8 μm.

Example 1

(47) Fabrication. We fabricated metasurfaces substantially as described above. FIG. 7 is a flowchart in which cross-sectional cartoon views illustrate several sequential stages in our fabrication process.

(48) As shown in the figure, a quantum heterostructure is grown epitaxially on an indium phosphide (InP) substrate (701). The substrate is flipped, and the heterostructure is adhered, topside-down, to a sapphire substrate using an epoxy adhesive (702). The InP substrate is removed by optional lapping, followed by wet etching. The cylindrical Mie resonators are then patterned using electron beam lithography (703) and defined by a combination of inductively coupled plasma etching (ICP) and reactive ion etching (RIE) (704).

(49) More specifically, the heterostructure used in this study was grown on a semi-insulating InP substrate using molecular beam epitaxy. It comprises multiple repetitions of the asymmetric coupled In.sub.0.53Ga.sub.0.47As QWs sandwiched between Al.sub.0.52In.sub.0.48As barrier layers as shown in FIG. 3. The total thickness of the multi-QW stack was 852 nm. Since the minimum height of the resonator needs to be about 1.5 μm to support the lowest order MD mode, additional Al.sub.0.52In.sub.0.48As layers were grown symmetrically on top and bottom of the multi-QW stack to make the total epilayer thickness equal to 1.5 μm. The epilayer thickness determined the height (h) of the resonators used in this study to be 1.5 μm.

(50) After growth, the wafer was flipped and adhesively bonded to a sapphire substrate using benzocyclobutene. After bonding, the InP substrate was selectively removed to leave the epilayer on top of the sapphire substrate.

(51) The resonators were then defined via e-beam lithography. Hydrogen silsesquioxane was used as the e-beam resist which was converted to SiO.sub.2 on e-beam exposure. After development, inductively-coupled plasma reactive ion etching was used for defining and fabricating the cylinders using the SiO.sub.2 layer on top of the cylinders as a mask. The SiO.sub.2 was not removed after the etching and together with the sapphire substrate on the bottom, it provided the refractive index contrast needed for the resonators to support the Mie modes.

(52) The Mie modes were spectrally scaled by changing the radius (R) of the resonators. In all the metasurfaces, the periodicity (p) of the array of the resonators along x and y was fixed to be 3R.

(53) Experimental Results. FIG. 8 shows the experimentally measured optical power of the SH signal reflected from the Mie metasurfaces as a function of the pump wavelength. Three metasurfaces are represented in the figure, with respective cylinder radii of 1.2 μm (curve 801), 1.3 μm (curve 802), and 1.4 μm (curve 803). The cylinder height in each case was 1.5 μm.

(54) The designed resonant wavelength of the IST in these metasurfaces was 7.8 μm. However, it is evident from the figure that the maximum SH generation efficiency occurred at an aspect ratio R/h of 0.93 and a pump wavelength of only about 7.62 μm. That is, the maximum SH generation efficiency is observed at a pump wavelength slightly shorter than the designed resonant wavelength of the IST. Without wishing to be limited by mechanism, we propose that this is because the field enhancement inside the resonators is larger at the shorter wavelengths.

(55) FIG. 9 is a graph of the experimentally measured linear reflectance spectrum of a metasurface with cylinders of R=1.4 μm. Rabi splitting of the MD resonance can be seen, indicating strong light-matter coupling between the MD resonance and the IST, leading to formation of intersubband polaritons. MD, ED, MQ, and EQ correspond to magnetic dipole, electric dipole, magnetic quadrupole, and electric quadrupole modes, respectively.

(56) FIGS. 10, 11, and 12 illustrate the interplay between the nonlinearity and the field enhancement.

(57) The SHG efficiency is determined by the product of the second-order nonlinear susceptibility times the field enhancement. The numerically calculated susceptibility is plotted versus wavelength in FIG. 10. The numerically calculated field enhancement is plotted versus wavelength in FIG. 11. In FIGS. 10 and 11, a vertical black broken line indicates the designed resonant wavelength of the IST.

(58) Although the susceptibility is maximum at 7.8 μm, the field enhancement is larger at shorter wavelengths. As a result, the maximum SHG efficiency will be shifted to a wavelength shorter than the wavelength of maximum susceptibility.

(59) FIG. 12 is a graph of the experimentally measured reflected SH power conversion efficiency as a function of pump wavelength. The measurements were for normal pump incidence a metasurface with cylinders of R=1.4 μm. The efficiency is plotted in normalized units obtained by dividing the reflected SH power by the square of the incident pump power. Two peaks at about 7.6 μm and about 8 μm are observed, corresponding to the two polariton branches. The maximum SHG efficiency is observed at a pump wavelength of about 7.65 μm.

(60) In FIG. 13, we have plotted experimental measurements of the peak reflected SH power optical power versus the square of the pump optical power for h=1.5 μm, R=1.4 μm, and a wavelength of 7.65 μm.

(61) The second-harmonic nonlinear conversion factor can be determined from the slope of a linear fit 1300 to the curve of FIG. 13. As seen in the figure, this factor reached values as high as 0.5 mW/W.sup.2.

(62) Another figure of merit is the second-harmonic power conversion efficiency, which is defined as the ratio of second-harmonic power to pump power.

(63) FIG. 14 is a graph of the experimentally measured conversion efficiency versus incident peak pump power for the metasurface of FIG. 13. The inset shows the corresponding experimentally-measured spectrum of SH radiation. The single spectral peak at one-half the pump wavelength, confirms the SHG process.

(64) The measurements plotted in FIGS. 13 and 14 were corrected by a factor of 4 to compensate for limited collection efficiency of our optics.

(65) It can be seen from FIG. 14 that the maximum conversion efficiency is 0.013% at a peak pump intensity of 11 kW/cm.sup.2.

(66) As FIG. 14 shows, the effectiveness of ISTs can be limited by the saturation of the IST transition at high pump powers. However, the pump intensity at saturation can be significantly increased by detuning one of the resonances by a few percent, rather than using equally spaced electronic energy levels as in the present example.

Example 2

(67) We performed numerical calculations of the SH generation efficiency for an array of resonators of several finite sizes. Our results are shown in FIG. 15.

(68) The figures shows the conversion efficiency, in normalized units, for normal incidence on the metasurface as a function of pump wavelength. The radius R of the cylindrical resonator elements is fixed at R=1.4 μm. Each array of size Nis a square N×N array with a period of 3R.

(69) It can be seen from the figure that the spectral dependence of normalized SHG efficiencies is similar for different array sizes. We confirmed that there was a good correspondence to the experimental data.

(70) It can also be seen from the figure that the absolute magnitude of the SHG efficiency initially increases rapidly as the number N of resonator elements in the array (i.e., the array size) increases from 1 to about 25. As N continues to increase above 25, the efficiency converges toward a value of about 4 mW/W.sup.2. The field enhancements that we have computed for finite arrays are always smaller than those for infinite arrays. Consequently, the finite-sized arrays always have lower SHG efficiencies than a comparable infinite array.

(71) Normal emission of second-harmonic light. As noted above, normal emission of second-harmonic light is suppressed when the pump beam is at normal incidence. This is a consequence of the in-plane inversion symmetry of the cylindrical resonator structures, together with the presence of only a single element χ.sub.zzz.sup.(2) in the second-order nonlinear susceptibility tensor. Because of these factors, the nonlinear polarization at the SH frequency is induced only in the z direction. This, in turn, results in the SH radiation being emitted only at off-normal angles, with zero SH power radiated along the normal direction.

(72) Measures can be taken, however, to break the in-plane inversion symmetry and thereby permit normal emission of SH radiation. One possible approach, for example, is to reduce the symmetry of the individual resonator elements. Another possible approach is to modify the array so that it has a unit cell that lacks inversion symmetry. That could be achieved, for example, with a unit cell that is a cluster of two or more resonator elements in an arrangement without inversion symmetry.

(73) The simplest approach, however, would be to tilt the direction of the pump beam so that the irradiation is off-normal. A suitable tilt angle would suffice to break the inversion symmetry and permit SH power to be radiated in the normal direction.

(74) We performed numerical simulations to test this approach. Our results are shown in FIG. 16, which is an intensity map in which the quantities plotted on the respective axes are pump angle of incidence and pump wavelength, and the intensity represents the SH power conversion efficiency. It will be evident from the figure that near the design pump wavelength, the SH power conversion efficiency is greatest at tilt angles in the range 10°-20°.

(75) Controlling the sign of the susceptibility. As mentioned above, it is possible to control both the magnitude and sign of χ.sup.(2) along the growth direction. The way this would be done is analogous to periodic poling of nonlinear crystals. Such control is possible because χ.sup.(2) of the IST systems is proportional to the product of the three transition dipole moments z.sub.12, z.sub.23, and z.sub.13. The respective signs of these moments can be flipped by reversing the order of the growth of the quantum wells in the asymmetric quantum well pairs.

(76) The SH generation efficiency in each resonator is proportional to the product of χ.sup.(2) and the second-harmonic electric field profile. Thus, it would be possible, in principle, to design Mie resonators in which the sign of χ.sup.(2) is engineered to locally match the polarity of the SH field and thereby enhance the SH generation efficiency.

(77) Other nonlinear optical phenomena. IST-based all-dielectric Mie metasurfaces can also be designed for other nonlinear optical phenomena such as sum frequency generation, third-harmonic generation, etc. They therefore offer an exciting and fundamentally new platform for studying nonlinear optics at the nanoscale. Such polaritonic all-dielectric metasurfaces thus offer a paradigm shift in the field of nonlinear metasurfaces that can offer new opportunities for developing new classical or quantum light sources for various applications ranging from spectroscopy, imaging, to sensing.

(78) By way of example, FIG. 17 is a diagram showing the growth sequence for a four-level IST system capable of generating both second-harmonic and third-harmonic radiation. As will be evident in the figure, there are now three asymmetric potential wells 1701, 1702, 1703. The III-V material system is the same as in the preceding examples. The widths of the respective potential barriers between wells are 1.6 nm The widths of the wells, in sequence from left to right in the view of the figure, are 4.3 nm, 2 nm, and 1.8 nm

(79) Alternative material systems. The example embodiments described above relate to quantum-well heterostructures fabricated in the InGaAs/AlInAs material system. Those examples should be understood as illustrative and non-limiting. Other material systems are also under contemplation as possible host systems for devices according to the principles described here and therefore falling within the scope and spirit of the present invention.

(80) For example, alternate material systems may be based on Group III antimonide compounds. In one such system, the quantum well heterostructures may be constituted by aluminum antimonide (AlSb) and indium arsenide (InAs).

(81) By way of illustration, FIG. 18 is an energy-level diagram of an AlSb/InAs coupled well heterostructure. Wavefunctions 1801, 1802, 1803 for the lowest three electronic states are notionally represented in the figure, with vertical separations corresponding to their separations in energy. The energy separations that are shown could potentially support the generation of second-harmonic light at a wavelength of about 1.5 μm from pump light at a wavelength of about 3 μm.