Gain-Assisted Giant Third-Order Nonlinearity of Epsilon-Near-Zero Multilayered Metamaterials

We investigate the third-order nonlinear optical properties of epsilon-near-zero (ENZ) Au/dye-doped fused silica multilayered metamaterials in the visible spectral range for TM incident by using nonlocal effective medium theory at different incidence angles. The nonlocal response affects the permittivity of anisotropic metamaterials when the thickness of the layer cannot be much smaller than the incident wavelength. By doping pump dye gain material within the dielectric layer to compensate for the metal loss, the imaginary part of the effective permittivity is reduced to 10−4, and the optical nonlinear refractive index and nonlinear absorption coefficient are enhanced. The real and imaginary parts of the permittivity are simultaneously minimized when the central emission wavelength of the gain material is close to the ENZ wavelength, and the nonlinear refraction coefficient reaches the order of 10−5 cm2/W, which is five orders of magnitude larger than that of the nonlinear response of the metamaterial without the gain medium. Our results demonstrate that a smaller imaginary part of the permittivity can be obtained by doping gain materials within the dielectric layer; it offers the promise of designing metamaterials with large nonlinearity at arbitrary wavelengths.


Introduction
Recently, epsilon-near-zero (ENZ) materials are attracting great interest due to their optical tunneling effect [1], tailoring of radiation patterns [2], enhanced nonlinearity [3][4][5][6][7][8][9][10][11][12][13][14][15], broadband [16] and high-speed tunable devices [17][18][19][20]. Currently, there are two main mechanisms to realize ENZ with different mechanisms and materials: (1) natural materials near the bulk plasma and phonon resonances; (2) anisotropic metamaterials with alternating stacks of negative permittivity metals, positive permittivity dielectrics [21][22][23][24] and metal pillars embedded in dielectrics [25,26]. The ENZ condition can be achieved at any wavelength of interest by reasonably adjusting the thickness of the metal and dielectric layers and the diameter of the metal pillars in the anisotropic metamaterial. The effective medium theory is valid only when the layer thickness of the anisotropic metamaterial is much smaller than that of the wavelength. However, the layer thickness is too thin; in a few nanometers, it leads to the appearance of unwanted effects in materials, such as: dislocations, island formation and size effects [27]. The effect of spatial dispersion, i.e., nonlocal, becomes increasingly important under the ENZ condition of the anisotropic metamaterial; the appearance of strong nonlocal properties significantly alters the optical properties of the material and can lead to the appearance of additional transverse or longitudinal waves [28,29]. Therefore, it is necessary to design the structure with a more suitable nonlocal effective medium theory [30].
ENZ metamaterials vanish a real part of the permittivity, but metals with a large imaginary part of the permittivity have large optical losses, which seriously affects practical applications and is not conducive to improving the nonlinear response. To circumvent these constraints, optical losses can be compensated for by a gain-doped medium with optically pumped dye molecules [31][32][33][34][35][36] or quantum dot [37][38][39][40][41]. Gain materials have little effect on the real part of the permittivity and greatly reduce the imaginary part of the permittivity at the ENZ wavelength. Due to its negative imaginary part and the essentially invariant real part of the effective permittivity, it improves the Kerr-like optical nonlinear response.
In this paper, the linear and nonlinear optical properties of the ENZ anisotropic multilayered metamaterial with alternating Au/dye-doped fused silica at different incidence angles in the visible range for TM light are investigated by considering spatial dispersion using a more comprehensive theoretical approach, i.e., nonlocal effective medium theory. The gain-doped medium reduces the imaginary part of the transverse permittivity to 10 −4 when the emission wavelength of the gain material is close to the ENZ wavelength. Compared with the metamaterial without gain medium, it has a narrower resonance bandwidth and a 4.92-fold increase in the electric field squared at the ENZ wavelength; the nonlinear coefficient is improved by 5 orders of magnitude.

Theory and Method
The multilayered metamaterial is composed of alternating Au/dye-doped fused silica with six cycles arranged on the SiO 2 substrate, as shown in Figure 1. The thickness of gold is 14.9 nm and dye-doped fused silica is 116.2 nm. The optical loss of the metamaterial is compensated for by gain-doped medium with pump dye molecule rhodamine R800. The relative permittivity of this gain-doped medium is a simple Lorentz model [42].
where Γ pump = δ abs I pump h f 30 = 6.5 × 10 9 s −1 is the pumping rate, δ abs = 3.14 × 10 −16 cm 2 is the absorption cross section, f 30 = 441 THz is the absorption frequency. ε 0 is the vacuum permittivity, ε 1 = 2.25 is the dielectric permittivity, δ a is the coupling strength of the polarization density in the emission frequency, λ a = 711 nm , ω a = 422 THz is the center emission frequency, ∆ω a is the frequency line-width, τ 21 = 500 ps, τ 32 = τ 10 =100 fs, N 0 is concentration of the R800 and N 0 = 1.98 × 10 18 cm −3 corresponds to about 3.3 mM. Figure 2 shows the variation in the relative permittivity for the gain-doped medium with respect to the wavelength, the vertical dashed line is the central emission wavelength and the negative imaginary part of the permittivity clearly shows the loss compensated in the gain-doped medium.    The loss is compensated for when the central emission wavelength of the gain medium is close to the ENZ wavelength of the ENZ multilayered metamaterial. The effective medium theory [43,44] is valid only when the layer thickness is small, and the nonlocal effect is not negligible at larger layer thicknesses and under the ENZ condition, which qualitatively changes the optical properties of the structure. The complex effective permittivity in-plane ∥ = and out-of-plane = is calculated by using the following nonlocal effective medium theory [30].
For the metamaterial without gain medium, where α, β, γ, , are: The loss is compensated for when the central emission wavelength of the gain medium is close to the ENZ wavelength of the ENZ multilayered metamaterial. The effective medium theory [43,44] is valid only when the layer thickness is small, and the nonlocal effect is not negligible at larger layer thicknesses and under the ENZ condition, which qualitatively changes the optical properties of the structure. The complex effective permittivity in-plane ε = ε nloc xx or ε nloc yy and out-of-plane ε ⊥ = ε nloc zz is calculated by using the following nonlocal effective medium theory [30].
ε nloc yy = ε loc yy 1 + For the metamaterial without gain medium, where α, β, γ, ε loc xx , ε loc zz are: For the metamaterial with gain medium, where α, β, γ, ε loc xx , ε loc zz are: The permittivity of effective medium theory in-plane and out-of-plane is ε loc = ε loc xx = ε loc yy = ε loc xy and ε loc ⊥ = ε loc zz , respectively. k 0 = 2π λ is the incident wave vector, k x = k 0 sin θ i , k z = k 0 cos θ i , θ i is the incident angle, ε d = 2.25 is the medium relative permittivity, the metal permittivity ε m can be calculated by the Drude-Lorentz model [45], f m = a 1 a is the metal filling ratio, f d = a 2 a is the medium filling ratio, a 1 is the metal thickness, a 2 is the medium thickness and a = a 1 + a 2 .
The permittivity of TE and TM waves propagating at the angle of θ in the metamaterial is as follows [25]: It is worth noting that nonlocal parallel permittivities vary from each other in terms of magnitude, as shown in Figure 3a. This means that nonlocality induces an effective biaxiality of the ENZ multilayered metamaterial, even though the structure is spatially uniform in the XY-plane directions. Compared with the local effective medium theory, the ENZ wavelength occurring within the visible spectral range is red-shifted, and the imaginary part is reduced. The ENZ wavelength in the x-direction is 711 nm, and the imaginary part is 0.1207. The nonlocality has a great influence on the out-of-plane permittivity, as shown in Figure 3b. Since we are particularly interested in how the gain medium affects the nonlinearity, only ε = ε nloc xx is considered later. Figure 4 shows the relative permittivity versus wavelength after adding the gain medium; the gain medium has little effect on the real part of the permittivity and a great effect on the imaginary part, and the imaginary part of the permittivity is reduced to 5.36 × 10 −4 at the ENZ wavelength 711 nm.
the nonlinearity, only ∥ = is considered later. Figure 4 shows the relative permittivity versus wavelength after adding the gain medium; the gain medium has little effect on the real part of the permittivity and a great effect on the imaginary part, and the imaginary part of the permittivity is reduced to 5.36 × 10 −4 at the ENZ wavelength 711 nm.

Linear Response
The optical properties are affected not only by changing the layer thickness and permittivity, but also by different incidence angles. The characteristics of reflectance R, transmittance T and absorption A of the metamaterial without gain medium for TM wave at different incidence angles are shown in Figure 5, absorption A = 1 -T − R, and the dashed line is the ENZ wavelength 711 nm. The plasmon-polaritonic coupling of the metal-dielectric interface greatly confines the electric field intensity in the dielectric layer [46][47][48][49], leading to multiple distinct resonances of reflectance, transmittance and absorption between wavelengths of 450-800 nm. The metal-dielectric-metal sandwich structure in the the nonlinearity, only ∥ = is considered later. Figure 4 shows the relative permittivity versus wavelength after adding the gain medium; the gain medium has little effect on the real part of the permittivity and a great effect on the imaginary part, and the imaginary part of the permittivity is reduced to 5.36 × 10 −4 at the ENZ wavelength 711 nm.

Linear Response
The optical properties are affected not only by changing the layer thickness and permittivity, but also by different incidence angles. The characteristics of reflectance R, transmittance T and absorption A of the metamaterial without gain medium for TM wave at different incidence angles are shown in Figure 5, absorption A = 1 -T − R, and the dashed line is the ENZ wavelength 711 nm. The plasmon-polaritonic coupling of the metal-dielectric interface greatly confines the electric field intensity in the dielectric layer [46][47][48][49], leading to multiple distinct resonances of reflectance, transmittance and absorption between wavelengths of 450-800 nm. The metal-dielectric-metal sandwich structure in the

Linear Response
The optical properties are affected not only by changing the layer thickness and permittivity, but also by different incidence angles. The characteristics of reflectance R, transmittance T and absorption A of the metamaterial without gain medium for TM wave at different incidence angles are shown in Figure 5, absorption A = 1 -T − R, and the dashed line is the ENZ wavelength 711 nm. The plasmon-polaritonic coupling of the metaldielectric interface greatly confines the electric field intensity in the dielectric layer [46][47][48][49], leading to multiple distinct resonances of reflectance, transmittance and absorption between wavelengths of 450-800 nm. The metal-dielectric-metal sandwich structure in the multilayered metamaterial is similar to an F-P resonant cavity, which enhances the electric field intensity and, therefore, produces transmittance windows, absorption peaks and reflectance dips near the ENZ wavelength. The resonances of transmittance, reflectance and absorption exhibit a similar blueshift phenomenon with increasing incident angle, and the transmittance of the multilayered metamaterial increases with increasing incident angle, but the reflectance and absorption conversely decrease. multilayered metamaterial is similar to an F-P resonant cavity, which enhances the electric field intensity and, therefore, produces transmittance windows, absorption peaks and reflectance dips near the ENZ wavelength. The resonances of transmittance, reflectance and absorption exhibit a similar blueshift phenomenon with increasing incident angle, and the transmittance of the multilayered metamaterial increases with increasing incident angle, but the reflectance and absorption conversely decrease. In Figure 6, we reveal transmittance characteristics of the structure with gain medium at different incident angles. The transmittance resonance peak at 711 nm in Figure 6 is narrower than that in Figure 5a because of the reduced permittivity after adding the gain medium. The structure also has absorption at 711 nm, but the absorption is smaller than the gain, so there is a narrow-band amplification of the transmittance peak, and transmittance resonance peaks at the remaining wavelengths are similar to Figure 5a, because the permittivity is the same as that of the structure without gain medium. According to the continuity of the electric displacement field for the boundary conditions at the interface In Figure 6, we reveal transmittance characteristics of the structure with gain medium at different incident angles. The transmittance resonance peak at 711 nm in Figure 6 is narrower than that in Figure 5a because of the reduced permittivity after adding the gain medium. The structure also has absorption at 711 nm, but the absorption is smaller than the gain, so there is a narrow-band amplification of the transmittance peak, and transmittance resonance peaks at the remaining wavelengths are similar to Figure 5a, because the permittivity is the same as that of the structure without gain medium. According to the continuity of the electric displacement field for the boundary conditions at the interface for ENZ metamaterials, that is, |E| ∝ ε −1 |E 0 |, ENZ metamaterials with the gain medium have a strong squared of electric field (|E| 2 ) enhancement due to the small magnitude of the permittivity (ε), which is 4.92 times larger than in the case without the gain medium at normal incidence. The electric field intensity at the ENZ wavelength also decreases with increasing angle due to the blueshift of the transmittance resonance wavelength and the decrease in the transmittance peak, as shown in Figure 7d-f. for ENZ metamaterials, that is, | | ∝ | |, ENZ metamaterials with the gain medium have a strong squared of electric field (| | ) enhancement due to the small magnitude of the permittivity ( ), which is 4.92 times larger than in the case without the gain medium at normal incidence. The electric field intensity at the ENZ wavelength also decreases with increasing angle due to the blueshift of the transmittance resonance wavelength and the decrease in the transmittance peak, as shown in Figure 7d-f.   for ENZ metamaterials, that is, | | ∝ | |, ENZ metamaterials with the gain medium have a strong squared of electric field (| | ) enhancement due to the small magnitude of the permittivity ( ), which is 4.92 times larger than in the case without the gain medium at normal incidence. The electric field intensity at the ENZ wavelength also decreases with increasing angle due to the blueshift of the transmittance resonance wavelength and the decrease in the transmittance peak, as shown in Figure 7d-f.

Nonlinear Response
The nonlinear refractive index and nonlinear absorption coefficient of the ENZ multilayered metamaterial were investigated using nonlinear effective medium theory [46]. Here, the effective nonlinear third-order susceptibility χ 3 e f f of the structure is a weighted average of the constituent materials; since the nonlinear susceptibility of dielectric χ (3) d is much smaller than that of the Au χ  [50][51][52][53]. However, we are interested in the effect of the gain medium on the third-order nonlinearity rather than optimizing the optimal third-order nonlinearity. In what follows, we assume that Au is dispersion-less [23] and that the third-order susceptibility of Au is χ [54]. The third-order nonlinear susceptibility of the ENZ multilayered metamaterial parallel to the interface and perpendicular to the interface is as follows [49].
where ε nloc ⊥ ε m displays the local field enhancement factor in the metallic component. For the TM wave at all incidence angles, the third-order susceptibility is derived as follows: where ε NL and ε NL ⊥ are the in-plane and out-of-plane nonlinear permittivity, respectively [23]: By replacing ε NL and ε NL ⊥ in Equation (19) with Equations (20-21), the nonlinearity in the denominator is too small compared with linearity, so the effects of the nonlinearity in the denominator and the electric field |E| 4 can be neglected.
To make a further approximation, assuming that the nonlinear component is much smaller than the linear component, Snell's law can be written as ε NL TM (θ) sin 2 (θ)=ε TM (θ) sin 2 (θ)=ε i sin 2 (θ i ). θ i is the angle of incidence and ε i is the permittivity of the air medium, therefore, Equation (22) can be rewritten as: The complex nonlinear response is: where ε 0 is the permittivity in a vacuum, n 0 is the complex effective linear refractive index, Re(n 0 ) is the real part of the complex effective linear refractive index and c is the speed of light. The nonlinear refraction coefficient and nonlinear absorption coefficient are as follows: We investigated the intensity-dependent refraction index n of the multilayered metamaterial. The nonlinear refractive index change is greater than the linear refractive index change under ENZ conditions, which violates the nonlinear perturbation theory [3]. Therefore, the conventional intensity-dependent refractive index, n(I) = n 0 + n 2 I, is not applicable; we directly use the susceptibility to calculate the intensity-dependent refractive index, i.e., n(I) = ε (1) + 3χ (3) |E| 2 . The incident light intensity I directly affects the permittivity of the gain medium as in Equation (1). The permittivity imaginary Im(ε ) is closest to 0 as the incident light intensity increases to 6.15 MW/cm 2 , and the refractive index Re(n) = 5.57 is at its maximum at I = 6.15 MW/cm 2 , as shown in Figure 8b, but the nonlinear bandwidth becomes narrower than I = 6.05 MW/cm 2 . When I > 6.15 MW/cm 2 , Re(n) decreases due to the increase in Im(ε ). We use I = 6.05 MW/cm 2 in the following after considering the nonlinear coefficient and the bandwidth size. For metamaterials without the gain medium, Re(n) in Figure 8a is much smaller than that in Figure 8b with the gain medium. The above results show that multilayered metamaterials with gain media exhibit an order of magnitude of nonlinear refractive index change under ENZ conditions. Nanomaterials 2022, 12, x FOR PEER REVIEW 9 of 16 are as follows: We investigated the intensity-dependent refraction index n of the multilayered metamaterial. The nonlinear refractive index change is greater than the linear refractive index change under ENZ conditions, which violates the nonlinear perturbation theory [3]. Therefore, the conventional intensity-dependent refractive index, ( ) = + , is not applicable; we directly use the susceptibility to calculate the intensity-dependent refractive index, i.e., ( ) = ( ) + 3 ( ) | | . The incident light intensity I directly affects the permittivity of the gain medium as in Equation (1). The permittivity imaginary ( ∥ ) is closest to 0 as the incident light intensity increases to 6.15 MW/cm 2 , and the refractive index Re(n) = 5.57 is at its maximum at I = 6.15 MW/cm 2 , as shown in Figure 8b, but the nonlinear bandwidth becomes narrower than I = 6.05 MW/cm 2 . When I > 6.15 MW/cm 2 , Re(n) decreases due to the increase in ( ∥ ). We use I = 6.05 MW/cm 2 in the following after considering the nonlinear coefficient and the bandwidth size. For metamaterials without the gain medium, Re(n) in Figure 8a is much smaller than that in Figure 8b with the gain medium. The above results show that multilayered metamaterials with gain media exhibit an order of magnitude of nonlinear refractive index change under ENZ conditions. The metamaterial displays a strong wavelength-dependent nonlinear response; the conversion between self-focusing to self-defocusing and saturated absorption and antisaturated absorption can be achieved by controlling the wavelength. For the metamaterial without gain medium, the | | is the maximum and reaches the order of 10 −10 cm 2 /W at the wavelength of 717.8 nm, as shown in Figure 9a. For the metamaterial with gain medium, as shown in Figure 10a, the complex linear refractive index is well reduced and the third-order nonlinear susceptibility is increased due to the reduction in the imaginary part of the dielectric constant at 710.8 nm, thus, the | | is the maximum and reaches the order of 10 −5 cm 2 /W, which is an increase of 5 orders of magnitude compared with the nonlinear refractive index for normal incidence in Figure 9a. Similarly, for the metamaterial without gain medium, the |β| at 707.3 nm is the maximum and reaches the order of 10 −5 cm/W, as shown in Figure 9b. For the metamaterial with gain medium, as shown in Figure 10b, the The metamaterial displays a strong wavelength-dependent nonlinear response; the conversion between self-focusing to self-defocusing and saturated absorption and antisaturated absorption can be achieved by controlling the wavelength. For the metamaterial without gain medium, the |n 2 | is the maximum and reaches the order of 10 −10 cm 2 /W at the wavelength of 717.8 nm, as shown in Figure 9a. For the metamaterial with gain medium, as shown in Figure 10a, the complex linear refractive index is well reduced and the third-order nonlinear susceptibility is increased due to the reduction in the imaginary part of the dielectric constant at 710.8 nm, thus, the |n 2 | is the maximum and reaches the order of 10 −5 cm 2 /W, which is an increase of 5 orders of magnitude compared with the nonlinear refractive index for normal incidence in Figure 9a. Similarly, for the metamaterial without gain medium, the |β| at 707.3 nm is the maximum and reaches the order of 10 −5 cm/W, as shown in Figure 9b. For the metamaterial with gain medium, as shown in Figure 10b, the |β| is the maximum at the wavelength of 710.8 nm and reaches the order of 1 cm/W, which is an increase of 5 orders of magnitude in the nonlinear absorption coefficient compared with Figure 9b. Table 1 shows the comparisons of nonlinear coefficients (nonlinear refraction coefficient and nonlinear absorption coefficient) for similar systems, and the results of this article are shown in bold. The nonlinear refraction coefficient and nonlinear absorption coefficient reach the order of 10 −5 cm 2 /W and 1 cm/W, respectively, which are several orders of magnitude larger than the results from other literature. Therefore, the gain-doped medium improves the nonlinear response, and the dominant mechanism for improving the nonlinearity is the factor 1/Re(n 0 )n 0 . |β| is the maximum at the wavelength of 710.8 nm and reaches the order of 1 cm/W, which is an increase of 5 orders of magnitude in the nonlinear absorption coefficient compared with Figure 9b. Table 1 shows the comparisons of nonlinear coefficients (nonlinear refraction coefficient and nonlinear absorption coefficient) for similar systems, and the results of this article are shown in bold. The nonlinear refraction coefficient and nonlinear absorption coefficient reach the order of 10 −5 cm 2 /W and 1 cm/W, respectively, which are several orders of magnitude larger than the results from other literature. Therefore, the gain-doped medium improves the nonlinear response, and the dominant mechanism for improving the nonlinearity is the factor 1 Re( ) ⁄ .
(a) (b)    [4] AZO film 1310/10 −13 1340/10 −8 [5] Metasurface 1240/10 −9 1260/10 −5 [10] Metasurface 1250/10 −9 1300/10 −4 [23] Multilayered metamaterial 470/10 −8 470/10 −3 [24] Multilayered metamaterial 500/10 −9 500/10 −4 Finally, we compare nonlinear coefficients of gain media with different incidence angles and different concentrations: 0 mM in Figure 11, 2.5 mM in Figure 12, 3.3 mM in Figure 13. Since the imaginary part of the permittivity increases and the electric field intensity coupled in the medium decreases with the increasing angle, the absolute value of the nonlinear coefficient decreases over the entire wavelength range. We consider concentrations of R800 dye molecules in the dielectric layer from 0 mM to 3.3 mM; on the one hand, the reduction in the imaginary part of the effective permittivity in the plane, Im ε ≈ 0.1207 for 0 mM to Im ε ≈ 10 −4 for 3.3 mM, i.e., the loss is well compensated. On the other hand the narrower resonance bandwidth leads to an increase in the intensity of the coupled electric field in the dielectric layer with increasing concentration of the gain medium. The nonlinear coefficients, near the ENZ wavelength, are increased by 1 order and 5 orders of magnitude for concentrations of 2.5 mM and 3.3 mM, respectively. Further, a broadband nonlinear response is shown for 0 mM concentration. When the gain medium concentration increases and loss is compensated, the nonlinear displays an angle/spectrum narrowband characteristic. Therefore, the nonlinearity of the ENZ multilayered metamaterial must be improved with an appropriate concentration of the gain medium. Table 2 shows the comparisons of nonlinear coefficients (nonlinear refraction coefficient, nonlinear absorption coefficient) at different doping concentrations. and 5 orders of magnitude for concentrations of 2.5 mM and 3.3 mM, respectively. Further, a broadband nonlinear response is shown for 0 mM concentration. When the gain medium concentration increases and loss is compensated, the nonlinear displays an angle/spectrum narrowband characteristic. Therefore, the nonlinearity of the ENZ multilayered metamaterial must be improved with an appropriate concentration of the gain medium. Table  2 shows the comparisons of nonlinear coefficients (nonlinear refraction coefficient, nonlinear absorption coefficient) at different doping concentrations.

Conclusions
In conclusion, the linear and nonlinear properties of the ENZ multilayered amaterial are investigated based on the nonlocal effective medium theory, and the wavelength of interest can be adjusted arbitrarily by adjusting the thickness of the d tric and metal layers. The metal-dielectric-metal in the ENZ multilayered metamate similar to the F-P cavity; the electric field intensity can be coupled in the dielectric but metals have large optical losses, leading to a large imaginary part of the permit

Conclusions
In conclusion, the linear and nonlinear properties of the ENZ multilayered metamaterial are investigated based on the nonlocal effective medium theory, and the ENZ wavelength of interest can be adjusted arbitrarily by adjusting the thickness of the dielectric and metal layers. The metal-dielectric-metal in the ENZ multilayered metamaterial is similar to the F-P cavity; the electric field intensity can be coupled in the dielectric layer, but metals have large optical losses, leading to a large imaginary part of the permittivity, which seriously affects practical applications and is not conducive to improving the nonlinear response. Therefore, the doping of R800 dye molecules with a concentration of 3.3 mM in the dielectric layer can be used to reduce the imaginary part of the permittivity in the ENZ multilayered metamaterial to 5.36 × 10 −4 . Compared with the metamaterial without gain medium, the squared electric field is enhanced by 4.92 times, and the nonlinear coefficient is increased by 5 orders of magnitude. Moreover, the conversion of self-focusing and self-defocusing as well as saturable absorption and anti-saturable absorption can be achieved by controlling the wavelength, and the absolute value of the nonlinear coefficient decreases with increasing angle. We believe that our investigation of the gain-doped ENZ multilayered metamaterial can lead to important applications in nonlinear optics and quantum optics.