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

Graphene-based metamaterial transmitarray antenna design for the generation of tunable orbital angular momentum vortex electromagnetic waves

Open Access Open Access

Abstract

A graphene-based metamaterial transmitarray antenna is proposed to generate tunable orbital angular momentum (OAM) vortex waves in terahertz. Theoretical design of the transmitarray antenna has been developed by using the transmission line network model, and a multilayer graphene-based metamaterial element has been designed. By changing the chemical potentials of the graphene sheets, the 360° transmission phase range of the element is achieved in a broad band from 4.2 THz to 5.6 THz. By arranging the metamaterial element into a transmitarray, the OAM waves with tunable modes including l = 0, ±1, ±2 and the mode purity greater than 0.96 are generated. Simulation results are given to demonstrate good performance of the proposed design, which provides a feasible way to the efficient generation and manipulation of the OAM vortex waves in terahertz.

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

1. Introduction

Orbital angular momentum (OAM) vortex waves have attracted increasing attentions in recent years due to its promising potential ability in wireless communication. Different from the conventional plane waves, the OAM vortex beam has an azimuthal phase of exp(jlφ) and its Poynting vector abides by a spiral trajectory around the optical axis, where l is the topological charge [1]. With the theoretically unlimited orthogonal eigenstates, the OAM beam has been considered to boost communication capacity.

There have been many attempts to generate the OAM beams including spiral phase plate [2], holographic diffraction gratings [3], spiral reflectors [4], and antenna arrays [5]. Recent progress in metamaterials has enabled new ways to manipulate the OAM beams in microwave and optics [6]. By properly designing sub-wavelength topology element, various the OAM waves including multiple beams and different modes [7] and non-diffractive beams [8] are generated. However, in these reported designs the OAM modes cannot be changed, once the metamaterial-based structures are designed. In recent years, graphene has caused tremendous interests owing to unique electrical and mechanical properties. The graphene has the highest mobility in known materials. By changing electrochemical potential, the conductivity of the graphene can be flexibly tuned. With the tunable feature, the graphene becomes a natural candidate for the designs of the tunable devices including the absorbers [9], reflectarrays [10], lens [11], modulation [12], etc.

The contribution of the work is to develop a graphene-based metamaterial transmitarray antenna to flexibly generate the OAM vortex beam with the tunable modes. A design method based on a transmission line network model has been proposed to design a stacked graphene-based metamaterial unit cell. By applying different external gate voltages to the graphene sheets, the transmission phase range of the designed unit cell can cover 360°, which distinguishes the proposed design from the conventional metamaterial/metasurface based designs in which the phase variation is achieved by either adjusting the dimensions of the unit cell [7] or rotating the unit cell [8]. By arranging the unit cell into an array along the azimuthal direction, the resultant transmitarray has a helical phase profile of exp(jlφ). Illuminated by a horn antenna, the transmission wave through the transmitarray becomes an OAM vortex wave. As the chemical potentials of the graphene sheets vary with the change of the applied gate voltages, the OAM wave with different modes including 0, ±1, and ± 2 can be dynamically obtained in a wide frequency band from 4.2 THz to 5.6 THz.

2. Design of metamaterial transmitarray

Figure 1 shows the proposed metamaterial unit cell and the corresponding transmission line network model. The metamaterial unit cell is a four-layer configuration, each of the four layers is same and consists of a graphene-quartz-graphene sandwich structure. The side length of the graphene sheet is 20 µm, and the thickness of the quartz hd is 5 µm. An air layer with a thickness of h0=10 µm is placed between two adjacent sandwich structures. The relative permittivity of the quartz is 3.78. Considering that the thickness of the monolayer graphene is about 0.335 nm, the monolayer graphene is regarded as an infinitesimally thin conductive sheet. The corresponding conductivity can be analytically approximated as according to Kubo formula [1315]

$$\sigma = {\sigma _{{\mathop{\textrm {int}}} ra}} + {\sigma _{{\mathop{\textrm {int}}} er}},$$
in which intra-band conductivity σintra and inter-band conductivity σinter are given as, respectively
$$\left\{ \begin{array}{l} {\sigma_{{\mathop{\textrm {int}}} ra}} \approx i\frac{{{e^2}{k_B}T}}{{\pi {\hbar^2}({\omega + i2\Gamma } )}}\left[ {\frac{{{\mu_c}}}{{{k_B}T}} + 2\ln \left( {{e^{ - \frac{{{\mu_c}}}{{{k_B}T}}}} + 1} \right)} \right],\\ {\sigma_{{\mathop{\textrm {int}}} er}} \approx i\frac{{{e^2}}}{{4\pi \hbar }}\ln \left[ {\frac{{2|{{\mu_c}} |- ({\omega + i2\Gamma } )\hbar }}{{2|{{\mu_c}} |+ ({\omega + i2\Gamma } )\hbar }}} \right]. \end{array} \right.$$

 figure: Fig. 1.

Fig. 1. The designed metamaterial unit cell and the corresponding TL network model.

Download Full Size | PDF

Here e is electron charge, kB is Boltzmann constant, ℏ is reduced Planck’s constant, ω is radian frequency, T = 300 K is room temperature, μc is defined as chemical potential of graphene, and Γ=1/(2τ) is charged particle scattering rate with the momentum relaxation time τ = 2ps.

With the designed unit cell, the corresponding TL network is shown in Fig. 1. Each graphene sheet is considered as a shunt admittance Ys(i)=1/σ(i) (i = 1, 2). The quartz and the air layers, respectively, are equivalent to the TLs with the phase constants of βd and β0 and the lengths of hd and h0. Each part in the unit cell are cascaded together to achieve the whole TL network model. The transfer matrix of the whole network [A] [16] can be written as

$$[A ]= \left[ {\begin{array}{{cc}} A&B\\ C&D \end{array}} \right] = [{{A_1}} ]\cdot [{{A_2}} ]\cdot [{{A_1}} ]\cdot [{{A_2}} ]\cdot [{{A_1}} ]\cdot [{{A_2}} ]\cdot [{{A_1}} ],$$
where
$$[{{A_1}} ]= \left[ {\begin{array}{{cc}} 1&0\\ {{Y_{s(1)}}}&1 \end{array}} \right]\left[ {\begin{array}{{cc}} {\cos {\beta_d}{l_d}}&{j{{\sin {\beta_d}{l_d}} \mathord{\left/ {\vphantom {{\sin {\beta_d}{l_d}} {{Y_d}}}} \right.} {{Y_d}}}}\\ {j{Y_d}\sin {\beta_d}{l_d}}&{\cos {\beta_d}{l_d}} \end{array}} \right]\left[ {\begin{array}{{cc}} 1&0\\ {{Y_{s(2)}}}&1 \end{array}} \right],$$
$$[{{A_2}} ]= \left[ {\begin{array}{{cc}} {\cos {\beta_0}{l_0}}&{j{{\sin {\beta_0}{l_0}} \mathord{\left/ {\vphantom {{\sin {\beta_0}{l_0}} {{Y_0}}}} \right.} {{Y_0}}}}\\ {j{Y_0}\sin {\beta_0}{l_0}}&{\cos {\beta_0}{l_0}} \end{array}} \right].$$
Here Yd and Y0 are wave admittances of the quartz and the air, respectively. Therefore, the transmission coefficient S21 of the network can be obtained as [17]
$${S_{21}} = \frac{2}{{A + B \cdot {Y_0} + {C \mathord{\left/ {\vphantom {C {{Y_0} + D}}} \right.} {{Y_0} + D}}}}.$$
In order to verify the effectiveness of the proposed TL network, the comparison of S21 between the TL network and full-wave (FW) simulation for different chemical potentials of the graphene sheets is given in Fig. 2. A good agreement between them is observed for different μc1 and μc2. Here μc1 and μc2 are the chemical potentials of the top and bottom graphene sheets in each layer, respectively. Note that in the FW simulation of the unit cell, an infinite periodic model including floquet port and periodic boundary is used, as shown in Fig. 1. With the TL network model, the transmission phases range of the unit cell are calculated in a wide frequency band from 4.2 THz to 5.6 THz when the chemical potentials of the graphene sheets μc1c2 vary from 0 eV to 1 eV, as shown in Fig. 3. It can be observed that in the frequency band of 4.2∼5.6 THz, the 360° transmission phase range of the proposed unit cell is achieved.

 figure: Fig. 2.

Fig. 2. Comparison of S21 between the FW simulation and the TL network model. (a). Magnitude of S21. (b) Phase of S21.

Download Full Size | PDF

 figure: Fig. 3.

Fig. 3. Transmission phase of the unit cell for possible μc1 and μc2 in a wide band.

Download Full Size | PDF

The proposed metamaterial unit cell is arranged into an array along the azimuthal direction, as shown in Fig. 4. The OAM vortex wave can be generated by the transmitarray as the variation of the transmission phase with the azimuthal angle φ continuously covers 360°. However, it is infeasible for a continuous phase variation. A good approximation is to discretize the transmission phase of 360° into N piecewise segments, each of which covers the phase of 360°/N. Here 8 segments are used to approximately cover the transmission phase of 360°, as shown in Fig. 4. The dimensions of the transmitarray are as follows: L1 = 1400 µm, L2 = 1410 µm, and Lg = 0.5 µm. Table 1 gives the chemical potentials of 8 segments and the corresponding S21 at 4.9 THz. It can be observed that the transmission phase difference between two adjacent segments is 45° and the total phase range of the 8 segments covers 360°. Moreover, the magnitude of S21 is better than −2 dB, which means the small transmission loss of the designed transmitarray. A standard circular horn antenna as the feeding source is placed at the top of the transmitarray with a focal length to diameter ratio of 0.29. A square observation plane with a side length of 3000 µm is located at the bottom of the transmitarray with a distance of 3150 µm, on which the phase, the magnitude, the complex magnitude of the electric field are solved. The phase distribution corresponding to the OAM wave with the mode of l=+1 is observed on the observation plane at 4.9 THz, as shown in Fig. 5.

 figure: Fig. 4.

Fig. 4. The graphene-based metamaterial transmitarray.

Download Full Size | PDF

 figure: Fig. 5.

Fig. 5. OAM vortex wave with tunable modes at 4.9 THz.

Download Full Size | PDF

Tables Icon

Table 1. Chemical potentials of the graphene in 8 segments and the corresponding S21.

To tune the mode of the OAM wave, the variable gate voltages Vd are applied between each graphene sheet of the transmitarray and the ground. Here the gold as a ground is placed around the transmitarray, as shown in Fig. 4. Moreover, a small gap of Ly = 0.5 µm is introduced between two adjacent graphene sheets located in the same plane to isolate different gate voltages. It is worthwhile pointing out that the use of the small gap does not affect the performance of the transmitarray. With the voltage, the chemical potential μc of the graphene can be estimated by the following formula [18]:

$$|{{\mu_c}} |\approx \hbar {v_F}{\{{\pi {a_0}|{{V_d} - {V_{Dirac}}} |} \}^{1/2}},$$
in which, vF9×105 m/s represents the Fermi velocity of the Dirac fermions, a0 9×1016 m−2V1 is estimated through a single capacitor model, and VDirac=0.8 V denotes the Dirac voltage offset caused by the natural doping.

When the chemical potentials in the segments 2, 4, 6, and 8 are changed to be same as those in the segments 1, 3, 5, and 7, respectively, the OAM wave with l = 2 is generated, as shown in Fig. 5. If the chemical potentials in the 8 segments are same, the OAM wave with l = 0, i.e., the plane wave, is obtained. If we tune the gate voltages to exchange the chemical potentials of the segments 2, 3, and 4 with those of the segments 8, 7, and 6, respectively, the OAM wave with l=−1 is generated, as shown in Fig. 5. In addition, a doughnut-shaped radiation patterns characterized by the OAM wave with different modes can be observed in Fig. 5. In order to evaluate the performance of the OAM wave produced by the transmitarray, a spectral analysis of Fourier transform [19,20] is implemented. The magnitude and phase of the OAM wave are sampled along a circle with a radius of 1000 µm, as shown in the inset of Fig. 6. The mode purities for the OAM wave with l = 1 and l = 2 are larger than 96%. Furthermore, as shown in Fig. 3, the 360° transmission phase in a wide frequency band means that the OAM wave can be generated in a broad band. As the chemical potentials of the graphene sheet are chosen in the frequency band from 4.2 THz to 5.6 THz according to Fig. 3, the OAM wave with l = 1 can be achieved, as shown in Fig. 7. Similar to the case at 4.9 THz, the OAM wave with other modes can be generated in the wideband of 4.2∼5.6 THz.

 figure: Fig. 6.

Fig. 6. Mode purity of the OAM wave.

Download Full Size | PDF

 figure: Fig. 7.

Fig. 7. The OAM wave with l = 1 in a wideband.

Download Full Size | PDF

3. Conclusion

In this paper, a graphene-based metamaterial transmitarray is designed to generate the tunable OAM vortex wave in a wide band covering 4.2∼5.6 THz. With the transmission line network analysis, a multilayer metamaterial element is designed and the 360° transmission phase range can be achieved in a broad band. With the proposed multilayer element, the transmitarray is designed to produce the high-purity OAM wave with the tunable modes of l = 0, ±1, ±2 in a broad band. The simulated results are given to verify the proposed design.

Funding

National Natural Science Foundation of China (61771359); Natural Science Basic Research Plan in Shaanxi Province (2018JM6006).

References

1. L. Allen, M. W. Beijersbergen, R. J. C. Spreeuw, and J. P. Woerdman, “Orbital angular-momentum of light and the transformation of Laguerre–Gaussian laser modes,” Phys. Rev. A 45(11), 8185–8189 (1992). [CrossRef]  

2. R. Niemiec, C. Brousseau, K. Mahdjoubi, O. Emile, and A. Menard, “Characterization of an OAM flat-plate antenna in the millimeter frequency band,” Antennas Wirel. Propag. Lett. 13, 1011–1014 (2014). [CrossRef]  

3. F. Tamburini, E. Mari, T. Bo, C. Barbieri, and F. Romanato, “Experimental verification of photon angular momentum and vorticity with radio techniques,” Appl. Phys. Lett. 99(20), 204102 (2011). [CrossRef]  

4. F. Capasso, P. Genevet, J. Lin, and M. A. Kats, “Holographic detection of the orbital angular momentum of light with plasmonic photodiodes,” Nat. Commun. 3(1), 1278 (2012). [CrossRef]  

5. B. Thidé, H. Then, J. Sjöholm, K. Palmer, J. Bergman, T. D. Carozzi, Y. N. Istomin, N. H. Ibragimon, and R. Khamitova, “Utilization of photon orbital angular momentum in the low-frequency radio domain,” Phys. Rev. Lett. 99(8), 087701 (2007). [CrossRef]  

6. S. Yu, L. Li, G. Shi, C. Zhu, X. Zhou, and Y. Shi, “Design, fabrication, and measurement of reflective metasurface for orbital angular momentum vortex wave in radio frequency domain,” Appl. Phys. Lett. 108(12), 121903 (2016). [CrossRef]  

7. S. Yu, L. Li, G. Shi, C. Zhu, and Y. Shi, “Generating multiple orbital angular momentum vortex beams using a metasurface in radio frequency domain,” Appl. Phys. Lett. 108(24), 241901 (2016). [CrossRef]  

8. K. Zhang, Y. Y. Yuan, D. W. Zhang, X. M. Ding, B. Ratni, S. N. Burokur, M. J. Lu, K. Tang, and Q. Wu, “Phase-engineered metalenses to generate converging and non-diffractive vortex beam carrying orbital angular momentum in microwave region,” Opt. Express 26(2), 1351 (2018). [CrossRef]  

9. Y. Zhang, Y. Shi, and C. H. Liang, “A broadband tunable graphene-based metamaterial absorber,” Opt. Mater. Express 6(9), 3036 (2016). [CrossRef]  

10. Y. Shi and Y. Zhang, “Generation of wideband tunable orbital angular momentum vortex waves using graphene metamaterial reflectarray,” IEEE Access 6, 5341–5347 (2018). [CrossRef]  

11. A. Andryieuski, A. V. Lavrinenko, and D. N. Chigrin, “Graphene hyperlens for terahertz radiation,” Phys. Rev. B 86(12), 121108 (2012). [CrossRef]  

12. Y. Ding, X. Zhu, S. Xiao, H. Hu, L. H. Frandsen, N. A. Mortensen, and K. Yvind, “Effective electro-optical modulation with high extinction ratio by a graphene–silicon microring resonator,” Nano Lett. 15(7), 4393–4400 (2015). [CrossRef]  

13. V. P. Gusynin, S. G. Sharapov, and J. P. Carbotte, “Magneto-optical conductivity in graphene,” J. Phys.: Condens. Matter 19(2), 026222 (2007). [CrossRef]  

14. V. P. Gusynin, S. G. Sharapov, and J. P. Carbotte, “Sum rules for the optical and hall conductivity in graphene,” Phys. Rev. B: Condens. Matter Mater. Phys. 75(16), 165407 (2007). [CrossRef]  

15. G. W. Hanson, “Dyadic greens functions and guided surface waves for a surface conductivity model of graphene,” J. Appl. Phys. 103(6), 064302 (2008). [CrossRef]  

16. D. M. Pozar, Microwave Engineering (Wiley, 2012).

17. J. S. Hong and M. J. Lancaster, Microstrip Filters For RF/Microwave Applications (Wiley, 2001).

18. A. Vakil and N. Engheta, “Transformation optics using graphene,” Science 332(6035), 1291–1294 (2011). [CrossRef]  

19. M. J. Strain, X. L. Cai, J. W. Wang, J. B. Zhu, D. B. Phillips, L. F. Chen, M. Lopez-Garcia, J. L. O’Brien, M. G. Thompson, M. Sorel, and S. Yu, “Fast electrical switching of orbital angular momentum modes using ultra-compact integrated vortex emitters,” Nat. Commun. 5(1), 4856 (2014). [CrossRef]  

20. B. Jack, M. J. Padgett, and S. F. Arnold, “Angular diffraction,” New J. Phys. 10(10), 103013 (2008). [CrossRef]  

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 (7)

Fig. 1.
Fig. 1. The designed metamaterial unit cell and the corresponding TL network model.
Fig. 2.
Fig. 2. Comparison of S21 between the FW simulation and the TL network model. (a). Magnitude of S21. (b) Phase of S21.
Fig. 3.
Fig. 3. Transmission phase of the unit cell for possible μc1 and μc2 in a wide band.
Fig. 4.
Fig. 4. The graphene-based metamaterial transmitarray.
Fig. 5.
Fig. 5. OAM vortex wave with tunable modes at 4.9 THz.
Fig. 6.
Fig. 6. Mode purity of the OAM wave.
Fig. 7.
Fig. 7. The OAM wave with l = 1 in a wideband.

Tables (1)

Tables Icon

Table 1. Chemical potentials of the graphene in 8 segments and the corresponding S21.

Equations (7)

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

σ = σ int r a + σ int e r ,
{ σ int r a i e 2 k B T π 2 ( ω + i 2 Γ ) [ μ c k B T + 2 ln ( e μ c k B T + 1 ) ] , σ int e r i e 2 4 π ln [ 2 | μ c | ( ω + i 2 Γ ) 2 | μ c | + ( ω + i 2 Γ ) ] .
[ A ] = [ A B C D ] = [ A 1 ] [ A 2 ] [ A 1 ] [ A 2 ] [ A 1 ] [ A 2 ] [ A 1 ] ,
[ A 1 ] = [ 1 0 Y s ( 1 ) 1 ] [ cos β d l d j sin β d l d / sin β d l d Y d Y d j Y d sin β d l d cos β d l d ] [ 1 0 Y s ( 2 ) 1 ] ,
[ A 2 ] = [ cos β 0 l 0 j sin β 0 l 0 / sin β 0 l 0 Y 0 Y 0 j Y 0 sin β 0 l 0 cos β 0 l 0 ] .
S 21 = 2 A + B Y 0 + C / C Y 0 + D Y 0 + D .
| μ c | v F { π a 0 | V d V D i r a c | } 1 / 2 ,
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.