Wide-Band Circularly Polarized ReflectarrayUsing Graphene-Based Pancharatnam-Berry Phase Unit-Cells for Terahertz Communication

A wide-band and high gain circularly polarized (CP) graphene-based reflectarray operating in the THz regime is proposed and theoretically investigated in this paper. The proposed reflectarray consists of a THz CP source and several graphene-based unit-cells. Taking advantages of the Pancharatnam Berry (PB) phase principle, the graphene-based unit-cell is capable of realizing a tunable phase range of 360° in a wide-band (1.4–1.7 THz) by unit-cell rotating, overcoming the restriction of intrinsic narrow-band resonance in graphene. Therefore, this graphene-based unit-cell exhibits superior bandwidth and phase tunability to its previous counterparts. To demonstrate this, a wide-band (1.4–1.7 THz) focusing metasurface based on the proposed unit-cell that exhibits excellent focusing effect was designed. Then, according to the reversibility of the optical path, a CP reflectarray was realized by placing a wide-band CP THz source at the focal point of the metasurface. Numerical simulation demonstrates that this reflectarray can achieve a stable high gain up to 15 dBic and an axial ratio around 2.1 dB over the 1.4–1.7 THz band. The good radiation performance of the proposed CP reflectarray, as demonstrated, underlines its suitability for the THz communication applications. Moreover, the design principle of this graphene-based reflectarray with a full 360° phase range tunable unit-cells provides a new pathway to design high-performance CP reflectarray in the THz regime.


Introduction
Terahertz (THz) science and technology that is undergoing an unprecedented rapid development has attracted considerable attention in the scientific community. As a promising technology, it has enormous potential in applications such as communication [1], sensing [2], imaging [3], detection [4] and spectroscopy [5]. In terms of communication applications, the THz communication system is particularly promising as it has a sufficient bandwidth that is capable of handling high data rate up to 100 Gb/s [6]. Consequently, new requirements for THz antennas have been put forward. Generally, high gain THz antennas are indispensable due to the intrinsic high path loss characteristic of THz waves in free space. Besides, considering the performance deterioration of linearly polarized antennas when they are polarization misaligned, circularly polarized (CP) antennas are often required in scenarios where multi-path effect are prominent. Therefore, high gain CP THz antennas are highly desirable.
By controlling the phase distribution on the array surface, the far-field radiation properties can be easily manipulated without bulky and lossy beamforming networks. In terms of CP reflectarray, several CP reflectarrays with excellent performances have been designed at microwave or millimeter wave frequencies [15][16][17][18][19]. It seems that reflectarrays can be chosen as a competitive candidate for high gain and CP THz applications. Nevertheless, most of these reflectarrays rely on judiciously designed metallic structures, which cannot simply scale to the THz frequencies because of the prominent skin effect of the conventional metal.The ohm loss of the metallic will immensely affect the array gain and radiation efficiency. Toavoid these drawbacks brought by the skin effect of metals at THz frequencies, there is a great demand for new types of THz antenna. Fortunately, graphene has several advantages over convention metals at the THz band. It has a high electron mobility up to 230,000 cm 2 /Vs at room temperature [20], and a low electrical resistivity about 10 −6 Ω·cm [21,22] in THz band, demonstrating lower loss than conventional metals. In addition, the surface plasmons resonant frequencies of graphene are quite lower than that in metals, which are often in optical frequencies.
Meanwhile, graphene surface plasmons exhibit extremely small wavelengths (λ/10-λ/100) and tight field confinement on the graphene sheet [23,24], while maintaining reasonably small losses in the THz band. Furthermore, the imaginary conductivity of graphene is highly tunable via chemical doping or electrical gating [25][26][27][28][29][30][31] at THz frequencies, which is impossible or inefficient if metals are used. Based on excellent physical properties of graphene at THz frequencies, several graphene-based THz antennas have been reported in recent years [32][33][34][35][36][37][38][39][40][41]. It evaluates the feasibility of a fixed beam reflectarray antenna at THz based on graphene and compares its performance to a similar implementation using gold for the first time [42]. Soon after, diverse graphene-based reflectarrays operating at THz frequencies have been proposed [43][44][45][46]. As is well known, tunable unit-cell with a full 360 • reflected phase coverage is crucial for realizing a high-performance reflectarray. A small reflected phase tunable range of the unit-cell often leads to deteriorative radiation performance, limiting the function of the whole reflectarray. However, phase tunable ranges of unit-cells in previously reported graphene-based reflectarrays or even metasurfaces [47] are essentially realized by tuning physical parameters of the graphene-based structures. Thus, due to the intrinsic resonant properties and finite losses (such asa damped oscillator) of the graphene-based material, only a narrow band phase tunable range of 300 • has been presented [20,[42][43][44], which restrict their applications in scenarios where high radiation performances are required. Especially, to the best of our knowledge, the CP graphene-based reflectarrays have also been seldom reported in previous literature, although there is a great demand on them at THz communications. One of the reason is that it is hard to achieve an excellent phase tunability by using graphene-based unit-cells. Therefore, designing a unit-cell with a wide-band tunable phase range of 360 • is a meaningful and challenging work for the graphene-based reflectarray design.
In this article, we propose a graphene-based wide-band CP reflectarray operating in the THz regime. The proposed reflectarray consists of a THz CP source and several graphene-based unit-cells. Based on PB phase principle [48], the graphene-based unit-cell is designed to obtain a wide-band (1.4-1.7 THz) tunable phase range of 360 • , which overcomes the restriction of intrinsic narrow band resonant in graphene. Therefore, the graphene-based unit-cell exhibits superior bandwidth and phase tunability to its previous counterparts. Based on the proposed graphene-based unit-cell, a focusing metasurface operating in 1.4-1.7 THz bands is demonstrated firstly, which exhibits prominent focusing effect. Then, according to the reversibility of the optical path, we realize a CP reflectarray by placing a wide-band CP THz source at the focal point of the metasurface. Numerical simulation demonstrates that this reflectarray can achieve a stable high gain up to 15 dBic and an axial ratio around 2.1 dB over the 1.4-1.7 THz bands. The good radiation performance of the proposed CP reflectarray, as demonstrated, underlines its suitability for the THz communication applications.

PB Principle
To have a clear physical insight of the proposed CP reflectarray, firstly, we have deduced the principle of PB phase, obtaining the necessary phase requirement of the unit cell. As shown in Figure 1, we apply a rectangle to represent a PB unit cell. In Figure 1a, a right-handed CP (RHCP) incident wave normally (along the −z direction) illuminates the surface, with an incident electric field expression and the electric field of reflected wave will be where → k denotes the wavenumber in free space. r x and r y represent the reflectivity of x-polarized and y-polarized waves, respectively. Similarly, ϕ x and ϕ y are the reflected phases of x-polarized and y-polarized waves, respectively. To have a clear physical insight of the proposed CP reflectarray, firstly, we have deduced the principle of PB phase, obtaining the necessary phase requirement of the unit cell. As shown in Figure 1, we apply a rectangle to represent a PB unit cell. In Figure 1a, a right-handed CP (RHCP) incident wave normally (along the −z direction) illuminates the surface, with an incident electric field expression inc jkz E ( x jy )e =+ (1) and the electric field of reflected wave will be ( r xe jr ye )e − =+ (2) where k denotes the wavenumber in free space. rx andry represent the reflectivity of x-polarized and y-polarized waves, respectively. Similarly, φx and φy are the reflected phases of x-polarized and y-polarized waves, respectively.  From Figure 1, we can easily obtain that the reflectivity ru and reflected phase φu in the uvz coordinate are equal to the reflectivity rx and reflected phase φx in the xyz coordinate, respectively. In a similar way, the reflectivity rv and reflected phase φv in the uvz coordinate are equal to the reflectivity ry and reflected phase φy in the xyz coordinate, respectively. Therefore, ru = rx, φu = φx, rv = ry, and φv = φy. Then, Equation (5)   As shown in Figure 1b, when the unit cell is anticlockwise rotated with an angle of θ, the relationship between the uv coordinate and the xy coordinate can be represented by Thus, in the uvz coordinate, we substitute Equation (3) into Equation (1), and the incident waves can be expressed as and, corresponding reflected waves are From Figure 1, we can easily obtain that the reflectivity r u and reflected phase ϕ u in the uvz coordinate are equal to the reflectivity r x and reflected phase ϕ x in the xyz coordinate, respectively. In a similar way, the reflectivity r v and reflected phase ϕ v in the uvz coordinate are equal to the reflectivity r y and reflected phase ϕ y in the xyz coordinate, respectively. Therefore, r u = r x , ϕ u = ϕ x , r v = r y , and ϕ v = ϕ y . Then, Equation (5) can be expressed as From Equation (6), if r x = r y , it is found that the reflected waves consist of two components, RHCP and left hand CP (LHCP), respectively.
and when |ϕ x − ϕ y | = π, Equations (7) and (8) change to Obviously, there is only an RHCP component inside the reflected wave, and the phase variation is 2θ, equal to two times of the rotation angle. Contrary to the metallic plate, the reflected waves of the PB unit-cell based metasurface only contain the co-polarized components. It is worth noting that the conditions of |ϕ x − ϕ y | ≈ π and r x ≈ r y must be fulfilled.

Graphene Based PB Unit-Cell
Graphene can strongly interact with electromagnetic waves in THz regime through plasmonic resonance [23,24]. However, for practical applications, wave-graphene interactions have to be further improved. In this paper, a graphene-based PB phase unit-cell is proposed, as shown in Figure 2b. The structure of the unit-cell consists of a top layer rectangular graphene patch and a square grounded quartz glass (SiO 2 ) substrate. When we extend these unit-cells periodically along both the x and y directions, a 2-D graphene-based metasurface can be obtained, as shown in Figure 2a. Incident terahertz waves can excite the plasmonic resonances of the graphene patches, and can be totally reflected by the bottom metallic ground. The top layer graphene-patches work as a partially reflecting mirror, and the bottom metallic ground operates as a fully reflecting mirror, respectively. Because of the rectangular shape of the graphene patch, the phases of the reflected waves can be independently manipulated by changing w or l, as shown in Figure 2d. The reason is that the E x component of the incident wave can only excite the plasmonic resonance in the x direction, and E y component of the incident wave can only excite the plasmonic resonance in the y direction [49]. Based on this characteristic, we can easily design a unit-cell fulfilling the conditions of |ϕ x −ϕ y | ≈ π and r x ≈ r y .
The modeling method of the graphene-based unit-cell is given in the following section. In the terahertz frequencies, the complex surface conductivity of graphene is determined by intraband transition. It can be approximated by a semi-classical Drude model [50] where e is the elementary charge, k B is the Boltzmann's constant, is the reduced Plank's constant, T is temperature, τ is the relaxation time, ω is the radian frequency, and E f is the Fermi energy. When E f is much larger than the thermal energy k B T, the complex conductivity of graphene can be simplified to The electron relaxation time τ is the function of the carrier mobility µ, the Fermi energy E f , and Fermi velocity v f . It can be expressed as τ = E f µ ev f µ will reduce the loss and enhance the efficiency of the device. The carrier mobility µ often has a variation range from~1000 cm 2 /Vs to~230,000 cm 2 /Vs with varied fabricated technologies [20]. In our simulation, according to the experimental results in Ref. [51], it is reasonable to assume that the Fermi energy E f = 0.64 eV (corresponding to electron concentration of n = 3 × 10 13 cm −2 in Ref. [51]), carrier mobility µ = 10,000 cm 2 /Vs, temperature T = 300 K and Fermi velocity v f = 10 6 m/s.  The reflectivity and phases of the proposed unit-cell were full-wave simulated using Ansoft HFSS 2017 software (Ansoft, Pittsburgh, PA, USA). The unit-cell depicted in Figure 2b,c is modeled with a graphene patch on the top layer, then it was deposited on a square grounded SiO2 substrate. The SiO2 substrate has a relative permittivity of εr = 3.75, and a loss tangent tanδ = 0.0184. The parameter w and l are the width and length of the graphene patch in the x and y directions, respectively, and p = 15 µm denotes unit-cell side-length which also equal to a periodicity to form the 2-D metasurface. Besides, t = 10 nm and d = 26 µm are the thickness of the metallic ground and quartz glass (SiO2) substrate, respectively. In our simulation, the master and slaver boundary condition was added to the unit-cell to modeling an infinite array. Meanwhile, Floquet port is placed at z = 6d and utilized to interact with the periodic unit-cell structure. In addition, the reflectivity and phases are obtained by the parametric sweep module of HFSS solver. The center frequency is set to 1.52 THz.
Such a PB unit cell can be looked as an anisotropic scatterer, which has a polarization dependent response. In normally incident case, we optimize the geometric parameters of the unit-cell, and choose the geometric parameters w = 13.39 μm, and l = 3.2 μm. Figure 3a displays the reflectivity and reflected phases of the x-and y-polarized waves, respectively. The reflectivity of xand y-polarized waves are almost equal, with values near 1, indicating high efficiency. A nearly constant 180° reflected phase difference between the x-and y-polarized incident waves can be achieved in a wide frequency band (1.4-1.7 THz), realizing exactly 180° at 1.52 THz, as shown in Figure 3b. The reflectivity and phases of the proposed unit-cell were full-wave simulated using Ansoft HFSS 2017 software (Ansoft, Pittsburgh, PA, USA). The unit-cell depicted in Figure 2b,c is modeled with a graphene patch on the top layer, then it was deposited on a square grounded SiO 2 substrate. The SiO 2 substrate has a relative permittivity of ε r = 3.75, and a loss tangent tanδ = 0.0184. The parameter w and l are the width and length of the graphene patch in the x and y directions, respectively, and p = 15 µm denotes unit-cell side-length which also equal to a periodicity to form the 2-D metasurface. Besides, t = 10 nm and d = 26 µm are the thickness of the metallic ground and quartz glass (SiO 2 ) substrate, respectively. In our simulation, the master and slaver boundary condition was added to the unit-cell to modeling an infinite array. Meanwhile, Floquet port is placed at z = 6d and utilized to interact with the periodic unit-cell structure. In addition, the reflectivity and phases are obtained by the parametric sweep module of HFSS solver. The center frequency is set to 1.52 THz.
Such a PB unit cell can be looked as an anisotropic scatterer, which has a polarization dependent response. In normally incident case, we optimize the geometric parameters of the unit-cell, and choose the geometric parameters w = 13.39 µm, and l = 3.2 µm. Figure 3a displays the reflectivity and reflected phases of the xand y-polarized waves, respectively. The reflectivity of xand y-polarized waves are almost equal, with values near 1, indicating high efficiency. A nearly constant 180 • reflected phase difference between the xand y-polarized incident waves can be achieved in a wide frequency band (1.4-1.7 THz), realizing exactly 180 • at 1.52 THz, as shown in Figure 3b.
According to the principle of PB phase, the proposed graphene-based unit-cell can realize co-polarized conversion, which can convert the illuminated CP wave into a co-polarized CP wave efficiently. Figure 3c gives the co-polarized and cross-polarized conversion ratios when there is a normally RHCP incident wave, verifying the co-polarized conversion. In 1.4-1.7 THz bands, high-efficiency co-polarized transformation has been achieved. Parameters r RR and r LR indicate the co-polarized conversion ratio and cross-polarized conversion ratio, respectively. According to the principle of PB phase, the proposed graphene-based unit-cell can realize co-polarized conversion, which can convert the illuminated CP wave into a co-polarized CP wave efficiently. Figure 3c gives the co-polarized and cross-polarized conversion ratios when there is a normally RHCP incident wave, verifying the co-polarized conversion. In 1.4-1.7 THz bands, high-efficiency co-polarized transformation has been achieved. Parameters rRR and rLR indicate the co-polarized conversion ratio and cross-polarized conversion ratio, respectively.

Graphene Metasurface for Focusing
After verifying the wide-band high efficiency co-polarized converting capability of the proposed PB unit-cell, we can compensate arbitrary phase in a range of 0°-360° by element rotating.
To focus the reflected CP wave of a normally incident plane wave, a phase distribution must be fulfilled, where f represents focal length, λ0 represents wavelength in free space, φ0 represents initial phase at the original point. Following Equation (12), we design a focusing graphene metasurface with 51 × 51 unit-cells at 1.52 THz, and the focal length is set to 190 μm. The phase distribution is illustrated in Figure 4a, and corresponding unit-cell distribution is plotted in Figure  4b. In the numerical simulation, a normally RHCP incident wave is placed along the −z direction, obtaining energy distribution in the xoz plane at 1.4, 1.52, 1.6 and 1.7 THz, respectively, as depicted in Figure 5a-d. From these results, the prominent focusing effect can be observed, verifying our design. Furthermore, as the frequency increases, the focal length increases proportional, in accord with Equation (12). It is worth noting that a full 360° reflected phase range is quite important for generating such focusing metasurface. If narrower phase range is applied, as presented by previous studies [20,42,43], focusing performances will be degraded because of out of phase values in some positions of the metasurface.

Graphene Metasurface for Focusing
After verifying the wide-band high efficiency co-polarized converting capability of the proposed PB unit-cell, we can compensate arbitrary phase in a range of 0 • -360 • by element rotating.
To focus the reflected CP wave of a normally incident plane wave, a phase distribution must be fulfilled, where f represents focal length, λ 0 represents wavelength in free space, ϕ 0 represents initial phase at the original point. Following Equation (12), we design a focusing graphene metasurface with 51 × 51 unit-cells at 1.52 THz, and the focal length is set to 190 µm. The phase distribution is illustrated in Figure 4a, and corresponding unit-cell distribution is plotted in Figure 4b. In the numerical simulation, a normally RHCP incident wave is placed along the −z direction, obtaining energy distribution in the xoz plane at 1.4, 1.52, 1.6 and 1.7 THz, respectively, as depicted in Figure 5a-d. From these results, the prominent focusing effect can be observed, verifying our design. Furthermore, as the frequency increases, the focal length increases proportional, in accord with Equation (12). It is worth noting that a full 360 • reflected phase range is quite important for generating such focusing metasurface. If narrower phase range is applied, as presented by previous studies [20,42,43], focusing performances will be degraded because of out of phase values in some positions of the metasurface.

High Gain CP Graphene-Based Reflectarray
Previously, we have proposed a graphene metasurface which has prominent focusing capability in a wide-band of 1.4-1.7 THz. Then, according to the principle of optical path reversibility [52], we can place a wide-band CP THz source at the focal point of the metasurface. Thus, the incident spherical wave from the THz source can be converted to reflected plane wave with enhanced gain. In our simulation, we choose a THz CP horn antenna as the feeding sourceand optimize its radiation pattern (−10 dB beam width) for effective illumination. Simulated radiation patterns and axial ratios of the proposed feeding source at 1.4, 1.52, 1.6, and 1.7 THz are demonstrated in Figures 6 and 7, respectively. Obviously, a stable −10 dB RHCP beam-width around 120° can be achieved, as shown in the figures. According to the calculated radiation angle, center feeding distance 190 μm is chosen, in keeping with the focal length f = 190 μm of the proposed metasurface.

High Gain CP Graphene-Based Reflectarray
Previously, we have proposed a graphene metasurface which has prominent focusing capability in a wide-band of 1.4-1.7 THz. Then, according to the principle of optical path reversibility [52], we can place a wide-band CP THz source at the focal point of the metasurface. Thus, the incident spherical wave from the THz source can be converted to reflected plane wave with enhanced gain. In our simulation, we choose a THz CP horn antenna as the feeding sourceand optimize its radiation pattern (−10 dB beam width) for effective illumination. Simulated radiation patterns and axial ratios of the proposed feeding source at 1.4, 1.52, 1.6, and 1.7 THz are demonstrated in Figures 6 and 7, respectively. Obviously, a stable −10 dB RHCP beam-width around 120° can be achieved, as shown in the figures. According to the calculated radiation angle, center feeding distance 190 μm is chosen, in keeping with the focal length f = 190 μm of the proposed metasurface.

High Gain CP Graphene-Based Reflectarray
Previously, we have proposed a graphene metasurface which has prominent focusing capability in a wide-band of 1.4-1.7 THz. Then, according to the principle of optical path reversibility [52], we can place a wide-band CP THz source at the focal point of the metasurface. Thus, the incident spherical wave from the THz source can be converted to reflected plane wave with enhanced gain. In our simulation, we choose a THz CP horn antenna as the feeding sourceand optimize its radiation pattern (−10 dB beam width) for effective illumination. Simulated radiation patterns and axial ratios of the proposed feeding source at 1.4, 1.52, 1.6, and 1.7 THz are demonstrated in Figures 6 and 7, respectively. Obviously, a stable −10 dB RHCP beam-width around 120 • can be achieved, as shown in the figures. According to the calculated radiation angle, center feeding distance 190 µm is chosen, in keeping with the focal length f = 190 µm of the proposed metasurface.   It is worth noting that the co-polarization radiation patterns have little asymmetric feature. Due to the center-feeding regime [53,54] of the proposed reflectarray, the CP feeding horn will inevitably interact with the reflected waves, causing distortion of the   It is worth noting that the co-polarization radiation patterns have little asymmetric feature. Due to the center-feeding regime [53,54] of the proposed reflectarray, the CP feeding horn will inevitably interact with the reflected waves, causing distortion of the  It is worth noting that the co-polarization radiation patterns have little asymmetric feature. Due to the center-feeding regime [53,54] of the proposed reflectarray, the CP feeding horn will inevitably interact with the reflected waves, causing distortion of the radiation pattern. Thus, we have to optimize the design parameters to maximally reduce the blockage effect of the feeding structure. Besides, we can alleviate such interaction by off-axial feeding regime. Furthermore, wide-band gain and axial ratio results are displayed in Figure 9, stable gain around 15 dBic and axial ratio around 2.1 dB can be achieved in the whole 1.4-1.7 THz bands, exhibiting both high gain and CP capability, simultaneously.
Materials 2018, 11, x FOR PEER REVIEW 9 of 13 radiation pattern. Thus, we have to optimize the design parameters to maximally reduce the blockage effect of the feeding structure. Besides, we can alleviate such interaction by off-axial feeding regime. Furthermore, wide-band gain and axial ratio results are displayed in Figure 9, stable gain around 15 dBic and axial ratio around 2.1 dB can be achieved in the whole 1.4-1.7 THz bands, exhibiting both high gain and CP capability, simultaneously.   . Simulated RHCP gain and axial ratio of the proposed reflectarray. The red and the blue curves represent the RHCP gain and axial ratio, respectively.

Conclusions
We have systematically investigated the spectral responses of graphene-based PB phase unit-cell and demonstrated a practical implementation of wide-band high gain CP reflectarray based on this configuration in the THz regime. Based on PB phase principle, an arbitrary reflected phase controlling can be realized, extremely extending the phase modulated ability of conventional graphene-based unit-cells that only have a phase controlling range of 300°. Finally, a focusing graphene-based metasurface and a reflectarray were designed, simulated and optimized. Simulation results exhibit excellent performances as theoretical expectations. The proposed reflectarray has a stable gain of 15 dBic and axial ratio of 2 dB over the 1.4-1.7 THz bands, demonstrating high gain and CP characteristics, simultaneously. This reflectarray may have a promising application in future THz communications. Meanwhile, for practical fabrication, we can design the graphene-based unit-cell which is consist of five layers, such as graphene layer, alumina layer, polysilicon layer, quartz glass layer, and ground. The polysilicon can be applied as an electrode. The Fermi energy is related to the conductivity of graphene, and can be dynamically tuned by varying the DC voltage (VDC) between the graphene and the polysilicon. Detailed techniques to fabricate these kinds of graphene reflectarray can be found in Reference [55], supporting the feasibility of our design. In addition, it is worth noting that the general technical procedure formulated herein facilitates further production of such graphene-based devices for various applications.
Author Contributions: L.D. proposed the idea and designed the structure; Y.Z. and J.Z. analyzed the data; and L.D. and C.Z. wrote the paper.

Conclusions
We have systematically investigated the spectral responses of graphene-based PB phase unit-cell and demonstrated a practical implementation of wide-band high gain CP reflectarray based on this configuration in the THz regime. Based on PB phase principle, an arbitrary reflected phase controlling can be realized, extremely extending the phase modulated ability of conventional graphene-based unit-cells that only have a phase controlling range of 300 • . Finally, a focusing graphene-based metasurface and a reflectarray were designed, simulated and optimized. Simulation results exhibit excellent performances as theoretical expectations. The proposed reflectarray has a stable gain of 15 dBic and axial ratio of 2 dB over the 1.4-1.7 THz bands, demonstrating high gain and CP characteristics, simultaneously. This reflectarray may have a promising application in future THz communications. Meanwhile, for practical fabrication, we can design the graphene-based unit-cell which is consist of five layers, such as graphene layer, alumina layer, polysilicon layer, quartz glass layer, and ground. The polysilicon can be applied as an electrode. The Fermi energy is related to the conductivity of graphene, and can be dynamically tuned by varying the DC voltage (VDC) between the graphene and the polysilicon. Detailed techniques to fabricate these kinds of graphene reflectarray can be found in Reference [55], supporting the feasibility of our design. In addition, it is worth noting that the general technical procedure formulated herein facilitates further production of such graphene-based devices for various applications.