Design and Optimization on Successive RIS-assisted Multi-hop Wireless Communications ^†^†thanks: This work is supported by the National Natural Science Foundation of China (NSFC) under Grants NO.12141107, and the Interdisciplinary Research Program of HUST (2023JCYJ012).

Rujing Xiong, Jialong Lu, Jianan Zhang, Minggang Liu, Xuehui Dong, Tiebin Mi, Robert Caiming Qiu, R. Xiong et al. are with the School of Electronic Information and Communications, Huazhong University of Science and Technology, Wuhan 430074, China (e-mail: rujing@hust.edu.cn).

Abstract

As an emerging wireless communication technology, reconfigurable intelligent surface (RIS) has become a basic choice for providing signal coverage services in scenarios with dense obstacles or long tunnels through multi-hop configurations. Conventional works of literature mainly focus on alternating optimization or single-beam calculation in RIS phase configuration, which is limited in considering energy efficiency, and often suffers from inaccurate channel state information (CSI), poor convergence, and high computational complexity. This paper addresses the design and optimization challenges for successive RIS-assisted multi-hop systems. Specifically, we establish a general model for multi-hop communication based on the relationship between the input and output electric fields within each RIS. Meanwhile, we derive the half-power beamwidth of the RIS-reflected beams, considering the beam direction. Leveraging these models and derivations, we propose deployment optimization and beam optimization strategies for multi-hop systems, which feature high aperture efficiency and significant gains in signal power. Simulation and prototype experiment results validate the effectiveness and superiority of the proposed systems and methods.

Index Terms:

reconfigurable intelligent surface (RIS), multi-hop, aperture efficiency, half-power beamwidth, deployment optimization, beam optimization, prototype experiment.

I Introduction

The Reconfigurable intelligent surface (RIS) has emerged as a promising technique for wireless communication networks [1, 2, 3]. By dynamically adjusting the reflection phase shifts of numerous passive units, RIS enables flexible wireless channel control and configuration, thereby significantly improving wireless signal transmission rates and reliability. In particular, RIS is a digitally controlled metasurface composed of a large number of cost-effective, well-designed passive reflecting units, each capable of independently manipulating the characteristics of incident electromagnetic waves, including phase, amplitude, and polarization. Due to its great potential, RIS has been widely investigated in the scope of wireless communications [4].

Despite the extensive literature on the design and optimization of various RIS-assisted wireless systems, prior works have primarily focused on enhancing wireless links through single signal reflection using one or multiple RISs [5, 6, 7, 8, 9]. This approach may be inadequate to improve wireless signal quality under some adverse propagation conditions, such as indoor/outdoor environments with dense blockages or obstructions. Additionally, the studies indicated that an RIS should be deployed near the base station (BS) or user equipment (UE) for enhanced communication performance [10, 11, 12]. This requirement further complicates the deployment of RIS.

To address the above challenges, cooperative multi-hop reflections are developed [13, 14, 15, 16, 17]. In [13], the authors studied the multi-reflection between a BS and multiple users, jointly designing the active and passive beamforming at BS and RISs based on deep reinforcement learning. The authors in [14] investigated a multi-route multi-hop cascaded RIS-aided system, with the objective of maximizing the sum rate for multiple users. In [18, 15], the optimal reflection path for multi-RIS beam routing problem is studied. Additionally, the authors in [17] optimized the multi-reflection paths for interference management and throughput performance boosting in multi-hop multi-RIS cooperative communications. [16] explored the integration of multi-hop RISs and open radio access networks, in which the best RIS-assisted path was selected for the overall network sum-rate maximization.

Multi-hop RIS plays a crucial role in facilitating communications within complex terrain and extensive urban environments [18]. Compared to the single-reflection link, the multireflection link in general provides more degrees of freedom to bypass dense and scattered obstacles in complex environments. Moreover, the multireflection links offer higher beamforming gains, which can counteract the product-distance path loss that also increases with the number of RIS reflections [19].

In the current literature on phase configuration within multi-hop communications, the methods can be primarily classified into two categories. The first considers the channels between multiple RISs in a cascaded manner, and achieving the optimization objective through alternating optimizations [20, 21, 22, 23]. These methods have two significant drawbacks. One is that RISs are passive reflectors with a large number of units, rendering the acquisition of channel state information between RISs particularly challenging. Additionally, as the number of RISs increases, the expression to solve optimization problems for each RIS becomes increasingly complex, the iterations are highly prone to converging to local suboptimal solutions, and the computational complexity escalates sharply.

Another feasible method involves considering RIS one by one, in which the phase configuration of each RIS can be determined efficiently according to the channel state information of inner RIS [24, 18, 25], or the relative positions of devices in the communication link [17]. This phase configuration method is frequently employed in optimal routing-path selection within multi-RIS networks.

For RIS-assisted communications, different feed source positions always lead to variations in the field distribution on the RIS surfaces, which in turn cause changes in directional gain. Our previous experimental results have also observed this phenomenon [26]. According to radiation theory in array antennas [27, 28], this variability can be attributed to the aperture efficiency of RIS. As a passive reflector array, the magnitude and uniformity of the electromagnetic wave incident on the RIS’s aperture determine the efficiency and directional power gain it can provide. These characteristics have mostly been ignored in existing studies on RIS phase configuration.

In this paper, we take into account the aperture efficiency and beamwidth of the reflecting beam and derive a closed-form expression for beamwidth. Based on the derivations, we and propose strategies for the deployment and beam optimization in multi-hop RIS-assisted communications. In beam radiation design, we employ our previously proposed MA algorithm to achieve multiple beams, whose computational efficiency, optimality, and robustness have been demonstrated in our earlier work [9]. The proposed multi-hop system features high aperture efficiency and significant power gains.

I-A Contributions

The main contributions of this paper can be summarized as follows:

•

Generalized signal model for multi-hop RIS-assisted Communication. We construct a system model that considers the propagation of signals between the base station (BS) and user equipment (UE) through successive multi-hop RIS configurations. The multiple RIS-reflection channel components among the propagation are categorized into two classes. We characterize the relationship between the input and output electric fields among RISs within each class and integrate them to establish a general multi-hop signal model. With this model, the electric field at each point can be expressed analytically.
•

Derivation of the closed-form expression for the beamwidth of RIS-reflected beams. As a reflectarray antenna, we derive the array factor incorporating the beam direction of RIS reflected. Based on that, we further present the closed-form expression for the half-power beamwidth (HPBW) of the reflecting beams, which closely correlates with the number of units (or aperture size) of RIS and the reflecting beam direction.
•

Optimization strategy on RIS deployment and beam design. Leveraging the derived beamwidth expressions and signal model, we propose optimization methods for RIS deployment and beam design, which eliminates the need for CSI estimation and prioritize high aperture efficiency. The derived results and proposed optimization strategies are validated through simulation and prototype experiments. In the prototype test, we demonstrate that employing optimized multi-beam can increase signal power gain by approximately 10 dB compared to widely used single-beam approaches.

I-B Outline

The remainder of the paper is organized as follows. In Section II, we present the modeling of successive RIS-assisted multi-hop communications. Section III analyzes the radiation of RIS, mainly on the aperture efficiency and array factor. Further, the closed-form expression of beamwidth for RIS-reflected beam is derived. In Section IV, we proposed the deployment and beam optimization strategies for enhanced multi-hop RIS-assisted communication. Section V is dedicated to the performance evaluations of the derived results and proposed methods through numerical and prototype experiments. Finally, this paper is concluded in Section VI.

I-C Notations

Notations in this paper are defined as follows. The imaginary unit is indicated by $j$ . The magnitude and complex components of a complex number are represented by $|\cdot|$ and $\mathcal{\Im}(\cdot)$ , respectively. ${\rm arg}(\cdot)$ stands for the angle function. Unless explicitly specified, lower and upper case bold letters denote vectors and matrices. The conjugate transpose and transpose of $\mathbf{A}$ are represented as $\mathbf{A}^{H}$ and $\mathbf{A}^{T}$ , respectively.

II System Model

In this section, we consider the scenario of the multi-hop RIS-assisted downlink communications, as depicted in Fig. 1. For modeling, we categorize the propagations in multi-hop communications into two types: BS-RIS-UE, where the transmitter is the BS and considered a point source; and RIS-RIS-UE, where the signal source is the reflected RIS. These two types of signal models can be linked through the relationship between the incident field and the reflected field on RISs, thereby we can establish a general model for multi-hop RIS-assisted communication.

II-A BS-RIS-UE

Refer to caption — Figure 1: RIS-assistanted multi-hop wireless communications.

For our analysis, we first concentrate on an arbitrary RIS situated on the $xoy$ -plane, comprising multiple units located at $\mathbf{p}_{n}=[x_{n},\ y_{n},\ z_{n}]^{T}$ , $n=1,\ldots,N$ . The RIS is illuminated by a single incident electromagnetic wave (EM) wave originating from a point source location $(r^{\text{i}},\theta^{\text{i}},\phi^{\text{i}})$ . For clarity, $r^{\text{i}}$ denotes the distance from the source to the anchor of the RIS, such as the geometrical center for a planar array or the left end for a uniform linear array, while $r^{\text{i}}(n)$ represents the distance from the source to the $n$ -th unit at $\mathbf{p}_{n}$ , $(\theta^{\text{i}},\phi^{\text{i}})$ denotes the elevation and azimuth angles, as shown in Fig. 2

In the case of a point source, radiation propagates radially. As the EM wave propagates toward the unit at $\mathbf{p}_{n}$ , the attenuation behavior is characterized by the factor $e^{-j2\pi r^{\text{i}}(n)/\lambda}/r^{\text{i}}(n)$ . Suppose the original electric field from the point source is $E^{\text{i}}(r^{\text{i}},\theta^{\text{i}},\phi^{\text{i}})$ . The incident electric field at $\mathbf{p}_{n}$ could be represented as

E^{\text{i}}(n)=E^{\text{i}}(r^{\text{i}},\theta^{\text{i}},\phi^{\text{i}})e^% {-j2\pi r^{\text{i}}(n)/\lambda}/r^{\text{i}}(n).

(1)

Similarly, if $E^{\text{s}}(n)$ denotes the scattered electric field of the unit at $\mathbf{p}_{n}$ , then the electric field at the observation point is $E^{\text{s}}(n)e^{-j2\pi r^{\text{s}}(n)/\lambda}/r^{\text{s}}(n)$ .

As indicated in our previous study [9], the electric field at the observation point along the propagation path through the $n$ -th unit at $\mathbf{p}_{n}$ can be obtained by combining the above components, as

E^{\text{s}}_{n}(r^{\text{s}}_{n},\theta^{\text{s}},\phi^{\text{s}})=\tau(% \theta^{\text{s}}(n),\phi^{\text{s}}(n);\theta^{\text{i}}(n),\phi^{\text{i}}(n% ))e^{j\omega_{n}}\\ \frac{e^{-j2\pi r^{\text{i}}(n)/\lambda}e^{-j2\pi r^{\text{s}}(n)/\lambda}}{r^% {\text{i}}(n)r^{\text{s}}(n)}E^{\text{i}}(r^{\text{i}},\theta^{\text{i}},\phi^% {\text{i}}).

(2)

Here, $\tau(\theta^{\text{s}},\phi^{\text{s}};\theta^{\text{i}},\phi^{\text{i}})$ is the scattering pattern of an isolated unit, which is dependent on both incident angle $(\theta^{\text{s}},\phi^{\text{s}})$ and scattered angle $(\theta^{\text{s}},\phi^{\text{s}})$ . $\omega_{n}$ represents the phase configuration of the $n$ -th unit within the RIS.

In isotropic scattering, incident EM waves uniformly scatter in all directions over the hemisphere of reflection, irrespective of the angle of incidence or observation. Consequently, the unit’s scattering pattern is simplified to $\tau_{n}(\theta^{\text{s}},\phi^{\text{s}};\theta^{\text{i}},\phi^{\text{i}})=\tau$ . The electric field at the observation point is determined by the superposition of individual fields scattered by $N$ units. As

E^{\text{s}}(r^{\text{s}},\theta^{\text{s}},\phi^{\text{s}})\\ =E^{\text{i}}(r^{\text{i}},\theta^{\text{i}},\phi^{\text{i}})\sum_{n=1}^{N}% \tau e^{j\omega_{n}}\frac{e^{-j2\pi r^{\text{i}}(n)/\lambda}e^{-j2\pi r^{\text% {s}}(n)/\lambda}}{r^{\text{i}}(n)r^{\text{s}}(n)}.

(3)

This expression can be written in matrix forms, as

	$\displaystyle E^{\text{s}}(r^{\text{s}},\theta^{\text{s}},\phi^{\text{s}})$	$\displaystyle=\tau\begin{bmatrix}\frac{e^{-j2\pi r^{\text{s}}(1)/\lambda}}{r^{% \text{s}}(1)}&\cdots&\frac{e^{-j2\pi r^{\text{s}}(N)/\lambda}}{r^{\text{s}}(N)% }\end{bmatrix}$		(4)
		$\displaystyle\begin{bmatrix}e^{j\omega_{1}}&&0\\ &\ddots&\\ 0&&e^{j\omega_{N}}\end{bmatrix}\begin{bmatrix}\frac{e^{-j2\pi r^{\text{i}}(1)/% \lambda}}{r^{\text{i}}(1)}\\ \vdots\\ \frac{e^{-j2\pi r^{\text{i}}(N)/\lambda}}{r^{\text{i}}(N)}\end{bmatrix}E^{% \text{i}}(r^{\text{i}},\theta^{\text{i}},\phi^{\text{i}}).$		(4)

This signal model characterizes the input-output behavior of the RIS, enabling us to calculate the required RIS phase configuration based on the desired electric field. In other words, for a given RIS phase configuration, one only needs to determine the elevation and azimuth angles of the observation point relative to the RIS to calculate the electric field at an observation location.

II-B RIS-RIS-UE

A notable feature in multi-hop RIS setups is that, except for the first RIS, all other RIS nodes are illuminated by the reflected signal from its former RIS. This means that these illuminations are often directional and non-uniform, especially in the near-fields. In this section, we consider the field distribution on RIS to model the signal propagation.

Suppose totaly $K$ RISs are employed, the $(k-1)$ -th RIS (denoted as RIS $\{k-1\}$ ) employs $N$ effective units reflecting the signal beam to illuminate the $k$ -th RIS (denoted as RIS $\{k\}$ , with $M$ units), where $k\geq 2$ . For RIS $\{k\}$ , the incident electric field at the $m$ -th unit (coordinates $\mathbf{p}_{m}=[x_{m},\ y_{m},\ z_{m}]^{T}$ , $m=1,\ldots,M$ ) can be calculate through the scattered field of RIS $\{k-1\}$ at $\mathbf{p}_{m}$ , which can be expressed as

E^{\text{i}}_{k}(m)=E^{\text{s}}_{k-1}(r^{\text{s}}_{k-1}(m),\theta^{\text{s}}% _{k-1},\phi^{\text{s}}_{k-1}),

(5)

where $r^{\text{s}}_{k-1}(m)$ represent the distance from the $m$ -th unit on RIS $\{k\}$ to the reference center point on RIS $\{k-1\}$ .

Furthermore, the electric field observed at a specific point after reflection by the hop-node RIS $\{k\}$ can be written as

		$\displaystyle E_{k}^{\mathrm{s}}(r^{\text{s}}_{k},\theta^{\text{s}}_{k},\phi^{% \text{s}}_{k})$		(6)
		$\displaystyle=\sum_{m=1}^{M}\tau e^{j\omega_{m}}E^{\text{i}}_{k}(m)\frac{e^{-j% 2\pi r^{\text{s}}_{k}(m)/\lambda}}{r^{\text{s}}_{k}(m)}$
		$\displaystyle=\sum_{m=1}^{M}\tau e^{j\omega_{m}}E^{\text{s}}_{k-1}(r^{\text{s}% }_{k-1}(m),\theta^{\text{s}}_{k-1},\phi^{\text{s}}_{k-1})\frac{e^{-j2\pi r^{% \text{s}}_{k}(m)/\lambda}}{r^{\text{s}}_{k}(m)}.$

It should be noted that this expression is also an iterative calculation about $E^{\text{s}}_{k-1}(r^{\text{s}}_{k-1}(m),\theta^{\text{s}}_{k-1},\phi^{\text{s% }}_{k-1})$ , with $2\leq k\leq K$ . The observed electric field at a certain point depends on the spatial scattered field of the all $K$ RISs.

II-C General Model

While $k=1$ , for each $m$ , $E^{\text{s}}_{k-1}(r^{\text{s}}_{k-1}(m),\theta^{\text{s}}_{k-1},\phi^{\text{s% }}_{k-1})$ can be obtained through (3). Thus, by combining equation (3), (5), and (6), we can construct a general cascaded model for signal transmission from the BS to a certain UE, with $K$ multi-hop RISs employed.

For simplicity and clarity, we consider a classic two-hop signal model in analysis to reveal this propagation of wireless signals along multi-hop RISs. As illustrated in Fig, 3, the signal propagates along the BS-RIS $\{1\}$ -RIS $\{2\}$ -UE path, and no direct links between the BS and RIS $\{2\}$ . The number of units on RIS $\{1\}$ and RIS $\{2\}$ is $N$ and $M$ , respectively.

After two hops, the scattered field observed at the UE point can be expressed as

$\displaystyle E^{\text{s}}(p)$		(7)
$\displaystyle=$	$\displaystyle\sum_{m=1}^{M}\tau e^{j\omega_{m}}E^{\text{i}}_{2}(m)\frac{e^{-j2% \pi r^{\text{s}}_{2}(m)/\lambda}}{r^{\text{s}}_{2}(m)}$
$\displaystyle=$	$\displaystyle\sum_{m=1}^{M}\tau e^{j\omega_{m}}E^{\text{s}}(r^{\text{s}}(m),% \theta^{\text{s}},\phi^{\text{s}})\frac{e^{-j2\pi r^{\text{s}}_{2}(m)/\lambda}% }{r^{\text{s}}_{2}(m)}$

For each $1\leq m\leq M$ , $E^{\text{s}}_{k-1}(r^{\text{s}}_{k-1}(m),\theta^{\text{s}}_{k-1},\phi^{\text{s% }}_{k-1})$ can be calculated through (3). Thus, the equation holds

E^{\text{s}}(p)=\sum_{m=1}^{M}\tau e^{j\omega_{m}}\frac{e^{-j2\pi r^{\text{s}}% _{2}(m)/\lambda}}{r^{\text{s}}_{2}(m)}E^{\text{i}}(r^{\text{i}},\theta^{\text{% i}},\phi^{\text{i}})\sum_{n=1}^{N}\tau e^{j\omega_{n}}\\ \frac{e^{-j2\pi r^{\text{i}}(n)/\lambda}e^{-j2\pi r^{\text{s}}(m,n)/\lambda}}{% r^{\text{i}}(n)r^{\text{s}}(m,n)}.

(8)

Here, $r^{\text{s}}(m,n)$ represents the distance between the $m$ -th unit on RIS $\{2\}$ and the $n$ -th unit on RIS $\{1\}$ .

This equation can be written in matrix forms as in (9). We have the canonical linear representation to describe the input/output behaviors of multi-hop RIS-assisted communications. The advantage lies in employing a simple system of linear equations to describe input/output field, making it suitable for analyzing and optimizing the performance of RIS-aided systems.

	$\displaystyle E^{\text{s}}(p)=$	$\displaystyle\tau^{2}\begin{bmatrix}\frac{e^{-j2\pi r^{\text{s}}_{2}(1)/% \lambda}}{r^{\text{s}}_{2}(1)}&\cdots&\frac{e^{-j2\pi r^{\text{s}}_{2}(M)/% \lambda}}{r^{\text{s}}_{2}(M)}\end{bmatrix}\begin{bmatrix}e^{j\omega_{1}}&&0\\ &\ddots&\\ 0&&e^{j\omega_{M}}\end{bmatrix}\begin{bmatrix}\frac{e^{-j2\pi r^{\text{s}}(1,1% )/\lambda}}{r^{\text{s}}(1,1)}&\cdots&\frac{e^{-j2\pi r^{\text{s}}(1,N)/% \lambda}}{r^{\text{s}}(1,N)}\\ \vdots&\ddots&\vdots\\ \frac{e^{-j2\pi r^{\text{s}}(M,1)/\lambda}}{r^{\text{s}}(M,1)}&\cdots&\frac{e^% {-j2\pi r^{\text{s}}(M,N)/\lambda}}{r^{\text{s}}(M,N)}\\ \end{bmatrix}\begin{bmatrix}e^{j\omega_{1}}&&0\\ &\ddots&\\ 0&&e^{j\omega_{N}}\end{bmatrix}$		(9)
		$\displaystyle\begin{bmatrix}\frac{e^{-j2\pi r^{\text{i}}(1)/\lambda}}{r^{\text% {i}}(1)}\\ \vdots\\ \frac{e^{-j2\pi r^{\text{i}}(N)/\lambda}}{r^{\text{i}}(N)}\end{bmatrix}E^{% \text{i}}(r^{\text{i}},\theta^{\text{i}},\phi^{\text{i}}).$		(9)

III Radiation Analysis of the RIS

One of the primary objectives of successive multi-hop RIS is to enhance the received signal power at the UEs. Traditional multi-hop communications have almost overlooked the transmission efficiency between RISs, leading to exponential signal energy loss during the multi-hop process. In this section, we take into account the aperture efficiency and half-power beamwidth (HPBW) analysis for the design and optimization of multi-hop systems, which are demonstrated to be crucial for signal improvement.

III-A Aperture Efficiency

We first consider the aperture efficiency, a parameter that measures the ability of an antenna’s physical aperture to convert available input power into radiated, significantly impacting the overall performance and signal strength of the antenna. In aperture-type antennas such as RISs, the physical aperture area $A_{\text{p}}$ determines the theoretical maximum directivity of the antenna, which is

D_{\text{m}ax}=\frac{4\pi}{\lambda^{2}}A_{\text{p}}.

(10)

The achievable gain of RIS can be computed as

G=\varepsilon_{\mathrm{ap}}D_{\text{m}ax}=\varepsilon_{\mathrm{ap}}\frac{4\pi}% {\lambda^{2}}A_{\text{p}},

(11)

where $\varepsilon_{\text{ap }}$ is the aperture efficiency of the array and $0\leq\varepsilon_{\text{ap }}\leq 1$ , $A_{\text{p}}$ is the physical aperture area of RIS. Since the wavelength and physical aperture area are easily determined, the study of gain reduces to an analysis of aperture efficiency, which can be expressed as a product of sub-efficiencies. [29]:

\varepsilon_{\text{ap }}=e_{\text{r}}\varepsilon_{\text{s }}.

(12)

Here $e_{\text{r}}=$ illumination efficiency, $\varepsilon_{\text{s }}=$ spillover efficiency.

Regarding a RIS comprising $N$ units with $A$ rows and $B$ columns, the radiation efficiency $e_{\text{r}}$ is associated with the electric field distribution that the feed source irradiates onto the reflecting surface (RIS), which can be calculated through the following equation [27, 30].

e_{\text{r}}=\frac{\left|\sum_{a=1}^{A}\sum_{b=1}^{B}\left|E_{ab}\right|\right% |^{2}}{AB\sum_{a=1}^{A}\sum_{b=1}^{B}\left|E_{ab}\right|^{2}},

(13)

It can be observed that $e_{\text{r}}$ is maximized when the electric field $E_{ab}$ at each unit on the RIS is equal. In other words, the aperture field on the RIS is uniformly distributed.

Spillover efficiency $\varepsilon_{\text{s }}$ measures that portion of the feed pattern that is intercepted by the effective aperture of the RIS plate (and redirected through the aperture into the main beam) relative to the total energy[29]. As

\varepsilon_{s}=\frac{\int_{0}^{2\pi}\int_{0}^{\theta_{0}}\left|F_{t}\left(% \theta_{t},\phi_{t}\right)\right|^{2}\sin\theta_{t}d\theta_{t}d\phi_{t}}{\int_% {0}^{2\pi}\int_{0}^{\pi/2}\left|F_{t}\left(\theta_{t},\phi_{t}\right)\right|^{% 2}\sin\theta_{t}d\theta_{t}d\phi_{t}},

(14)

where $\theta_{0}$ is the angular aperture of the reflector (angle between the line from the feed source to the edge of the antenna aperture and the central axis of the antenna aperture), $F_{t}\left(\theta_{t},\phi_{t}\right)$ represents the radiation pattern function of the feed source. The spillover efficiency $\varepsilon_{s}$ is maximized when the numerator $\int_{0}^{2\pi}\int_{0}^{\theta_{0}}\left|F_{t}\right|^{2}\sin\theta_{t}d% \theta_{t}d\phi_{t}$ and denominator $\int_{0}^{2\pi}\int_{0}^{\pi/2}\left|F_{t}\right|^{2}\sin\theta_{t}d\theta_{t}% d\phi_{t}$ are equal, that is when the radiation beam illuminate solely on the RIS board.

Recall that $G$ represents the directional gain achievable by the RIS acting as an aperture antenna. It’s maximum is achieved when the aperture efficiency $\varepsilon_{\mathrm{ap}}$ peaks. Fig.II illustrates the relationship between the distance from the feed source to the RIS and the efficiency and directional gain when a horn antenna is used as the feed source. It can be observed that as the distance increases, the illumination efficiency gradually increases, while the spillover efficiency gradually decreases. Meanwhile, the aperture efficiency, which is the product of these two factors, and the directionally correlated gain exhibit a trend of initially increasing and then decreasing.

Based on the equation(13), (14), and the antenna radiation theory [28, 31], in communication links where RIS serves as a reflective feed source, one crucial approach to maximizing antenna gain $G$ can be designing a precise beam that uniformly covers the RIS surface area. As illustrated in Fig. 5.

III-B The Half-Power Beamwidth

To quantify the signal coverage performance, the half-power beamwidth (HPBW) is one of the most important metrics [28]. HPBW is conventionally defined as the beamwidth within which the signals’s normalized radiated power is higher than -3 dB in a narrowband system. In this paper, we consider that there is a uniform electric field distribution within the 3-dB beamwidth coverage area.

With ideal aperture efficiency, the RIS can be modeled as a uniform amplitude and spacing array. The array factor can be derived by treating the units to be point sources. According to the pattern multiplication rule described in antenna theory [28], for an array composed of identical units, the total field can be synthesized by multiplying the array factor of the isotropic sources by the field of a single unit.

Assuming the position vector as follows,

the transmitter (Tx):

\mathbf{r}_{1}=\left[\cos\phi_{1}\sin\theta_{1},\sin\phi_{1}\sin\theta_{1},% \cos\theta_{1}\right]^{T},

the receiver (Rx):

\mathbf{r}_{2}=\left[\cos\phi_{2}\sin\theta_{2},\sin\phi_{2}\sin\theta_{2},% \cos\theta_{2}\right]^{T},

the observation direction:

\mathbf{r}=[\cos\phi\sin\theta,\sin\phi\sin\theta,\cos\theta]^{T},

the unit in the $a$ -th row and $b$ -th column:

\mathbf{r}_{ab}=[(b-1)d,(a-1)d,0]^{T}.

Denote the phase of the ( $a,b$ )-unit on RIS as $\phi_{ab}$ , as the main beam direction is achieved at ( $\theta_{2},\phi_{2}$ ), the phase configuration is

\phi_{ab}=-k\left(\mathbf{r}_{1}*\mathbf{r}_{ab}^{T}+\mathbf{r}_{2}*\mathbf{r}% _{ab}^{T}\right),

(15)

where $k=2\pi/\lambda$ is the wave number. The array factor of RIS can be given by

\mathrm{AF}(\theta,\phi)=\sum_{a=1}^{A}\sum_{b=1}^{B}e^{jk\left(\mathbf{r}*% \mathbf{r}_{ab}^{T}-\mathbf{r}_{2}*\mathbf{r}_{ab}^{T}\right)},

(16)

which can also be rewritten in (17).

	$\displaystyle\mathrm{AF}(\theta,\phi)$	$\displaystyle=\sum_{a=1}^{A}\sum_{b=1}^{B}e^{jkd\left[((b-1)\cos\phi\sin\theta% +(a-1)\sin\phi\sin\theta)-\left((b-1)\cos\phi_{2}\sin\theta_{2}+(a-1)\sin\phi_% {2}\sin\theta_{2}\right)\right]}$		(17)
		$\displaystyle=\sum_{a=1}^{A}e^{jkd(a-1)\left(\sin\phi\sin\theta-\sin\phi_{2}% \sin\theta_{2}\right)}*\sum_{b=1}^{B}e^{jkd(b-1)\left(\cos\phi\sin\theta-\cos% \phi_{2}\sin\theta_{2}\right)}$		(17)

Denote $\Psi_{1}=kd\left(\sin\phi\sin\theta-\sin\phi_{2}\sin\theta_{2}\right)$ , and $\Psi_{2}=kd\left(\cos\phi\sin\theta-\cos\phi_{2}\sin\theta_{2}\right)$ , it holds

\mathrm{AF}(\theta,\phi)=\frac{\sin\left(\frac{A}{2}\Psi_{1}\right)}{\sin\left% (\frac{1}{2}\Psi_{1}\right)}\frac{\sin\left(\frac{B}{2}\Psi_{2}\right)}{\sin% \left(\frac{1}{2}\Psi_{2}\right)}e^{j\left(\frac{A-1}{2}\Psi_{1}+\frac{B-1}{2}% \Psi_{2}\right)}

(18)

If the reference point is the physical center of the array, the array factor of (18) reduces to

\mathrm{AF}(\theta,\phi)=\frac{\sin\left(\frac{A}{2}\Psi_{1}\right)}{\sin\left% (\frac{1}{2}\Psi_{1}\right)}\frac{\sin\left(\frac{B}{2}\Psi_{2}\right)}{\sin% \left(\frac{1}{2}\Psi_{2}\right)}.

(19)

It can be normalized as

\mathrm{AF}(\theta,\phi)=\frac{1}{A}\frac{1}{B}\frac{\sin\left(\frac{A}{2}\Psi% _{1}\right)}{\sin\left(\frac{1}{2}\Psi_{1}\right)}\frac{\sin\left(\frac{B}{2}% \Psi_{2}\right)}{\sin\left(\frac{1}{2}\Psi_{2}\right)}.

(20)

The above expression can be approximated by

\mathrm{AF}(\theta,\phi)\simeq\left(\frac{\sin\left(\frac{A}{2}\Psi_{1}\right)% }{\frac{A}{2}\Psi_{1}}\frac{\sin\left(\frac{B}{2}\Psi_{2}\right)}{\frac{B}{2}% \Psi_{2}}\right)

(21)

With the array factor, we proceed to investigate the HPBW of the directional main beam. Considering the rotational symmetry of the beam, we take the plane containing the azimuth angle of the beam direction as the reference plane, i.e., $\phi=\phi_{2}$ .

Substitute $\phi=\phi_{2}$ into $\Psi_{1}$ and $\Psi_{2}$ , we hold the equations

	$\displaystyle\Psi_{1}$	$\displaystyle=kd\sin\phi_{2}(\sin\theta-\sin\theta_{2}),$		(22)
	$\displaystyle\Psi_{2}$	$\displaystyle=kd\cos\phi_{2}(\sin\theta-\sin\theta_{2}).$		(22)

Denote $\Psi_{3}=\frac{A}{2}\Psi_{1}$ , and $\Psi_{4}=\frac{B}{2}\Psi_{2}$ , the factor array can be written as

\mathrm{AF}(\theta,\phi)\simeq\frac{\sin\left(\Psi_{3}\right)}{\Psi_{3}}\frac{% \sin\left(\Psi_{4}\right)}{\Psi_{4}}

(23)

For the main beam in the plane at any azimuth angle $\phi_{2}$ , the HPBW depends on the elevation angle $\theta$ . To simplify and clarify, we illustrate this with $\phi_{2}=0$ , where it holds

\mathrm{AF}(\theta,\phi)\simeq\frac{\sin\left(\Psi_{4}\right)}{\Psi_{4}}.

(24)

The HPBW can be obtained through

\mathrm{HP}=|\theta^{*}_{1}-\theta^{*}_{2}|,

(25)

where

	$\displaystyle\theta^{*}_{1}$	$\displaystyle=\sin^{-1}\left(1.391\times\frac{2}{Bkd}+\sin\theta_{2}\right),$		(26)
	$\displaystyle\theta^{*}_{2}$	$\displaystyle=\sin^{-1}\left(-1.391\times\frac{2}{Bkd}+\sin\theta_{2}\right).$		(26)

Theorem 1.

Given a reflecting angle $\theta_{2}$ of RIS, the half-power beamwidth can be obtained through $\mathrm{HP}=|\theta^{*}_{1}-\theta^{*}_{2}|$ , where $\theta^{*}_{1}=\sin^{-1}\left(1.391\times\frac{2}{Bkd}+\sin\theta_{2}\right)$ , and $\theta^{*}_{2}=\sin^{-1}\left(-1.391\times\frac{2}{Bkd}+\sin\theta_{2}\right)$ .

Proof.

According to the AF approximation in (24), the 3-dB point for the array occurs when

\Psi_{4}=\frac{B}{2}kd\left(\sin\theta^{*}-\sin\theta_{2}\right)=\pm 1.391,

we obtain

\theta^{*}=\sin^{-1}\left(\pm 1.391\times\frac{2}{Bkd}+\sin\theta_{2}\right).

Then

\mathrm{HP}=|\theta^{*}_{1}-\theta^{*}_{2}|.

∎

It is noteworthy that when the reflecting elevation angle is relatively large, the AF approximation function (21) tends to underestimate the actual value, resulting in a computed HPBW that is smaller than the actual one.

While the main beam direction is arbitrary, the HPBW can also be obtained by constructing and solving an equation through the 3-dB value point (see Appendix B).

For an $N$ -unit uniformly distributed linear array with spacing $d$ . The array factor expression in (16) can be reduced as

\displaystyle\mathrm{AF}(\theta)=\sum_{n=1}^{N}e^{j(n-1)kd(\sin\theta-\sin% \theta_{2})}

(27)

which can be normalized and approximated as (see Appendix A)

\mathrm{AF}(\theta)\simeq\left[\frac{\sin\left(\frac{N}{2}\Psi_{1}\right)}{% \frac{N}{2}\Psi_{1}}\right]

(28)

where $\Psi_{1}=kd(\sin\theta-\sin\theta_{2})$ .

The 3-dB point and beamwidth can also be obtained through (25) and (26).

The closed-form expressions (25) and (26) characterize the relationship between the main beam’s beamwidth in a reflective array and the number of units, as well as the reflection direction, for the first time. In the subsequent section, we will optimize the deployment and beam design of RIS based on these findings.

IV Optimization Strategy

The width of the beam reflected by the RIS is highly dependent on the location of the receiver. In this section, we optimize RIS-assisted communications from both deployment and beamforming perspectives. By leveraging the aperture efficiency theory and half-power beamwidth derived in the previous section, our objective is to enhance the performance of the multi-hop communications through the homogenization of aperture field distribution.

IV-A Deployment Optimization

In successive multi-hop RIS-assisted communication, signals are typically point-to-point. In such single-input single-output scenarios, the phase configuration of the RIS can be directly and easily computed with the location of the transmitter and receiver [29, 9]. One of the choices to maximize aperture efficiency and directional gains is to optimize the deployment of RISs according to the beamwidth. Specifically, one can select an appropriate installation distance based on a fixed RIS aperture, or determine the required RIS aperture (the number of units) based on a fixed distance. For the sake of clarity, we employ a uniform linear array with $N$ units for analysis.

As shown in Fig. 6, RIS $\{1\}$ reflecting the incident signals to RIS $\{2\}$ . The effective aperture of RIS $\{2\}$ is $L\cos\alpha$ .

IV-A1 RIS position optimization

For a fixed-aperture RIS, selecting an appropriate installation position to optimize the distance $r$ can effectively optimize the aperture field distribution and improve the aperture efficiency. This optimization can be formulated as

\displaystyle\text{(P1)}\ \min_{r}\ |2r\sin(\frac{\theta_{h}}{2})-L\cos\alpha|,

(29)

where $\theta_{h}$ is denoted as the HPBW. Recall that $\theta_{h}=|\theta^{*}_{1}-\theta^{*}_{2}|$ , we have the optimum $r_{\rm opt}$ as

r_{\rm opt}=\frac{L\cos\alpha}{2\sin(\frac{|\theta^{*}_{1}-\theta^{*}_{2}|}{2}% )},

(30)

where $\theta^{*}_{1}=\sin^{-1}\left(1.391\times\frac{2}{Nkd}+\sin\theta^{\text{s}}\right)$ , and $\theta^{*}_{2}=\sin^{-1}\left(-1.391\times\frac{2}{Nkd}+\sin\theta^{\text{s}}\right)$ .

IV-A2 RIS aperture optimization

Based on the relationship between beamwidth, the number of array units, and beam direction, we can also determine the required RIS aperture size (i.e., the number of units needed in a uniform array) for signal relays at a fixed location.

Assuming the size of each unit within the uniform linear array is $\lambda/2$ , the relationship between the number of units $N$ and the RIS aperture $L$ is given by $L=N\cdot\lambda/2$ . We can formulate the optimization problem as

\text{(P2)}\ \min_{N\in\mathbb{Z}^{+}}\ |2r\sin(\frac{\theta_{h}}{2})-L\cos% \alpha|.

(31)

The optimum can be achieved by

N_{\rm opt}={\rm round}(\frac{4r\cdot\sin(\frac{|\theta^{*}_{1}-\theta^{*}_{2}% |}{2})}{\lambda\cdot\cos\alpha}).

(32)

Here, ${\rm round(\cdot)}$ denotes the operator of mapping to the nearest whole number, $\theta^{*}_{1}=\sin^{-1}\left(1.391\times\frac{2}{Nkd}+\sin\theta^{\text{s}}\right)$ , and $\theta^{*}_{2}=\sin^{-1}\left(-1.391\times\frac{2}{Nkd}+\sin\theta^{\text{s}}\right)$ .

Optimizing the position and aperture size of the RIS can effectively focus the beam and ensure uniform illumination of the RIS’s effective aperture, thereby achieving high aperture efficiency and directional gains. We will verify this through simulations in Section V.

IV-B Beam Optimization

For a large-aperture RIS, if the position is fixed, a single beam is often insufficient to cover the entire aperture area (as shown in Fig. 5), especially in the near field. It becomes necessary to design the reflecting beam through optimization methods. In this context, traditional methods often struggle to achieve the desired beam coverage and efficiency due to issues such as side-lobe interference and complex implementation requirements.

It has been demonstrated that the phase optimization methods to achieve more beam numbers may reduce the antenna directivity but has no significant impact on the width of the beam width[27, 9]. To design desired radiation patterns for RISs, we propose a modified method based on the MA algorithm, whose convergence and optimality have been proven in our previous work [9]. For simplicity and clarity, we assume all the RISs have the same number of units.

We adopt the scheme to optimize phase configuration RIS by RIS. Specifically, for RIS $\{1\}$ , we optimize its phase configuration to achieve uniform coverage of the reflected beam across the aperture of RIS $\{2\}$ , aligning the coverage area with that of RIS $\{2\}$ , based on the signal expression (4). For RIS $\{k\}$ ( $k\geq 2$ ), the analysis approach remains consistent, transitioning to the signal model as depicted in (6).

Similarly, initially taking the two-hop RIS scenario as an example, the beam design can be classified into two types. If the width of a single beam is over or exact illumination on RIS $\{2\}$ , as depicted in Fig. 5 and Fig. 5, we set the center of RIS $\{2\}$ as the target beam direction, and directly obtain the required phase configuration for RIS $\{1\}$ by

\omega_{n,{\rm opt}}=2\pi(r^{\text{i}}(n)+r^{\text{s}}(n))/\lambda,\ n=1,% \cdots,N.

(33)

Otherwise, for partial illuminations in Fig. 5, we put forward the following optimization steps for RIS beam design, termed the multi-beam approach.

The optimization begins with determining the sampling points on RIS $\{2\}$ , which will serve as the target directions for RIS $\{1\}$ to generate multiple beams, thereby achieving a uniform illumination on RIS $\{2\}$ . To avoid excessive computational complexity and potential unsolvable cases caused by too many sampling points, as well as overly narrow illuminated areas due to too few sampling points, we determine the points according to the relationship between the beamwidth and RIS aperture. The number of points can be determined through

Z=2\times\lceil\frac{L\cos\alpha}{2r\sin(\frac{\theta_{h}}{2})}\rceil-1,

(34)

where $\lceil\cdot\rceil$ denotes the ceiling function.

Set the $Z$ points on RIS $\{2\}$ as the hypothetic UEs, we explore a multi-beam design strategy with the objective of maximizing the minimum power. Formally, this strategy can be expressed as

\max_{\omega_{1},\dots,\omega_{N}\in[0,2\pi)}\min_{z}\ \left\{P^{\text{s}}(r^{% \text{s}}_{z},\theta^{\text{s}},\phi^{\text{s}}),z=1,\ldots,Z\right\}

(35)

Here $P^{\text{s}}(\cdot)=\lvert E^{\text{s}}(\cdot)\rvert^{2}$ represents the power at a specific location.

Recall the equation in (4), if we define $\mathbf{w}=[e^{j\omega_{1}},\cdots,e^{j\omega_{N}}]^{H}$ , and $\mathbf{h}=[h_{1},\cdots,h_{N}]^{T}$ with $h_{n}=\frac{e^{-j2\pi r^{\text{i}}(n)/\lambda}e^{-j2\pi r^{\text{s}}(n)/% \lambda}}{r^{\text{i}}(n)r^{\text{s}}(n)}$ , the optimization can be reformulated as

\displaystyle\text{(P3)}\ \max_{\omega_{1},\dots,\omega_{N}\in[0,2\pi)}\min_{z% }\ \left\{\mathbf{w}^{H}\mathbf{H}_{z}\mathbf{w},z=1,\ldots,Z\right\},

(36)

where $\mathbf{H}=\mathbf{h}\mathbf{h}^{H}$ .

For RIS $\{k\}$ with $2\leq k\leq K-1$ , this optimization problem can also be formulated, with the distinction that $h_{n}=E^{\text{s}}_{k-1}(r^{\text{s}}_{k-1}(n),\theta^{\text{s}}_{k-1},\phi^{% \text{s}}_{k-1})\frac{e^{-j2\pi r^{\text{s}}_{k}(n)/\lambda}}{r^{\text{s}}_{k}% (n)}$ , as in (6).

Problem (P3) belongs to the class of semi-infinite max-min problems, which can be optimally solved by our previously proposed MA algorithm [9].

Algorithm 1 multi-beam approach for RIS phase configuration

1:procedure OP(optimization on RIS phases for

k\in[1,K]

)

2: if

k\leq K-1

then

3: calculate the samp point

Z

on RIS

\{k+1\}

by (34)

4: solve (36) for

\mathbf{w}_{\rm opt}^{k}

by MA method [9]

5: else

6: calculate

\omega_{n,\rm{opt}}

by (37)

7: calculate

\mathbf{w}_{\rm opt}^{k}

with

\omega_{n,\rm{opt}}

8: end if

9: return

\{\mathbf{w}^{1}_{\rm opt},\mathbf{w}^{2}_{\rm opt}\cdots,\mathbf{w}^{K}_{\rm opt}\}

10:end procedure

After $K$ hops, a single UE receives the final signal. As for the last RIS $\{K\}$ , and the phase configuration can be obtained as

\omega_{n,\rm{opt}}^{k}=\arg\{E^{\text{s}}_{k-1}(r^{\text{s}}_{k-1}(n),\theta^% {\text{s}}_{k-1},\phi^{\text{s}}_{k-1})\}+2\pi r^{\text{s}}_{k}(n)/\lambda

(37)

Thus, the phase of all $K$ RISs are all optimally configured. The multi-beam approach for partial illuminations can be summarized in Algorithm 1

The proposed beam optimization method features high flexibility in the multi-hop transmission of the wireless signal. The multi-beam approach leverages not only the positions of the transmitter-receiver and RISs but also utilizes channel state information if available. The primary distinction between them lies in the elements within matrix $\mathbf{H}$ in optimization problem (P3).

V Simulation and Prototype Experiments

We now investigate the beamforming capabilities of RISs to demonstrate the effectiveness of the derived methods and optimization strategies in multi-hop RIS-assisted communications. Initially, we validate the relationship between beamwidth, beam direction, and the number of units based on the derived HPBW expressions. Subsequently, we analyzed the aperture efficiency under different RIS deployment strategies, confirming that the optimized deployments are effective.

Furthermore, we evaluate the impact of different aperture sizes and distances between RISs on the received signal power by employing the proposed multi-beam method and traditional single-beam method. The simulation and prototype-experiment outcomes highlight the importance of the proposed strategies for enhancing communication quality in multi-hop RIS-assisted communications. In all conducted simulations, the coordinate of the Tx remains fixed at (0, 0, 3).

V-A Beamwidth of the Reflecting Wave

In Fig 7, we compare the radiation patterns derived from array factor theory with those obtained through simulation tests. We employed a linear RIS array with 16 units. The results in Fig. 7 indicate that at a beam direction of 45°, the derived array factor closely matches the actual beam radiation pattern. Additionally, when using the MA algorithm to generate multiple beams, the main beams at $-30^{\circ}$ , $0^{\circ}$ , and $30^{\circ}$ all align with the derived, as shown in Fig. 7.

To further validate the relationship between the derived beamwidth, beam direction, and the number of units, we conduct experiments with different parameter settings and compare the theoretical beamwidth calculations with the actual beamwidths. The outcomes are illustrated in Table I. It can be observed that the HPBW becomes smaller with an increasing number of units. Additionally, the beam direction influences the reflected beamwidth, with beams at smaller elevation angles exhibiting narrower widths. In summary, the derived beamwidth in this work aligns closely with the actual measurements. When the number of units exceeds 16, the discrepancy can be less than $0.02^{\circ}$ .

TABLE I: Half-power beamwidth correspond to uniform linear RISs.

RIS units	Beam direction $\theta$	Theoretical HPBW	Measured HPBW	Gap
8	$30^{\circ}$	$14.73^{\circ}$	$14.84^{\circ}$	$0.11^{\circ}$
8	$45^{\circ}$	$18.25^{\circ}$	$18.38^{\circ}$	$0.13^{\circ}$
8	$60^{\circ}$	$28.56^{\circ}$	$28.85^{\circ}$	$0.29^{\circ}$
16	$30^{\circ}$	$7.33^{\circ}$	$7.2^{\circ}$	$0.13^{\circ}$
16	$45^{\circ}$	$9.01^{\circ}$	$9.02^{\circ}$	$0.01^{\circ}$
16	$60^{\circ}$	$12.97^{\circ}$	$12.99^{\circ}$	$0.02^{\circ}$
32	$30^{\circ}$	$3.66^{\circ}$	$3.67^{\circ}$	$0.01^{\circ}$
32	$45^{\circ}$	$4.49^{\circ}$	$4.49^{\circ}$	$0^{\circ}$
32	$60^{\circ}$	$6.38^{\circ}$	$6.38^{\circ}$	$0^{\circ}$
64	$30^{\circ}$	$1.83^{\circ}$	$1.83^{\circ}$	$0^{\circ}$
64	$45^{\circ}$	$2.24^{\circ}$	$2.24^{\circ}$	$0^{\circ}$
64	$60^{\circ}$	$3.18^{\circ}$	$3.18^{\circ}$	$0^{\circ}$
64	$75^{\circ}$	$6.26^{\circ}$	$6.27^{\circ}$	$0.01^{\circ}$

V-B Aperture Efficiency in Deployment Optimization

Deployment optimization in this work primarily targets the enhancement of aperture efficiency. By adjusting the current RIS’s aperture size (number of units) or its distance from the preceding RIS based on the width of the beam reflected by the preceding RIS, higher aperture efficiency can be achieved in multi-hop setups. For comparison, we defined the ratio of RIS effective aperture to beam illumination (EA-B ratio) in either linear or planar array as

{ABR}=\frac{L\cos\alpha}{2r\sin(\frac{\theta_{h}}{2})}.

(38)

When the EA-B ratio nears $\mathbf{1}$ , the theoretical aperture efficiency and achievable directional gain are maximized. We employ two successive RISs for illustration.

V-B1 RIS position optimization for the problem (P1)

We initially simulate a scenario where the number of RIS units is fixed at $32\times 32$ , maintaining a constant reflection beam angle between RIS $\{1\}$ and RIS $\{2\}$ to ensure an unchanged beamwidth. By adjusting the coordinates of RIS $\{2\}$ , we observe and record its aperture efficiency.

The results in Fig. 8 show that as the distance between the two RISs increases, RIS $\{2\}$ ’s aperture efficiency initially rises and then declines. The specific values in Table. II indicate that when RIS $\{2\}$ is positioned at (9, 0, 9), the aperture efficiency is maximized, with the corresponding EA-B ratio being $1.0013$ . Compared to its position at (19, 0, 19), RIS $\{2\}$ achieves an additional power gain of nearly 2 dB.

V-B2 RIS aperture optimization for the problem (P2)

In addition to the positions, we also evaluate the impact of RIS aperture size on the aperture efficiency. In this test, we fix the positions of RIS $\{1\}$ and RIS $\{2\}$ at (0, 0, 0) and (5, 0, 5), respectively, and measure the changes in aperture efficiency by simultaneously increasing the number of rows and columns (aperture size) of RISs. As shown in Fig. 9, the results indicate that when the number of RIS units is $24\times 24$ , RIS $\{2\}$ achieves a maximum aperture efficiency of $0.6062$ , where the corresponding EA-B ratio equals $1.0132$ .

The two experiments on RIS deployment demonstrate that changing the position and aperture size of RISs significantly impacts the aperture efficiency and directional gain of RISs. This underscores the practical significance of optimizing RIS deployment in real-world scenarios.

TABLE II: Aperture efficiency and Direction Gain of RIS

\{2\}

, while RIS

\{1\}

located at

(0,0,0)

RIS $\{2\}$ coordinate	Inner-RIS Distance (m)	EA-B Ratio	Aperture efficiency	Direction gain (dBi)
(5, 0, 5)	7.07	$1.8023$	0.5097	32.16
(6, 0, 6)	8.49	$1.5019$	0.5409	32.41
(7, 0, 7)	9.90	$1.2874$	0.5598	32.56
(8, 0, 8)	11.31	1.1264	0.5902	32.79
$\mathbf{(9,\ 0,\ 9)}$	$\mathbf{12.73}$	$\mathbf{1.0013}$	$\mathbf{0.6086}$	$\mathbf{32.93}$
(10, 0, 10)	14.14	$0.9011$	0.6019	32.88
(11, 0, 11)	15.56	$0.8192$	0.5842	32.75
(13, 0, 13)	18.38	0.6932	0.5346	32.36
(15, 0, 15)	21.21	0.6008	0.4828	31.92
(17,0,17)	24.04	0.5301	0.4367	31.48
(19,0,19)	26.87	0.4743	0.3959	31.06

V-C Aperture Field Distribution in Beam Optimization

In practical scenarios where the aperture size and position of RISs are fixed and immutable, optimizing the beams remains crucial for achieving higher gains on received power, particularly in near-field environments. Based on the proposed signal model, we compare the effects on received signal power between traditional direct calculation (single-beam) methods [17, 18] and the proposed multi-beam approach in partial illumination scenarios.

As aforementioned, the advantage of the multi-beam strategy lies in its ability to homogenize the aperture field distribution of RISs. To illustrate this, we generate heat maps illustrating the field distributions corresponding to both single-beam and multi-beam methods at two distinct positions of RIS $\{2\}$ . Figs. 10 and 10 demonstrate that using traditional single-beam methods leads to signal energy over-concentration on the central portion of the RIS board, intensifying as the distance between RISs decreases. In contrast, the multi-beam approach effectively alleviates this issue, by concurrently generating beam coverage across various regions of RIS $\{2\}$ using RIS $\{1\}$ , as indicated in Figs. 10 and 10.

V-D Different Aperture Field Distribution Lead to Various Received Power

To further investigate the impact of field distribution variations on received signal power, we conduct experiments comparing the multi-beam and single-beam approaches in two-hop RIS-assisted communication. The experiments consist of two primary sets: the first set observes the effect of beam optimization on received signal power with varying numbers of RIS units, while the second set examines the received signal power at different distances between RISs.

We first position Tx, Rx, and two RISs at (0, 0, 3), (5, 0, 0), (0, 0, 0), and (5, 0, 5) respectively, and varied the number of rows (columns) of the RIS. As illustrated in Fig. 11, when the number of rows (columns) is relatively small, the received signal power corresponding to both the single-beam and multi-beam approaches is comparable. As the number increases, the single-beam shows an advantage within a certain range. However, when the number exceeds a certain threshold (e.g., $32\times 32$ ), the advantage of the multi-beam approach becomes increasingly evident.

These findings can be attributed to the fact that when the number of RIS units is small, the beam’s minimum resolution is insufficient, leading to over-illumination of the RIS aperture (as in Fig. 5). This results in similar aperture field distributions with the multi-beam and single-beam approaches. As the number of units gradually increases, the beam width and aperture size become closer (EA-B ratio approaches 1), the single-beam approach produces more regular and uniform beams, thereby slightly outperforming the multi-beam approach. However, with a further increase in the aperture size of the RIS, the multi-beam approach exhibits advantages in field distribution and demonstrates increasingly superior performance.

In another distance experiment, we maintain the positions of Tx and RIS $\{1\}$ constant and changed the position of Rx to (10, 0, 0). By adjusting the coordinates of RIS $\{2\}$ and its distance from RIS $\{1\}$ , we observed the received signal power. The results in Fig. 12 indicate that when the distance between RISs is relatively short, the multi-beam approach exhibits superior performance regarding received signal power. For instance, at a distance of 1.41 meters between RISs, this gain exceeds 5 dB.

From these two sets of experiments, it can be concluded that the multi-beam approach yields significant performance gains in partial illumination scenarios, validating the effectiveness of our proposed beam optimization strategy. Additionally, the results from Figs. 11 and 12 also indicate that selecting the appropriate number of units and RIS positions is crucial for enhancing the received signal power, as highlighted in the deployment optimization.

V-E Experimental Results with Prototype System

To assess the effectiveness of the proposed beam optimization method in practical settings, we conduct experiments with an operational two-hop RIS-assisted communication system. Within the framework of the present model, the phase disparities among different paths are computed directly based on the positions and configurations of the units. This characteristic facilitates beamforming without the necessity for explicit channel estimation procedures.

The prototype validations are conducted utilizing a two-hop 2-bit RIS-assisted wireless system operating at the central frequency of 3.4 GHz, as illustrated in Fig. 13. Each RIS is composed of $16\times 16$ units, with each unit carrying two PIN diodes, corresponding four different response combination when biased in forward and reverse respectively. The source signal with baseband frequency at 100 kHz is generated and modulated using linear frequency modulation (LFM) by employing the USRP 2954R. The signal is transmitted through a transmitting horn antenna (Tx), reflected by the two-hop RISs, and ultimately received by a receiving horn (Rx).

We initially measure the received signal power with no specific phase configuration. The Tx and Rx antenna is positioned in front of the RIS $\{1\}$ and RIS $\{2\}$ with a distance of $1.50$ m and $1.92$ m, respectively. The distance between the two RIS centers is set at $3.74$ meters, with a reflection angle of $32.5^{\circ}$ , as shown in Fig. 14. All the devices are at the same horizontal height. The observed power spectrum is illustrated in Fig. 16, it shows that the received power level is approximately -84.90 dBm.

For each of the comparative objectives, we perform different methods to obtain the phase configurations on RIS $\{1\}$ and RIS $\{2\}$ based on the proposed signal model. These configurations are then transformed into control codebooks. As shown in Fig.15. Four different colors represent four distinct phase configurations with intervals of $\pi/2$ . The codebooks in Fig. 15 and Fig. 15 are obtained through phase compensation methods (belongs to single-beam method) according to the locations of the equipment [29], while Fig. 15 and 15 are from the proposed Algorithm 1.

The experimental results in Fig.16 indicate that, although we can achieve the optimal phase configuration for the two RIS units in the single-beam scheme, the resulting signal gain is not substantial, only around 10 dB, due to the beam illuminating an area far smaller than the RIS aperture. However, when considering the aperture efficiency of the RIS and employing a multi-beam approach to increase the illuminated area of the RIS, the achieved signal power gain is further enhanced, surpassing 19 dB compared to the unconfigured RIS phase. These outcomes underscore the effectiveness of the proposed beam optimization strategy in multi-hop RIS-assisted communications.

VI Conclusion and Future Directions

This paper focuses on the design and optimization of multi-RIS-assisted communications, particularly in multi-hop scenarios. As a reflective aperture antenna, the aperture efficiency of RIS determines its maximum directional gain, a factor often overlooked in current beamforming schemes. This study examines the field distribution on the RIS aperture, derives a closed-form expression for the half-power beamwidth of beams reflected by RIS, and elucidates the relationship between beamwidth, RIS aperture size (number of units), and beam direction. Drawing from beamwidth and aperture considerations, the paper proposes deployment and beam optimization strategies to enhance the performance of multi-hop RIS-assisted communications without the requirement of any CSI estimation. Simulation and prototype experiments validate the accuracy of the derived results and the effectiveness of the proposed strategies. These optimization frameworks and methods are applicable across a broad spectrum of multi-hop RIS networks, including communications in tunnels and urban environments with dense blockages/obstructions.

Appendix A Proof of Equations

For the equation

\displaystyle\mathrm{AF}(\theta)=\sum_{n=1}^{N}e^{j(n-1)\psi},

Multiplying both sides by $e^{j\psi}$ , it can be written as

(\mathrm{AF}(\theta))e^{j\psi}=e^{j\psi}+e^{j2\psi}+e^{j3\psi}+\cdots+e^{j(N-1% )\psi}+e^{jN\psi}

Subtracting the last from the former reduces to

\mathrm{AF}(\theta)\left(e^{j\psi}-1\right)=\left(-1+e^{jN\psi}\right)

(39)

which can also be written as

	$\displaystyle\mathrm{AF}(\theta)$	$\displaystyle=\left[\frac{e^{jN\psi}-1}{e^{j\psi}-1}\right]=e^{j((N-1)/2]\psi}% \left[\frac{e^{j(N/2)\psi}-e^{-j(N/2)\psi}}{e^{j(1/2)\psi}-e^{-j(1/2)\psi}}\right]$		(40)
		$\displaystyle=e^{j[(N-1)/2]\psi}\left[\frac{\sin\left(\frac{N}{2}\psi\right)}{% \sin\left(\frac{1}{2}\psi\right)}\right]$		(40)

In fact, if the reference point is the physical center of the array, the array factor reduces to

\mathrm{AF}(\theta)=\left[\frac{\sin\left(\frac{N}{2}\psi\right)}{\sin\left(% \frac{1}{2}\psi\right)}\right]

(41)

which can be normalized and approximated as

\mathrm{AF}(\theta)\simeq\left[\frac{\sin\left(\frac{N}{2}\psi\right)}{\frac{N% }{2}\psi}\right]

Appendix B

Plot of function $f(x)=\frac{\sin x}{x}$ , and $f(x,y)=\frac{\sin x}{x}\frac{\sin y}{y}$ .

References

[1] T. J. Cui, M. Q. Qi, X. Wan, J. Zhao, and Q. Cheng, “Coding metamaterials, digital metamaterials and programmable metamaterials,” Light-Sci. Appl., vol. 3, no. 10, pp. e218–e218, Oct. 2014.
[2] E. Basar, M. Di Renzo, J. De Rosny, M. Debbah, M.-S. Alouini, and R. Zhang, “Wireless communications through reconfigurable intelligent surfaces,” IEEE Access, vol. 7, pp. 116 753–116 773, Aug. 2019.
[3] M. Di Renzo et al., “Smart radio environments empowered by reconfigurable intelligent surfaces: How it works, state of research, and the road ahead,” IEEE J. Sel. Areas Commun., vol. 38, no. 11, pp. 2450–2525, Jul. 2020.
[4] Q. Wu et al., “Intelligent surfaces empowered wireless network: Recent advances and the road to 6G,” Proc. IEEE, Jun. 2024, early access, doi: 10.1109/JPROC.2024.3397910.
[5] C. You, B. Zheng, and R. Zhang, “Channel estimation and passive beamforming for intelligent reflecting surface: Discrete phase shift and progressive refinement,” IEEE J. Sel. Areas Commun., vol. 38, no. 11, pp. 2604–2620, Nov. 2020.
[6] C. Pan et al., “Multicell MIMO communications relying on intelligent reflecting surfaces,” IEEE Trans. Wirel. Commun., vol. 19, no. 8, pp. 5218–5233, May 2020.
[7] Z. Yang et al., “Energy-efficient wireless communications with distributed reconfigurable intelligent surfaces,” IEEE Trans. Wirel. Commun., vol. 21, no. 1, pp. 665–679, Jul. 2021.
[8] T. V. Nguyen, D. N. Nguyen, M. D. Renzo, and R. Zhang, “Leveraging secondary reflections and mitigating interference in multi-IRS/RIS aided wireless networks,” IEEE Trans. Wirel. Commun., Jan. 2023.
[9] R. Xiong, K. Yin, T. Mi, J. Lu, K. Wan, and R. C. Qiu, “Fair beam allocations through reconfigurable intelligent surfaces,” Dec. 2023. [Online]. Available: https://arxiv.org/abs/2310.15911
[10] Q. Wu and R. Zhang, “Intelligent reflecting surface enhanced wireless network via joint active and passive beamforming,” IEEE Trans. Wirel. Commun., vol. 18, no. 11, pp. 5394–5409, Aug. 2019.
[11] S. Zhang and R. Zhang, “Intelligent reflecting surface aided multi-user communication: Capacity region and deployment strategy,” IEEE Trans. Commun., vol. 69, no. 9, pp. 5790–5806, Sep. 2021.
[12] Q. Wu, S. Zhang, B. Zheng, C. You, and R. Zhang, “Intelligent reflecting surface-aided wireless communications: A tutorial,” IEEE Trans. Commun., vol. 69, no. 5, pp. 3313–3351, May 2021.
[13] C. Huang et al., “Multi-hop RIS-empowered terahertz communications: A DRL-based hybrid beamforming design,” IEEE J. Sel. Areas Commun., vol. 39, no. 6, pp. 1663–1677, Apr. 2021.
[14] R. Liang, J. Fan, J. Luo, and J. Joung, “A multi-route cascaded IRSs beamforming scheme for mmwave communication systems,” IEEE Trans. Veh. Technol., vol. 71, no. 6, pp. 6813–6818, Jun. 2022.
[15] W. Mei and R. Zhang, “Multi-beam multi-hop routing for intelligent reflecting surfaces aided massive MIMO,” IEEE Trans. Wirel. Commun., vol. 21, no. 3, pp. 1897–1912, Mar. 2022.
[16] M. Abou Houran et al., “Intelligent reflecting surfaces assisted cellular V2X based open RAN communications,” IEEE Trans. Veh. Technol., Jan. 2024, early access, doi: 10.1109/TVT.2024.3356977.
[17] X. Ma, H. Zhang, X. Chen, Y. Fang, and D. Yuan, “Multi-hop multi-RIS wireless communication systems: Multi-reflection path scheduling and beamforming,” IEEE Trans. Wirel. Commun., vol. 23, no. 7, pp. 6778–6792, Jul. 2024.
[18] W. Mei and R. Zhang, “Cooperative beam routing for multi-IRS aided communication,” IEEE Wirel. Commun. Lett., vol. 10, no. 2, pp. 426–430, Feb. 2021.
[19] W. Mei, B. Zheng, C. You, and R. Zhang, “Intelligent reflecting surface-aided wireless networks: From single-reflection to multireflection design and optimization,” Proceedings of the IEEE, vol. 110, no. 9, pp. 1380–1400, Sep. 2022.
[20] Z. Zhang and Z. Zhao, “Rate maximizations for reconfigurable intelligent surface-aided wireless networks: A unified framework via block minorization-maximization,” Dec. 2021. [Online]. Available: https://arxiv.org/abs/2105.02395
[21] D. Shen, Z. Zhang, and L. Dai, “Joint beamforming design for RIS-assisted cell-free network with multi-hop transmissions,” Tsinghua Sci. Technol., vol. 28, no. 6, pp. 1115–1127, Dec. 2023.
[22] R. Liang and J. Fan, “Energy-efficient mmwave IoT communications with multihop IRS-assisted systems,” IEEE Internet Things J., vol. 10, no. 21, pp. 19 344–19 355, Nov. 2023.
[23] H. Niu, Z. Chu, F. Zhou, C. Pan, D. W. K. Ng, and H. X. Nguyen, “Double intelligent reflecting surface-assisted multi-user MIMO mmwave systems with hybrid precoding,” IEEE Trans. Veh. Technol., vol. 71, no. 2, pp. 1575–1587, Feb. 2022.
[24] Y. Han, S. Zhang, L. Duan, and R. Zhang, “Cooperative double-IRS aided communication: Beamforming design and power scaling,” IEEE Wirel. Commun. Lett., vol. 9, no. 8, pp. 1206–1210, Aug. 2020.
[25] W. Mei and R. Zhang, “Cooperative multi-beam routing for multi-irs aided massive MIMO,” in Proc. IEEE Int. Conf. Commun. IEEE, Jun. 2021, pp. 1–6.
[26] R. Xiong et al., “Multi-ris-aided wireless communications in real-world: Prototyping and field trials,” Jul. 2023. [Online]. Available: https://arxiv.org/abs/2303.03287
[27] P. Nayeri, F. Yang, and A. Z. Elsherbeni, Reflectarray antennas: theory, designs, and applications. Hoboken, NJ, USA: Wiley, 2018.
[28] C. A. Balanis, Antenna theory: analysis and design. Hoboken, NJ, USA: Wiley, 2016.
[29] W. L. Stutzman and G. A. Thiele, Antenna theory and design. Hoboken, NJ, USA: Wiley, 2012.
[30] Z. Zhang, LinJin, and X. Shu, Radar antenna technology. Beijing, China: PHEI, 2005.
[31] W. Tang et al., “Wireless communications with reconfigurable intelligent surface: Path loss modeling and experimental measurement,” IEEE Trans. Wirel. Commun., vol. 20, no. 1, pp. 421–439, Jan. 2020.

Design and Optimization on Successive RIS-assisted Multi-hop Wireless Communications ††thanks: This work is supported by the National Natural Science Foundation of China (NSFC) under Grants NO.12141107, and the Interdisciplinary Research Program of HUST (2023JCYJ012).