Figure 1(a) shows the fundamental structure (structure A) for theoretical study with a spherical nano-cavity of R in radius embedded in an infinitely large material of n2 in refractive index. The refractive index of the medium filling the nano-cavity is n1. Here, the origin of the coordinate system is set at the center of the spherical cavity, at which an x- or z-oriented source dipole is located. In Fig. 1(b) we illustrate the structure (structure B) for studying the FRET behavior when both energy donor and acceptor are located inside a spherical nano-cavity at the coordinate origin and (0, 0. d) respectively. Structures C-E illustrated in Figs. 1(c)-1(e) are used for investigating the nanoscale-cavity effect on SP coupling. Figure 1(f) illustrates the reference structure (structure R) in which a radiating dipole is situated in a homogeneous medium of n1 in refractive index.
To study the nanoscale-cavity effect, we first derive the electromagnetic field distribution produced by an x-oriented dipole with the intrinsic dipole moments \({\vec {p}_0}={p_0}\hat {x}\) in a spherical nano-cavity structure (structure A) based on the classical electromagnetic theory. It is noted that the electromagnetic derivation does not include the Purcell effect, which describes the emission behavior change of a source dipole caused by the scattered field from its surrounding structure. In other words, the scattered field of the source dipole can change its own dipole moment. In this situation, the derived results are still valid if the source dipole moment \({\vec {p}_0}={p_0}\hat {x}\) is replaced by a modified dipole moment \({\vec {p}_m}={p_m}\hat {x}\). The modified dipole moment can be obtained by assuming that the source dipole acts as a two-level oscillation system and solving its optical Bloch equations under the influence of the scattered field [23, 30]. Generally, the modified dipole moment \({\vec {p}_m}={\vec {p}_0}+\Delta \vec {p}\) can be derived from the steady-state solution of the optical Bloch equations to give [30]
$$\Delta \vec {p}= - \frac{{i{T_2}\left( {\rho _{{11}}^{{(0)}} - \rho _{{00}}^{{(0)}}} \right)}}{\hbar }\frac{{{{\left| {\frac{{{{\vec {p}}_0} \cdot \vec {E}}}{2}} \right|}^2}}}{{{{\left| {\vec {E}} \right|}^2}}}\frac{{\vec {E}}}{{1+\frac{{{T_1}{T_2}}}{\hbar }{{\left| {\frac{{{{\vec {p}}_0} \cdot \vec {E}}}{2}} \right|}^2}}}$$
1
.
Here, \(\vec {E}\) is the scattered electric field at the dipole position, \(\rho _{{11}}^{{(0)}}\) and \(\rho _{{00}}^{{(0)}}\) are the density matrix elements at zero field of the upper and lower states in the two-level system, respectively. Also, T1 and T2 are the population relaxation time and the dipole moment dephasing time, respectively. We can set T2 = 2T1 in the absence of collision broadening [31]. According to the Wigner-Weisskopf theory of spontaneous emission, T1 is equal to the reciprocal of Einstein’s A coefficient as [32]
$${T_1}=\frac{{6{\varepsilon _0}h{c^3}}}{{{\omega ^3}{n_b}{{\left| {{{\vec {p}}_0}} \right|}^2}}}$$
2
Here, h is the Planck constant, ω is the angular frequency of the electromagnetic field, nb is the background refractive index, ε0 and c are the dielectric permittivity and light speed, respectively, in vacuum. From the last factor of Eq. (1) we can see the saturation effect, indicating that \(\Delta \vec {p}\) is small when the scattered electric field is very strong. The radiation of the modified source dipole produces a new scattered field, which leads to a new solution of the optical Bloch equations and hence a new modified dipole moment. Such a procedure is iterated until a steady state is reached for giving the final modified dipole moment. The numerical implementation of the Purcell effect is called the feedback process [23]. In numerical computations, we choose \(\rho _{{11}}^{{(0)}} - \rho _{{00}}^{{(0)}}=1\) and \(\left| {{{\vec {p}}_0}} \right|\)= 2 x 10− 10/ω (Ω-m). The choice for the \(\left| {{{\vec {p}}_0}} \right|\) value is actually immaterial because it will be normalized later in the numerical results shown below. Depending on the surrounding structure, the orientation of the source dipole can also be changed after the feedback process is applied. However, in a cavity structure of circular symmetry with respect to the dipole orientation, such as structure A illustrated in Fig. 1(a) the feedback process does not change the dipole orientation. Because of the strong scattering by the nano-cavity structure, the Purcell effect or feedback process is important in the current study. In particular, with the nanoscale dimension, the scattering of near field for producing the Purcell effect can provide us with novel emission behaviors for useful applications.
In a spherical nano-cavity structure, the radiated fields, including near and far fields, from a dipole at an arbitrary position inside the cavity can be analytically derived. However, to simplify the demonstration of formulations, we only show the derivations in a simple case, in which the source dipole is placed at the center of the spherical cavity, as illustrated in Fig. 1(a). In terms of the spherical coordinates (r, θ, φ) the total fields at \(\overrightarrow{r}\) inside the spherical cavity generated by an x-dipole with the modified dipole moment \({\vec {p}_m}={p_m}\hat {x}\) are derived to give
$$\begin{gathered} E_{r}^{{in}}=\frac{{i{k^2}{p_m}}}{{2\pi {\varepsilon _1}{\varepsilon _0}r}}\left\{ {h_{1}^{{(1)}}(kr) - \alpha {j_1}\left( {kr} \right)} \right\}\sin \theta \cos \phi \\ =\frac{{i{k^3}{p_m}}}{{2\pi {\varepsilon _1}{\varepsilon _0}}}\left\{ {\left( {2\alpha - 1} \right)\frac{{\cos (kr)}}{{{{(kr)}^2}}} - (2\alpha +1)\frac{{\sin (kr)}}{{{{(kr{\text{)}}}^3}}} - i\left( {\frac{{\sin (kr)}}{{{{(kr)}^2}}} - \frac{{\cos (kr)}}{{{{(kr{\text{)}}}^3}}}} \right)} \right\}\sin \theta \cos \phi \\ \end{gathered}$$
3
$$\begin{gathered} E_{\theta }^{{in}}=\frac{{i{k^2}{p_m}}}{{4\pi {\varepsilon _1}{\varepsilon _0}r}}\left\{ {\frac{{d\left[ {rh_{1}^{{(1)}}(kr)} \right]}}{{dr}} - \alpha \frac{{d\left[ {r{j_1}(kr)} \right]}}{{dr}}} \right\}\cos \theta \cos \phi \\ = - \frac{{i{k^3}{p_m}}}{{4\pi {\varepsilon _1}{\varepsilon _0}}}\left\{ {(\alpha - 1)\left( {\frac{{\sin (kr)}}{{kr}}+\frac{{\cos (kr)}}{{{{(kr)}^2}}} - \frac{{\sin (kr)}}{{{{(kr)}^3}}}} \right)+i\left( {\frac{{\cos (kr)}}{{kr}} - \frac{{\sin (kr)}}{{{{(kr)}^2}}} - \frac{{\cos (kr)}}{{{{(kr)}^3}}}} \right)} \right\}\cos \theta \cos \phi \\ \end{gathered}$$
4
$$\begin{gathered} E_{\phi }^{{in}}= - \frac{{i{k^2}{p_m}}}{{4\pi {\varepsilon _1}{\varepsilon _0}r}}\left\{ {\frac{{d\left[ {rh_{1}^{{(1)}}(kr)} \right]}}{{dr}} - \alpha \frac{{d\left[ {r{j_1}(kr)} \right]}}{{dr}}} \right\}\sin \phi \\ =\frac{{i{k^3}{p_m}}}{{4\pi {\varepsilon _1}{\varepsilon _0}}}\left\{ {(\alpha - 1)\left( {\frac{{\sin (kr)}}{{kr}}+\frac{{\cos (kr)}}{{{{(kr)}^2}}} - \frac{{\sin (kr)}}{{{{(kr)}^3}}}} \right)+i\left( {\frac{{\cos (kr)}}{{kr}} - \frac{{\sin (kr)}}{{{{(kr)}^2}}} - \frac{{\cos (kr)}}{{{{(kr)}^3}}}} \right)} \right\}\sin \phi \\ \end{gathered}$$
5
with
$$\alpha =\frac{{{\zeta ^2}h_{1}^{{(1)}}(\kappa R){{\left\{ {\frac{{d\left[ {rh_{1}^{{(1)}}(kr)} \right]}}{{dr}}} \right\}}_{r=R}} - h_{1}^{{(1)}}(kR){{\left\{ {\frac{{d\left[ {rh_{1}^{{(1)}}(\kappa r)} \right]}}{{dr}}} \right\}}_{r=R}}}}{{{\zeta ^2}h_{1}^{{(1)}}(\kappa R){{\left\{ {\frac{{d\left[ {r{j_1}(kr)} \right]}}{{dr}}} \right\}}_{r=R}} - {j_1}(kR){{\left\{ {\frac{{d\left[ {rh_{1}^{{(1)}}(\kappa r)} \right]}}{{dr}}} \right\}}_{r=R}}}}$$
6
.
Here, \(h_{1}^{{(1)}}(x)\) is a spherical Hankel function of the first kind with order 1 and \({j_1}(x)\) is a spherical Bessel function of the first kind with order 1. Also, \({\varepsilon _1}=n_{1}^{2}\) and \({\varepsilon _2}=n_{2}^{2}\) (k = 2πn1/λ and κ = 2πn2/λ) are the dielectric constants (propagation constants) inside and outside the cavity, respectively. Meanwhile, ζ = n2/n1. In addition, λ is the wavelength in vacuum. The electric fields outside the cavity are derived to give
$$E_{r}^{{out}}=\frac{{i{k^2}\beta {p_m}}}{{2\pi {\varepsilon _1}{\varepsilon _0}r}}h_{1}^{{(1)}}\left( {\kappa r} \right)\sin \theta \cos \phi = - 2\frac{{{k^3}\zeta {p_m}}}{{4\pi {\varepsilon _1}{\varepsilon _0}}}\beta \left( {\frac{i}{{{{\left( {\kappa r} \right)}^2}}} - \frac{1}{{{{\left( {\kappa r} \right)}^3}}}} \right){e^{i\kappa r}}\sin \theta \cos \phi$$
7
$$E_{\theta }^{{out}}=\frac{{i{k^2}\beta {p_m}}}{{4\pi {\varepsilon _1}{\varepsilon _0}r}}\frac{{dh_{1}^{{(1)}}(\kappa r)}}{{dr}}\cos \theta \cos \phi =\frac{{{k^3}\zeta {p_m}}}{{4\pi {\varepsilon _1}{\varepsilon _0}}}\beta \left( {\frac{1}{{\kappa r}}+\frac{i}{{{{\left( {\kappa r} \right)}^2}}} - \frac{1}{{{{\left( {\kappa r} \right)}^3}}}} \right){e^{i\kappa r}}\cos \theta \cos \phi$$
8
$$E_{\phi }^{{out}}= - \frac{{i{k^2}\beta {p_m}}}{{4\pi {\varepsilon _1}{\varepsilon _0}r}}\frac{{d\left[ {rh_{1}^{{(1)}}(\kappa r)} \right]}}{{dr}}\sin \phi = - \frac{{{k^3}\zeta {p_m}}}{{4\pi {\varepsilon _1}{\varepsilon _0}}}\beta \left( {\frac{1}{{\kappa r}}+\frac{i}{{{{\left( {\kappa r} \right)}^2}}} - \frac{1}{{{{\left( {\kappa r} \right)}^3}}}} \right){e^{i\kappa r}}\sin \phi$$
9
with
$$\begin{gathered} \beta =\frac{{h_{1}^{{(1)}}(kR){{\left\{ {\frac{{d[r{j_1}(kr)]}}{{dr}}} \right\}}_{r=R}} - {j_1}(kR){{\left\{ {\frac{{d[rh_{1}^{{(1)}}(kr)]}}{{dr}}} \right\}}_{r=R}}}}{{{\zeta ^2}h_{1}^{{(1)}}(\kappa R){{\left\{ {\frac{{d[r{j_1}(kr)]}}{{dr}}} \right\}}_{r=R}} - {j_1}(kR){{\left\{ {\frac{{d[rh_{1}^{{(1)}}(\kappa r)]}}{{dr}}} \right\}}_{r=R}}}} \\ = - i{e^{ - i\kappa R}}{\left\{ \begin{gathered} \frac{1}{\zeta }\left[ { - {\zeta ^2}\left( {\sin (kR)+\frac{{\cos (kR)}}{{kR}} - \frac{{\sin (kR)}}{{{{(kR)}^2}}}} \right) - \left( { - \frac{{\cos (kR)}}{{kR}}+\frac{{\sin (kR)}}{{{{(kR)}^2}}}} \right)} \right] \hfill \\ - i\left[ {{\zeta ^2}\left( {\sin (kR)+\frac{{\cos (kR)}}{{kR}} - \frac{{\sin (kR)}}{{{{(kR)}^2}}}} \right)\frac{{kR}}{{{{(\kappa R)}^2}}}+\left( {\cos (kR) - \frac{{\sin (kR)}}{{kR}}} \right)\left( {1 - \frac{1}{{{{(\kappa R)}^2}}}} \right)} \right] \hfill \\ \end{gathered} \right\}^{ - 1}} \\ \end{gathered}$$
10
.
For obtaining the expressions after the second equality signs in Eqs. (3)-(5) and (7)-(10) we have used the expansions of the two special functions as
$$h_{1}^{{(1)}}(kr)= - \frac{{{e^{ikr}}}}{{kr}}\left( {1+\frac{i}{{kr}}} \right)$$
11
$${j_1}(kr)= - \frac{{\cos (kr)}}{{kr}}+\frac{{\sin (kr)}}{{{{(kr)}^2}}}$$
12
.
In the equations above, the nanoscale-cavity effect plays its roles through α and β besides the modified dipole moment, pm. If the cavity does not exist, i.e., \({n_1}={n_2}\), \({\varepsilon _1}={\varepsilon _2}\), and k = κ, α = 0 and β = 1. In this situation, Eqs. (3)-(5) and (7)-(9) are reduced to the well-known near-field expressions produced by a radiating dipole in a homogeneous medium [33].
The FRET efficiency is determined by the near-field energy of the donor absorbed by the acceptor. This absorbed donor energy is proportional to the total field intensity at the position of the acceptor produced by the donor if the acceptor absorption coefficient is not affected by the donor intensity. In structure B illustrated in Fig. 1(b) in which the donor and acceptor are located at the coordinates of (0, 0, 0) and (0, 0, d) respectively, the enhancement factor of the FRET efficiency caused by the nanoscale-cavity effect is proportional to the normalized field intensity at the position of the acceptor. With an x-dipole, the normalized field intensity, In(x) is given by
$${I_{n(x)}}=\frac{{{{\left| {E_{\theta }^{{in}}} \right|}^2}}}{{{{\left| {E_{\theta }^{0}} \right|}^2}}}=\frac{{{{\left| {\frac{{i{k^2}{p_m}}}{{4\pi {\varepsilon _1}{\varepsilon _0}}}\frac{1}{r}\left\{ { - \alpha \frac{{dr{j_1}(kr)}}{{dr}}+\frac{{drh_{1}^{{(1)}}(kr)}}{{dr}}} \right\}} \right|}^2}}}{{{{\left| {i\frac{{{k^2}{p_0}}}{{4\pi {\varepsilon _1}{\varepsilon _0}}}\frac{1}{r}\left( {\frac{{drh_{1}^{{(1)}}(kr)}}{{dr}}} \right)} \right|}^2}}}={\left| {\frac{{{p_m}}}{{{p_0}}}} \right|^2}\left\{ {1+{{\left| {\frac{{\frac{{dr{j_1}(kr)}}{{dr}}}}{{\frac{{drh_{1}^{{(1)}}(kr)}}{{dr}}}}} \right|}^2}{{\left| \alpha \right|}^2} - \operatorname{Re} \left( {\frac{{\frac{{dr{j_1}(kr)}}{{dr}}}}{{\frac{{drh_{1}^{{(1)}}(kr)}}{{dr}}}}\alpha } \right)} \right\}$$
13
.
Here, \(E_{\theta }^{0}\) is the θ-component of the electric field with the dipole source in the homogeneous medium of n1 in refractive index [structure R in Fig. 1(f)]. It can be obtained from \(E_{\theta }^{{in}}\) by setting α = 0 and pm = p0. From Eq. (6) when kR < < 1 and κR < < 1, we can obtain
$$\alpha \sim \frac{{3i\left( {{\zeta ^2} - 1} \right)}}{{2{\zeta ^2}+1}}\frac{1}{{{{\left( {kR} \right)}^3}}}\left\{ {1+\frac{3}{5}\left( {\frac{{4{\zeta ^2}+1}}{{2{\zeta ^2}+1}}} \right){{\left( {kR} \right)}^2}+O\left[ {{{\left( {kR} \right)}^3}} \right]} \right\}$$
14
.
We can see that α diverges when kR→0. Under the condition of kR < < 1 and κR < < 1, the normalized field intensity in Eq. (13) is reduced to
$${I_{n(x)}}\sim {\left| {\frac{{{p_m}}}{{{p_0}}}} \right|^2}\left\{ {{{\left[ {1 - \frac{{2\left( {{\zeta ^2} - 1} \right)}}{{2{\zeta ^2}+1}}{\eta ^3}} \right]}^2} - \frac{6}{5}\left( {\frac{{{\zeta ^2} - 1}}{{2{\zeta ^2}+1}}} \right)\left[ {1 - \frac{{2\left( {{\zeta ^2} - 1} \right)}}{{2{\zeta ^2}+1}}{\eta ^3}} \right]\left[ {{\eta ^2}+2\left( {\frac{{4{\zeta ^2}+1}}{{2{\zeta ^2}+1}}} \right)} \right]{\eta ^3}{{\left( {kR} \right)}^2}} \right\}$$
15
.
Here, η = r/R with r < R. When kR→0, Eq. (15) is reduced to
$${I_{n(x)}}\sim {\left| {\frac{{{p_m}}}{{{p_0}}}} \right|^2}{\left[ {1 - \frac{{2\left( {{\zeta ^2} - 1} \right)}}{{2{\zeta ^2}+1}}{\eta ^3}} \right]^2}$$
16
.
In Eq. (15) we can see that the normalized field intensity depends on η. When η << 1, i.e., near the source dipole, the normalized field intensity is equal to \({\left| {{p_m}/{p_0}} \right|^2}\). In other words, only the Purcell effect can contribute to the intensity enhancement. As a matter of fact, when η << 1, pm→p0, i.e., the Purcell effect is negligibly weak. This is so because near the source dipole the intensity is dominated by the direct radiation field. The scattered field makes an insignificant contribution to the total field near the source dipole.
When a z-oriented source dipole at the cavity center is considered, we can also derive the electric field distribution and hence the normalized intensity. Under the condition of kR < < 1, we can obtain the normalized intensity, In(z) as
$${I_{n(z)}}\sim {\left| {\frac{{{p_m}}}{{{p_0}}}} \right|^2}\left\{ {{{\left[ {1+\left( {\frac{{{\zeta ^2} - 1}}{{2{\zeta ^2}+1}}} \right){\eta ^3}} \right]}^2} - \frac{6}{5}\left( {\frac{{{\zeta ^2} - 1}}{{2{\zeta ^2}+1}}} \right)\left[ {1+\left( {\frac{{{\zeta ^2} - 1}}{{2{\zeta ^2}+1}}} \right){\eta ^3}} \right]\left[ {{\eta ^2} - \left( {\frac{{4{\zeta ^2}+1}}{{2{\zeta ^2}+1}}} \right)} \right]{\eta ^3}{{\left( {kR} \right)}^2}} \right\}$$
17
.
When kR→0, Eq. (17) is reduced to
$${I_{n(z)}}\sim {\left| {\frac{{{p_m}}}{{{p_0}}}} \right|^2}{\left[ {1+\left( {\frac{{{\zeta ^2} - 1}}{{2{\zeta ^2}+1}}} \right){\eta ^3}} \right]^2}$$
18
.
The limiting result in Eq. (18) for z-dipole is different from that in Eq. (16) for x-dipole. In Eq. (16) besides the factor of dipole moment ratio, the normalized intensity is proportional to a factor smaller than unity. However, in Eq. (18) this proportion factor is always larger than unity. This comparison result implies that the nanoscale-cavity effect on FRET efficiency depends on the alignment of the donor and acceptor positions with respect to the donor dipole orientation. When the donor and acceptor positions are aligned with the donor dipole orientation, FRET efficiency can more likely be enhanced through the nanoscale-cavity effect.
In the far field limit, i.e., κr > > 1, the electrical fields outside the cavity shown in Eqs. (7)-(9) can be simplified to \(E_{r}^{{out}}\sim 0\),
$$E_{\theta }^{{out}}\sim \frac{{{k^3}\zeta {p_m}}}{{4\pi {\varepsilon _1}{\varepsilon _0}}}\beta \frac{{{e^{i\kappa r}}\cos \theta \cos \phi }}{{\kappa r}}$$
19
$$E_{\phi }^{{out}}\sim - \frac{{{k^3}\zeta {p_m}}}{{4\pi {\varepsilon _1}{\varepsilon _0}}}\beta \frac{{{e^{i\kappa r}}\sin \phi }}{{\kappa r}}$$
20
Equations (19) and (20) are used for evaluating the total radiated power, P, of a dipole inside the cavity to give
$$\begin{gathered} P=\frac{1}{2}\operatorname{Re} \left\{ {\int_{0}^{{2\pi }} {\int_{0}^{\pi } {\left( {{{{\mathbf{\vec {E}}}}^{out}} \times {{{\mathbf{\vec {H}}}}^{out*}}} \right) \cdot {\mathbf{\hat {r}}}{r^2}\sin \theta d\theta d\phi } } } \right\}=\frac{{{n_2}}}{{2{Z_0}}}\int_{0}^{{2\pi }} {\int_{0}^{\pi } {\left( {{{\left| {E_{\theta }^{{out}}} \right|}^2}+{{\left| {E_{\phi }^{{out}}} \right|}^2}} \right){r^2}\sin \theta d\theta d\phi } } \\ =\frac{{k_{0}^{4}{{\left| {{p_m}} \right|}^2}{n_2}}}{{12\pi \varepsilon _{0}^{2}{Z_0}}}{\left| \beta \right|^2} \\ \end{gathered}$$
21
Here, Z0 is the intrinsic impedance in vacuum. Without the Purcell effect or feedback process, the total radiated power, P0, of a dipole \({\vec {p}_0}={p_0}\hat {x}\) in a homogeneous medium of n1 in refractive index is given by
$${P_0}=\frac{{k_{0}^{4}p_{0}^{2}{n_1}}}{{12\pi \varepsilon _{0}^{2}{Z_0}}}$$
22
.
Therefore, the normalized radiated power Pn is given by
$${P_n}=\frac{P}{{{P_0}}}=\zeta {\left| \beta \right|^2}{\left| {\frac{{{p_m}}}{{{p_0}}}} \right|^2}$$
23
.
Here, we can see that the radiated power enhancement caused by the nanoscale-cavity effect is determined by two factors, including \({\left| {{p_m}/{p_0}} \right|^2}\) and \(\zeta {\left| \beta \right|^2}\). The former (latter) is caused by the Purcell effect (classical scattering). In the limiting case of kR > > 1,
$${\left| \beta \right|^2}\sim \frac{1}{{\left( {{\zeta ^2} - 1} \right){{\sin }^2}(kR)+1}}$$
24
.
Therefore, in the far-field limit, when R is large, the normalized radiated power shows a periodical function of photon energy with hc/(2n1R) in period. By using n1 = 1.577 and R = 250 (200) nm, the variation period is 1.576 (1.970) eV. In another limiting case of kR < < 1,
$${\left| \beta \right|^2}\sim \frac{{9{\zeta ^4}}}{{{{\left( {2{\zeta ^2}+1} \right)}^2}}}\left\{ {1 - \frac{1}{5}\left( {\frac{{10{\zeta ^4} - 9{\zeta ^2}+1}}{{2{\zeta ^2}+1}}} \right){{(kR)}^2}} \right\}+O({(kR)^4})$$
25
.
When kR→0, we can obtain
\({\left| \beta \right|^2} \to {\left( {\frac{{3{\zeta ^2}}}{{2{\zeta ^2}+1}}} \right)^2}\) , \({P_n} \to \frac{{9{\zeta ^5}}}{{{{\left( {2{\zeta ^2}+1} \right)}^2}}}{\left| {\frac{{{p_m}}}{{{p_0}}}} \right|^2}\). (26)
It is noted that when kR→0, \({\left| {{p_m}/{p_0}} \right|^2}\)→1 because insignificant Purcell effect can be produced under this condition. Therefore, by assuming that n1 = 1.577 and n2 = 2.399, \({\left| \beta \right|^2}\)→1.523, and Pn →2.317 when kR→0. It is also noted that when the source dipole is placed at the center of a spherical cavity, the normalized radiated power is independent of the dipole orientation.