Multistability at arbitrary low optical intensities in a metallo-dielectric layered structure

We show that a nonlinear metallo-dielectric layered slab of subwavelength thickness and very small average dielectric permittivity displays optical multistable behavior at arbitrary low optical intensities. This is due to the fact that, in the presence of the small linear permittivity, one of the multiple electromagnetic slab states exists no matter how small is the transmitted optical intensity. We prove that multiple states at ultra-low optical intensities can be reached only by simultaneously operating on the incident optical intensity and incidence angle. By performing full wave simulations, we prove that the predicted phenomenology is feasible and very robust.

We show that a nonlinear metallo-dielectric layered slab of subwavelength thickness and very small average dielectric permittivity displays optical multistable behavior at arbitrary low optical intensities. This is due to the fact that, in the presence of the small linear permittivity, one of the multiple electromagnetic slab states exists no matter how small is the transmitted optical intensity. We prove that multiple states at ultra-low optical intensities can be reached only by simultaneously operating on the incident optical intensity and incidence angle. By performing full wave simulations, we prove that the predicted phenomenology is feasible and very robust. Hysteresis is probably one of the most intriguing feature a nonlinear system can exhibit both theoretically and for its main application, i.e. designing of memory devices. In optics, bistability and related hysteresis have attracted a large research interest in the last decades [1,2] since they are the basic ingredients for devising optical memories, logic gates and optical computing devices [3]. However, the actual exploitation of the hysteresis behavior has been hampered by the fact that the electromagnetic bistable phenomenology generally occurs at high optical intensities. The research effort aimed at reducing the bistability threshold is considerable and almost all considered strategies rely on the enhancement of nonlinear effects [4][5][6][7]. Metallo-dielectric multilayer structures have been shown to host considerable nonlinearity enhancement [8,9] and this has lead to the prediction of bistability [10,11] characterized by a low intensity threshold [12]. Optical bistability at relatively low optical intensities has also been predicted in the presence of surface plasmon polaritons [13] where the low bistability threshold is due to the strong plasmon field enhancement which in turn produces an efficient enhancement of the nonlinear effects. Following a very different route, it has recently been proposed that an extremely marked nonlinear behavior can be observed, instead of by enhancing the nonlinear response, by substantially reducing the medium linear dielectric permittivity [14] and this strategy has allowed to predict low-intensity transmissivity directional hysteresis exhibited by nonlinear metamaterial slabs with very small linear permittivity [15].
In this Letter, we show that a nonlinear metallodielectric layered slab of subwavelength thickness and very small average dielectric permittivity is characterized by a multistable behavior which is accessible at arbitrary low optical intensities. More precisely, we prove that the slab transmissivity, in addition to the angular multistable behavior occurring at low (but not arbitrary small) optical intensities (as discussed in Ref. [15]), is a multi-valued function of the incident optical intensity at each fixed incidence angle. The nonlinear electromagnetic matching at the output free-space slab interface combined with the small linear permittivity produces the multiplicity of electromagnetic states, one of which has a large electric field component normal to the slab even if the free-space transmitted intensity approaches zero. As a consequence, the lowest intensity at which transmissivity is multi-valued turns out to vanish for normal incidence and this implies that, at a very small incidence angle, multistability occurs at ultra-low optical intensities. The overall slab transmissivity turns out to be a rather articulated multi-valued function of both incident optical intensity and incidence angle and we show that ultra-low intensity multistable states can be reached only by simultaneously varying the incident optical intensity and incidence angle along suitable paths. We perform full wave simulations to investigate both the nontrivial metamaterial homogenization in the strong nonlinear regime and the role played by plasmonic resonances. We show that the induced plasmon polaritons excited at the metaldielectric interfaces affect the numerical value of the slab transmissivity solely in proximity of the hysteresis jumps without altering the overall hysteresis behavior which is therefore proven to be a feasible and robust phenomenology.
Consider the metallo-dielectric layered slab reported in Fig.1, embedded in vacuum, of thickness L = 258 nm and consisting of alternating, along the y-axis, two isotropic and non-magnetic nonlinear Kerr layers of thicknesses d 1 and d 2 . We choose to operate with a monochromatic radiation of wavelength λ = 810 nm and, if the layers thicknesses are such that d 1 + d 2 ≪ λ, the slab behaves as an homogeneous medium which, in the presence of a transverse magnetic field (TM), is characterized by the standard Kerr constitutive relations . The effective dielectric permittivity ǫ and nonlinear susceptibility χ are the average values of their underlying constituents and we choose here ǫ = −0.054 and χ = 8.72 × 10 −17 m 2 /V 2 (see below for the actual feasibility of these values). Note that the negative and very small value of the effective dielectric permittivity is achieved here by averaging the negative and positive permittivities characterizing metal and dielectric layers, respectively. As reported in Fig.1, a TM plane wave of amplitude E i is made to impinge onto the slab interface at z = 0 with incidence angle θ thus producing a TM reflected plane wave of amplitude E r and a transmitted TM plane wave of amplitude E t . We have numerically evaluated the slab transmissivity T = |E t /E i | 2 for a number of incident optical intensities and incidence angles and in Fig.2 we report T as a function of the input intensity I in = (1/2) ǫ 0 /µ 0 |E i | 2 for the incidence angles θ = 0.028 rad (panel (a)) and θ = 0.170 rad (panel (b)). The reported slab transmissivities are multi-valued functions of the incident optical intensity. The multiplicity (at last three-fold) of different slab electromagnetic states, corresponding to the same transmissivity, physically arises (as discussed in Ref. [15]) from the impact of the nonlinearity on the electromagnetic matching condition at the slab interfaces, an effect made possible by the very small value of |ǫ|. Specifically, the continuity, at z = L, of the displacement field component normal to the output slab interface, which is a cubic equation that, for a given E t and θ yields at last three different values of the real amplitude A z = (E t /|E t |) exp(−i cos θ)E z (L) characterizing the longitudinal electric field component at z = L. In Fig. 2(c) we report the surface comprising all the possible values of |χ| 1/2 A z together with values of |χ| 1/2 A z (solid lines) corresponding to the transmissivities of Fig.2(a) and (b), the capital letters labelling some reference states. Note that the incident intensities at which multistability occurs (of the order of 0.15 M W/cm 2 and 4 M W/cm 2 for the transmissivities of Fig.2(a) and 2(b), respectively) are smaller than those normally required for observing standard optical bistability. However it is evident that the smaller |ǫ| the smaller the intensity threshold required for observing the considered multistability so that very smaller intensity thresholds (of the order of W/cm 2 ) are very likely to be attained.
Multistable transmissivities of Fig.2(a) and 2(b) show two main mutual differences which deserve to be discussed. First, multi-valued transmissivity of Fig.2(a) has, for I in > I N , three different branches that does not break whereas multi-valued transmissivity of Fig.2(b) has, for I in > I E , an upper branch that does not break and two lower branches characterized by the two breaking states B and H. The origin of this difference is easily understood by considering the surface of Fig.2(c) from which it is evident that the breaking of a transmissivity branch occurs for states close to critical points at which surface folding takes place (states B and H of Fig.2(c)). Second, even though in both cases there exists a lower intensity above which the transmissivity is multi-valued (see the state N and E of Fig.2(a) and 2(b), respectively), transmissivities of Fig.2(a) and 2(b) are multistable for I in > I N = 0.04 M W/cm 2 and I in > I E = 1.92 M W/cm 2 , respectively so that I N ≪ I E or, in other words the smaller the incidence angle the smaller the incident optical intensity at which multistable states are accessible. In order to closer investigate this remarkable point we have plotted in Fig.3 the slab transmissivity as a function of both the incident optical intensity and incidence angle thus providing a complete description of the slab multistability. The transmissivity of Fig.3 is a multi-sheet surface along which we have also plotted, for clarity pur-  Fig.2(a) and 2(b) and the capital letters label the same reference states as in Fig.2

. In the inset the plane of independent parameters
Iin and θ is reported together with the region at which slab multistability occurs (shaded region).
poses, the transmissivities of Fig.2(a) and 2(b) (dashed lines) together with the very same reference states labelled with capital letters. The intersection between the upper and intermediate sheets yields the lower intensity at which multistability occurs and it turns out that such intensity is a monothonically increasing function of incidence angle θ and, most importantly, that this intensity vanishes if the angle vanishes. To underline this result, in the inset of Fig.3 we have reported the region of the independent external parameter plane (I in , θ) where multistability occurs (shaded region) and, since this region contains the excitation state I in = 0 θ = 0, it is evident that multistable electromagnetic slab states are available at arbitrary small incident optical intensities (and very small incidence angles). In order to physically grasp the origin of such an unusual situation, consider the states providing multistability and belonging to the lower sheet of the surface reported in Fig.2(c). For |E t | = 0 we have A z = 0 and A z = ± −2ǫ/(3χ) = ±0.19/ √ χ (see the above mentioned matching condition at z = L) so that the surface exists even at arbitrary small values of |E t | (see Fig.2(c)) and multistable states exist even in the limiting case where there is no transmitted waves. Since the input intensity is smaller than or equal to the output one, it is evident that the described mechanism provides multistable states even at ultra-low input optical intensities.
The transmissivity reported in Fig.2(b) clearly shows hysteresis behavior since the state B is a breaking point of its transmissivity branch. Suppose to switch on the incident plane wave with θ = 0.17 rad and I in = I A = 0.035M W/cm 2 so that the state A of Fig.2(b), 2(c) and Fig.3 is excited since before the illumination the slab was not polarized. By increasing the input intensity and holding the incidence angle fixed, the transmissivity follows the curve joining the point A and B of Fig.2(b) while the electromagnetic state continuously varies as in Fig. 2(c). From the state B, if the intensity is further increased, it is evident that the electromagnetic state undergoes a sudden jump to the state C (see Fig.2(c)) on the lower surface sheet since there is no allowed state continuously joined to B. As a consequence, the transmissivity undergoes a sudden jump to a higher value (see Fig.2(b)). If now, starting from the state C, the input intensity is decreased, the electromagnetic state varies along the curve from C to D of Fig.2(c) whereas the transmissivity in Fig.2(b) assumes the values along the curve from C to D which are different from those attained along the forward path (hysteresis). At the state D, if the intensity is further decreased, the state undergoes a jump to the state E of Fig.2(c), belonging to the central surface sheet, since the states from D to F of Fig.2(c) have, according to Fig.2(b), higher input intensities.
The transmissivity reported if Fig.2(a), although multi-valued, does not exhibit hysteresis since its branches do not present breaking points. However the low intensity multistable states can always be reached by selecting appropriate paths of the excitation plane (I in , θ). Suppose, for example, the multistability responsible state L of Fig.2(a), 2(c) and Fig.3 (of intensity I L = 0.08 M W/cm 2 ) is required to be reached. A possible operational way consists in exciting the state A and reaching the state C as discussed in the above hysteresis example. Now if from the state C the intensity is hold fixed and the angle is decreased, from Fig.3 it is evident that the state I is reached since both the states C and I belongs to the same surface sheet of Fig.2(c).
From the state C one can now fix the angle and decrease the incident optical intensity thus eventually reaching the prescribed target state L. From L one can decrease the intensity thus reaching the state M from which, a further intensity decrement, produces (in analogy to the above described hysteresis loop) a sudden jump of the trasmittivity to N of Fig.2(a), the corresponding state N of Fig.2(c) belonging to the central surface sheet. The very same procedure can even be used to excite the above discussed multistable states at much smaller incident optical intensities (and incidence angle).
In order to discuss the feasibility of the above predicted slab multistability we have performed 3D full-wave finiteelement simulations [16] for evaluating the transmissivity of the slab in the presence of a TM radiation of freespace wavelength λ = 810 nm. We have chosen d 1 = 2 nm, d 2 = 5.25 nm, L = 258 nm, ǫ 1 = −28.81 + 10i, ǫ 2 = 10.9 − 3.8i, χ 1 = 3.16 × 10 −16 m 2 /V 2 and χ 2 = 0. The parameters of medium 1 coincide with those of silver [17], characterized by a very large nonlinear susceptibility [18], with the imaginary part of the permittivity corrected by the layer size effect (since d 1 = 2 nm) [19] whereas medium 2 is a linear dielectric with gain (to compensate the metal losses). The effective permittivity and nonlinear susceptibility of the considered sample are ǫ = (d 1 ǫ 1 + d 2 ǫ 2 )/(d 1 + d 2 ) = −0.054 + 0.007i and  Fig.4(a) we report the comparison between the transmissivity evaluated through full-wave simulations for θ = 0.17 rad (dotted line) and the transmissivity of Fig.2(b). We note that good agreement exists between the results of the two kind of simulations and, most importantly, that finite-element simulations still predicts the above discussed multistability thus proving its robustness. The origin of the discrepancies between the transmissivities lies in the fact that a surface plasmon resonance occurs when the TM wave impinges on the slab as reported in Fig.1. Therefore an electric field y-component (E y ) of plasmonic origin arises mainly at the edges of the layers inside the medium and it is characterized by a subwavelength varying profile and an evanescent field tail in vacuum which does not contribute to the power flow. Evidently, if E y is much smaller than E x and E z , the homogenization theory correctly describes the slab nonlinear behavior. In Fig.4(b) we report the maximum values (within the medium) attained by the three field components for θ = 0.17 rad as a function of the incident optical intensity corresponding to the evaluated full-wave transmissivities reported in Fig.4(a). It is evident that, outside the shaded region, E y (dot-dashed line) is much smaller than both E x and E z (solid and dashed lines respectively) and in the corresponding regions of Fig.4(a) the agreement with the homogenization approach is very satisfactory. Within the shaded region of Fig.4(b), E y is not negligible and this partially breaks the validity of the homogenization approach thus leading to the discrepancies of the transmissivities in the range 1 M W/cm 2 < I in < 4 M W/cm 2 as reported in Fig.4(a).
In conclusion we have shown that a nonlinear metallodielectric layered slab of subwavelength thickness with very small average permittivity is characterized by a very low bistability threshold and, most importantly, that it provides multistable states at arbitrary low optical intensities. Such ultra-low intensity multistable states are always accessible through a suitable excitation procedure involving both the incident intensity and the incidence angle. Since the main difficulty of usefully exploiting standard optical bistability is related to the required high values of the optical intensity, we believe that our findings can have a large impact on the devising of a novel class of nanophotonic devices with complex memory functionalities and operating at very low optical intensities (as compared to those generally required in standard nonlinear optics).