We demonstrate that when the frequency of the external field differs from the lasing frequency of an autonomous spaser, the spaser exhibits stochastic oscillations at low field intensity. The plasmon oscillations lock to the frequency of the external field only when the field amplitude exceeds a threshold value. We find a region of values of the external field amplitude and the frequency detuning (the Arnold tongue) for which the spaser synchronizes with the external wave.
© 2011 OSA
Joule (ohmic) losses in metallic nano-inclusions embedded in a host metamaterial medium are the main obstacle to realizing a practical “perfect” lens with subwavelength resolution as well as to other applications of metamaterials . Losses can be compensated for, however, in amplifying media. The authors of [2,3] proposed making a stackable lens with alternating metamaterial and gain layers This can be accomplished by directly introducing gain inclusions such as molecules or quantum dots (QDs) into the matrix [4–6]. The latter scheme has recently been realized experimentally [7–10] demonstrating that it is possible to compensate for losses.
The combination of a gain medium and metallic nanoparticles (NPs) results in the emergence of a spaser (Surface Plasmon Amplification by Stimulated Emission of Radiation) first suggested by Bergman and Stockman  and realized experimentally in . Spasers, as well as nanoplasmonics overall, is a new fast growing field with tremendous promise. Schematically, the spaser is a system of inversely excited two-level QDs surrounding metal nanoparticles [11,13]. Its principles of operation are analogous to those for a laser with the role of photons played by surface plasmons (SPs) localized at a NP that serves as the resonator [11,14,15]. In other words, in a spaser, near-fields of the NP are generated and amplified. The amplification of the SPs occurs due to non-radiative energy transfer from QDs. This process takes place due to the dipole-dipole (or any other near-field ) interaction between the QD and the plasmon NP. The efficiency of such a mechanism depends on the probability of the non-radiative excitation of the surface plasmon, which is greater than the radiation of the photon , where is the distance between the centers of the QD and the NP and k is the photon wavenumber. Stimulated radiation from the QD into the plasmon mode results in spasing. The excitation of the plasmon mode is carried out via the excitation of the QD.
In theoretical studies of a metamaterial with compensated losses, it is assumed that electromagnetic wave propagation can be described by the Maxwell equations with real valued effective electric permittivity and magnetic permeability. In its turn, the use of effective permittivities and refraction coefficients assumes that the dipole moment of a NP oscillates with the frequency of the external field and the amplitude of this oscillation is determined by the external field. In the absence of the external field, there should be no oscillations of a dipole moment. The classical linear description treats the gain medium as a system with a negative imaginary part of the permittivity [18–21].
However, when spasers are used as gain inclusions, their dipole moments are excited not only by the external field but also by the radiation produced by QDs, which actually tends to compensate for loss. Unlike the classical linear description [18–21], the semi-classical analysis shows that the spaser in the presence of pumping is a self-oscillating system. The dipole moment of the spaser’s NP oscillates autonomically even in the absence of the external field. The autonomic frequency of this self-oscillation is determined by the plasmon frequency, transition frequency of the gain inclusion and characteristic times of relaxation in NP and excitation of gain inclusions . Therefore, in order to develop a correct description of metamaterials with spasers, it is necessary to study the interaction of the spaser with the external field in details.
In this paper, we study analytically and numerically the operation of a spaser driven by an external optical wave. We demonstrate that the pumping drastically changes the spaser’s behavior in the optical field in comparison with the behavior of a passive QD-NP pair. In particular, when the frequency of the external field is detuned from the autonomous frequency of the spaser, even infinitesimally weak field drives the spaser to stochastic oscillations. The spaser can be synchronized to the external field only when the field amplitude exceeds a threshold value.
2. Equations of motion of the spaser
The simplest model of the spaser consists of a two-level QD of size which is positioned at a distance r from a metallic NP of size . The whole system is immersed into a solid dielectric or semiconductor matrix with dielectric permittivity . The dipole moment of a typical QD of size 10 nm is 20 Debye . For a NP with Ag-core and SiO2-shell with radius , the classical dipole moment near the frequency of the plasmon resonance, ω, can be estimated to be . The energy of the NP-QD quasistatic dipole-dipole interaction is estimated to be which gives Thus, the Rabi frequency (the coupling constant) is two orders of magnitude smaller than the generation frequency. This allows for the use of the slowly varying envelope approximation in the investigation of the dynamics of the spaser.
The transition processes of the NP and two-level QD interacting in a spaser can be described by the model Hamiltonian11,23,24], the operator determines their interaction. To take into account the interaction of the spaser, which is an open system, with surroundings (pumping, phonons in metal, etc.) following  we introduce the non-Hermitian operator that describes relaxation and pumping processes. Here is the annihilation operator of the dipole SP, is the operator for the dipole moment of the QD, is the transition operator between ground and excited states of the QD, is the QD dipole moment matrix element. Quantization of the plasmon field without taking into account losses is carried out in the standard way . Thus, we obtain [11,15]11]. In the dipole approximation, in Eq. (4) we retain only one dipole mode () with the frequency. For a spherical NP this dipole mode is uniform inside the particle and has the form of the field of the dipole with a unitary dipole moment, , so that and outside and inside of the NP, respectively. Note, that even though is a unit vector, it is a dimensional quantity, therefore, below we will show explicitly.
In the absence of the external field, for the NP near field we have . The near field is not equal to zero at the resonance frequency only. However, in this case, the field vanishes at the infinity and . The latter integral is equal to zero because both inside and outside the NP. Furthermore, outside of the particle. This leads to the following result
Assuming that the frequencies of the QD transition and the frequency of the dipole SP are close, , we will look for the time dependencies of and in the form and , where and are the slowly changing amplitudes and is the autonomous frequency of the spaser which we seek. In the rotating wave approximation  we can neglect rapidly oscillating terms, , and obtain the interaction operator in the form of the Jaynes–Cummings Hamiltonian 11,28]).
The commutation relations for operators and are standard: and , where the operator describes the inversion of the occupancies of the ground and excited states of the QD,, and are the occupancy operators for the ground and excited states, . Using the Hamiltonian, Eqs. (1)-(3), and the commutation relations for operators and we obtain the Heisenberg equations of motion for the operators , , and [28,29]27]. The operator for the occupancy inversion describes the pumping [24,27]. We now neglect quantum fluctuations and correlations and consider , and as complex quantities (c-numbers), so that we can use complex conjugation instead of Hermitian conjugation [11,28,30,31]. The difference between the occupancies of the upper and lower levels D(t) must be a real valued quantity because the respective operator is Hermitian. The quantities σ(t) and a(t) are the complex amplitudes of the dipole oscillations of the QD and SP, respectively.
For the system of Eqs. (8)-(10) has a trivial stationary solution only, while for there are two stationary solutions. The trivial unstable solution corresponds to the absence of SPs, while the stable one corresponds to laser generation:Fig. 1 ). The stationary values of a, σ, and D shown in Fig. 1.
3. The spaser in the external electromagnetic field
Let us consider the dynamics of the NP and QD in the field of the external optical wave, . Assuming that the external field is classical and taking into account the dipole interaction only, we can write the Hamiltonian in the formEq. (2), , and . As before, we use the Heisenberg equations for operators , , and to obtain equations of motion for “slowly varying” amplitudes,30]).
The quantities defined by Eqs. (13)-(15) differ from those in Eqs. (9)-(11). Without the external field the frequency in the substitutions and is found from the solution of stationary Eqs. (8)-(10). When the external field is present, we use the representations and . We look for a stationary solution of the system of Eqs. (13)-(15). The stationary values of D satisfy a cubic equation and are shown in Fig. 2 as a function of the amplitude of the external field E. For (Fig. 2а) in the absence of the external field and two points, stable and unstable, correspond to the points marked by dots in Fig. 1.
4. Spaser synchronization
Only the lower branch of the synchronization region corresponds to a stable solution. The stationary solutions at zero external field and with no detuning () are shown in Fig. 1. When detuning is present (we still assume that but ), the stable stationary solution exists for only (Fig. 2b). This result is confirmed by our numerical simulation shown in Fig. 3 . Such behavior of a self-oscillating system is referred to as synchronization by an external periodic influence . The region in which the synchronization takes place is known as the Arnold tongue. To obtain an analytical estimate of the Arnold tongue boundary, we find the solution of Eqs. (13)-(15) in the first approximation with respect to the field. For the case , using the substitution and in Eq. (15) and equating its imaginary parts, we obtain
Since , represents the phase shift between the external field and spaser oscillations. For , , and we can use their values obtained from Eq. (11) in the zeroth order with respect to the field approximation. As the result, we haveEq. (18) takes the form of Adler’s equation :
The phase dynamics can be viewed as the motion of a particle sliding along the potential profile (Fig. 4 ) in a viscous liquid. For a small field and/or large detuning, , the phase difference of the system and the external field increases monotonically. For , the “particle” should be trapped in one of the minima of the potential function. In this regime the oscillation of the system is synchronized and the phase difference is time independent. In the approximation considered, the width of the synchronization region is proportional to the field amplitude. Within the Arnold tongue, if the energy of the interaction of the NP with the external field is smaller than the QD-NP interaction, the amplitude of the auto-oscillations depends weakly upon the amplitude of the external field, which in this case plays the role of a synchronizer.
When the threshold field value is determined by the equation For large values of the detuning, becomes independent of (see Fig. 3). Taking into account that and for plasmonic NPs, one can estimate the asymptotic value of the Arnold tongue boundary, . Our numerical calculations show that tends to a plateau for , which gives . This agrees with calculations of [3,10] in which wave propagation in the system of spasers was considered. In [3,10] the amplitude of the incident wave was a few orders of magnitude larger than . Such a field is comparable with the near field inside the spaser. In this case, spaser becomes synchronized with the external field for any value of detuning; it ceases to be an autonomous system and responds linearly to the external field as can be seen in Fig. 3. For such strong fields, losses cannot be compensated exactly. Moreover, our numerical calculations show that for a very strong field, the population inversion of the QD decreases, the excitation of SPs by QDs is inhibited, and the spaser becomes just a lossy NP. Thus, when detuning is substantial, spasers cannot be used as active inclusions for the perfect lens because both loss and amplification destroy the perfect image, so that the exact compensation of loss is required.
For applications it is important to know not only a response of a single spaser to an external electromagnetic wave but also the behavior of a system of interacting spasers randomly or regularly distributed in the dielectric matrix. The interaction between spasers may lead to a cooperative behavior of spasers. This problem is a subject of our future research.
The authors are indebted to Yu. E. Lozovik and C. Z. Ning for useful discussions. This work was supported by RFBR Grants Nos. 10-02-91750 and 11-02-92475 and by a PSC-CUNY grant.
References and links
1. W. Cai and V. Shalaev, Optical Metamaterials (Springer, 2010).
2. S. A. Ramakrishna and J. B. Pendry, “Removal of absorption and increase in resolution in a near-field lens via optical gain,” Phys. Rev. B 67(20), 201101 (2003). [CrossRef]
3. A. Fang, T. Koschny, and C. M. Soukoulis, “Lasing in metamaterial nanostructures,” J. Opt. 12(2), 024013 (2010). [CrossRef]
4. M. A. Noginov, G. Zhu, M. Bahoura, J. Adegoke, C. E. Small, B. A. Ritzo, V. P. Drachev, and V. M. Shalaev, “Enhancement of surface plasmons in an Ag aggregate by optical gain in a dielectric medium,” Opt. Lett. 31(20), 3022–3024 (2006). [CrossRef]
6. A. K. Sarychev and G. Tartakovsky, “Magnetic plasmonic metamaterials in actively pumped host medium and plasmonic nanolaser,” Phys. Rev. B 75(8), 085436 (2007). [CrossRef]
7. M. A. Noginov, G. Zhu, V. P. Drachev, and V. M. Shalaev, “Surface plasmons and gain media,” in Nanophotonics with Surface Plasmons, V. M. Shalaev and S. Kawata, eds. (Elsevier, 2007).
8. N. I. Zheludev, S. L. Prosvirnin, N. Papasimakis, and V. A. Fedotov, “Lasing spaser,” Nat. Photonics 2(6), 351–354 (2008). [CrossRef]
9. E. Plum, V. A. Fedotov, P. Kuo, D. P. Tsai, and N. I. Zheludev, “Towards the lasing spaser: controlling metamaterial optical response with semiconductor quantum dots,” Opt. Express 17(10), 8548–8551 (2009), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-17-10-8548. [CrossRef]
10. N. Meinzer, M. Ruther, S. Linden, C. M. Soukoulis, G. Khitrova, J. Hendrickson, J. D. Olitzky, H. M. Gibbs, and M. Wegener, “Arrays of Ag split-ring resonators coupled to InGaAs single-quantum-well gain,” Opt. Express 18(23), 24140–24151 (2010), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-18-23-24140. [CrossRef]
11. D. J. Bergman and M. I. Stockman, “Surface plasmon amplification by stimulated emission of radiation: quantum generation of coherent surface plasmons in nanosystems,” Phys. Rev. Lett. 90(2), 027402 (2003). [CrossRef]
12. M. A. Noginov, G. Zhu, A. M. Belgrave, R. Bakker, V. M. Shalaev, E. E. Narimanov, S. Stout, E. Herz, T. Suteewong, and U. Wiesner, “Demonstration of a spaser-based nanolaser,” Nature 460(7259), 1110–1112 (2009). [CrossRef]
13. M. I. Stockman, “Spasers explained,” Nat. Photonics 2(6), 327–329 (2008). [CrossRef]
14. M. I. Stockman, “Ultrafast nanoplasmonics under coherent control,” New J. Phys. 10(2), 025031 (2008). [CrossRef]
15. M. I. Stockman, “The spaser as a nanoscale quantum generator and ultrafast amplifier,” J. Opt. 12(2), 024004 (2010). [CrossRef]
16. I. R. Gabitov, B. Kennedy, and A. I. Maimistov, “Coherent amplification of optical pulses in metamaterials,” IEEE J. Sel. Top. Quantum Electron. 16(2), 401–409 (2010). [CrossRef]
17. V. V. Klimov, “Nanoplasmonics,” Phys. Usp. 51(8), 839–844 (2008). [CrossRef]
18. L.-M. Si, T. Jiang, K. Chang, T.-C. Chen, X. Lv, L. Ran, and H. Xin, “Active microwave metamaterials incorporating ideal gain devices,” Materials 4(1), 73–83 (2011). [CrossRef]
19. A. D. Boardman, Y. G. Rapoport, N. King, and V. N. Malnev, “Creating stable gain in active metamaterials,” J. Opt. Soc. Am. B 24(10), A53–A61 (2007). [CrossRef]
20. Y.-Y. Yu, D.-Z. Lin, L.-S. Huang, and C.-K. Lee, “Effect of subwavelength annular aperture diameter on the nondiffracting region of generated Bessel beams,” Opt. Express 17(4), 2707–2713 (2009), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-17-4-2707. [CrossRef]
21. A. Govyadinov, V. Podolskiy, and M. A. Noginov, “Active metamaterials: Sign of refractive index and gain-assisted dispersion management,” Appl. Phys. Lett. 91(19), 191103 (2007). [CrossRef]
22. A. Muller, Q. Q. Wang, P. Bianucci, C. K. Shih, and Q. K. Xue, “Determination of anisotropic dipole moments in self-assembled quantum dots using Rabi oscillations,” Appl. Phys. Lett. 84(6), 981–983 (2004). [CrossRef]
23. L. Novotny and B. Hecht, Principles of Nano-Optics (Cambridge University Press, 2006).
24. R. H. Pantell and H. E. Puthoff, Fundamentals of Quantum Electronics (Wiley, 1969).
25. P. W. Milonni, The Quantum Vacuum: An Introduction to Quantum Electrodynamics (Academic Press, 1994).
27. M. O. Scully and M. S. Zubairy, Quantum Optics (Cambridge University Press, Cambridge, 1997).
28. I. E. Protsenko, A. V. Uskov, O. A. Zaimidoroga, V. N. Samoilov, and E. P. O’Reilly, “Dipole nanolaser,” Phys. Rev. A 71(6), 063812 (2005). [CrossRef]
29. K. Li, X. Li, M. I. Stockman, and D. J. Bergman, “Surface plasmon amplification by stimulated emission in nanolenses,” Phys. Rev. B 71(11), 115409 (2005). [CrossRef]
30. A. S. Rosenthal and T. Ghannam, “Dipole nanolasers: A study of their quantum properties,” Phys. Rev. A 79(4), 043824 (2009). [CrossRef]
31. A. N. Lagarkov, A. K. Sarychev, V. N. Kissel, and G. Tartakovsky, “Superresolution and enhancement in metamaterials,” Phys. Usp. 52(9), 959–967 (2009). [CrossRef]
32. A. Pikovsky, M. Rosenblum, and J. Kurths, Synchronization. A Universal Concept in Nonlinear Sciences (Cambridge University Press, 2001).