6G

IEEE websites place cookies on your device to give you the best user experience. By using our websites, you agree to the placement of these cookies. To learn more, read our Privacy Policy. Accept & Close Skip to Main Content INSTITUTIONAL SUBSCRIBERS: Are you having difficulty accessing IEEE Xplore when working remotely? Contact IEEE Now IEEE.org IEEE Xplore Digital Library IEEE-SA IEEE Spectrum More Sites Cart Create Account Personal Sign In Institutional Sign In Browse My Settings Get Help Subscribe Enter keywords or phrases (Note: Searches metadata only by default. A search for 'smart grid' = 'smart AND grid') Advanced Search Other Search Options Journals & Magazines > IEEE Access > Volume: 7 Wireless Communications and Applications Above 100 GHz: Opportunities and Challenges for 6G and Beyond Publisher: IEEE Cite This PDF 8 Author(s) Theodore S. Rappaport ; Yunchou Xing ; Ojas Kanhere ; Shihao Ju ; Arjuna Madanayake ; Soumyajit Mandal ; Ahmed Alkhateeb ; Georgios C. Trichopoulos View All Authors 35 Paper Citations 17223 Full Text Views Open Access Comment(s) Abstract Document Sections I. Introduction II. mmWave and THz Applications III. Multi-Port Array Transceiver Architectures Based on Space-Time Sigma-Delta Theory IV. Novel Wideband Multi-Beamformers at Low Size, Weight, and Power (SWaP) V. Moving to Beyond 5G and Above 100 GHz – A New Type of Channel Authors Figures References Citations Keywords Metrics Abstract: Frequencies from 100 GHz to 3 THz are promising bands for the next generation of wireless communication systems because of the wide swaths of unused and unexplored spectrum. These frequencies also offer the potential for revolutionary applications that will be made possible by new thinking, and advances in devices, circuits, software, signal processing, and systems. This paper describes many of the technical challenges and opportunities for wireless communication and sensing applications above 100 GHz, and presents a number of promising discoveries, novel approaches, and recent results that will aid in the development and implementation of the sixth generation (6G) of wireless networks, and beyond. This paper shows recent regulatory and standard body rulings that are anticipating wireless products and services above 100 GHz and illustrates the viability of wireless cognition, hyper-accurate position location, sensing, and imaging. This paper also presents approaches and results that show how long distance mobile communications will be supported to above 800 GHz since the antenna gains are able to overcome air-induced attenuation, and present methods that reduce the computational complexity and simplify the signal processing used in adaptive antenna arrays, by exploiting the Special Theory of Relativity to create a cone of silence in over-sampled antenna arrays that improve performance for digital phased array antennas. Also, new results that give insights into power efficient beam steering algorithms, and new propagation and partition loss models above 100 GHz are given, and promising imaging, array processing, and position location results are presented. The implementation of spatial consistency at THz frequencies, an important component of channel modeling that considers minute changes and correlations over space, is also discussed. This paper offers the first in-depth look at the vast applications of THz wireless products and applications and provides approaches ... (View more) Topic: Millimeter-wave and Terahertz Propagation, Channel Modeling and Applications Published in: IEEE Access ( Volume: 7 ) Page(s): 78729 - 78757 Date of Publication: 06 June 2019 Electronic ISSN: 2169-3536 INSPEC Accession Number: 18770124 DOI: 10.1109/ACCESS.2019.2921522 Publisher: IEEE Funding Agency: Prof. Theodore S. Rappaport gave a talk based on this paper about Wireless Communication and Applications above 100 GHz, with a focus on opportunities and challenges for ... View more Hide Full Abstract SECTION I.Introduction The tremendous funding and research efforts invested in millimeter wave (mmWave) wireless communications, and the early success of 5G trials and testbeds across the world, ensure that commercial wide-spread 5G wireless networks will be realized by 2020 [1]. The use of mmWave in 5G wireless communication will solve the spectrum shortage in current 4G cellular communication systems that operate at frequencies below 6 GHz [2]. However, the increasing number of new applications such as virtual/augmented reality (VR/AR), autonomous driving, Internet of Things (IoT), and wireless backhaul (as a replacement for labor-intensive installation of optical fiber) [3], [4], as well as newer applications that have not been conceived yet, will need even greater data rates and less latency than what 5G networks will offer. Today, within the global unlicensed wireless mmWave band of 60 GHz, there is ~7 GHz of bandwidth available [5], and in such a wide bandwidth, data rates on the order of 100 Gigabits per second (Gbps) can only be achieved with transmission schemes having a spectral efficiency of at least 14 bit/s/Hz, which requires symbol fidelity that is not feasible using currently known digital modulation techniques or transceiver components [6]–[7][8]. Therefore, data rates on the order of 100 Gbps or more will flourish at frequencies above 100 GHz, where the available spectrum is massively abundant [9]. Fig. 1 illustrates the applications and range of frequencies available from the sub-THz regime up through and beyond the optical spectrum, and shows how mmWaves and THz frequencies are three and two orders of magnitude, respectively, below the frequencies of visible light. At optical and infrared frequencies, issues like the impact of atmospheric and water absorption on the signal propagation, ambient sunlight, required low transmission power budget due to eye-safety limits, and high diffusion losses on rough surfaces limit their use for wireless communication systems [10]. Ionizing radiation, which includes ultraviolet, x-rays, galactic radiation, and gamma-rays, is dangerous since it is known to have sufficiently high particle energy to dislodge electrons and create free-radicals that can lead to cancer [11], [12] and is believed to be a major health risk for interplanetary travel [12], [13]. The adverse health effects of ionizing radiation may be negligible, however, if used with care [14]. Ionizing radiation can be used for gauging the thickness of metals, Roentgen Stereophotogrammetry, astronomy, nuclear medicine, sterilizing medical equipment, and pasteurizing certain foods and spices [15]. Unlike ionizing radiation, mmWave and THz radiation are nonionizing because the photon energy is not nearly sufficient (0.1 to 12.4 meV, which is more than three orders of magnitude weaker than ionizing photon energy levels) to release an electron from an atom or a molecule, where typically 12 eV is required for ionization [11], [12], [16]. Since ionizing radiation is not known to be a concern at mmWave and THz frequencies, and heating is believed to be the only primary cancer risk [11],[12] the Federal Communications Commission (FCC) and International Commission on Non-Ionizing Radiation Protection (ICNIRP) standards [17], [18] are designed principally to protect against thermal hazards, particularly for the eyes and skin where these tissues are most sensitive to heat from radiation due to lack of blood flow. However, we must point out that with the likelihood of THz sources becoming more widely available, there should be careful work done to understand the biological and molecular impact of THz radiation on human health [12], since, even though THz is more than two orders of magnitude lower in frequency than ionizing radiation, it would be prudent to know with certainty that heating is the only health concern at THz [11]. FIGURE 1. - The electromagnetic spectrum, and various applications as a function of frequency. FIGURE 1. The electromagnetic spectrum, and various applications as a function of frequency. View All FIGURE 2. - Surface scattering across the EM spectrum a) Most surfaces appear smooth at microwave frequencies (specular scattering), b) Same surfaces exhibit significant considerable roughness in optical spectrum (diffuse scattering) c) in THz regime most building surfaces exhibit significant diffuse scattering and strong specular reflections. FIGURE 2. Surface scattering across the EM spectrum a) Most surfaces appear smooth at microwave frequencies (specular scattering), b) Same surfaces exhibit significant considerable roughness in optical spectrum (diffuse scattering) c) in THz regime most building surfaces exhibit significant diffuse scattering and strong specular reflections. View All FIGURE 3. - Leveraging mmWave imaging and communications can position devices with centimeter accuracy, even in NLOS scenarios. Figure (a) shows the experimental setup adopted by [50] to illustrate the capability of mmWave imaging-based positioning. In Figure (b) the user location is projected on the constructed mmWave image. FIGURE 3. Leveraging mmWave imaging and communications can position devices with centimeter accuracy, even in NLOS scenarios. Figure (a) shows the experimental setup adopted by [50] to illustrate the capability of mmWave imaging-based positioning. In Figure (b) the user location is projected on the constructed mmWave image. View All FIGURE 4. - (a) A spatially-oversampled $N_{x}\times N_{y}$ uniform rectangular antenna array (URA). The antenna spacing is $\lambda _{m}/\left ({2K_{x}}\right)$ and $\lambda _{m}/\left ({2K_{y}}\right)$ along the $x$ -and $y$ -axes, respectively. Here $\lambda _{m}$ is the electromagnetic wavelength, $K_{x}\geq 1$ and $K_{y}\geq 1$ are the spatial oversampling factors along $x$ and $y$ , respectively, and assume $K_{x}=K_{y}=K_{u}$ for convenience. (b) The ROS of waves received by the 2-D URA (green) consists of a narrow light cone that is sparse in the 3-D spacetime frequency domain. FIGURE 4. (a) A spatially-oversampled Nx×Ny uniform rectangular antenna array (URA). The antenna spacing is λm/(2Kx) and λm/(2Ky) along the x -and y -axes, respectively. Here λm is the electromagnetic wavelength, Kx≥1 and Ky≥1 are the spatial oversampling factors along x and y , respectively, and assume Kx=Ky=Ku for convenience. (b) The ROS of waves received by the 2-D URA (green) consists of a narrow light cone that is sparse in the 3-D spacetime frequency domain. View All FIGURE 5. - An overview of a space-time $\Delta -\Sigma $ ADC where spatio-temporal noise-shaping occurs along both spatial- and temporal-frequency domains. The analog filtering module inside the multi-port ADC couples the inputs and output ports in such a manner that the quantization noise arising from the 1-bit quantizers are shaped to be outside the $K_{u}$ -compressed “cone of silence” of the array. The plane-wave signals of interest, on the other hand, now lie inside the “cone of silence” and can be extracted by a sharp digital beamformer whose minor lobes will significantly attenuate the quantization noise of the system. FIGURE 5. An overview of a space-time Δ−Σ ADC where spatio-temporal noise-shaping occurs along both spatial- and temporal-frequency domains. The analog filtering module inside the multi-port ADC couples the inputs and output ports in such a manner that the quantization noise arising from the 1-bit quantizers are shaped to be outside the Ku -compressed “cone of silence” of the array. The plane-wave signals of interest, on the other hand, now lie inside the “cone of silence” and can be extracted by a sharp digital beamformer whose minor lobes will significantly attenuate the quantization noise of the system. View All While 5G, IEEE 802.11ay, and 802.15.3d [19], [20] are being built out for the mmWave spectrum and promise data rates up to 100 Gbps, future 6G networks and wireless applications are probably a decade away from implementation, and are sure to benefit from operation in the 100 GHz to 1 THz frequency bands where even greater data rates will be possible [3], [7], [10]. The short wavelengths at mmWave and THz will allow massive spatial multiplexing in hub and backhaul communications, as well as incredibly accurate sensing, imaging, spectroscopy, and other applications described subsequently in this paper [21]–[22][23][24]. The THz band, which we shall describe as being from 100 GHz through 3 THz, can also enable secure communications over highly sensitive links, such as in the military due to the fact that extremely small wavelengths (orders of microns) enable extremely high gain antennas to be made in extremely small physical dimensions [25]. Although we note the formal definition of the THz region is 300 GHz through 3 THz, some have begun to use the terms “sub-THz” or “sub-mmWave” (e.g. using frequency or wavelength) to define the 100–300 GHz spectrum. There are tremendous challenges ahead for creating commercial transceivers at THz frequencies, but global research is addressing the challenges. For example, the DARPA T-MUSIC program is investigating SiGe HBT, CMOS/SOI and BiCMOS circuit integration, in hopes of achieving power amplifier threshold frequencies ft of 500–750 GHz [26]. A survey of power amplifier capabilities since the year 2000 is given in [27]. It should be clear that the semiconductor industry will solve these challenges, although new architectures for highly dense antenna arrays will be needed, due to the small wavelengths and physical size of RF transistors in relation to element spacing in THz arrays. Section III provides some promising design approaches for future digital arrays. Since there is very high atmospheric attenuation at THz band frequencies, especially at frequencies above 800 GHz (see Fig. 6), highly directional “pencil” beam antennas (antenna arrays) will be used to compensate for the increased path loss due to the fact that the gain and directivity increase by the square of the frequency for a fixed physical antenna aperture size [6], [28]. This feature makes THz signals exceedingly difficult to intercept or eavesdrop [4], [10], [25], [29]. However, a narrow pencil-like beam does not guarantee immunity from eavesdropping, and physical-layer security in THz wireless networks and transceiver designs that incorporate new counter-measures for eavesdropping will be needed [30]. FIGURE 6. - Atmospheric absorption beyond the natural Friis free space loss of electromagnetic waves at sea level versus frequency under different humidity conditions [101], [104]. FIGURE 6. Atmospheric absorption beyond the natural Friis free space loss of electromagnetic waves at sea level versus frequency under different humidity conditions [101], [104]. View All Energy efficiency is always important for communication systems, especially as circuitry moves up to above 100 GHz, and a theoretical framework to quantify energy consumption in the presence of vital device, system, and network trade-offs was presented in [31], [32]. The theory, called the consumption factor theory (CF, with a metric measured in bps/W), provides a means for enabling quantitative analysis and design approaches for understanding power trade-offs in any communication system. It was shown in [31], [32] that the efficiency of components of a transmitter closest to the output, such as the antenna, have the largest impact on CF [31]. The power efficiency increases with increasing bandwidth when most of the power used by components that are “off”, e.g., ancillary, to the signal path (e.g., the baseband processor, oscillator, or a display) is much greater than the power consumed by the components that are in line with the transmission signal path (e.g., power amplifier, mixer, antenna) [31], [32]. For a very simple radio transmitter, such as one that might be used in low cost IoT or “smart dust”applications where the power required by the ancillary baseband processor and oscillator is small compared to the delivered radiated power, the power efficiency is independent of the bandwidth [31]. Thus, contrary to conventional wisdom, the CF theory proves that for antennas with a fixed physical aperture, it is more energy efficient to move up to mmWave and THz frequencies which yield much wider bandwidth and better power efficiency on a bits per second per watt (bps/W) basis, as compared to the current, sub-6 GHz communication networks. Global regulatory bodies and standard agencies such as the FCC [33], the European Telecommunication Standards Institute (ETSI) [34], and the International Telecommunication Union (ITU) [35], are seeking comments to allocate frequency bands above 95 GHz for point-to-point use, broadcasting services, and other wireless transmission applications and use cases [36]–[37][38][39]. In fact, in March 2019, the FCC voted to open up spectrum above 95 GHz for the first time ever in the USA, and provided 21.2 GHz of spectrum for unlicensed use shown in Table 1, and permitted experimental licensing up to 3 THz [40]. The mmWave coalition [41], which is a group of innovative companies and universities united in the objective of removing regulatory barriers to technologies using frequencies ranging from 95 GHz to 275 GHz in the USA, submitted comments to the FCC and to the National Telecommunications and Information Administration (NTIA) for developing a sustainable spectrum strategy for America’s future, and urged NTIA to facilitate greater access to spectrum above 95 GHz for non-Federal use in January 2019 [41]. The Institute of Electrical and Electronics Engineers (IEEE) formed the IEEE 802.15.3d [20] task force in 2017 for global Wi-Fi use at frequencies from 252 GHz to 325 GHz, creating the first worldwide wireless communications standard for the 250–350 GHz frequency range, with a nominal PHY data rate of 100 Gbps and channel bandwidths from 2 GHz to 70 GHz [20]. The use cases for IEEE 802.15.3d include kiosk downloading (dubbed the “Information Shower” by an author of this paper) [42], intra-device radio communication [43], connectivity in data centers, and wireless fiber for fronthaul and backhaul [20], [41], [44]. As indicated in [41], an interesting issue is the system simulation and fabrication of antennas to ensure that ground-based fixed and mobile sub-THz systems do not interfere with satellites and space-based sensors operating in the same sub-THz bands above earth. TABLE 1 Unlicensed Spectrum Proposed by FCC [33] Table 1 This paper is organized as follow: Section II provides insights into mmWave and THz applications. Section III explores the possibility of using array signal processing techniques for improved front-end performance at frequencies above 100 GHz. Section IV demonstrates viability for true time delay (TTD) beamformers for wideband digital arrays. Section V reviews key wireless propagation fundamentals and presents measurements and research above 100 GHz which could be used for many novel applications, including future 6G communications. Channel sounding systems and measured results at above 100 GHz are given in Section VI, respectively. 28, 73 and 142 GHz channels are compared in Section VII. Scattering mechanism at mmWave and THz frequencies is explained in Section VIII. Section IX investigates positioning approaches (e.g., position location that can be used in navigation) at THz frequencies that promise unprecedented accuracy. Section X discusses the implementation of spatial consistency, an important component of channel modeling at THz frequencies, and Section XI concludes the paper and summarizes key areas that warrant future research. SECTION II.mmWave and THz Applications The ultra-high data rates facilitated by mmWave and THz wireless local area and cellular networks will enable super-fast download speeds for computer communication, autonomous vehicles, robotic controls, the information shower [42], high-definition holographic gaming, entertainment, video conferencing, and high-speed wireless data distribution in data centers [10]. In addition to the extremely high data rates, there are promising applications for future mmWave and THz systems that are likely to evolve in 6G networks, and beyond. These applications can be categorized into main areas such as wireless cognition, sensing, imaging, wireless communication, and position location/THz navigation (also called localization, or positioning), as summarized in Table 2. TABLE 2 Promising Applications at mmWave and THz Table 2 A. Wireless Cognition Wireless cognition is the concept of providing a communication link that enables massive computations to be conducted remotely from the device or machine that is doing real time action [46]. For example, a lightweight drone fleet may not have the power or weight budget to conduct massive computations on board the apparatus, but with a wide enough channel bandwidth and sufficiently fast data rate, real time computations for extremely complex tasks, such as contextual awareness, vision, and perception may be carried out at a fixed base station or edge server that is in wireless connection and supporting real time cognition for the drone fleet. Robots, autonomous vehicles, and other machines may be similarly designed to exploit cognitive processing performed remotely from the machine using wireless, with the ability to perform tasks without the benefit of local cognition on the platform [45], [46]. When one considers the growth of computational power provided by Moore’s law, it can be seen that the modest price of 1000 USD (the cost of today’s smart phone) will likely be able to purchase a computer with computational capabilities that are on the order of the human brain by the year 2036. This observation stems from results in [54], [55] which extrapolate the increase in computational capabilities over time since 1965 [56]. The following analysis shows that Terahertz frequencies will likely be the first wireless spectrum that can provide the real time computations needed for wireless remoting of human cognition [52]. There are about 100 billion (1011) neurons in the human brain, each of which can fire 200 times per second (5 ms update rate), and each neuron is connected to about 1000 others, resulting in a computation speed of 20×1015 floating-point operations per second (flops) [54], [55], which, if each operation is assumed to be binary, will require a data rate of 20,000 Tbps: ===Human Brain flops (Computation Speed)1011 neurons×200 flop/sec×103/neuron20×1015 flop/sec=20 petaflops/sec×1 bit/flop20,000 Tbps.(1) View SourceRight-click on figure for MathML and additional features. Each neuron has write access to 1000 bytes resulting in a memory size of the human brain of 100 Terabytes [54], [55]: Storage==1011 neurons×103 bytes/neuron1014 bytes=100 TB.(2) View SourceRight-click on figure for MathML and additional features. Today’s state-of-the art 1000 USD computer technology performs one trillion (e.g., 1012) computations/sec, which is four orders of magnitude less than the speed of the human brain. Future wireless generations (e.g., 6G or 7G) are likely to allocate up to 10 GHz RF channels for each user in the THz regime, and by assuming that each user is able to exploit 10 bits/symbol modulation methods and 1000 times increase in channel capacity using yet-to-be-invented concepts beyond cooperative multipoint (CoMP) and Massive-MIMO, it is readily seen in (3) that data rates of 100 Terabytes/sec will be achieved. R==10 GHz channel×10 bits/(sec cdot Hz)×103100 Tbps.(3) View SourceRight-click on figure for MathML and additional features.From (1) and (2), it is clear that a 100 Tbps link is plausible in a 10 GHz channel bandwidth, providing 0.5% of real time human computational power. Ambitiously, if 100 GHz channel bandwidths are used, 1 Petabits/sec of information, or 5% of the real time computational power of the human brain, could be carried over wireless [52]. B. Sensing Sensing applications enabled by mmWave and THz frequencies exploit the sub-mm wavelength and frequency selectivity of the measured environment to gain knowledge about the environment based on the observed signature of the propagated signal. Sensing applications may exploit both the vastly wider channel bandwidths at above 100 GHz, as well as the frequency-selective resonances and absorption of various materials. The ability to implement very high gain antennas in a small physical form factor also enables sensing applications that are very directional, and as frequencies increase, the spatial resolution becomes much finer as a function of wavelength, thereby enabling sub-millimeter spatial differentiation as frequencies exceed 300 GHz. Through beam scanning, it will become possible to create images of physical spaces, through the systematic monitoring of received signal signatures at a wide array of different angles. Because electrical beam steering algorithms can be implemented in real time (sub microsecond) and radio propagation distances are small (e.g., on the order of meters in a room) thus lending to propagation times that are less than 10 ns, it shall be feasible to measure the properties of a room, an office, or a complex environment in a matter of seconds or less. This ability opens up a new dimension of wireless, that enables future wireless devices to do “wireless reality sensing” and gather a map or view of any location, leading to detailed 3-D maps of the world created on the fly and uploaded/shared on the cloud by future devices [57], [58]. Also, since certain materials and gasses have vibrational absorption (e.g., resonances) at particular frequencies throughout the THz band, it becomes possible to detect the presence of certain items based on frequency scanning spectroscopy. The presence of certain chemicals or allergens in food or drink or air, or other defects in the world around us, may be sensed based on spectroscopy that could be enabled in smart phones for 6G and beyond. THz will enable new sensing applications such as miniaturized radars for gesture detection and touchless smartphones, spectrometers for explosive detection and gas sensing [49], THz security body scanning, air quality detection [24], personal health monitoring systems [48], precision time/frequency transfer, and wireless synchronization [10], [22], [59]. By building real time maps of any environment, it may also be possible to predict channel characteristics at a mobile device, to aid in the alignment of directional antennas, to provide on-the-fly position location, as well as adaptation of wireless capabilities [53]. This capability could also be fed to the cloud to enable a real time collection capability for mapping and sensing the world, which then could be used in commercial applications for transportation, shopping, and other retail uses. C. Imaging Radar at mmWave and THz frequencies are more effective than light or infrared-based imaging such as Light Detection and Ranging (LIDAR), due to the smaller impact that weather and ambient light have on the THz channel. It is worth noting that although LIDAR can provide higher resolution, LIDAR cannot work when it is foggy, raining, or cloudy [60]. However, mmWave and THz radar can be used for assisting driving or flying in foul weather, as well as in military and national security [10], [21], [25]. High-definition video resolution radars that operate at several hundred gigahertz will be sufficient to provide a TV-like picture quality and will complement radars at lower frequencies (below 12.5 GHz) that provide longer range detection but with poor resolution [21], [23]. Dual-frequency radar systems will enable driving or flying in very heavy fog or rain [21]. THz waves can augment human and computer vision to see around corners and to “view” NLOS objects, which enables unique capabilities in rescue and surveillance, autonomous navigation, and localization. A building surface (e.g., wall, floor, door) typically behaves to a first order as mirrors (e.g., perfect reflectors of THz energy), thus allowing THz imaging to see around corners and behind walls if there are sufficient reflection or scattering paths (see Section VIII) [61], [62]. NLOS imaging methods based on visible and infrared light have also been proposed in [63]–[64][65]. However, the optical wavelength is smaller than the surface roughness of most surfaces, thus optical NLOS imaging requires complex hardware and computationally expensive reconstruction algorithms, while exhibiting short imaging distances (< 5 m). In general, the practical deployment of visible light systems is still not clearly evolved due to weakly scattered signals, small field of view, and long integration times. Alternatively, lower frequency (< 10 GHz) NLOS radar systems suffer less loss and objects appear relatively much smoother. However, in the lower frequency spectrum, materials are semitransparent, edge diffraction become much stronger, and images are easily cluttered due to strong multi-reflection propagation [66]–[67][68][69][70]. Additionally, radar systems require accurate knowledge of the static geometry and are limited to object detection rather that detailed images of the hidden scene. THz waves combine many of the advantageous traits of both microwaves and visible light. Namely, they feature small wavelengths and wide bandwidths that allow high spatial resolution images with moderately sized imaging systems [71], [72]. Additionally, THz scattering exhibits both specular and diffuse scattering from most building surfaces, as depicted in Fig. 2 (see Section VIII for further details). The strong specular component turns the surfaces into something close to “electrical mirrors”, thus allowing to image objects around obstacles, while maintaining spatial coherency (narrow beams) and high spatial resolution [62]. Radar imaging systems illuminate the scene with THz waves and generate 3D images by computing the time of flight (ToF) of the backscattered signals. When the path of the scattered signal involves multiple bounces on surrounding surfaces, the resulting 3D image appears distorted. If the LOS surfaces act as mirrors due to the strong specular reflections, a corrected image of the NLOS objects can be reconstructed by applying relatively simple mirroring transformations [61]. D. Precise Positioning Leveraging mmWave imaging and communications for centimeter-level localization is an exciting feature that will likely be incorporated in future portable devices operating above 100 GHz [52]. A simple example of what is possible was illustrated in [50] using an experimental setup, as shown in Fig. 3(a). The experiment demonstrates excellent accuracy when a LOS does not exist, and illustrates how positioning and sensing will become part of communication devices. Figure 3(a) shows the setup, which consists of (i) a synthetic aperture radar (SAR) emulating a 13-cm aperture linear antenna array and operating at 220–300 GHz, (ii) two drywall sheets, and (iii) a single-antenna user that is in the NLOS of the SAR. First, a mmWave image (e.g., map) of the surrounding environment, including both LOS and NLOS objects, is constructed by rapid beam steering and analyzing the radar returns at very narrow beam angle separations. Second, the user sends an uplink pilot that is used by the antenna array, for example at a base station or an access point, to estimate the angle and time of arrival (AOA and TOA). Finally, the AOA and TOA are used to back-track the path of the user signal on the constructed mmWave image, as shown in Fig. 3(b), to determine the user location (see Section II-D for more details). Considering Fig. 3(a), this mmWave imaging/communication-based method was able to localize the user within a 2 cm radius over a distance of 2.8 m, and as shown in Section IX, excellent accuracy are possible over much greater distances. Leveraging mmWave and THz imaging for localization has interesting and unique benefits compared to other methods. The mmWave imaging/communication methods proposed in [50], [53], [58] can localize users in the NLOS areas, even if their travel paths to the base station/access points experience more than one reflection (e.g., multiple bounces). For example, the experiment in Fig. 3(a) was able to localize the user even though its path experiences two bounces. Further, different from classical simultaneous localization and mapping (SLAM) methods, which require prior knowledge and calibration of the environment, the proposed mmWave imaging/communication-based technique in [50] does not require any prior knowledge of the environment. By building or downloading the map of the environment, the mobile device will be able to exploit many other capabilities, such as predicting the signal level, using real time site-specific prediction [57], or uploading the map to the cloud that compiles physical maps, or which uses such maps for mobile applications [50], [53], [58]. Finally, due to leveraging the large bandwidth available at frequencies above 100 GHz, the LOS and NLOS users can be localized with centimeter accuracy. These promising user features highlight the potential of merging mmWave imaging and communications for precise positioning in 6G systems and beyond. Using mmWave or THz imaging (e.g., a mmWave camera as discussed in Section II-C) to reconstruct 3D maps of the surroundings in unknown environments, it will be possible to merge sensing and imaging and position location all at the same time. As shown subsequently, mmWave and THz signals reflect strongly from most building materials which enable imaging of hidden objects (NLOS imaging) [50], [53], [58], although as shown in Section VIII, scattering can also be well modeled and predicted. Based on the 3D maps of the physical surroundings, and the time and angular information from a mobile (TOA and AOA), centimeter level localization and mapping are achievable with the massive bandwidth and large antenna arrays at mmWave and THz frequencies [50], [53], [58]. SECTION III.Multi-Port Array Transceiver Architectures Based on Space-Time Sigma-Delta Theory The adoption of mmWave and THz frequency bands above 100 GHz may require new approaches to transceiver front-end design, due to the challenges of power efficiency and the vastly increased number of antennas that will be possible in a small physical area. In fact, the physical size of RF amplifiers may not be small enough to be physically located behind individual antenna elements at THz frequencies, possibly necessitating hybrid beamforming approaches for practical implementations [73], [74]. Key transceiver performance metrics such as i) noise figure (NF) and linearity of the low-noise amplifier (LNA), and ii) the peak output power Psat and power-added efficiency (PAE) of the power amplifier, degrade rapidly with increasing operating frequency f0 [6], [27], [75]. For example, a survey of the recent literature shows the product of Psat , gain, and PAE scales as 1/f20 [27]. Phase noise also increases with frequency, creating increased noise levels for practical coherent demodulation. Thus, it becomes challenging to provide sufficient performance at these frequencies to facilitate higher-order modulations, such as those adopted by modern 5G wireless and Wi-Fi systems. The larger number of antennas used at mmWave and THz frequencies begs consideration of array signal processing techniques that can lead to improved antenna array performance in the face of noise and linearity limitations for mmWave and THz transceivers. In this section, we illustrate a promising new technique of spatially-oversampled antennas, and present new phased array device architectures and compact computational approaches that can provide accurate beam steering with exponentially less power and physical chip real estate than today’s methods of phased array beamforming. A mmWave digital phased array transmits and receives far-field propagating electromagnetic plane-waves. One promising viewpoint is to consider the spatio-temporal frequency-domain region of support (ROS) of these plane-waves as a region of causality that contains all possible signals of interest seen by the array, while ignoring all other components outside of the ROS [76]–[77][78][79][80]. This region of causality, which is a frequency-domain version of the light-cone that is well-described by the Special Theory of Relativity, is a conical region of the spatio-temporal frequency domain inside of which all plane-waves for a given space-time beam pattern are constrained to exist (see Fig. 4). For a Nyquist-sampled planar antenna array in the x−y plane, the cone-shaped frequency domain region defined by the surface ω2x+ω2y=ω2ct contains all propagating RF beams used for wireless communications, where ωx and ωy are spatial frequencies for the wave along the x and y axes, while ωct is its temporal frequency in units such that the speed of light is unity. When the antenna arrays are spatially-oversampled (see Fig. 4(a)), the cone-shaped ROS becomes correspondingly narrower, with a surface given by ω2x+ω2y=ω2ct/K2u where Ku∈ℤ+ is the spatial over-sampling factor (see Fig. 4(b)). Specifically, the opening angle of the cone is given by θ′=tan−1(1/Ku) (see Fig. 4(b)), which decreases as Ku increases. Thus, the fractional volume of 3-D spacetime occupied by the cone-shaped ROS of all propagating waves decreases with Ku . In other words, spatial oversampling of the antenna array makes the cone increasingly sparse within 3-D spacetime. As a simple numerical example, consider a carrier frequency of 300 GHz, corresponding to a free-space wavelength of λ=1 mm. A conventional Nyquist-sampled antenna array would have an antenna spacing of λ/2=0.5 mm, resulting in an open cone angle θ′=45∘ . By contrast, a 4× spatially-oversampled array would have an oversampled antenna spacing of λ/8=0.125 mm, resulting in θ′=14.0∘ . Such sparsity does not allow the array to generate sharper beams (since beam width is ultimately limited by the total aperture), but does result in other benefits in terms of reduced noise figure (NF), increased linearity, and a relaxing of the requirement of analog-to-digital converters or digital-to-analog converters (ADCs/DACs) resolution, as described in more detail in the next sub-section. A. Cones of Silence in Antenna Arrays For a given spatial over-sampling factor Ku , all possible electromagnetic plane-waves and their beam-shaped wideband spectral ROSs across an antenna array lie inside the cone-shaped region of causality. Because of relativistic constraints, no propagating waves can exist outside this region. This is a direct result of the nature of the wave equation as observed in the frequency domain. Thus, the rest of the spectral domain is devoid of any propagating waves. Hence we denote the cone-shaped ROS defined by a particular antenna array’s over-sampling factor as an array’s “cone of silence”, and our design mission becomes one of shifting all possible noise and distortion to outside of this cone. Such cones directly arise from fundamental physics, namely the causality of electromagnetic plane-waves in the spacetime frequency domain (ωx,ωy,ωct)∈ℝ3 [76]–[77][78][79][80], and are thus present for all uniformly-spaced rectangular arrays. We now show that these naturally occurring “cones of silence” can be quite effective in improving the performance of antenna array electronics as the industry moves to THz bands. B. Moving From N 2-Port Tansceivers to One Multi-Port Transceiver for Array Processing In conventional digital phased array transceivers, a single transceiver is replicated N times to form an N -element array. In the proposed approach, these N individual transceivers are replaced by a single multi-port transceiver. For example, N conventional multi-bit ADCs in an N element array are replaced by a single KuN -port ADC which uses 1-bit resolution. The multi-port ADC takes in KuN number of analog inputs and produces KuN number of 1-bit digital outputs. The algorithm, implemented in circuit form, shapes noise and distortion in both space- and time-frequency domains, such that the dominant parts of noise and distortion energy lies outside the “cone of silence” of the signals of interest. The fact that the noise and distortion has been shaped to be outside the “cone of silence”, while the signal of interest itself lies inside it, allows lower-resolution ADC outputs to be employed in an KuN element digital beamformer that generates outputs with a higher digital resolution than in conventional N -element array receivers. Such multi-port flash ADCs have the potential to achieve ≈2Ku (i.e., exponential) savings in number of comparators at a Ku (i.e. linear) cost in number of antennas. The idea of noise shaping is motivated by the fact that, at sub-THz carrier frequencies with basebands spanning tens of GHz, fully-digital arrays having KuN number of “crude 1-bit ADCs making up the multi-port ADC” is possible, unlike building N high-resolution ADCs for conventional arrays, which is more extensive from a power and part count perspective. The noise-shaping approach allows a spatially over-sampled dense 1-bit receiver array to have the same quantization noise performance as a 3–4 bit resolution array which is Nyquist spaced. C. Spatially Noise-Shaped Circuits for Antenna Arrays The fact that “cones of silence” exist for all phased arrays and become increasingly sparse in 3-D spacetime as Ku increases was previously exploited to improve array performance, but offers even greater promise for mmWave and Terahertz phased array antennas. As was first shown in [81] and later expanded upon in [82]–[83][84][85][86], it is possible to exploit such multi-dimensional spatio-temporal feedback between the individual transceiver channels, i.e., along both spatial as well as temporal dimensions. These feedback loops emulate the performance of Δ−Σ algorithms, which are widely used to improve the resolution of ADCs and digital-to-analog converters (DACs) [6], [87]. Fig. 5 shows an example Δ−Σ ADC and LNA compared with the topology of a conventional digital phased array. A critical and valuable practical result of using the phased array architecture in Fig. 5 is that it becomes possible to use RF circuitry (which may be physically large compared to antenna element spacing) behind numerous antenna elements (either in a true digital or hybrid beamforming implementation). Conventional array designs simply may not be able to accommodate, physically, the RF transistors in a phased array package. A conventional Δ−Σ ADC is designed to separate desired signals (i.e., the inputs to be digitized) from undesired ones (i.e., quantization noise) along the temporal frequency axis, thus increasing the ADC’s precision beyond that of its quantizer. In fact, many Δ−Σ ADCs use a 1-bit quantizer (i.e., a single comparator) but nevertheless provide several bits of effective resolution (i.e., effective number of bits or ENOB≫1 ). Such improved performance arises by combining two key principles: i) temporal oversampling to expand the frequency range over which the ADC’s quantization noise is spread, thus reducing its power spectrum density (PSD); and ii) noise shaping to actively filter quantization noise out of the small frequency range (at low frequencies) within which the input signals of interest must be confined (due to oversampling). Note that a Δ−Σ DAC operates in a fundamentally similar way, but with the roles of the analog input and digital output signals reversed. The key point exploited in [81]–[82][83][84][85][86] is that the potential benefits of i) spatially oversampling an antenna array, and ii) shaping unwanted thermal and quantization noise of the array electronics away from the “cone of silence” occupied by propagating electromagnetic waves, are exactly analogous to those provided by temporal oversampling and noise shaping in a conventional Δ−Σ ADC or DAC. In particular, it becomes clear that it is possible to add appropriate spatial-domain, and spatio-temporal domain feedback loops (i.e., carefully-designed electrical coupling networks between the amplifiers and data converters present in adjacent transceiver channels as shown in Fig. 5) that mimic the noise shaping provided by a Δ−Σ converter, but in the spatio-temporal (rather than temporal) frequency domain. As a result, the thermal noise, quantization noise, and nonlinear distortion added by low-noise amplifiers (LNAs), ADCs, and DACs within the array can be actively shaped to lie outside the array’s “cone of silence”, effectively separating them from the propagating electromagnetic waves received or transmitted by the antennas. As a result, spatial noise shaping improves all the important performance metrics for a receive array: it reduces NF, greatly increases linearity, and significantly increases ENOB [85], [86]. These improvements enable receiver power consumption to be reduced for a given set of performance metrics. For example, one can use individually noisy and nonlinear amplifiers and low-resolution quantizers to emulate high-performance versions, thus greatly reducing hardware complexity and overall power consumption. As a result, the design and cooling requirements of dense multi-beam arrays can be greatly simplified. In practice, the available performance improvements due to spatial and spatio-temporal noise-shaping are eventually limited by the increase in mutual coupling between the antenna elements as the spatial oversampling factor Ku increases. Intuitively, a highly-oversampled array approaches a continuous aperture, implying that no further improvements in NF, linearity, or ENOB are possible [85], [86]. However, both simulations and experiments confirm that significant performance benefits are achievable using common antenna elements (e.g., patches and dipoles) for moderate oversampling ratios (Ku≤4 ). SECTION IV.Novel Wideband Multi-Beamformers at Low Size, Weight, and Power (SWaP) We now demonstrate promising new approaches in computation for obtaining beamforming array patterns with much less computation power and time than conventional methods. An N -element uniformly-spaced linear antenna array can be used to realize an N -beam multi-beamformer by computing N -element discrete Fourier transforms (DFTs) along the spatial dimension. These computations can be efficiently performed using a spatial Fast Fourier transform (FFT), which can be implemented using both analog as well as digital approaches [88] and reduces the brute force computational complexity of the DFT from (N2) down to (NlogN) . Moreover, approximate versions of the DFT can be computed with even lower complexity than the FFT, as shown by Cintra, Bayer, and their collaborators [89]–[90][91][92][93]. The Cintra-Bayer approach to approximate-DFT computation allows significantly smaller arithmetic complexity over FFT-based approaches (in some cases, reducing the computation burden by an order of magnitude for large N ) at the cost of only ~2 dB penalty in worst-case sidelobe level [89]–[90][91][92][93]. Note that while today’s mobile devices do not concern themselves with the power savings due to FFT computations (because fidelity of the modulation is vital for bit error rate), the same may not be true for future wireless devices at mm-wave and THz. In fact, the availability of massive channel bandwidths and likely increase in sensitivity to power consumption (e.g., prolonged battery life) may remove the present-day need for perfect fidelity in DFT/FFT processing and making approximate versions attractive [93]. However, both DFT/FFT and approximate-DFT approaches are narrowband digital beamformers and are therefore not suitable for emerging wideband 5G/6G digital arrays [94]. In fact, true time delay (TTD) beamformers are required for achieving electronically-steered wideband RF beams. Unlike phased-array (i.e., DFT-based) beamformers, which are limited to a relatively small fractional bandwidth because of “beam squint” [95], a TTD beamformer achieves frequency independent wideband RF beams as demanded by 5G and impending 6G applications. However, the digital arithmetic complexity required to generate a massive number of TTD beams can be quite staggering compared to the traditional narrowband methods used in digital phased-arrays. Thus, while the use case for generating a large number of TTD digital beams in emerging wireless systems is a strong one, the resulting arithmetic circuit complexity and power consumption often becomes a limiting factor. A novel approach to multi-beam TTD beamforming has been explored, where the delay Vandemonde Matrix (DVM) and its factorization [96], [97], has the potential to significantly lower the net arithmetic complexity of the underlying digital system, well below the brute-force matrix computations required for the direct-computed DVM approach [95]. Extensions to FFT-like algorithms for the case of the multi-beam DVM-based TTD beamformer are the subject of ongoing research. In particular, fast algorithms for DVM factorization have already been found for relatively small N , while more general versions applicable to larger values of N are being actively developed [95]–[96][97]. What is fascinating to note is that the use of TTD approaches has also found promise in mmWave and THz circulators and duplexers, which are non-linear devices that allow one-way transmission for coupling multiple signals to a single RF port [98]. Thus, it is clear that TTD approaches for both antennas and circulators/duplexers offer great promise in reducing noise and increasing isolation and linearity for future mm-wave and THz electronics as compared to the electronics used in present-day wireless systems [99]. SECTION V.Moving to Beyond 5G and Above 100 GHz – A New Type of Channel The importance of atmosphere absorption effects on the transmission of high frequency signals has long been recognized [100], [101]. At lower frequencies (below 6 GHz), the attenuation of a propagating wave is mainly caused by molecular absorption in free space, which is minimal [102], but at higher frequencies, as the wavelength approaches the size of dust, rain, snow, or hail, the effects of Mie scattering become more severe [6], [103], [104]. The various resonances of oxygen, hydrogen, and other gases in the atmosphere, however, cause certain frequency bands to suffer from significant signal absorption. Fig. 6 illustrates how the frequency bands of 183 GHz, 325 GHz, 380 GHz, 450 GHz, 550 GHz, and 760 GHz suffer much greater attenuation over distance in typical air due to the atmospheric absorption, which is in addition to the natural Friis’ free space loss, making these particular bands well suited for very short range and secure communications such as “whisper radio” applications, where massive bandwidth channels will attenuate very rapidly out to a few tens of meters, meters, or fractions of a meter [2], [4], [6], [25], [104]. Fig. 6 also shows many mmWave and THz bands suffer surprising little loss as compared to sub 6-GHz bands, perhaps less than 10 dB/km of additional loss than caused by free space propagation in air up to 300 GHz. These frequency bands could easily be used for high speed 6G mobile wireless networks with up to km size coverage range, and perhaps even up to 10 km or beyond in fixed applications [4], [25], [105]–[106][107]. In fact, Fig. 6 shows that much of the spectrum between 600 and 800 GHz suffers 100 to 200 dB/km attenuation, which is only 10 to 20 dB over a 100 m distance (the typical radius of a small cell). As shown subsequently, high gain electronically steered directional antennas, that are physically small, will easily overcome this atmospheric attenuation, meaning that the mobile industry will be able to work well up to 800 GHz in the future using the small-cell architectures envisioned for 5G if line-of-sight or strong reflectors exist. Work in [108] showed that particles such as rain drops, snow, and hail caused substantial attenuation at frequencies above 10 GHz. 73 GHz signals attenuate at 10 dB/km when the rain rate is 50 mm/h [2], [6], [109]. Interestingly, as shown in [2], [101], [110] rain attenuation flattens out from 100 GHz to 500 GHz, meaning that rain will not cause any additional attenuation at operating frequencies above 100 GHz (see Fig. 7). For all mmWave frequencies where the urban cell sizes are on the order of 200 m, rain or snow attenuation can clearly be overcome with additional antenna gain (obtainable by switching in more antenna array elements) [2], [4], [6]. FIGURE 7. Rain attenuation flattens out above 100 GHz. The rain attenuation at 1 THz is 10 dB/km for moderate rainfall of 25 mm/h, only 4 dB/km more than the rain attenuation experienced at 28 GHz [2], [101], [110]. View All Work in [102], [111]–[112][113][114][115] showed effects of rain, fog, dust, and air turbulence at THz and infrared frequencies based on measurements with less than 10 m links, which indicated a relatively small change (a few dB) in the received power compared to the power received in clear weather. THz outdoor measurements during a snow fall at 200 GHz were presented in [103], and a 2 dB decrease in the received power was observed across an 8 m LOS link indicating that a higher transmitted power was required during a snow fall to maintain the same data rate. As shown in Fig. 7, mmWave and THz wireless channels are expected to experience ~30 dB/km of maximum rain attenuation above 100 GHz. Since the first Fresnel zone is tiny at THz frequencies, the antenna beams of wireless devices are easily blocked [116], [117], which requires rapidly steerable phased array transceivers to support both adequate transmission range and to provide adaptive beam steering [2], [8], [60], [110]. The array angular resolution approximates λ/L , where λ is the wavelength at the carrier frequency and L is the aperture of the array. For example, consider a fixed wireless system or mobile base station using a 340 GHz carrier with a 35 cm (14 inch) aperture – this permits a sharp 0.14 ° beam resolution [21]. Three system examples – a 140 GHz, 10 Gbps picocell backhaul unit using 64-element arrays; a 340 GHz, 8-channel 160 Gbps MIMO backhaul link with each channel comprised of a 64-element subarray module; and a heads-up-display 400 GHz automotive imaging radar using a linear 64×512 array and frequency scanned beam steering were presented in [60], thus proving the capability of building THz systems. There is still a recurring misconception among researchers, that the wireless channels at higher frequencies would experience more loss, since they have been taught to only consider the free space path loss (FSPL) with omnidirectional antennas. However, it is important to note that at higher frequencies, the antennas (for a given physical size) will be more directional and have more gain. Friis free space equation (4) with antenna gain (5) readily demonstrate that higher frequency links are viable with less loss, not more [6], [28]. With directional antennas at both link ends, the received power is given by (6). Theoretically, the path loss in free space decreases quadratically as frequency increases, so long as the physical size of the antenna (effective aperture) is kept constant over frequency at both link ends [4], [6], [28], [118]. The astounding result of improved coverage at higher frequencies for a given antenna aperture at both ends of the link is clear in Fig. 8, where the received power at 140 GHz free space is 5.7 dB and 14 dB greater than the power received at 73 GHz and at 28 GHz, respectively, for the same TX output power and for identical physical antenna areas at all frequencies (see Ch.3 in [6]). The better wireless links at higher frequencies enable wider bandwidths to be used while keeping the same signal to noise ratio as at lower frequencies, implying that highly directional steerable antennas will enable mobile systems to overcome air-induced attenuation well into the THz region. FIGURE 8. Received power vs. distance with (i) TX and RX are both directional, (ii) TX is directional but RX is omnidirectional, and (iii) both TX and RX are omnidirectional at 28 GHz, 73 GHz, and 140 GHz. Directional antennas with equal effective aperture (Ae=2.9cm2 ) at both TX and RX have much less path loss at higher frequencies (see Ch.3 in [6] and [28]). View All FIGURE 9. Block diagram of 140 GHz broadband channel sounder system with 4 GHz 3 dB RF bandwidth used at NYU WIRELESS. View All SECTION VI.Channel Sounding and Channel Properties Above 100 GHz There are three main channel sounding techniques in the THz range, namely, THz Time Domain Spectroscopy (THz-TDS), Vector Network Analyzer (VNA) based channel sounding, and correlation based sounding [6], [28], [118]–[119][120][121]. THz-TDS is based on sending ultrashort pulsed laser light from a common source to the transmitter (TX) and the receiver (RX). The TX converts the ultrashort light pulse to the THz range and the detector at the RX transforms the received field strength of the THz impulse into an electrical signal when the optical impulse hits the detector [119], [120]. The short THz-TDS pulses cover a huge bandwidth and are excellent for estimating electrical and scattering parameters of sample materials. However, due to the large size of the spectrometer and the limited output power, it is not suitable to be used over a wide range of indoor or outdoor scenarios or for measuring the wireless channels at more than a few meters of distance. Four-port VNAs are commonly used for THz range channel sounding, where the two additional ports (compared to the traditional two-port VNAs used at lower frequencies) are used to generate a local oscillator for the mixer in the frequency extensions, that are used to elevate the VNA stock frequency range to much higher frequencies through heterodyning [120], [121]. A VNA sweeps discrete narrowband frequency tones across the bandwidth of interest to measure the S21 parameter of the wireless channel. Due to the long sweep time across a broad spectrum which can exceed the channel coherence time, VNA-based channel sounders are typically used in a static environment and require a cable that can be a tripping hazard over tens or hundreds of meters [28], [118], [121]. Correlation based channel sounder systems transmit a known wideband pseudo-random sequence. At the RX, the received signal is cross-correlated with an identical but slightly delayed pseudo-random sequence, providing autocorrelation gain at the expense of a slightly longer acquisition time (on the order of tens of ms) [118], [121], [122]. Sliding correlator chips have recently been produced that offer a 1 Gbps baseband spread spectrum sequence [123], and sliding correlators generally enable cable-free operation over useful mobile communication distances of up to 200 m at sub-THz frequencies, depending on transmit power, bandwidth, and antenna gain [28], [121]. A 140 GHz channel sounder with 1 GHz RF bandwidths has been used to measure indoor channels at NYU [28], with new measurements results given in Section VII. A correlation-based channel sounder at 300 GHz with 8 GHz bandwidth was presented and evaluated in [120] with the same wired clock source being connected to both the TX and RX, which use a subsampling technique to avoid the expense of high-speed A/D converters. A 12th order M-sequence was used with a subsampling factor of 128, and the theoretically maximum measurable Doppler frequency was 8.8 kHz, equivalent to a velocity of 31.7 km/h at 300 GHz [120]. Measurements at 100, 200, 300, and 400 GHz with a 1 GHz RF bandwidth THz-TDS channel sounder showed that both indoor LOS and NLOS (specular reflection from interior building walls) links could provide a data rate of 1 Gbps [104]. Propagation loss measurements for estimating the performance of a communication link in the 350 GHz frequency band were presented in [124], where a VNA based system was used with 26 dBi gain co-polarized horn antennas at both the TX and RX. The presence of water absorption lines in the spectra at 380 GHz and 448 GHz was very evident. Data rates of 1 Gbps for a 8.5 m link and 100 Gbps for a 1 m T-R separation distances were shown to be possible via wireless communication links at 350 GHz [124]. Channel and propagation measurements at 300 GHz were presented in [125], [126], where a VNA based channel system with 26 dBi gain horn antennas at both TX and RX was used to analyze the channel characteristics from 300–310 GHz with an IF frequency bandwidth of 10 kHz. Maximum transmission rates of several tens of Gbps for LOS and several Gbps for NLOS paths were shown to be achievable [126]. THz band indoor propagation measurements were conducted in [127] using a VNA based system covering a frequency range from 260 GHz to 400 GHz with 25 dBi gain horn antennas at both TX and RX within a T-R separation range of 0.95 m. Measurement results showed that Tbps throughput was achievable in the THz band. However, robust beam forming algorithms will be required in THz band communications. Acoustic ceiling panels, which were shown to be good reflectors at THz band, could be used as low cost components to support NLOS links [127]. Absorption coefficients of glass, medium density fiberboard, and plexiglass were measured where the S21 parameter was recorded at a T-R separation distance of 10 cm with and without the obstructions brought into the ray path, showing that a 2.5 cm thick window and a 3.5 cm thick fiberboard door would induce absorption attenuation of approximately 65.5 dB and 86.7 dB at 300 GHz [125], respectively. Second generation transceivers demos for D-band radars with Doppler, loop-back, and 4 Gbps NLOS radio links, over the air and at distances exceeding one meter were presented in [128]. The transceiver chip-set operated at 122 GHz based on low-voltage typologies with a 7.5 dB noise figure and −95 dBc/Hz phase noise while consuming only 220 mW. SECTION VII.Comparing 28, 73, and 142 GHz Channels Propagation measurements in the 140 GHz band were conducted in a shopping mall [129], [130] using a VNA based channel sounder with a 19 dBi horn antenna at the RX, and a 2 dBi bicone antenna at the TX. It was shown that the numbers of clusters and multipath components (MPCs) in each cluster in the 140 GHz band, an average of 5.9 clusters and 3.8 MPCs/cluster, were fewer compared to the 28 GHz band, which had an average of 7.9 clusters and 5.4 MPCs/cluster [129]. Work in [131], [132] presented D-band propagation measurements in a very close-in environment around a personal computer using a VNA based sounder. Indoor directional path losses at 30 GHz, 140 GHz, and 300 GHz were compared using different path loss models in [132]. Although the LOS path loss models predicted path loss exponents (PLE) close to 2.0, the multi-frequency close-in free space reference distance model with a frequency-dependent term (CIF) and alpha-beta-gamma (ABG) model had slightly better PLE and standard deviation stability for the indoor environments than the single-frequency CI and Floating Intercept (FI) models (although a sensitivity analysis such as provided in [133] has not yet been performed for indoor sub-THz channels) [4], [132], [134], [135]. Indoor channel measurements in the 126–156 GHz frequency band were given in [136] using a VNA-based channel sounder and 20 dBi gain horn antennas at both TX and RX within a T-R separation distance of 10.6 m. The PLE for the measured D band signal was close to 2 (free space) since the LOS component is much stronger than the secondary paths. The root mean square delay spread was found to be 15 ns [136]. Channel modeling and capacity analysis for THz band electromagnetic (EM) wireless nano-networks were explored in [137], where a propagation model accounting for the total path loss and the molecular absorption for EM communications in the THz bands (0.1-10 THz) was developed based on radiative transfer theory. Channel capacity of the THz band was investigated by using the proposed path loss model for different power allocation schemes, showing that for several tens of millimeters propagation distance, the THz channel supports very large bit-rates (up to few terabits per second), which enables a radically different communication paradigm for nano-networks [137]. A. Calibration, Free Space Path Loss (FSPL), and Antenna Cross-Polarization Discrimination A key to all measurements is using a standard approach for calibration, that assures repeatable measurements by any research team at any frequency [138]. 140 GHz FSPL verification measurements were conducted at T-R separation distances of 1, 2, 3, 4, and 5 m using the standard calibration and verification method taught in [138], and the results are shown after removal of antenna gains in Fig. 10, together with 28 GHz and 73 GHz path loss measurement data [2], [139]. The measured path loss at 140 GHz agrees with Friis FSPL equation [140]. The CI path loss model with 1 m reference distance [2] fits perfectly to the measured data, showing that the CI model is viable well above 100 GHz. FIGURE 10. 28, 73 and 140 GHz free space path loss (after subtracting out all antenna gains) verification measurements at distances of 1, 2, 3, 4, and 5 m. View All As shown in Fig. 10, the difference in the path loss measured at 73 GHz and 140 GHz in far field distances within a few meters is 5.85 dB, which is extremely close to the theoretical value calculated by Friis FSPL equation (20×10log1014073=5.66 dB), indicating the high accuracy and proper calibration of the channel sounder system. The difference in the path loss between the 28 and 73 GHz measurements at the same distance is 8.45 dB, which is virtually identical to Friis FSPL equation (20×10log107328=8.32 dB) [138], where calibration and path loss values were calculated using (4)–(6). Pr=G=Pr=PtGtGrλ2(4πd)2Ae4πλ2PtAetAerd2λ2(4)(5)(6) View Source Polarization diversity can provide up to twice the data rates in dual-polarized communication systems, as compared to single co-polarized systems without any extra bandwidth or spatial separation between antennas [141], [142]. Measurements and calibration for antenna cross-polarization discrimination (XPD) at 28, 73, and 140 GHz were conducted to analyze the polarization characterization of the horn antennas [138], [143]. The XPD values also are required to analyze the polarization effects of partitions at different frequencies. The architectures of the channel sounders used in the measurements given here have been extensively described in [28], [121]. The LO, IF frequencies, and antenna characteristics are summarized in Table 3. The output power from the TX waveguide (with the antenna removed) was measured using a PM5-VDI/Erickson Power Meter having a highly accurate power resolution of 10−4 mW. TABLE 3 Summary of NYU Channel Sounder Systems and Antennas Used in Measurements at 28 GHz, 73 GHz and 142 GHz [2], [28], [121] To properly calibrate our measurements, FSPL measurements at 28, 73, and 140 GHz were firstly conducted at 3, 3.5, 4, 4.5, and 5 m, according to the reference procedures introduced in [138], with vertically polarized 27 dBi gain horn antennas at both the TX and RX. The FSPL measurement results after removing the antenna gains are shown in Fig. 8, indicating a negligible difference (less than 0.5 dB) from the theoretical Friis FSPL equation over all distances, which validates the accuracy of the channel sounder. After the free space power measurements with co-polarized antennas, measurements were conducted at the same distances but with cross-polarized antennas (e.g., V-H). Cross-polarization was realized by using a waveguide twist which rotates the antenna by 90°. The insertion loss caused by the twist was measured and calibrated out. The received power for the cross-polarized antennas is shown in Fig. 11. The XPD was calculated by taking the difference between the path losses of the V-V and V-H antenna configurations at same distances. See [138] for a detailed description of the XPD verification and calibration procedure used, that can be universally applied. Note that at a fixed TR separation distance, the free space received powers for the H-H and H-V configurations were within 1 dB of the V-V and V-H received powers, respectively, showing reciprocity with cross-polarization measurements. FIGURE 11. Measured free space antenna XPD at 28, 73, and 140 GHz. The solid lines and the dash lines represent the path loss measured with co-polarized and cross-polarized antennas, respectively. The XPD values calculated across five distances are within 1 dB at each frequency, which verifies the XPD measurement procedure is correct. View All B. mmWave and THz Partition Loss Measurements Wideband mmWave and THz networks, as well as precise ray-tracer algorithms [153], [154], will require accurate channel models that accurately represent the partition losses induced by common building objects [53], [138]. Therefore, partition loss of common building materials needs to be extensively investigated for 5G mmWave wireless systems and future THz wireless communications in and around buildings. Partition loss measurements at 2.5 and 60 GHz of a variety of materials like drywall, office whiteboard, clear glass, mesh glass, using a wideband channel sounder with vertically polarized antennas at both TX and RX were presented in [146], where root mean square (RMS) delay spreads were found to be much lower at 60 GHz than at 2.5 GHz. Based on the measured data, a partition-based path loss model which provided fast and accurate link budget predictions in a multipath-rich environment was developed in [145], [146]. Partition loss measurements at 73 GHz for V-V and V-H polarization configurations for glass doors, drywall, steel doors, and clear glass were conducted in [147]. The co-polarized partition loss for glass doors and windows was found to be 5–7 dB at 73 GHz, and steel doors induced a partition loss as large as 40–50 dB, showing that different materials could be properly chosen for propagation or interference isolation between neighboring rooms [28], [138], [147]. The partition loss of common outdoor and indoor materials are summarized in Table 4 and 5, respectively, and reveal that not much is known above 100 GHz. TABLE 4 Partition Loss of Common Outdoor Materials TABLE 5 Partition Loss of Common Indoor Materials Effective attenuation of typical building materials such as concrete bricks, wood, tiles, and a gypsum plate were measured at 100 GHz with the TX and RX boresight-aligned [152]. The co-polarization situations (both TX and RX were vertically or horizontally polarized) were analyzed, and the effective attenuation of most of the building materials was observed to be polarization sensitive [152]. Work in [149] presented penetration loss measurements of concrete walls and plasterboard walls over the frequency range of 900 MHz to 18 GHz, as shown in Table 5, and the penetration loss was shown to not necessarily increase monotonically with respect to frequency. A comparison of measured partition loss versus frequency for a variety of different common building materials was shown in [150]. As measured in [150], an 8-inch concrete wall had a typical partition loss of 2 to 4 dB at 900 MHz for both V-V and H-H antenna configurations. The attenuation through solid wood, having a thickness of 2 cm, was measured to be about 5 dB and 8 dB at 45 GHz for both H-H and V-V antenna configurations, respectively, while the attenuation through a concrete slab having a thickness of 3 cm was measured to be about 13 dB at 45 GHz for both H-H and V-V antenna configurations [151]. Penetration loss measurements in THz band (0.1 - 10 THz) of plastic board, paper, and glass which had losses of 12.47 dB/cm, 15.82 dB/cm, and 35.99 dB/cm were presented in [148]. However, the measurements in [148] were conducted in a very short range (less than 10 cm) and the thicknesses of sample materials were less than 1 cm which may cause large uncertainty or error. Absorption coefficients of glass, medium density fiberboard, and plexiglass were measured where the S21 parameter was recorded at a T-R separation distance of 10 cm with and without the obstructions brought into the ray path, showing that a 2.5 cm thick window and a 3.5 cm thick fiberboard door would induce absorption attenuation of approximately 65.5 dB and 86.7 dB at 300 GHz [125], respectively. Partition loss measurements at 28, 73, and 140 GHz were conducted at the NYU WIRELESS research center. T-R separation distances of 3, 3.5, 4, 4.5, and 5 m were used and the TX/RX antenna heights were 1.6 m (see Fig. 3 in [138]). These distances are greater than 5×Df , where Df is the Fraunhofer distance to ensure a plane wave is incident upon the measured material, and the dimensions of the measured materials are large enough such that the radiating wavefront from the TX antenna is illuminated on the material without exceeding the projected HPBW angle spread from the TX antenna [138]. At each distance, five measurements/power delay profiles (PDP) were recorded with slight movements on the order of half a wavelength, taking the average of the power in the first arriving multipath component of the recorded PDPs, to exclude the multipath constructive or destructive effects. With 4 GHz RF bandwidth, the channel sounder systems used in this paper have a time resolution of 0.5 ns, which means any multipath components spaced more than 0.15 m in propagation distance can be resolved [28]. The partition loss results of common building construction materials, drywall and clear glass, at 28, 73 and 140 GHz, were listed in Table 6 and Table 7. The mean partition loss of clear glass at 28 GHz for co-polarized situation (V-V and H-H), see Table 6, was about 1.5 dB with a standard deviation (std) of 0.5 dB. But the same glass had greater partition losses of 7.17 dB and 10.22 dB at 73 GHz and 140 GHz, respectively, for the co-polarization situation. TABLE 6 Partition Loss Measurement Results at 28, 73, and 140 GHz for Clear Glass TABLE 7 Partition Loss Measurement Results at 28, 73, and 140 GHz for Drywall It is worth noting that the antenna free space XPD has not been subtracted from the cross-polarized partition loss measurements shown in Fig. 12, Fig. 13, Table 6 and Table 7. Subtracting the XPD results in a negative value of partition loss due to the polarization coupling effects (depolarization) of the building materials. It is shown in Fig. 12 and 13 that the partition loss tends to increase with frequencies for both co-polarized and cross-polarized antenna configurations. FIGURE 12. Partition loss measurement results of clear glass with a thickness of 1.2 cm at 28, 73, and 140 GHz. Partition loss trends to increase with frequencies for both co-polarized and cross-polarized antenna configurations. View All FIGURE 13. Partition loss measurement results of drywall with a thickness of 14 cm at 28, 73, and 140 GHz. View All FIGURE 14. Setup for scattering measurements. θi is the angle the incident ray makes with the normal and θs is the scattered angle. View All FIGURE 15. Photograph of the scattering measurement setup. The scattered power of drywall was measured at a distance 1.5 m away from the wall, in angular increments of 10°. View All SECTION VIII.Scattering Above 100 GHz Scattering at mmWave and THz frequencies is an important propagation mechanism [6], [62]. The surfaces of buildings, terrain, walls, and ceilings have usually been assumed to be electrically smooth since their surface height variations are small compared to the carrier wavelength at lower frequencies (e.g., greater than 5 cm at below 6 GHz) [155]. At today’s sub-6 GHz cellular and WLAN systems, the reflection process is dominated by a strong specular path at an angle of reflection equal to the angle of incidence, and scattering, being a weaker propagation phenomenon, is negligible. However, in mmWave bands and above, the roughness of surfaces are comparable to the carrier wavelength and the illuminated scatterers may actually create signal paths that are as substantial as, or even stronger than reflected paths, depending on the incident angle [6], [62], [155]. A. Directive Scattering Model The directive scattering (DS) model is widely used to predict the scattered power in optics [156]. The single lobe DS model [157], [158] assumes that the main scattering lobe is steered in the direction of the specular reflection (θr shown in Fig. 16). The DS model has been used in [158] to model the RF propagation environment of a hospital room at 60 GHz. The PDP of the environment agreed well with simulations using the DS model, up to an excess delay of 30 ns. The DS model has also been tested at 1.296 GHz in [157], where the DS model agreed with the scattering from rural and suburban buildings. When an electromagnetic wave impinges upon a surface at an incident angle θi , the scattered electric field at any particular scattering angle θs can be calculated using the DS model. The DS scattered electric field in the incident plane is given by: |Es|2==|Es0|2⋅(1+cos(Ψ)2)αR(SKdtdr)2lcosθiFαR⋅(1+cos(Ψ)2)αR(7) View SourceRight-click on figure for MathML and additional features.where Es is the scattered electric field at the scattering angle Ψ . Es0 is the maximum scattered electric field, which is adopted from an effective roughness model [157]. K is given by 60 PtGt‾‾‾‾‾‾‾√ which is a constant depending on the transmitted power and the transmitter antenna gain [159]. dt and dr are the distances between the scatterer and the transmitter and receiver, respectively. l is the length of the scattering object. Ψ is the angle between the reflected wave and the scattered wave, as shown in Fig. 16. αR determines the width of the scattering lobe with higher values of αR implying a narrower scattering lobe. FIGURE 16. - A radio wave incident at an angle $\theta _{i}$ with respect to the normal to the targeted rough surface. $\theta _{r}$ and $\theta _{s}$ are reflected and scattered angle, respectively. From Snell’s law, reflections obey $\theta _{i} = \theta _{r}$ . $\Psi $ is the angle between reflected and scattered waves. FIGURE 16. A radio wave incident at an angle θi with respect to the normal to the targeted rough surface. θr and θs are reflected and scattered angle, respectively. From Snell’s law, reflections obey θi=θr . Ψ is the angle between reflected and scattered waves. View All The received power at the receiver can then be calculated as [118]: Pr=PdAe=|Es|2120π⋅Grλ24π=|Es|2Grλ2480π2(8) View SourceRight-click on figure for MathML and additional features.where Pd is the power flux density of the scattered wave, and Ae is receiver antenna aperture [118]. Gr is the receiver antenna gain, and λ is the wavelength of the radio wave. The scattered power is calculated using (7)–(8), and the reflected power is calculated using: Γrough=ρs⋅Γsmooth,(9) View SourceRight-click on figure for MathML and additional features.where ρs is the scattering loss factor and ρs=exp[−8(πhrmscosθiλ)2]I0[8(πhrmscosθiλ)],(10) View SourceRight-click on figure for MathML and additional features.where I0 is the zero-order Bessel function of the first kind. When θi=1∘ , 30°, and 45° (small incident angles), it was shown in [62] that the maximum scattered power (in the reflection direction as in Snell’s law) computed by the DS model is stronger than the reflected power obtained from the rough surface reflection model for a rough surface. Simulations were performed using three materials with incident angle θi ranging from 10° to 90° and over frequencies from 1 GHz to 1 THz [62]. It was shown that as the frequency increases, the received scattered power increases, which can be expected since the surfaces tend to be rough as the frequency increases and rougher surfaces cause greater scattered power (even greater than the reflected power). Moreover, the maximum scattered power is received when the incident wave impinges upon the surface along the normal direction. The scattered power falls off sharply when the incident wave becomes grazing, and most of the incident power is reflected [62]. The single lobe DS model may be modified to incorporate power an additional back-scatter lobe [157]. The DS scattered electric field taking into account for backscattering is given by: |Es|2=|Es0|2[Λ⋅(1+cos(Ψ)2)αR+(1−Λ)⋅(1+cos(Ψi)2)αR](11) View Sourcewhere Ψi is the angle between the scattered ray and the incident ray and Λ determines the relative strength of the back scattered lobe with respect to the main scattered lobe. B. Scattering Measurements and Results A smooth metal plate or a conformal metal foil attached on the wall (having a size much larger than the radius of first Fresnel zone for frequencies of 100, 200, 300, and 400 GHz) was shown to provide 6–10 dB more power from specular reflection (at an incident angle of 50°) than the bare painted cinderblock wall did [104], indicating that the effect of scattering from the rough surface was significantly smaller than the effect of absorption from the bare painted cinderblock wall, even at 400 GHz. Specular reflection loss with varied incident angles measurements at 400 GHz were also conducted in [104] with 5 different incident angles ranging from 20° to 60° which showed that when the incident angle is large, reflection loss is small (absorption loss plus the scattering loss) and the scattering loss is negligible. However, when the incident angle is small (e.g., incident ray impinges the wall perpendicularly), the scattering loss is not negligible (4 dB difference than predicted result when scattering losses are neglected) [104]. Scattering measurements at 60 GHz were presented in [155] using a 2 GHz wide baseband sounder to measure the non-specular diffuse scattering with an incident angle of 15°, 30°, and 45°. Two types of building materials were used: a red stone wall (a rough wall surface) and a concrete pillar (a smooth wall surface). The received power was measured in the range of 0° to 90° from the normal to the wall and the term of power concentration was defined as the angular span corresponding to 90% of the received power. It was shown that the angular span is centered at the incident angle and higher incident angles lead to smaller angular spread [155]. The power concentration of the rough wall surface was shown to be significantly larger than the case of a smooth wall surface and the coupling effects (depolarization) is more severe when reflecting/scattering from rough wall surface [155]. Scattering measurements of a drywall were conducted at 142 GHz using the channel sounder system summarized in Table 3, and the measurement setups and photos are shown in Fig. 14 and 15. During the scattering measurements, both the heights of TX and RX were set at 1.2 m on an arc with a radius of 1.5 m. Incident angles of 10°, 30°, 60°, 80° were chosen to study the scattering performance from a small incident angle to a large angle with respect to the normal. The RX power was measured in the same plane of the incident ray and the normal of the surface from 10° to 170° in 10°/step (0° and 180° are not able to measure due to the physical size of the channel sounder system). The path loss is calculated by: PL[dB]=PTX[dBm]−PRX[dBm]+GTX[dBi]+GRX[dBi],(12) View Sourcewhere PTX [dBm] is the true transmit power in dBm, PRX [dBm] is the record received power by the channel sounder system, GTX [dBi] and GRX [dBi] are the TX and RX antenna gain in dB, respectively. The scattered loss (SL) was then normalized with respect to the FSPL: SL=FSPL(f,d)[dB]=PL[dB]−FSPL(f,d)[dB],32.4+20log10(f)+20log10d,(13)(14) View SourceRight-click on figure for MathML and additional features.where f is the carrier frequency in GHz and d is the TR separation distance in meter. Measured scatter patterns of different incident angles at 142 GHz are shown in Fig. 17. The peak scattered power (scattered plus reflected) was observed at the specular reflected angle (Snell’s law), and peak scattered power was greater at large incident angles than at small incident angles (9.4 dB difference between 80° and 10°, where most of the energy is due to reflection and not scattering). At all angles of incidence, measured power is within 10 dB below the peak power in a ±10° angle range of the specular reflection angle, likely a function of antenna patterns. In addition, back scattered are also observed (negative RX angles) when the incident angle is small (e.g., 30° and 10°) but was more than 20 dB below the peak received power,which means that the surface of drywall can be modeled as smooth surface at 140 GHz and the reflected power can be used to estimate the received power in the specular reflection direction, especially close to grazing. FIGURE 17. - Polar plots of the scattered power off drywall, measured at four incident angles. Power at scattered angles is measured in 10° increments. FIGURE 17. Polar plots of the scattered power off drywall, measured at four incident angles. Power at scattered angles is measured in 10° increments. View All A comparison between measurements and predictions made by the dual-lobe DS model (11) with TX incident angle θi=10∘ , 30°, 60°, and 80° is shown in Fig. 18. It is shown that the peak received power of the simulation at the specular reflection angle agrees well with the measured data (within 2 dB). The received power (measured and simulated) changes versus the RX pointing angle (θr as shown in Fig. 14). A back-scattered peak is observed in both the measurement data and the dual-lobe DS model prediction for all the incident angles, and the back-scattered power could be used for imaging around the corner [50], [58]. FIGURE 18. Comparison between measured data and the dual-lobe Directive Scattering model with TX incident angle θi=80∘,60∘,30∘ , and 10° at 142 GHz. (a) θi=80∘ . (b) θi=60∘ . (c) θi=30∘ . (d) θi=10∘ . View All SECTION IX.User Localization in Millimeter Wave and Terahertz Bands Positioning is the determination of the location of a user that is fixed or moving, based on the known locations of other base stations (BSs). Although users can be localized to within 4.9 m under open sky [160] by the global positioning system (GPS), in obstructed environments such as urban canyons and indoors, the positioning accuracy of GPS is poor. In such environments, 3D maps of the environment along with accurate ToF and angle of arrival (AoA) measurements could be leveraged for precise positioning. A. Time of Flight (ToF) and Phase-Based Positioning In LOS environments, the ToF of a signal can be used to estimate the TR separation distance d (since d=c⋅t where c is the speed of light and t is the ToF). The RX’s position may then be estimated by trilateration. In NLOS environments, ranging introduces a positive bias in the position estimates since the path length of reflected multipath rays is longer than the true distance between the user and the BS. The ultra-wide bandwidths available at THz frequencies, allow the RXs to resolve finely spaced multipath components and accurately measure the signal’s ToF. In [161], the phase of the received signal in LOS was used to estimate the TR separation distance at 300 GHz. By tracking the baseband modulated signal’s phase at the RX and manually correcting for phase ambiguity, decimeter-level accuracy was achieved up to a distance of 40 m [161]. To get the absolute ToF, synchronization is required between the TX and RX, which is not easy in a large network of nodes widely separated. A timing error of 1 ns may result in a positioning error of 30 cm, if no correction for clock bias or clock drift is applied. It is easier for the RX to measure the time-difference-of-arrival (TDoA) of signals arriving from multiple synchronized TXs. In order to determine the relative time delays between the signals, cross-correlation is done. The peaks in the cross-correlation function correspond to the signal propagation time delay. B. Angle Based Positioning The line of bearing to the BS is given by the AoA of the signal at the user. The user is localized to the point of intersection of the lines of bearing from several BSs, as seen in Fig. 19. In NLOS environments, due to specular reflections off walls and metallic surfaces, rays do not arrive from the direction of the BS, leading to accuracy penalties [53]. FIGURE 19. - The UE measures the relative AoAs, $\theta _{1} $ and $\theta _{2} $ from $BS_{1}$ , $BS_{2}$ , and $BS_{3} $ . Narrow antenna beamwidths at THz frequencies facilitate accurate AoA measurements. FIGURE 19. The UE measures the relative AoAs, θ1 and θ2 from BS1 , BS2 , and BS3 . Narrow antenna beamwidths at THz frequencies facilitate accurate AoA measurements. View All FIGURE 20. - Three multipath components arrive at the user (RX) shown above - one LOS component (in green) and two NLOS components (in blue and black). Of the six candidate locations for the user (RX), based on AoD and ToF measurements at the BS (TX) ( $\text {CL}_{1} - \text {CL}_{6} $ ), three candidate locations $(\text {CL}_{3}, \text {CL}_{5}, \text {CL}_{6}) $ correspond to the actual location of the user. The position of the user is estimated to be the modal candidate location (i.e. $\text {CL}_{3}, \text {CL}_{5}, \text {CL}_{6} $ ). FIGURE 20. Three multipath components arrive at the user (RX) shown above - one LOS component (in green) and two NLOS components (in blue and black). Of the six candidate locations for the user (RX), based on AoD and ToF measurements at the BS (TX) (CL1−CL6 ), three candidate locations (CL3,CL5,CL6) correspond to the actual location of the user. The position of the user is estimated to be the modal candidate location (i.e. CL3,CL5,CL6 ). View All FIGURE 21. - The CDF of the positioning error, when an RX is localized using a single TX. The median positioning error is 6.7 cm when $\sigma _{ToF} $ and $\sigma _{AoD} $ were 0.25 ns and 0.5° respectively. FIGURE 21. The CDF of the positioning error, when an RX is localized using a single TX. The median positioning error is 6.7 cm when σToF and σAoD were 0.25 ns and 0.5° respectively. View All The AoA of the signal is found from detecting the maximum amplitude of the received signal as a function of angular direction. Many methods, such as a sum and difference method used in radars or optical receivers [162] or peak power method may be used. In [53], [163], all possible directions were scanned by a rotating RX antenna and the angle at which maximal power was received was determined to be the AoA. Alternatively, an antenna array may be used to estimate the AoA, by beamforming and subspace based methods [164]. The Cramer-Rao lower bound (CRLB) is a bound on the minimum variance of an unbiased estimate of a parameter. An increase in the number of antenna array elements leads to greater angular resolution and thus a reduction in the CRLB for estimating AoA [165]. C. NYURay – A 3D mmWave Ray Tracer Ray tracing can be done to simulate channel characteristics at different locations. By using a ray tracer that replicates measured data with high accuracy, costly and time-intensive measurement campaigns can be avoided. Ray tracing allows researchers to test positioning algorithms. In addition, ray tracers can be used to create a radio map of an environment, which could replace the off-line stage of fingerprinting algorithms. Building on the 2-D mmWave ray tracer developed in [53], NYURay, a 3-D mmWave ray tracer has been developed. NYURay is a hybrid ray tracer which combines shooting-bouncing ray (SBR) based ray tracing with geometry based ray tracing [58]. D. Map-Based Localization Work in [50] demonstrated how ToF, AoD measurements, along with a map of the obstacles encountered by the signal could be used to localize a user as described in Section II-D. However, [50] assumes that all materials act like mirrors. The measurements described in Section VII illustrate that at 140 GHz, the partition loss through glass and dry wall is close to 10dB, indicating that materials do not act like perfect reflectors at mmWave and THz frequencies. The authors of [166] used ToF and an indoor map for localization, using the principal of “virtual access points” to get additional anchor points to localize users that are in NLOS. However [166] does not use AoA information, which will be available with high accuracy using directional antennas at THz frequencies. In [167], authors used backwards ray tracing using AoA and power observations to estimate the TX location. The transmitted signal had a power of 30 dBm and a bandwidth of 300 MHz. A thermal noise floor of −85 dBm was assumed. Assuming perfect AoA resolution, TXs were localized to within 5 m, 90 % of the time, using six RX locations. E. Map-Assisted Localization With a Single BS Map-assisted positioning with angle and time (MAP-AT) is a novel map-assisted localization technique that shall now be presented. MAP-AT may be used when at least two multipath components arrive at the user. The user need not be in LOS of the BS. The user may either be the TX or the RX of a known radio signal. If the user is a TX, ToF and AoA information is required. If the user is an RX, ToF and AoD information is required. In combination with a 3-D map of the area, the location of the user can be determined. First consider the case where the user is an RX. If the ToF and AoD of the ray that reaches the user is known and if the ray reaches the user either directly from the BS or through one reflection or by passing through one obstruction, there are two possible locations of the user. If the ray reached the user after one reflection, the BS and user must lie on the same side of the reflecting object. If the ray reached the user directly from the BS, or through one obstruction, the BS and user must lie on a straight line joining the BS and user. The possible locations of the user, based on ToF and AoD measurements at the BS shall henceforth be referred candidate locations. The process of finding candidate locations is repeated for all multipath components. The actual user location corresponds to the candidate location identified by the maximum number of multipath components. Errors in ToF measurements cause the BS to incorrectly estimate the path length to the user. Due to inaccurate AoD measurements the BS incorrectly estimate the user’s bearing. As a result, a candidate locations estimated by MAP-AT using imprecise ToF and AoD information will not coincide with the user’s true location. However, it is likely that the candidate locations will be close to the user’s true location, and hence close to one-another. MAP-AT is modified to first group the candidate locations that are close to one another (at a distance of up to 40 cm) and then estimate the user position to be the centroid of the candidate location group with maximal members. If the user is a TX, a similar technique can be used for user localization, using AoA information at the BS and doing backwards ray tracing instead of forwards ray tracing. 1) Simulation Environment Simulations of localization using candidate locations were done by synthesizing ToF and AoD measurements at 30 BS-user combinations using NYURay, of which 18 were in NLOS and 12 were in LOS. The BS and user locations were taken from previous indoor propagation measurement campaigns conducted in the NYU WIRELESS research center on the 9th floor of 2 MetroTech Center in downtown Brooklyn, New York [168]. The research center has an open-office architecture, with cubicles, walls made of drywalls and windows. The entire floor has dimensions of 35 m by 65.5 m as shown in Fig. 22. FIGURE 22. - The three error circles centered at each user location correspond to the rms error over 100 simulations for three noise levels: $\sigma _{ToF} = 0.25$ , 0.5 and 1 ns with $\sigma _{AoD} = 0.5^\circ $ in the azimuth and elevation plane. A single BS was used to localize each of the 30 user locations. FIGURE 22. The three error circles centered at each user location correspond to the rms error over 100 simulations for three noise levels: σToF=0.25 , 0.5 and 1 ns with σAoD=0.5∘ in the azimuth and elevation plane. A single BS was used to localize each of the 30 user locations. View All The position of each user was determined using a single BS. To make the simulations realistic, zero mean Gaussian noise with standard deviation σToF and σAoD was added to the ToF and AoD measurements. The positioning error for each RX was defined to be equal to the Euclidean distance between the position estimate and the true position of the RX. 2) Localization Error Analysis The cumulative distribution function (CDF) for the root mean squared (rms) positioning error in presence of Gaussian noise over 100 simulation runs at each user location is plotted in Fig. 21. When σAoD=0.5∘ and σToF=0.25 ns, the median positioning error was 6.7 cm. Three concentric error circles of radii equal to the rms positioning error corresponding to the three levels of noise added to ToF measurements have been drawn centered at each user location in Fig. 22. The error circles for two outlier user locations having errors of 4.0 m and 4.9 m have not been plotted. The localization error of the outliers can be reduced by using additional BSs for positioning the users [58]. Table 8 illustrates how the localization error varies with TR separation distance in LOS and NLOS environments. TABLE 8 Performance of the Map-Assisted Localization Algorithm for Different TR Separation Distances in LOS and NLOS Environments at NYU Wireless Table 8 SECTION X.Spatial Statistical Channel Modeling THz frequency poses new challenges for channel modeling. Diffraction effect at THz frequencies will not be as prominent as it is at microwave frequencies. Scattering, instead, is expected to be more notable since the wavelength at THz frequencies are comparable to the size of dust, rain, snow, and roughness of walls [113]. Thus, the THz channel will be even more sparse than the mmWave channels, and LOS and specular reflection may be the only two reliable propagation paths. A. Spatial Consistency Spatial consistency, an important channel modeling component, has been acknowledged by 3GPP standards [169]. Most of current statistical channel models are drop-based, which are only able to generate channel impulse responses (CIRs) for a particular user at a randomly picked location in a distance range [170]. There is no spatial correlation between two consecutive simulation runs. Thus, it is impossible to generate dynamic CIRs with spatial correlation based on a user’s motion within a local area [73], [171], [172]. Spatial consistency enables channel models to provide spatially consistent and smoothly time-evolving CIRs for closely spaced locations. B. Early Research on Spatial Consistency Even though spatial consistency is a relatively new channel modeling component, some valuable works have been done not only for channel models but also for channel simulators. 3GPP Release 14 proposed a three-dimensional (3D) geometry-based stochastic channel model (GSCM) and also provided alternative spatial consistency procedures to enable spatially consistent mobility simulations [169]. 5GCM proposed three spatial consistency approaches, which the first one uses spatially correlated random variables as small-scale parameters, and the other two are based on the geometry and fixed locations of scatterers [173]. COST 2100 introduced spatial consistency based on the visible region over which a group of MPCs is visible at a user [174]. QuaDRiGa proposed a time-evolving channel model for user trajectory based on the WINNER II model using geometrical approaches [175]. NYUSIM channel model, a spatial statistical channel model, has also been enabled spatial consistency by generating spatial-correlated large-scale parameters and time-variant small-scale parameters [170]. Some local area measurements in a street canyon scenario have been conducted, and the results indicated that the correlated distance of the large-scale parameters such as the number of time clusters is about 10–15 meters [176]. C. Challenges for Spatial Consistency at mmWAVE and THz Frequencies The Doppler shift will be more severe at THz frequencies, which is 90 kHz at 1 THz with a speed of 60 mph. Further, as the channel coherence time decreases, the channel status will change fast and needs to be updated more frequently. Thus, it is critical to provide consecutive and realistic CIRs as time evolves, which can be provided by spatial consistency. For super resolution imaging, understanding the minor differences of the channel over a few wavelengths, e.g., over hundreds of wavelengths (hundred of micrometers), will require new measurements and models to determine the tiny variations. Statistical modeling for sensing and imaging at THz will require much greater spatial resolution than ever used before in channel modeling. Most current research on spatial consistency is based on the assumption of omnidirectional antennas at the transmitter and receiver [175], [177]. However, considering that high-gain steerable antennas and beamforming will be used at mmWave and THz frequencies, spatial consistency must be studied and implemented with directionality. When a narrowbeam antenna points to a certain direction, just a few scatterers are illuminated and can generate MPCs. Thus, when the user is moving, the part of environment that can be “seen” will change fast. Spatial consistency needs to be maintained with advanced beam tracking schemes to provide accurate CIRs with time evolution. Many challenges have not been resolved yet before the actual implementation of spatial consistency for narrowbeam directional antennas. The correlation distances of different parameters are abstracted from field measurements, but most measurements only provide reference correlation distance of spatial consistency for omnidirectional antennas. The antenna beamwidth and the ability of beam steering will have an impact on the value of correlation distance, which has not been investigated but will be at much smaller distance scales than ever used before. In addition, more user cases (scenarios) are introduced in the latest 3GPP Release 15 [5] for 5G than the previous generations of cellular systems, and even more user cases will surely evolve into 6G and beyond. Therefore, scenario-specific values of several key channel model parameters will be needed in the future. SECTION XI.Conclusion This paper has provided a broad and deep look at fundamental opportunities, challenges, and approaches to creating future wireless, sensing, and position location systems for spectrum above 100 GHz, which are likely to be a part of the 6G era in the 2025–2035 time frame. The paper described recent regulatory and standard body activities aimed at fostering future wireless systems that will use multi-GHz bandwidth channels at above 100 GHz, and which may support data rates in excess of 100 Gbps. We demonstrated a large number of promising applications that future THz frequencies will support, as computing power concurrently grows to approach the processing power of the human brain. Many applications of THz wireless will enable novel cognition, sensing, imaging, communications, and positioning capabilities that will be used by automated machinery, autonomous cars, and new human interfaces, all enabled by the ultra-wide bandwidth and ultra-short wavelength at THz which appears to be a promising spectrum for future wireless communications beyond the mmWave regime. We showed how directional steerable antennas will enable mobile communications well into the THz band, where the combined antenna gains can overcome the atmospheric and weather losses that have been believed to be prohibitive in the past. Antenna array technologies will need to exploit new approaches and physical architectures, such as spatial noise-shaping, beamforming, hybrid beamforming, and the cones of silence, as these were shown to provide significant performance benefits and will solve design constraints such as the physical size of RF components with a vast number of antenna elements. We also indicated how simultaneous imaging and sensing with location capabilities may be enabled by the move to above 100 GHz. This paper surveyed past work and presented new propagation measurements above 100 GHz, including a comparison of cross-polarization discrimination (XPD) and partition loss measurements for common building materials at 28, 73, and 140 GHz. We also derived scattering theory results that are useful for understanding how the frequency and surface roughness influence scattering behavior across the mmWave and THz bands, and provided reflection/scattering measurements at 140 GHz that validated the Directive Scattering (DS) theory. We noted that extensive measurements in urban, suburban, and rural areas are necessary to study the channel characteristics using statistical and deterministic approaches at frequencies above 100 GHz. It is clear from this work that THz positioning will support centimeter level accuracy and may also support imaging, even in non line of sight (NLOS) environments. Challenges such as power-efficient devices, cost-effective integrated circuit solutions, and practical phased arrays that may be interconnected with minimal loss loom as impediments to 6G and THz product development, and offer open research and development problems that are being investigated by DARPA and other global research agencies. Other challenges remain, such as how to efficiently construct and code maps of the environment, how to model the THz channel for spatial consistency, and how to reduce computational complexity in spatial multiplexing and beam codebooks for massive numbers of antenna elements, and this paper gave some insights into possible solutions in these areas. Certainly not least, as the world begins to look to 6G and beyond, at frequencies above 100 GHz, there should be careful studies to determine the impact that THz radiation has on the biological effects of humans and animals, and to see if such radiation could be used as a safe and effective replacement for more dangerous ionizing imaging methods such as X-rays and CT/PET scans, provided losses due to water and tissue are not too great. Authors Figures Citations Keywords Metrics More Like This A hybrid adaptive antenna array IEEE Transactions on Wireless Communications Published: 2010 Design of an adaptive antenna array for tracking the source of maximum power and its application to CDMA mobile communications IEEE Transactions on Antennas and Propagation Published: 1997 View More Top Organizations with Patents on Technologies Mentioned in This Article References References is not available for this document. IEEE Personal Account CHANGE USERNAME/PASSWORD Purchase Details PAYMENT OPTIONS VIEW PURCHASED DOCUMENTS Profile Information COMMUNICATIONS PREFERENCES PROFESSION AND EDUCATION TECHNICAL INTERESTS Need Help? US & CANADA: +1 800 678 4333 WORLDWIDE: +1 732 981 0060 CONTACT & SUPPORT Follow About IEEE Xplore | Contact Us | Help | Accessibility | Terms of Use | Nondiscrimination Policy | Sitemap | Privacy & Opting Out of Cookies A not-for-profit organization, IEEE is the world's largest technical professional organization dedicated to advancing technology for the benefit of humanity. © Copyright 2020 IEEE - All rights reserved. Use of this web site signifies your agreement to the terms and conditions.