Gain-assisted extraordinary optical transmission through periodic arrays of subwavelength apertures

We theoretically investigate the amplification of extraordinary optical transmission (EOT) phenomena in periodic arrays of subwavelength apertures incorporating gain media. In particular, we consider a realistic structure consisting of an opaque silver film perforated by a periodic array of slits and clad on each side by an optically pumped dielectric thin film containing rhodamine dye molecules. By solving the semiclassical electronic rate equations coupled to rigorous finite-element simulations of the electromagnetic fields, we show how the resonant electric-field enhancement associated with EOT properties enables complete ohmic loss compensation at moderate pump intensity levels. Furthermore, our calculations show that, as a consequence of the strong spatial hole-burning effects displayed by the considered structures, three separate regimes of operation arise: the system can behave as an absorber, an optical amplifier or a laser, depending on the value of the pump intensity. A discussion on the feasibility of reaching the lasing regime in the considered class of structures is also presented.


Introduction
The resonant transmission of light through subwavelength apertures in metal films has been the subject of considerable research activity over the last 10 years [1,2]. The physics underlying extraordinary optical transmission (EOT) properties relies essentially on the ability of surface plasmons (SPs) to create significant electric field enhancements at the interfaces of the metallic film. This field enhancement assists externally incident light to be funneled through subwavelength apertures, yielding much larger transmission coefficients than those expected from the classical theory of diffraction [3,4]. The key role played by SPs in this process makes, on the one hand, EOT phenomena natural candidates as the basis for a number of applications in which the subwavelength confinement of the electromagnetic (EM) field is a crucial ingredient. However, on the other hand, the large ohmic losses exhibited by SPs in the optical regime (even in the case of noble metals) represent a fundamental limitation that reduces drastically the practical applicability of EOT properties.
The above-mentioned drawback can actually be viewed as a more general challenge faced by any nanophotonic phenomena relying on the subwavelength confinement of an electric field close to a metallic surface. Recently, the incorporation of optically or electrically pumped gain media (dye molecules, quantum wells or dots, etc) was proposed as a possible route to solve the limitation introduced by ohmic losses in the performance of nanoplasmonic systems, including plasmonic waveguides and metamaterials [5][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20][21]. In this approach, the external EM cavity present in most conventional active devices (which provides the optical feedback needed for amplification or lasing action in those systems [22]) is replaced by the strong EM field confinement properties associated with SPs. More specifically, the large local electric field enhancement associated with the plasmonic resonances supported by the system induces a substantial increase of the interaction strength between the external EM radiation and the gain medium. As a consequence, the effective gain coefficient of the system is dramatically boosted (compared to the gain coefficient displayed by the same active material in bulk) to the extent that full compensation of ohmic losses was reported recently in fishnet-like metamaterials [17]. In this context, despite its significant fundamental and applied interest, a detailed analysis of the generalization of the above-described ideas to loss compensation in EOT phenomena is lacking.
The aim of this paper is to theoretically study optical amplification of EOT properties in periodic arrays of subwavelength apertures incorporating optically pumped gain media. Specifically, we consider a realistic structure formed by a periodic array of subwavelength slits perforated in an optically thick silver metallic film and sandwiched between two thin layers of a dielectric containing rhodamine dye molecules. We assume that such a structure is illuminated by means of a canonical pump-probe configuration. The theoretical framework Figure 1. Sketch of the active structure analyzed in this work. The inset shows the energy level diagram, and the corresponding electronic transitions and lifetimes, of the four-level model used to describe the interaction between external radiation and the active regions of the system. applied in this paper consists essentially in extending the traditional frequency-domain analysis of gain amplification and lasing action [22,23] to structures displaying a highly non-uniform gain distribution at the subwavelength scale. In order to do that, we generalize recent works in this context [7,11,14] and solve the electronic rate equations of a four-level model of the active medium coupled to rigorous electrodynamical calculations of both the pump and probe obtained by means of the finite-element method (FEM). Using this theoretical framework, the response of the considered class of structures is probed by analyzing how its steady-state spectral response depends on the pump beam intensity. The corresponding results are then explained in terms of the evolution with the pump intensity of the spatial distribution of the electronic population inversion in the active regions of the structure. Our numerical results show that, for realistic configurations, full compensation of ohmic losses can be achieved by using moderate pump intensities below the damage threshold of the system. Our simulations also allow us to distinguish three different regimes of operation that arise as the pump intensity is increased: the whole structure can behave as an absorber, an optical amplifier or a laser. In this paper, we also provide a discussion on how lasing action could be achieved in realistic configurations.
This paper is organized as follows. In section 2, a description of the investigated structures is presented, together with an outline of the theoretical model used in our numerical experiments. Section 3 is devoted to the analysis of loss compensation and optical amplification in the considered class of structures. Final conclusions and remarks are presented in section 4.

Theoretical framework
We consider the active structure shown in figure 1. It consists of a periodic array (period p) of subwavelength slits of width a milled on a silver metallic film of thickness w. The slit array is 4 sandwiched between two dielectric thin films of width h d containing rhodamine dye molecules, which, when optically pumped, act as a gain medium. Air regions extend infinitely above and below the whole structure. We also assume that the considered system is illuminated by means of a generic pump-probe configuration. In this configuration, part of the EM energy of a pump field (specifically, the energy carried by one of its Fourier components at a frequency ω = ω p ) is transferred to the Fourier component of a probe field at a different frequency ω = ω s (usually, ω p > ω s ) via the variation of atomic potential energy occurring inside a suitably chosen gain medium.
To accurately model the interplay between the different physical mechanisms involved in this problem, we employ a theoretical framework based on a frequency-domain approach. This kind of theoretical analysis has been extensively applied to describe optically pumped gain media both in the context of traditional active dielectric structures [22,23], as well as, more recently, in the context of gain and lasing action in plasmonic waveguides and metamaterials [7,11,14,17]. In what follows we summarize the particular numerical implementation of this framework employed in this paper.
Firstly, to describe the electronic dynamics of the rhodamine dye molecules, we assume that their corresponding electronic energy levels and transition rates can be accounted for by a four-level system (see schematics in the inset of figure 1). This approximation has shown its accuracy to model optical gain induced by the considered class of active dyes [24]. Under this assumption, the whole EM response of the structure under continuous-wave (cw) excitations can be obtained by solving the following set of linearly coupled partial differential equations, which govern the spatio-temporal dynamics of the electronic population densities (N i (r, t), i = 0, 1, 2, 3) of each of the energy levels displayed in the inset of figure 1 [22]: where τ i j (i, j = 0, 1, 2, 3) are the lifetimes of the non-radiative transitions between the ith and the jth energy levels. W a (r) and W e (r) control the radiative stimulated absorption and emission transitions in the system (i.e. they determine the exchange of energy between the gain medium and its EM surroundings) at frequencies ω a and ω e , respectively (see the definition of these frequencies in the inset of figure 1). For the pump-probe configuration considered in this work, these functions can be expressed as W a,e (r) = σ a,e (ω p,s ) 0 c n p,s 2h ω a,e E(r; ω p,s ) 2 (5) with n p,s being the value of the refractive index of the gain composite (i.e. the dielectric host and the dye molecules embedded in it) at the pump (probe) frequency n p,s = n(ω p,s ). E(r; ω p,s ) stand for the spatial profile of the Fourier component of the electric field in the structure at the pump and probe frequencies. 0 and c are the free space permittivity and the speed of light in vacuum, respectively. Finally, the functions σ a,e (ω) represent the absorption and emission cross sections of the dye. These magnitudes can be written in terms of the corresponding normalized lineshapes as σ a,e (ω) = σ (0) a,e L a,e (ω). The constants σ (0) a,e control the overall quantum yield q of the gain process q = σ (0) e /σ (0) a . Note that, owing to the fact that these lineshape functions L a,e (ω) peak at ω = ω a,e , the efficiency of the transfer of energy between the pump and probe pulses is maximal when the pump and probe central frequencies are tuned to the absorption and emission frequencies, respectively.
As described above, a crucial ingredient of the present formalism is the calculation of the EM field profile for both pump and probe wavelengths, E(r; ω p,s ). Here those fields are obtained by means of the FEM (Comsol Multiphysics package), which models Maxwell's equations exactly, except for the discretization. In these simulations, the dielectric constant used to describe the metallic regions of the system, m (ω), is given by a canonical Drude response with a Lorentzian resonance, where ∞ is the non-resonant permittivity. ω p and γ are the plasma frequency and the relaxation rate of the metal under consideration, respectively. The parameters of the Lorentz oscillator are the oscillator strength L , the resonant frequency ω L and the damping constant γ L . We now turn to the analysis of the steady-state solutions of equations (1)-(4). In general, a time-domain approach would be necessary to track the whole spatio-temporal dynamics described by that set of equations [13,18,[25][26][27][28][29]. However, if the duration of the external field pulses is much longer than the time scale needed to access the steady-state values of the electronic population densities, a cw theoretical approach can be employed to describe the optical response of the analyzed system with accurate precision. In addition, in the operation regimes considered in this paper, the physics underlying optical amplification can be fully accounted for in terms of just the steady-state behavior of the system. This also allows us to increase significantly the computational efficiency of our calculations with respect to those performed in a time-domain analysis. Specifically, by taking the limit ∂ N i (r, t)/∂t → 0 in the set of equations (1)-(4), after some straightforward algebra, one can find that the steady-state population inversion N (r) of the stimulated emission electronic transition (the one between energy levels 2 and 1) can be expressed as where N TOT is the total density of dye molecules present in the system. The functions a (r) and e (r) are defined as a (r) As can be deduced from equation (7), the dependence of N (r) on both the pump and probe E-field profiles enters through the dependence of W a (r) and W e (r) on the corresponding field intensities, |E(r; ω p )| 2 and |E(r; ω s )| 2 , respectively (see equation (5)). Now, in most active materials that can be described by a four-level system (including the active dye considered in this paper), the stimulated emission transition is orders of magnitude slower than the other electronic transitions occurring in the system, i.e. {τ 32 , τ 10 } τ 21 . This fact allows us to rewrite equation (7) in a more compact way, namely Once N (r) is computed, the complex permittivity of the active material g (r) at a frequency ω s can be obtained from the following expression [11,14,22]: where g (ω) is the relative permittivity of the considered gain composite (i.e. the dielectric material hosting the molecules and the molecules themselves). α(r, ω s ) = σ a (ω s )N 0 (r) and g(r, ω s ) = σ e (ω s )N 2 (r) are the local absorption and gain coefficients, respectively. As can be deduced from equation (9), in general, a self-consistent approach would be needed to compute the interaction of the gain medium and the pump-probe illumination. The optical pumping increases the gain coefficient of the whole structure, which varies the complex permittivity of the system. In turn, this change in g (r, ω s ) modifies the interaction between the pump and the gain material, which can modify the original assumption on how the pump and the probe were interacting with the dye molecules. However, although this is the most general scenario, for all the cases considered in this paper, this self-consistent approach can be simplified to a simpler two-step problem [11,14]. Specifically, as we show in this work, for the class of structures considered here, and for pump and probe intensities below the damage threshold of the system, one finds that σ (0) a E(r; ω p ) σ (0) e E(r; ω s ). By applying this condition to equation (8), it is straightforward to conclude that the effect of the probe beam on the population inversion N (r) is negligible, and that, consequently, the two-step approach considered here gives virtually the same results as those obtained by using a full self-consistent calculation.

Numerical results and discussion
The active medium considered in our simulations consists of a dielectric host material (of refractive index n h = 1.62) doped with rhodamine 800 (Rh800) dye molecules. The total density of dye molecules embedded in the host dielectric is assumed to be N TOT = 6 × 10 18 cm −3 [17,18]. The parameters defining the response of the dye molecules are taken from available experimental data [17,18]. Specifically, the absorption and emission cross sections of each molecule have a Lorentzian shape defined by the same half-widths, a = e = 1/(20 fs), and are centered at λ a = 680 nm and λ e = 710 nm, respectively. The values of the maxima absorption and emission cross sections are σ (0) a = 3.14 × 10 −16 cm 2 and σ (0) e = 2.43 × 10 −16 cm 2 (see the definition of these magnitudes in section 2). The stimulated emission lifetime is τ 21 = 500 ps. This value for τ 21 corresponds to a bulk configuration. The calculation of the possible position-dependent variation of τ 21 due to the presence of the metallic regions forming the considered structure represents a computational tour-de-force that is outside the scope of the present simulations. On the other hand, the metallic regions of the considered system are assumed to be made of silver. Thus, the following values are used in the expression for dispersion relation given in equation (6): ∞ = 4.6, ω p = 9.0 eV, γ = 0.07 eV, ω L = 4.9 eV, γ L = 1.2 eV and L = 1.10 [18,30]. In the rest of this paper, we refer to the structure in which the total density of dye molecules embedded in the host dielectric layers is negligible as the passive structure.
The computational domain used in our FEM simulations is bounded by perfectly matched layers along the z-direction (see the definition of axes in figure 1), whereas periodic boundary conditions are applied along the x-direction. The defined computational window is then sampled by a non-uniform mesh whose maximum element size is more than 15 times smaller than the operation wavelength. In addition, we note that the transmission and reflection spectra shown in this paper were obtained by integrating the flow of the Poynting vector through line monitors properly placed below and above the air-dielectric interfaces and parallel to them. The corresponding absorption spectrum is then obtained just by subtracting, at each wavelength, transmission and reflection from unity. Figure 2 renders the results from a first set of simulations in which the spectral response of the passive case is computed for several values of the thickness h d of the dielectric layers cladding the metallic array. The geometrical parameters defining the metallic slit array considered in these simulations are w = 150 nm, p = 520 nm and a = 100 nm (see the definition of these parameters in figure 1). As we show below, these geometrical parameters are optimized so that the main SP resonance supported by slit array is tuned to the stimulated emission wavelength (λ e = 710 nm). In addition, to maximize the whole optical amplification process, all simulations are performed assuming that the pump and probe signals are centered at λ a and λ e , respectively. This optimization of the geometrical parameters of the structure can be viewed as a particular implementation of the above-mentioned general approach to loss compensation in nanophotonics, consisting in controlling the optical response of a structure by tuning its plasmonic resonance to an electronic radiative transition of a suitably chosen class of molecules embedded in the system.
As shown in figure 2, the spectra computed for each value of h d are dominated by the presence of three main resonances. The physical origin of these resonant features emerges  [31]. In contrast, as revealed by figure 3(b), the resonance located at λ ≈ 710 nm confines most of the field near the metal-dielectric interfaces, and thereby its physical origin can be linked to the presence of SPs decorating the interfaces of the slit array [31]. This fact makes this latter resonance the optimal candidate for enhancing the interaction between the external radiation and the dye molecules embedded in the considered class of structures. Moreover, as observed in figure 2, the spectral location of this resonance is rather sensitive to small variations of h d , which allows us to easily tune the resonance to the stimulated electronic transition of the gain medium. For the specific structure and gain medium considered in this paper, the optimal thickness of the dielectric layers is h d = 200 nm (see figure 2). Thus, in the rest of the paper we assume this value for h d . We point out that, although it is outside the framework of the fundamental analysis presented here, further enhancement of the interaction between external radiation and the active medium could be achieved by optimizing the structure so that λ a and λ e are simultaneously tuned to two SP resonances of the system. Figure 4 summarizes our numerical results on how the spectral response of the considered structure is modified when optically pumped dye molecules are embedded in the top and bottom dielectric layers sketched in figure 1. In these simulations, we assume the same geometrical parameters as those given above for the passive system. Note also that all the electric-field magnitudes are normalized to the magnitude of the saturation electric field, defined as |E th | = 2hω a /( 0 cn τ 21 σ (0) a ) (whereh is the reduced Planck constant and n is the refractive index of gain composite, n = ( g ) 1/2 ). For a given set of parameters describing the electronic four-level system, the value of |E th | corresponds to the magnitude of the pump electric field (|E in |) at which the electronic population of the ground is depleted to one half of its initial value (i.e. at which N 0 = 0.5 × N TOT ). Figures 4(a)-(c) clearly demonstrate that the system displays two notably different responses as a function of the optical pumping strength. For values of |E in | such that |E in | < |E th |, the whole structure behaves essentially as an optical absorber centered at the stimulated emission wavelength. This is clearly visualized in the absorption peak centered at λ e observed in figure 4(c) for the case |E in |/|E th | = 0.1. Note that the resonant peak located at λ ≈ 620 nm is barely influenced by the optical pumping, essentially because the corresponding wavelength is shifted away from the center of the gain profile. As the value of |E in | is increased, the absorption peak starts smearing out (see figure 4(c)), until the characteristics of the transmission and reflection spectra of the passive system are fully recovered by the active one. We have found numerically that, in the considered structure, this occurs when |E in |/|E th | ≈ 0.9 (see the insets of figures 4(a) and (c)). Now, as shown in figure 4(a), if the value of |E in | is further increased, the EM power generated through the pump-probe frequency conversion taking place inside the gain medium starts compensating for the intrinsic ohmic losses present in the system (i.e. the absorption losses due to the penetration of the EM field into the metallic regions). The loss-compensation process flattens the absorption peak centered at λ e (see the main and inset panels in figure 4(c)), and leads to values of the maximum transmission at resonances larger than those found in the passive structure (see the inset of figure 4(a)). This dependence of the spectral response on |E in | continuously evolves as the pump field is further increased, up to a value of |E in | at which the optical pumping is large enough to completely compensate for the ohmic losses present in the structure. In the analyzed system, this condition is reached when |E in |/|E th | ≈ 1.9.
Importantly, when complete ohmic loss compensation occurs, 100% transmission is achieved at λ = λ e . This result is linked to the full transmission properties observed in EOT phenomena occurring in metals that behave as perfect metals (i.e. in which Im( m ) = 0) [4,32]. To investigate this point in more detail, we have compared the numerical results obtained for the active structure operating at |E in |/|E th | ≈ 1.9 with those corresponding to a passive structure in which the condition Im( m ) = 0 is artificially imposed. Figure 5 shows the comparison between the transmission and reflection spectra obtained for each case. The excellent agreement found between both calculations demonstrates that, owing to the action of the gain medium, the system does effectively act as a passive system in which the silver regions are replaced by perfect metal areas.
For values of the pump field magnitude larger than those corresponding to complete ohmic loss compensation, the action of the gain medium induces the considered system to effectively behave as an optical amplifier [23]. This is clearly illustrated in figure 4(c), which shows that for |E in | > 1.9|E th |, negative values of absorption are obtained in a wavelength interval centered at λ = λ e . Note that, due to the physical limit set by the damage threshold of the structure, values of |E in | larger than 5|E th | would lead to unrealistic predictions [14]. As displayed in figure 4(a), in this optical amplification regime, the values of maximum transmission at resonance can be significantly larger than unity. As discussed below, this fact opens up the possibility of inducing lasing action in the system. We also note that, by comparing figures 4(a) and (b), it can be deduced that, within this optical amplification regime, most of the optical power generated in the gain medium ends up being redirected to the transmission region rather than to the reflection region. This could also be an important factor that facilitates the practical implementation of the results discussed in this paper.
Further physical insight into this problem can be obtained by analyzing how the evolution with the pump intensity of the effective response of the system can be accounted for in terms of the corresponding variation of the population inversion spatial profile N (r). Figure 6 renders the results for spatial distribution of the population inversion, normalized to N TOT (panels (a), (c) and (e)), and the corresponding imaginary part of the dielectric constant inside the dielectric layers (right column, panels (b), (d) and (f)). Three different values of the pump electricfield magnitude are considered in this figure. Figures 6(a) and (b) correspond to |E in |/|E th | = 0.1, figures 6(c) and (d) display the results for |E in |/|E th | = 1, while in figures 6(e) and (f) |E in |/|E th | = 5 is assumed. These values for |E in |/|E th | correspond to three representative spectral responses of the system (see the red, green and blue lines in the spectra shown in figure 4). All the results displayed in figure 6 have been computed at the stimulated emission wavelength, λ e = 710 nm. Figure 6(a) shows how, for the lowest magnitude of the incident electric field, the density of inverted dye molecules is negligible in nearly every position inside the active regions of the structure. As a consequence, these active regions act effectively as an almost uniform absorbing material with its resonance centered at λ = λ e (see figure 6(b)). In contrast, for |E in |/|E th | = 1, the response of the system changes substantially, and the active dielectric regions display some hot-spots in which a significant density of inverted molecules is found (see figure 6(c)). These regions are mainly located near the metallic interfaces of the slit array, basically because it is precisely in those regions where significant values of electric-field enhancement are reached. The corresponding spatial distribution of g shows that in these regions of maximal population inversion, the sign of the absorption coefficient changes from positive to negative (see figure 6(d)), i.e. the active medium acts locally as the gain medium in those regions. The effective response of active dielectric layers is then determined by the spatial average (weighted by the electric-field intensity) of the contribution to the absorption coefficient of the absorbing and gain regions. In particular, for |E in |/|E th | = 1 this average is approximately zero, which explains why in this case the spectral response of the structures coincides with that of the passive system (see figure 4). For the largest value of |E in | considered here (|E in |/|E th | = 5), the optical pumping is so strong that the regions displaying significant population inversion occupy almost completely the whole volume defined by the top and bottom dielectric layers. Thus, as seen in figure 6(e), in this operation regime, these layers can be effectively considered as formed by a uniform optical amplifying medium. Importantly, note the crucial role played by the high non-uniformity of EM field distributions in determining the effective response of the system. Specifically, as discussed above, this non-uniformity induces the presence of spatial holes in the gain medium. In these spatial holes, the local population inversion is not saturated, or alternatively, the population inversion density is smaller than in the regions of maximum electric field intensity. This occurs much in the same way as in the well-known spatial holeburning processes observed in traditional active and lasing systems [22,23]. Thus, we believe that the results reported in this paper can be considered as a novel instance of how spatial holeburning effects at the subwavelength scale can be responsible for the effective response of an active nanophotonic structure.
Finally, to quantitatively characterize the different operation regimes described above, we analyze how the total decay rate of the EM energy stored in the considered active structure can be decomposed into the different absorption and gain channels governing the effective response of the system. As we discuss below, this analysis also allows us to investigate the conditions under which the studied systems can display lasing action at the subwavelength scale. Specifically, the total decay rate of the structure ( T ) can be expressed as a function of the pump E-field magnitude |E in | as where rad and Ohm are the decay rates corresponding to the intrinsic radiation and ohmic losses present in the system, respectively. Rh800 is the decay rate associated with the absorption of the dye molecules in the absence of external pumping. Note that ohm , rad and Rh800 are only determined, respectively, by material properties of the slit array (ohmic losses), by the geometry of the system (radiation losses) and by intrinsic properties of the rhodamine molecules (the linewidth and maximum value of the corresponding absorption cross section). Therefore, these three magnitudes are independent of |E in |. On the other hand, the function g (|E in |) determines the rate at which the external EM energy is pumped into the system (thereby it has the opposite sign from rad , ohm and Rh800 in equation (10)); this magnitude obviously satisfies g (|E in | = 0) = 0. The value of rad can be easily obtained from the linewidth of the resonant transmission peak observed at λ = λ e in the case of the passive structure in which the condition Im( m ) = 0 is imposed (see figure 5). Once the value of rad is known, the calculation of Ohm is also straightforward from the linewidth of the resonant transmission peak computed for the passive structure, now assuming a realistic m in the metallic regions of the system (see the black line in figure 4). Finally, for fixed values of rad and Ohm , the value of Rh800 can be readily deduced from the linewidth of the active system assuming no optical pumping (i.e. assuming |E in | = 0). Thus, once the values of rad , Ohm and Rh800 have been computed, the dependence of g on |E in | can be obtained just by tracking how the linewidth of transmission resonance observed in the active structure evolves as |E in | is increased. Figure 7 renders the corresponding results for the described absorption and gain rates. In this figure, the decay rates are normalized to the total decay rate of the active system with no external pumping T (|E in | = 0), which is given by T (|E in | = 0) = rad + Ohm + Rh800 (see equation (10)). The horizontal red line in figure 7 marks the value of the normalized molecular absorption losses, Rh800 / (|E in | = 0), while the green line corresponds to the normalized sum of all intrinsic absorption losses of the system ( Ohm + Rh800 )/ (|E in = 0|). Blue dots represent the calculated values of g / (|E in = 0|) for the corresponding values of |E in |/|E th | displayed in the figure.
As illustrated in figure 7, three different intervals of operation (labeled as A, B and C) emerge from this analysis. For the values of |E in |/|E th | < 0.8 (yellow region in figure 7), we obtain g < Rh800 . This interval corresponds to the above-described low-pumping regime, in which the gain medium effectively behaves as an absorbing material. For values of |E in | such that 0.8 < |E in |/|E th | < 1.6 (cyan region in figure 7), g becomes larger than Rh800 and, consequently, it contributes to reduce the total linewidth of the considered transmission resonance (see equation (10) and the inset of figure 7). From a physical point of view, these results suggest that, for low pump intensities, the EM energy from the active medium is first invested in compensating for the absorption losses associated with the dye molecules (interval A in figure 7). Once those losses have been compensated for, the pump-probe frequency conversion process taking place inside the gain medium starts compensating for the ohmic losses stemming from the penetration of the EM in the metallic regions of the system (this regime is labeled B in figure 7). Following this reasoning, one expects that once all intrinsic absorption losses have been compensated for, further increase of |E in | would lead to a new regime in which the energy provided by the gain medium starts compensating for the leakage of energy due to the radiation losses. This is precisely what is observed in the interval C of figure 7. We note that the values of |E in |/|E th | at which the crossover between the different regimes of operation takes place agree well with the predictions for the same values obtained previously from the analysis of the evolution with |E in |/|E th | of the maximum transmission at resonance (see figure 4). From figure 7, it is also clear that the value of g / (|E in | = 0) saturates far from unity. This is essentially due to the fact that the saturation of the gain occurs before all the losses present in the system are fully compensated. Importantly, this full compensation of losses actually determines the onset of lasing action in the considered structure (see the thick light blue line in figure 7), similarly to what occurs in the different plasmonic structures displaying lasing action reported in the recent literature [33][34][35][36]. In this context, in order to access the lasing regime, one can envision a system in which the radiation losses are dramatically decreased by, for instance, reducing the width of the apertures. Alternatively, the lasing region could be accessible to a configuration with a similar geometry as that shown in this paper, but in which the damage threshold is significantly increased and, in which, in addition, the dye molecules have been replaced by a different gain medium showing notably larger saturation electric field intensities.

Conclusions
In conclusion, we have theoretically studied loss compensation in EOT phenomena enabled by optically pumped active media. Specifically, we have analyzed the spectral response of an active structure formed by a periodic array of subwavelength slits milled in an optically thick silver metallic film and sandwiched between two thin dielectric layers including rhodamine dye molecules. The corresponding results have been accounted for in terms of the evolution with the pump field magnitude of the spatial distribution of the density of electronically inverted molecules present in the system. Our calculations have shown that, for a conventional pump-probe illumination, moderate values of the magnitude of the pump electric field can lead to full compensation of the absorption ohmic losses induced by the penetration of the field in the metallic regions. We have also reported on the existence of three main operation regimes for the analyzed system. For low values of pump field, the whole structure effectively behaves as an optical absorber centered at the stimulated emission wavelength. As the pump field is increased, the system starts behaving as a special class of optical amplifiers in which the EM energy provided by the gain medium is first invested in compensating for the ohmic losses. Once the ohmic losses are fully compensated, the energy provided by the gain-assisted pump-probe frequency conversion process allows compensating also the radiation losses present in the structure. We have also quantitatively characterized the conditions that would lead to full compensation of both ohmic and radiation losses and therefore to the onset of lasing action in the considered class of structures.