# The Terahertz Controlled Duplex Isolator: Physical Grounds and Numerical Experiment

Authors: Konstantin Vytovtov; Said Zouhdi; Rostislav Dubrovka; Volodymyr Hnatushenko
Journal: International Journal of Microwave Science and Technology
 (2016)

---

## Abstract

Electromagnetic properties of an anisotropic stratified slab with an arbitrary orientation of the anisotropy axis under an oblique incidence of a plane harmonic wave are studied. The dependence of the eigenwave wavenumbers and the reflection coefficient on an anisotropy axis orientation and frequency is investigated. For the first time, the expression for the translation matrix is obtained in the compact analytical form. The controlled two-way dual-frequency (duplex) isolator based on the above described slab is presented for the first time. It is based on the properties of the anisotropic structure described here but not on the Faraday effect.


---

## Body

## 1. Introduction

The terahertz domain actively developed in the last decades [1–5]. Various devices have been created and numerous electromagnetic structures have been used in this frequency range. One of the most interesting classes of terahertz devices is the class of nonreciprocal devices [1–10]. Usually, the Faraday rotation effect itself is exploited naturally to achieve desired nonreciprocal performance [1]. The electrically tunable Faraday effect in aHgTe thin film has been reported in [2]. Faraday rotation has also been studied for pump pulse of terahertz radiation [3]. In [4], the heterostructures based on the piezoelectric and semiconductor layers have been proposed for radio frequency applications. Other examples of nonreciprocal properties of nonlinear devices in terahertz range have been presented in [5] where a Faraday isolator operating on the coupled microwave resonators circuit has been proposed. A theoretical study of the reflection of infrared radiation from antiferromagnets and M-type barium hexagonal ferrite using an attenuated total reflection geometry has been presented in [6]. Nonreciprocal devices using attenuated total reflection and thin film magnetic layered structures have been described in [7]. Thus, we can clearly see that this topic is quite relevant in modern microwave science.In this paper, we consider nonreciprocal properties of an anisotropic structure based on dependence of the reflection coefficient on an incident wave direction and an anisotropy axis orientation, but not on the Faraday effect. The physical grounds and numerical simulations of a controlled bidirectional dual-frequency isolator not based on the Faraday effect are presented. For this, the analytical investigation and numerical calculations of a stratified anisotropic structure are carried out. At first, the elements of the translation matrix are written in the compact analytical form (4) for the general case of an anisotropic medium for the first time. It is important that these elements are expressed in elementary algebraic functions. The obtained form is very useful in subsequent numerical experiments.In this work, it also is found that the reflection coefficients are different for different incidence anglesα and -
α. Thus, a structure shows nonreciprocal properties.These properties allow us to build nonreciprocal elements in the terahertz domain. The physical principles of construction of a two-way double-frequency (duplex) isolator are described here. These principles are based on dependence of reflection coefficient on an incidence wave direction and an orientation of an anisotropy axis. It should be noted that the Faraday effect is also present in these structures and also defines nonreciprocal properties of anisotropic structures. Today, a wide variety of isolators is described in contemporary scientific literature [7–9]. A feature of the presented device is the fact that it is a two-way dual-frequency one. Switching an isolator (changing a direction and a frequency) is implemented by two solenoids. Its tuning is done mechanically by changing an orientation of external magnetic field. Note that here a very complicated geometry of the problem is used to create the isolator. But only this geometry allows us to obtain the pronounced minimums of reflectance for a given incidence angle, as well as a wide range of total reflection.


## 2. Description of the Isolator Scheme

A presentation of the offered duplex isolator is shown in Figure1(a). This device contains two transceiver antennas, two solenoids generating a magnetizing field, and a stratified anisotropic slab (a nonreciprocal element).Figure 1
(a) Appearance of the controlled double-frequency bidirectional isolator; (b) geometry of an anisotropy layer.
(a)
(b)An orientation of the axes within a homogeneous layer of a stratified anisotropic slab is presented in Figure1(b). Here, z
′ is an anisotropy axis, z is a normal to a structure, k
1 and k
2 are wave vectors of forward refracted waves within an anisotropic layer, β is an inclination angle, ψ is an angle between an incidence plane and a plane containing an anisotropy axis, and α is an incidence angle. The choice of such a complicated geometry is due to the need to obtain an asymmetric (nonreciprocal) structure, the need to obtain wide range of total reflection, and the need to control the structure. Of course nonreciprocal properties can be observed in the simplest case of orientation of an anisotropy axis along interfaces [10, 11], but in that case total reflection regions are not sufficiently wide and a structure has low selective properties in comparison with the given case. Additionally, here mechanical tuning of an isolator is assumed. And using such a complicate geometry gives us the additional freedom to control the device. Moreover, properties of a structure having such a geometry are poorly understood. Therefore, the results of these investigations have both practical and theoretical significance.The layers of a stratified structure are described by a scalar permittivity and a dyadic permeability in the gyrotropic form:(1)
μ
=
μ
x
x
j
μ
x
y
0
-
j
μ
x
y
μ
x
x
0
0
0
μ
z
zwith(2)
μ
x
x
=
1
+
ω
M
ω
0
+
j
/
T
ω
0
+
j
/
T
2
-
ω
2
,
μ
x
y
=
ω
M
ω
ω
0
+
j
/
T
2
-
ω
2
.Here, ω
M
=
4
π
γ
M, ω
0
=
γ
H
0, γ is the gyromagnetic relation, H
0
=
0.1
T is an external magnetic field, and T is a relaxation time.In our numerical investigation, we take into account losses and frequency dispersion of the ferrite material and the elements in (1) are written in the Lanfau-Lifshits form [12]. The dependence of the permeability elements on a frequency for the first and second layers of a period is given in Figure 2. Here, μ
x
x
′ and μ
x
y
′ are the real parts of the elements, and μ
x
x
′
′ and μ
x
y
′
′ are the imaginary parts of the elements. The resonance frequencies are 1.1 GHz and 1.8 GHz for the first and second layers accordingly.Figure 2
The dependence of permeability elements on frequency for the first (a) and second (b) layers of a period.
(a)
(b)


## 3. Method of Calculation

The well-known translation matrix method and the reflection matrix method [13, 14] are used in this study of a stratified anisotropic slab. Both of these methods are accurate and analytical and these can be used in any frequency range for both lossless and lossy media without restrictions. According to the definition, a translation matrix relates tangential field components on both boundaries of a plane-parallel structure:(3)
E
x
z
E
y
z
H
x
z
H
y
z
=
M
z
E
x
0
E
y
0
H
x
0
H
y
0
.Here, M
(
z
) is a 2
×
2 matrix for an isotropic medium or a 4
×
4 translation matrix for an anisotropic medium.The method is described in the literature in detail for both isotropic [13] and anisotropic [14] media. The matrices have been found for a normal orientation of an anisotropy axis, for a tangential orientation of an anisotropy axis, for arbitrary orientation of an anisotropy axis, and for a normal incidence of a plane harmonic wave. However, the expressions for the most general case of an anisotropic medium with an arbitrary orientation of the anisotropy axes under an oblique incidence of a wave have not been obtained in a compact analytical form. In this paper, the translation matrix elements of a homogeneous layer are presented for the first time as (4)
M
m
l
=
∑
i
=
1
4
-
1
m
+
i
γ
m
i
det
⁡
G
l
i
det
⁡
G
exp
⁡
-
j
k
i
z
.Here, G
l
i is the minor of the element γ
m
i of the 4 × 4-matrix G. Matrix G is defined by the expression (5)
E
x
E
y
H
x
H
y
=
γ
11
γ
12
γ
13
γ
14
γ
21
γ
22
γ
23
γ
24
γ
31
γ
32
γ
33
γ
34
γ
41
γ
42
γ
43
γ
44
A
1
e
x
p
⁡
-
j
k
1
z
A
2
e
x
p
⁡
-
j
k
2
z
A
3
e
x
p
⁡
-
j
k
3
z
A
4
e
x
p
⁡
-
j
k
4
z
=
G
A
1
e
x
p
⁡
-
j
k
1
z
A
2
e
x
p
⁡
-
j
k
2
z
A
3
e
x
p
⁡
-
j
k
3
z
A
4
e
x
p
⁡
-
j
k
4
z
,where A
i is a constant. Here alsok
iis a wavenumber defined by well-known fourth-order dispersion relation [12–16]:(6)
k
4
+
a
3
k
3
+
a
2
k
2
+
a
1
k
=
0
.The coefficients a
i in (6) must be obtained directly from Maxwell equations.The expression of matrixG in (5) can be obtained directly from Maxwell equations. Actually, expression (5) is a general solution of Maxwell’s equations for this case [12–14] (see Appendix). The achieved form of the translation matrix (4) essentially simplifies the solution of the programming problem and improves accuracy of numerical calculations. Moreover, no analytical expression of the translation matrix element has been presented in scientific literature for this complicate case.The reflection matrix relates the tangential components of incident and reflected waves [15]:(7)
E
x
r
e
f
l
E
y
r
e
f
l
=
R
E
x
i
n
c
E
y
i
n
c
.Here E
x
i
n
c and E
y
i
n
c are the components of an incident field, E
x
r
e
f
l and E
y
r
e
f
l are the components of a reflected field, and R is a 2 × 2 reflection matrix.


## 4. Nonreciprocal Properties of Anisotropic Structure

The main element of the isolator is a stratified anisotropic slab (Figure1) possessing nonreciprocal properties. In this section, we consider the main physical features of a slab, which ensure its isolator properties. An orientation of the axes for the studied case is shown in Figure 1(b). Angle ψ is given in Figure 1(b). The presented complicated geometry of the problem gives us the interesting features of electromagnetic field in the structure. Such properties cannot appear in a case of more simple geometry. In particular, study of wavenumbers in a homogeneous anisotropic layer suggests that the wavenumbers can be real, imaginary, or complex even in a lossless medium. In the considered case, imaginary parts of wavenumbers indicate existence of a complex wave [16] but not in presence of losses in an anisotropic structure. Indeed a resulting structure of electromagnetic field in a layer is determined by a superposition of all four eigenwaves only, but not each eigenwave separately. In addition, these numbers are not the same for the incidence angles α and −α (Figure 3). The numerical calculations analytical investigations give us the following result (Figure 4): k
1
α
=
-
k
4
(
-
α
), k
2
α
=
-
k
3
(
-
α
), k
3
α
=
-
k
2
(
-
α
), and k
4
α
=
-
k
1
(
-
α
). Thus, we can see nonreciprocal properties of the medium in the considered case.Figure 3
The dependence of the wavenumbers on an incidence angle (β
=
60
°, ψ
=
30
°, f
=
10 THz, μ
x
x
=
4.9, μ
x
y
=
2.9, and μ
z
z
=
0.99).
(a)
(b)
(c)
(d)Figure 4
The dependence of the reflection on a frequency and an incidence angle (β = 30°, ψ = 50°), d
1
=
0.78 mm, and d
1
/
d
2 = 1.3125.It also is clear (Figure3) that two of the four wavenumbers can be real and two others can be complex for the certain parameters of the structure (10
°
<
α
<
16
°). It should also be noted that all four wave numbers have positive (negative) real part (Figure 3) in the certain range of the incidence angles (α
>
16
°). Such unusual combinations of wavenumbers are possible only in the case of this complex geometry of the problem. Similarly, it also is verified that the field components are different for the incidence angles α and −α.Obviously, this kind of wavenumber affects reflection properties of a medium. So now we analyze the dependence of the reflection coefficient on a frequency and an anisotropy axis orientation. One of the possible variants of this dependence is presented in Figure4. The study shows that there is no regularity in the reflection coefficient dependence on a frequency and an orientation of the anisotropy axis. Here a very wide range of total reflection even for small incidence angles can be seen. Their existence is determined by multiple reflection phenomenon, existence of four eigenwaves with different wavenumbers and propagation angles in each layer, presence of complex waves, and the losses in the medium. The resonant character of the reflection coefficient is determined by the effects described above also. Note that there are no sharp edges of the reflection coefficient in the neighborhood of the unity. It also is shown (Figure 5) that the dependence of the reflection coefficient of a frequency is not symmetric with respect to α = 0. This is due to the asymmetric dependence of wave numbers and the field components on the incidence angle.Figure 5
The dependence of the reflection coefficient on the incidence angle forf
=
6.35
·
10
12, ψ = 50°, d
1
=
0.78 mm, and d
1
/
d
2 = 1.3125; (1) β = 30° solid line and (2) β = −30° dash line.


## 5. Operation Principle of the Isolator

The operation principle of the device is based on an angle dependence of medium properties (in particular, the reflection coefficient) but not on the Faraday effect. Wherein for the certain structure parameters and an orientation of the anisotropy axis (an external magnetic field direction) it is possible to obtain the total reflection at an incidence angleα and the minimal reflection at an incidence angle -
α. Analogously, it is possible to choose an orientation of the external magnetic field for the given structure parameters such that it is possible to obtain the total reflection for -
α and the minimal reflection for α. Thus, changing a direction of wave propagation (α or -
α) can be accomplished by changing a direction of an external magnetic field. Analogously it is possible to find the two frequencies that give such an effect.Now let us describe, for clarity, the operation principle of a bidirectional double-frequency isolator shown in Figure1. If a voltage is applied to solenoid 1, a signal transmits from antenna 1 to antenna 2 without attenuation at a frequency f
1 as the reflection coefficient is equal to unite, a signal cannot propagate from antenna 2 to antenna 1 at this frequency as the reflection coefficient is about zero, and a signal transmits from antenna 2 to antenna 1 without attenuation at a frequency f
2; however, inverse propagation is not possible at this frequency. When voltage is applied to solenoid 2, a signal transmits from antenna 2 to antenna 1 substantially without attenuation at a frequency f
1, and it transmits from antenna 1 to antenna 2 without attenuation at a frequency f
2. Circular aperture of the antennas is used because the Faraday effect (rotation of a polarization plane) is taken into account.


## 6. The Parameters of the Isolator

A very important characteristic of an isolator is its amplitude response. During the researches, numerous calculations have been carried out and the optimal structure for isolator is chosen. This structure includes the 12 double-layered periods. The first layer of the structure is FeF2 and the second one is MnO. The antiferromagnetic resonance frequencies in the considered case are 1.1 THz and 1.8 THz for the first and second layers accordingly. The dependence of the element on a frequency is given in Figure 2. Here the operating frequency f
=
6.35 THz and the angles β = ±30° and ψ = 50° are also taken. In our problem, the operating frequency is essentially higher than the resonant one. The losses in FeF2 and MnO can be neglected at this frequency (Figure 2). Here we consider an above resonance frequency and we have μ
x
x = 1.415 and μ
x
x = −0.42 for the first layer of the period at this frequency and μ
x
x = 1.2 and μ
x
x = −0.63 for the second layer. In the given work we consider a TE incident wave. The dependence of the reflection coefficient on an incidence angle for this structure (the amplitude characteristic) is presented in Figure 5. It is seen that the reflection coefficient is minimal (R
=
0.05) at α
=
76.4
° and the one is equal to unit at α = −76.4° for β = 30° (the solid line).Thus, this structure has isolator properties and it passes a signal only in a forward direction and does not pass a signal in an opposite direction. The stopband for an opposite direction is 9.6° at the level 0.707, 15.8 at the level 0.5, and 27.7 at the level 0.1. Now, let us choseβ
=
-
30
° (the dash line). For this case, the minimum of the reflection coefficient (R
=
0.03) is at α
=
-
76.4
°.The frequency dependence of the reflection coefficient is presented in Figure6. This graph is presented for the incidence angle α = 76.4°. It is obvious that this structure has bad frequency selective properties and low squareness coefficient. The bandwidth is 1.88 THz (29.38%) at the 0.707 level (a half power level), and the bandwidth is 3.13 THz (48.9%) at the 0.5 level of the maximum. The squareness is 1.66 at levels 0.5 and 0.707. It also is verified that changing orientation of magnetic field allows us to shift a reflection minimum in a small range. And changing value of this field does not give us any useful effect.Figure 6
The dependence of reflection coefficient on a frequency atα = 76.4°, d
1
=
0.78 mm, and d
1
/
d
2 = 1.3125.The results demonstrating the possibility of mechanical tuning of the isolator by changing an angle of an anisotropy axis inclination are shown in Figure7. The minimum is shifted from α
=
72
° to α
=
82
° if an inclination angle is changed from β
=
25
° to β
=
40
°. Simultaneously, the angle bandwidth is narrowed and it is equal to 28°, 10°, and 6° correspondingly. Decreasing the inclination angle less than β
=
25
° and increasing it more than β
=
50
° leads to disappearance of the angle selective properties of the structure.Figure 7
Dependence of the reflection coefficient on an incidence angle for difference inclination anglesβ, d
1
=
0.78 mm, d
1
/
d
2
=
1.3125, and f
=
6.35
·
10
12.Figure8 shows the results of calculations, allowing us to improve the described isolator. It is seen in Figure 8 that the reflection minimum is observed at α
=
76.4
° when the frequency is 6.35 GHz and the reflection minimum corresponds to α
=
-
76.4
° at the frequency 4.05 GHz. Thus, the structure can be used in duplex terahertz systems as a two-way isolator. Moreover, when the direction of the magnetic field is switched from ψ
=
50
° to ψ
=
-
50
° the above minima are reversed. Therefore, if voltage is applied to solenoid 1 then a signal transmits from antenna 1 to antenna 2 substantially without attenuation at frequency 6.35 GHz, and a signal transmits from antenna 2 to antenna 1 without attenuation at frequency 4.05 GHz. If voltage is applied to solenoid 2 then a signal transmits from antenna 2 to antenna 1 substantially without attenuation at frequency 6.35 GHz, and a signal transmits from antenna 1 to antenna 2 without attenuation at a frequency 4.05 GHz.Figure 8
The dependence of the reflection coefficient on an incident angle ford
1
=
0.78 mm,  d
1
/
d
2 = 1.3125, β
=
30
°; (1) ψ
=
50
° and (2)  ψ
=
-
50
°.


## 7. Conclusions

Electromagnetic properties of a stratified anisotropic slab with an arbitrary orientation of the anisotropy axis are studied analytically and numerically. The translation and reflection matrices methods are used. The analytical expression of the translation matrix is obtained in the convenient compact form for the first time. The dependence of the eigenwave wavenumbers in an anisotropic slab and reflection coefficient on an anisotropy axis orientation and an incidence angle is investigated. It is shown that the structure shows nonreciprocal properties. It is important that these properties are not based on the Faraday effect. Although this effect is still inherently present in an anisotropic structure, the controlled two-way dual-frequency (duplex) isolator is presented here for the first time. Its operation principle is based on the properties of an anisotropic structure described here for the first time. The very complicated geometry of an anisotropic slab is used in this isolator. But this geometry (an arbitrary orientation of magnetic field) allows us to obtain the necessary properties of the isolator and to control a resonant frequency and an angle of a wave transmission mechanically.



---

*Source: 1468508-2016-01-17.xml*