Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Far-zone analytical expressions for electromagnetic waves produced by a uniform-current loop placed in a uniaxial dielectric-magnetic material

Open Access Open Access

Abstract

We investigate the radiation generated by a current loop situated within an unbounded uniaxial dielectric-magnetic material characterized by a uniform current distribution. We employ dyadic Green functions in the frequency domain to derive closed-form expressions for the far-field radiation. Our analytical findings distinguish between two scenarios: one in which the loop’s axis aligns parallel to the optic axis and another where it is perpendicular to the optic axis. In cases where the loop’s axis parallel’s the optic axis, only magnetically extraordinary wave is emitted as a consequence. However, when the loop’s axis is perpendicular to the optic axis, both electrically and magnetically extraordinary waves are emitted. Our results demonstrate a pronounced dependence of the radiation pattern on the loop’s size, as observed for various loop radii.

© 2024 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Uniaxial materials are anomalous anisotropic materials possessing a distinctive symmetry axis [13]. Electromagnetic waves can propagate through these materials with two mutually orthogonal states of polarization and with different phase velocities [4]. Uniaxial materials naturally exist in the world, e.g., calcite, zircon and rutile [5]. Uniaxial materials can also be assembled artificially, e.g., thin electrical sheets of various isotropic materials arranged as a layered stack [6,7]. Uniaxial materials have attracted peculiar attention due to their application in several aberrant phenomenon, e.g., negative-refraction and sub-wavelength imaging [8]. Although uniaxial material are prevalently employed in fabrication of optical devices, they are also utilized in nano imaging [9].

In past various authors have been studied radiation emitted from a point magnetic-dipole laying in uniaxial materials [4]. Salem and Caloz have analyzed the electromagnetic radiations generated from an arbitrary current carrying loop antenna embedded in the isotropic material [10]. Pozar has studied the radiation and scattering of electromagnetic waves exuding from a microstrip antenna laying in uniaxial-dielectric substrate [11]. Li et al., has computed the electromagnetic fields radiated by a thin circular loop-antenna in both near and far zones by adopting the dyadic Greens functions technique [12]. Ghaffar et al., has studied the focusing of electromagnetic waves into a uniaxial crystal such that the waves are refracted by an inhomogeneous slab [13]. Moreover, Hayat and Faryad have studied the radiation emitted from a Hertzian dipole having finite length placed in a uniaxial dielectric material. They have adopted the dyadic Green function approach in frequency domain to evaluate radiation fields generated by a current loop in near and far fields [14]. Recently, closed-form expressions for electromagnetic waves generated by a current loop in a uniaxial dielectric medium in the far zone were studied by us [15].

To effectively model real radiation sources in uniaxial materials such as quantum dots, wires, and sheets, it is necessary to consider magnetic dipoles. Consequently, we aim to establish the radiation characteristics of an electric current loop, which serves as an effective representation of a magnetic dipole. This paper assumes an arbitrary radius for the loop, making the problem suitable for simulating real extended radiation sources. Additionally, employing a current loop to model the magnetic dipole not only facilitates its direct realization but also extends the field generalization from an isotropic medium to the simplest form of anisotropic medium, namely, a uniaxial dielectric-magnetic medium. In the current study, we establish closed-form expressions for electromagnetic fields emitted from a limited sized loop, where the axis of symmetry of loop is either aligned with or orthogonal to the optic axis of the material.

The article’s overview is as follows. The far-zone approximated dyadic Green functions for the uniaxial dielectric-magnetic material are simulated in Section 2. The analytical results for the fields generated by the current loop when its axis of symmetry is aligned with the optic axis of the uniaxial dielectric-magnetic material, are evaluated in Section 3. The same analytical findings for the instance when loop’s axis (it stands for the axis normal to the disk enclosed by the loop, which is the $z$ axis) of symmetry is orthogonal to the optic axis are presented in Section 4. In Section 5 a comparison for the radiation pattern of a current loop against that of point magnetic dipole has been presented. The final outcomes and conclusion are presented in Sections 6. and 7 respectively.

The time dependence $\exp (i\omega t)$ is implied in this article, where $i=\sqrt {-1}$ and $\omega$ represents angular frequency. Furthermore, $\mu _o$, $\epsilon _o$ and $k_o=\omega \sqrt {\epsilon _o \mu _o}$ symbolize the vacuum permeability, vacuum permittivity and wave number respectively. Vectors are represented by letters with arrow heads and doubly underlined symbols represent the dyadics. The identity dyadic is denoted by $\underline {\underline {I}}$.

2. Dyadic Green functions

The relevance between a vector source and corresponding vector field is explicated in terms of dyadic Green function and can be represented as [4]

$$\vec{E}(\vec{r})=\iiint_{V'} \underline{\underline{G}}^{ee}\vec{(R)}\cdot{\vec{J}}_e{(\vec{r}')}d^3{r'},$$
$${\vec{H}(\vec{r})}=\iiint_{V'} {\underline{\underline{G}}^{me}\vec{(R)}}.{{\vec{J}}_e{(\vec{r}')}} d^3{r'}\,.$$

Here volume covered by current density $\vec {J}_e{(\vec {r}')}$ is denoted by $V'$ with

$$\vec{R} = {\vec{r}}- {\vec{r}'}.$$

The materials whose permittivity and permeability dyadics have a common uniaxiality axis are known as uniaxial dielectric-magnetic materials [4]. The permittivity and permeability dyadics of the uniaxial dielectric-magnetic materials are given as [16]

$$\underline{\underline{\epsilon}}=\epsilon_{o}\underline{\underline{\epsilon}}_{r}=\epsilon_{o}\left[\epsilon^{{\perp}}\underline{\underline{I}}+\left(\epsilon^{{\parallel}}-\epsilon^{{\perp}}\right)\hat{\mathbf{c}}\hat{\mathbf{c}}\right],$$
and
$$\underline{\underline{\mu}}=\mu_{o}\underline{\underline{\mu}}_{r}=\mu_{o}\left[\mu^{{\perp}}\underline{\underline{I}}+\left(\mu^{{\parallel}}-\mu^{{\perp}}\right)\hat{\mathbf{c}}\hat{\mathbf{c}}\right],$$
where $\epsilon ^{\parallel }$, $\epsilon ^{\perp }$, $\mu ^{\parallel }$, and $\mu ^{\perp }$ can be complex values, whereas $\bf {\hat {c}}$ represent the direction of the optic axis of the material.

Since we are concerned with the fields that a current loop exudes in the far-zone, so here the far-zone approximation of the dyadic Green functions are relevant only, and their exact expressions are given in Refs. [4,17]. The approximated dyadic Green functions in the far-zone are given as [4,17]

$$\underline{\underline{G}}^{ee}\vec{(R)} \approx i\omega\mu^{{\perp}}\mu_{o} \Bigg\{ g_{e}({\vec{R}}) \frac{\epsilon^{{\parallel}}}{\epsilon^{{\perp}}}\frac{\left[\vec{R}\times\left(\vec{R}\times{\bf{\hat{c}}}\right)\right]\left[\vec{R}\times\left(\vec{R}\times{\bf{\hat{c}}}\right)\right]}{ R_e^2 \mid{\vec{R}}\times{\bf{\hat{c}}}\mid^2} +\frac{\mu^{{\parallel}}}{\mu^{{\perp}}}g_m(\vec{R})\frac{\left({\vec{R}}\times{\bf{\hat{c}}}\right)\left({\vec{R}}\times{\bf{\hat{c}}}\right)}{\mid{{\vec{R}}}\times{\bf{\hat{c}}}\mid^2}\Bigg\},$$
$${{\underline{\underline{G}}^{me}}{\vec{(R)}}}\approx i k_o n_o \Bigg\{-\frac{\epsilon^{{\parallel}}}{\epsilon^{{\perp}}}g_e({\vec{R}})\frac{\left({\vec{R}}\times{\bf{\hat{c}}}\right)\left[{\vec{R}}\times\left({\vec{R}}\times{\bf{\hat{c}}}\right)\right]}{ {R_e} \mid{\vec{R}}\times{\bf{\hat{c}}}\mid^2} + \frac{\mu^{{\parallel}}}{\mu^{{\perp}}} g_m({\vec{R}}) \frac{\left[{\vec{R}}\times\left({\vec{R}}\times{\bf{\hat{c}}}\right)\right]\left({\vec{R}}\times{\bf{\hat{c}}}\right)}{ {R_m} \mid{\vec{R}}\times{\bf{\hat{c}}}\mid^2} \Bigg\},$$
where
$$n_o=\sqrt{\epsilon^{{\perp}} \mu^{{\perp}}}, \quad k_o=\omega\sqrt{\epsilon_o\mu_o},$$
and the scalar Green functions representing both of the extraordinary waves can be written as [4]
$${g_e(\vec{R})}=\frac{\exp\left(i k_o n_o R_e\right)}{4\pi R_e},\quad {g_m(\vec{R})}=\frac{\exp\left(i k_o n_o R_m\right)}{4\pi R_m},$$
where $R_e$ and $R_m$ are expressed as
$$R_e=\sqrt{\epsilon^{\delta} \left({\vec{R}}\times{\bf{\hat{c}}}\right)^2 + \left({\vec{R}}\cdot{\bf{\hat{c}}}\right)^2} ,$$
and
$$R_m=\sqrt{\mu^{\delta}\left({\vec{R}}\times{\bf{\hat{c}}}\right)^2 + \left({\vec{R}}\cdot{\bf{\hat{c}}}\right)^2},$$
the ratio of the permittivities and the permeabilities of the uniaxial dielectric-magnetic material are given as
$$\epsilon^{\delta}= \frac{\epsilon^{{\parallel}}}{\epsilon^{{\perp}}}, \quad \mu^{\delta}=\frac{\mu^{{\parallel}}}{\mu^{{\perp}}}.$$

3. Optic axis and loop axis are parallel

Let us suppose a current loop having radius $a$ carrying constant current $I_o$ laying in xy plane as illustrated in Fig. 1. In that case the loop axis and the optic axis are directed parallel to the $z$-axis.

 figure: Fig. 1.

Fig. 1. Illustration of a current loop having radius $a$, with a constant current flowing through it and laying in a uniaxial dielectric-magnetic material, such that loop axis and the optic axis are aligned in the same direction.

Download Full Size | PDF

In spherical coordinates electric current density phasor can be expressed as

$${\vec{J}}_e{(\vec{r})} = \frac{I_o}{a} \delta \left(r-a\right)\delta\left(\theta-\frac{\pi}{2}\right) {\hat{\phi}}.$$

The electric field radiated from a current loop in the far-zone can be evaluated using Eqs. (6) and (13) in Eq. (1), so that we get

$$\scalebox{0.9}{$\displaystyle{\vec{E}(\vec{r})}= i \omega \mu_o \mu^{{\perp}} I_o a \int_0^{2\pi} \Bigg\{ \epsilon^{\delta} g_e(\check{\vec{R}})\frac{\left[\check{\vec{R}}\times\left(\check{\vec{R}}\times\hat{\bf{z}}\right)\right]\left[\check{\vec{R}}\times\left(\check{\vec{R}}\times\hat{\bf{z}}\right)\right]}{ {R_e^2} \mid\check{\vec{R}}\times{\bf{\hat{z}}}\mid^2} + \mu^{\delta} g_m(\check{\vec{R}} ) \frac{\left(\check{\vec{R}}\times{\bf{\hat{z}}}\right)\left(\check{\vec{R}}\times{\bf{\hat{z}}}\right)}{\mid\check{\vec{R}}\times{\bf{\hat{z}}}\mid^2}\Bigg\} \cdot {\hat{\phi}^\prime} d\,\phi^\prime,$}$$
where, using Eq. (3) we can write
$$\check{\vec{R}}=(x-x^\prime){\bf{\hat{x}}}+(y-y^\prime){\bf{\hat{y}}}+z\bf{\hat{z}},$$
whereas the ${\hat {\phi }^\prime }$ unit vector in spherical coordinates is given by
$${\hat{\phi}^\prime} ={-}\sin\phi^\prime{\bf{\hat{x}}}+\cos\phi^\prime\bf{\hat{y}}.$$

Since spherical coordinates are related with the Cartesian coordinates as follows

$$x=r\sin{\theta} \cos{\phi}, \quad y=r\sin{\theta} \sin{\phi}, \quad z=r\cos{\theta},$$

As the net electric field and magnetic field are given by:

$$\left\{\begin{array}{l} \vec{E}=\vec{E}_{e}+\vec{E}_{m} \\ \vec{H}=\vec{H}_{e}+\vec{H}_{m} \end{array}\right.$$

Here we are not providing the detailed calculation because these expressions can be obtained by the same method that we already provided in Ref. [15]. Following the same procedure given in Ref. [15] the final expression for electrically extraordinary wave that can be extracted from Eq. (14) is given by

$${\vec{E}}_e(\vec{r}) = \frac{i \mu_o\mu^{{\perp}} I_o c a}{2r^2 n_o \Theta_e^2} \exp(i k_o n_o r_{ez}) \frac{\cos^2{\theta}}{\sin^2{\theta}} \times \textit{J}_1\bigg(\frac{k_o n_o a \epsilon^{\delta} \sin{\theta}}{\Theta_e}\bigg) {\bf{\hat{y}}}.$$
where $J_1 \bigg (\frac {k_o n_o a \epsilon ^{\delta } \sin {\theta }}{\Theta _e}\bigg )$ is the first kind of Bessel function, $r_{ez}=r\Theta _e$ and $\Theta _e = \sqrt {\epsilon ^{\delta }\sin ^2{\theta } + \cos ^2{\theta }}$. As from the start we are retained the terms proportional to $1/r$ in the far-zone, however, the above equation contains $1/r^2$ so under this condition the above equation redue to
$${\vec{E}}_e(\vec{r})\approx 0.$$

It shows that the extraordinary electric field emitted by the current loop has a smaller amplitude, approximated to zero, if the loop axis is aligned with the optic axis of the material.

Similarly the final expression for magnetic part of electric field that can be separated from Eq. (14) is given as

$${\vec{E}}_m(\vec{r}) = \frac{\omega\mu_o \mu^{{\parallel}} I_o a \exp(i k_o n_o r \Theta_m)}{2 r \Theta_m} \textit{J}_1 \left(\frac{k_o n_o a \mu^{\delta} \sin{\theta}}{\Theta_m}\right) {\bf{\hat{y}}}.$$
where $r_{mz} = r\Theta _m$, and $\Theta _m = \sqrt {\mu ^{\delta } \sin ^2{\theta }+cos^2{\theta }}$.

Repeating the same procedure as explained earlier for electric field and also given in Ref. [15] the expression for $\vec {H}_{e}$ is given by

$${\vec{H}_e(\vec{r})}={-} \frac{i I_o a\exp(i k_o n_o r_{ez}) \cos{\theta} }{4\pi r^2 \Theta_e \sin^2{\theta}} \textit{J}_1 \bigg( \frac{\epsilon^{\delta} k_o n_o a\sin{\theta}}{\Theta_e} \bigg) {\bf{\hat{x}}}.$$

Since the final extraordinary magnetic field contains term proportional to $1/r^2$ only, so we can write

$${\vec{H}_e(\vec{r})} \approx 0,$$
which shows that in the far-zone magnetically extraordinary wave is missing from radiation when the loop axis and the optic axis are parallel to the $z$ axis. Moreover, the expression for magnetic field is
$${\vec{H}_m(\vec{r})}=\frac{k_o n_o I_o \mu^{\delta} a \exp(i k_o n_o r \Theta_m)}{2 r \Theta_m^2} {\left(\sin\theta\mathbf{ \hat{z}}-\cos{\color{red}{\theta}}\mathbf{\hat{x}} \right)}\times {\textit{J}_1{\bigg(\frac{k_o n_o \mu^{\delta}a \sin{\theta}}{\Theta_m} \bigg)}}.$$

Hence, Eqs. (21) and (24) reveal that only magnetic part of the extraordinary waves is exuded in the far-zone when the loop axis and the optic axis are aligned in the same direction, i.e., parallel to $z$-axis. Moreover, the emanated fields satisfy the orthogonality condition, similar to the behaviour of the ordinary electromagnetic radiations.

The time averaged-power, emitted from the current loop in the far-zone is given as [18]

$${\frac{d\,P_m}{d\,\Omega}}= \frac{1}{2} {\bf{\hat{r}}} \cdot Re{({\vec{E}}_m\times{\vec{H}_m}^*)}r^2.$$

Using Eqs. (21) and (24) in Eq. (24), and simplifying the expression we can write

$${\frac{d\,P_m}{d\,\Omega}}=\frac{k_o^2 c n_o \mu_o \mu^{{\parallel}} \mu^{\delta} I_o^2 a^2}{8 \Theta_m^3} {\textit{J}}_1^2\bigg( \frac{k_o n_o a \mu^{\delta} \sin{\theta}}{\Theta_m}\bigg) .$$

The final result of magnetic power emanated from a current loop in far-zone do not include $\epsilon ^{\delta }$, i.e., magnetic waves are polarized in a direction orthogonal to the optic axis.

4. Optic axis and the loop axis are mutually perpendicular

Now, let us assume that the optic axis of the material is turned along $x$-axis while axis of the current loop is still in $z$-direction, i.e., both are perpendicular to each other, the schematic diagram of this configuration is depicted in Fig. (2). Here we are going to use the same electric and magnetic current densities described in the Eq. (13), while we choose ${\bf {\hat {c}}}={\bf {\hat {x}}}$ for a uniaxial dielectric-magnetic material.

 figure: Fig. 2.

Fig. 2. Illustration of a current loop having radius $a$, with a constant current flowing through it and laying in a uniaxial dielectric-magnetic material, such that loop axis and the optic axis are mutually perpendicular.

Download Full Size | PDF

Our main goal is to compute the waves emitted from that current loop at a very far away from point P. Since the procedure for computing the expressions for waves and time-averaged was already discussed in Ref. [15] so we will not provide the detailed calculations. We will only present the final result and skip the details. Following the same procedure as followed earlier for Sec. 3 the expressions for $\vec {E}_{e}$, $\vec {E}_{m}$, $\vec {H}_{e}$, and $\vec {H}_{m}$ are given as follows:

$${\vec{E}_e{(\vec{r})}}=\frac{\omega \mu_o\mu^{{\perp}} I_o a \epsilon^{\delta} y \exp(i k_o n_o r \Phi_e)}{2(y^2+z^2) r^2\Phi_e \sin{\theta} \rho_d} \bigg[-(y^2+z^2){\bf{\hat{x}}}+xy{\bf{\hat{y}}}+xz{\bf{\hat{z}}}\bigg]\times \textit{J}_1 \bigg( \frac{k_o n_o a \rho_d \sin{\theta}}{\Phi_e}\bigg),$$
$${\vec{E}_m{(\vec{r})}}=\frac{\omega \mu_o \mu^{{\parallel}} I_o a \exp(i k_o n_o r_{mx}) \cos{\theta} \cos{\phi}}{2\Phi_m (y^2+z^2) \sigma_d} (z{\bf{\hat{y}}}-y{\bf{\hat{z}}}) \times \textit{J}_1 \bigg(\frac{k_o n_o a \sin{\theta} \sigma_d}{\Phi_m}\bigg) ,$$
$${\vec{H}_e(\vec{r})}=\frac{k_o n_o I_o a \epsilon^{\delta} y \exp(i k_o n_o r \Phi_e)}{2 (y^2+z^2) r \rho_d \sin{\theta}} (y{\bf{\hat{z}}}-z{\bf{\hat{y}}}) \textit{J}_1 \bigg(\frac{k_o n_o a \rho_d \sin{\theta}}{\Phi_e}\bigg) ,$$
$${\vec{H}_m(\vec{r})}=\frac{k_o n_o a I_o \mu^{\delta} \exp\left(i k_o n_o r_{mx} \right)}{2 r \Phi_m^2 (y^2+z^2) \sigma_d} \cos{\theta} \cos{\phi} \bigg[-(y^2+z^2){\bf{\hat{x}}} +xy{\bf{\hat{y}}}+xz{\bf{\hat{z}}}\bigg] \textit{J}_1 \left(\frac{k_o n_o a \sigma_d \sin{\theta}}{\Phi_m}\right) ,$$
where different factors used in above equations are given as follows:
$$r_{(e,m)x}= r\Phi_{(e,m)} ,$$
$$\Phi_{(e,m)}= \sqrt{(\epsilon^{\delta}, \mu^{\delta})(\sin^2{\theta}\sin^2{\phi}+\cos^2{\theta})+\sin^2{\theta}\cos^2{\phi}},$$
$$(\rho_d,\sigma_d) = \sqrt{\cos^2{\phi}+({\epsilon^{\delta}}^2,{\mu^{\delta}}^2)\sin^2{\phi}}.$$

Furthermore, it is observed that in the far-zone, the radiated extraordinary wave consists of both electrically extraordinary wave and magnetically extraordinary wave parts for the instance when axis of the loop is orthogonal to the optic axis of the material. We can add the powers of both parts of the extraordinary wave to get the total power radiated per unit time per unit solid angle, given as

$$\frac{d\,P}{d\,\Omega}=\frac{d\,P_e}{d\,\Omega} +\frac{d\,P_m}{d\,\Omega},$$
where for electric part we have
$$\frac{d\,P_e}{d\,\Omega}=\frac{1}{2} {\bf{\hat{r}}} \cdot Re{({\vec{E}}_e\times{\vec{H}_e}^*)}r^2,$$
and for magnetic part
$$\frac{d\,P_m}{d\,\Omega}=\frac{1}{2} {\bf{\hat{r}}} \cdot Re{({\vec{E}}_m\times{\vec{H}_m}^*)}r^2.$$

In order to calculate time-averaged power of the extraordinary waves emitted by the current loop in the far-zone using Eqs. (26) and (27) in Eq. (32) and switching the variables into spherical coordinates we can write

$$\frac{d\,P_e}{d\,\Omega}=\frac{ck_o^2 n_o \mu_o \mu^{{\perp}} I_o^2 a^2 {\epsilon^{\delta}}^2 }{8 \Phi_e \rho_d^2} \left( \frac{ \sin^2{\phi}}{\sin^2{\theta}\sin^2{\phi}+\cos^2{\theta}} \right) \times\textit{J}_1^2\left(\frac{k_o n_o a \rho_d \sin{\theta}}{\Phi_e}\right).$$

The time-averaged power of the magnetically extraordinary wave emitted by the current loop in the far-zone can be computed by using Eqs. (27) and (28) within Eq. (33) we can obtain

$$\frac{d\,P_m}{d\,\Omega}=\frac{c k_o^2 n_o \mu_o \mu^{{\parallel}} \mu^{\delta} I_o^2 a^2 } {8 \Phi_m^3\sigma_d^2 } \left(\frac{\cos^2{\phi}\cos^2{\theta}}{\sin^2{\theta}\sin^2{\phi}+\cos^2{\theta}} \right)\times \textit{J}_1^2\left(\frac{k_o n_o a \sigma_d \sin{\theta}}{ \Phi_m}\right).$$

Now, the total power emitted by the current loop per unit time per unit solid angle can be computed by using Eqs. (34) and (35) in Eq. (31) we obtain

$$\begin{aligned} \frac{d\,P}{d\,\Omega} =& \frac{k_o^2 c n_o \mu_o I_o^2 a^2}{8 (\sin^2{\theta}\sin^2{\phi}+\cos^2{\theta})} \Bigg[ \frac{\mu^{{\perp}} {\epsilon^{\delta}}^2 \sin^2{\phi}}{\Phi_e \rho_d^2} \textit{J}_1^2\left(\frac{k_o n_o a \rho_d \sin{\theta}}{\Phi_e}\right)\\ &+ \frac{\mu^{{\parallel}} \mu^{\delta} \cos^2{\theta}\cos^2{\phi}}{\Phi_m^3 \sigma_d^2} \textit{J}_1^2 \left(\frac{k_o n_o a \sigma_d \sin{\theta}}{\Phi_m}\right) \Bigg] . \end{aligned}$$

On substituting $\epsilon ^{\delta }=1$ (which gives $\Phi _e=1$ and $\rho _d=1$) and $\mu ^{\delta }=1$ (which means that $\Phi _m=1$ and $\sigma _d=1$ ) in the Eq. (36) and by substituting $\Theta _m =1$ and $\mu ^{\delta }=1$ into Eq. (25) we observe that both equations are equivalent, displaying same result as for the power radiated by the current loop in an isotropic material [10].

5. Validation of work

A point magnetic dipole can be effectively emulated by a small electric current loop when the loop’s radius $a$ is very short indicating that $|k_o n_o a|<<1$. Under this approximation, it becomes possible to replicate the outcomes described in Secs. 3 and 4 for a point magnetic dipole. These achieved results can then be compared to those presented in Chapter 3, Section 3.5.3 of Ref. [4]. Additionally, by selecting values such as $\epsilon ^{\delta }=1$, $\mu ^{\delta }=1$, $\Theta _m=1$, $\Phi _e=1$ and $\Phi _m=1$, a point magnetic dipole placed within an isotropic material can be obtained. Notably, these findings align with those detailed in Section 2.1.4 of Chapter 2 in Ref. [4].

6. Numerical results and discussion

The main outcomes of this article were the analytical expressions of the electromagnetic waves emitted by a current loop in the far-zone, represented in Eqs. (25) and (34)–(36). The far-zone magnetically extraordinary wave emanated by a current loop when the loop axis is aligned in the direction of the optic axis of the material is illustrated in Fig. (3) for $a=0.1\lambda _o$, and for $a=0.3\lambda _o$. To demonstrate their applications, we give illustrative numerical results for radiations emitted by different current loops having various radii, while they are placed in a uniaxial dielectric-magnetic material having parameters as follows $\mu ^{\parallel }=1.5$, $\mu ^{\perp }=1.2$, $\epsilon ^{\parallel }=1.8$, $\epsilon ^{\perp }=2.5$, $I_o=0.1 A$ and $\lambda _o=0.584$ $\mu m$ [19] for the current loops with radii $a=0.1\lambda _o$ and $a=0.3\lambda _o$.

 figure: Fig. 3.

Fig. 3. Far-zone radiation patterns for the time averaged power of the magnetically extraordinary wave emitted by a current loop that is computed in Eq. (25), when the loop axis is aligned with optic axis of the material. The closed-patterns are sketched for parameters $\epsilon ^{\parallel }=1.8$, $\epsilon ^{\perp }=2.5$, $\mu ^{\perp }=1.2$, $\mu ^{\parallel }=1.5$ and $I_o=0.1 A$ with wavelength of $\lambda _o=0.584$ $\mu m$ in angle ranges $\theta \in [0,\pi ]$ and $\phi \in [\pi /2 , 3\pi /2]$ for current loop having radius $a=0.1\lambda _o$ (left) and $a=0.3 \lambda _o$ (right).

Download Full Size | PDF

For more prominent visualization we sketched Eq. (25) for ranges $\theta \in [0,\pi ]$ and $\phi \in [\pi /2 , 3\pi /2]$. The plots display that wave patterns are analogous with the plot for radiations emitted by a point magnetic-dipole, when laying in an isotropic material.

Moreover, with the increase in the loop size, i.e., $a=0.3 \lambda _o$ emitted radiations are suppressed along the loop and spread along axis of the current loop. Also no radiations are emitted along the optic axis.

When the optic axis of the material is at right angles to the axis of the current loop, the extraordinary electric part emission pattern are illustrated in Fig. 4. Since the pattern is symmetric about $yz$-plan so we are sketching Eq. (34) for ranges $\theta \in [0,\pi ]$ and $\phi \in [\pi /2 , 3\pi /2]$ for two different radii of the current loop, i.e., for $a=0.1\lambda _o$ and $a=0.3 \lambda _o$.

 figure: Fig. 4.

Fig. 4. Far-zone radiation patterns for the time averaged power of the electrically extraordinary wave emitted by a current loop that is computed in Eq. (34), when the loop axis is orthogonal to the optic axis of the material. The closed-patterns are sketched for parameters $\epsilon ^{\parallel }=1.8$, $\epsilon ^{\perp }=2.5$, $\mu ^{\perp }=1.2$ and $I_o=0.1 A$ with wavelength of $\lambda _o=0.584$ $\mu m$ in angle ranges $\theta \in [0,\pi ]$ and $\phi \in [\pi /2 , 3\pi /2]$ for current loop having radius $a=0.1\lambda _o$ (left) and $a=0.3 \lambda _o$ (right).

Download Full Size | PDF

These emission patterns reveal that for small loop size $a=0.1 \lambda _o$ radiations are emitted in a direction that is perpendicular to both loop axis and optic axis as plotted in Fig. 4(a). Whereas Fig. 4(b) depicts emission pattern of a large sized loop $a=0.3 \lambda _o$, illustrating that emitted radiations are entirely suppressed along loop axis and the plane perpendicular to it, with the maximum suppression in a direction normal to both the loop axis and the optic axis.

The emission pattern of magnetically extraordinary wave, evaluated in Eq. (35), exuded by the current loop when the loop axis and optic axis of the material are orthogonal to each other are illustrated in Fig. 5 for $a=0.1\lambda _o$ and $a=0.3 \lambda _o$.

 figure: Fig. 5.

Fig. 5. Identical to Fig. 4, but the pattern is for the extraordinary wave’s magnetic part as evaluated in Eq. (35).

Download Full Size | PDF

The radiations are pointed perpendicular to the loop’s axis as well as the optic axis of the material, when the loop size is minimal. However, when the loop’s size increases, the radiations are wholly suppressed along the loop’s axis and the plane which is orthogonal to it, with the sharpest suppression being in a direction that is orthogonal to both of optic axis and the loop’s axis.

7. Conclusions

Analytical results were computed to comprehend the patterns of radiations exuded in the far-zone from a uniform current loop laying in a uniaxial dielectric-magnetic material. The corresponding numerical findings were reported for a selected uniaxial dielectric-magnetic material. Similar far-field results of emitted fields and power were evaluated for a point magnetic-dipole. The loop axis was assumed to be identical or orthogonal to optic axis of the material. When loop axis and optic axis were aligned in same direction only magnetically extraordinary wave were emanated in far-zone. Both electrically extraordinary and magnetically extraordinary waves were exuded from the current loop in far-zone when its axis is laying normal to the direction of optic axis. Although both electrically extraordinary wave and magnetically extraordinary waves were squashed in the direction of optic axis of material. In all of these instances it is obvious that the enlargement of loop size appreciably increased the directionality of emitted radiations patterns. For limiting case, i.e., for loop of very small radius exuded radiation patterns were similar with that of a point magnetic-dipole laying in same uniaxial dielectric-magnetic material. It is evidently noticed that radiation patterns of a current loop have an accurate resemblance with that of a point magnetic-dipole, solely if the loop size is quite tiny.

Funding

Ministry of Education – Kingdom of Saudi Arabia (IFKSUOR3–611–2); Zhejiang Normal University (YS304023905).

Acknowledgments

The authors extend their appreciation to the Deputyship for Research and Innovation, "Ministry of Education" in Saudi Arabia for funding this research (IFKSUOR3-611-2). The authors also acknowledge, M. Shah the postdoctoral fellowship supported by Zhejiang Normal University, China who have equally contributed to this work by checking the different sections of the manuscript.

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. C. D. Gribble, Optical Mineralogy: Principles And Practice (CRC Press, 1993).

2. O. Wiener, “Die theorie des mischkorpers fur das feld der stationaren stromung,” Abhandlungen der Sachsischen Gesellschaft der Akademischen Wissenschaften in Mathematik und Physik 32, 507–604 (1912).

3. A. Lakhtakia, “Constraints on effective constitutive parameters of certain bianisotropic laminated composite materials,” Electromagnetics 29(6), 508–514 (2009). [CrossRef]  

4. M. Faryad and A. Lakhtakia, Infinite-space dyadic Green functions in electromagnetism, vol. 2 (Morgan Claypool Press, 2018).

5. G. Jellison, F. Modine, and L. Boatner, “Measurement of the optical functions of uniaxial materials by two-modulator generalized ellipsometry: rutile (tio 2),” Opt. Lett. 22(23), 1808–1810 (1997). [CrossRef]  

6. M. Born and E. Wolf, Principles of optics: electromagnetic theory of propagation, interference and diffraction of light (Elsevier, 2013).

7. D. R. Smith, W. J. Padilla, and D. Vier, “Composite medium with simultaneously negative permeability and permittivity,” Phys. Rev. Lett. 84(18), 4184–4187 (2000). [CrossRef]  

8. R. A. Depine and A. Lakhtakia, “Diffraction by a grating made of a uniaxial dielectric–magnetic medium exhibiting negative refraction,” New J. Phys. 7, 158 (2005). [CrossRef]  

9. K. G. Balmain, A. Luttgen, and P. C. Kremer, “Resonance cone formation, reflection, refraction, and focusing in a planar anisotropic metamaterial,” Antennas Wirel. Propag. Lett. 1, 146–149 (2002). [CrossRef]  

10. M. A. Salem and C. Caloz, “Electromagnetic fields radiated by a circular loop with arbitrary current,” IEEE Trans. Antennas Propag. 63(1), 442–446 (2015). [CrossRef]  

11. D. Pozar, “Radiation and scattering from a microstrip patch on a uniaxial substrate,” IEEE Trans. Antennas Propag. 35(6), 613–621 (1987). [CrossRef]  

12. L.-W. Li, M.-S. Leong, P.-S. Kooi, et al., “Exact solutions of electromagnetic fields in both near and far zones radiated by thin circular-loop antennas: A general representation,” IEEE Trans. Antennas Propag. 45(12), 1741–1748 (1997). [CrossRef]  

13. A. Ghaffar, S. Latif, and Q. Naqvi, “Study of the focusing of a field refracted by an inhomogeneous slab into a uniaxial crystal by using maslov’s method,” J. Mod. Opt. 55(21), 3513–3528 (2008). [CrossRef]  

14. A. Hayat and M. Faryad, “On the radiation from a hertzian dipole of a finite length in the uniaxial dielectric medium,” OSA Continuum 2(4), 1411–1429 (2019). [CrossRef]  

15. A. Hayat and M. Faryad, “Closed-form expressions for electromagnetic waves generated by a current loop in a uniaxial dielectric medium in the far zone,” J. Opt. Soc. Am. B 36(8), F9–F17 (2019). [CrossRef]  

16. A. Lakhtakia, V. V. Varadan, and V. K. Varadan, “Time-harmonic and time-dependent dyadic green’s functions for some uniaxial gyro-electromagnetic media,” Appl. Opt. 28(6), 1049–1052 (1989). [CrossRef]  

17. W. S. Weiglhofer, “Dyadic green’s functions for general uniaxial media,” in IEE Proceedings H-Microwaves, Antennas and Propagation, vol. 137 (IET, 1990), pp. 5–10.

18. J. D. Jackson, Classical Electrodynamics (John Wiley & Sons, 1999).

19. A. Yariv and P. Yeh, Optical waves in crystals, vol. 5 (Wiley, 1984).

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (5)

Fig. 1.
Fig. 1. Illustration of a current loop having radius $a$, with a constant current flowing through it and laying in a uniaxial dielectric-magnetic material, such that loop axis and the optic axis are aligned in the same direction.
Fig. 2.
Fig. 2. Illustration of a current loop having radius $a$, with a constant current flowing through it and laying in a uniaxial dielectric-magnetic material, such that loop axis and the optic axis are mutually perpendicular.
Fig. 3.
Fig. 3. Far-zone radiation patterns for the time averaged power of the magnetically extraordinary wave emitted by a current loop that is computed in Eq. (25), when the loop axis is aligned with optic axis of the material. The closed-patterns are sketched for parameters $\epsilon ^{\parallel }=1.8$, $\epsilon ^{\perp }=2.5$, $\mu ^{\perp }=1.2$, $\mu ^{\parallel }=1.5$ and $I_o=0.1 A$ with wavelength of $\lambda _o=0.584$ $\mu m$ in angle ranges $\theta \in [0,\pi ]$ and $\phi \in [\pi /2 , 3\pi /2]$ for current loop having radius $a=0.1\lambda _o$ (left) and $a=0.3 \lambda _o$ (right).
Fig. 4.
Fig. 4. Far-zone radiation patterns for the time averaged power of the electrically extraordinary wave emitted by a current loop that is computed in Eq. (34), when the loop axis is orthogonal to the optic axis of the material. The closed-patterns are sketched for parameters $\epsilon ^{\parallel }=1.8$, $\epsilon ^{\perp }=2.5$, $\mu ^{\perp }=1.2$ and $I_o=0.1 A$ with wavelength of $\lambda _o=0.584$ $\mu m$ in angle ranges $\theta \in [0,\pi ]$ and $\phi \in [\pi /2 , 3\pi /2]$ for current loop having radius $a=0.1\lambda _o$ (left) and $a=0.3 \lambda _o$ (right).
Fig. 5.
Fig. 5. Identical to Fig. 4, but the pattern is for the extraordinary wave’s magnetic part as evaluated in Eq. (35).

Equations (39)

Equations on this page are rendered with MathJax. Learn more.

$$\vec{E}(\vec{r})=\iiint_{V'} \underline{\underline{G}}^{ee}\vec{(R)}\cdot{\vec{J}}_e{(\vec{r}')}d^3{r'},$$
$${\vec{H}(\vec{r})}=\iiint_{V'} {\underline{\underline{G}}^{me}\vec{(R)}}.{{\vec{J}}_e{(\vec{r}')}} d^3{r'}\,.$$
$$\vec{R} = {\vec{r}}- {\vec{r}'}.$$
$$\underline{\underline{\epsilon}}=\epsilon_{o}\underline{\underline{\epsilon}}_{r}=\epsilon_{o}\left[\epsilon^{{\perp}}\underline{\underline{I}}+\left(\epsilon^{{\parallel}}-\epsilon^{{\perp}}\right)\hat{\mathbf{c}}\hat{\mathbf{c}}\right],$$
$$\underline{\underline{\mu}}=\mu_{o}\underline{\underline{\mu}}_{r}=\mu_{o}\left[\mu^{{\perp}}\underline{\underline{I}}+\left(\mu^{{\parallel}}-\mu^{{\perp}}\right)\hat{\mathbf{c}}\hat{\mathbf{c}}\right],$$
$$\underline{\underline{G}}^{ee}\vec{(R)} \approx i\omega\mu^{{\perp}}\mu_{o} \Bigg\{ g_{e}({\vec{R}}) \frac{\epsilon^{{\parallel}}}{\epsilon^{{\perp}}}\frac{\left[\vec{R}\times\left(\vec{R}\times{\bf{\hat{c}}}\right)\right]\left[\vec{R}\times\left(\vec{R}\times{\bf{\hat{c}}}\right)\right]}{ R_e^2 \mid{\vec{R}}\times{\bf{\hat{c}}}\mid^2} +\frac{\mu^{{\parallel}}}{\mu^{{\perp}}}g_m(\vec{R})\frac{\left({\vec{R}}\times{\bf{\hat{c}}}\right)\left({\vec{R}}\times{\bf{\hat{c}}}\right)}{\mid{{\vec{R}}}\times{\bf{\hat{c}}}\mid^2}\Bigg\},$$
$${{\underline{\underline{G}}^{me}}{\vec{(R)}}}\approx i k_o n_o \Bigg\{-\frac{\epsilon^{{\parallel}}}{\epsilon^{{\perp}}}g_e({\vec{R}})\frac{\left({\vec{R}}\times{\bf{\hat{c}}}\right)\left[{\vec{R}}\times\left({\vec{R}}\times{\bf{\hat{c}}}\right)\right]}{ {R_e} \mid{\vec{R}}\times{\bf{\hat{c}}}\mid^2} + \frac{\mu^{{\parallel}}}{\mu^{{\perp}}} g_m({\vec{R}}) \frac{\left[{\vec{R}}\times\left({\vec{R}}\times{\bf{\hat{c}}}\right)\right]\left({\vec{R}}\times{\bf{\hat{c}}}\right)}{ {R_m} \mid{\vec{R}}\times{\bf{\hat{c}}}\mid^2} \Bigg\},$$
$$n_o=\sqrt{\epsilon^{{\perp}} \mu^{{\perp}}}, \quad k_o=\omega\sqrt{\epsilon_o\mu_o},$$
$${g_e(\vec{R})}=\frac{\exp\left(i k_o n_o R_e\right)}{4\pi R_e},\quad {g_m(\vec{R})}=\frac{\exp\left(i k_o n_o R_m\right)}{4\pi R_m},$$
$$R_e=\sqrt{\epsilon^{\delta} \left({\vec{R}}\times{\bf{\hat{c}}}\right)^2 + \left({\vec{R}}\cdot{\bf{\hat{c}}}\right)^2} ,$$
$$R_m=\sqrt{\mu^{\delta}\left({\vec{R}}\times{\bf{\hat{c}}}\right)^2 + \left({\vec{R}}\cdot{\bf{\hat{c}}}\right)^2},$$
$$\epsilon^{\delta}= \frac{\epsilon^{{\parallel}}}{\epsilon^{{\perp}}}, \quad \mu^{\delta}=\frac{\mu^{{\parallel}}}{\mu^{{\perp}}}.$$
$${\vec{J}}_e{(\vec{r})} = \frac{I_o}{a} \delta \left(r-a\right)\delta\left(\theta-\frac{\pi}{2}\right) {\hat{\phi}}.$$
$$\scalebox{0.9}{$\displaystyle{\vec{E}(\vec{r})}= i \omega \mu_o \mu^{{\perp}} I_o a \int_0^{2\pi} \Bigg\{ \epsilon^{\delta} g_e(\check{\vec{R}})\frac{\left[\check{\vec{R}}\times\left(\check{\vec{R}}\times\hat{\bf{z}}\right)\right]\left[\check{\vec{R}}\times\left(\check{\vec{R}}\times\hat{\bf{z}}\right)\right]}{ {R_e^2} \mid\check{\vec{R}}\times{\bf{\hat{z}}}\mid^2} + \mu^{\delta} g_m(\check{\vec{R}} ) \frac{\left(\check{\vec{R}}\times{\bf{\hat{z}}}\right)\left(\check{\vec{R}}\times{\bf{\hat{z}}}\right)}{\mid\check{\vec{R}}\times{\bf{\hat{z}}}\mid^2}\Bigg\} \cdot {\hat{\phi}^\prime} d\,\phi^\prime,$}$$
$$\check{\vec{R}}=(x-x^\prime){\bf{\hat{x}}}+(y-y^\prime){\bf{\hat{y}}}+z\bf{\hat{z}},$$
$${\hat{\phi}^\prime} ={-}\sin\phi^\prime{\bf{\hat{x}}}+\cos\phi^\prime\bf{\hat{y}}.$$
$$x=r\sin{\theta} \cos{\phi}, \quad y=r\sin{\theta} \sin{\phi}, \quad z=r\cos{\theta},$$
$$\left\{\begin{array}{l} \vec{E}=\vec{E}_{e}+\vec{E}_{m} \\ \vec{H}=\vec{H}_{e}+\vec{H}_{m} \end{array}\right.$$
$${\vec{E}}_e(\vec{r}) = \frac{i \mu_o\mu^{{\perp}} I_o c a}{2r^2 n_o \Theta_e^2} \exp(i k_o n_o r_{ez}) \frac{\cos^2{\theta}}{\sin^2{\theta}} \times \textit{J}_1\bigg(\frac{k_o n_o a \epsilon^{\delta} \sin{\theta}}{\Theta_e}\bigg) {\bf{\hat{y}}}.$$
$${\vec{E}}_e(\vec{r})\approx 0.$$
$${\vec{E}}_m(\vec{r}) = \frac{\omega\mu_o \mu^{{\parallel}} I_o a \exp(i k_o n_o r \Theta_m)}{2 r \Theta_m} \textit{J}_1 \left(\frac{k_o n_o a \mu^{\delta} \sin{\theta}}{\Theta_m}\right) {\bf{\hat{y}}}.$$
$${\vec{H}_e(\vec{r})}={-} \frac{i I_o a\exp(i k_o n_o r_{ez}) \cos{\theta} }{4\pi r^2 \Theta_e \sin^2{\theta}} \textit{J}_1 \bigg( \frac{\epsilon^{\delta} k_o n_o a\sin{\theta}}{\Theta_e} \bigg) {\bf{\hat{x}}}.$$
$${\vec{H}_e(\vec{r})} \approx 0,$$
$${\vec{H}_m(\vec{r})}=\frac{k_o n_o I_o \mu^{\delta} a \exp(i k_o n_o r \Theta_m)}{2 r \Theta_m^2} {\left(\sin\theta\mathbf{ \hat{z}}-\cos{\color{red}{\theta}}\mathbf{\hat{x}} \right)}\times {\textit{J}_1{\bigg(\frac{k_o n_o \mu^{\delta}a \sin{\theta}}{\Theta_m} \bigg)}}.$$
$${\frac{d\,P_m}{d\,\Omega}}= \frac{1}{2} {\bf{\hat{r}}} \cdot Re{({\vec{E}}_m\times{\vec{H}_m}^*)}r^2.$$
$${\frac{d\,P_m}{d\,\Omega}}=\frac{k_o^2 c n_o \mu_o \mu^{{\parallel}} \mu^{\delta} I_o^2 a^2}{8 \Theta_m^3} {\textit{J}}_1^2\bigg( \frac{k_o n_o a \mu^{\delta} \sin{\theta}}{\Theta_m}\bigg) .$$
$${\vec{E}_e{(\vec{r})}}=\frac{\omega \mu_o\mu^{{\perp}} I_o a \epsilon^{\delta} y \exp(i k_o n_o r \Phi_e)}{2(y^2+z^2) r^2\Phi_e \sin{\theta} \rho_d} \bigg[-(y^2+z^2){\bf{\hat{x}}}+xy{\bf{\hat{y}}}+xz{\bf{\hat{z}}}\bigg]\times \textit{J}_1 \bigg( \frac{k_o n_o a \rho_d \sin{\theta}}{\Phi_e}\bigg),$$
$${\vec{E}_m{(\vec{r})}}=\frac{\omega \mu_o \mu^{{\parallel}} I_o a \exp(i k_o n_o r_{mx}) \cos{\theta} \cos{\phi}}{2\Phi_m (y^2+z^2) \sigma_d} (z{\bf{\hat{y}}}-y{\bf{\hat{z}}}) \times \textit{J}_1 \bigg(\frac{k_o n_o a \sin{\theta} \sigma_d}{\Phi_m}\bigg) ,$$
$${\vec{H}_e(\vec{r})}=\frac{k_o n_o I_o a \epsilon^{\delta} y \exp(i k_o n_o r \Phi_e)}{2 (y^2+z^2) r \rho_d \sin{\theta}} (y{\bf{\hat{z}}}-z{\bf{\hat{y}}}) \textit{J}_1 \bigg(\frac{k_o n_o a \rho_d \sin{\theta}}{\Phi_e}\bigg) ,$$
$${\vec{H}_m(\vec{r})}=\frac{k_o n_o a I_o \mu^{\delta} \exp\left(i k_o n_o r_{mx} \right)}{2 r \Phi_m^2 (y^2+z^2) \sigma_d} \cos{\theta} \cos{\phi} \bigg[-(y^2+z^2){\bf{\hat{x}}} +xy{\bf{\hat{y}}}+xz{\bf{\hat{z}}}\bigg] \textit{J}_1 \left(\frac{k_o n_o a \sigma_d \sin{\theta}}{\Phi_m}\right) ,$$
$$r_{(e,m)x}= r\Phi_{(e,m)} ,$$
$$\Phi_{(e,m)}= \sqrt{(\epsilon^{\delta}, \mu^{\delta})(\sin^2{\theta}\sin^2{\phi}+\cos^2{\theta})+\sin^2{\theta}\cos^2{\phi}},$$
$$(\rho_d,\sigma_d) = \sqrt{\cos^2{\phi}+({\epsilon^{\delta}}^2,{\mu^{\delta}}^2)\sin^2{\phi}}.$$
$$\frac{d\,P}{d\,\Omega}=\frac{d\,P_e}{d\,\Omega} +\frac{d\,P_m}{d\,\Omega},$$
$$\frac{d\,P_e}{d\,\Omega}=\frac{1}{2} {\bf{\hat{r}}} \cdot Re{({\vec{E}}_e\times{\vec{H}_e}^*)}r^2,$$
$$\frac{d\,P_m}{d\,\Omega}=\frac{1}{2} {\bf{\hat{r}}} \cdot Re{({\vec{E}}_m\times{\vec{H}_m}^*)}r^2.$$
$$\frac{d\,P_e}{d\,\Omega}=\frac{ck_o^2 n_o \mu_o \mu^{{\perp}} I_o^2 a^2 {\epsilon^{\delta}}^2 }{8 \Phi_e \rho_d^2} \left( \frac{ \sin^2{\phi}}{\sin^2{\theta}\sin^2{\phi}+\cos^2{\theta}} \right) \times\textit{J}_1^2\left(\frac{k_o n_o a \rho_d \sin{\theta}}{\Phi_e}\right).$$
$$\frac{d\,P_m}{d\,\Omega}=\frac{c k_o^2 n_o \mu_o \mu^{{\parallel}} \mu^{\delta} I_o^2 a^2 } {8 \Phi_m^3\sigma_d^2 } \left(\frac{\cos^2{\phi}\cos^2{\theta}}{\sin^2{\theta}\sin^2{\phi}+\cos^2{\theta}} \right)\times \textit{J}_1^2\left(\frac{k_o n_o a \sigma_d \sin{\theta}}{ \Phi_m}\right).$$
$$\begin{aligned} \frac{d\,P}{d\,\Omega} =& \frac{k_o^2 c n_o \mu_o I_o^2 a^2}{8 (\sin^2{\theta}\sin^2{\phi}+\cos^2{\theta})} \Bigg[ \frac{\mu^{{\perp}} {\epsilon^{\delta}}^2 \sin^2{\phi}}{\Phi_e \rho_d^2} \textit{J}_1^2\left(\frac{k_o n_o a \rho_d \sin{\theta}}{\Phi_e}\right)\\ &+ \frac{\mu^{{\parallel}} \mu^{\delta} \cos^2{\theta}\cos^2{\phi}}{\Phi_m^3 \sigma_d^2} \textit{J}_1^2 \left(\frac{k_o n_o a \sigma_d \sin{\theta}}{\Phi_m}\right) \Bigg] . \end{aligned}$$
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.