{'fraction_non_alphanumeric': 0.03168274317962245, 'fraction_numerical': 0.011225444340505144, 'mean_word_length': 4.918497004901071, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Reconfigurable intelligent surfaces: Channel characterization and modeling

Reconfigurable intelligent surfaces (RISs) are two dimensional (2D) metasurfaces which can intelligently manipulate electromagnetic waves by low-cost near passive reflecting elements. RIS is viewed as a potential key technology for the sixth generation (6G) wireless communication systems mainly due to its advantages in tuning wireless signals, thus smartly controlling propagation environments. In this paper, we aim at addressing channel characterization and modeling issues of RIS-assisted wireless communication systems. At first, the concept, principle, and potential applications of RIS are given. An overview of RIS based channel measurements and experiments is presented by classifying frequency bands, scenarios, system configurations, RIS constructions, experiment purposes, and channel observations. Then, RIS based channel characteristics are studied, including reflection and transmission, Doppler effect and multipath fading mitigation, channel reciprocity, channel hardening, rank improvement, far field and near field, etc. RIS based channel modeling works are investigated, including largescale path loss models and small-scale multipath fading models. At last, future research directions related to RIS-assisted channels are also discussed.

objects such as buildings. In this case, an RIS can be employed to provide strong LOS links from Tx to RIS and from RIS to Rx, respectively, as shown in Fig. 1. A RIS can be easily deployed on internal walls of indoor environments, external facades of buildings, glasses of windows, etc. As the distance of the link via RIS is longer than the direct Tx-Rx link, a large number of RIS elements is essential to combat the higher path loss.
Starting from the basic RIS-assisted communication scenario, many new applications are promising and under developments, as illustrated in Fig. 2. The potential application scenarios include: 1) High frequency band communications: All the spectra will be fully explored in the 6G era, especially the high frequency bands. As the frequency band increases, the wireless channel shows higher path loss, which limits the communication distance. On one hand, the molecular absorption in high frequency bands aggravates the path loss. On the other hand, the high frequency signal will easily form an inaccessible shadow behind the obstacles due to the weak diffraction ability and serious penetration loss. RIS can be well combined with high frequency band communications by creating alternative reflecting path and employing massive elements to enable high gains and reconstruct propagation path towards the target user to avoid blockage. Initial works on RIS-assisted millimeter wave (mmWave), terahertz (THz) [24]- [28], and optical [29] wireless communications are promising. Since the optical wireless communication is mainly limited to LOS propagation and easily impaired by obstacles and geometric misalignment losses, RIS can help to confront non-LOS (NLOS) propagation by creating new links to improve the coverage range and reduce the outage probability. It showed that placing RIS closer to Rx can reduce outage probability when the jitter of Rx is more violent, while the middle of Tx and Rx is a better choice for other cases [30]. In addition, multi-branch RIS can further reduce outage probability compared with single-branch RIS [31].
2) Space-air-ground-sea integrated networks: While 5G is mainly for terrestrial communications, 6G will be spaceair-ground-sea integrated networks to enable the vision of global coverage. Satellite, unmanned aerial vehicle (UAV), and maritime communications will be integrated with terrestrial communications. The merits of RIS can be well exploited in the space-air-ground-sea integrated networks to achieve a wide range of signal coverage. As an example, for UAV communications, the trajectory mobility of UAV can create favorable propagation conditions, while the RIS can further improve the system performance by being deployed on the ground to reflect signals from ground base stations to UAV [32].
3) Non-orthogonal multiple access (NOMA): RIS-assisted NOMA communication systems can improve the spectral efficiency and user connectivity by opportunistically exploring the users' different channel conditions [16]. Compared to conventional MIMO-NOMA systems, RIS-NOMA systems can overcome the challenges of the random fluctuation of wireless channels, blockage, and user mobility in an energyefficient way [16]. 4) Simultaneous wireless information and power transfer (SWIPT): Since RIS is a passive device, it can act as a relaylike function to reflect signals mixed with information and energy. The drawback of low energy efficiency at the energy receivers can be solved with the aid of RIS. The wireless power transfer range can be enlarged and the number of required energy beams can be reduced, as shown in [16]. 5) Mobile edge computing (MEC): MEC is essential to offload the data traffic of wireless networks by employing the storage at edge servers to tackle the computation latency and guarantee the privacy and security of big data analysis [12]. Because RIS can control the environment by adjusting the reflection coefficients, a beamforming gain can be achieved. Using this feature, in the edge network, RIS can increase the success rate of offloading of edge devices, thereby improving the performance of the entire network and reducing the endto-end signal transmission delay. In addition, the adjustment efficiency of the RIS will further be improved with the utility of computing power of the edge server when RIS is deployed near it, thereby enhancing system coverage and transmission rate, as well as reducing edge network transmission and processing delays. 6) Physical layer security (PLS): As RIS can realize realtime reconfigurable electromagnetic environment and dynamic programmable wireless channel, it will further enhance the signal quality of intended users and decrease the interference signals from other users and then minimize the uncertainty of the environment to the utmost extent. Therefore, it is a potential key technology for PLS application [12]. 7) Sensing and localization: Radio localization includes three parts, i.e., measurements, a reference system, and inference algorithms [33]. RIS can be utilized to create a strong and consistent multipath in the propagation environment and thus can improve the accuracy of localization. For example, RIS can support localization in harsh factory by circumventing 8) Smart radio environment (SRE): Due to its ability to intelligently control the propagation environments, RIS can be deployed to benefit our daily life, including smart city/home/transportation/factory [9]. In general, these applications are essential parts of SREs.
We remark in each application scenario mentioned above, the specific RIS-assisted channel characteristics need to be considered in channel modeling. For example, in high frequency band communications, large bandwidth will bring frequency non-stationarity. In space-air-ground-sea integrated networks, temporal non-stationarity needs to be considered due to the movement of satellite and UAV, which will cause Doppler effect and multipath fast fading. When the size of RIS is large enough, the near field and spatial non-stationarity effects should also be considered. Current researches mainly concentrate on the theoretical analysis of performance improvements provided by RIS. However, the underlying wireless channel, which is crucial to the performance analysis, is not fully investigated. Many simplified assumptions in the theoretical analysis need to be validated by real RIS based channel measurements and experiments. Meanwhile, only a few of them take the RIS electromagnetic properties and their effects on the channel characteristics into account, such as the hardware impairments and mutual coupling between different elements. Moreover, most of the current works consider only the large-scale path loss, while the small-scale multipath fading is ignored or simplified to independent and identically distributed (i.i.d.) Rayleigh fading, which may lead to an overestimated performance gain. To fill these gaps, we focus on the detailed investigation of RIS based channel measurements and experiments, channel characteristics, and channel models, which are of importance to the design, optimization, and performance evaluation of RIS-assisted wireless communication systems. The novelty and main contributions of this paper are as follows.
1) The available RIS based channel measurements and experiments are comprehensively surveyed, unveiling the new applications and performance gains of RIS-assisted wireless communications. Many new RIS based channel observations have been thoroughly analyzed. 2) The various large-scale path loss models are investigated and compared, pointing out their applicabilities and limitations.
3) The available small-scale multipath fading models are reviewed. The features of each modeling approach are discussed in details. The remainder of this paper is organized as follows. Section II investigates the available RIS based channel measurements and experiments. In Section III, the important channel characteristics brought by RIS are comprehensively surveyed. In Section IV, the state-of-the-art channel modeling methods for RIS based channel are presented, including large-scale path loss models and small-scale multipath fading models. In Section V, future research directions on RIS based channel measurements and modeling are discussed. Conclusions are drawn in Section VI.

II. RIS CHANNEL MEASUREMENTS AND EXPERIMENTS
Recently, many types of RISs have been designed and fabricated for different purposes related to 6G research. Up to now, most of the RIS-assisted channel measurements and experiments in the literatures are conducted in relatively simple indoor environments at sub-6 GHz and mmWave bands. Although the prototypes and demonstrations are not enough, they are essential for the validation, theoretical analysis, and potential insights.
A key question is how to realize the RIS-assisted wireless communication systems. In [34], a RIS-assisted wireless communication system prototype was presented. By changing the biasing voltage, the capacitance value of the varactor diode was adjusted and thus the reflection coefficient of each element was controlled, while the phase response was limited to 255 • . The system was based on national instruments (NI) universal software radio peripheral (USRP) platform. At the Tx side, the USRP generated a single tone incident electromagnetic wave to programmable RIS as the carrier signal. At the Rx side, the USRP down-mixed the received modulated radio frequency (RF) signal and sent the baseband signal to the host computer for synchronization and demodulation processing. Real-time video streaming experiment was conducted at 4 GHz in a typical indoor environment with RIS-Rx distance of 4 m, achieving 2.048 Mbps transmission rate. Experiment results showed that the QPSK constellation diagram after equalization was clear and stable, and the video stream can be transmitted smoothly and clearly even in the absence of channel coding. Then, the authors improved the design of RIS with a full 360 • phase response [35]. Real-time video streaming experiment was conducted at 4.25 GHz in an indoor environment with RIS-Rx distance of 3 m, achieving 6.144 Mbps transmission rate by using 8PSK modulation. Furthermore, the real-time RIS-assisted 2×2 MIMO-16QAM wireless communication experiment was conducted at 4.25 GHz in an indoor environment with RIS-Rx distance of 1.5 m [36]- [39]. The 16QAM constellation diagrams of the two bitstreams recovered by the Rx were clear and dense, achieving 20 Mbps transmission rate. Since the reflection amplitude and phase responses of an element were coupled, a non-linear modulation technique was introduced to enable RIS-based QAM modulation under constant envelope constraint. In addition, the discrete phase shift was shown to be negligible when the quantization step was larger than 8.
As a comparison, two RISs working at 2.3 GHz and 28.5 GHz bands with 16×16 elements were utilized for transmission demonstration on NI USRP platform in [40]. Measured in a compact anechoic chamber, a 21.7 dBi antenna gain was obtained by the RIS at 2.3 GHz, while a 19.1 dBi antenna gain was achieved at 28.5 GHz. A 20 m communication link was formed in an indoor environment, and a variety of modulation modes, including QPSK, 16QAM, 64QAM, etc, were employed, achieving 28.34 Mbps data transmission rate at 2.3 GHz.
In [41], a RIS-assisted wireless communication prototype was presented at 3.6 GHz. The RIS was capable of reflection and transmission. Each element was composed of two mirror symmetric layers. The NI USRP and directional double-ridged horn antennas were employed to transmit and receive signals at both sides. One Rx antenna was deployed at the same side, while the other Rx antenna was deployed on the opposite side in an indoor environment, where the walls were covered with wave absorbing materials to avoid the signals scattered from the environment. The results proved that the proposed beamforming scheme can generate directional beam towards the intended users on both sides of the RIS.
In [42], path loss channel measurements were conducted in an anechoic chamber to validate the proposed path loss channel model. The Keysight E8267D signal generator, N9010A signal analyzer, two horn antennas, and three different RISs were utilized. The distances from Tx to RIS and from Rx to RIS were tuned to compare the measurement results with the proposed path loss model for the near field broadcasting and far field beamforming cases. In addition, new refined path loss models were proposed in [43] and corresponding mmWave channel measurements were conducted in an anechoic chamber. Specifically, two RISs were employed. The first one worked at 27 GHz, where the second one worked at 33 GHz. Two measurement systems were built. The first one was in an anechoic chamber to measure the impact of angle of reception on the power distribution of the signal reflected from the RIS. The second one was in a large empty room to measure the impacts of the distances of Tx-RIS and RIS-Rx. The new path loss models were validated. The authors also compared the received power of the 27 GHz RIS and a metal plate with the same shape and size in the anechoic chamber. The results showed that the RIS with a uniform coding state and a high reflection efficiency had the same reflection characteristics as a metal plate of the same shape and size.
In [44], the aforementioned RISs at 4.25 GHz and 27 GHz were employed to study the channel reciprocity. The RIS at 4.25 GHz was made with varactor diodes and its phase can be continuously manipulated when the voltage of the RIS controller varied from 0 V to 21 V, whereas the RIS at 27 GHz was made with PIN diodes and its phase was only binary programmable. The uplink and downlink received complex signals were recorded in an indoor environment under different distances and angles from the Tx and Rx to the RIS, respectively. The RIS was also configured with identical, gradient, and stripe patterns. Good agreement between the uplink and downlink wireless channels was obtained, which validated the channel reciprocity of RIS-assisted time division duplex (TDD) wireless communication systems.
In [45], [46], a RIS at 5.8 GHz was employed for indoor and outdoor trials based on NI USRP platform. A simple channel reciprocity test was conducted by connecting two horn antennas to the two port R&S vector network analyzer (VNA). Two horn antennas were placed close in front of the RIS to minimize the effects of multipath components (MPCs) other than the one reflected by the RIS. The measured S 12 and S 21 parameters represented the uplink and downlink channels, respectively. Similar to [44], the amplitudes and phases of the two channels showed consistent, which validated the channel reciprocity. In addition, initial indoor and outdoor over the air (OTA) tests were carried out. In the indoor test, about 26 dBi gain was achieved by using RIS. In outdoor tests, about 27 dBi and 14 dBi gains were obtained for 50 m and 500 m links, respectively. The larger distance had performance degradation as the wireless channel became more complicated.
In [47], RIS-assisted mmWave channel measurements were conducted in indoor laboratory and office environments. Different from the conventional case where the RIS is between the Tx and Rx, the RIS was utilized as Tx. It was based on linearly polarized binary-phase unit cells with 180 • phaseshifting capability. The phase shift provided by each unit cell through a proper bias of the PIN diodes was computed using an optimization tool. The channel sounder consists of a four-ports VNA, the RIS, one monopole antenna, and an antenna positioner. Channel measurements were conducted at 28 GHz with 4 GHz bandwidth. The results showed that antenna beamforming can largely reduce the path loss and delay spread in the main beam directions, which led to an improvement of the link budget.
In [48], [49], experiments were conducted to assess the RISassisted ambient backscatter communications. The source, the tag, and the reader were all equipped with dipole antennas. The USRP was used at both sides. The transmitted signal was at 5.4 GHz with 500 MHz bandwidth. The reader received the combined signals of the source, the RIS, and the tag. As the RIS controlled the reflected signal, a hot spot or coherent spots on the tag and the reader can be created to assist the communication between the tag and the reader.
In [50], the concept of active RIS was proposed to break the double fading effect of passive RIS. Different from the existing passive RIS that reflects signals passively without amplification, active RIS can actively amplify the reflected signals. Experiments were conducted by using an active RIS at 2.36 GHz. The incident signal and the pump input were coupled in a varactor-diode-based reflection-type amplifier to generate the reflected signal with amplification.
In [51], a RIS was fabricated at 28 GHz based on a binaryphase tunable metasurface. Experiments were conducted in near field and far field cases. In near field case, the RIS was designed to be a typical reflectarray antenna and proved to be able to beamform mmWave signals very efficiently, achieving directivity around 30 dBi and reaching angles as high as 60 • as predicted by the analytical model. In far field case, two horn antennas were connected to a VNA without LOS propagation. The RIS was proved to be able to create a LOS propagation, which provided a 25 dBi gain over the received power. Based on the design described in [51], a RIS was fabricated at 28 GHz in [52]. Four RISs with unit cells of 20×20 were assembled to achieve 40×40 unit cells. Indoor measurements were conducted to demonstrate the functions of RIS. The RIS was placed on a wall and two horn antennas were used at the Tx and Rx sides to emulate the base station and user equipment, respectively. The Tx-RIS and RIS-Rx distances were 5.5 m and 2 m, respectively. To create a RIS-assisted LOS path, the positions of Tx and Rx were estimated in the beamforming algorithm. A 30 dBi gain was achieved in a wide frequency range of 3 GHz. In addition, NI USRP platform was utilized for QPSK signal transmission, where the four modulated symbols were well distinguished.
In [53], a prototype of the binary-phase metasurface was designed and tested. The simulation showed that the reflection amplitudes had little changes but the reflection phases were opposite at 28 GHz when the PIN diode switches. The functions of steerable sum and difference beams were realized. In addition, the beemsteering performance was analyzed in [54]- [57], which can achieve precise tracking with an angle mismatch of 0.5 • -2.5 • and minimizes the beam training overhead through flexible beam refinement.
In [58], experiments were conducted for RIS-assisted wireless communications in rich scattering environments. Two cases were studied. First, an RIS was deployed to shape the channel impulse response (CIR) to enable higher achievable communication rates. Second, the RIS-tunable propagation environment was leveraged as an analog multiplexer to localize non-cooperative objects using wave fingerprints, even when they were outside the LOS propagation region.
In [59], measurements were conducted to show that RIS can enable attackers to perform complex signal manipulation attacks on the physical layer in real time with low-cost. A laptop was connected to the network to measure the end-to-end connection quality. The proposed environment reconfiguration attack scheme was able to largely reduce the throughput.
In [60], a RIS at 28 GHz was designed to perform as reflectarray and transmitarray on each side. Five functions were achieved, including reflection magnitude, transmission magnitude, reflection phase difference, transmission phase difference, and polarization conversion ratio. The channel reciprocity was also validated by connecting two horn antennas to a VNA and put the RIS between them.
A summary of the above-mentioned RIS channel measurements and experiments is given in Table I, classifying the frequency bands, scenarios, system configurations, RIS constructions, experiment results, and channel observations. As can be seen, initial channel measurements and experiments have unveiled various characteristics of RIS-assisted wireless communication systems, such as the transmission rate, received signal power, path loss, and channel reciprocity. Meanwhile, many new applications are proved to be possible, including RIS-assisted ambient backscatter communications, active RIS, RIS as a transmitarray/reflectarray, RIS for physical layer attack, RIS-assisted rich scattering environment, etc.

III. RIS CHANNEL CHARACTERISTICS
As RIS has the ability to tune the amplitude and/or the phase of incident signals, many new channel characteristics will appear, which are substantially different from conventional wireless channels. Some of the specific channel characteristics have been validated by channel measurements and experiments, while others are theoretically analyzed. Here, we mainly introduce channel characteristics of reflection and transmission, Doppler effect and multipath fading mitigation, channel reciprocity, channel hardening, rank improvement, near field and far field, and half-power beamwidth (HPBW) and aperture efficiency. Other new channel characteristics include: multipath labeling [61], which intended to inject a label on propagation paths through the RISs; OTA equalization [62], which adjusted the phases of RIS elements according to the incoming signals to maximize the magnitude of the first tap and made the inter-symbol interference negligible; electromagnetic interference [63], which had a non-negligible impact on communication performance, especially when the size of the RIS grew large; etc.

A. Reflection and transmission
In traditional wireless communication systems, transmission signals will have a large penetration loss caused by objects, which can be up to several tens of dB, especially for mmWave and THz bands. Thus, the wireless signals are mainly received by the Rx antenna via reflection, scattering, and diffraction.
For RIS-assisted communications, the RIS can be configured in the reflecting mode, the transmitting mode, or simultaneously transmitting and reflecting (STAR) mode [15], [64]- [66]. In the reflecting mode, the users at the transmitter's same half-space are served, while the users at the other half-space are served in the transmission mode. In the simultaneously reflection and transmission mode, the users in the full space are able to be served. Preliminary simulation results showed that STAR-RISs can significantly reduce the base station power consumption compared to conventional reflecting/transmittingonly RISs [65], [66]. Up to now, most of the works on RISassisted communications mainly focus on the reflecting mode, while the other two modes need further investigation.

B. Doppler effect and multipath fading mitigation
When the user equipment is moving, it will cause Doppler effect and fast multipath fading. As the RIS can tune the phase of incident signals, it can align the phases of MPCs arriving from different directions, i.e., the rapid fluctuations in the received signal strength due to the Doppler effect can be effectively reduced by using the real-time tunable RIS [67]. The simulation results proved that the multipath fading effect can be totally eliminated when all reflectors in a propagation environment are coated with RIS, while even a few RIS can significantly reduce the Doppler spread as well as the deep fading in the received signal for general propagation environments with several interacting objects [67]. In [68], the authors showed that the received power can be maximized without adding Doppler spread to the system.
In [69], a RIS-assisted high-mobility vehicular communication system was considered. The RIS is fixed between the roadside BS and high-mobility vehicle. The predictable information of vehicle, such as speed, trajectory, and timetable was utilized to design a sub-optimal phase shift set of RIS, which made the RIS capable of maximizing instantaneous signal-to-noise ratio (SNR), minimizing Doppler spread, and keeping low delay spread simultaneously.
Different from [69], a RIS-assisted high-mobility vehicular communication system was studied in [70], [71], where the RIS was mounted on the top of high-mobility vehicle. An efficient two-stage transmission protocol was proposed to achieve RIS channel estimation and refraction optimization for data transmission by exploiting the quasi-static channel between the high-mobility RIS and user as well as the LoS dominant channel between the BS and RIS. Simulation results showed that the proposed protocol can achieve the full RIS passive beamforming gain and convert the overall BS-user channel from fast to slow fading to enable reliable transmission.

C. Channel reciprocity
Channel reciprocity means that the wireless channel from the Tx side to the Rx side is the same with the wireless channel from the Rx side to the Tx side [72]. This phenomenon comes from the symmetry of Maxwell's equations with respect to time [44]. In general, when the control signals applied to the unit cells remain unchanged, commonly designed and fabricated RISs inherently obey the reciprocity theorem. This is because the materials that constitute RISs, such as metal patches, dielectric layers, and electronic components, often conform to the Rayleigh-Carson reciprocity theorem [42], [44]. The experiments at 4.25 GHz and 27 GHz bands demonstrate channel reciprocity, which show consistent between uplink and downlink channels. In addition, by using active nonreciprocal circuits, performing time-varying controls, and employing nonlinearities and structural asymmetries, the RIS channel reciprocity can be broken. An example of reciprocity breaking by RIS was shown in [17]. The channel reciprocity was also validated in [45] and [60].
An angle-dependent phase shifter model for varactor-based RISs was proposed based on the equivalent circuit in [73]. The simulation results indicated that the angle-reciprocity of varactor-based RIS only held under small incident angles of both forward and reverse incident electromagnetic waves, thus limited the channel reciprocity in RIS-assisted TDD systems.
In [74], a RIS-assisted pilot decontamination scheme was proposed for TDD massive MIMO system by breaking the channel reciprocity. Specifically, the phases of RIS elements were configured differently during the uplink transmission and downlink transmission. In this case, the inter-cell interference due to the RIS-assisted link can be mitigated while the non-RIS link interference still existed. If there was only the RIS-assisted link between the Tx and Rx, by increasing the RIS elements and randomly configuring the phases, the pilot contamination can be mitigated asymptotically. The system performance was then proved to be significantly improved.

D. Channel hardening
MIMO channels can provide spatial diversity to reduce SNR variations. In particular, i.i.d. Rayleigh fading channels give rise to channel hardening, where the SNR variations average out as the number of antennas increase. In [75], the authors provided a new definition of channel hardening in RIS-assisted communications, which stated that the ratio between the system SNR and the square of the number of RIS elements tends to be a constant value when the number of RIS elements tends to infinity. The simulation results revealed that the random SNR variations in the case of random phases remain large since there is no hardening. When the direct path is present, the SNR increases very slowly with the number of RIS elements. Hence, an RIS should be properly configured to benefit from the channel hardening and SNR gain.
In [76], the channel hardening effect was investigated for RIS-assisted massive MIMO system. Once the channel hardening effect is achieved, the effective channel gain can be replaced by its mean value. The base station can replace the instantaneous channel gain by its mean value for signal detection in the uplink. In the downlink, each user can treat the mean value of the effective channel gain as the true one to detect the desired signal and no downlink pilot overhead is required for channel estimation.
In [77], the channel hardening effect was theoretically analyzed and the performance bound was asymptotically derived for the RIS-assisted system. The derived ergodic rates showed that hardware impairments, noise, and interference from estimation errors and the NLOS path became negligible as the number of antennas increases.

E. Rank improvement
In point-to-point MIMO communications, large spatial multiplexing gains can be achieved at the high SNR region, while high SNR mainly appears in LOS scenarios where the channel matrix has low rank and therefore does not support spatial multiplexing [78]. In this case, RIS can enrich the propagation environment by adding MPCs with distinctively different spatial angles, thus a multiplexing gain is achieved even when the direct path has low rank. The simulation results in [78] revealed that the performance greatly depends on the path loss and deployment angles.
The rank deficiency can be analyzed in the form of the eigenvalues of the covariance matrix [76]. The RIS was proved to be able to compensate for the rank deficiency appearing in massive MIMO communications.
To solve the problem of rank deficiency due to the unobstructed LOS in UAV communications, RIS was utilized to enrich the channel in [79]. An RIS placement approach was proposed to improve the spatial multiplexing gains and maximize the average channel capacity in a predefined drone corridor. The RIS-assisted channel was shown to have similar performance at sub-6 GHz and mmWave bands.

F. Near field and far field
In the far field, the total channel gain grows linearly with the number of antenna elements N , while it grows as N 2 for RIS [19], [80], [81]. However, the growth rate eventually tapers off when entering the near field. In general, the Rayleigh distance criterion is applied to distinguish the near field and far field cases. In near field, the spherical wavefront effect needs to be modeled, while the plane wavefront can be assumed in far field. Meanwhile, the distances and geometry relationship for each pair of Tx/Rx and RIS elements will influence the electromagnetic response of RIS-assisted channel directly. An exact near field channel model is needed to evaluate the system performance whenever the distance between the Tx/Rx and RIS is comparable with the dimension of RIS [82].

G. HPBW and aperture efficiency
The HPBW of RIS was investigated in uniform linear array (ULA) and uniform rectangular array (URA) configurations under the far field condition without considering path loss and fading [83]. It showed that RIS's HPBW is equal to antenna array's HPBW at the same reflect/transmit angle with maximal ratio combining applied, while it is greater than or equal to antenna array's HPBW when all weights are equal. Another feature is that RIS's HPBW will change with the incident angles. In addition, the aperture efficiency of RIS was studied in [84].

IV. RIS CHANNEL MODELING A. Large-scale path loss models
As the RIS-assisted channel is the cascade of three subchannels, the large-scale path loss models are affected by the electromagnetic properties of RIS and show great differences with traditional channels. Several path loss models have been proposed for different scenarios and taken different parameters into account. A comprehensive survey of the path loss models are given below.
In [3], the authors illustrated the mechanism of RIS by revisiting the two-ray path loss model, which consisted of a LOS ray and a ground-reflected ray. The RIS with Q elements was assumed to be laid on the ground to assist the communications. The received power was approximated as where P t is the transmitted power, and d tr is the Tx-Rx distance. The received power is proportional to Q 2 and decay with the inverse of (d tr ) 2 . By increasing the RIS elements, the path loss is able to be combated. In [85], the authors proposed a far field path loss model for RIS communications based on physical optics technique. It also explained why the surface consists of many elements that individually act as diffuse scatterers but can jointly beamform the signal in a desired direction with a certain beamwidth. Specifically, the far field path loss was expressed as where sinc(x) = sin(x)/x is the sampling function, a and b are the sizes of the rectangular RIS. G t and G r are the Tx and Rx antenna gains, respectively. d 1 and d 2 are the distances from the Tx and Rx to the RIS, respectively. θ i , θ r , and θ s are the incident angle, desired reflection angle, and observed reflection angle, respectively. When we observe at the desired reflection angle, i.e., θ s = θ r , the path loss model simplifies to The received power is proportional to the square of the RIS area and 1/(d 1 d 2 ) 2 , which was different from the conjecture in [3] that the received power would be proportional to 1/(d 1 + d 2 ) 2 . That conjecture might hold for an infinitely large RIS or in the near field, if the RIS is configured to act as a mirror, but provably not in the far field setup.
In [81], the authors compared conventional massive MIMO with RIS-assisted communications and proved that the RISassisted communications can never provide higher information rate or SNR. Numerical results showed that thousands of reflecting elements are needed as the RIS can not amplify the transmitted signal.
A general RIS-assisted path loss model was proposed in [42], which was expressed as where d x and d y are the size of each element along x axis and y axis, respectively. M and N are the number of elements along x axis and y axis, respectively. G is the gain of one element, F combine n,m denotes the effect of the normalized power radiation pattern of each element on the received signal power, Γ n,m is the reflection coefficient of each element with amplitude of A and phase of φ, and r t n,m and r r n,m denote the distances from the Tx and Rx to each element. As we can see, the received power is related to the Tx/Rx antenna gains, the gain, size, number, and specific electromagnetic properties of elements, and the travel distances. In addition, the general path loss model was simplified to three cases, i.e., far field beamforming (both Tx and Rx are in the far field), near field beamforming (at least one in the near field), and near field broadcasting (at least one in the near field).
In far field beamforming case, the path loss was expressed as where F (θ t , φ t ) and F (θ r , φ r ) are the normalized radiation patterns of each element at the incident direction (θ t ,φ t ) and the reflection direction (θ r ,φ r ).
In near field beamforming case, the distances of Tx and Rx to each element is different and the corresponding path loss was expressed as In the near field broadcasting case, the path loss was approximated as Note that for each case, the phase shift of each element was differently designed. The proposed path loss models were also validated by channel measurements at 4.25 GHz and 10.5 GHz in an anechoic chamber.
In [43], the authors further refined their proposed path loss models in [42] to make it easier to use. The impact of the radiation patterns of the antennas and RIS elements was considered as an angle-dependent loss factor. In the far field beamforming case, the path loss model was refined as The near field beamforming path loss was refined as In addition, the path loss model for one element was proposed as The refined path loss models were then validated by channel measurements at 27 GHz and 33 GHz bands in both anechoic chamber and indoor environments.
In [86], a path loss model for RIS communications was proposed by considering the number of elements, the physical dimensions of the RIS, and the radiation pattern and spacing of constituent elements. It was expressed as where η is the power efficiency of RIS, which is less than 1 for passive RIS. A n and φ n are the amplitude and phase of the nth element, respectively. d 1n and d 2n are the distances from the Tx and Rx to the nth element, respectively. The path loss expression for far field scenario was also derived. The proposed path loss model was then applied to a reflectarraytype RIS, yielding insights into performance as a function of the size of the RIS, proximity of the RIS to the Tx and Rx, and the criteria used to control the elements.
In [87], the authors unified the opposite behavior of RIS as a scatterer and as a mirror based on the free-space Green function. The derived path loss model is similar to [86]. The results showed that depending on its size and distance, the RIS can be observed as a zero, one, or two dimensional object whose radiated power exposes a dependence with the fourth, third, or second power of the distance, respectively.
In [88], a simple path loss model was proposed for RIS communications based on the general scalar theory of diffraction and the Huygens-Fresnel principle. The RIS was modeled as a sheet of electromagnetic material of negligible thickness. The proposed approach can identify the conditions under which an RIS of finite size can be approximated as an anomalous mirror or as a scatterer.
In [89], a physics-consistent analytical free-space path loss model for RIS was proposed based on the vector generalization of Green's theorem. The path loss model can be applied to 2D homogenized metasurfaces. The path loss is formulated in terms of a computable integral that depends on the transmission distances, the polarization of the radio waves, the size of the surface, and the desired surface transformation. Closedform expressions were obtained in far field and near field deployments.
In [90], a path loss model was proposed for RIS-assisted THz communication systems, taking into account the THz band and RIS characteristics. Specifically, the proposed path loss model revealed the relationships between the RIS specifications, such as size, number of RIS elements, element's size, reflection coefficient, radiation pattern, phase shift, and the transmission parameters, such as frequency band, Tx-Rx distance, THz specific parameters, and so on.
In [46], a radar cross section (RCS) based path loss model for RIS communications was proposed. The received power was related to the distances from the Tx/Rx to the RIS, the angles in the Tx-RIS-Rx triangle, and the effective area and the reflection coefficient of each element. In addition, channel measurements were conducted in different scenarios to validate the proposed model in both near field and far field conditions.
In [91], [92], the RIS was partitioned into tiles, where each tile consisted of several unit cells. The far field tile response function g is obtained based on the concept of RCS. Then, the free-space path loss for RIS based link is given as where P L t = ( λ 4πd1 ) 2 and P L r = ( λ 4πd2 ) 2 are the freespace path losses for Tx-RIS and RIS-Rx links, respectively. In addition, the continuous and discrete tile response functions were derived for a given transmission mode.
In [97], a similar path loss model was proposed consistent with the path loss model in [42] by using the concept of RCS, which is more common for scattering problems. The maximum RCS of the unit cell can be approximated as the gain of the unit cell multiplied by its area.
To illustrate the far field path loss, Fig. 3 shows an RISassisted communication setup. The Tx and Rx are on the focal of an ellipse with 2D coordinates of (-c, 0) and (c, 0), respectively. The Tx-Rx distance is d = 2c = 200 m. The long semi-axis of the ellipse is a = 200 m. Thus, b = 173 m and d 1 + d 2 = 400 m. The Rayleigh distance for a 1 m×1.02 m RIS at 10.5 GHz in [42] is about 140 m. We simulate the path loss when the RIS moves from d 1 = 140 m to d 1 = 200 m. The path loss is shown in Fig. 4. As can be seen, the path loss decreases from 73 dB to 69 dB. In case of free space path loss with a distance of d 1 + d 2 = 400 m, the path loss at 10.5 GHz is 105 dB. If the 21 dB antenna gains at both sides are excluded, the multiplied path losses of d 1 and d 2 would be larger than the free space path loss. When the size of the RIS reduces from 102×100 to 51×50, the path loss increases about 10 dB, which indicates that a very large size up to ten thousands of RIS elements are needed to combat the path loss.
In Fig. 5, we compare the condition of RIS with and without phase control capability. The general physical model in [92] and the path loss model in [43] are consistent. Note that the received power is relative to a transmitted power of 0 dBm. The simulated frequency band is 5.4 GHz, M = N = 60, d x = d y = 0.33λ, d 1 = 20 m. The incident angle is 60 • and the reflection angle is 45 • . For the programmable RIS, each element is set to a different phase as the way in [43]. It is clear that the received power is largely improved after the phase control of RIS.

B. Small-scale multipath fading models
In the literatures, some simple end-to-end fading channel models were used to assess the RIS-assisted wireless communication system performance. The Rician fading model was used for RIS-assisted base station to UAV channel and UAV to ground channel in [93] and for RIS-assisted downlink NOMA cellular system in [94]. In [95], the performance of RIS-assisted communication system was analyzed under Nakagami-m fading model. In [96], a generalized fading model was applied to RIS-assisted communication system performance analysis in the presence of phase noise, where the Fox's H model was utilized as a general fading distribution with Rayleigh and Rician as special cases.   A few RIS-assisted small-scale channel models have been proposed recently, including phase-shift model, transmissionmode model, physics based channel model, beamspace channel model, and geometry based stochastic model (GBSM), which are summarized in Table II and thereby overviewed  below. 1) Phase-shift model & transmission-mode model: The phase-shift model is widely used in the literature which means the cascaded three sub-channels, i.e., the direct Tx-Rx channel, Tx-RIS channel, and the RIS-Rx channel. However, the phaseshift model has two drawbacks [91], [92]. The first one is that the ranks of Tx-RIS matrix and RIS-Rx matrix are limited by the number of scatterers, which is small compared to the number of RIS elements. The second one is that the physical properties of RIS, the incident and reflected angles, and polarization are not considered. Thus, a transmission-mode model was proposed in [92] by three key ideas. Firstly, it partitioned the large RIS into tiles and calculated tile response function. Then, the low rank channel matrices were decomposed. At last, different transmission modes were selected for each tile, where each transmission mode effectively corresponded to a given configuration of the phase shifts that the unit cells of the tile apply to an impinging electromagnetic wave.
2) Physics based channel model: In [98], [99], a physics based channel model was proposed for RIS communications. It considers the RIS and the scattering environment as a whole by studying the signal's multipath propagation. The joint channel contains a specular component by the RIS LOS link and a scatter component by the NLOS direct link. Channel statistics were analyzed by considering M RIS MPCs and N MPCs for the scattering environment. The model suggests that the RIS-aided wireless channel can be approximated by a Rician distribution.
In [100], the Tx-RIS channel was modeled by the summation of impinging plane waves. Then, the joint spatialtemporal correlation functions for RIS under both isotropic and non-isotropic scattering were derived. A closed-form four dimensional (4D) sinc function is derived to describe the joint spatial-temporal correlation in the isotropic scattering environment. For the more practical non-isotropic scattering, a realistic 3D angular distribution extracted from mmWave channel measurements in the urban micro scenario was utilized.
In [101]- [103], a spatially-stationary channel model was proposed for holographic MIMO small-scale fading in the far filed case. It modeled the array with a massive number of antennas in a compact space, which is the case of RIS. The small-scale fading was modeled as a zero-mean, spatiallystationary, and correlated Gaussian scalar random field satisfying the Helmholtz equation in the frequency domain, which was equivalent to the scalar wave equation in the time domain.
In [104], an analytical channel model was proposed for RIS-assisted free space optical (FSO) systems based on the Huygens-Fresnel principle. Several effects, such as the nonuniform power distribution of Gaussian beams, the impact of RIS size, the positions of the laser source, the lens, and the phase-shift configuration of the RIS were considered. The simulation results validated the accuracy of the proposed analytical model.
3) Beamspace channel model: In [105], a beamspace channel model was proposed for RIS communications, as shown in Fig. 6. Different from the phase-shift model, the proposed model viewed RIS as a controlled scattering cluster reflecting its own MPCs. The model is based on the extended Saleh Valenzuela (SV) channel model. An antenna segmentation method was proposed to fit RIS in the developed far field model and showed the characteristics of its MPCs and how they can be controlled.

4) GBSM:
The GBSM approach has been widely used in 5G and 6G channel modeling, such as the proposed channel models in [106]- [111]. In [112]- [115], an open-source SimRIS channel simulator was presented for RIS-assisted communications, as shown in Fig. 7. The proposed channel model was applicable to various indoor and outdoor scenarios and various frequency bands, including mmWave band. The model is similar to 3GPP cluster channel model. It incorporated the LOS probabilities between terminals, array responses of RISs and Tx/Rx units, RIS element gains, realistic path loss and shadowing models, and environmental characteristics in different scenarios and frequency bands. However, it is only applied to far field case. In [116], the SimRIS channel simulator was adopted to provide performance insights for design of RIS-assisted smart radio environments. The impact of positioning/placement of RISs, the number of reflecting elements and the tilt/rotation of RISs in both single and multi-RIS environments were investigated. Results revealed that these aspects are crucial to achieving capacity and reliability enhancements in smart radio environments.
In [117], a GBSM was proposed for RIS communications, as shown in Fig. 8, where the three ellipses were used to model Tx-RIS, RIS-Rx, and Tx-Rx channels, respectively. The Tx and Rx were assumed to be equipped with a single antenna, whereas the RIS was a rectangular array. As the proposed model assumed that the Tx, RIS, and Rx were placed in the x-y plane of a Cartesian coordinate system, it is a 2D GBSM. The spatial cross-correlation functions (CCFs) of sub-channels associated with different RIS elements were considered.
In [118], the 2D GBSM proposed in [117] was extended to 3D wideband case by adopting three ellipsoids, as shown in Fig. 9. The Tx and Rx were equipped with multiple antennas. The non-stationary property was considered. Statistical properties, such as temporal autocorrelation function (ACF), spatial CCF, and frequency correlation function (FCF) were analyzed.
In [119], we proposed a 3D non-stationary RIS-assisted MIMO GBSM. The twin cluster concept was applied. The channel model supports movements of Tx, Rx, and clusters. The RIS-assisted channel is divided into three sub-channels, i.e., Tx-RIS channel, RIS-Rx, and Tx-Rx channel. The path loss model in [44] was incorporated in the channel coefficients. The temporal ACF and spatial CCF were analyzed.
We further take the practical discrete RIS phase shifts into account in [120]. The cluster evolution in the space domain and practical discrete phase shifts are considered. The statistical properties, including temporal ACF, spatial CCF, local Doppler spread, and RMS delay spread are investigated. We find that RIS plays a role in separating the whole channel. We take the temporal ACF as an example to compare channel statistical properties with and without RIS, as shown in Fig.  10. It can be seen that the absolute value of temporal ACF employing RIS is much larger than the situation without RIS. In fact, temporal ACF reflects the correlation among different time instances. A larger value indicates that channel is more stable with RIS. The phenomenon proves that RIS can improve the stability of channel. This is because RIS can eliminate the phase differences among LOS components and the strength of the received signal is always strong at different time instances.
In [121], a RIS-assisted channel model was proposed for THz band. The Tx and Rx were assumed to be equipped with a single antenna. The spherical wavefront, i.e., the phase difference between adjacent RIS elements was considered. The path loss model for THz bands was employed, which relied on the plate scattering paradigm and molecular absorption. The simulation results showed that a RIS can significantly improve the energy efficiency of MIMO system in near field and beamfocusing.
5) Keyhole channel model: In the concept of keyhole MIMO channel, the channel was considered as a dyad with one degree of freedom, where each RIS element was a perfect example of keyhole. In [130], the RIS-MIMO channel was modeled as a keyhole MIMO channel. The direct Tx-Rx channel was ignored. The joint Tx-RIS-Rx channel was modeled as the sum of the contributions from each RIS element. Based on that, a channel estimation framework was proposed using single value decomposition (SVD) to separate the cascaded channel links and estimate each link separately. Compared to the phase shift model of RIS, this method considers the amplitude response of RIS element, making the channel model more accurate. With the number of RIS elements increasing, the dimension of Tx-RIS channel coefficient matrix and RIS-Rx channel coefficient matrix will increase linearly, and the Fig. 6. A beamspace channel model for RIS-assisted communication scenario [105].
complexity of this channel model also increases. 6) ML-based channel model: ML-based approach has been a new method in channel modeling [131]. Compared to the traditional channel modeling methods, the most important two advantages are that it can predict channel characteristics by training channel measurement data and can explore the complex nonlinear relationships between channel characteristics, frequency band, scenario, and RIS-assisted system configurations. In [132], the unsupervised machine learning (ML) algorithm, i.e., expectation-maximization (EM) algorithm was applied to RIS-assisted channel modeling. Several aspects were considered, including beamforming, spatial channel correlation, phase-shift errors, arbitrary fading conditions, and coexistence of direct and RIS channels. With approximation of a mixture of two Nakagami-m distributions, the EM algorithm was used to estimate the mixture weight, mean powers, and fading parameters, which can characterize the RIS's composite channels as correlated/independent over conventional/generalized fading paths in the presence of phase errors.

A. RIS channel measurements and parameters estimation
Current experiments are mainly conducted to demonstrate RIS-assisted wireless communication systems and transceiver designs. Most of the experiments are conducted in simple environments with a short distance, which is far away from reality. Only path loss channel measurements were conducted in an anechoic chamber, while the multipath fading channel measurements are missing. The physical properties of RIS will have great effects on the incident signals. However, RIS can only tune part of the wireless channel, i.e., the MPCs which interact with RIS. Other MPCs, such as the LOS component and reflected MPCs by other scatterers in the environments will have no change. Although channel characteristics have been theoretically analyzed, to further validate the theoretical analysis, RIS experiments and channel measurements are indispensable.
In RIS-assisted communication scenario, the MPCs which interact with RIS will be affected by the physical properties of RIS. Other MPCs reflected by other scatterers will be independent with RIS. We designed an algorithm to identify whether a MPC interacts with RIS is carried out by investigating this MPC's responses under the regulation of different RIS transmission modes. Similar to the space-alternating generalized expectation-maximization (SAGE) algorithm, the maximum likelihood estimation and grouped coordinate ascent method are used to estimate MPCs' general parameters, i.e., the amplitude, delay, angle of arrival (AoA), angle of departure (AoD), and Doppler frequency. For MPCs that interact with RIS, an extra and novel maximum likelihood estimation for the incident angle and reflect angle at the RIS is considered in each iteration of the algorithm. This extra estimation is realized based on these MPCs' amplitudes under the regulations of different RIS transmission modes [134].
The root mean square estimation error (RMSEE) for different SNRs and for different number of RIS transmission modes  is shown in Fig. 11. It can be seen that RMSEE reduces with an increase in the number of RIS transmission modes, meaning that more RIS transmission modes lead to higher algorithm's accuracy.

B. A pervasive RIS-assisted channel model
None of the mentioned RIS-assisted channel models are general enough to take all channel characteristics into account and can be applied to various frequency bands and scenarios. Different channel models will lead to underestimates or overestimates of RIS-assisted communication system performance evaluation. This situation makes it difficult to acquire a standardized RIS-assisted channel model which is pervasive and can be an add-on component. A pervasive RIS-assisted channel model can apply unified channel modeling framework Fig. 9. A 3D ellipsoid GBSM for RIS-assisted communication scenario [118].
such as the GBSM method, use the unified expression of CIR, and integrate different channel characteristics brought by RIS. By adjusting channel model parameters, the pervasive RISassisted channel model is able to adapt to various frequency bands, scenarios, and RIS configurations. It is also vital for further performance analysis. The new proposed circuit-based channel model considering antenna size and mutual coupling can be a promising method for RIS-assisted channel modeling [133].

C. AI enabled RIS channel modeling
Artificial intelligence (AI) has now been widely applied to wireless communications and channel modeling, including MPC clustering, scenario classification, and channel statistical properties prediction [2]. However, AI enabled wireless  channel modeling is still not fully explored, especially AI enabled RIS channel modeling. Recent works have shown that generative adversarial network (GAN) is able to learn the end-to-end wireless channels in some simple scenarios [135]- [138]. Different from the traditional channel, RISassisted channel is a cascade of three sub-channels, i.e., Tx-Rx, Tx-RIS, and RIS-Rx channels. The inherent nonideal channel characteristics depend on the specific analysis of the three sub-channels. However, as RIS is passive, it is not able to obtain the solely responses of Tx-RIS and RIS-Rx channels from traditional channel measurements. Currently, there are no open RIS-assisted channel datasets, which would limit the application of AI enabled RIS channel modeling.

D. RIS and electromagnetic information theory
In principle, the current channel models are for propagation channel, while the radio channel model is lacking. The radio channel takes into account the effects of Tx and Rx antennas. A novel circuit-based radio channel model was proposed in [133], which can be applied to RIS-assisted radio channel modeling. Due to its unique characteristics, RIS can be a bridge between electromagnetic theory and information theory, building the foundation of electromagnetic information theory. Here, some challenges are pointed out.
Firstly, the essence of RIS amplitude and phase control is actually realized by the circuit behind it, so the RF link may need more consideration. The certain circuit parameters may be included when obtaining channel characteristics or calculating the system capacity. Secondly, there is mutual coupling between a large number of RIS sub-wavelength units. Therefore, the entire antenna array can be regarded as continuous when compactly arranged (such as holographic MIMO). This requires starting from Maxwell' s electromagnetic theory to determine the assumptions and applications of the unit from discrete to continuous. An accurate and easy-to-handle model that can describe the function of the electromagnetic characteristics of the RIS is needed. Thirdly, analysis and optimization of point-to-point RIS wireless networks is an open task, especially for optimizing the phase of the surface and system capacity under different configurations.

E. Applications of RIS channel research to RIS-assisted communications
The RIS channel research will be essential for RIS-assisted communications, since the performance of RIS-assisted communication systems are largely affected by the underlying channel characteristics, which is summarized in Table III.
Take the channel capacity as an example. The channel capacity is a strict upper limit of the rate at which information can be reliably transmitted on the communication channel. RIS can improve the channel capacity by enhancing the received signal power. The ergodic capacity of RIS-assisted single-input and single-output (SISO) system was analyzed in [139]. In [140] and [141], the channel capacity of RIS-assisted MIMO system in the unicast scenario and multicast scenario were studied, respectively. The alternating optimization algorithm was used to jointly optimize the phase shift matrix of RIS and covariance matrix of the transmitted signal vector. The numerical results illustrated that the channel capacity increased logarithmically with the number of antennas and the square of the number of reflecting elements and it tended to be a constant when both of them went to infinity at a fixed ratio. However, it is inversely proportion to the user number. Degree of freedom (DoF) offered by RIS will be limited by the number of elements of RIS, since it is impractical to create a RIS containing infinite reflecting elements. To solve this gap, authors in [142] proposed irregular RIS, which added a selection matrix in the channel transform function. Current analysis of RIS-assisted communications is mainly based on the simplification of RIS and the wireless channel, thus enabling closed-form expressions. More accurate RIS based channel characteristics and channel models can help to improve the design and performance analysis of RIS-assisted wireless communication systems.

F. Advanced hardware of RISs
In most RISs reported above, their main functions are focused on reconfiguring the beam directions and coverages. Recently, some advanced metasurfaces have been developed, such as digital coding and programmable metasurfaces [22], [143]- [149], which can not only manipulate the beam directions and coverages, but also control the number of beams [22], [143], [144], polarization states [145], and waveforms (e.g. cone beams [146] and orbital-angular-momentum vortex beams [147]- [149]) in real time and programmable way. The introduction of time-domain coding increases an additional degree of freedom to control the frequency spectra of electromagnetic waves [150], [151]. Hence space-time-coding digital metasurfaces can manipulate both spatial waves and their spectra simultaneously and independently [152], [153]. More importantly, the digital coding metasurfaces make a bridge between the digital world and physical world, and hence they can be used to process the digital information directly, resulting in information metasurfaces [154]- [156]. The information metasurfaces have been applied to build new-architecture wireless communication transmitter systems [157]- [163]. Based on the information metasurfaces, new electromagnetic information theory have also been presented by considering both electromagnetic waves and bit-stream information simultaneously [164], [165]. The digital coding and programmable metasurfaces, space-time-coding digital metasurfaces, and information metasurfaces can be regarded as the advanced RISs. In the future, investigating the channel characteristics and channel models of such advanced RISs may produce revolutionary developments of the wireless communications.
VI. CONCLUSIONS In this paper, we have provided a comprehensive survey on state-of-the-art RIS experiments and channel measurements, channel characteristics, and large-scale path loss models and small-scale multipath fading channel models. Although various RIS prototypes have been made, the experiments and especially channel measurements are still in the infancy. RIS based communications bring several merits to wireless channel characteristics, including simultaneous reflection and transmission, Doppler effect and multipath fading mitigation, channel reciprocity, channel hardening, rank improvement, etc. The large-scale path loss models and small-scale multipath fading models have been compared. In addition, we have pointed out future directions which need further investigation, i.e., RIS prototypes, experiments, and channel measurements, pervasive RIS channel models, AI enabled RIS channel modeling, RIS and electromagnetic information theory, applications of RIS channel research to RIS-assisted communications, and advanced hardware of RISs. He has participated in ITU-R WP 3J/3K/3M and SG3 meetings regularly since 2014. He has made considerable contributions to the work of several ITU-R Recommendations. He also participates regularly in ITU-R WP 5D meetings and contributes actively to spectrum sharing studies. In 3GPP, he has contributed to RAN1 channel model standards in 3GPP TR38901, relating to V2V, NTN, IIOT, etc.
His current interests include research into 3GPP/ITU radio propagation characteristics at frequencies that may be considered for systems beyond 5G/6G, fixed links and satellite communication systems, the use of reconfigurable intelligent surface (RIS), and the modelling of radio propagation in spectrum sharing studies.
Zhimeng Zhong received the B.E., M.S., and Ph.D. degrees from Xi'an Jiaotong University, China, in 2002China, in , 2005China, in , and 2008, respectively, all in Electronic Engineering. He joined Huawei Technologies CO., LTD from 2009, and about 11 years experiences on wireless communication system research and development. Dr. Zhong is a senior member of the IEEE. In these years, he mainly focused on wireless channel measurement, channel modeling, and MIMO algorithm & system design. He is the team leader of wireless channel research in Huawei. Especially, he was in charge of the 3GPP 5G channel model standardization work in Huawei, and gave contributions on 3GPP mmWave channel model, V2V channel model, and IIOT channel model standardization. Recently, he also focused on massive MIMO algorithm design based on some special channel characteristics and machine learning technologies.