Abstract

In this paper, the impact of the intelligent reflecting surface (IRS) on the propagation characteristics is analyzed from the arguments of the physical and electromagnetic characteristics of the channel. Analytical expressions for free-space pathloss are derived based on the terahertz multiple-input multiple-output (MIMO) channel model in different communication scenarios. This calculation method is broadly applicable and is used to characterize pathloss as a function of the configuration of our system model. Numerical results show that setting the phase matrix of IRS based only on the angle of arrival (AoA) and angle of departure (AoD) can significantly reduce the pathloss when the far-field condition is met. However, when the distance between the transceiver and the IRS is the same order of magnitude or less than the side length of the IRS, both direction and distance information are required to configure the IRS. This model provides an insightful consideration in the design of IRS-assisted wireless system.

1. Introduction

To meet the requirement of ultrawideband communication for over 100 billion user equipment in the forthcoming sixth-generation (6G) network, higher frequency communication and a variety of wireless technologies have been proposed [1]. Control of the propagation channel is a promising method to break the channel capacity boundary in the future [2]. Enhancing the complexity of physical layer technology is a commonly employed strategy to achieve higher data rates [3]. As an emerging technology, the IRS changes the channel characteristics during propagation directly, improving the efficiency and bandwidth of the communication system.

The IRS consists of a large number of low-cost passive components that can be artificially designed and programmed to change the phase of the received signal. By changing the phase matrix, beamforming and interference suppression at specific receivers can be realized [4, 5]. Therefore, the IRS technology has generated great interest of academia and industry, the RIS TECH Alliance (RISTA) was formed in April 2022, and the RIS-related discussions are also underway in the 3rd Generation Partnership Project (3GPP) Rel-18 [6].

Many important aspects of IRS-assisted wireless communication systems are investigated, including electromagnetic characteristics [7], channel estimation [8], modeling [9], and joint precoding [10]. However, there is only a limited number of papers that consider the design and analysis of IRS-assisted communication in the terahertz (THz) band [11]. The strong atmosphere molecule absorption attenuation and penetration loss during propagation make the implementation of the IRS necessary.

The realistic and pervasive three-dimensional (3D) system model is crucial for the analysis of IRS [12]. However, the standardized channel model has not been developed presently, and many proposed models are specific to limited communication scenarios.

To the best of the authors’ knowledge, the study of IRS-assisted MIMO communication system models in THz band is still in its infancy [1319]. In [13], the latest research on reconfigurable intelligent surface- (RIS-) assisted wireless networks is summarized and the differences between the principles of RIS technology and similar technologies are analyzed. In [14], a channel model for double-IRS-assisted unmanned aerial vehicle (UAV) communication operated at sub-6G band is proposed. In [15], two power allocation methods that can achieve energy savings are proposed, one using the gradient descent method to calculate the phase coefficient of RIS before power allocation and the other using sequential fractional planning to optimize the phase shift, and it is found that RIS-based power allocation can improve the energy efficiency by three times compared to using conventional amplify-and-forward (AF) relaying. In addition, since decode-and-forward (DF) relaying outperforms AF relaying, the papers [16, 17] further compare the performance of DF relaying with IRS. In [18], although it analyzed in detail the coverage enhancement of cellular networks with the assistance of RIS, it mainly focused on the millimeter wave band and did not further consider the terahertz range. In [19], the IRS-assisted free-space pathloss model for wireless communication was established for several scenarios and actual measurement data were provided for verification, but the paper only studied the single-input single-output (SISO) case at 10 GHz.

Against this background, the main contributions of this paper are listed as follows: (1) consider the MIMO transmission framework to improve the single channel transceiver antenna gain to solve the problem of large attenuation of terahertz signal atmospheric molecular absorption and establish a model for the application of IRS in the MIMO channel and (2) the IRS-assisted wireless communication channel model is introduced from the perspective of electromagnetic theory. A generalized free-space pathloss (FSPL) model is derived, and the specific forms of the pathloss in three practical scenarios are further deduced. The effect of the distance between the transceiver and the IRS, the transmission direction, and the IRS size on the FSPL is elucidated.

The following sections of the paper are organized as follows: Section 2 presents the system model, and the general problem formulation is given. Section 3 presents the analytical expressions for the pathloss of the model in different scenarios. The numerical results are analyzed and discussed in Section 4. Finally, conclusions are given in Section 5.

Notations: the bold letters denote matrices and vectors. The superscript is the Hermitian transpose, the expectation is , and is the complex Gaussian distribution. A diagonal matrix is with in the diagonal. Finally, is the identity matrix.

2. System Model

In the 3D system model we proposed, the base station (BS) is represented by S and the user equipment (UE) is represented by D, as shown in Figure 1. When the line-of-sight (LoS) path is obscured by obstacles such as buildings, we assume that the LoS path does not exist due to the strong directivity of the THz wave.

As shown in Figure 2, the S with M antennas transmit the signal to the D with N antennas. The IRS is located in the - plane in a Cartesian coordinate system, and the geometric center of the IRS is located at the origin of the coordinate system. The IRS consists of phase-shifted units, and the side length of each unit is . The main symbols in the model and their description are shown in Table 1.

In particular, , , and are used to denote the Euclidean distances of BS-IRS, IRS-UE, and BS-UE. , , and stand for the complex CIR of the BS-IRS, IRS-UE, and BS-UE, respectively. is the vector expression of the antenna unit at the source; is the vector expression of the antenna unit at the destination. is the unit direction vector of the distance ; is the unit direction vector of the distance .

The phase shift matrix introduced by the IRS is denoted as , where . The signal is transmitted from the and satisfies . Therefore, the signal received by D can be expressed as where is the transmit power, is the additive white Gaussian noise, and .

3. The Pathloss Model

It is assumed that the polarization of the transceiver is always matched in the transmission process. According to the Friis formula, the received signal power of IRS-assisted wireless communication is as follows:

Proof. The proof is provided in the appendix.

In formula (2), , , and denote the radiation pattern of the IRS, the antenna of S, and the antenna of D, respectively. and denote the angle between the transmitting and receiving antenna unit and the peak of the radiation pattern of the BS. indicates the phase shift generated during the propagation. denotes the terahertz wave absorption loss coefficient, and is obtained from [20]. is relative to the specific communication scenario, and is phase of the unit cell configured by the IRS.

Formula (2) indicates that the transceiver is interchangeable, which has important implications for the uplink and downlink design. According to formula (2), the pathloss can be obtained as where due to the strong directionality of the terahertz antenna, we use the normalized power radiation model expressed as formula (4) [21] to represent the radiation pattern of the transceiver antenna and the reflection units of the IRS.

The directions of peak value in radiation pattern of the transmitting antenna and the receiving antenna are both pointing at the geometric center of the IRS. In the model we proposed, we assume that the phase shift is continuous. Therefore, the phase configuration of the IRS is separated into the following three scenarios.

3.1. General Scenario

In the general scenario, the phase shift . Therefore, . The pathloss can be written as

3.2. Far-Field and Near-Field PiBF Scenarios

In the PiBF scenario, we take the exact value of which relate to , , and . To obtain the minimum pathloss, the phase shift of the unit cell is . Therefore, .

When the distance among the transceiver and the IRS meets the far-field condition which means the and , the distance difference among the transceiver and different unit cells on the IRS can be ignored. relates to and . Therefore, the pathloss can be written as

The distance only meets requirement of near field, including (1) and , (2) and , and (3) and . relates to , , and . The pathloss can be written as

3.3. PiBS Scenario

In this scenario, we set the phase of to minimize the pathloss of the direction of S to the IRS and IRS to the D. is related to and , which is the unit vector of distances and . Therefore, can be calculated by where is the vector representing the position of the center of unit cell on the IRS pointing from the origin, as shown in Figure 3.

The pathloss model for the PiBS scenario can be obtained from formulas (3) and (8): where and . denotes the phase shift processing by the IRS, denotes the total equivalent distance, and denotes the equivalent distance on the IRS due to the different positions of the reflecting unit cells.

Although the pathloss obtained by the PiBS is not optimal, it has practical significance. In the real communication scenarios, it is easier to determine the orientation of the transceiver than the exact propagation distance.

4. Performance Analysis and Simulation Results

In this section, we analyze the performance of the proposed channel model through numerical results. The simulations of the pathloss are conducted based on the deformation expressions of the initial expression (3) in the general scenario (5), far-field PiBF (6), near-field PiBF (7), and PiBS scenario (9). In the numerical simulation, the considered system parameters are shown in Table 2.

It is assumed that the horizontal locations of the S, D, and IRS are fixed as the - plane. Figure 4(a) shows the pathloss variation with IRS size when  m; it can be observed that both the PiBS and the PiBF successfully transmit the signal relative to the severe pathloss of the general scenario, and the pathloss is only related to the radiation pattern of the IRS units. According to the simulation results, the pathloss decreases as the number of units on the IRS increases in the PiBS and PiBF scenarios. However, in the general scenario, the pathloss shows a downward trend with great fluctuation. Figure 4(b) shows the pathloss variation with IRS size when  m. Compared to Figure 4(a), the pathloss changes drastically in the PiBS scenario. In this case, the equivalent distance cannot substitute the actual distance. In addition, the near-field PiBF and far-field PiBF no longer coincide entirely when the distance is short, which is in agreement with our expected results.

Figure 5 shows the variation of the pathloss with the distance and . Simulation is based on setting , , and the size of the IRS . In Figure 5(a), compared with the general scenario, the PiBF scenario has a significant improvement of about 34.9 dB, and for the PiBS is about 32 dB. As a result, although the impact of phase on the received power is reduced, it is not completely eliminated. In Figure 5(b), the simulation results are obtained near the threshold between the near and far field. Effective performance can still be obtained in this case, while Figure 5(c) shows the simulation results of the near-field case. According to the simulation result of the PiBS scenario, simply representing the actual distance by the equivalent distance will result in an error in the calculation when the distance between the transceiver and IRS is close, similar to the size of the IRS. Therefore, an insightful conclusion can be obtained that when the far-field condition is met, the pathloss during transmission can be significantly reduced based on the orientation information. While in the near-field condition, information of the precise distance is essential.

The signal transmitted between different height of the S and D through the IRS in the four scenarios is shown in Figure 6. is the vertical distance between the geometric center of S and the - plane. is the vertical distance between the geometric center of D and the - plane. The sign of and is the same as the -axis. Compared to the traditional situation, the pathloss decreases apparently and the fluctuation reduced efficiently in Figures 6(a)6(d). The maximum loss in the PiBS, far-, and near-field PiBF scenarios are 63%, 59%, and 58% of the general scenario, respectively. We set the minimum value to 160 dB in Figure 6(a) for a more distinct comparison.

In Figure 7, the simulation is based on setting  m and the size of the IRS . As shown in Figure 7(a), the pathloss is relatively low only when and depends on the angle in orientation dramatically. Figures 7(b)7(d) show that after the phase shift of the unit cells on the IRS, the phase dependence and pathloss at the same AoA and AoD as the general scenario are reduced largely. Compared with Figures 7(c) and 7(d), Figure 7(b) reflects the asymmetric changes of AoA and AoD since the peculiarities of the phase processing of the PiBS scenario.

Figure 8 uses the same parameters as Figure 7 and depicts the snapshot at . In Figure 8, the performance when the direction of the destination is unknown for the IRS is presented. For the general scenario, pathloss reaches the minimum only when . However, for the PiBS and PiBF scenarios, there is a significant improvement on the robust of pathloss in the dimension of azimuth.

5. Conclusion

In this paper, the FSPL models in different scenarios for the IRS-assisted THz MIMO wireless communication system based on the physical and electromagnetic properties of the IRS are investigated. According to the initial form of pathloss, the FSPL of general scenario, far- and near-field PiBF scenarios, and PiBS scenario are proposed. Numerical results show that IRS can significantly reduce the pathloss during the propagation in THz band. Moreover, when the far-field condition is met, a similar effect of the PiBS and the PiBF scenario can be achieved. The results shown that the IRS-assisted channel propagation characteristics are related to the IRS size and the direction and distance of the transceivers. This model validates that the communication with IRS can be considered an effective technology to meet the demands of 6G wireless networks.

Appendix

The power received by the unit cell on the IRS that emitted by the antenna unit of S is where . Assuming equal weights for each channel, we distribute the total transmit power equally.

The power reflected by the unit cell on the IRS and received by the antenna unit of D is

The efficiency factor is the ratio of the received power to the reflected power on the unit cell on the IRS and represents the insertion loss required to achieve phase change on the unit cell as well as the finite unit cell transmission efficiency. because the IRS is passive. is the atmospheric molecular absorption loss during THz wave transmission. The transmit power can be obtained by formulas (A.1) and (A.2):

The signal received by the receiving antenna at the destination is

Substituting formulas (A.3) and (A.4) into the receiving power

Formula (2) can be derived by the process above.

Data Availability

The data used to support the findings of this study is confidential.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the Key-Area Research and Development Program of Guangdong Province, China (2019B010157001), and the National Natural Science Foundation of China (no. 61821001).