1. Introduction
Planar photonic crystal structures with various combinations of their constituent layers have attracted close attention of researchers due to their capabilities for effective control of optical radiation [
1,
2,
3,
4,
5]. Control of reflection and transmission spectra for such structures, as a rule, is achieved by including layers in their composition, the material parameters of which depend on easily variable external factors (electric and magnetic fields, angle of incidence and polarization of incident radiation) [
6,
7,
8,
9,
10]. Photonic crystal structures are used to create a wide class of radiation control devices (filters, switches, modulators), performing the function of filtering and amplifying various types of propagating waves [
11,
12,
13,
14].
The presence of the photonic band gaps in the photon spectrum is characteristic of one-dimensional layered periodic structures. Violation of periodicity in the structure leads to the appearance of a narrow transmission band in the photonic band gap. Microcavities are similar structures, the formation of which requires two Bragg mirrors (BMs) with a change in the order of layers in one of them. The working cavity (the area between the mirrors), as a rule, is filled with an active medium. Standing waves are established in the working cavity, and transmission resonances (defect modes) arise in the photonic band gap if the distance between the BMs is a multiple of half the propagating radiation wavelength.
The contrast of the photonic band gap and the defect mode in the spectrum significantly depends on the number of periods in the BMs; with their increase, the contrast increases. In this regard, it is of interest to consider the properties of a photonic cell (PC) with a minimum number of periods in BMs, which can also be considered as a microcavity with fairly effective control of the radiation passing through it (or reflected). Graphene and various planar structures based on it can be considered the most promising photonic materials. Graphene is a 2D structure and has a number of unique properties: zero band gap, high electron mobility, high optical transparency [
15,
16,
17,
18,
19]. Graphene may have the properties of an amplifying medium when excited and creates an inversion in its energy structure [
20,
21,
22]. The active state of graphene can be obtained through the drift of charge carriers as a result of current pumping. Doped and inverted graphene is discussed in paper [
21,
22], decorated graphene is discussed in paper [
23]. In recent years, structures called “photonic graphene” have been of great interest. These are two-dimensional metasurfaces made of nano- or microparticles located in space according to a certain law at distances shorter than the wavelength of propagating radiation. An important property of such gratings is tunability under external influence, the role of which can be played by the intensity of the light wave, polarization, magnetic field, temperature [
24,
25].
In this work, the photonic properties of one of the possible symmetrical graphene-containing photonic cells are studied based on a numerical analysis of the general relationships for the transmission (reflection) coefficients and polarization characteristics of a planar layered magnetically active structure, in which the active state of graphene is achieved due to the drift of charge carriers as a result of current pumping. It has been shown that it is possible to effectively control the transmission and reflection spectra, as well as the polarization of transmitted radiation, by changing the Fermi energy of graphene, by changing the external magnetic field, by changing the thickness of the dielectric (the gap between the graphene layers), and also by changing the orientation of the polarization plane incident radiation.
2. Geometry and Material Parameters of a Photonic Cell
From the many possible types of PCs, we study a cell
, consisting of two layers of a magnetic semiconductor
, two graphene monolayers
, separated by a dielectric layer
, the thickness of which can be changed during the experiment (such a dielectric can be an air gap). A PC is shown schematically in
Figure 1.
As a magnetically active material, we have chosen a well-studied magnetic semiconductor
, which combines the properties of a semiconductor (sufficiently high conductivity, can be
p- and
n-type depending on the production method) and a magnetic (sufficiently high magnetization) [
26]. This material is characterized by cyclotron and ferromagnetic resonances with a narrow width of the resonance curve, which significantly expands the possibilities of its practical application [
27].
The permittivities and permeabilities of such a material are tensor quantities:
The components of the dielectric permittivities tensor are given by the relations [
28]:
where plasma and cyclotron frequencies are introduced
and
. For numerical modeling, we will use the parameters of the above
p-type magnetic semiconductor:
is the lattice part of the permittivity,
and
is the concentration of carriers
cm
, their effective mass
,
, collision frequency
[
26,
28]. The effective permittivity is
taking into account the transverse orientation of the external magnetic field with respect to the wave vector
. Note the resonant nature of the dependencies
. For each value
, in a narrow region near the frequency
there is a sharp increase in the imaginary part of the effective permittivity, which should lead to sharp dips in the transmission spectrum at these frequencies.
The non-zero components of the permeability tensor
,
,
have the form [
29]:
where
,
is the saturation magnetization,
,
is the magnetic resonance linewidth,
is the gyromagnetic ratio. Effective permeability for the considered geometry
.
Figure 2 shows the frequency dependences of the imaginary part of the effective permittivity and permeability obtained at
kOe (black, red, blue, green curves), which correspond to absorption lines at cyclotron and ferromagnetic resonances. The maxima of these dependencies lie in the low-frequency terahertz region for permittivities and in the microwave region for permeability. The maxima of these dependences correspond to cyclotron and magnetic resonances. It is possible to control the material parameters of an effective medium using a magnetic field; this follows from the given dependencies. In this case, an increase in the external magnetic field shifts the resonance curves to higher frequencies.
The terahertz response of a graphene monolayer carrying a direct current for frequencies below the electron-electron collision frequency is determined by its surface conductivity
. Due to the strong electron-electron interaction in graphene, it is correct to use the hydrodynamic approach at the indicated frequencies. The expression for the surface conductivity of graphene, valid for the low frequency region
THz and obtained in the hydrodynamic approximation [
30,
31,
32,
33], taking into account the normal incidence of the wave on the layer, has the form
Here the surface conductivity of graphene is introduced
,
is the Fermi energy in the absence of carrier drift,
is a parameter,
and
u are the Fermi speed and the drift speed, which is determined by the magnitude of the applied static electric field (
is the potential difference between the graphene layers). The real part of the graphene conductivity becomes negative at drift speeds greater than the phase speed of the electromagnetic wave, and due to the Vavilov-Cherenkov effect in the terahertz range, amplification is achieved [
31,
32].
Figure 3 shows the dependence of the conductivity real part of graphene on the frequency (a), obtained on the basis of Equation (
4) for the values
,
meV (red, green, blue curves) and
m/s,
ps.
Figure 3b,c shows the conductivity maps of the dependences
and
, which make it possible to determine the conductivity value for any values of the specified parameters. Dotted lines separate positive and negative values
. Thus, the real part of the conductivity
takes negative values and graphene becomes an amplifying medium in the frequency range
THz at non-zero Fermi energy values. The combination of the frequency ranges of graphene amplification and the resonances of the material parameters of a magnetic semiconductor indicates that using an external magnetic field it will be possible to control the optical characteristics of a photonic cell.
For the permittivity of graphene, we can write , since high-frequency conductivity and permittivity are related to each other by Maxwell’s equations. This means that the imaginary part of the permittivity of graphene at may turn out to be negative, which indicates the possibility of implementing amplification with a graphene monolayer in the specified range. Graphene layers can be in both unexcited and active states. In unexcited graphene, the Fermi energy (chemical potential) is at the Dirac point and is equal to zero. In this case, the valence band is completely filled, the conduction band is completely free, and there is no band gap. It is possible to shift the chemical potential into the conduction band or valence band by applying a voltage of different polarity to a graphene sheet, and thereby change the surface conductivity of the graphene. A planar magnetic field applied in the plane of monolayers does not affect the electronic subsystem of graphene due to its 2D geometry.
3. Amplitude Characteristics of the Photonic Cell
Consider the transformation of transmission spectra with changes in the Fermi energy and magnetic field to identify the spectral features that arise when linearly polarized radiation passes through a PC. We assume that the incident wave falls on the PC normally to the plane of its layers and propagates along the symmetry axis of the structure. In this case, the inclination angle of the polarization plane, measured from the axis, can vary from 0 to . For numerical analysis, we select the geometric parameters of the cell: the thickness of each of the side layers of the magnetic semiconductor is assumed to be equal to m, the thickness of the graphene monolayer is assumed to be nm, and the thickness of the dielectric (or air gap) is further considered to be a parameter. The wavelength m corresponds to the operating frequency THz. For the air gap , the optical and geometric thickness of the working cavity are the same.
The transfer matrix connecting the wave field at the input
and output
of the photonic cell has the form
. Here
are the transfer matrices of individual layers, which in the plane wave approximation (valid for the above transverse dimensions of the layers) can be represented as follows:
where
,
are the propagation constants in the corresponding layers,
,
and
c are the frequency and speed of the wave in vacuum. For dielectric layers
,
, for graphene layers
,
, for semiconductor layers
,
in the case of a TE wave and
,
in the case of a TM wave.
The amplitude transmission and reflection coefficients for the entire structure are determined through the matrix elements of the transfer matrix
G [
23]:
If the polarization plane of the incident wave is oriented at an angle
to the direction of the external magnetic field (
), we represent the electric field of the wave incident on the photonic cell as follows:
At the exit from the photonic cell, the electric field is determined by the expression:
Taking into account the obtained expressions, the energy transmission and reflection coefficients can be presented as follows:
When taking into account absorption (amplification) in the layers, the fraction of energy absorbed by the structure is determined by the quantity. If there is only absorption in the structure, if there is both absorption and amplification, the integral absorption coefficient can be either greater or less than zero.
Figure 4 for two photonic cells with a working cavity filled with a dielectric SiO
2 (a–c) and air (d–f), shows the frequency dependences of the transmission coefficients for incident waves with the orientation of the plane of polarization
(green, blue, red curves), which become in the cell
and
waves of elliptical polarization, respectively. We chose silicon dioxide (SiO
2) as a dielectric, which is an optically isotropic material and has high optical transparency.
The optical thicknesses of the working cavity in both photonic cells are the same and equal
, where
m corresponds to the operating frequency
THz, values
,
. For dielectric SiO
2 in the terahertz range
[
34,
35], the geometric thickness of the dielectric must be equal to
so that the optical thicknesses of the dielectric and the air gap are equal. The optical anisotropy of a PC is determined by the presence of magnetically active layers
with a magnetization orientation directed along the
axis (as in the applied magnetic field). It can be seen that significant changes occur in the spectrum with changes in the optical thickness of the dielectric (air gap). First of all, with decreasing
, the number of spectral peaks in the frequency range under consideration also decreases. With
, a band gap appears in the spectrum, similar to the well-developed photonic band gap in structures with a large number of periods. At
, a defect mode is observed at the center of the band gap (at frequency
THz), which is absent at thickness
. Thus, at
, an additional defect mode is also observed, which is split off from the right edge of the band gap. The indicated transformation of the spectrum is associated with a change in the phase relationships of the reflected waves in the photonic cell.
The nature of the given spectra also indicates that the strongest magneto-optical interaction with the structure is experienced by waves with the orientation of the polarization plane at the entrance to the photonic cell , the weakest for waves with . This is due to the fact that the first of these waves inside a photonic cell is a TM wave, the second is a TE wave. The most noticeable changes associated with the presence of graphene layers occur for waves with in the low-frequency region of the spectrum.
The spectrum of radiation passing through a photonic cell can be effectively controlled not only by the thickness of the air gap, but also by changing the energy state of graphene (i.e., its Fermi energy). Fermi energy control of graphene is based on its dependence on the external electric field strength
, which is given by [
20,
32]:
where
is the Fermi-Dirac distribution,
is the energy of the electron,
is the Fermi speed, where
eV,
b is the distance between neighboring atoms in the graphene structure,
is the dielectric constant of the substrate.
Figure 5 shows the dependence of the Fermi energy of graphene on the external constant electric field strength for
nm and the permittivity
of the substrate, respectively, for vacuum, SiO
2 and
(blue, red, green curves). It can be seen that with an increase in the external static field
applied to the graphene layer, the value of the Fermi energy increases the faster, the greater the permittivity of the substrate
. In further numerical analysis, we assumed to use as a substrate for graphene SiO
2 (green curve).
Figure 6 shows the frequency dependence of the transmission and reflection coefficients at the angles of polarization plane orientation of the incident wave
(green, red, blue curves), which were obtained for a photonic cell with
,
meV and
kOe.
It can be seen that the transformation of the spectrum occurs with increasing Fermi energy, which is expressed primarily in the appearance of amplification in the low-frequency region. In this region, the transmission coefficient, as well as the reflection coefficient, become significantly greater than unity. The amplification region shifts from the region of low to the region of higher frequencies with increasing Fermi energy, and the magnitude of the amplification decreases, which is consistent with the dependence . In this case, a significant deformation of the band gap also occurs, associated with a shift of the low-frequency edge of the amplification to the region of higher frequencies. A slight increase in the amplitude of the defect mode, which also becomes greater than unity, is observed in the center of the band gap. It is characteristic that at meV the transmission coefficients in the low-frequency region do not exceed unity, while the defect mode and the value of T at the right edge of the band gap become significantly greater than unity. The reflectance spectra in the low-frequency region also experience significant changes depending on the magnitude of the magnetic field. In the high-frequency region, the spectrum changes are insignificant. As the magnetic field increases, the most significant transformation is experienced by the spectra for waves with , which become TM modes inside the photonic cell.
Figure 7 shows a map of the transmission coefficient when a linearly polarized wave is incident on a PC for three orientations of the polarization plane
. The thickness of the air gap of the cell is taken equal to
. These maps determine the dependence of the coefficient
T on frequency and Fermi energy at external magnetic field values
kOe. For a clearer image of the amplification regions, we selected the low-frequency part of the spectrum and the transmission coefficient scale was limited to the value
. Representing the transmission spectrum with tone diagrams makes it possible to convey the continuous dependence of the coefficient
T on two parameters at once: frequency and Fermi energy. In particular, for the given values of
and
the regions of the parameters
f and
are visible, where the amplification of the radiation propagating in the photonic cell is most effectively manifested.
For example, one of the amplification regions for almost all values of has clear shapes for the parameters , THz and meV. For an angle of , the amplification regions are the geometric sum of the amplification regions for waves with polarization angles and , which is expressed in an increase in the number of amplification regions.
4. Polarization Characteristics of a Photonic Cell
To determine the polarization properties of a PC, we will assume that the plane of polarization of the incident wave is at an angle
to the direction of the magnetizing field
and is determined by the direction of the electric vector of the wave in the plane
. In the geometry under consideration, the eigenwaves for the photonic cell under study are waves of the TE and TM type with wave field components
and
, respectively. At the exit from the photonic cell, the components of the wave electric field are determined, according to Equation (
8), by the following expressions:
where the corresponding amplitudes are:
,
. To describe the polarization state of a wave passing through a photonic cell, we introduce a complex polarization variable [
36,
37]:
where
is the phase mismatch of eigenwaves when passing through a photonic cell. The main polarization characteristics of the wave are the angle of inclination
of the major axis of the ellipse to the OX axis (azimuth of the polarization ellipse) and ellipticity
, which is equal to the ratio of the semi-axes of the polarization ellipse and is determined by the ellipticity angle
[
36,
37]:
Since the transmission coefficients of TE and TM waves depend on the parameters that control the state of the structure, the polarization of the wave transmitted through the photonic cell can be controlled by changing the control parameters.
Figure 8a,b shows polarization ellipses corresponding to the magnetic field values
kOe (red, green, blue, black lines) and: polarization angles of the incident wave
, Fermi energy
meV, frequency
THz (
Figure 8a); polarization angle
and different values of frequency and Fermi energy (indicated in
Figure 8b).
It can be seen that at angles and the polarization at the exit from the photonic cell does not change, i.e., the transmitted wave remains linearly polarized with the same orientation of the polarization plane. With increasing angle , the polarization of the transmitted wave becomes elliptical with ellipticity parameters depending on frequency, Fermi energy, external magnetic field and angle . The semimajor axis of the transmitted wave ellipse is greater than unity, since the selected frequency and Fermi energy value fall into the amplification region.
The polarization of the wave passing through the photonic cell can also be controlled at a fixed
by changing other parameters, for example, the Fermi energy.
Figure 8b shows the polarization ellipses of the transmitted wave at
and various values of frequency and Fermi energy, which are indicated in the figures. The arrows in the figures indicate the direction in which the polarization ellipse is traversed as the wave propagates. Blue and black ellipses at the indicated frequencies correspond to left-handed polarized waves (the electric field vector rotates clockwise), red and green ellipses correspond to right-handed polarized waves (the electric field vector rotates counterclockwise).
Figure 9 shows contour plots that determine the dependence of ellipticity as a function of frequency and Fermi energy for two values of the thickness of the air gap between the graphene layers
. Constant values of the magnetic field
are marked in the figure, as well as lines of equal ellipticity levels. It can be seen that with an increase in the magnetic field, the transmitted right-handed polarized radiation changes to left-handed polarized; centers are also visible in which the ellipticity reaches maximum values
, which corresponds to the degeneration of the ellipse into a circle with the corresponding direction of traversal of the electric field vector.
From a comparison of the given contour plots it follows that the position of the areas with the maximum and minimum ellipticity values, as well as changes in the direction of rotation of the electric vector, significantly depends on the thickness of the air gap .
Thus, by changing the thickness of the air gap, it is also possible to control not only the parameters of the polarization ellipse, but also the rotation direction of the electric vector, i.e., type of transmitted wave.
5. Discussion
The study of layered periodic structures, which are used as functional elements, plays an important role in the development of current areas of science and technology. Therefore, close attention has been paid to the analysis of dispersion properties, reflectivity and transmission capacity of planar photonic crystal structures consisting of a large number of periods and various materials (passive and active) for many years. In this work, we study the optical properties of a photonic cell, which is a structure with a minimum number of periods, which makes it possible to implement a photonic band gap in its transmission spectrum and control its configuration and parameters using external electric and magnetic fields. For these purposes, the work considers one of the possible resonator-type photonic cells, consisting of two pairs of “magnetic semiconductor-graphene” layers and a dielectric layer between them.
The use of graphene in the structure is associated with the possibility of implementing amplification in the terahertz range, which can be achieved through several mechanisms. Thus, graphene with an inverted population of charge carriers, which can arise in graphene due to optical, injection or diffusion pumping, can act as an active medium for amplifiers and generators of THz radiation. However, for the near THz range it is promising to use graphene with charge carrier drift created by passing a high-density current through graphene without destroying it, due to technical and fundamental difficulties in creating inversion. Thus, to achieve a drift velocity
, it is necessary to create an electric field with a strength of 15 V/
m, while the current density in graphene is 173 A/cm at a Fermi energy of 200 meV [
32].
The influence of a magnetic field on the properties of a magnetic semiconductor is associated with the resonant dependence of the components of its permittivity and permeability on the magnetic field strength, which is included in the expressions for the cyclotron frequency and the ferromagnetic resonance frequency. Free electrons of a magnetic semiconductor begin to rotate around the magnetic field lines at a frequency close to the cyclotron resonance frequency, which interferes with the propagation of an electromagnetic wave inside the material. Spin oscillations will be excited in a magnetic semiconductor, which also impede the penetration of an electromagnetic wave into the medium at frequencies close to the ferromagnetic resonance frequency (microwave range), which should lead to a decrease in the transmission coefficient.
A magnetic field applied in the plane of graphene layers, due to their transverse monoatomicity, does not affect the dynamics of electrons. Since the radius of the cyclotron orbit , even at T, is much greater than the thickness of graphene, in this case the magneto-optical effect in graphene is strongly suppressed, which does not lead to the influence of the magnetic field on the transmission (reflection) coefficients.
Thus, the transformation of the spectra can be achieved both by changing the Fermi energy of graphene and its conductivity under the influence of an external electric field, and by changing the permittivity and permeability of the magnetic semiconductor layers when exposed to an external magnetic field.