## Abstract

We consider emitting nanoparticles in dielectric nanocomposites with varying refractive index contrast and geometry. For that we develop a simple and universal method to calculate the emission enhancement in nanocomposites that employs solely the calculation of the effective refractive index and electric field distributions from three quasistatic calculations with orthogonal polarizations. The method is exemplified for dilute nanocomposites without electromagnetic interaction between emitting particles as well as for dense nanocomposites with strong particle interaction. We show that the radiative decay in dielectric nanocomposites is greatly affected by the shape and arrangement of its constituents and give guidelines for larger enhancement.

© 2018 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

## 1. Introduction

The radiative lifetime of the emitter influences the efficiency of a light source as there are always competing non-radiative processes present which dissipate energy [1]. The shorter the radiative lifetime is compared to other decay processes, the more likely the emitter’s energy is transformed into light and not into heat. Besides the improvement of the efficiency of light sources, a fast radiative decay is also important for fast switching as the radiative decay rate sets the limit for the maximal modulation rate of a light source [2].

The radiative decay rate of an emitter, such as an erbium atom, is not only dependent on the emitter itself but can also be changed by modifying its environment. Strongly increased decay rates originate from a strong electromagnetic interaction between the luminescent medium and the optical field. Emission enhancement can be achieved in optical devices with local field enhancement such as resonators [3], slow light structures [4], or by focusing the optical fields into a small volume [5]. Of particular interest is the emission enhancement in nanocomposite media employing emitting nanoparticles and a host medium [6,7]. Typical considerations of this enhancement are limited to diluted spherical and elliptical inclusions [7,8] or flat slots [9] of an emitting medium in a high index host. There is no simple description available for dense nanocomposites of arbitrary geometry.

Commonly, the optical properties of dense dielectric nanocomposites are described by effective medium models [10]. These models fall short if one tries to calculate the optical emission as they can only describe the average electromagnetic fields. However, for emission processes local fields cannot be neglected. Here, we derive an equation for the average emission enhancement in dielectric nanocomposites based on the reciprocity relation and a quasistatic approximation. We discuss the enhancement for two shapes of luminescent particles in diluted composites as limiting cases and present enhancement calculations for two representative arrangements of the luminescent medium in dense composites using our approach. .

## 2. Calculating the emission of a dipole via reciprocity relation

Although the emission of a photon is a quantum mechanical process, we can accurately model the change of emission properties upon change of the environment by a classical dipole. The lifetime of a quantum emitter placed in vacuum divided by the lifetime inside the environment, ${f}_{e}=\frac{{\tau}_{vac}}{{\tau}_{enviroment}}$, is identical to the ratio of the power a continuously oscillating dipole emits inside the environment to the power emitted in vacuum, ${f}_{e}=\frac{{P}_{enviroment}}{{P}_{vac}}$ [11]. Thus, we will use the latter ratio as our metric for the emission enhancement.

The methods available to describe the emission of dipoles in structured media are usually based on Green’s function approach [12]. For a composite, this approach requires complex numerical calculations for every point in the volume. In nanocomposite the structure parameters are usually much smaller than the optical wavelength. Here, we make use of this fact and utilize a simple quasistatic approximation in combination with the reciprocity relation. This consideration allows us to obtain the average emission enhancement for a distributed luminescent medium from three quasistatic calculations with orthogonal polarizations.

The investigated structure is sketched in Fig. 1(a). It is a nanocomposite composed of an emitting material (drawn in red) inside a host material. We start by considering only a single emitter inside this nanocomposite. We model this emitter as an electric dipole. The power emitted by the dipole can be evaluated by integrating the intensity on the closed surface in the far field.

To calculate the far field emission of the source, we can simplify the problem by dividing the emission process into two parts: Initially, the dipole emits a wave and subsequently the wave propagates to our defined surface in the far field region. The wavelength of the propagating field is assumed to be much larger than the characteristic length of our nanocomposite and we can describe the wave propagation in the far field accurately with effective material parameters. We cannot do that for the dipole emission as this process is governed by the fields exactly at the position of the dipole. Therefore, we artificially divide the composite in two regions: In the outer region, we consider a macroscopically isotropic medium. In the inner region, we model the exact shape of the nanocomposite to obtain the electric field at the dipole position. The placement of the boundary was chosen such that it does not affect the fields at the emitter. The region outside has the same refractive index as the average refractive index of the nanocomposite, thus there will be no reflections at this boundary. Our goal is now to obtain the far fields emitted by the dipole such that we can quantify the emitted power.

We start by calculating the field, ${E}_{1}$, that the dipole produces far away from the source (see Fig. 1). An easy way to calculate the far field is to place a test dipole at the same positon and with the same orientation as the field we are interested in (Fig. 1(c)) [13]. This dipole will emit light as well and produce a field, ${E}_{2}$, at the initial dipole. We can connect the two cases in Fig. 1(b) and (c) with the help of the Rayleigh-Carson reciprocity theorem:

With this relation in place we can use the second case to calculate the fields in the first case. For simplicity, we set ${J}_{1}=\left|{J}_{1}\right|{e}_{J1}\delta \left(r-{r}_{1}\right)$ and ${J}_{2}=\left|{J}_{2}\right|{e}_{J2}\delta \left(r-{r}_{2}\right)$ with $\left|{J}_{1}\right|=\left|{J}_{2}\right|$ and both current densities in phase. With this Eq. (1) becomes a simple scalar product

when ${J}_{2}$ is parallel to ${E}_{1}\left({r}_{2}\right)$. We use the convention of vectors being represented by bold letters and scalar amplitudes being represented by light letters. Accordingly, if we know ${E}_{2}\left({r}_{1}\right)$ we also know ${E}_{1}\left({r}_{2}\right)$. In Fig. 1(c), ${J}_{2}$ creates a wave that has an amplitude of ${E}_{2}^{\text{'}}=\frac{1}{4\pi {\u03f5}_{0}}{k}_{0}^{2}\frac{{e}^{jk\left({r}_{1}-{r}_{2}\right)}}{\left|{r}_{1}-{r}_{2}\right|}\frac{\left|{J}_{1}\right|}{\omega}$, where ${k}_{0}=\omega /c$, once it reaches the encircled region [14]. Far away from the source the curvature of the wave fronts can be neglected and the wave can be locally approximated by a plane wave with corresponding amplitude. Thus, to find the electric field at the position of the initial dipole we have to illuminate the encircled region with a plane wave and find the relation between ${E}_{2}\left({r}_{1}\right)$ and ${E}_{2}^{\text{'}}$. The ratio of the amplitudes of the two electric fields is the local field enhancement factor. It can be obtained by a full optical simulation. But nanostructured media have a characteristic length much smaller than the wavelength. Thus the plane wave excitation can be simplified to a quasistatic calculation with a constant bias field along the optical polarization. However the ratio between the outside and the inside field is obtained, we find thatThe intensity radiated in this direction is $I\left({r}_{2}\right)=\frac{{n}_{eff}}{2{Z}_{0}}{\left|{E}_{1}\left({r}_{2}\right)\right|}^{2}$, where ${n}_{eff}$ is the effective refractive index of the composite medium. Note that although we mentioned in the introduction that an effective refractive index is not sufficient to calculate the optical emission, we can use it here to model the wave propagation in the far field.

From the electric field in Eq. (3) we can determine the intensity at position ${r}_{2}$. We can calculate and sequentially integrate the intensity over a closed surface to obtain the emitted power, if we know ${E}_{1}$ at all positions on the surface. This would involve separate backward excitation for each direction. However, we show in the appendix that it is sufficient to illuminate the encircled region in Fig. 1 with three plane waves polarized in the three orthogonal orientations, ${E}_{2,x}^{\text{'}}$, ${E}_{2,y}^{\text{'}}$ and ${E}_{2,z}^{\text{'}}$. Consequently, the power emitted by the dipole amounts to:

*i*as a superscript in ${E}_{2}^{i}$.

To obtain the emission enhancement factor, we need to divide this power by the power the emitter would emit into vacuum which is

Equation (7) is useful in case we want to determine the emission of a dipole with a certain orientation. However, in many cases emitters are unpolarized. As we show in the appendix, the average emission enhancement of an unpolarized emitter still depends on the local field enhancement averaged for three orthogonally polarized incident fields ${E}_{2,x}^{\text{'}}$, ${E}_{2,y}^{\text{'}}$ and ${E}_{2,z}^{\text{'}}$ and looks similar to the polarized case:

Equations (7) and (8) describe a convenient way to calculate the emission enhancement of nanocomposites. In summary, the following two steps are needed: Firstly, one obtains the effective refractive index from Maxwell Garnett theory [10], a numerical simulation or any other suited method. Secondly, one chooses a representative volume around the emitter and calculates the electric field at the position of the emitter due to a static electric bias along the three principal directions. (If the volume features rotational symmetries, calculations for only one or two directions are sufficient.) The fields and the effective index are subsequently plugged into Eq. (7) in case of a polarized emitter or into Eq. (8) in case of an unpolarized emitter. In the following we give a few application examples for this method.

## 3. Diluted nanocomposites

As a first example, we calculate the emission enhancement of a diluted nanocomposite. The term *diluted* means two things in this example. Firstly, the emitting nanoparticles are so far away from each other such that they do not interact by their near fields. Secondly, the volume fraction of the nanoparticles is so small that the effective refractive index of the composite is equal to that of the host material.

To illustrate the importance of the nanoparticle shape on the emission we compare a sphere and a platelet with large aspect ratio. To calculate the emission enhancement factor of a dipole inside a sphere we use Eq. (7) for which we need to know the field enhancement factor. Since we consider a sphere which is much smaller than the wavelength, the field is practically in phase over the whole volume and the problem is equivalent to the case of a sphere biased with a static electric field. In this quasistatic case the field inside the sphere is related to the outside biasing field by [15],

where ${n}_{host}^{}$ is the refractive index of the host material and ${n}_{emitter}$ that of the particle. The field is constant inside the sphere. Inserting Eq. (9) back into Eq. (7) gives the emission enhancement factorDue to the reduced symmetry, the platelet consideration is slightly more complicated. We assume that the platelet has a large aspect ratio such that we can neglect fringe fields at the edges and model it as an infinitely extended slab. The field inside the platelet is then defined by the continuity condition at the platelet-host interfaces. For an electrical field perpendicular to the interfaces we have

and for fields parallel to the interfaces we haveBecause the dipole sees only the fraction of the field that is parallel to it (Eq. (3)), a dipole oriented perpendicular to the interfaces has an emission enhancement factor ofOn the other hand, a dipole parallel to the interfaces has no advantage to an unstructured material of refractive index ${n}_{host}$:

The platelet and the sphere show a strongly different dependence on the host material refractive index. The emission enhancement factor of the spherical nanoparticle scales linearly for a host material index much higher than the particle index, because the field enhancement converges to ${(3/2)}^{2}$. The emission enhancement factor scales linearly for emitters orientated parallel to the interface of the platelet. But the emission of perpendicular emitters scales with the fifth power of the host refractive index. We exemplify this shape dependence for erbium doped CaF_{2} nanocrystals. These crystals emit at around 1.55µm optical wavelength [16] and have a refractive index of 1.43 [17,18]. We investigate how the emission enhancement factor changes if we place the particles inside host materials of varying refractive index. We vary the refractive index from 1 which would be a very porous material like an aerogel [19] to 3.5 which corresponds to silicon.

The resulting curves are shown in Fig. 2. For a spherical particle we see emission enhancement factors of up to 6.7 which is almost 5 times more than what can be achieved in CaF_{2} alone. It is perhaps more practical to compare the emission in the nanostructured medium to the emission in amedium consisting purely out of the bulk host materialinstead of comparing to vacuum. This factor is ${f}_{e,bulk}={\text{f}}_{e}/{n}_{host}$, i.e. in this case ${f}_{e,bulk}={f}_{e}/1.43=4.7$. In stark contrast to the spherical particle the platelet can show an enhancement factor of ${f}_{e}=125({f}_{e,bulk}=87)$. This value is reduced if one considers an unpolarized emitter as the parallel dipoles show only an enhancement scaling with ${n}_{host}$. Nonetheless, with an averaged enhancement of ${f}_{e}=44({f}_{e,bulk}=31)$ the platelet emits 6.5 times stronger than the spherical particle. And the emission enhancement continuously grows proportionally to ${n}_{host}^{5}$.

## 4. Dense nanocomposites

We now turn to dense nanocomposites. In the case of a dense composite the volume fraction of the emitting particles is high enough such that they affect each other’s radiation properties. Further, the emitting material cannot be neglected for the calculation of the effective refractive index of the nanocomposite. To describe emission in these structures we can use the same approach we used for the diluted nanocomposites but for the dense composites it is not justified to neglect fringe fields which are not aligned with the biasing electric field. We can use Eq. (8) and calculate the local fields for every position in the luminescent medium.

As an example for a dense composite, we investigate the two periodic composites shown in Fig. 3. The structure is in both cases a cubic crystal. whose unit cell consists of a cube surrounded by an other material. In Fig. 3(a) the cubes are the emitting material and the surrounding medium is passive and in Fig. 3(b) we consider the reversed case. Again, the emitting material is erbium doped CaF_{2} and the host material’s refractive index is varied from 1 to 3.5. We chose a volume ratio of emitter material to host material of 50:50. To calculate the emission enhancement factor of the system we need to know the effective index of the system and the ratio of the local field to the biasing field. We obtain both via an eigenmode calculation with CST Microwave Studio. A detailed description of the simulation is given in the appendix. In contrast to the previous case, the field enhancement factor now spatially varys inside the emitting material and thus we calculate a volume averaged emission enhancement factor. Because the volume is much smaller than wavelength, the three static solutions yield all the possible optical states in the system. From which follows that our approach is, in this case, completely equivalent to calculating the emission enhancement factor via the local density of states [20].

One might expect that due the same amount of host and emitting material the average emission enhancement factor of the two arrangements in Fig. 3 is the same. As one can see in Fig. 4(a) this is not the case. The composite with emitting material inside the cubes scales linearly and thus resembles the behavior of the isolated sphere in Fig. 2(a). In contrast to that, the material with the emitter outside the cube shows a nonlinear dependence on the host refractive index. The reason for this behavior is similar to the emission enhancement of the platelet and can be seen in Fig. 4(b). If the emitting material is in between the cubes it forms platelet-like structures. The biasing field normal to them causes a field enhancement. The effect is not as strong as for the platelets because the fields are significantly reduced close to the edges of the cubes. Further, in contrast to the diluted composite, the effective refractive index of the dense composite is influenced by the emitting material. For a structure with ${n}_{host}=3.5$, we obtain ${n}_{eff}=2.2$, if the host medium is inside the cubes and ${n}_{eff}=2.5$, if it is outside. The lower ${n}_{eff}$ in the first case is disadvantageous for the emission enhancement (see Eq. (8)), but from Fig. 4(a) we see that the local field enhancement exceeds the effect of the smaller ${n}_{eff}$ leading to a net enhancement.

The effect of the fringe fields can be significantly reduced if the volume fraction of the luminescent material is reduced. For a 10% volume fraction the emission enhancement for a host with refractive index of 3.5 can reach ${f}_{e}=28({f}_{e,bulk}=20)$ if the host is concentrated in the cubes.

## 5. Conclusion

In conclusion, we presented a method to calculate the emission enhancement in nanocomposites and applied this method to diluted and dense composites. While common effective medium models are not applicable and direct Green’s function simulations are cumbersome, simple quasi static simulations can be used to calculate the average emission enhancement. Using the derived equations, we show that the emission can be greatly altered by the shape of the constituents and offers a route towards enhancing the performance of this material class.

Platelet based nanocomposites in high index materials like silicon can provide emission enhancement factors (compared to vacuum emission) in the order of 40, which is of the order of the enhancement factors achieved with high Q dielectric cavities [21,22]. Though higher emission enhancements can be obtained in small gap metal-dielectric resonators [5] and with metal-dielectric metamaterials [23], these structures as well as dielectric cavities require deliberated patterning on the scale of the wavelength or below and cannot easily be scaled up to yield large emitting volumes. The use of metals is also strongly connected to absorption losses. On contrary the presented novel approach can show large emission enhancement in bulk dielectric media even for luminescent material volume fraction up to 10%.

## Appendix 1 Derivation of average emission enhancement factor

Dipole 1 in Fig. 5(a) produces a far field in the direction $\left(\theta ,\phi \right)$ of ${E}_{1}={E}_{1,\theta}+{E}_{1,\phi}$. It produces no fields in the radial direction as only transverse waves propagate in the farfield.

We can calculate the fields in the two polarizations separately by placing test dipoles with respective orientation as indicated in Fig. 5(b) and (c). Like we did in the main text, we can relate the fields by

We assumed that $\left|{J}_{1}\right|=\left|{J}_{2,\theta}\right|=\left|{J}_{2,\phi}\right|$ and that all currents are in phase. We can again calculate the far fields ${E}_{2,\theta}^{\text{'}}$ and ${E}_{2,\phi}^{\text{'}}$ and obtain ${E}_{2}^{\theta}\left({r}_{1}\right)$ and ${E}_{2}^{\phi}\left({r}_{1}\right)$ from a quasistatic consideration. This approach would require a separate calculation for each direction $\left(\theta ,\phi \right)$. We can reduce the number of necessary calculations by expanding the fields at ${r}_{1}$ in three excitations with orthogonal polarizations:

Here, we utilize the superposition principle i.e. that in the quasistatic regime we can decompose the fields ${E}_{2}^{\theta}\left({r}_{1}\right)$ and ${E}_{2}^{\phi}\left({r}_{1}\right)$ into three fields each caused by a biasing field in direction of one of the Cartesian axes. Note that ${E}_{2}^{x}$, ${E}_{2}^{y}$ and ${E}_{2}^{z}$ represent local fields excited by three orthogonal polarizations, and thus they do not have to be orthogonal to each other. We can obtain ${E}_{1,\theta}\left({r}_{2}\right)$ and ${E}_{1,\phi}\left({r}_{2}\right)$ by inserting Eqs. (18) and (19) into (16) and (17).

## Appendix 2 Calculating the emission enhancement factor with an eigenmode solver

We can find the emission enhancement factor with a numerical eigenmode calculation. For this purpose we simulate a unit cell of our periodic structure with boundary conditions as indicated in Fig. 6. Since we are in the quasistatic regime and the mode is practically in phase over the whole unit cell, we set the phase difference $\text{\Delta}\varphi $ to a small number, e.g. 0.1°. Slight changes of this value will proportionally change the mode frequency which is of no relevance since in the quasistatic regime the modal fields are frequency independent. We use a small non zero phase to obtain the local field distribution and the ${n}_{eff}$ from the same simulation.

We use the eigenmode solver of CST Microwave Studio. The solver returns the modes frequency and its fields normalized to carry 1J of electromagnetic energy in the simulation volume. To calculate the emission enhancement factor we need to obtain the effective index and the ratio between local field and biasing field.

We can calculate the effective index from the k-vector:

where $a$ is the volume length, ${c}_{0}$ is the vacuum speed of light and $\omega $ is the modes radial frequency. The biasing field is the average electric field of the wave traveling through the volume. We can calculate it from the power the mode carries:The power can be either obtained by integrating the Poynting flux in or out of the volume or we can use the fact that the mode is normalized to 1J:where ${v}_{g}$ is the group velocity equal to phase velocity in the effective medium. This simplification can be made for periodic structures with periods much smaller than the wavelength if material dispersion can be neglected [24].## Funding

German Federal Ministry of Education and Research (BMBF) (01GR0468); 382-PiGnano project of ERA.Net RUS Plus 2013-2018 initiative under Consortium Agreement with Swiss Federal Laboratories for Materials Science and Technology (Switzerland) and Hamburg University of Technology (Germany); Deutsche Forschungsgemeinschaft (392323616); Hamburg University of Technology.

## Acknowledgments

The authors acknowledge the support of CST, Darmstadt, Germany with their Microwave Studio software.

## References and links

**1. **P. Anger, P. Bharadwaj, and L. Novotny, “Enhancement and Quenching of Single-Molecule Fluorescence,” Phys. Rev. Lett. **96**(11), 113002 (2006). [CrossRef] [PubMed]

**2. **K. L. Tsakmakidis, R. W. Boyd, E. Yablonovitch, and X. Zhang, “Large spontaneous-emission enhancements in metallic nanostructures: towards LEDs faster than lasers,” Opt. Express **24**(16), 17916–17927 (2016). [CrossRef] [PubMed]

**3. **E. M. Purcell, “Spontaneous emission probabilities at radio frequencies,” Phys. Rev. **69**, 681 (1946).

**4. **S. Ek, P. Lunnemann, Y. Chen, E. Semenova, K. Yvind, and J. Mork, “Slow-light-enhanced gain in active photonic crystal waveguides,” Nat. Commun. **5**, 5039 (2014).

**5. **G. M. Akselrod, C. Argyropoulos, T. B. Hoang, C. Ciracì, C. Fang, J. Huang, D. R. Smith, and M. H. Mikkelsen, “Probing the mechanisms of large Purcell enhancement in plasmonic nanoantennas,” Nat. Photonics **8**(11), 835–840 (2014). [CrossRef]

**6. **T. Senden, F. T. Rabouw, and A. Meijerink, “Photonic Effects on the Radiative Decay Rate and Luminescence Quantum Yield of Doped Nanocrystals,” ACS Nano **9**(2), 1801–1808 (2015). [CrossRef] [PubMed]

**7. **K. Dolgaleva, R. W. Boyd, and P. W. Milonni, “Influence of local-field effects on the radiative lifetime of liquid suspensions of Nd:YAG nanoparticles,” J. Opt. Soc. Am. B **24**(3), 516–521 (2007). [CrossRef]

**8. **D. Toptygin, “Effects of the Solvent Refractive Index and Its Dispersion on the Radiative Decay Rate and Extinction Coefficient of a Fluorescent Solute,” J. Fluoresc. **13**(3), 201–219 (2003). [CrossRef]

**9. **L. Rogobete, H. Schniepp, V. Sandoghdar, and C. Henkel, “Spontaneous emission in nanoscopic dielectric particles,” Opt. Lett. **28**(19), 1736–1738 (2003). [CrossRef] [PubMed]

**10. **A. H. Sihvola, *Electromagnetic Mixing Formulas and Applications* (IET, 1999).

**11. **A. E. Krasnok, A. P. Slobozhanyuk, C. R. Simovski, S. A. Tretyakov, A. N. Poddubny, A. E. Miroshnichenko, Y. S. Kivshar, and P. A. Belov, “An antenna model for the Purcell effect,” Sci. Rep. **5**(1), 12956 (2015). [CrossRef] [PubMed]

**12. **A. N. Poddubny, P. A. Belov, and Y. S. Kivshar, “Spontaneous radiation of a finite-size dipole emitter in hyperbolic media,” Phys. Rev. A **84**(2), 23807 (2011). [CrossRef]

**13. **P. Bharadwaj, B. Deutsch, and L. Novotny, “Optical Antennas,” Adv. Opt. Photonics **1**(3), 438–483 (2009). [CrossRef]

**14. **J. D. Jackson, *Classical Electrodynamics*, 3rd ed. (Wiley, 1999).

**15. **D. J. Griffiths, *Introduction to Electrodynamics* (Pearson Education, 2014).

**16. **G. A. Kumar, C. W. Chen, J. Ballato, and R. E. Riman, “Optical Characterization of Infrared Emitting Rare-Earth-Doped Fluoride Nanocrystals and Their Transparent Nanocomposites,” Chem. Mater. **19**(6), 1523–1528 (2007). [CrossRef]

**17. **H. H. Li, “Refractive index of alkaline earth halides and its wavelength and temperature derivatives,” J. Phys. Chem. Ref. Data **9**(1), 161–290 (1980). [CrossRef]

**18. **M. Daimon and A. Masumura, “High-accuracy measurements of the refractive index and its temperature coefficient of calcium fluoride in a wide wavelength range from 138 to 2326 nm,” Appl. Opt. **41**(25), 5275–5281 (2002). [CrossRef] [PubMed]

**19. **T. Bellunato, M. Calvi, C. Matteuzzi, M. Musy, D. L. Perego, and B. Storaci, “Refractive index of silica aerogel: Uniformity and dispersion law,” RICH **2007**(595), 183–186 (2008).

**20. **A. Taflove, A. Oskooi, and S. G. Johnson, *Advances in FDTD Computational Electrodynamics: Photonics and Nanotechnology* (Artech House, 2013).

**21. **M. Pelton, J. Vukovic, G. S. Solomon, A. Scherer, and Y. Yamamoto, “Three-dimensionally confined modes in micropost microcavities: quality factors and Purcell factors,” IEEE J. Quantum Electron. **38**(2), 170–177 (2002). [CrossRef]

**22. **T. J. Kippenberg, A. L. Tchebotareva, J. Kalkman, A. Polman, and K. J. Vahala, “Purcell-Factor-Enhanced Scattering from Si Nanocrystals in an Optical Microcavity,” Phys. Rev. Lett. **103**(2), 027406 (2009). [CrossRef] [PubMed]

**23. **I. Iorsh, A. Poddubny, A. Orlov, P. Belov, and Y. S. Kivshar, “Spontaneous emission enhancement in metal–dielectric metamaterials,” Phys. Lett. A **376**(3), 185–187 (2012). [CrossRef]

**24. **A. A. Krokhin, P. Halevi, and J. Arriaga, “Long-wavelength limit (homogenization) for two-dimensional photonic crystals,” Phys. Rev. B **65**(11), 115208 (2002). [CrossRef]