Hexa-X-D2 3 v1.0
Hexa-X-D2 3 v1.0
Hexa-X-D2 3 v1.0
Project Name:
A flagship for B5G/6G vision and intelligent fabric of technology enablers connecting human,
physical, and digital worlds
Hexa-X
Deliverable D2.3
Radio models and enabling techniques
towards ultra-high data rate links and
capacity in 6G
Date of delivery: 31/03/2023 Version: 1.0
Start date of project: 01/01/2021 Duration: 30 months
Document properties:
Document Number: D2.3
Document Title: Radio models and enabling techniques towards ultra-high data
rate links and capacity in 6G
Editor(s): Ahmad Nimr (TUD)
Authors: Ahmad Nimr (TUD), Marko E. Leinonen (OUL), Muris
Sarajlic (EAB), Tommy Svensson (CHA), Katsuyuki Haneda
(AAU), Andreas Wolfgang (QRT), Han Yu (CHA), Charitha
Madapatha (CHA), Hao Guo (CHA), Pål Frenger (EAB),
Rafael Puerta (EAB), Behrooz Makki (EAB), Yinggang Li
(EAB), Omer Haliloglu (EBY), Kilian Roth (INT), Leonardo
Gomes Baltar (INT), Hardy Halbauer (NOG), Le Hang Nguyen
(NOG), Mohamed Abdelbasset Aliouane (ORA), Jean-Marc
Conrat (ORA), Pekka Kyösti (OUL), Nuutti Tervo (OUL),
Aarno Pärssinen (OUL), Xinxin Yang (QRT), Meik
Dörpinghaus (TUD), Stephan Zeitz (TUD), Roberto Bomfin
(TUD)
Contractual Date of Delivery: 31/03/2023
Dissemination level: PU1/
Status: Final
Version: 1.0
File Name: Hexa-X D2.3
Revision History
This report presents the findings of studies and performance analysis on 6G radio system for achieving
ultra-high data rate links and capacity in the context of Hexa-X work package 2: “Novel radio access
technologies towards 6G”. The report begins by presenting an overview of use cases relevant to sub-
THz communication, defining representative scenarios, technical requirements, and performance
metrics. Then, technical studies follow in the areas of channel modelling, radio architecture including
transceiver design and distributed MIMO, in addition to signal processing and transmission schemes.
In particular, radio transceiver design approaches are explored considering hardware non-idealities,
channel properties, and waveform impacts. The report introduces description and analysis based on
measurement for different hardware components and channel model using measurement data, which
are used in the design and analysis of the physical layer modulation and waveforms. Furthermore, D-
MIMO architecture options and corresponding transmission and processing schemes are investigated
with performance analysis focusing on coverage and throughput. Finally, system-level simulation is
conducted to illustrate the impact of different deployment options and design parameters on the
performance.
Keywords
100 – 300 GHz, hardware model, beamforming, channel measurements, distributed MIMO, physical
layer, radio model, simulation model
Disclaimer
The information and views set out in this deliverable are those of the author(s) and do not necessarily
reflect views of the whole Hexa-X Consortium, nor the official opinion of the European Union.
Neither the European Union institutions and bodies nor any person acting on their behalf may be held
responsible for the use which may be made of the information contained therein
This project has received funding from the European Union’s Horizon 2020 research and
innovation programme under grant agreement No 101015956.
Hexa-X Deliverable D2.3
Executive Summary
This report is the third deliverable of Hexa-X project work package 2 (WP2): “Novel radio access
technologies towards 6G”. It focuses on 6G radio system design aspects and enabling techniques for
fulfilling the requirements of emerging use cases, especially, link data rate of 100 Gbps. These include
radio frequency (RF) transceiver architectures and hardware (HW) components models for frequency
above 100 GHz, channel measurements and modelling, deployment scenarios with cellular-based and
distributed MIMO (D-MIMO) approaches, and signal processing topics on waveform, beam
management, and D-MIMO transmission schemes. Moreover, performance analysis is presented with
different objectives, such as link-level evaluation of modulation under HW nonidealities, and system-
level analysis to assess the performance of different design and deployment options.
The report starts with an overview of a set of use cases enabled by sub-THz communication in the band
(100 GHz-300 GHz) introduced in deliverable [HEX21-D21] with key performance indicators (KPIs)
requirements for data rate and latency. Based on range requirements, three scenarios (mid-range, short
range, and very short range) are defined with detailed specifications including the positioning
requirements from WP3 report [HEX21-D31]. For each scenario, several radio options are foreseen
depending on the operating band and device class (user equipment (UE), or access point (AP)). To
narrow the design space, system performance metrics are used to evaluate the design and deployment
scenarios, including coverage, energy performance, link reliability and latency, beam failure rate, and
positioning metrics, for the. Moreover, design considers the impacts of radio channel, including path
loss, fading, and opportunities for spatial beams. This is based on the results of channel models for
material interaction in the range (2-260 GHz), and a using the stored channel model data at 140 GHz.
The report then presents technical guidelines for RF transceiver design to achieve link data rate at a
particular range with detailed discussions and evaluation examples for each step. Various parameters
and degrees of freedom are discussed, including bandwidth, signal-to-noise ratio (SNR), analogue-to-
digital convertors (ADC), HW nonidealities, and waveform impact on the power amplifier (PA)
backoff. Moreover, it discusses RF implementation aspects and strategies to optimize the system such
as channelization and carrier aggregation (CA). Based on the analysis, it is concluded that sub-array-
based RF transceiver is a practically feasible architecture for above 100 GHz, which implies analogue
beam steering per subarray and digital precoding over multiple RF chains. Different configurations for
the RF chains are possible; e.g., spatial multiplexing, or implementation of ultra-wideband waveforms
following CA approach. However, conventional initial beam access is impractical due to the large
number of beams to be probed, and it is necessary to consider side information for beam management.
The document describes several hardware models for local oscillator (LO) phase noise and PA, and
their accuracy is justified based on measurements and circuit analysis. These models are considered in
the signal processing analysis to design waveforms and modulation schemes that mitigate HW non-
idealities, or to develop solutions for estimation and compensation. Accordingly, performance analysis
is presented for one or more HW impairments, with special focus on DFTS-OFDM and single carrier
(SC), as they show tolerance to PA nonidealities compared to OFDM [HEX21-D22], and zero-crossing
modulation (ZXM) with 1-bit quantization as a potential energy-efficient modulation scheme candidate.
With respect to other performance metrics, a study on RAN latency to achieve 0.1 ms is conducted,
showing that processing latency with high throughput is the most critical factor. This needs to be
considered in the design of processing architecture considering I/O speeds, memory access speed, and
power consumption. The report also provides several performance analyses in the context of D-MIMO
architecture, focusing on spectral efficiency, throughput, and coverage for different transmission
schemes and techniques to cope with blockage. Finally, a system-level simulation example is studied
to illustrate the impact of radio configuration and deployment options on the performance in terms of
area throughput and energy consumption, showing the potential for optimizing radio design in terms of
transceiver architecture and deployment scenario to achieve KPIs at lower energy cost.
Table of Contents
List of Figures...............................................................................................................................7
List of Tables ..............................................................................................................................11
List of Acronyms ........................................................................................................................13
1 Introduction .........................................................................................................................18
1.1 Objective of the document............................................................................................18
1.2 Structure of the document ............................................................................................18
2 Scenarios and technical requirements ...............................................................................20
2.1 Representative communication scenarios for sub-THz ................................................20
2.1.1 Relevant use cases ...................................................................................................20
2.1.2 Sub-THz scenarios ..................................................................................................21
2.2 Flexible radio hardware design.....................................................................................22
2.2.1 Radio hardware options ...........................................................................................22
2.2.2 Radio hardware design methodology ......................................................................23
2.3 Performance metrics .....................................................................................................25
3 Channel model .....................................................................................................................26
3.1.1 Theoretical model and ITU model ..........................................................................26
3.1.2 Measurement description and data processing ........................................................27
3.1.3 Conductivity/permittivity estimation.......................................................................28
3.1.4 Comparison between measurement and the ITUF model .......................................30
3.1.5 ITU model validation above 100 GHz ....................................................................32
3.2 Stored channel model at 140 GHz ................................................................................33
3.2.1 Principle...................................................................................................................33
3.2.2 Stored channels from measurements .......................................................................34
3.2.3 Channel generation ..................................................................................................35
4 Radio architecture and models...........................................................................................37
4.1 RF transceiver architecture modelling methods ...........................................................37
4.1.1 Modelling from signals to RF performance ............................................................39
4.1.1.1 Data rate .............................................................................................................40
4.1.1.2 Sampling rate and resolution requirements of data converters ADC & DAC....42
4.1.1.3 ADC analysis based on the requirements ...........................................................43
4.1.1.4 Output power vs. Psat (backoff) and waveform .................................................44
4.1.1.5 Minimum SNR/EVM for radio link ...................................................................45
4.1.2 Link range analysis using 100 Gbps as an example ................................................45
4.1.2.1 General considerations on the use of the model .................................................45
4.1.2.2 Modelling example for 100Gbps communications ............................................47
4.2 Hardware models ..........................................................................................................49
4.2.1 Phase noise ..............................................................................................................49
4.2.1.1 Millimetre-wave and sub-THz phase noise for frequency-multiplied LO chains49
4.2.1.2 System phase noise behaviour and limitations on the achievable performance and
bandwidth ...........................................................................................................51
4.2.2 Power amplifier modelling ......................................................................................55
4.2.2.1 Frequency scalable saturated power models with different semiconductor
technologies........................................................................................................55
4.2.2.2 Simplified memoryless amplifier models at 300 GHz .......................................57
4.2.2.3 Behavioural PA modelling with memory ...........................................................60
4.3 D-MIMO radio architecture..........................................................................................63
4.3.1 Introduction .............................................................................................................63
4.3.2 Architectures ...........................................................................................................65
List of Figures
Figure 2-1: Generic radio architecture .................................................................................................. 22
Figure 2-2: Deployment scenarios of different radio hardware options. .............................................. 23
Figure 2-3 Radio design flow chart....................................................................................................... 24
Figure 3-1: Multipaths inside a material slab........................................................................................ 27
Figure 3-2: pp_HrITUf and pp_HtITUf for glass .................................................................................. 27
Figure 3-3: Measurement system. ......................................................................................................... 28
Figure 3-4: pp_hrMUTt for plexiglass. ................................................................................................. 29
Figure 3-5: Permittivity estimation method. ......................................................................................... 29
Figure 3-6: Estimation method comparison. ......................................................................................... 29
Figure 3-7: pp_htMUTt for mortar. ...................................................................................................... 29
Figure 3-8: Material permittivity. ......................................................................................................... 30
Figure 3-9: pp_HrMUT(f) and pp_HtMUTf compared with ITU and ITUF models. Plexiglass (a1,b1),
Mortar (a2,b2), Melamine chipboard (a3,b3). ...................................................................................... 31
Figure 3-10: (a) Band- and aperture-limited channel response from a measurement and (b) its band- and
aperture-unlimited model as propagation paths. Data are from a shopping mall measurement at 140
GHz [NJK+18]. ..................................................................................................................................... 34
Figure 3-11: Power delay profiles of three example links of the stored channel model. ...................... 36
Figure 3-12: Time variant attenuation pattern of frontal and lateral human blockage events. ............. 36
Figure 3-13: Angles and path gains of three example links of the stored channel model..................... 36
Figure 4-1: System model for 6G radio analysis from end-to-end link perspective. Beamforming,
MIMO, and distributed radio network need to consider the whole signal path and limited isolation (i.e.
coupling) between orthogonal channels. [HEX21-D21]. ...................................................................... 37
Figure 4-2: RF transceiver reference architecture for hybrid beamforming based on sub-array per beam
principle. Key parameters for model abstraction are named. ................................................................ 38
Figure 4-3 Required relative BWs for fixed data rate with different modulations and coding rates. ... 41
Figure 4-4: ADC dynamic range and sampling rate requirements for various combinations resulting to
100Gbps data rate. ................................................................................................................................ 43
Figure 4-5: ADC dynamic range (SNDR) and power consumption (P) is compared to ADC dynamic
range requirements for the OFDM and CW waveforms (solid oval). Some commercial state-of-the-art
ADCs with high resolution and sampling rate (dashed ovals) have been added to complement the data
from [Mur22]. ....................................................................................................................................... 43
Figure 4-6: EVM for 5GNR and SC (r=0.35) 64QAM modulations from the 1 dB compression point of
the receiver (or transmitter) where non-linearity is dominated by the 3rd order term. ......................... 44
Figure 4-7: Probability of the peaks above average signal power. ....................................................... 45
Figure 4-8: Number of antennas, EIRP and relative bandwidth as a function of frequency for 100 Gbps
target data rate with different modulation orders and link distance. Legends and axis to which each of
the curves points are shown in figure (a). ............................................................................................. 48
Figure 4-9: Number of antennas, EIRP and relative bandwidth as a function of frequency for 100 Gbps
considering 10 dB and 20 dB to describe pathloss of NLOS radio channel. ........................................ 48
Figure 4-10: LO chain with frequency multipliers. .............................................................................. 49
Figure 5-8: Processing time scaling for SCS of 1.92 MHz. .................................................................. 76
Figure 5-9: Comparison of SC-FDE and DFT-s-OFDM with phase noise and non-linear PA distortion.
.............................................................................................................................................................. 77
Figure 5-10: Dependency on the deployed number of quantization bits of SC-FDE (upper row) and
DFT-S-OFDM (lower row) for QPSK, 16QAM and 64QAM. ............................................................ 79
Figure 5-11: 4-bits constellation pattern for standard NR 16QAM (left) and PA friendly 16APSK
(right). ................................................................................................................................................... 80
Figure 5-12: 4-bits modulation: CCDF of the PAPR (left) and BLER as a function of required Eb/No
(right) of SC-FDE waveform using NR 16QAM and 16APSK constellation pattern. ......................... 81
Figure 5-13: 6-bits modulation: CCDF of the PAPR (left) and BLER as a function of required Eb/No
(right) of SC-FDE waveform using NR 16QAM and 16APSK constellation pattern. Inset: applied 4
rings constellation pattern for 64APSK. ............................................................................................... 81
Figure 5-14: Complex baseband of a DFTS-OFDM transmitter. ......................................................... 82
Figure 5-15: Complex baseband of a FDF-DFTS-OFDM transmitter.................................................. 83
Figure 5-16: Principle of power boosting. ............................................................................................ 83
Figure 5-17: Block PAPR for the set of chosen waveforms. ................................................................ 84
Figure 5-18: link budget comparison of two systems using different waveforms. ............................... 84
Figure 5-19: Minimum OBO satisfying EVM and ACLR requirements. ............................................. 85
Figure 5-20: BLER comparison of DFTS-OFDM and FDF-DFTS-OFDM under the influence of PA
nonlinearity. .......................................................................................................................................... 86
Figure 5-21: MCL difference vs throughput. ........................................................................................ 86
Figure 5-22: Block Error Rate (BLER) for 16-QAM, R = ¾, subcarrier spacing of 960, 1920 and 3840
kHz, with Ng=2,4,8,16,32 and 64 pilot groups. .................................................................................... 88
Figure 5-23: Model of the time frame structure, taken from [GSD+22]. ............................................. 90
Figure 5-24: MSE of the estimators and BCRB over the sample index k within the observation frame,
taken from [ZGD+22]. .......................................................................................................................... 92
Figure 5-25: Illustration of typical out-phasing concept. ...................................................................... 93
Figure 5-26: Illustration of over-the-air combined outphasing concept. .............................................. 93
Figure 5-27: Simulated data rates at different bandwidth and link distance in four scenarios (a)-(d)
derived in Table 14. .............................................................................................................................. 96
Figure 5-28: Beam management finite state machine. .......................................................................... 97
Figure 5-29: L1 mobility in D-MIMO with spatial repetition transmission. ........................................ 99
Figure 5-30: SE performance of CZF and IADZF with CJT for different number of serving RUs in a
scenario with 16 RUs. ......................................................................................................................... 100
Figure 5-31: SE performance of CZF with CJT for different serving RU subset periodicities, e.g., 100
ms, 500 ms and 1 s. ............................................................................................................................. 100
Figure 5-32: SE performance for CJT and NCJT schemes for different number of serving RUs in the
scenario of 16 UEs when serving RU subset update periodicity is 100 ms. ....................................... 100
Figure 5-33: SE behaviour evaluated in 10s-time frame for CJT and NCJT schemes with serving subset
size of 4 RUs and different subset update periodicities. ..................................................................... 100
Figure 5-34: Modelling the topology of a dynamic D-MIMO network as a bipartite graph with partial
connectivity. ........................................................................................................................................ 101
Figure 5-35: CDF of the downlink achievable rate per user under the dynamic scenario. ................. 103
Figure 5-36: Rate versus pilot dimension under the dynamic scenario. ............................................. 103
Figure 5-37: Experimental setup block diagram. ................................................................................ 104
Figure 5-38: NCR-based DL system. .................................................................................................. 105
Figure 5-39: An illustration of how NCR helps the network in the presence of blockage. ................ 106
Figure 5-40: The required number of RIS elements to achieve the same performance as in the NCR-
assisted network .................................................................................................................................. 106
Figure 5-41: Generalized mesh front/backhaul problem formulation................................................. 106
Figure 5-42: Mesh resource optimization options. ............................................................................. 107
Figure 5-43: Outage for Single and Dual connectivity in the UE-RU link. Outage is caused by
shadowing/blocking which is a function of distance. Link-loss in mesh backhaul caused by interference.
The number of APs with core connections was 9. Network simulation parameters are summarized in
Table 5-17. .......................................................................................................................................... 108
Figure 5-44: Distribution of number of hops for a different of core-DU/RUs and a different number of
RUss. The algorithm used for calculating routes through the mesh was an iteratively applied Dijkstra
optimizing for link SINR. ................................................................................................................... 109
Figure 5-45: IAB network with deployment constraints..................................................................... 109
Figure 5-46: Service coverage as a function of the minimum distance constraint between the nodes for
different UE data rate thresholds. ....................................................................................................... 110
Figure 5-47: Symmetric IAB setup with donor at the centre .............................................................. 111
Figure 5-48: Service coverage as a function of the minimum distance constraint between the nodes in
different UE densities for blocking aware and blocking unaware optimization scenarios. ................ 111
Figure 5-49: Service coverage as a function of the radius of the constrained areas for different
deployment scenarios. ......................................................................................................................... 111
Figure 6-1: Deployment scenario with 2 BSs (left) and 4 (right) BSs illuminating the area from the
ceiling.................................................................................................................................................. 113
Figure 6-2: Subpanel configurations with 1, 2 and 4 subpanels. ........................................................ 114
Figure 6-3: Area throughput for different number of simultaneously served UEs. ............................ 114
Figure 6-4: Area throughput for different number of simultaneously served UEs. ............................ 115
Figure 6-5: UE Throughput CDF for 2 and 4 BS over the whole deployment area, 1 subpanel and
different number of simultaneously served UEs. ................................................................................ 115
Figure 6-6: Exemplary power consumption for 2- and 4-BS scenarios and different ADC resolution.
............................................................................................................................................................ 116
Figure 6-7: Path loss models. .............................................................................................................. 117
Figure 6-8 Measured path powers, the beam power and found independent beam azimuth directions of
an example link. .................................................................................................................................. 118
Figure 6-9 Number of independent beams in 132 indoor links. ......................................................... 118
Figure 6-10 CDF of the number of independent beams in 132 indoor and 157 outdoor links using 10,
and 20 dB dynamic range. .................................................................................................................. 118
List of Tables
Table 2-1: sub-THz communication scenarios requirements................................................................ 21
Table 3-1: Estimated and ITU permittivity and conductivity. .............................................................. 32
Table 3-2: 140 GHz spatio-temporal channel sounding. ...................................................................... 34
Table 4-1: Common RF and waveform parameters (assumptions) used in the example analysis. ....... 41
Table 4-2: Required bandwidth for fixed data rate using different modulations with guard bands...... 41
Table 4-3: ADC sampling rate (Gsps) for 100 Gbps data rate for different modulations and coding rates.
.............................................................................................................................................................. 42
Table 4-4: ADC dynamic range and sampling rate requirements for various combinations resulting to
100 Gbps data rate with 5/6 coding rate. .............................................................................................. 42
Table 4-5: Minimum SNR/EVM and PAPR (i.e. required BO in this case) for different modulations in
the analysis. ........................................................................................................................................... 45
Table 4-6: HW models used in the analysis. ......................................................................................... 47
Table 4-7: Model parametrization for phase noise model 1 given in (4-4)........................................... 51
Table 4-8: Model parametrization for phase noise model 2 given in (4-4)........................................... 51
Table 4-9: Simulation parameters for link phase noise analysis. .......................................................... 53
Table 4-10: Common memoryless nonlinear models of RF amplifiers. ............................................... 58
Table 4-11: Parameter set for common memoryless amplifiers models fitted for circuit level simulations
(sim) and measurements (meas) of 290 GHz SiGe amplifier. .............................................................. 60
Table 5-1: Qualitative assessment of selected waveforms. ................................................................... 71
Table 5-2: Extrapolation of NR numerology for FFT size 4096. ......................................................... 72
Table 5-3: Extrapolation of NR SCS for FFT size of 1024 and 2048................................................... 72
Table 5-4: Terms for user plane latency evaluation. ............................................................................. 73
Table 5-5: Assumed values of parameters for latency calculation, processing time tied to slot duration.
.............................................................................................................................................................. 74
Table 5-6: Assumed values of parameters for latency calculation, processing time independent of slot
duration. ................................................................................................................................................ 75
Table 5-7: System parameters for DFTS-OFDM and SC-FDE comparison study. .............................. 77
Table 5-8: SNR in dB at which a BLER of 10-1 is reached dependent on the OPBO........................... 78
Table 5-9: Parameter setting. ................................................................................................................ 78
Table 5-10: EVM requirements. ........................................................................................................... 85
Table 5-11: Indices of the time-domain pilots for PN estimation ......................................................... 87
Table 5-12: Loss in terms of Eb/N0 for 16QAM pilot schemes at BLER = 10-2. NA = not achieved. 88
Table 5-13: Loss in terms of Eb/N0 for 64QAM pilot schemes at BLER = 10-2. NA = not achieved. 88
Table 5-14 Simulation parameters. ....................................................................................................... 92
Table 5-15: Simulation parameters for data-rate analysis with different link distances and signal
bandwidth.............................................................................................................................................. 95
Table 5-16: Simulation Model Specifications .................................................................................... 100
Table 5-17: Basic network parameters................................................................................................ 108
Dissemination level: public Page 11 / 127
Hexa-X Deliverable D2.3
List of Acronyms
Term Description
3D 3 Dimensional
3GPP 3rd Generation Partnership Project
5G 5th Generation mobile communication
5G-PPP 5G Public Private Partnership
6G 6th Generation mobile communication
ABG Alpha-Beta-Gamma
ACLR Adjacent Channel Leakage Power
ADC Analogue-to-Digital Converter
AGC Automatic Gain Control
AM-AM Amplitude-to-Amplitude
AM-PM Amplitude-to-Phase
AoA Angle-of-Arrival
AoD Angle-of-Departure
AP Access Point
APSK Amplitude Phase Shift Keying
AR Augmented Reality
ARoF Analog Radio over Fibre
AWG Arbitrary Waveform Generator
AWGN Additive White Gaussian Noise
B5G Beyond 5G
BB Baseband
BCRB Bayesian Cramér-Rao bound.
BER Bit Error Rate
BFW Beamforming Weights
BiCMOS Bipolar CMOS
BLER Block Error Rate
BM Beam Management
BS Base Station
BW Bandwidth
CA Carrier Aggregation
CC Component Carriers
CCDF Complementary Cumulative Density Function
CI Close-In
CJT Coherent Joint Transmission
CMOS Complementary Metal Oxide Semiconductor
CPE Common Phase Error
HW Hardware
HWI Hardware Impairment
i.i.d. independent and identically distributed
IAB Integrated Access and Backhaul
IBKM Information-Based k-Means
IC Integrated Circuit
ICI Inter-Carrier Interference
ICP Input Compression Point
IF Intermediate Frequency
IFFT Inverse Fast Fourier Transform
InP Indium Phosphide
IQ In-phase and Quadrature
ISI Inter Symbol Interference
KPI Key Performance Indicator
LCA Life Cycle Assessment
LDPC Low-Density Parity-Check
LINC Linear Amplification with Nonlinear Components
LLS Lower Layer Split
LNA Low Noise Amplifier
LO Local Oscillator
LOS Line-of-Sight
MAC Medium Access Control
MCL Maximum Coupling Loss
MDF Medium Density Fibreboard
MDS Maximum Distance Separable
MIMO Multiple Input Multiple Output
MMSE Minimum Mean Square Error
mmWave Millimetre-wave
MR Mixed Reality
MSE Mean-Squared Error
MT Mobile Termination
MTC Machine Type Communication
MZM Mach-Zehnder Modulator
NCJT Non-coherent Joint Transmission
NCR Network-Controlled Repeater
NF Noise Figure
NLOS Non-Line-of-Sight
NR New Radio
NRZI Non-Return-to-Zero Inverse
1 Introduction
Hexa-X is one of the 5G-PPP projects under the EU Horizon 2020 framework. It is a flagship project
that develops a Beyond 5G (B5G)/6G vision and an intelligent fabric of technology enablers connecting
human, physical and digital worlds.
This document is the third deliverable (D2.3) of Work Package 2 (WP2) - “Novel radio access
technologies towards 6G”. The work in WP2 focuses on radio design aspects for future mobile
communication systems, including radio frequency (RF) transceiver implementation and hardware
(HW) modelling in sub-THz bands (100 – 300 GHz), distributed multiple-input multiple-output (D-
MIMO) architectures, and signal processing and transmission schemes. The research focuses on the
following key aspects:
• Radio enablers and technology roadmap
• Radio and antenna implementation aspects, hardware component models and architecture
• Hardware-aware waveform and modulation design
• Hardware-aware beamforming design
• Distributed large MIMO systems for beyond 5G and 6G
• Channel measurement and modelling above 100 GHz
The first deliverable of WP2 (D2.1) [HEX21-D21] presents use cases related to Tbps communications,
and identifies gaps in the current and future technologies, which steer the work in WP2. The second
deliverable (D2.2) [HEX21-D21] provides initial analysis of radio models and performance analysis at
the radio system level, modelling of non-linearities, channel measurement and modelling, and
waveform analysis with the modelled radio impairments.
In Section 5, signal processing topics related to the design of high-throughput 6G systems are presented.
First, an overview is given of selected waveform design topics and related digital transceiver design
guidelines. Then, guidelines are given on how the beam management techniques need to be adapted to
the large number of beams needed in sub-THz systems. Thereafter, a look into a selection of topics in
the area of D-MIMO and integrated access and backhaul is provided.
Section 6 is dedicated for system performance analysis, where the impact of deployment scenarios on
the power consumption is evaluated. In addition, insights are provided on the influence of the
propagation channel on the link quality based on channel measurements in several environments.
The document is concluded in Section 7.
malls and subways, or outdoor streets and parks. While considering multiple users, the deployment
option should be high density cells. The involved devices include augmented reality (AR) glasses and
tactile gloves that are locally connected to a terminal that is connected to an access point (AP). As the
devices are held, or worn, by a human, the mobility is at walking speed and can be stationary. Depending
on the application, the data rate ranges from 1 Gbps to 100 Gbps with E2E latency less than 20 ms. The
positioning accuracy should be on the order of 1 cm and the latency of less than 100 ms.
Fixed wireless access and fronthaul: In this case, sub-THz can be used to provide directional wireless
link between stationary nodes with high data rate (100 Gbps) and low latency (0.1-100 ms) as alternative
to optical fibres for fronthaul links for interconnecting small cell APs, connecting remote digital twins,
and providing fixed wireless access for residential users. Such link is commonly deployed in outdoor
environment between stationary APs, but can also be relevant to indoor. For supporting the sensor
infrastructure web, location accuracies around 0.1 - 1 meter with latencies of 10-100 ms are required.
Mid-range wireless access: this scenario can be encountered in use cases that require fixed assess with
no mobility. The targeted data rate is 100 Gbps at a maximum range of 200 m, which can be realized
Dissemination level: public Page 21 / 127
Hexa-X Deliverable D2.3
by means of high gain antenna arrays. The communicating nodes such as APs have the same
capabilities, and therefore, the same radio design can serve for both end nodes.
Short-range wireless access: this scenario requires lower peak data rate, namely 10 Gbps, but mobility
should be considered, which requires efforts to cope with the channel variation via beam tracking.
When the mobility is controllable in terms of location and trajectory and/or other contextual information
(e.g., device position and orientation) is available, such information can be used to predict the beam
switching. The maximum range to be covered is 10-100 m in indoor or outdoor environment, with the
possibility to exploit spatial multiplexing. Two radio HW needs to be designed differently for the
infrastructure AP and the mobile terminal considering the constraints on size, cost, and energy
consumption.
Very short-range wireless access: this scenario is defined by 100 Gbps peak data rate and a very short
distance - up to 10 m – for providing connectivity to mixed reality and holographic communications
with E2E latency less than 20 ms depending on the application. The corresponding channel is indoor
and outdoor, with mobility at walking speed. For different constraints on the AP and user equipment,
as in the short-range scenario, two different radio HW are needed.
the radio options in RAN is shown in Figure 2-2, which highlights the three scenarios and related use
cases.
error rate (SER) computed by the Q-function. In this step, the impact of the number of possible
spatial streams is provided.
2- Link model: this consists of the transmitter and receiver SNR models, considering the HW
impairments that depend on the signal, such as PA and LNA non-linearities, phase noise,
quantization noise, and linear distortions such as IQ imbalance and filtering, in addition to the
additive noise, which depends on the bandwidth and noise figure. Moreover, the impact of DSP,
waveform, and other communication overhead should be considered. The receive signal power
depends on the transmit power and the link budget. The link budget is a function of range,
carrier frequency, and antenna array and antenna element gains, and it is affected by the path
losses of the propagation environment. Regulatory requirements may limit the maximum total
transmit power and the spectral mask of the transmitted signal that needs to be considered. In
addition, it is necessary to include margins to account for unforeseen or unlikely effects during
the system design. This model involves many parameters, and some of them are bounded by
the HW limitations, such as the PA transmit power, the antenna array size and number of
elements constraints. These constraints need to be considered in limiting the number of possible
solutions. Another factor to be considered is to minimize the power consumption and the
effective isotropic radiated power (EIRP) radiation, in addition to reducing the cost.
3- Evaluation under radio realistic channel: the parameters obtained need to be evaluated using
channel measurements to determine the possible number of independent spatial directions. In
theory, this evaluation requires accessing full channel measurements between all pairs of
antennas to study the channel rank seen by the antenna array. Each singular value represents
the path gain of a spatial direction/beam. However, in a practical sense, obtaining the spatial
information can be achieve by assessing several directions, each served by a subarray. High
number of independent directions/beams allows the adaptation of transmit power or reducing
the bandwidth. On the contrary, in the case of a single non-line-of-sight (NLoS) beam, the link
budget may deteriorate, and thus, an update of the design is required. Therefore, an extensive
simulation covering different locations determines whether the design requires changes or it
can be used in an adaptive way.
4- Evaluation of DSP algorithms: in order to validate the overall air interface design, the selected
radio architecture should be used to evaluate the selected waveforms and verify the ability of
achieving reliable transmission within the latency constraints of the use cases. In addition,
beamforming procedures for initial access and management should be evaluated on top of the
selected architectures for the corresponding scenarios.
3 Channel model
This section describes channel models and their input parameters that were populated in the project. In
order to support various types of channel models ranging from fully physical model, e.g., ray-optics
simulations, to fully stochastic channel model based on the geometry, two basis channel models are
provided: 1) material parameters for 2-260 GHz and 2) stored channel model at 140 GHz that uses
processed measured channel responses as a seed to reproduce channels based on a Monte Carlo
simulation.
Air Air
Material
t
Incident wave Path 1T
Path 1R
t Transmitted
paths
Path 2R
Reflected
paths
…..
Path 2T
…..
…..
The ITU recommendation [ITU P.2040-2] proposes a constant permittivity value regardless the
frequency and expresses the frequency-dependent conductivity 𝝈 as 𝝈 = 𝒄𝒇𝒅 . 𝜺𝒓 , c and d are given for
some materials such as glass. 𝐻𝑟𝐼𝑇𝑈 (𝑓) and 𝐻𝑡𝐼𝑇𝑈 (𝑓) are 𝐻𝑟𝑀𝑜𝑑 (𝑓) and 𝐻𝑡𝑀𝑜𝑑 (𝑓) applying the ITU
model. 𝐻𝑟𝐼𝑇𝑈 (𝑓) and 𝐻𝑡𝐼𝑇𝑈 (𝑓) are illustrated in Figure 3-2 for the glass (𝑑 = 8 mm, 𝜀𝑟 = 6.25, 𝜎 =
0.0043𝑓 1.19 ).
When analysing propagation channels in indoor or outdoor environments, the large-scale effects are
defined as referring to the environment or Tx/Rx distance effects, while the small-scale effects refer to
the multipath complex combination. It can be observed in Figure 3-2 similar effects with a large-scale
frequency variation that depends on the permittivity and conductivity. The higher the permittivity, the
lower the reflection losses. The higher the conductivity, the higher the penetration losses. Small-scale
effects are created by the interferences between the multipath created inside the slab. Due to the
penetration loss increase with the frequency, small-scale effects tend to become less relevant at high
frequencies.
As a first approach, the discussion about the ITU model validity above 100 GHz may be divided in two
sub-discussions. The first sub-discussion is about the validity of the ITU theoretical framework. For
instance, does the permittivity remain constant and is the conductivity modelled by 𝑐𝑓 𝑑 ? Are the
conductivity and permittivity still relevant parameters to compute reflection and transmission losses?
The second sub-discussion is about the numerical values of 𝜀𝑟 , c and d. Shall they be updated to
correctly represent reflection and transmission losses? An ITU-like model was defined to help the
analysis related to the first sub-discussion. This model is called ITUF and keeps the same theoretical
approach of the ITU model but parameters 𝜀𝑟 , c and d are fitted on measurement. The notations
𝐻𝑟𝐼𝑇𝑈𝐹 (𝑓) and 𝐻𝑡𝐼𝑇𝑈𝐹 (𝑓) correspond to the transfer functions 𝐻𝑟𝑀𝑜𝑑 (𝑓) and 𝐻𝑡𝑀𝑜𝑑 (𝑓) when applying
the ITUF model. A pp prefix added to H functions indicates the power profile of H functions with
𝑝𝑝_𝐻(𝑓) = 10 ∗ log10 (|𝐻(𝑓)|2 ).
for reflection and transmission measurements, respectively. Moreover, 𝐻𝑟 (𝑓) and 𝐻𝑡𝑀𝑈𝑇 (𝑓)
𝑀𝑈𝑇
represent the intrinsic material transmission and reflection frequency response. They are computed
according to (3-1),
𝑀𝑒𝑎𝑠
𝑇𝐺[𝑆11 𝐹𝑆
(𝑓) − 𝑆11 (𝑓)] 𝑀𝑒𝑎𝑠
𝑇𝐺[𝑆21 (𝑓)] (3-1)
𝐻𝑟𝑀𝑈𝑇 (𝑓) = 𝑀𝑒𝑡 𝐹𝑆 𝐻𝑡𝑀𝑈𝑇 (𝑓) = 𝐹𝑆 ,
𝑇𝐺[𝑆11 (𝑓) − 𝑆11 (𝑓)] 𝑇𝐺[𝑆21 (𝑓)]
where 𝑇𝐺[. ] is the time gating operator. It is a temporal filter applied on the inverse fast Fourier
transform (IFFT) of S parameters that rejects all undesirable components such as wall reflections or
material edge diffractions. Here, ℎ𝑟𝑀𝑈𝑇 (𝑡) and ℎ𝑡𝑀𝑈𝑇 (𝑡) denote the IFFT of 𝐻𝑟𝑀𝑈𝑇 (𝑓) and 𝐻𝑡𝑀𝑈𝑇 (𝑓),
respectively. In the following sections, the data analysis and model comparison will be mainly based
on power profile. A pp_ prefix added to a complex vector indicates the vector power profile. For
instance, 𝑝𝑝_ℎ𝑟𝑀𝑈𝑇 (𝑡) is the power profile of ℎ𝑟𝑀𝑈𝑇 (𝑡) and is equal to 10 ∗ log10 (|ℎ𝑟𝑀𝑈𝑇 (𝑡)|2 ).
Path 1R Path 1T
Consequently, separate estimation methods are applied for the permittivity and conductivity. First, the
permittivity is estimated by two different methods based either on the delay of transmitted paths or on
the amplitude of reflected paths. Only path 1R amplitude and path 1T delay are considered, as
illustrated in Figure 3-5, because second order paths are often below the noise level at the highest
frequencies. The initial measurement sub-bands were slightly rearranged in order to define nine new
sub-bands of about 30 GHz bandwidth (2-30 GHz, 30-50 GHz, 50-75 GHz, 75-110 GHz, 110-140 GHz,
140-170 GHz, 170-200 GHz, 200-230 GHz, 230-260 GHz). An analysed bandwidth of roughly 30 GHz
is a good trade-off to keep a good resolution for the path detection, to consider a constant value for the
permittivity and to assess the large-scale frequency dependency of permittivity.
Method 1 is based on the delay. Let 𝑡𝑃1𝑇 be the delay of path 1T. The wave speed inside the material is
𝑐
equal to 𝜀𝑣 , 𝑐𝑣 being the wave speed in vacuum. Then, the relative permittivity 𝜀𝑟 is calculated as
√ 𝑟
expressed in (3-2).
𝑡1𝑇 ∗ 𝑐𝑣 2 (3-2)
𝜀𝑟 = ( + 1) .
𝑑
√ 𝜀𝑟 − 1
Method 2 is based on the amplitude. Let’s define ℎ1𝑅
𝑟 the amplitude of path 1R. ℎ𝑟 =
1𝑅
, then, 𝜀𝑟
√𝜀𝑟 +1
can be calculated as in (3-3).
1 + ℎ1𝑅
2 (3-3)
𝑟
𝜀𝑟 = ( ) .
1 − ℎ1𝑅
𝑟
Figure 3-6 shows the permittivity estimated for the plexiglass, mortar, and chipboard. The permittivity
estimation is sometimes unreliable for rough surface materials or inhomogeneous materials, or multi-
layer materials because 𝑝𝑝_ℎ𝑟𝑀𝑈𝑇 (𝑡) or 𝑝𝑝_ℎ𝑡𝑀𝑈𝑇 (𝑡) do not show sufficient distinction of path 1R or
path 1T. For instance, as illustrated in Figure 3-6, the mortar permittivity was not estimated by method
1 above 140 GHz. Permittivity values estimated by method 1 and method 2 are very similar for most of
the materials but diverge for some materials when the frequency increase. It can be observed in Figure
3-6 that the permittivity estimated by method 2 for the mortar decreases artificially due to the scattering
from rough surfaces when the wavelength is similar to the surface small irregularities. This point is
discussed more in details in the next paragraph. It can also be seen that the chipboard permittivity
estimated by method 2 strongly increases with frequency, whereas the permittivity estimated by method
1 remains constant. The permittivity increases by method 2 due to the melamine at the chipboard surface
or due to the chipboard composite structure. The small wood particles may become more efficient
reflectors when the wavelength decreases.
Even if the permittivity is mainly used to compute reflection losses in propagation channel simulation,
the final permittivity values are estimated by method 1, i.e., based on the delay of transmission
measurement, because this method seems to be more representative of the material electrical properties
and less dependent of the material surface state.
Figure 3-8 presents the permittivity for various materials. It confirms that the permittivity is independent
of the frequency even at frequency above 100 GHz. The permittivity is averaged over the different
frequency bands, and the mean value is reported in Table 3-1
more due to the complex recombination between path 1R and path 2R but is due to the complex
recombination of scattered components. The scattering increases the reflection loss average value and
makes the ITUF model inaccurate or not valid as the difference model/measurement may be higher than
10 dBs. It can be seen that, for some materials (MDF, PFB, glasswool, etc), 𝑝𝑝_𝐻𝑟𝑀𝑈𝑇 (𝑓) increases
with frequency as illustrated by the melamine chipboard example in Figure 3-9. The 𝑝𝑝_𝐻𝑟𝑀𝑈𝑇 (𝑓)
increase compared to 𝑝𝑝_𝐻𝑟𝐼𝑇𝑈𝐹 (𝑓) is generally limited to a few dBs.
(a1) (b1)
(a2) (b2)
(a3) (b3)
Figure 3-9: 𝒑𝒑_𝑯𝑴𝑼𝑻
𝒓 (𝒇) and 𝒑𝒑_𝑯𝑴𝑼𝑻
𝒕 (𝒇) compared with ITU and ITUF models. Plexiglass (a1,b1),
Mortar (a2,b2), Melamine chipboard (a3,b3).
The analysis of transmission results leads to similar conclusions. The frequency-dependent conductivity
can be well modelled by the ITUF model “on average”. The absolute difference between 𝑝𝑝_𝐻𝑡𝐼𝑇𝑈𝐹 (𝑓)
and 𝑝𝑝_𝐻𝑟𝑀𝑈𝑇 (𝑓) is generally about a few dB for frequencies less than 60 GHz but may exceed 10 dB
for frequencies above 100 GHz. For instance, the ITUF model agrees with mortar or chipboard results
up to 60 GHz, but some variations appear above 60 GHz. Above 100 GHz, the transmission loss
differences between the ITUF model and measurements can exceed 10 dB, and the three points results
are uncorrelated. Figure 3-7 shows 𝑝𝑝_ℎ𝑡𝑀𝑢𝑡 (𝑡) for the mortar at different sub-bands. Path 1T is highly
attenuated, and some delayed components appear at frequencies above 100 GHz. It suggests that the
volumic diffuse scattered component inside the material due to the material inhomogeneity becomes
dominant over the specular transmission component.
frequencies, the ITU model could be improved by adding to 𝐻𝑡𝐼𝑇𝑈 (𝑓) a statistical variable but as
indicated hereabove for the reflection loss, this approach requires additional measurements
Finally, it can be concluded that the ITU model needs to be improved for frequencies above 100 GHz
to represent the measured reality of reflection and transmission well, while it could be used as it is in
many simulations related to 6G sub-THz scenarios. Reflection loss errors due to rough surfaces may
not be a concern in office environment, shopping mall, airport, etc as most of the material are quite
smooth. Transmission loss errors due to the material inhomogeneity may not be a concern as they
concern high transmission losses at high frequencies. Simulating a transmission error of 20 dB instead
of 30 dB may not significantly impact system simulation if a blockage is defined by a transmission loss
higher than 20 dB. Future work will be focused on multi-frequency measurement in real environments
to better evaluate the propagation channel differences between cmW and sub-THz frequencies.
(a)
(b)
Figure 3-10: (a) Band- and aperture-limited channel response from a measurement and (b) its band- and
aperture-unlimited model as propagation paths. Data are from a shopping mall measurement at 140 GHz
[NJK+18].
were obtained for each radio link; 𝑞 is the link index, 𝐿𝑞 is the number of paths, δ(⋅) is the delta
function, 𝑃𝑞,𝑙 , φrx
𝑞,𝑙 , and τ𝑞,𝑙 are the power (squared magnitude of path gain), the azimuth direction of
arrival, and the propagation delay of the 𝑙th path, respectively. The missing angular knowledge of
multipath components, e.g., azimuth angles on the Tx side and elevation angles at the Rx, was
supplemented by fusing the available multipath parameter estimates, i.e.,(3-4), with the available
detailed geometric database of the measurement environment [DKH21]. The resulting channel data,
which are published in [DHK23], consist of double-directional multipath components denoted as
𝐿𝑞
where Ωrx
𝑞,𝑙 and Ωtx
𝑞,𝑙are the direction of arrival and departure, including azimuth and elevation angles.
Current data are measured with vertically polarized antennas; hence the definitions are restricted here
to a single polarized case only, neglecting the polarization characteristics.
Table 3-2: 140 GHz spatio-temporal channel sounding.
Entrance Hall Suburban Residential City Centre
Number of LOS Links 5 5 5 5
Number of NLOS Links 12 32 13 19
Link Distance Range (m) 3-66 2-172 20-175 10-178
Random snapshots of frequency response matrices can be generated by introducing random initial phase
𝜑𝑞,𝑙 for each multi-path, where phase terms are drawn from the uniform distribution in [0,2𝜋] .
Moreover, the temporal dimension and time variability can be included by introducing small Doppler
frequencies 𝜈𝑙,𝑞 for each path. This models a small-scale virtual motion, where only phases of path
component change over time, but other propagation parameters remain constant. The resulting
snapshot/time variant frequency response matrix is
𝐿𝑞
𝑇
𝐇𝑞 (𝑡, 𝑓) = ∑ 𝐠 rx (Ωrx
𝑞,𝑙 ) √𝑃𝑞,𝑙 𝑒
𝑗(𝜑𝑞,𝑙 +2𝜋𝜈𝑞,𝑙 𝑡) −𝑗2𝜋𝑓𝜏𝑞,𝑙
𝑒 𝐠 tx (Ωtx 𝑀×𝑁 (3-7)
𝑞,𝑙 ) ∈ ℂ .
𝑙=1
Examples of band-/aperture-unlimited double directional power angular delay profiles (PADPs) are
shown in Figure 3-11 and Figure 3-13. Three colours indicate three propagation paths of three Tx/Rx
locations measured at 140 GHz in an indoor environment. Figure 3-11 illustrates power delay profiles
(PDP) of measured propagation channel on three different Tx/Rx locations. One can observe that within
examples, the maximum excess delay varies between 60 and 200 ns. Moreover, e.g., the link number 3
(blue) has 12 multipath within 20 dB dynamic range. Figure 3-13 depicts power-angular information of
the same paths. The top left figure indicates azimuth arrival/departure angles, top right elevation
arrival/departure angles, bottom left arrival azimuth/elevation angles, and bottom right departure
arrival/departure angles. In the illustrated examples, the elevation and azimuth angles are within
approximately 40 deg and 220 deg range, respectively.
Another source of time variability is the blockage of a propagation path by a moving obstacle, e.g., a
human body. Such events are probable, especially in short link distance indoor environments, when
antennas are not elevated high and persons are moving in the environment. Short wavelengths at sub-
THz radio frequencies make Fresnel ellipsoids of each multipath component small and a human body
can severely block the signal up to 40 dB attenuation. The human blockage effect was measured at 140
GHz frequency at 3.5 m link distance with different persons and orientations of persons in [ZBK+23].
Blockage events are modelled as supplemental time variability in (3-7) by introducing attenuation term
𝛼(𝑡) such that the resulting transfer function becomes
𝐿𝑞
𝑇
𝐇𝑞 (𝑡, 𝑓) = ∑ 𝐠 rx (Ωrx
𝑞,𝑙 ) √𝛼𝑞,𝑙 (𝑡)𝑃𝑞,𝑙 𝑒
𝑗(𝜑𝑞,𝑙 +2𝜋𝜈𝑞,𝑙 𝑡) −𝑗2𝜋𝑓𝜏𝑞,𝑙
𝑒 𝐠 tx (Ωtx 𝑀×𝑁 (3-8)
𝑞,𝑙 ) ∈ ℂ .
𝑙=1
The attenuation term is taken directly from the measurement reported in [ZBK+23]. Figure 3-12 shows
two example attenuation patterns, one with lateral and one with frontal crossing of a human body across
a multipath component. In the channel model one can choose which path or paths are blocked, the initial
time instance of each blockage event, and the type of blockage, i.e., lateral or frontal crossing of a
human body. These choices can be made also random. If path 𝑙 experiences blockage, the term 𝛼𝑞,𝑙 (𝑡)
gets such values as in Figure 3-12 in linear units, starting with the selected time instant. Term 𝛼𝑞,𝑙 (𝑡) =
1 ∀𝑡 are used for those paths without blockage event. A few remarks can be made from Figure 3-12
about the human shadowing at short link distances. The maximum attenuation reaches approximately
40 dB. The frontal blocking has shorter duration and slightly less severe attenuation, as expected. There
are diffraction patterns just before after the actual blockage. Potentially these power variations could be
used to predict an upcoming blockage event.
Figure 3-11: Power delay profiles of three Figure 3-12: Time variant attenuation pattern of
example links of the stored channel model. frontal and lateral human blockage events.
Figure 3-13: Angles and path gains of three example links of the stored channel model.
3.3 Summary
This section presented a Hexa-X channel model, consisting of 1) a list of complex permittivity estimates
for materials that make up our living environments for a frequency range from 2 to 260 GHz and 2)
wideband MIMO channel model at 140 GHz based on measured spatio-temporal channel sounding data
in indoor and outdoor scenarios. The latter is implemented in MATLAB and is published in Hexa-X
Zenodo community [DHK23].
Figure 4-1: System model for 6G radio analysis from end-to-end link perspective. Beamforming, MIMO,
and distributed radio network need to consider the whole signal path and limited isolation (i.e. coupling)
between orthogonal channels. [HEX21-D21].
complexity especially at high data rates and bandwidths. Moreover, already in FR2 bands of 5G NR
most transceivers moved to sub-array-based RF architecture, where each sub-array consisting of
multiple phase and amplitude steered antennas was pointed towards beam containing one data stream
(or sometimes multiple streams). With array gain and spatial filtering, the content of each stream was
possible to digitize and process in data converters and digital size with reasonable sampling rate,
resolution, and power consumption. With advances in technology the border frequency between fully
digital and hybrid approach will gradually rise. However, the technology evolution in the area of RF
and mixed-signal processing is moderate (or one could say sometimes even slow) compared to pure
digital (Moore’s law) and therefore when moving towards higher data rates and thus higher carrier
frequencies it is anticipated that sub-array-based RF transceiver architectures are likely to be de facto
also in extreme data rates in 6G.
The sub-array-based approach is sometimes called as a hybrid-MIMO architecture. However, it is only
one specific (but the most straightforward and many times the most practical) case of that. The approach
has been theoretically generalized in [ARA+14]. However, the generalization will lead in large antenna
arrays to implementation complexity especially in RF wiring and interconnects that is not possible to
realize. Only limited scale trials have been demonstrated for this architecture summing signals to
multiple paths only from few antennas [MSH+18]. One must also note that in any hybrid architecture
dynamic range requirements of all RF and analogue blocks change compared to fully digital MIMO.
(NOTE: If multiple signals are transmitted from the same antenna power amplifiers (PAs) and any other
components in the signal chain must be dimensioned taking this additional dynamic range into account
[TAT+16]. The same applies for blocks in the receiver after any signal combining. All signals within
the sub-array are amplified to achieve higher SNR but also at the same time larger dynamic range from
the following components.)
For the reasons above, for any 6G transceivers the simplest sub-array-based hybrid architecture is the
best reference for the architectural studies. The same approach reduces to fully digital approach if the
number of antennas in each sub-array is only one. Naturally, in that case it makes only sense to perform
phase shifting in digital domain removing any analogue phase shifter from the RF parts. The same
model is also applicable for lens antenna-based solutions as the amplification comes from the lens and
opportunities to further steer the beam are limited to activate the appropriate antenna feed for specific
direction [AAK16]. As long as the sub-arrays are combined into one digitized stream in one or multiple
phases the number of antennas is not fundamentally limited to any specific value. Of course, it has direct
impact to the performance partitioning of the receiver or the transmitter impacting power consumption
etc. But combining is feasible and recent implementation shows 384 element phased array consisting
of multiple Radio Frequency Integrated Circuit (RFICs) and modules. Finally, the inverse relation
between antenna array/element gain and beamwidth as shown for example in [RKL+20] equally valid
for antenna arrays and lenses should be noted.
Gant,TX Gant,RX
fs,,DAC DaTX DATX Psat,PA NFLNA DARX DaRX fs,,ADC
Nb,DAC Nb,ADC
(SNDR,DAC) BTX1 BRX1 (SNDR,ADC)
I DF D A DF I
E A
Garr,TX Garr,RX D E
STX1 SRX1
Q DF nant,TX nant,RX DF Q
H(s)fRF,BW
D A
E A D E
BWRF
BTX2 BRX2
I DF D A DF I
STX2
E A
BWRF,ch D E
SRX2
Q DF D A DF Q
E A D E
Figure 4-2: RF transceiver reference architecture for hybrid beamforming based on sub-array per beam
principle. Key parameters for model abstraction are named.
A block diagram of the reference transceiver architecture for 6G analysed in this document is given in
Figure 4-1. It contains key parameters that are needed in the analysis. The highest level of abstraction
provides means to analyse coarsely the link performance, architectural constraints and boundaries of
Dissemination level: public Page 38 / 127
Hexa-X Deliverable D2.3
various technologies using scalar performance values from different entities. The basics of the approach
is described in detail in [TTP17]. The key restriction of the method is ignorance of the frequency
responses in the signal chain and other components. Therefore, all impacts from frequency selective
filtering, beam squint, NLOS channel, etc., should be modelled using implementation margins that are
related to specific implementation or physical boundaries in the particular waveform, HW realization
or both at the same time. However, this is a conventional approach when defining specification limit
for example for minimum signal-to-noise ratio (SNR) for the receiver/detection. However, over the
different generations of wireless systems this approach has proven to be accurate enough for
architecture and coarse performance tuning of the RF transceivers.
The scalar method uses still true bandwidth estimates for the signal (BW RF,ch), radio band (BWRF) and
ADC or digital-to-analogue converter (DAC) sampling rates (fs,ADC & fs,DAC). In the first order, this
could be understood as ideal brick-wall filters with linear phase response, possibly with some scaling.
Fortunately, when designing with actual (complex) waveforms (Ex. in Matlab or other waveform-based
test benches), more elaborate channel models and realistic filters, the same performance parameters
used in the simplified model are applicable as it to represent HW performance in nominal, and possibly
also in other environmental conditions. Of course, those would require more detailed HW simulations
in the background. Just the methods to calculate signal quality using error vector magnitude (EVM) or
any bit/symbol/packet level error rates requires more complex end-to-end simulations using typically
baseband equivalent modelling. Many of the scalar parameters are applicable to those models as well.
If more precision is needed, signal levels including waveform, noise, nonlinearity etc. could be refined
with their frequency responses. However, then the accuracy of the model needs to be significantly better
and for example memory effects of PAs will require additional details of the model and in most cases
also the dedicated HW implemented for the function. The detailed limits of HW become then
complicated and in many cases relationship to system analysis blurry. The most important aspect
however is to consider the absolute power levels when operating the models in link level analysis or
transmitter or receiver specifications as was done in [TTP17]. The impact of the absolute signal level is
introduced for example in Section 3.2 of [HEX21-D22]. The dynamic range depends on many aspects
like signal (i.e., thermal noise) bandwidth, link range, gain control both in transmitter and receiver,
number of antennas and its relation to linearity.
The lack of proper modelling of signal level dependence is one of the key reasons for misunderstanding
the RF constraints and opportunities. This is natural as making a detailed and accurate model requires
detailed understanding of the HW implementation boundaries and their relations to HW architecture.
Modelling those properly is a complex and tedious task and requires a lot of expertise. Therefore,
models are often developed only to certain simulation cases with a smaller validity range. This is
acceptable if the validity of the models is thoroughly understood and obeyed in the analysis.
For additional details of implementation and technology aspects reader is recommended to familiarize
with SoA publications from IEEE Solid-State Circuits, Microwave Theory and Technology, and
Antennas and Propagation Societies. Overviews can be found for example from [PAB+20] and
[GFF+21].
Table 4-1: Common RF and waveform parameters (assumptions) used in the example analysis.
Parameter symbol example alternative
User data rate with protocol overhead 𝑅𝑢 100 Gbps 1Tbps
coding rate 𝑅𝑐 5/6 uncoded
carrier frequency 𝑓𝑐 140 GHz 300 GHz
In this example, two basic waveforms are used i.e., single carrier QAM and OFDM with QAM
subcarriers. As protocol overheads are not counted for, the numbers cannot be compared in the case of
OFDM directly to achievable 5G NR data rates where overhead is quite significant.
First the bandwidth (BW) and the (data converter) sampling rate vs. data rate will be defined for various
modulations. As BW and sampling rate are somewhat different, they need to be treated separately.
Required RF bandwidth reserves some guard bands to enable channel selection and defines then channel
raster. Generally, it can be defined for any data rate in relative scale. For single carrier (SC) modulations
pulse shaping filter defines the additional BW. A typical number of roll-off r=0.35 is used here. In case
of OFDM the filtering shape of the modulation is sharp but for analogue/RF filtering typical 10% guard
band as in 3GPP is adopted. When calculating total spectral efficiency this must be considered. Thus,
any fixed data rate can be defined for the required bandwidth vs data-rate (defined here as bandwidth
relative to the user data-rate, 𝐵𝑊/𝑅𝑢 ) with different modulations and coding rates. For SC waveforms,
the required relative BW can be written as
1 (4-1)
𝐵𝑊𝑅 (𝑆𝐶) = ,
(1 + 𝑟)𝑅𝑐 𝑀
where 𝑟 is the factor for raised cosine filtering, 𝑅𝑐 coding rate and 𝑀 modulation index. Similarly, for
OFDM the relative BW is
1 (4-2)
𝐵𝑊𝑅 (𝑂𝐹𝐷𝑀) = ,
(1 + 𝐺)𝑅𝑐 𝑀
where 𝐺 is the relative (nominal) guard band between the channels. Based on the values given above,
the relative BWs for fixed data rate are defined in Table 4-2, whose content is illustrated in Figure 4-3.
Table 4-2: Required bandwidth for fixed data rate using different modulations with guard bands.
𝑀 𝐵𝑊𝑅 (OFDM) 𝐵𝑊𝑅 (SC)
Code rate 𝑅𝑐 5/6 1/2 1/3 5/6 1/2 1/3
BPSK 1 1.32 2.20 3.30 1.62 2.70 4.05
QPSK 2 0.66 1.10 1.65 0.81 1.35 2.03
16-QAM 4 0.33 0.55 0.83 0.41 0.68 1.01
32-QAM 5 0.26 0.44 0.66 0.32 0.54 0.81
64-QAM 6 0.22 0.37 0.55 0.27 0.45 0.68
128-QAM 7 0.19 0.31 0.47 0.23 0.39 0.58
256-QAM 8 0.17 0.28 0.41 0.20 0.34 0.51
4.1.1.2 Sampling rate and resolution requirements of data converters ADC & DAC
Nyquist rate for baseband (BB) in-phase and quadrature (IQ) signals is two times the baseband BW i.e.,
equivalent to RF BB. To recover also clock, the phase information and sample rate twice at the RF BW
(four times BB BW in both I & Q channels) are required. Therefore, ADC sampling rate is defined
simply as twice the coded data rate using coding rate and modulation index as
2 (4-3)
𝑓𝑠 = 𝑅 .
𝑅𝑐 𝑀 𝑢
For 100 Gbps data rate the sampling rate (Gsps) can be defined for different modulations and coding
rates according to Table 4-3. Note that the sampling rate requirement in Table 4-3 are with assumption
that only the actual signal is quantized/sampled. However, in many cases the requirement is higher, for
example, when using correction techniques such as digital predistortion (DPD) to correct transmitter
nonlinearity. In case of transmitter linearization, so called feedback or measurement, receiver requires
often higher bandwidth (1.5-3 times) [LWZ20] to provide enough bandwidth for linearization,
depending on the target linearization performance.
Table 4-3: ADC sampling rate (Gsps) for 100 Gbps data rate for different modulations and coding rates.
Code rate 𝑅𝑐 5/6 1/2 1/3
BPSK 240 400 600
QPSK 120 200 300
16-QAM 60 100 150
32-QAM 48 80 120
64-QAM 40 66.7 100
128-QAM 34.3 57.1 85.7
256-QAM 30 50 75
Moreover, it would be desirable to use the same ADC resources for both observation receiver and actual
data receiver, which expands the sampling requirements [Ter22]. DAC sampling rates can be defined
equivalently. Note that similar bandwidth extension is needed also in the entire RF signal path. This is
because the predistorted signal bandwidth is larger than the actual signal bandwidth before the PA.
Table 4-4: ADC dynamic range and sampling rate requirements for various combinations resulting to 100
Gbps data rate with 5/6 coding rate.
OFDM SC (r =0,35)
RF RF
BW fs SNDR BW fs SNDR
BW BW
GHz GHz ADC GHz GHz ADC
factor factor
16-QAM 0.33 33 60 44.6 16-QAM 0.41 40.5 60 41.7
64-QAM 0.22 22 40 48.4 64-QAM 0.27 27 40 45.7
256-QAM 0.17 16.5 30 55.7 256-QAM 0.20 20.25 30 53.1
To evaluate ADC requirements and power consumption also required resolution for different
waveforms needs to be evaluated based on their minimum SNR requirements also considering some
margin for clipping and the fact that ADC cannot dominate the total noise figure (NF) of the receiver
and is preferably almost negligible. Therefore, ADC quantization noise must be roughly 15 dB below
the amplifier input referred noise from RF and analogue components giving only 0.1dB penalty to input
referred noise level. However, this is one factor that significantly increases ADC resolution requirement
compared to classical digital only noise estimates for digital modulations. Also, some margin needs to
be reserved for gain control etc. Details of these analysis are given in public deliverable [HEX21-D12].
Based on that resolution and sampling rate for different modulations are defined for small coding rate
or 5/6 i.e., some coding will be used but as little as possible to minimize bandwidth overhead. The
results for selected modulations are given in Table 4-4 and illustrated in Figure 4-4. It is assumed that
transmitter and receiver split the EVM budget equally as is quite common practice in very high data
rate scenarios in 5G NR.
Figure 4-4: ADC dynamic range and sampling rate requirements for various combinations resulting to
100Gbps data rate.
Figure 4-5: ADC dynamic range (SNDR) and power consumption (P) is compared to ADC dynamic
range requirements for the OFDM and CW waveforms (solid oval). Some commercial state-of-the-art
ADCs with high resolution and sampling rate (dashed ovals) have been added to complement the data
from [Mur22].
ADC models based on data in [Mur22] with selected commercial SoA converters can be further refined
based on the formulas presented by Walden [WAL99] and Schreier [ST05] with boundaries defined in
[Mur22]. However, drawing conclusion simply by using the formulas is not straightforward and
dependencies between dynamic range and power consumption cannot be drawn just based on
boundaries. Data in [Mur22] is captured from various scientific publications that are prototypes to show
feasibility or new ideas and they may not necessarily include all the components comprehensively like
clocking needed in complete ADC sub-systems. Therefore, comparison requires also perspectives to
data available publicly from commercial ADCs. It is well-known that ADCs are somewhat more
complex circuits compared to corresponding DACs having similar sampling rate and resolution.
Therefore, DACs are known to consume less power and if ADC is feasible, it should be possible to have
also corresponding DAC available. As DAC power consumption and performance depends also on
(analogue) load impedance it is driving. Therefore, it is more difficult to make similar predictions of
figure-of-merits (FoMs) for DACs as presented for ADCs in [Mur22].
Figure 4-6: EVM for 5GNR and SC (r=0.35) 64QAM modulations from the 1 dB compression point of the
receiver (or transmitter) where non-linearity is dominated by the 3rd order term.
For a transmitter, definition of nonlinearity requirements depends not only on the crest factor but more
complex non-linear behaviour (depending on the PA class, technology and design aspects), opportunity
to clip the highest peaks somewhat even in digital side reducing PAR, and possible linearization. Those
have direct impact how close to saturate output power (Psat) one can drive the PA. In order to make first-
order estimate on the relation of Psat to a class-A type of PA without linearization or digital clipping one
can estimate necessary backoff from probability of the peaks. In this case, BO for different modulations
is determined based on 0.1% peak probability as shown in Figure 4-7 leading to 5.6-8.5 dB BO
depending on the modulation. This is the major factor for the link range when evaluating the PA output
power for the same PA. Values are summarized in 4.1.1.5.
• Front-end losses between antenna elements and active circuitry (PA & LNA) have been
neglected (‘zero loss’).
• Tx and Rx are assumed to have separate antenna arrays i.e. no transmit-receiver switch (TRx
switch) is included that has inevitable loss both in transmit and receive directions.
• Analysis is based on scalar values and flat filtering responses that are not feasible in real
implementations. Possible effects of physical bandpass filters and resonators and beam squint
due to phase shifter in typical arrays have been ignored.
• Perfect beam alignment is assumed (no beam pointing error estimated).
• Only LoS channel model is adopted here but more sophisticated models can be applied later.
• Non-linearity aspects would require more elaborate models and analysis both in transmit and
receive ends of the link.
• Only single-carrier link is analysed here but mechanisms to enhance this to multi-carrier /
carrier aggregation schemes is quite straightforward.
• Power consumption and therefore also heat, heat management, and form factor may and likely
will be a limiting factor in cases when a large number of antennas is required. It is possible to
enhance the analysis to that direction but it is not done yet.
• This analysis gives guidelines for more detailed power consumption estimates but details are
not given in this section. Some modelling aspects with a concrete example have been given in
[HEX21-D21] and in [HEX21-D22].
• Models are specifically done for 100-300 GHz range but they give also guidelines over a
broader frequency range. However, especially, frequencies below 10 GHz should not be
directly compared using the models before further characterization of their validity at those
frequencies.
• Model validity needs to be understood and acknowledged.
Therefore, all results given in the following can be considered as fundamental technological
boundaries that cannot be exceeded but definitely there will be additional sources of loss either
deteriorating the communications range or causing a need to add parallelism (i.e. number of antenna
elements) to compensate the impact. Analysis is made for phased arrays but can be quite
straightforwardly enhanced to other antenna types and solutions.
Despite of all disclaimers as above the model and results are valid in the relative scale and they give
also a coarse ballpark for absolute values that is a core requirement for any implementation and
architectural estimate to get start with. The proposed method gives opportunities to narrow the gap
between conventional communications system design and HW architecture and implementation design
process at an early stage. However, it will not remove the need to mutually optimize and refine the
system based on the same principles when standardization and product development are considered.
The proposed model and the method give means to:
• model system boundaries over a large range in RF spectrum based on the data rate requirement.
• model with key parameters including noise and output power over broad range of radio
spectrum and get guidelines for preferred range of spectrum for particular application with
certain data rate requirement.
• utilize highly systematic and solid approach over large frequency range.
• enhance modelling to different sub-topics and bring more accurate models for analysis.
• give concrete inputs to model system with more accurate models in waveform simulations (Ex.
MATLAB), RF/mixed-mode hardware simulations (Ex. Cadence, ADS, CST, HFSS) and
digital hardware simulations (Ex. System-C, Verilog, VHDL).
• bring abstracted results of the refined simulations back to link level analysis from all domains
involved.
• co-design system at different hierarchical levels.
• facilitate and maintain consistent understanding of the key aspects need to be understood from
the system.
The results in the following examples are compared to each other using the required number of antenna
elements in a phased array into account. It assumes that PA output power, Rx NF or any other parameter
will not change when parallelism in HW is added/enhanced. That is not generally valid and further
refinements are needed based on the studies on antenna isolation, its impact to power amplifiers,
examples of complete phased array transceiver implementations etc. are somewhat included. However,
such examples are always based mostly on specific scenario/HW involved to analysis and their
generalization is not easy. Therefore, those can and should be considered only when architectural design
on a specific HW platform will be started. One just needs to understand that some implementation
margin needs to be accepted at a later phase.
OFDM and SC modulations for the same bandwidth might sound small but that is an indicative number
for differences on RF power consumption including PAs. As energy efficiency is critical now and it
will be even more for 6G this is a significant difference. Saving 10% by some other means would be
very difficult if not possible in apples-to-apples analysis.
Detailed conclusions will be left for the future work, but charts given here with explanations in the
charts themselves and figure captions should give good idea on some of the core aspects here.
The number of antennas, as shown in Figure 4-9, naturally needs to be increased as loss is increasing.
Light blue box indicates the range when relative bandwidth is 20% or below. Also, 75dBm line to show
EIRP limit is shown. It will not be a problem in this case.
Figure 4-8: Number of antennas, EIRP and relative bandwidth as a function of frequency for 100 Gbps
target data rate with different modulation orders and link distance. Legends and axis to which each of the
curves points are shown in figure (a).
75 dBm
Figure 4-9: Number of antennas, EIRP and relative bandwidth as a function of frequency for 100 Gbps
considering 10 dB and 20 dB to describe pathloss of NLOS radio channel.
𝑁 𝑓𝑜 𝛼𝑧,𝑛
𝑁𝑟𝑒𝑓 ∏ 2𝑛=1 (1 + ( ) )
𝑓𝑐 𝑓𝑧,𝑛 (4-4)
𝑆(𝑓𝑜 ) = 10 10 ( ) ,
𝑓𝑝𝑙𝑙 𝑀 𝑓𝑜 𝛼𝑝,𝑚
∏ (1 + ( ) )
( 𝑚=1 𝑓𝑝,𝑚 )
where 𝑓𝑜 is the offset frequency from the carrier frequency 𝑓𝑐 , 𝑁𝑟𝑒𝑓 is the reference phase noise level at
PLL operating frequency 𝑓𝑝𝑙𝑙 , fz,1 .. fz,N are the corner frequencies of the zeros, fp,1 .. fp,M the pole
frequencies, αz,1 .. αz,N are the orders of the zeros, and αp,1 .. αp,M the orders of the poles, respectively.
2
𝑓
The frequency multiplication part is described by term (𝑓 𝑐 ) that shifts the whole phase noise curve
𝑝𝑙𝑙
up 20 dB / decade. In the expression (4-4), the flat part of the phase noise is done by the last zero, that
turns the phase noise spectra into flat. Because the flat part is the most important phase noise contributor
in very wideband systems, that 6G systems generally are, (4-4) can be also rewritten in another form as
𝑓𝑜 𝛼𝑧,𝑛
2
𝑁𝑟𝑒𝑓 ∏𝑁
𝑛=1 (1 + (
) ) 𝑁𝑡ℎ,𝑝𝑙𝑙
𝑓𝑐 𝑓𝑧,𝑛 (4-5)
𝑆(𝑓𝑜 ) = ( ) 10 10 + 10 10 ,
𝑓𝑝𝑙𝑙 𝑓𝑜 𝛼𝑝,𝑚
∏𝑀
𝑚=1 (1 + (𝑓 ) )
( 𝑝,𝑚 )
where the additional term 𝑁𝑡ℎ,𝑝𝑙𝑙 describes the (white) noise level at the PLL output. Using (4-5) is
more systematic as it helps to derive noise calculations for the LO chain with potential buffer amplifiers
and frequency multipliers. It also gives easy parameter for analysing the phase noise performance with
very wideband systems as single parameter tells the noise level that is then scaled up by 20 dB/decade
(memory-dependent filter parts from equation (4-5) can be then neglected with analysing maximum
achievable performance with perfect cancellation of the memory-dependent phase noise terms). Also,
(4-5) needs generally less parameters as the flat part are generated with one additive term instead of
separate zero expression.
(a) the example TI LMX2596 PLL with different (b) the PLL phase noise spectra at 15 GHz
PLL output frequencies all scaled to 150 GHz and compared to the derived models 1 and 2.
An example parametrization for both described models is given as in Table 4-7 and Table 4-8 and the
corresponding curves at 15 GHz in Figure 4-11 (b). As one can see, both models give the same curve
and their behaviour approximate also well the original PLL data given in the datasheet.
Table 4-7: Model parametrization for phase noise Table 4-8: Model parametrization for phase noise
model 1 given in (4-4). model 2 given in (4-4).
fvco 15 GHz
fvco 15 GHz
𝑁𝑟𝑒𝑓 −62 dB
𝑁𝑟𝑒𝑓 −62 dB
𝑁𝑡ℎ,𝑝𝑙𝑙 −149.8 dB
𝑛, 𝑚 𝑓𝑧,𝑛 𝛼𝑧,𝑛 𝑓p,𝑚 𝛼p,𝑚
𝑛, 𝑚 𝑓𝑧,𝑛 𝛼𝑧,𝑛 𝑓p,𝑚 𝛼p,𝑚
1 25 k 1.3 1 1.1
1 25 k 1.3 1 1.1
2 21 M 2.8 650 k 3
2 N/A N/A 650 k 3
For fair comparison, the comparison among between the proposed model parametrization to the existing
upscaled 3GPP reference models is illustrated in Figure 4-12 (a) 150 GHz and (b) 300 GHz. As one can
observe, different models have rather significant differences in the overall characteristics.
Figure 4-12: Comparison of the derived phase noise models against 3GPP. The models 1 and 2 given in
this document gives the equivalent phase noise spectra.
4.2.1.2 System phase noise behaviour and limitations on the achievable performance
and bandwidth
Figure 4-13 illustrate the SSB phase noise spectra together with different system level aspects that have
impact on how different regions matter on the signal phase variations from the system perspective. Note
that the spectrum is given in logarithmic frequency axis, as typically made, to highlight different
regions. Waveform selections and frame structure has impact on the overall performance of the system
under phase noise. In particular, the phase noise impact on different waveforms depends on how the
phase is synchronized in the receiver and how the data and pilot sequences are distributed across the
time domain signal with respect to the phase jitter characteristics.
In multi-carrier waveforms, the lower offset frequency components of the phase noise cause so-called
common phase error (CPE) due to long symbols that can be compensated in the receiver by proper
phase noise tracking reference signals. The higher offset frequencies and the flat part of the phase noise
then causes inter-symbol interference that generally sets limitation on the achievable SNR. In single
carrier waveforms, the symbols are short in time domain which means that the potential phase drift over
a single symbol is mostly dominated by the flat phase noise if fast symbol rate is assumed. The potential
impact of the lower offset frequency components of the phase noise depends also on the signal length,
i.e., the speed of the phase synchronization in the receiver. These two (minimum frequencies that
matters for the system and the maximum frequency, i.e., symbol rate), which are defined as 𝑓1 and 𝑓2
for calculating the overall jitter from the phase noise in [HEX21-D22, Eq. (3-14)]. In addition to
waveform level parameters, also the environment has impact on the phase variations of the received
signal. For example, mechanical vibration causes random phase variations for the link in terms of
Doppler effect, but also part of the vibration may be conducted to the synthesizer itself and cause phase
to vary [HNH+07, HNH09]. These may be important especially in industrial applications, for example
in factory automation, drones, etc. The impact of the mechanical vibration on the phase noise spectra is
generally experienced in lower offset frequencies, but in the very high centre frequencies, their
harmonic impacts may be seen in even tens of MHz offset frequencies.
Figure 4-13: Illustration on the different aspects in system level that have impact on the importance of
different regions on the phase noise characteristics in different offset frequencies.
In the literature [DPS20, CHK+17], phase noise with high offset frequencies (wideband phase noise) is
identified as one of the most critical performance bottlenecks especially with very wide signal
bandwidths versioned in 6G. However, the impact of the operating frequency has not been analysed in
detail in the open literature in a systematic way. In the following, the limitation of the phase noise on
the the achievable bandwidth and SNR in different centre frequencies is analysed. The analysis assumes
that the flat part of the phase noise dominates. This is the case if the signal structure is generated in such
a way that the occurrence frequency of the phase noise tracking signals is more than the corner
frequency where the phase noise spectra turn to flat. In practice, this should be at least tens of MHz that
is certainly possible at least with very wideband single-carrier waveforms. This basically means that
when 𝑓𝑝𝑡𝑟𝑠 > 2𝑓𝑧,𝑁 in model 1(4-4) is chosen, the memory-dependent part of the phase noise is
assumed to be compensated. In model 2 (4-5) this then means that the remaining part of the phase noise
is only the flat term of the phase noise. By doing this, the remaining standard deviation of the phase
noise can be written as
2 𝑁𝑡ℎ,𝑝𝑙𝑙
fc (4-6)
𝜎𝑝𝑛 = √𝐵 ⋅ ( ) 10 10 ,
f𝑝𝑙𝑙
where 𝐵 is the bandwidth of the signal (symbol rate). Thus, the phase noise-limited SNR can be written
as
2 𝑁𝑡ℎ,𝑝𝑙𝑙
fc
𝑆𝑁𝑅𝑃𝑁 = −20 log10 (√𝐵 ⋅ ( ) 10 10 ) (4-7)
f𝑝𝑙𝑙
= −𝑁𝑡ℎ,𝑝𝑙𝑙 − 10 log10 (𝐵) − 20 log10 (𝑓𝑐 ) + 20 log10(𝑓𝑝𝑙𝑙 ).
Note that 𝑆𝑁𝑅𝑃𝑁 in the equation mimics the SNR budget given only for the phase noise. In addition to
this, the overall SNR is naturally also limited by the rest of the RF nonidealities that set their own limits.
In the end this means that the achieved overall SNR is always less than the SNR budget given for the
phase noise. The equation can be also defined to calculate maximum achievable bandwidth for a given
centre frequency and SNR budget as
𝑓𝑣𝑐𝑜 2 (−𝑁𝑡ℎ,𝑝𝑙𝑙10
−𝑆𝑁𝑅𝑃𝑁
) (4-8)
𝐵=( ) 10 .
𝑓𝑐
An important note from the expressions (4-7) and (4-8) is that the bandwidth dependency is 10
dB/decade, and the carrier frequency impact is 20 dB/decade. This generally means that it is doubling
the centre frequency for a given fixed SNR target, decreases the achievable bandwidth to be four times
less. As these both (bandwidth and centre frequency) increases when going towards 6G systems, this
result is of a particular importance.
A single radio link generally includes at least a contribution from two LOs, i.e., phase noise from
transmitter and receiver. Usually, in cellular links, it is often assumed that the mobile terminal has worse
performance. The overall link phase noise can be calculated as the root mean square (RMS) sum of the
noise contributors (Tx and Rx) in the best possible case if receiver own thermal noise and multipath
propagation would not limit the performance of the phase noise compensation. To get the values of the
previous analysis into context, four downlink scenarios in different carrier frequencies and bandwidth
are analysed. The scenarios with the parameters are given in Table 4-9.
Table 4-9: Simulation parameters for link phase noise analysis.
𝑓𝑝𝑙𝑙 15 GHz
The results in all four scenarios are illustrated as contour plot in Figure 4-14 (a)-(c). The results are
illustrated as contour-plots where each curve in the figure describes the bandwidth-carrier frequency
pair for a given phase noise limited SNR. One should note that in lower frequency communication
systems with less than 1 GHz of signal bandwidth, the remaining phase noise is rarely a performance
limiting factor for the system and hence the SNR budget given for the phase noise is usually in a range
of 30-40 dB at minimum and hence the overall system performance is limited by other factors. To give
few numerical examples in the figure: at 𝑓𝑐 = 30 GHz, the maximum bandwidth for 30 dB 𝑆𝑁𝑅𝑝𝑛 in
the worse scenario (10 dB worse mobile terminal PLL) is 24 GHz, meaning that phase noise is not a
limiting factor. However, for the same scenario the maximum bandwidth would be around only 950
MHz at 𝑓𝑐 = 150 GHz and 240 MHz at 𝑓𝑐 =300 GHz. These numbers are serious limiting factors. Same
scenarios for 20 and 25 dB 𝑆𝑁𝑅𝑝𝑛 targets would be 9.5 and 2.4 GHz for 𝑓𝑐 =150 GHz and be 2.4 GHz
and 750 MHz at 𝑓𝑐 =150 GHz. Moreover, even with symmetric link, assuming that both link ends have
extremely good PLL performance, the phase noise is very drastically limiting the achievable signal
bandwidth in high carrier frequencies.
(a) Single link end (b) Whole link with equal PLL performances in
Tx and Rx
(c) Whole link with 6 dB worse (d) Whole link with 10 dB worse PLL
performance
Figure 4-14: Phase noise-limited SNR. The different curves in the figure indicate the different phase noise
limited SNRs.
Note that the analysis above assumed perfect cancellation of the lower frequency components of the
phase noise. Hence, as some room must be left also for the compensation performance, the performance
is decreasing further. Moreover, also the other RF nonidealities are contributing to the overall effective
received SNR, the numbers given as an example gives rather pessimistic view for achieving very wide
signal bandwidth in high frequencies at least with higher order modulations. However, the authors think
that the problem can be solved in multiple domains and gives good research directions for further
analysis and future work. First, it is evident that more SNR budget must be given for the phase noise at
sub-THz region than in the lower frequencies. This is also natural as the target spectral efficiencies and
lower order modulations are often versioned to be used that can survive with less signal to noise ratio.
In waveform design, using analogue carrier-aggregation schemes with proper filtering could help to
reduce the impact. There the trade-off is then to choose wide enough carriers to coop with lower offset
frequency components of the phase noise but narrow enough carriers to tolerate the flat part of the phase
noise. Such a channelized processing is also useful to generate and analyse signals by using multiple
lower-bandwidth DACs and ADCs. Furthermore, different filtering schemes in OFDM-like waveforms
may help to reduce the inter-symbol-interference caused by the wideband phase noise. Also, using
analogue filtering to reduce phase noise in the PLL output could help to reduce the phase noise in some
scenarios, but the filtering must be designed together with desired LO tuning range. Especially in those
scenarios, understanding the noise budget in the multiplier chain and the overall noise level of the
multipliers themselves would be critical.
Figure 4-15: Saturated output power models based on published PAs on linear frequency scale for
different semiconductor technologies.
The InP-based power amplifier models are shown in Figure 4-15 (d) and Figure 4-16 (d). The dominant
process node for the best-performing InP-based power amplifiers is 250 nm. The 130 nm node is used
at the highest frequency PAs. The modelling of the saturated power of InP-based PAs requires two
linear equations on both linear and logarithm scales. The threshold frequency of the modelling lines is
230 GHz since the slope changes from a flat (-1) to a very steep (-111) declining performance on the
logarithm scale. The InP power modelling has similar curves regardless of the frequency scale (linear
or logarithm).
Figure 4-16: Saturated output power models based on published PAs on logarithm frequency scale for
different semiconductor technologies.
Figure 4-18: Performance of the 290 GHz SiGe amplifier, together with common memoryless amplifier
models with fitted parametrization.
For each measured power level, an individual two-port calibration was performed. The simulated and
measured power compression curve of the amplifier is shown in Figure 4-18 (a) and (c), and the phase
plot over power levels is shown in Figure 4-18 (b) and (d) with black dots, respectively. The simulations
are circuit-level simulations of the amplifier design. The measured AM-AM curve clearly follows a
typical amplifier's compression curve. The variation can be seen from the AM-PM curve in the range
of 20 degrees which corresponds to 55 um distance at 300 GHz. This variation is difficult to avoid when
multiple mechanical probe adjustments are performed between calibrations and measurements at the
300 GHz range. Hence, the presented AM-PM here is having some uncertainty that cannot be avoided.
Based on the measured and simulated AM-AM and AM-PM at 290 GHz, all four typically used PA
compression models were fitted to the datasets. The founded parameters are gathered to Table 4-11 and
the corresponding curves are plotted to Figure 4-18. Note that the measured plots are purposely drawn
over the fitting region to understand the behaviour outside the measurement range. In the models,
amplitude is given normalized voltage (amplitude) against 1 Ω reference (can be used to calculate power
directly as |𝑥|2 ), while AMPM parameters are given in degree scale. Figure 4-18 (c), Figure 4-18 (d)
also contain the relative phase shift of the amplifier, which have been normalized away in the
parametrization given in Table 4-11. Among all the models, Rapp model with gives the best
compromise since has meaningful behaviour also outside the measurement range. This is important
thing to keep in mind when performing system level simulations, as the model validity can be usually
guaranteed only within the region to which it has been fitted. Both parametrizations, simulated and
measured, for the Rapp model are given and may be used for system level simulations. Note that the
parametrizations given here are not normalized and hence they also include the output power delivery
capability of the amplifier under test.
Table 4-11: Parameter set for common memoryless amplifiers models fitted for circuit level simulations
(sim) and measurements (meas) of 290 GHz SiGe amplifier.
AM-PM
Model (par). AM-AM (par) AM-AM (numeric)
(par)
AM-PM (numeric)
Rapp (sim.) {gr,xsat,s} {4.708,0.663,1.603} {α,q1,β,q2} {−740.2,1.945,0.298,1.797}
Rapp (meas.) {gr,xsat,s} {2.55,0.605,1.296} {α,q1,β,q2} {−2.84 · 103,3.00,0.239,2.640}
Ghorbani (sim.) {a1,a2,a3,a4} {29.07,1.454,37.22,−0.1057} {b1,b2,b3,b4} {1.235,0.0018,−0.461,−111.856}
Ghorbani (meas.) {a1,a2,a3,a4} {7.628,1.296,8.742,−0.2037} {b1,b2,b3,b4} {0.4163,0.0013,−0.994,−157.4}
White (sim.) {aw,bw,cw,dw} {1.008,6.604,−0.733,0.760} N/A N/A
White (meas.) {aw,bw,cw,dw} {0.146,−3.456,20.159,1.156} N/A N/A
Saleh (sim.) {as,bs} {4.58,10.32} {cs,ds} {−132.85,0.673}
Seleh (meas.) {as,bs} {2.52,5.16} {cs,ds} {−79.85,−1.37}
base-band output signal. In this simulation approach, input and output mismatch network effects in the
frequency domain are included.
The thermal noise, drain-current non-linearity and memory effects are introduced in the intrinsic circuit
model, as listed below.
• Thermal noise: Random thermal agitation of charges inside a conductor produces thermal noise
which is proportional to the temperature. This current voltage is given by Planck’s blackbody
4ℎ𝑓𝐵
radiation law [POZ11], 𝑖𝑛 = √ ℎ𝑓 , where h is Planck’s constant, k is Boltzmann’s
𝑅(𝑒 𝑘𝑇 −1)
constant, T is the temperature, B is the bandwidth of the system, f is the centre frequency of the
bandwidth, and R is the device resistance.
• Drain current non-linearity: Non-linearities from drain current and gate–source capacitance can
give rise to undesirable effects such as gain compression and the generation of spurious signals.
Since the non-linearity effects of the drain current are much larger than those from the gate–
source capacitance [VRD03]. Here, only the drain current non-linearity is considered in the PA
simulation. The nonlinear current source controlled by the intrinsic gate voltage (𝑣𝑔𝑠𝑖 }) can be
modelled by the hyperbolic tangent (tanh) function [AZR92]: 𝑖𝑑 = 𝑖0 (1 + tanh(𝑝1 (𝑣𝑔𝑠𝑖 −
2
𝑣𝑝𝑘 ) + 𝑝2 (𝑣𝑔𝑠𝑖 − 𝑣𝑝𝑘 ) + 𝑝3 (𝑣𝑔𝑠𝑖 − 𝑣𝑝𝑘 )3 )), where 𝑖0 scales the current. 𝑣𝑝𝑘 and𝑝1,2,3 can be
extracted by fitting 𝑣𝑔𝑠𝑖 -𝑖𝑑 curves.
• Memory effects: The output of a PA with memory effects is not only determined from the
instantaneous input, but also the input past and/or the system state [OML09]. Here only the
memory effects from a linear gate–source capacitance are considered. In this case, the gate
1 𝑡
voltage is given by [VRD03], 𝑣𝑔𝑠𝑖 (𝑡) =
𝐶
∫−∞ 𝑖𝑔𝑠 (𝑡 ′ )𝑑𝑡 ′ ), where the voltage at time t is
𝑔𝑠
proportional to all prior gate current values, not just to the instantaneous value. Self-heating
and trap-related mechanisms can also introduce memory effects. However, those are not
considered in this work since some published measurement results show their minor
dependence on memory effects [BG03, BP04].
Figure 4-21: Measured and circuit simulation AM-AM (a, c and e) and AM-PM (b, d and f) results of the
GaN with different bandwidths (10 MHz, 40 MHz and 160 MHz).
Figure 4-23: AM-AM (a, c, e) and AM-PM (b, d, f) behavioural modelling results based on measured and
simulation results with different bandwidths (10 MHz, 40 MHz and 160 MHz) of the GaN.
In summary, RF equivalent circuit simulation method is demonstrated for generating PA baseband input
output signal. The method including both memory effects and nonlinear effects is validated by GaN PA
measured results. The NMSE between the measured results and circuit simulation results of the GaN
PA is lower than -15 dB with bandwidth less than 160 MHz at 2.14 GHz. It is important to note that
this method is not limited to GaN PA but also can be applied to any kinds of PAs. Similar PA modelling
can be done for higher frequencies but needs appropriate measurement data for validation. More
experimental data is needed to further validate the presented method for generating input-output
baseband data.
solution together with reuse of already existing deployments will probably lead to that the interaction
between legacy systems such that dense D-MIMO deployments can co-exist with legacy macro sites.
In dense urban area scenarios that requires high user experience data rates and high traffic area capacity,
e.g., public squares, stadiums, airports, etc., D-MIMO networks can enhance performance with
capacity-based features, e.g., higher spectral efficiency at frequency range 1 (FR1) bands; in
conjunction with dense deployment on areas where capacity cannot be met by macro sites in an efficient
way. It will provide robust deployments that carry decent traffic regardless of site location; enable
efficient multi-user MIMO with coherent joint transmission with interference suppression capability.
One can apply flexible total radiated power (TRP) selection/clustering based on the scenario
requirements. Last but not least, outdoor D-MIMO installations need to have a low visual footprint due
to invisible, and aesthetic impacts. With radio stripes, it is possible to hide the installation in existing
construction elements in the environment.
For the indoor deployment scenarios Figure 4-24 e.g., factories, warehouses, offices focusing on dense
machine-type communication (MTC) and extended reality (XR), while D-MIMO not only provides
traffic area capacity, but also enhances reliability and resilience aspects. The line-of-sight probability
for D-MIMO is very high which makes it suitable for deployment on very high frequency bands where
radio propagation makes it challenging to provide a robust access link for mobile users. Phase coherent
joint transmission is hard to realize at frequency range 2 (FR2) bands, but non-coherent joint
transmission can help creating reliable access link due to macro diversity while it can still provide high
rates thanks to high bandwidth despite its low spectral efficiency performance. Layer 1 (L1) mobility,
e.g., combined cell type of radio units (RUs), which could be transparent to UE, can provide reliable
and resilient mobility without frequent handovers. Besides, distributed environment inherently
enhances network’s sensing and positioning capabilities. In some scenarios, coverage from existing
macro deployments may not exist. In these cases, D-MIMO networks can be used for data boosting but
may require support for standalone D-MIMO operation.
The introduction of D-MIMO to a network means, that i) the network will have to support multi-
connectivity of a user, ii) the network deployment will be densified and, iii) a much higher percentage
of fronthaul links will be used in the network. These three aspects will not only affect deployment
scenarios, signal processing, scheduling etc., but also affect user- and control-plane security. Multi-
connectivity will need an arrangement for authentication and key agreement between the network and
the user, which supports the connection of a user to multiple radio units. Deriving multiple keys and
managing multiple secure connections may result in some extra overhead on the user side but in
principle can be realized with today’s techniques. D-MIMO schemes that require that the signal
transmitted from multiple radio units combine coherently at the user will need a common encryption
key to be used for the different signal paths. This applies to both user- and control plane traffic. Note
that in the control plane, not only the higher layers but also the MAC layer will have to be secured, so
depending on the realization of the medium access control (MAC) layer in D-MIMO (how centralized
it is) it might be necessary, that the network supports multiple security associations of a user, even if
the higher layer control processing is centralized. To support the higher number of network nodes
implied by D-MIMO, as well as multi-connectivity, it will be favourable to terminate user plane security
more centrally in the network, for example where the application resides. This will reduce the number
of exposed network points and will allow for using common user plane keys. The densification of the
network using D-MIMO can be realized in different ways depending on how the functional split of
radio unit (RU), distributed unit (DU), and centralized unit (CU) is implemented. In either case, it will
be necessary that the nodes closer to the user, will be as cheap as possible (many nodes), and easy to
deploy. Easy to deploy means the fronthaul nodes should support automatic certificate enrolment
procedures such as specified for gNB for example. Based on this automatic enrolment, secure
communication between the fronthaul nodes can be established.
If cheap and easy to deploy is the goal, this typically means that there must be a trade-off with other
metrics. In this case, this means that cheap and easy to deploy nodes will not be as trustworthy as
traditional eNBs or gNBs and not as resistant against physical attacks. The nodes will also be closer to
the user, which makes them easier accessible for potential attackers. The network should not assume
these entities are very trustworthy, and measures should be taken to detect compromise and
malfunctioning and contain the impact of a compromised entity. In this context artificial intelligence
(AI) may potentially be used in the future for monitoring and anomaly detection.
4.3.2 Architectures
D-MIMO systems can be deployed with different optimization goals in mind. At sub-6GHz, the reason
for deploying advanced antenna solutions such as D-MIMO is mainly driven by the need for high-
spectral efficiency. At very high frequencies, the need for distributed antenna systems is driven by the
need to overcome the effects of a less reliable radio channel and to produce reliable communication
links to the user equipment (UEs). Distributed systems allow for an architecture where a serving antenna
is closer to a UE and thus can provide a more reliable link.
Figure 4-25: Illustration of architectural options Figure 4-26: Illustration of distributed MIMO with
for distributed MIMO. key design options.
Higher frequencies and larger bandwidth in the access put higher demands on the front/backhaul for the
distributed antennas as well as the processing of the high bandwidth signals. Processing and
transmission can be categorized in the following ways:
1. Transport media: The backhaul/fronthaul transport media can either be wired or wireless.
2. Signalling: The data transmitted over the media can further be distinguished in being digitally
encoded (e.g., common public radio interface (CPRI) like), or an analogue signal modulated
onto a carrier.
3. Processing: Processing (such as for example beamforming) in different nodes can either be
performed analogue or digitally. A further distinction can be made into centralized or
distributed processing.
4. Transmission: If multiple serving antennas serve a UE jointly, this can either happen through
coherent or non-coherent transmission by the serving antennas.
To enable scalable D-MIMO networks, there is a need to understand how much distribution of access
points is needed and how to address practical issues associated with distributed access points operating
at higher frequency bands. It is also necessary to understand the requirements on transport solutions and
how these requirements can be met. The optimum solution would be phase-coherent joint transmission
with centralized processing, but it is difficult to realize and might not meet the feasibility requirements.
The other extreme would be phase non-coherent transmission with duplicating data in each RU. While
this approach is technically feasible, it is inefficient. Furthermore, a combination of DMIMO and
repeaters (more details in Section 4.3.3.2) will be an interesting topic to study. Not in all cases, one may
require an RU, instead repeaters can extend the coverage in blind spots with relatively low cost.
The trade-off between localized (distributed) and centralized processing raises new challenges on the
fronthaul (FH). The increased number of antennas for massive MIMO drives up the required FH load
proportionally. There have been many FH challenges; some of them originated during 5G and some
new types of FH challenges have been identified for new generations: further increased system
bandwidth, e.g., at higher frequency bands, and multi-band support, more antennas, improved MU-
MIMO performance with more layers, baseband resource pooling and FH statistical multiplexing, Inter-
cell and interference coordination and joint processing, new D-MIMO schemes and deployments, e.g.,
star and serialized (cascaded) topologies. different lower layer split (LLS) options try to balance RU
and FH complexity, cost and capacity for different use cases, while achieving high performance. The
main differences can be summarized how to handle the following functionalities, e.g., where channel
estimation is performed, where the beamforming weights (BFWs) are calculated, what types of
beamforming types to support, which data, e.g., sounding reference signal (SRS), channel information,
BFWs, etc., is transported over the FH interface, and their latency requirements.
Figure 4-27: D-MIMO with wireless fronthaul operating at high bands while access links at low bands.
There are two different ways to provide fronthaul to different RUs: wired through fibre, or via wireless
links. An optimized wired fronthaul can be realized through so-called radio stripes, while wireless
fronthaul/backhaul can be realized efficiently by integration with the access network. Cascading or
serializing the fronthaul connections is also important for indoor installations, e.g., a factory
environment, one fronthaul cable can connect several RUs or Antenna Processing Units (APUs) in some
installations. Cable fronthauling is always preferable if it is easy, however wireless links are also an
option for fronthauling, see Figure 4-25. At higher frequencies (e.g., sub-THz), the spectrum is large
and there will not be much traffic in early phase of the utilization of these bands, even the UE is limited
in how much bandwidth it can process, i.e., there will be available spectrum that can be utilized for the
wireless fronthauling (WFH), Figure 4-27. Since RUs are generally not nomadic, robust fronthaul links
will be easy to maintain. Random blockage due to moving objects can be handled by redundant links.
Moreover, in dense urban areas the backhaul connection for the DU can be provided wirelessly by an
existing macro base station.
Figure 4-28: A wide bandwidth D-MIMO system can be implemented using ARoF and WDM.
structure model, two types of links are considered in this work: strong and weak links. It should be
noted that the strong or weak link can be defined by a threshold of channel strength, e.g., the distance,
or path loss between UE and RU. Strong links must be estimated, whereas weak links can be ignored,
but they have an impact on the pilot contamination. Figure 4-30 (a) indicates the t-th snapshot, while
Figure 4-30 (b) indicates the (t+1)-th after structural modification. In comparison to the previous
snapshot, in these figures the red and black lines represent new and lost connections, respectively.
(a)
(b)
Figure 4-31: An illustration of IAB network (a): A planned IAB network architecture based on
a hierarchy, (b) An IAB mesh network architecture.
“The time it takes to successfully deliver an application layer packet/message from the radio protocol
layer 2/3 SDU ingress point to the radio protocol layer 2/3 SDU egress point via the radio interface
in both uplink and downlink directions, where neither device nor Base Station reception is restricted
by DRX.”
NR frame structure extrapolation
Establishing the time structure of layer-1 transmissions is essential for solving the problem of latency
evaluation. The analysis assumes an NR-like time structure based on OFDM symbols with cyclic
prefix, arranged in slots. Time/frequency resources in a slot are assigned to both data and reference
signalling; hence, a slot is the basic unit of processing time. Slot-based time structure applies to OFDM
and OFDM-like waveforms, such as DFTS-OFDM; however, it can be assumed that other waveforms,
such as pure single carrier, can be represented by a slot-like structure, where data and reference signals
needed for correct reception of data are e.g. time-multiplexed.
As another set of working assumptions include
• Extrapolation of OFDM symbol time durations in NR for larger than NR sub-carrier spacing
(SCS) configurations. Higher SCS are tightly connected to large data bandwidths and high
sampling rates; therefore, the SCS can be seen as a measure of bandwidth and sampling rate
and can in this context be assumed to represent general waveforms, not only OFDM-like
waveforms.
• Adoption of the number of (OFDM) symbols in a slot as a parameter for the parametric latency
calculation. In the results, 1-symbol and 2-symbol slots are assumed. Note that 2-symbol mini
slots are also present in NR.
Numerologies/subcarrier spacings, corresponding OFDM symbol durations and corresponding
maximum supported BW (for standard FFT size of 4096 and not assuming frequency guard bands) are
shown in Table 5-2. As a clarification: a system under analysis can be based on OFDM-like waveforms
with symbols containing (4096 + 288) samples, with 1 or 2 such symbols forming a slot, or by
considering a system that usees the same sampling rate but a different waveform, where samples of this
waveforms form subgroups of (4096 + 288) samples each, and groups formed by 1 or 2 subgroups.
The latency analysis for a given numerology applies to both systems equally, given that they use the
same sampling rate, which is represented by the numerology.
Table 5-2: Extrapolation of NR numerology for FFT size 4096.
Maximum supported BW
OFDM symbol
sub-carrier spacing = sampling rate
duration + CP, in μs
(oversampling ratio = 1)
μ = 0 (SCS = 15 kHz) 71.4 61.44 MHz
μ = 1 (SCS = 30 kHz) 35.7 122.9 MHz
μ = 2 (SCS = 60 kHz) 17.85 245.8 MHz
μ = 3 (SCS = 120 kHz) 8.925 491.5 MHz
μ = 4 (SCS = 240 kHz) 4.4625 983 MHz
μ = 5 (SCS = 480 kHz) 2.2313 1.96 GHz
μ = 6 (SCS = 960 kHz) 1.1156 3.93 GHz
μ = 7 (SCS = 1920 kHz) 0.5578 7.86 GHz
μ = 8 (SCS = 3840 kHz) 0.2789 15.73 GHz
Another option is to further increase the SCS by limiting the FFT size to 1024 or 2048 sub-carriers. The
feasibility of further scaling is depending on the delay spread of the channel. Thus, current, and future
investigation of the channel characteristics with and without beamforming will show whether this is
feasible for an OFDM like system. To compensate for the effect of delay spread, it is also possible to
increase the cyclic prefix (CP) overhead. Current new radio (NR) systems (with normal CP) keep a
constant CP overhead of about 7%. The result of this extrapolation is shown in Table 5-3.
Table 5-3: Extrapolation of NR SCS for FFT size of 1024 and 2048.
sub-carrier spacing in OFDM symbol FFT size 1024 BW FFT size 2048 BW
kHz duration + CP, excluding guard band (~60 excluding guard band (~60
in μs % SCs used) % SCs used)
15 71.4 9.21 MHz 18.4 MHz
30 35.7 18.4 MHz 36.8 MHz
60 17.85 36.8 MHz 73.7 MHz
120 8.925 73.7 MHz 147 MHz
240 4.4625 147 MHz 295 MHz
480 2.2313 295 MHz 589 MHz
960 1.1156 589 MHz 1.18 GHz
1920 0.5578 1.18 GHz 2.36 GHz
3840 0.2789 2.36 GHz 4.72 GHz
7680 0.13945 4.72 GHz 9.43 GHz
Following the choice of RAN latency definition, the evaluation procedure for determining user plane
latency proposed by 3GPP for the NR standard [37.910, Section 5.7] may also be adopted. Under this
definition, the possibility of employing retransmissions to achieve a satisfactory performance level is
accepted. Evaluation of user plane latency between user equipment (UE) and base station (BS) is
illustrated in Figure 5-1, taken from [37.910, Section 5.7.1.1].
Figure 5-2: RAN latency for 0-3 retransmissions, Figure 5-3: RAN latency for 0-3 retransmissions,
slot size of 1 OFDM symbol and UL/DL BS/UE. slot size of 1 OFDM symbol and UL/DL BS/UE.
processing times of 1 slot (high processing processing times of 2 slots (low processing
capability). capability).
Figure 5-4: RAN latency for 0-3 retransmissions, Figure 5-5: RAN latency for 0-3 retransmissions,
slot size of 2 OFDM symbol and UL/DL BS/UE slot size of 2 OFDM symbol and UL/DL BS/UE
processing times of 1 slot (high processing processing times of 1 slot (high processing
capability). capability).
Figure 5-6: Processing time scaling for SCS of 480 Figure 5-7: Processing time scaling for SCS of
kHz. 960 kHz.
Discussion
From the results of the parametric analysis, it is evident that the RAN user plane latency of 0.1 μs is
generally achievable at numerologies 𝜇 = 5 and above (corresponding to bandwidths of 2 GHz and
above), if a) slot size is kept low (1 or 2 OFDM symbols) and b) processing times are on the order of
1-2 slots. In this case, retransmissions are possible. Latency, naturally, degrades with a larger number
of retransmissions and longer processing times (for slot = 2 symbols and processing time of 2 slots a
signal BW of 4 GHz is needed for meeting the latency requirement with 3 retransmissions).
Regarding the second part of the study, the results in in Figure 5-6 to Figure 5-8 show that for the used
slot structure the latency target of 100 μs can be met if the processing times at UE and BS scale down
by about an order of magnitude compared to state-of-the-art for SCS = 60 kHz. In nominal terms, this
means 10 – 50 𝜇𝑠 of processing time per slot at both UE and BS side. Also, it is worth noting that the
latency target in this case can be met with significantly higher processing time than 2 slots. This is
dependent on the SCS but assuming that a maximum of 1 retransmission is used to achieve the required
reliability a processing time of about 15 µs is sufficient.
For finding the ultimate answer on whether meeting the proposed RAN latency requirement is feasible,
it is necessary to find out reasonable UE and BS processing times, considering the projected
computational capability scaling, I/O speeds, memory access speeds and power consumption
constraints.
RF impairment PA nonlinearity
PA model memoryless GaN
PN model Model 2 described in Section 4.2.1.1 scaled to 140 GHz
Carrier frequency 140 GHz
Subcarrier spacing 3840 kHz
Channel BW (data subcarriers + guard bands) 2160 MHz
Number of BSs / numbers of UEs 1/1
Code Rates 0.75
SC-FDE role off 0.22
Modulation format QPSK, 16-QAM, 64-QAM
CPE removal per DFT-s-OFDM symbol, or per SC-FDE
PN compensation
block
Channel AWGN
receiver type linear (MMSE CE + MMSE equalization)
FEC NR low-density parity-check (LDPC)
PA output power back-off dB 3 to 10 dB
The results of this comparison are shown in Figure 5-9. In this figure, only the block error rate (BLER)
for an output back-off (OBO) of 10 dB is shown. It is obvious that in this case the performance is only
limited by the phase noise. The comparison shows that there is a limited difference between SC-FDE
and DFTS-OFDM with the given system configuration. However, it needs to be noted that in the shown
simulation the subcarrier spacing for the DFTS-OFDM system is relatively large and the phase noise
given by the model is not too large to limit system performance. Therefore, it is concluded that there is
a benefit of using SC-FDE over DFTS-OFDM if PN is the limiting factor of the system. For the case
that phase noise is not limiting system performance either SC-FDE or DFTS-OFDM can be chosen.
Figure 5-9: Comparison of SC-FDE and DFT-s-OFDM with phase noise and non-linear PA distortion.
Parameter Value
Carrier frequency (GHz) 120
FFT size 4096
Occ. BW (GHz) O:3.17 / D:295 / S:3.92
Coding LDPC 3/4
RRC β 0.3 (SC-FDE only)
Cyclic Prefix (%) 7
PT-RS (%) 3
Channel AWGN
Nr.of Tx, Rx Antenna 1,1
Details on the applied PN model can be found in earlier report [HEX21-D22] section 3.6.6.2. PN
contribution on the received signal is estimated and compensated according to simple mitigation
algorithms for both waveforms (cf. [HEX21-D22] section 4.3.1). Uniform quantizer with equal spaced
steps is used. For DFTS-OFDM a typical 4 clipping is applied at the quantizer input where σ2 is the
variance of the zero-mean gaussian distribution of the received signal amplitude. For the block
transmission SC-FDE signal the amplitude distribution is rather formed by the shape of the applied
filter. An adaptive quantizer is therefore deployed here which dynamically adjusts the dynamic range
and step size to the input amplitude on the symbol basis.
Figure 5-10: Dependency on the deployed number of quantization bits of SC-FDE (upper row) and DFT-
S-OFDM (lower row) for QPSK, 16QAM and 64QAM.
Figure 5-10 shows the dependency of BLER performance on the number of deployed quantization bits
of SC-FDE and DFTS-OFDM for QPSK, 16QAM and 64QAM. As can be seen from the plots, 4-5
quantization bits are in general sufficient for transmissions with 64QAM. Further SC-FDE shows a
more robust behaviour in case of very coarse quantization. Compare, for example, the 64QAM case,
optimal BLER performance can be reached with 4 quantization bits for SC-FDE transmission whereas
5 bits are needed for DFTS-OFDM to achieve the same BLER performance. Applying the power
consumption model described in [Mur22] this means that the ADC power consumption using SC-FDE
transmission waveform can be of a factor of 2 lower compared to the case of using DFTS-OFDM
waveform.
5.1.2.2.1 PA friendly constellation pattern and pulse shaping filter for SC-FDE
The signal amplitude envelope determines the peak to average power ratio (PAPR) level of the transmit
signal. Thanks to the simple transmit structure of SC-FDE design the modulation alphabet can influence
the envelope variation of the signal and hence reduce the PAPR. 3GPP specifies standard square QAM
as modulation alphabets for NR, as exemplary shown for 16QAM in Figure 5-11 (left plot). It’s
primarily designed for the subcarrier of an OFDM signal, where the impact on PAPR from individual
QAM schemes on the subcarrier is almost negligible. The impact of the QAM shape is relevant only
for single carrier signals (and depends there also further on other parameters like filter shape). Another
standard utilizing single-carrier waveform is the digital video broadcasting second generation (DVB-
S2) for satellite broadcasting [DVB-S2X]. As modulation constellation amplitude phase shift keying
(APSK) patterns are defined. Looking at the left plot of Figure 5-11, it’s intuitive that those corner
constellation points on the outer ring contributes the most to the signal peaks. Rearranging them to other
inner rings helps to reduce the variation in the signal envelope. On the other hand, reducing the
Euclidian distances between the constellation points increases the decoding error probability at the
receiver and hence the system performance. It is a trade-off between the PAPR reduction and the
tolerable packing density of the constellation points.
Figure 5-11: 4-bits constellation pattern for standard NR 16QAM (left) and PA friendly 16APSK (right).
A selected 2 rings constellation pattern is used now as an example [DVB-S2X]. The ratio of the ring
radii, 𝑟1 /𝑟2 , is 0.6. The number of constellation points on the outer and inner rings is 12 and 4,
respectively. Figure 5-12 shows the complementary cumulative density function (CCDF) distribution
of the PAPR (left plot) and the BLER performance (right plot) of the chosen constellation pattern. As
baseline PAPR and BLER performance of NR standard 16QAM constellation are used. Table 5-9
summarizes the main simulation parameters used to produce the results shown in Figure 5-12 & Figure
5-13.
It can clearly be seen from the plots that a large PAPR reduction of 1.8 dB can be achieved using
16APSK pattern compared to NR 16QAM case. It is, however, at cost of 0.3 dB loss in BLER
performance. When considering PN with simple estimation and compensation method as described in
[HEX21-D22], section 4.3.1, an additional 0.2 dB loss must be accounted for. In a scenario where PA
output power is the limiting factor, the PAPR reduction can be well harvested for coverage extension.
The loss in BLER performance can be compensated by e.g., a more robust channel coding rate i.e., at
some spectral efficiency cost. Higher modulation orders tolerate smaller decoding error probability.
Therefore, less PAPR reduction from reshaping constellation can be expected.
For the 6-bits modulation case, an exemplary APSK constellation with 8, 16, 20, 20 constellation points
on 4 rings is chosen (inset of the left plot) [DVB-S2X]. A code rate of ¾ was used and the radian ratios
are 0.27,0.59,0.97,1.40, respectively. As in case of 4-bits modulation above, Figure 5-13 depicts the
CCDF distribution of PAPR (left) and BLER performance (right) for the 4-ring APSK modulation
pattern (solid blue) and the baseline NR 64 QAM (solid black). It can be observed for the selected
APSK constellation pattern, a PAPR reduction by 0.6dB can be achieved, without any further penalty
in BLER performance for both cases, with and without PN consideration.
4-bits Modulation
Figure 5-12: 4-bits modulation: CCDF of the PAPR (left) and BLER as a function of required Eb/No
(right) of SC-FDE waveform using NR 16QAM and 16APSK constellation pattern.
6-bits Modulation
Figure 5-13: 6-bits modulation: CCDF of the PAPR (left) and BLER as a function of required Eb/No
(right) of SC-FDE waveform using NR 16QAM and 16APSK constellation pattern. Inset: applied 4 rings
constellation pattern for 64APSK.
As already mentioned above, pulse shaping filter is an additional design parameter for SC-FDE
transmissions which can be used to constrain the signal envelope and hence is more PA friendly and
efficiently extends the achievable transmission range. Standard root raised cosine (RRC) filter is
commonly applied to band limit a transmission. Explicitly designing the pulse shaping filter to minimize
the peak amplitude of the transmit signal allows to further reduce the PAPR level. The applied RRC
filter shown in Figure 5-12 and Figure 5-13 above uses roll-off factor of 0.3. The designed PAPR aware
filter, with similar excess bandwidth, achieves an additional 0.4 dB PAPR reduction, without further
penalty in BLER performance (green curves in Figure 5-13).
The study in this subsection shows the potential of SC-FDE transmissions to reduce the PAPR level by
selecting suitable constellation patterns and pulse shaping filter with moderate to no losses in the BLER
performance. This PA friendly behaviour makes this waveform a potential candidate for transmissions
in scenarios in which PA output power is a limiting factor and power/energy efficiency is of cardinal
relevance.
frequency, which in turn requires lower operational power to avoid dielectric breakdown. Decreased
power capability of PAs will naturally compromise coverage/performance at sub-THz, already
compromised by increased pathloss. One way of improving coverage consists of pushing the operating
point of PAs closer to saturation power level without increasing the level of nonlinear distortion, which
may be achieved by choosing waveforms with small envelope variations. One such waveform is DFT-
spread OFDM (DFTS-OFDM), which is akin to ordinary OFDM but with a smaller envelope variation.
In this contribution, ordinary DFTS-OFDM and a modification of DFTS-OFDM (employing frequency-
domain filtering) are compared in terms of potential coverage extension, under the influence of a
nonlinear PA.
DFTS-OFDM and FDF-DFTS-OFDM
0
.
.
.
0
Add Wideband
N-Point P/S
CP/GI/UW Filter
IFFT
QAM Nc-Point
symbols FFT
Nc assigned carriers
f
Total N subcarriers
from FFT
This bandwidth expansion is illustrated by the grey-coloured “tails” of the spectrum at the IFFT input
in Figure 5-15. On the other hand, if the occupied bandwidth is constrained to be exactly 𝐵𝑊𝑑𝑎𝑡𝑎 (e.g.
in order to satisfy coexistence requirements), then the number of data-bearing subcarriers has to be
reduced so that they occupy a reduced bandwidth
′
𝐵𝑊𝑑𝑎𝑡𝑎 (5-2)
𝐵𝑊𝑑𝑎𝑡𝑎 = .
(1 + 𝛽)
This case is denoted as FDF-DFTS-OFDM without bandwidth expansion. In order to ameliorate the
reduction of transmit power caused by reduced bandwidth, the pre-shaping signal power should be
boosted by a factor of 1 + 𝛽 (this will be referred to as power boosting in the follow-up). Moreover, the
loss of raw throughput induced by bandwidth reduction (2) can be compensated by increasing the code
rate by a factor of 1 + 𝛽. The principles of power boosting and of bandwidth expansion/no bandwidth
expansion are illustrated in Figure 5-16, where (1 − 𝛾) = 1/(1 + 𝛽).
FD
Repetition (Nyquist) Add Wideband
N-Point P/S
with Filter with CP/GI/UW Filter
IFFT
Period Nc roll-off
QAM Nc-Point factor b
symbols FFT
Figure 5-18: link budget comparison of two systems using different waveforms.
A companion illustration for the proposed system comparison methodology is given in Figure 5-18.
The link budget on the left-hand side corresponds to a hypothetical system using a waveform with high
envelope variation, requiring a large PA OBO to ensure low probability of clipping and consequently
low distortion. The link budget on the right-hand side may correspond to a system using a waveform
with low envelope variation but requiring a higher demodulation SNR to support the same performance
as system 1. There is a difference in the link budget headroom (labelled as maximum coupling loss
(MCL) in the figure) between the two systems, with system 2 having a larger MCL and therefore being
able to accommodate a larger propagation loss.
DFTS-OFDM vs FDF-DFTS-OFDM: comparison of minimum power backoff
Minimum power backoff of the PA is the difference (in dB) between the saturation power of the PA
and the mean output power, chosen to ensure that in-band and out-of-band performance requirements
are satisfied. In-band requirements are quantified by error vector magnitude (EVM) requirements, given
in Table 5-10.
OOB requirements are quantified by adjacent channel leakage ratio (ACLR), calculated as the ratio of
signal power in configured bandwidth 𝐵𝑊𝑐𝑜𝑛𝑓𝑖𝑔 centered at the carrier frequency 𝑓𝑐 (in-band power),
and signal power in configured bandwidth 𝐵𝑊𝑐𝑜𝑛𝑓𝑖𝑔 centered at frequency 𝑓𝑐 + 𝐵𝑊𝑐ℎ𝑎𝑛𝑛𝑒𝑙 , where
𝐵𝑊𝑐ℎ𝑎𝑛𝑛𝑒𝑙 ≥ 𝐵𝑊𝑐𝑜𝑛𝑓𝑖𝑔 . ACLR requirement for the evaluations in this contribution is set to 26 dB. The
PA model used in the evaluations is the memoryless GaN model described in [HEX21-D22, Section
3.6.6.1]. As for other parameters of the evaluation, subcarrier spacing of both waveforms is set to 960
kHz, 𝐵𝑊𝑐ℎ𝑎𝑛𝑛𝑒𝑙 = 2.16 GHz and 𝐵𝑊𝑐𝑜𝑛𝑓𝑖𝑔 = 1.88 GHz. Results are shown in Figure 5-19, and
indicate that use of frequency-domain shaping reduces the OBO by 0.5 – 1 dB, depending on the
modulation used.
DFTS-OFDM vs FDF-DFTS-OFDM: link level performance
In the second step in the analysis, link-level performance of the two waveforms under investigation is
assessed. The set of modulation and coding schemes (MCSs) selected for the analysis is taken from
[38.214, Table 6.1.4.1 – 1, q = 1]. The MCSs for FDF-DFTS-OFDM are the same as for DFTS-OFDM
but with code rate adjusted:
𝑟𝑐,𝐹𝐷𝐹−𝐷𝐹𝑇𝑆−𝑂𝐹𝐷𝑀 = 𝑟𝑐,𝐷𝐹𝑇𝑆−𝑂𝐹𝐷𝑀 (1 + 𝛽). (5-3)
A selection of link level performance results is shown in Figure 5-20, illustrating the differences in
BLER performance between DFTS-OFDM and FDF-DFTS-OFDM. MCSs selected for FDF-DFTS-
OFDM match those for DFTS-OFDM in throughput, as they effectively use less BW. For the
evaluations, practical channel estimation was used. Channel model is SISO CDL/TDL-A with delay
spread of 10 ns and UE speed of 3 km/h. Block fading was assumed in the evaluations, i.e. a random
channel impulse response is generated for each block, with block equal to one slot (14 DFTS-OFDM
symbols).
The results clearly show the negative impact on performance that the code rate increase at higher
modulation orders has on the performance of FDF-DFTS-OFDM, the effect which will eventually offset
the OBO gains when considered in the context of full link budget evaluation.
Figure 5-20: BLER comparison of DFTS-OFDM and FDF-DFTS-OFDM under the influence of PA
nonlinearity.
DFTS-OFDM vs FDF-DFTS-OFDM: total link budget comparison
Motivated by the illustration in Figure 5-18, let us define the maximum coupling loss (MCL) as
𝑀𝐶𝐿 = 𝑃𝑃𝐴,𝑚𝑎𝑥 − 𝑂𝐵𝑂 + 𝐺𝑎𝑛𝑡 − (𝑆𝑁𝑅 + 𝑃𝑛 ) [𝑑𝐵], (5-4)
where 𝑃𝑛 is thermal noise power and 𝐺𝑎𝑛𝑡 total antenna and array gain at Tx and Rx side. Let us also
define the information throughput as
𝑏𝑝𝑠 (5-5)
𝑇ℎ = 𝑆𝐸(1 − 𝐵𝐿𝐸𝑅) [ ],
𝐻𝑧
where SE denotes spectral efficiency. Now, between two transmission systems with same throughput
but different MCL supporting those throughputs, the system having a larger MCL achieves larger
coverage. Conversely, when comparing two systems operating at the same MCL, the system achieving
a larger throughput is “better” since it achieves a better throughput at the same coverage as the other
system.
What the results reveal is that FDF-DFTS-OFDM has advantage over DFTS-OFDM only at the lowest
spectral efficiencies, and that it provides at most 0.8 dB increase in MCL (corresponding to a coverage
extension of about 10% if free-space pathloss is assumed). Therefore, FDF-DFTS-OFDM may be
considered for supporting a moderate coverage extension for the users at the cell edge, whereas for
other users DFTS-OFDM is preferred.
5.1.2.3 Optimal pilot configuration for phase noise compensation at 140 GHz
This section presents a systematic numerical analysis on the pilot design for phase noise compensation
of DFTS-OFDM waveform at the 140 GHz band. The pilot scheme used in this analysis follow the
approach of the phase tracking reference signals (PTRS) of 3GPP [38.211], which consists of 2, 4 or 8
pilot groups with samples each. In addition, in order to have a more general set of pilots which is more
adequate to the sub-THz channel, an extension is considered with 16, 32 and 64 groups of pilots to
analyse the performance and data rate trade-off. In addition, each pilot group occupies 4 time-domain
samples to have a comparable setup with the PTRS of 3GPP. For compactness, only the indexes of the
first pilot of each group is displayed in Table 5-11 for 𝑁g =16, 32 and 64 groups, which were defined
by distributing the pilots as equally as possible along the time domain samples. For example, the pilots
for 𝑁g = 16 are [1 2 3 4 137 138 139 140 …], where the indexes [1 137 …] represent the first index of
each pilot group.
terms of energy per bit with the minimum performance loss of 0.9, 0.9 and 0.8 dB for the SCS of 960,
1920 and 3840 kHz, respectively. Increasing the number of pilot groups to 𝑁g = 64 increases the
performance loss to 1 dB because the amount of energy spent in the additional pilots does not
compensate the small performance gain in terms of SNR, which does not make it an interesting choice.
In addition, it is important to highlight that the pilot group number of 𝑁g = 16 is also an option to be
used, especially with SCS of 1920 and 3840 kHz, since the performance loss in relation to 𝑁g = 32 is
small. For SCS of 3840 kHz, the configuration 𝑁g = 8 is an option for the same reason.
Figure 5-22: Block Error Rate (BLER) for 16-QAM, R = ¾, subcarrier spacing of 960, 1920 and 3840
kHz, with Ng=2,4,8,16,32 and 64 pilot groups.
Table 5-12: Loss in terms of Eb/N0 for 16QAM pilot schemes at BLER = 10-2. NA = not achieved.
SCS 𝑁g = 2 𝑁g = 4 𝑁g = 8 𝑁g = 16 𝑁g = 32 𝑁g = 64
960 kHz NA NA 2.9 dB 1.2 dB 0.9 dB 1dB
1920 kHz NA 5.1 dB 1.6 dB 1 dB 0.8 dB 1dB
3840 kHz NA 2.7 dB 1.2 dB 0.85 dB 0.8 dB 1dB
The simulation has also been carried out for 64-QAM with 3/4 code rate NR LDPC code. The
performance loss results are shown in Table 5-13. As expected, for the higher MCS, the phase noise
has a higher impact in the performance loss. In addition, one important aspect to notice is that for SCS
of 960 kHz, all the PTRS of 3GPP [38.211] with of 𝑁g = 2, 4 and 8 pilot groups do not achieve the
BLER of 10-2, meaning that the extended set of pilots is necessary in this case. Also, the system with
𝑁g = 64 pilot groups provide the best trade-off for SCS = 960 and 1920 kHz, however, it is important
to notice that the performance gain of 𝑁g = 64 in relation to 𝑁g = 32 is relatively small, e.g., 0.4 and
0.1 for the SCS of 960 and 1920 respectively. For SCS = 3840 kHz, there is almost no gain of 𝑁g = 64
in relation to 𝑁g = 32. Thus, in practise 𝑁g = 32 can be a suitable choice because it leads to less pilot
overhead. Lastly, for the SCS of 3840 kHz, 𝑁g = 16 is also an option due to small loss with higher 𝑁g .
Table 5-13: Loss in terms of Eb/N0 for 64QAM pilot schemes at BLER = 10 -2. NA = not achieved.
SCS 𝑁g = 4 𝑁g = 8 𝑁g = 16 𝑁g = 32 𝑁g = 64
960 kHz NA NA 3.2 dB 1.5 dB 1.1dB
1920 kHz NA NA 1.7 dB 1 dB 0.9 dB
3840 kHz NA dB 2.3 dB 1.2 dB 0.9 dB 0.9 dB
Although these simulations considered 16-QAM and 64-QAM for FFT size of 2048, comments about
possible outcomes for other configurations can be discussed. In particular, it can be expected that the
FFT size does not play an important role in this case because it does not change the time duration of the
DFT-s symbol, which is solely dependent on the SCS choice. Also, as the MCS gets larger, e.g., 128/256
QAM, the required 𝑁g increases as observed from 16-QAM to 64-QAM. In addition, another case is
considered using a highly directive antenna where the multipath components (MPCs) are too small such
that the channel collapses into an AWGN. In general, if there are MPCs, it is expected that an extra
performance loss due to imperfect channel estimation and additional noise in the equalization process
at the receiver, i.e., the necessary values of 𝑁g found in Table 5-12 and Table 5-13 are expected to
increase, which increases the need of using the pilot schemes defined in Table 5-11.
For this investigation, the phase noise model from [KKP+14] is adapted, which models the phase noise
as a non-stationary process and that is commonly used to describe the phase noise behaviour of free
running oscillators. Here, the phase noise stems from two independent processes:
𝝓𝒌 = 𝝓𝟎,𝒌 + 𝝓𝟐,𝒌 ,
where 𝜙0,𝑘 is a white Gaussian phase noise process with variance 𝜎𝜙20 and 𝜙2,𝑘 is cumulative phase
noise processes, that stems from the integration of white noise inside the oscillators and can be
expressed in terms of an initial phase and a sum of phase increments:
𝜙2,𝑘 = 𝜙2,1 + ∑𝑘𝑖=1 𝜉2,𝑖
It can be shown that the phase increments are Gaussian distributed with variance 𝜎𝜉2 . Moreover, the
increments 𝜉2,𝑖 have no temporal correlation, i.e., 𝜙2,𝑘 is a Wiener process. In [KKP+14], there is an
additional third phase noise component, 𝜙3,𝑘 , that is left out here, as its energy is mostly concentrated
at very low offset frequencies, and which therefore can almost be treated as a constant phase offset.
Figure 5-23: Model of the time frame structure, taken from [GSD+22].
where 𝑥𝑘 and 𝑦𝑘 denote the samples from the noise free receive signal and the receive 1-bit quantized
samples, respectively. The index of the current pilot block within the observation frame is denoted by
𝑚. This estimate can be interpreted as an estimate of the average phase within the pilot block, i.e.,
(𝑚+1)𝑁
1 (5-7)
𝜙̂𝐿𝑆 (𝑚) = 𝜙̅(𝑚) = ∑ (𝜙0,𝑘 + 𝜙2,𝑘 ).
𝐿𝑝
𝑘=𝑚𝑁+1
In the post-processing stage, using knowledge on the statistics of the phase noise process, the estimates
of the pilot blocks are interpolated to acquire phase estimates for the data blocks. Here, two different
algorithms are used: a Kalman filter and the Rauch-Tung-Striebel (RTS) algorithm. The Kalman filter
works in two steps: a model-based prediction step is followed by a data-based correction step. It is
capable of using all observations from the past and is therefore a causal filter. The RTS algorithm forms
a non-causal extension to the Kalman filter and therefore requires a closed observation frame. It
performs a forward and a backward filtering sweep that both use the idea of the Kalman filter, but in
the backward sweep on the time-reversed sequence. Note that the latter approach leads to a lower MSE
at the cost of additional latency.
To set the results of the estimation algorithms into perspective, the optimal estimation performance that
can be achieved based on 1-bit quantized observations was also examined in [ZGD+22]. For this
purpose, the Bayesian Cramér-Rao bound (BCRB) was formulated. The BCRB is a lower bound for the
achievable Bayesian mean-squared error (MSE) of any estimator and is given by the inverse of the so-
called Bayesian information matrix that can be formulated as a sum of the expected Fisher information
and a term that reflects the a-priori knowledge that estimators can employ:
𝑩 = E𝝓 [𝑭(𝝓)] + E𝝓 [−∇2𝝓 log 𝑝𝝓 (𝝓)]. (5-8)
In the equation above, 𝝓 denotes the vector of estimation parameters, which in this case corresponds
the phase noise samples. Calculating the matrix for the system model presented here is not
straightforward and the details of the derivation can be found in [ZGD+22]. It is important to note that
the bound has been calculated only holds for estimators that create estimates solely based on
observations made during pilot observations and discard available information about the estimation
parameters that is contained in the observations that are made during data transmission. However, this
covers most practical estimation algorithms, as it allows the separation of synchronization and data
detection and, thus, enables low complexity receiver algorithms. The calculated bound is given by
−1
𝑐1 𝐸𝑠 −𝑐2 𝐸2𝑠 𝑐2 𝐸𝑠 𝑐2 𝐸𝑠 (5-9)
𝑩 −1
= ( 2 𝑒 𝜎𝑛 (𝐼0 ( 2 ) + 𝐼1 ( 2 )) 𝚫𝑃 + 𝑪−1
𝝓 ) ,
𝜋𝜎𝑛 𝜎𝑛 𝜎𝑛
where 𝜎𝑛2 is the variance of the additive noise in the channel, 𝐸𝑠 denotes the symbol energy (which for
this formulation is assumed to be constant for all symbols), 𝑐1 ≈ 4.0360 and 𝑐2 ≈ 0.3930 are
numerically determined approximation constants, and 𝐼𝑛 (𝑥) denotes the modified Bessel function of
the first kind and n-th order. The pilot-data structure is captured in the matrix 𝚫𝑃 that is a diagonal
matrix with ones on the main diagonal where pilot blocks are located within the observation frame and
zeros where data blocks are located. The matrix 𝑪𝝓 denotes the covariance matrix of the estimation
parameter vector 𝝓. As the RTS algorithm, the BCRB uses all available information that can be
extracted from the temporal correlation of the phase noise process and therefore, also requires a finite
observation frame.
In Figure 5-24 the Bayesian MSE of the estimation algorithms is shown, along with BCRB, i.e., the
diagonal elements of 𝑩−1 , over the time index 𝑘. The parameters used for this evaluation can be found
in Table 5-14. Naturally, the bound and the algorithms show the lowest MSE at the positions of the
pilot blocks, which explains the periodic minima in the MSE curves. These minima are very pronounced
for the Kalman filter but can be barely seen for the RTS algorithm and the BCRB, indicating that the
RTS algorithm could be used with even larger pilot block spacing. Moreover, the RTS algorithm
performs very close to the BCRB, and the MSE of the RTS algorithm and the BCRB exhibit a “bathtub”
shape. This is intuitive, as the missing observations before and after the observation frame explain the
rise of the MSE curve at the ends.
Figure 5-24: MSE of the estimators and BCRB over the sample index k within the observation frame,
taken from [ZGD+22].
is distributed into two components and their output phase is then varied to achieve the varying envelope.
Typically, the combining is performed using different types of conducted combiners in the circuit level
to feed the varying envelope signal to the antenna. Such combining shown in Figure 5-25. In principle,
any signal can be divided into two out-phase components that have constant amplitude. Let us simply
assume 𝑥(𝑡) to be the input signal of the system that has varying envelope. Now, 𝑥(𝑡) can be written
as a sum of two components 𝑥1 = exp(1𝑗𝜑1 (𝑡)) and 𝑥2 = exp(1𝑗𝜑2 (𝑡)) such as 𝑥(𝑡) = 𝑥𝑝ℎ𝑠 (𝑥1 (𝑡) +
𝑥2 (𝑡)). Note that |𝑥1 (𝑡)| = |𝑥2 (𝑡)| = |𝑥𝑝ℎ𝑠 (𝑡)| = 1. In out-phasing transmitter, the phases of the
signal components are selected to be 𝜑1 (𝑡)) = −𝜑2 (𝑡)). Now, by varying the phase, any value in the
real axis can be produced. This makes it possible to generate varying amplitude by two constant
envelope signals. The output phase of the signal is achieved by rotating the overall signal then also with
a phasor 𝑥𝑝ℎ𝑠 (𝑡) = exp(𝑗𝜑𝑝ℎ𝑠 ) . Hence, the overall signal consisting of amplitude and phase
components 𝑥(𝑡) = |𝑥(𝑡)|exp(𝑗 arg(𝑥(𝑡)) can be represented such that |𝑥(𝑡)| = |𝑥1 (𝑡) + 𝑥2 (𝑡)| and
arg(𝑥) = 𝜑𝑝ℎ𝑠 . This basically enables to have only phase-varying signal feeding each PA, while the
amplitude is generated as a sum of the components at the output. This enables to have low PAPR per
PA. In general, any type of waveform can be reproduced by using this technique.
(b)
(a)
Figure 5-25: Illustration of typical out-phasing concept.
Using the concept above in very high centre frequencies such as upper mmWave region is challenging.
In the circuit side, one of the challenges is the output-combiner that is typically rather lossy and hence
reduces the overall efficiency and reduces the maximum achievable output power. One way to solve
the problem is to replace the output combiner by simply individual antennas as proposed in [LR11].
This generally makes the output signal by combining two signals over-the-air. The interesting thing in
such concepts is that it eventually makes the varying amplitude by (digital) beamforming and nicely
combining waveform and modulation design, beamforming and RF design together. In large subarray-
based hybrid beamforming systems this concept can be utilized by combining multiple subarrays in
similar way. Such a concept that is shown in Figure 5-26.
be done in the digital side, while the over-the-air combining is then over two subarrays. In principle,
this is basically generating the overall amplitude-varying signal by hybrid beamforming over two
subarrays. Similar concepts are presented in the literature, for example, in [LXA+21] and [ASS+20],
where in the latter one the concept is called as “outspacing” transmitter.
The basic requirement for applying the concept in practice is to have a multi-antenna transmitter with
at least two digital signal chains. Applying the concept in 6G transmitter is possible as the architecture
has multiple subarrays. Hence, it can be implemented as an alternative method to generate the higher
PAPR signals. The main thing is that it does not require any extra RF hardware components. Added
complexity comes from the signal separation to amplitude and phase components and amplitude
separation to two out-phased components. This could be addressed in the waveform design side to
generate the waveform to directly provide these two components. For example, in typical single-carrier
modulations such as QAM this is straight forward [LR11]]) and performed in such a way that the typical
pulse shaping is then performed after the signal separation for both out-phasing components. However,
with typical multi-carrier transmitters, this is more challenging to address directly inside the waveform
generation. For generic 6G transmitter architecture, out-phasing by hybrid beamforming can be enabled
as a separate mode to use higher spectral efficiency modulation in some cases, while in other cases, the
individual subarrays can operate in more typical modes serving users at dedicated directions with their
own waveform or serving a single direction with one waveform same fed to both subarrays.
The main challenges of the technique are that it requires accurate calibration of the subarrays to make
the arrays exactly phase coherent. This is expected to be challenging with very wideband waveforms.
In addition, the other challenge is related to the fact that the quality (linearity) of the combined signal
is highly direction dependent that is also analysed in the original reference [LR11] by showing that the
achieved bit error rate is very sensitive for steering angle. The impact is similar than in directive
linearization schemes presented in [Ter22, TAT+17]. This further increases the calibration requirements
also to beam domain. On the other hand, direction dependent behaviour of the waveform can be also
beneficial as discussed in [Ter22] to, for example, make it even more challenging to detect the signal
outside the steering angle.
Note that the used technique does not, as such, increase the achieved 𝑃𝑠𝑎𝑡 in the steering angle as the
maximum amplitude is achieved when the out-phasing is selected to be zero. The overall amplitude
variations are, hence, achieved by using the array gain and the maximum array gain is achieved for the
peaks of the signal. Note that same maximum 𝑃𝑠𝑎𝑡 for the signal peaks would limit any type of
waveform eventually, but it is likely that in out-phasing mode it could be possible to go rather close to
𝑃𝑠𝑎𝑡 as the signal is per PA is expected to have low PAPR. Also, it enables to operate the PAs in more
efficient mode that, for example, could help the thermal management of the amplifiers. Moreover, the
PAs can be biased in more nonlinear mode for this case that can further improve the efficiency and
hence reduce the power consumption of the power amplifiers.
4.2 of this deliverable, also phase noise impact are dependent on the signal bandwidth in higher centre
frequencies. Hence, bandwidth also play a crucial role in the achievable link ranges.
Carrier aggregation can be also effectively used to scale bandwidth for different link scenarios. The
link-adaptation can be done basically in two domains: vary the modulation and coding scheme (MCS)
or vary the bandwidth. Varying MCS is simply changing the spectral efficiency. However, highest
spectral efficiency MCS is limited eventually at least by the EVM of the transmitter. The goal
envisioned in [HEX21-D21] was to reach 100 Gbit/s with 100-200 m link distance. This can mean 15
– 40 GHz of signal bandwidth to reach with realistic MCS selection in line-of-sight conditions. To get
the numbers into a context, the following simple analysis is conducted. Let us assume that the data rate
with a certain SNR and bandwidth can be calculated as
𝐷𝑅 = 𝐵[log 2 (1 + 𝛾(𝐵, 𝑑, 𝑓𝑐 )) − 𝑆𝐸𝑙𝑜𝑠𝑠 ], (5-10)
where γ is the effective SNR (caused by all RF nonidealities) in linear scale that depends on the signal
bandwidth B, centre frequency 𝑓𝑐 , and link distance 𝑑. The distance and frequency dependency are
caused by path loss, while the bandwidth dependency is caused by the noise bandwidth of the receiver.
The term 𝑆𝐸𝑙𝑜𝑠𝑠 represent the loss in the spectral efficiency with respect to the theoretical Shannon
capacity, which according to [Ter22: figure 2] is around 1.5-2bit/s/Hz for typical QAM modulations.
By following the methodology presented in [TTP17, Ter22, LTJ19, LJT19], the contributions of the Tx
EVM and Rx noise are combined to evaluate the effective signal to noise ratio achieved at different link
distances, bandwidth, and centre frequency.
A simulation example is conducted to highlight how the bandwidth vs distance relation can be seen in
different link scenarios. The simulation parameters are shown in Table 5-15. The TX power is taken as
a single PA power following the Silicon germanium (SiGe) process saturated power estimated derived
in Section 4.3 with roughly 9 dB output backoff from 𝑃𝑠𝑎𝑡 . Similarly, low-noise amplifier (LNA) noise
figure is following the estimates of the SiGe process derived in [HEX21-D21]. Used antenna element
gains are typical patch antenna gains and two different antenna configurations are selected. Note that in
practice the number of antennas is limited by power consumption, area as well as Tx power limitations
set by different authorities. The analysis is conducted for 150 GHz and 300 GHz centre frequencies
with free space path loss assumption. The results for four different scenarios are shown in Figure 5-23.
Each curve in the figure shows the fixed achieved data-rate with a certain link distance and bandwidth.
The data rate is derived based on Equation (5-10). Let us first look the scenario at 150 GHz at Figure
5-23 (a) – (b). The clear observation when comparing cases 1 and 2 at 150 GHz is that if the power
budget, mainly realized by higher number of antennas, satisfy well the SNR constraint of highest
achievable spectral efficiency (limited by Tx EVM), there is no significant link-distance dependency
on the data-rate in decent link distances up to 500 m. This is seen in case 2, where the curves for each
target data-rate are flat lines. However, in case 1, the link distance dependency in large bandwidths
becomes more evident. The scenario can be more clearly observed at cases 3-4 in 300 GHz shown in
Figure 5-23 (c) – (d). When looking this as a function of bandwidth, there seems to be “best bandwidth
selection” to achieve a certain data rate with a certain link distance. This dependency is mainly caused
by the fact that the different link scenarios in cases 3-4 clearly are in the boundary of Tx EVM limited
and Rx noise limited region. Hence, in these scenarios, varying both MCS and bandwidth, to serve a
user at certain link distance would be possible. To make it clear, highest MCS should be generally
always chosen if the SNR budget can afford it, but as the SNR depends also on the bandwidth, they
should be varied both for best compromise in the link adaptation. In practice, the best way the allocate
the bandwidth would be to vary the number of active CCs.
Table 5-15: Simulation parameters for data-rate analysis with different link distances and signal
bandwidth.
Case 𝒇𝒄 Tx EVM Ptx Ntx Gatx NF Garx Nrx
1 150 5 10 256 5 6 4 16
2 150 5 10 1024 5 6 4 64
3 300 5 -4.5 256 5 17 4 64
4 300 5 -4.5 1024 5 17 4 256
(a) (b)
(c) (d)
Figure 5-27: Simulated data rates at different bandwidth and link distance in four scenarios (a)-(d)
derived in Table 14.
The tasks performed in the state initial beam-alignment can consist of a number of procedures.
Combined they have the target that after the process is finalized communication via the interface
utilizing the beamforming is possible. In a standalone 3GPP system this process would consist of first
the initial beam-sweeping using the synchronization blocks with the UE measuring and selecting the
optimal beam. The mobile device would afterwards use the knowledge gained from receiving the
broadcast information to perform the random-access procedure. Afterwards a beam that can be utilized
to establish a communication between these points is available, thus the mobile device and the access
point have their beams connected. Note that it is also possible that the beam alignment fails, and the
system would transition back to the BM idle state.
connection no
longer needed
quality
criteria
initial beam- initial beam- satisfied
alignment alignment
BM parameters initial beam- success beams
idle alignment connected
initial beam-
alignment failed
beam recovery
success quality criteria
not satisfied
beam recovery
failed beam
recovery
beam pairs. Using the synchronization sequence with the broadcast channel as in NR for this beam
training this would result in 4.56 s of just transmitting reference signals to train the beams, even if the
subcarrier spacing is scaled o 3.84 MHz (28 * 15 kHz).
It is obvious that it will not be feasible to spend such a large amount of transmit power and time to align
the beams of a single mobile device. Note that how severe this problem is going to be depends on the
number of beams utilized by the access point as well as the mobile device. The number of beams depend
on the coverage requirement as well as the channel conditions. At this point in time these plus the
hardware capabilities of the future are not clear, thus accurate prediction and design of the related
procedure need to wait till more information is available. The same observations made for the initial
beam-alignment can be made for beam recovery. In this case due to the additional information of a past
beam that was fulfilling the quality criteria is available in addition to the information available for initial
beam alignment.
Different techniques to reduce the number of beams to be tested for initial beam-alignment as wells as
beam recovery need to be defined. All these can be categorized into utilizing any time of side
information available to the system. This side information can be historic observations, optical
observations data, measurements at another frequency, positioning information, and information about
the radio environment. Any combination of these is possible, however what is available for future
systems operating in the frequency range from 100 to 300 GHz is at this point unclear. But as the system
would be very inefficient with utilizing these, treating them as an integral part of the system design is
paramount. For the quality monitoring of the connection with aligned beams a reactive or proactive
approach is possible. In the reactive approach different quality criteria of the connections utilizing these
beams are constantly monitored. In the case that the quality falls below a (potentially dynamic)
threshold the network or the mobile device declares that communication with the current beams does
fail. In the proactive approach in addition to all aspects of the reactive approach the system tries to
anticipate a future degradation of the communication quality utilizing the current beam pair. This
prediction can be based on any side information available to the network.
Combining this discussion about beam management it is clear that future communication systems
utilizing the frequency range from 100 to 300 GHz will only work efficiently if the exploitation of side
information is included at an early stage of the system design. One example of such side information
are localization and sensing information as described in [HEX23-D33].
This section presents the following works: (a) Centralized and distributed precoding, i.e., zero-forcing
(ZF) in particular, are studied with L1 mobility, i.e., serving RU subset update, see Figure 5-29.
Distributed ZF done in clusters/subsets are interference-aware, i.e., consider minimizing the power
interfering to other clusters [HEX21-D22]. (b) NCJT using SRT among multiple RUs is explored in a
D-MIMO setup. (c) The impact of serving subset and precoder update periodicities are examined for
both CJT and NCJT. The precoder must be updated when serving RU subset has been updated. The RU
subsets, on the other hand, may not require as frequent updates as the precoder because the key factor
for RU-UE association is slow fading which is mainly affected by the distance to the RUs.
where 𝒒 is the data symbol vector, 𝐇 ̂ denotes the estimated channel matrix and 𝐇̃ is the bias matrix
incorporating the channel aging impact, and 𝐖 is calculated by using estimated channel matrix such
that 𝐇
̂ 𝐖=I where I is the identity matrix; however, the received signal quality is still affected by
channel aging. In case of centralized ZF (CZF), precoder for all UEs can be formulated as 𝐖 =
̂ 𝑯 (𝐇
𝐇 ̂𝐇̂ 𝑯 )−1. In case of interference-aware distributed ZF (IADZF), UE-centric precoders per each
subset can be found as 𝐖𝒌 = 𝐇 ̂ 𝒌 (𝐇 ̂ 𝒌𝐻 + 𝚲𝒌 )−1 , where 𝐇𝑘 = [𝒉′𝑘1 𝒉′𝑘2 ⋯ 𝒉′𝑘𝐿𝑘 ] is the channel
̂ 𝒌𝐇
gain matrix for the 𝑘-th UE with 𝐿𝑘 serving RUs for the k-th UE, and 𝚲𝒌 is the interference covariance
matrix [HEX21-D22]. The received signal is perturbed by a term of √𝑃𝐇 ̃ 𝐖𝒒 due to channel aging that
is scaled by precoder, will eventually degrade the performance.
Performance Evaluation:
Performance of the D-MIMO system with L1 mobility and NCJT operating at 28 GHz is analysed in a
DL indoor scenario, for different number of regularly deployed RUs, randomly distributed 𝐾=16 UEs
and 1000 blockers. It is assumed that each RU and UE are equipped with single antenna, but it is
straightforward to extend to the multi-antenna case. During L1 mobility, different serving subset update
periodicities, i.e., 100 ms and 500 ms, are considered to update the serving RU subset. Other simulation
parameters are given in Table 5-16. Perfect channel estimation and lossless, high-capacity FH is
assumed. RF imperfections, hardware impairments, phase noise, power amplifier nonlinearities are
omitted from the scope of this work.
Figure 5-30 shows spectral efficiency (SE) performance with different number of serving RUs with
100ms serving subset update periodicity. If more candidate RUs are deployed, better interference
management can be done with increased degrees of freedom. Larger serving subsets can attain high
performance for CZF in general and IADZF in case of dense deployments. Nonetheless, IADZF can
provide less performance for a smaller number of deployed RUs since more interference starts to exist
in the UE specific clusters. As serving subset enlarges, i.e., number of serving RUs increases,
performance difference between CZF and IADZF scales up.
Table 5-16: Simulation Model Specifications
Parameters Model Specification
Bandwidth (MHz) 200
RU max transmission power (dBm) 13
Propagation Model 3GPP InH
Duplexing TDD with 50% DL
Channel coherence time (ms) 100
UE speed (m/s) 5
Figure 5-30: SE performance of CZF and IADZF Figure 5-31: SE performance of CZF with CJT
with CJT for different number of serving RUs in for different serving RU subset periodicities, e.g.,
a scenario with 16 RUs. 100 ms, 500 ms and 1 s.
Figure 5-32: SE performance for CJT and NCJT Figure 5-33: SE behaviour evaluated in 10s-time
schemes for different number of serving RUs in frame for CJT and NCJT schemes with serving
the scenario of 16 UEs when serving RU subset subset size of 4 RUs and different subset update
update periodicity is 100 ms. periodicities.
Figure 5-31 demonstrates the performance degradation if longer periodicity is adopted for serving
subset update, where more performance loss has been observed for longer periodicities. And the
degradation is higher in case of larger serving subsets. Frequent update is important to keep high
performance for centralized precoding such that it loses much performance with slower updates. Figure
5-32 shows the impact of NCJT which attains less performance in multiple serving RU subsets relative
to the CJT due to the not-optimized parameters. Provided that distributed precoding already performs
worse, NCJT will cause relatively less performance loss compared with CZF. Figure 5-33 presents the
performance w.r.t. time during mobility in case of serving subsets of 4 RUs in a region of 81 deployed
RUs for both serving subset update periodicities. Herein CJT provides good performance only when
precoder updates, whereas SE fluctuates with serious degradation under longer update periodicity that
is like NCJT performance. For NCJT, regardless of update periodicity, performance is substantially
reduced. In short, frequent serving subset update is necessary for CJT, but not critical for NCJT.
Figure 5-34: Modelling the topology of a dynamic D-MIMO network as a bipartite graph with partial
connectivity.
To make the topology structure of dynamic D-MIMO more clearly, Figure 5-34 is an example of
dynamic D-MIMO process, employing sparse bipartite graphs 𝒢 = {𝒜, 𝒰, ℰ}, where 𝒜 =
{𝑎1 , 𝑎2 , 𝑎3 , . . . , 𝑎𝑀 }, 𝒰 = {𝑢1 , 𝑢2 , 𝑢3 , . . . , 𝑢𝐾 }, denote the RU nodes and users nodes, respectively. In
Figure 5-34 (a) and Figure 5-34 (b), the blue circle nodes represent RUs while the orange triangle nodes
represent users. The upper two figures indicate the bipartite graph of the t-th snapshot, while the bottom
one indicates the (𝑡 + 1)-th after structural modification. In comparison to the previous snapshot, in
these four figures the red and black lines represent new and lost connections, respectively. Here, the
edge set ℰ is defined as ℰ = {𝑒(𝑎𝑚 , 𝑢𝑘 ), 𝑚 ∈ [𝑀], 𝑘 ∈ [𝐾]}.
In this way, the received signal at m-th RU t-th snapshot during training phase can be written as
𝑡 𝑡
𝒓𝑡𝑚 = √𝜌𝑇𝑝 ∑ 𝑔𝑚𝑘 𝒔𝑘 + √𝜌𝑇𝑝 ∑ 𝑔𝑚𝑘 ′ 𝒔𝑘 ′ + 𝒘𝑚 , (5-12)
𝑘∈𝒦 𝑚 𝑘 ′ ∉𝒦 𝑚
where 𝜌 is the transmit power of each user, and 𝑤𝑚 is the i.i.d. additive AWGN at RU-m with 𝒩(0, 𝐈),
𝑡
𝑔𝑚𝑘 denotes the channel coefficient between RU-m and user-k at t-th snapshot, 𝒦 𝑚 denotes the
dominated users set connecting with the RU-m, and 𝒔𝑘 is the pilot sequence of user-k.
In the aforementioned definition of the topology structure of D-MIMO, the topology graph
modifications at the (𝑡 + 1)-th snapshot are discussed. This indicates that the connections between
users and RUs have varied, and the assigned pilot sequence needs to be revised. The robust pilot
assignment method indicates the specified pilot sequence strategy can be durable in a dynamic scenario;
that is, even if the topology structure has changed, the pilot sequence could continue to operate on the
updated graph if the distributed unit (DU) has received updated pathloss information.
MDS based channel estimation:
On the RU side, each RU receives a combination signal from at most 𝑇𝑝 dominant pilot signals,
allowing for the estimation of 𝑇𝑝 dominant channels. Regarding the maximum distance separable
(MDS) code generation process, it has been demonstrated that the Reed-Solomon code satisfies MDS
code requirements. As one type of linear codes, Reed-Solomon and its derivatives are perhaps the most
widely used codes due to their numerous advantages [Rot06]. The generator matrix is defined as
𝐶(𝐾, 𝑇𝑝 ) = [𝑐1 , 𝑐2 , . . . , 𝑐𝐾 ], and each element in this matrix is an integer number that demonstrates the
linear combination of pilots. In this way, the pilot sequence of the k-th user 𝒔𝑘 can be specified as
𝑇𝑝
𝒔𝑘 = ∑ 𝝓𝑗 𝑐𝑘𝑗 , (5-13)
𝑗=1
where 𝝓𝑗 is the orthogonal vector, and 𝑐𝑘𝑗 is the coefficient of j-th orthogonal vector. For example, in
Figure 5-34, it is consided that all dominating channels (links in the bipartite graph) must be estimated,
and the maximum dimension is same with the maximum number of edges connected to RU, which is 2
in the upper figure. A larger space than the initial snapshot is defined to enable the generated code to
be available for more snapshot, e.g., 𝑇𝑝 = 3. Then, generate the normalized MDS generated matrix as
𝒄1 1 0 0
𝒄 0 1 0 (5-14)
[𝒄2 ] = 0 0 1 .
3
7 1 6
𝒄4 [ 14 14 14]
Here, the 5th RU is demonstrated as the example. The received signal at the upper figure is specified as
𝑟5𝑢 = √𝜌𝑇𝑝 𝑔52
𝑢 𝑢
𝒄2 𝚽 + √𝜌𝑇𝑝 𝑔53 𝑢
𝒄3 𝚽 = √𝜌𝑇𝑝 𝑔52 𝑢
𝒄2 𝝓2 + √𝜌𝑇𝑝 𝑔53 𝒄3 𝝓3 , (5-15)
and once the topology changes, the received signal at the bottom figure becomes
𝑟5𝑏 = √𝜌𝑇𝑝 𝑔52
𝑏 𝑏
𝒄2 𝚽 + √𝜌𝑇𝑝 𝑔53 𝑏
𝒄3 𝚽 + √𝜌𝑇𝑝 𝑔54 𝒄4
7 1 6 (5-16)
𝑏 𝑏 𝑏
= √𝜌𝑇𝑝 𝑔52 𝒄2 𝝓2 + √𝜌𝑇𝑝 𝑔53 𝒄3 𝝓3 + √𝜌𝑇𝑝 𝑔54 ( 𝝓1 + 𝝓2 + 𝝓3 ),
14 14 14
As a result, both channels of upper and bottom figure at 5 RU can be resolved by using minimum
th
mean square error (MMSE) or least square (LS) channel estimation methods. The MDS-based pilot
assignment approach modifies the dynamic pilot assignment from dynamically assigning the pilot
sequence to the selected solvable channel for each RU. In other words, the RUs receive sufficient signal
to solve 𝑇𝑝 dominating channels, which may be predicted using a designed estimator according to the
graph structure e.g., MMSE. In fact, the criteria of the MDS-based pilot sequences are substantially less
stringent than those of general orthogonal connection-based pilot assignment techniques, yet the MDS-
based pilot sequence is remarkably robust. On the one hand, the dimension of the MDS code determines
the highest capability of estimating the channel from the received signal, however, the property of the
MDS has no bearing on the connection structure. Thus, a designed estimator may estimate the desired
channels. Moreover, this non-orthogonal pilot sequence eliminates the restriction on the number of
users in this cell. On the other hand, because each RU receives the non-orthogonal pilot sequence
concurrently from all users throughout all timeslots, the interference from all other weak links cannot
be disregarded. The detailed algorithm has been shown in here.
Simulation result
An open square situation is considered, where 𝑀 = 200 APs and 𝐾 = 100 users with a single
antenna are randomly distributed throughout a 1 km × 1 km square region. Modelling the large-scale
fading coefficient 𝛽𝑚𝑘 as [BD16], the MDS-based method is compared with existing methods, such as
semi-random [NAY+17], information-based k-means (IBKM) [CZE+21], and the sequential maximum
weighted induced matching (sMWIM) [YYC22], respectively. Regarding the dynamic system concept,
the topological structure is assumed to be dynamic. In this work, 𝑁 = 100 snapshots are used in which
the RU and user locations are randomly regenerated for each snapshot. Due to the sparse channel of the
mmWave model, it is assumed that 10% connection in sMWIM unless otherwise specified. The pilot
dimension 𝑇𝑝 of the semi-random method sets to 6, and for the RU selection of structured policies, each
RU maximizes the selection of 6 users. As there are more users in this cell, the value 𝜅 = 3 of sWMIW
is adopted in this context. To generate the bipartite graph and determine the MDS-based method
parameter, it is assumed that each user links 𝑁𝑖𝑛 = 8 APs with stronger connections in the bipartite
graph, and the pilot dimension for MDS is 𝑇𝑝 = 8.
Figure 5-35: CDF of the downlink achievable rate Figure 5-36: Rate versus pilot dimension under
per user under the dynamic scenario. the dynamic scenario.
To evaluate the fairness of the baselines, it is assumed that the users utilize the same pilot sequence
despite the changing topological structure, and that the DU modifies the pathloss information that is
used for MMSE channel estimation. Figure 5-35 depicts the cumulative distribution function (CDF) of
the downlink achievable rate per user of our MDS-based method in comparison to other benchmarks.
In this figure, the performance of the proposed MDS-based approach exceeds that of the considered
benchmarks, indicating that this method can sufficiently detect dynamic impacts. In contrast to the
outcome of a static condition, semi-random provides superior robustness performance to the sWMIW
and IBKM methods. This is owing to the fact that semi-random is not dependent on the system structure.
To illustrate the performance clearly, Figure 5-36 depicts the sum rate versus the pilot dimension 𝑇𝑝 for
all pilot assignment algorithms under dynamic conditions The MDS-based technique has the highest
sum rate when the pilot dimension is 𝑇𝑝 = 8, while the semi-random method has the best performance
when the pilot dimension is small, i.e., 𝑇𝑝 = 2,4. Unlike the performance under static conditions, the
robustness of other techniques that depend on the connection structure presents a worse performance.
and requirements for 5G NR BSs comprising the transmitter and receiver characteristics and
performance testing. Thus, a technology must fulfil all mandatory RF tests and requirements defined
by 3GPP to be considered NR compliant.
DSO
MZM
SSMF
PC VOA EA ESA
LS PD
AWG
where P is the transmit power at the BS and g is the amplification gain applied by the NCR. Also, in
such a DL system,
𝛾𝑖,𝑗 = 𝑟𝑗 𝐻𝑖,𝑗 𝑡𝑖 , (5-18)
represents the combined effective channel from node i (the BS or the NCR depending on the considered
link) to node j (the NCR or the UE depending on the considered link) with transmit and receive
beamforming, i.e., 𝑡𝑖 and 𝑟𝑗 , respectively. Also, 𝐻𝑖,𝑗 denotes the channel between transmitter i = {b, n},
and receiver j = {n, u}. Assuming perfect channel state information (CSI), the transmit and receive
beamformer can be obtained by, e.g., maximum-ratio transmission. The total noise power received by
the UE can be calculated as
𝑁𝑃 = 𝜎12 + 𝜎22 𝑔𝛾𝑁𝑈 , (5-19)
where 𝜎12 and 𝜎22 are the variance of the AWGN at the UE and NCR, respectively. In this way, the
received signal-to-noise ratio (SNR) at the UE is S/NP, and the achievable rate is given by
𝑅 = 𝐵 log 2 (1 + 𝑆𝑁𝑅), (5-20)
g = 100 dB
with HWI
without HWI
102
1 2 3 4 5 6 7 8
Number of NCR antennas
Figure 5-39: An illustration of how NCR helps Figure 5-40: The required number of RIS elements
the network in the presence of blockage. to achieve the same performance as in the NCR-
assisted network.
Figure 5-40 shows the required number of RIS elements to achieve the same throughput as in the NCR-
assisted setup. Here, the results are presented for both cases with ideal and HWI-affected RIS models.
Also, the transmit power at the BS and the maximum output power constraint at NCR are set to 43 dBm
and 40 dBm, respectively with 90 dB/100 dB NCR amplification gain. As demonstrated, to achieve the
same performance as in the NCR-assisted setup, a large number of RIS elements are required, and the
required numbers increase with the NCR amplification gain. For instance, compared to the cases with
7 antennas at the NCR with amplification gain of 100 dB, around 200 RIS elements are required to
obtain equivalent performance as in the NCR-assisted setup, if the RIS is ideal. With HWI, however,
the required number of RIS elements increases significantly, for instance, to more than 1000 RIS
elements in the considered example.
"In the study, the concept and structure of NCRs were introduced, along with conceptual simulations
and a comparison of the performance of NCR-assisted networks with cases using RISs.
• The path cost of a UE is determined by the maximum cost of any segment in the path between
UE and Core. The motivation for this choice is that the throughput of an end-to-end link
between a UE and the core is limited by the link segment with the lowest throughput.
The considered network settings are summarized in Table 5-17.
Figure 5-43: Outage for Single and Dual connectivity in the UE-RU link. Outage is caused by
shadowing/blocking which is a function of distance. Link-loss in mesh backhaul caused by interference.
The number of APs with core connections was 9. Network simulation parameters are summarized in
Table 5-17.
The graphs shown in Figure 5-41 show that for a low RU density, increasing the number of RUs results
in an increased number of outage rate is improved since the distance between RUs an UEs is reduced.
Allowing for two connections compared to only one from a UE, increases the robustness to shadowing
and thus reduces the outage. The effect of adding more RUs to the network saturates for large number
of RUs. Interference between the backhaul links starts to affect the reliability. Note that this simulation
was performed for a fixed number of core-connections, even when number of APs was increased.
Intuitively, adding core connections, is the same as reducing the number of mesh segments (‘cell size’).
In Figure 5-43, the number of hops required to reach the core from a UE serving RU is illustrated. For
a given number of core connections, increasing the number of RUs does not reduce the number of hops
further. For the chosen simulation setup moving with 4 Core connection RUs, increasing the number of
RUs beyond 25, does not reduce latency (number of hops) anymore. Interestingly the distribution of
hop-count starts to stabilize, which means that number of active links in the system starts to increase
linearly with the number of RUs despite them operating in a mesh configuration. Adding more RUs
mainly starts to improve reliability on the access links. Increasing the number of RUs with core-
connections to 16, decreases the number of hops to reach the core network even for a higher number of
RUs in the network. Compared to the same number of RUs and 4 core-connections, the number of hops
reduces, and thus inter-RU interference in the mesh.
Figure 5-44: Distribution of number of hops for a different of core-DU/RUs and a different number of
RUss. The algorithm used for calculating routes through the mesh was an iteratively applied Dijkstra
optimizing for link SINR.
where 𝑃𝑡 , ℎ𝑡,𝑟 , 𝐺𝑡,𝑟 , 𝐿𝑡,𝑟 denote the transmit power, small scale fading and pathloss according to 5GCM
Uma close-in model respectively. The antenna gain model represents the sectored-pattern antenna array
model described in [AMK+19], [MMC+20] where the main lobe antenna gain is significantly higher
than the side lobe gain. In IAB networks the UEs can be served either by the child IAB nodes or IAB
donor. In this context, following open access strategy and maximum received power rule the UEs are
associated with the serving nodes. In a similar manner, the backhaul link association is determined
based on the minimum pathloss rule. Moreover, the aggregated interference observed by the UE 𝑢 due
to the neighbouring interferers is given by
(5-22)
𝐼𝑢 = ∑ 𝑃𝑖 ℎ𝑖,𝑢 𝐺𝑖,𝑢 𝐿𝑖,𝑢 (||𝑥𝑖 − 𝑥𝑢 ||)−1 ,
𝑗∈𝑋𝑖,𝑢 {𝑤𝑢 }
where 𝑖 represents the nodes excluding the associated node 𝑤𝑢 of UE 𝑢. Similarly, the aggregated
interference on the backhaul links can be calculated as given in [MMM+21]. The achievable UE rates
are then calculated using Shannon’s capacity considering sufficiently long codewords, which is an
acceptable assumption in IAB networks [MMF+20]. Using the achievable UE data rates 𝑅𝑈 , the study
performs constrained deployment optimization such that service coverage given by
is maximized where 𝜌 is the defined threshold data rate. Here, the IAB placement with minimum inter-
IAB distance requirement and IAB placement in the presence of constrained areas requirement are
optimized using Algorithm 1 and Algorithm 2 respectively.
Here, greedy-based algorithms for optimizing the IAB placement with minimum inter-node distance
constraint and geographical constrains are proposed respectively. Thereby, using a rejection-sampling
method where multiple possible solutions are considered while satisfying the aforementioned
constraints, the service coverage probability is maximized.
Performance Evaluation:
Performance of the optimized IAB system with deployment constraints at 28 GHz is analysed in a DL
outdoor scenario, for different number of regularly deployed UEs, IAB child nodes and IAB donors.
The general simulation parameters are given in Table 5-18.
Table 5-18: Simulation parameters.
Parameters Model Specification
Career Frequency 28 GHz
Bandwidth 1 GHz
Blocking 500 km-2
Noise power 5 dB
Figure 5-46: Service coverage as a function of the minimum distance constraint between the nodes for
different UE data rate thresholds.
Figure 5-46 depicts a setup where the IAB donor is located at the centre and 3 child nodes are placed
symmetrically, equidistant from the donor. As demonstrated in Figure 5-46, in such setup, the service
coverage increases with the node separations, up to a point around 550 m. This is characterized by the
effect of decreased interference between the nodes while having proper coverage in the area. Thereafter,
since the UEs can be located in between too far nodes, the coverage starts to slightly drop. In this way,
there is an optimal distance between the nodes maximizing the coverage. Finally, the coverage
decreases significantly with increased UEs minimum rate requirements due to the increased requirement
of resources/nodes to serve such needs.
Figure 5-47 presents the results for the cases where the nodes locations are obtained only by considering
the minimum distance between them or when the blockages and the backhaul links’ qualities are also
considered in the optimization. As depicted, the service coverage drops when the constraint becomes
tighter, however, for all considered range of constraints, compared to hexagonal deployment,
constrained deployment optimization increases the network coverage significantly. Indeed, prior
knowledge of the blockages’ locations helps to improve the network deployment optimization,
especially when the UE density increases. Also, as observed, the effect of inter-node distance constraint
on the coverage increases with the UE density. Here, the results are presented for the cases where the
IAB donor has a main lobe antenna gain of 24 dBi and child IAB nodes have a main lobe antenna gain
of 18 dBi while the side lobe gains are significantly lower, i.e., -2 dBi.
Figure 5-48: Service coverage as a function of the Figure 5-49: Service coverage as a function of the
minimum distance constraint between the nodes in radius of the constrained areas for different
different UE densities for blocking aware and deployment scenarios.
blocking unaware optimization scenarios.
Figure 5-48 verifies the effect of geographical constraints, on the coverage of IAB networks.
Particularly, the coverage of the deployment optimized IAB networks in the cases where, the IAB nodes
are unable be placed in constrained areas either due to geographical or regulatory restrictions is studied.
Here, the results are presented for a network consisting of five circular constrained areas of radius c. As
demonstrated for low geographical constraints, network performance is not affected by the deployment
constraints. However, with the presence of large area constraints, the performance, i.e., service coverage
decreases. This is intuitively because there is an increased chance of low coverage for users within the
constrained areas when there are larger constraints for IAB placement. Also, since the IAB nodes get
packed outside the constrained areas, interference levels for the UEs outside the constrained areas
increases resulting in a further decrease in coverage. Finally, proper network planning boosts the
coverage significantly compared to random deployment. Also, compared to the case with child IAB
nodes distributed randomly in the unconstrained areas, the effect on the coverage is less severe when
the network is optimized for geographical constraints.
BS 0m BS BS 0m
0) (29 , 0, 3) (-15 , 0, 10) (0 , 0) (15 , 0, 10)
BS1 BS2
-15 m -15 m
y
30 m -30 m -60 deg 0m 30 m y -30 m 0m 30 m
et er = 60 x Widt h_met er = 60 Width_meter = 60
x
Figure 6-1: Deployment scenario with 2 BSs (left) and 4 (right) BSs illuminating the area from the ceiling.
The two scenarios are depicted in Figure 6-1. Each base station is equipped with a uniform planar array
(UPA) with 8 x 8 antenna elements. Analog beam forming is assumed. To enable more than one radio
frequency (RF) chains, the antenna can be divided into two or four subpanels with 8x4 and 4x4
elements, respectively. Each subpanel has one RF chain, so that 1, 2 or 4 RF chains are simultaneously
active, and the respective number of simultaneous users can be served. Each subpanel in each
configuration uses a grid-of-beams (GoB) with 9 x 9 beams covering an angle range of +/- 60° in both
x and y dimensions. The subpanel configurations are shown in Figure 6-2.
With a system simulation the total throughput achievable within the deployment area, assuming 4, 8,
16 and 32 simultaneously served UEs, has been compared for the 2-BS and the 4-BS scenario (Figure
6-3). A total transmission (TX) power per BS of 41 dBm for the 2-BS scenario and 35 dBm for the
4-BS scenario has been applied. Each subpanel is assumed to serve 1 UE. For comparison also, the
values for a fully digital GoB solution, which allows serving multiple UEs simultaneously, have been
simulated. However, it should be noted that this is a theoretical configuration, since the number of RF
chains for a full digital solution with large arrays might be prohibitively large. The purpose of this
analysis is to investigate the power consumption for serving a specific indoor scenario. Hence, the
Dissemination level: public Page 113 / 127
Hexa-X Deliverable D2.3
throughput in the total area is compared for the two configurations. Target is to find the most energy
efficient deployment for a specific area with different number of simultaneous UEs. Therefore, the
results are evaluated for the whole playground area and not “per cell”, because also the number of cells
is varying.
K=4
A P=4
Ms=N/4 4
N=no. of elements
K=2
2
P=2
C Ms=N/2
N=no. of elements
K=1
1
P=1
D Ms=N
N=no. of elements
300
250
Area TP [Gbit/s]
0
4 8 16 32
no. of simultaneously scheduled UEs in playground
Figure 6-3: Area throughput for different number of simultaneously served UEs.
In Figure 6-4 the 4-BS scenario is shown, and the number of subpanels is varying. With 1 subpanel up
to 4 UEs can be served simultaneously in the area (dark blue bar). With increasing number of UEs, the
throughput is distributed among the higher number of UEs. When increasing the number of subpanels,
an increasing number of UEs can be served, thus increasing the area throughput up to 8 UEs for 2
subpanels (light blue bars), and up to 16 UEs in case of 4 subpanels (magenta bars). The dark blue,
green and grey bars correspond with bars of same colours in Figure 6-3.
4 BS, Area TP, same no. of x-pol UE [Gbit/s], SBPL
350
300
250
Cell TP [Gbit/s]
0
4 8 16 32
no. of simultaneously scheduled UEs in playground
Figure 6-4: Area throughput for different number of simultaneously served UEs.
When looking at the UE throughput cumulative distribution function (CDF) (exemplary for 1 subpanel)
over the deployment area shown in Figure 6-5, it can be seen that for 4 BS the area throughput CDF
shows higher throughput than with 2 BS (red curves more to the right than blue curves).
Figure 6-5: UE Throughput CDF for 2 and 4 BS over the whole deployment area, 1 subpanel and
different number of simultaneously served UEs.
The next step is to assess the related power consumption of these configurations and bring it into relation
to the throughput, to come to energy efficiency comparison. Basis for power consumption analysis is
the power consumption model derived for antenna array hardware and high bandwidth data conversion
as described in [HW21] and [HEX21-D22The power consumption of the RF frontend, including data
conversion, is assessed. However, the power consumption model does not yet include baseband
processing and is subject to further study. In Figure 6-6 the results for 2 and 4 BSs and also different
ADC resolutions of 8 and 12 bits are shown.
It can be seen that for the 4 BS case, where the number of PAs is higher but the TX power can be lower,
the overall power consumption is reduced although the area throughput has increased. Thus, a further
increase of the number of BS could potentially further reduce the power consumption. But at the same
time the number of RF chains increases, which leads to a trade-off and related optimum depending on
the specific scenario.
Figure 6-6: Exemplary power consumption for 2- and 4-BS scenarios and different ADC resolution.
6.2.1 Pathloss, angular and delay dispersions of outdoor and indoor radio
links
The path loss, delay and angular spread models of a wireless channel are most influential characteristics
of channel to radio links and systems. These parameters are estimated for the line-of-sight (LOS) and
non-line-of-sight (NLOS) links in the entrance hall and three outdoor scenarios namely suburban,
residential and city centre, described in [HEX21-D22, Section 7]. Two path loss models such as the
close-in (CI) reference free space reference and the floating initial pathloss model, also called alpha-
beta-gamma (ABG) models, are considered here. The path loss of a signal with frequency 𝑓𝑐 at distance
𝑑 based on the CI model is expressed as
𝑑 (6-1)
𝑃𝐿CI = 𝐹𝑆𝑃𝐿(𝑓𝑐 , 1 m) + 10𝑛 log10 ( ) + 𝜒𝜎CI ,
𝑑0
where 𝐹𝑆𝑃𝐿(𝑓𝑐 , 1 m) is the free space path loss of the signal with frequency 𝑓𝑐 at 1 m distance, 𝑛 is the
path loss exponent, and 𝜒𝜎CI is the shadow fading for the CI model, and the ABG model is given by
𝑓 (6-2)
𝑃𝐿ABG = 10𝛼 log10(𝑑) + 𝛽 + 10𝛾 log10 ( ) + 𝜒𝜎ABG
1 𝐺𝐻𝑧
where 𝛼 is the distance-dependent loss coefficient, 𝛽 is an offset coefficient, 𝛾 is the frequency-
dependent loss coefficient, and 𝜒𝜎ABG is the shadow fading for the ABG model. The omnidirectional
path loss for the indoor and outdoor scenarios, together with the fitted path loss models are plotted in
Figure 6-7.
Note that the outdoor scenario is an ensemble of path loss estimates from the three mentioned outdoor
scenarios since there are only limited NLOS links that can provide sensible fitting to the model when
individual scenario is considered. It can be noticed that both the CI and ABG models provide equally
good fit for LOS and NLOS links in both scenarios.
The mean 𝜇 and standard deviation 𝜎 of angular and delay spread values are listed in Table 6-1: Angular
and delay spread statistics in indoor and outdoor scenarios.. The angular spread values are based on
azimuth and zenith angle of departure (AoD/ZoD) and azimuth and zenith angle of arrival (AoA/ZoA)
estimated from measurement-based ray-launcher presented in [DKH21]. Similar to the path loss
modelling, the statistics of various spread values for the outdoor case include those from the three
outdoor scenarios.
Figure 6-8 Measured path powers, the beam power and found independent beam azimuth directions of an
example link.
Measured PADPs were evaluated using 10° half power beam width (HPBW), 10 dB dynamic range
below the strongest beam, 2 GHz bandwidth (BW), and 0.5 correlation threshold between beams. An
example PADP, beam power and identified beam directions are illustrated by blue circles, red curve
and orange squares, respectively, in Figure 6-8. The identified independent number of beams in 132
measured indoor transmitter (Tx) and receiver (Rx) locations are shown in Figure 6-9. Figure 6-10
depicts empirical cumulative distribution functions (CDF) of the number of beams in 132 indoor and
157 outdoor locations using both 10, and 20 dB dynamic range. Median values of the number of beams
are two in both environments using dynamic range of 10 dB, and using 20 dB they are five and three in
indoor and outdoor environments, respectively.
7 Conclusions
Sub-THz radio access network has the potential to realize use cases with extreme performance
requirements in terms of data rate and latency. However, designing radio hardware at sub-THz upper
Millimetre-wave frequencies (100-300 GHz) is challenging due to the complexity of multi-antenna
implementations, hardware impairments, transmit power limitations, and propagation channel
properties. A single radio configuration may not be suitable for all use cases and deployment scenarios,
especially, when considering sustainability metrics. The radio design option depends on factors such as
required data rate, communication range, and device type. Moreover, the physical layer modulation
scheme and frame structure should fulfil the air interface latency requirements. Accordingly, the report
proposes an approach for analysing radio systems, which involves considering multiple radio options
based on different use cases. First, the requirements for each use case are determined, such as data rate,
range, and number of users, and then used for deriving high-level technical requirements and
communication scenarios. For each scenario, multiple radio options are proposed for different
frequencies and device class.
The report discusses cellular-based and D-MIMO deployment options and highlights the potential of
D-MIMO technology to optimize for coverage, enhance the capacity and reliability by cooperative
transmission over multi-links. However, high data rates can increase processing costs for fronthaul and
backhaul. The report explores various approaches to address this challenge.
A technical guideline for radio frequency (RF) transceiver design is presented with detailed discussions
and evaluation example in each step. First, the technical requirements to achieve the data rate are
presented in terms of bandwidth (BW), signal-to-noise ratio (SNR), and the challenges of analogue-to-
digital convertors (ADC). Then, it discusses the requirements to achieve the SNR for a give range
considering the RF nonidealities, and waveform impact on the power amplifier (PA) backoff. In
addition to this analysis, it provides insight on the RF implementation aspects and strategies such as
channelization and exploiting the concept of carrier aggregation (CA) to optimize the system. In each
of this step, there can be several options to select from. Thus, other performance metrics, such as energy
consumption and complexity, need to be considered to reduce the design space.
Based on the analysis, it is concluded that, for sub-THz frequencies, a hybrid architecture based on
subarrays is a practical choice. The number of subarrays controlled with a full RF chain can be
optimized depending on the number of simultaneous users and available spatial beams by the channel.
Moreover, it enables the implementation of wide-band waveforms by aggregating multiple subarrays
over the air. Furthermore, it provides control capability for power consumption by means of switching
off the unneeded chains. However, insight analysis on beam management reveals that conventional
beam tracking is inefficient due to the large number of beams to be probed with sub-THz narrow beams.
Therefore, new approaches are needed to exploit available side information.
The report emphasizes the importance of hardware and channel modelling in the design process to
consider non-idealities, compensate for path loss and exploit spatial beams. Several models for PA and
phase noise are described and their accuracy is evaluated using measurement of lab components. A
companion with state-of-the art models is introduced, proposing an extension of those models to cover
wider frequency band above 100 GHz. In addition, the report highlights that channel measurements
alone are insufficient to evaluate performance, as they depend on channel sounding and environmental
factors. Therefore, data processing and the derivation of an extended channel model are critical for
design. The modelling work in this report is complementary to the results presented in the second
deliverable [HEX21-D2.2].
After designing the radio, the physical layer signal processing is evaluated considering one or more
hardware models focusing on frame error rate (FRE). In the section of the waveforms, one approach is
based on designing hardware-friendly modulation schemes, and another approach is based on
developing digital techniques to estimate and mitigate the impact of hardware. For instance, OFDM
variants such as DFTS-OFDM and SC-FDE are seen as potential candidates due to their low peak to
average power ratio (PAPR) and resistance to phase noise, especially with optimized constellations
schemes. Moreover, these waveforms can also be easily integrated in the sub-array architecture based
on spatial combining or carrier aggregation. On the contrary, zero-crossing modulation (ZXM) with 1-
bit quantization, as an energy-efficient waveform, needs to cope with the phase noise by estimation and
compensation.
A study in the report, focusing on achieving air interface latency, demonstrates that achieving 0.1 ms
requires a short transmission time interval (TTI). However, the processing latency raises as a major
challenge, especially, with extremely high data rates. This needs to be considered in the processing
architecture design, considering I.O speeds, memory access speed, and power consumption.
Several performance examples demonstrate that optimizing deployment scenarios, transceiver
configurations, and ADC resolution can significantly impact the power consumption and improve
performance. Thus, the future work should focus on holistic optimization of radio design considering
different degrees of freedom not only in the infrastructure but also in the signal processing, to ensure
meeting the sustainability goals, while fulfilling the extreme requirements of emerging use cases.
References
[CZE+21] S. Chen, J. Zhang, E. Bjö rnson, J. Zhang, and B. Ai, “Structured massive access
for scalable cell-free massive MIMO systems,” IEEE J. Sel. Areas Commun.,
vol. 39, no. 4, pp. 1086-1100, Apr. 2021.
[DKH21] M. F. De Guzman, P. Koivumäki and K. Haneda, "Double-directional
Multipath Data at 140 GHz Derived from Measurement-based Ray-launcher,"
in Proc. 2022 IEEE 95th Vehicular Technology Conference (VTC2022-
Spring), 2022.
[DHK23] M. F. De Guzman, K. Haneda and P. Kyösti, “Measurement-based MIMO
channel model at 140GHz (Version 1),” available:
https://doi.org/10.5281/zenodo.7640353.
[DPS20] Erik Dahlman, Stefan Parkvall, Johan Skold, 5G NR - The Next Generation
Wireless Access Technology, 2nd Edition, 1. Chapter 26, §26.2: “mmW
LO generation and phase noise aspects” - September 18, 2020 ISBN:
9780128223208
[DVB-S2X] ETSI EN 302 307 V1.2.1; Digital Video Broadcasting (DVB); Second
generation framing structure, channel coding and modulation systems for
Broadcasting, Interactive Services, News Gathering and other broadband
satellite applications. DVB-S2 Extensions (DVB-S2X)
[FDB+19] G. Fettweis, M. Dörpinghaus, S. Bender, L. Landau, P. Neuhaus, and M.
Schlüter, “Zero crossing modulation for communication with temporally
oversampled 1-bit quantization,” in 2019 53rd Asilomar Conference on
Signals, Systems, and Computers, Pacific Grove, CA, USA, Nov. 2019, pp.
207–214.
[GDS+22] P. R. B Gomes, G T. de Araujo, B. Sokal, A L. F. de Almeida, B. Makki, G.
Fodor, “Channel estimation in RIS-assisted MIMO systems operating under
imperfections”, [Online]. Available at: https://arxiv.org/pdf/2207.02700.pdf.
[GFF+21] U. Gustavsson, P. Frenger, C. Fager, T. Eriksson, H. Zirath, F. Dielacher, C.
Studer, A. Pärssinen, R. Correia, J. Matos, D. Belo, N. Borges Carvalho,
“Implementation Challenges and Opportunities in Beyond-5G and 6G
Communication,” in IEEE Journal of Microwaves, vol. 1, no. 1, pp. 86-100,
winter 2021 (invited).
[GMA+22] H. Guo, B. Makki, M-S. Alouini and T. Svensson, "High-rate uninterrupted
internet of vehicle communications in highways: Dynamic blockage avoidance
and CSIT acquisition," IEEE Commum. Mag., vol. 60, no. 7, pp. 44-50, Jul.
2022.
[GMM+22] H. Guo, et al. "A Comparison between network-controlled repeaters and
reconfigurable intelligent surfaces,” Nov. 2013. [Online]. Available at:
https://arxiv.org/abs/2211.06974.
[GS91] A. Ghorbani and M. Sheikhan, “The effect of solid-state power amplifiers
(SSPAs) nonlinearities on MPSK and m-QAM signal transmission,” in6th Int.
Conf. on Digital Processing of Signals in Communications, Loughborough UK,
1991, p. 193–197.
[GSD+22] F. Gast, M. Schlüter, M. Dörpinghaus, H. Halbauer, and G. Fettweis, “Phase
noise tracking for receivers with 1-bit quantization and oversampling,” in ICC
2022 - IEEE International Conference on Communications, May 2022, pp.
1361–1366.
[HEX21-D11] Hexa-X project deliverable D1.1: 6G Vision, use cases and key societal values,
Eds. M. Hoffmann, M. Uusitalo, M.-H. Hamon, and B. Richerzhagen, Feb.
2021.
[HEX21-D12] Hexa-X project deliverable D1.2 Expanded 6G vision, use cases and societal
values, Eds. M. Hoffmann, M. Uusitalo, M.-H. Hamon, B. Richerzhagen, G.
D’Aria, A. Gati and D. Lopez, Apr. 2021.
[HEX21-D12] Hexa-X project deliverable D1.2 Expanded 6G vision, use cases and societal
values, Eds. M. Hoffmann, M. Uusitalo, M.-H. Hamon, B. Richerzhagen, G.
D’Aria, A. Gati and D. Lopez, Apr. 2021.
[HEX21-D13] Hexa-X project deliverable D1.3: Targets and requirements for 6G - initial E2E
architecture, Eds. B.M. Khorsandi, M. Hoffmann, M. Uusitalo, M.-H. Hamon,
B. Richerzhagen, G. D’Aria, A. Gati, and D. Lopez, Feb. 2022.
[HEX21-D21] Hexa-X project deliverable D2.1 Towards Tbps Communications in 6G: Use
Cases and Gap Analysis, Eds. K. Roth, N. Michailow, A. Pärssinen and H.
Halbauer, Jun. 2021.
[HEX21-D22] Hexa-X project deliverable D2.2 Initial radio models and analysis towards
ultra-high data rate links in 6G, Eds. Marko E. Leinonen, Dec. 2021.
[HEX21-D31] Hexa-X project deliverable D3.1: Localisation and sensing use cases and gap
analysis, Eds. H. Wymeersch, Dec. 2021.
[HEX23-D33] Hexa-X project deliverable D3.3: Final models and measurements for
localisation and sensing, Jul. 2023.
[HEX21-D51] Hexa-X project deliverable D5.1: Initial B5G/6G Architectural Components
and Enablers, Eds. M. Ericson, H. Flinck, P. Vlacheas, and S. Wänstedt, Dec.
2021.
[HH97] M. Honkanen and S.-G. Haggman, “New aspects on nonlinear power amplifier
modeling in radio communication system simulations,” in Proceedings of 8th
International Symposium on Personal, Indoor and Mobile Radio
Communications - PIMRC ’97, vol. 3, 1997, pp. 844-848 vol.3.
[HNH+07] A. Hati et al., "Vibration Sensitivity of Microwave Components," 2007 IEEE
International Frequency Control Symposium Joint with the 21st European
Frequency and Time Forum, Geneva, Switzerland, 2007, pp. 541-546, doi:
10.1109/FREQ.2007.4319131.
[HNH09] A. Hati, C. W. Nelson and D. A. Howe, "Vibration-induced PM and AM noise
in microwave components," in IEEE Transactions on Ultrasonics,
Ferroelectrics, and Frequency Control, vol. 56, no. 10, pp. 2050-2059, October
2009, doi: 10.1109/TUFFC.2009.1288.
[HW21] H. Halbauer, T. Wild, “Towards Power Efficient 6G Sub-THz Transmission,”
Joint EuCNC & 6G Summit, June 8-11 2021, Porto, Portugal
[IBS+98] M. Ibnkahla, N. Bershad, J. Sombrin, and F. Castanie, “Neural network
modeling and identification of nonlinear channels with memory: algorithms,
applications, and analytic models,” IEEE Transactions on Signal Processing,
vol. 46, no. 5, pp. 1208–1220, 1998.
[IEE09] IEEE, “60 GHz Impairment Modeling, Doc. 802.11-09/1213r1,” IEEE Std
802.11-09, pp. 1–29, Nov 2009.
[IEE16] IEEE, “IEEE TGay Evaluation Methodology, Doc. 802.11-15/0866r2,” IEEE
Std 802.11-15, pp. 1–27, Oct 2016.
[Imm04] K. A. Schouhamer Immink, Codes for Mass Data Storage Systems, 2. ed.,
Shannon Foundation Publishers, 2004.
[ITU P.2040-2] ITU-R recommendation P.2040-2, “Effects of building materials and structures
on radiowave propagation above about 100 MHz”, Jul. 2021.
[KGH+22] P. Kyösti, M. F. De Guzman, K. Haneda, N. Tervo, A. Pärssinen, ”How many
beams does sub-THz channel support?,” in IEEE Antennas and Wireless
Propagation Letters, vol. 21, no. 1, pp. 74-78, Jan. 2022.
[KKP+14] M. R. Khanzadi, D. Kuylenstierna, A. Panahi, T. Eriksson, and H. Zirath,
“Calculation of the performance of communication systems from measured
oscillator phase noise,” IEEE Trans. Circuits Syst. I, vol. 61, no. 5, pp. 1553–
1565, May 2014.
[LDF15] L. Landau, M. Dörpinghaus, and G. Fettweis, “Communications employing 1-
bit quantization and oversampling at the receiver: faster-than-Nyquist signaling
and sequence design,” International Conference on Ultra Wideband (ICUWB),
Montreal, Canada, Oct. 2015
[LJT+19] M. E. Leinonen, M. Jokinen, N. Tervo, O. Kursu and A. Pärssinen, "System
EVM Characterization and Coverage Area Estimation of 5G Directive mmW
Links," in IEEE Transactions on Microwave Theory and Techniques, vol. 67,
no. 12, pp. 5282-5295, Dec. 2019, doi: 10.1109/TMTT.2019.2951773.
[LMN+22] M. E. Leinonen, J.-P. Mäkelä, K. Nevala, N. Tervo, and A. Pärssinen,
“Repeatability of 220 - 330 GHz variable waveguide attenuator with frequency
extenders for 6G measurements,” in 2022 98th ARFTG Microwave
Measurement Conference (ARFTG), 2022, pp. 1–4
[LNT+21] M. E. Leinonen, K. Nevala, N. Tervo, and A. Pärssinen, “Linearity
measurement of 6g receiver with one transmission frequency extender
operating at 330 GHz,” in 2021 96th ARFTG Microwave Measurement
Conference (ARFTG), 2021, pp. 1–4.
[LR11] C. Liang and B. Razavi, "Transmitter Linearization by Beamforming," in IEEE
Journal of Solid-State Circuits, vol. 46, no. 9, pp. 1956-1969, Sept. 2011, doi:
10.1109/JSSC.2011.2148530.
[LTJ+19] M. E. Leinonen, N. Tervo, M. Jokinen, O. Kursu and A. Pärssinen, "5G mm-
Wave Link Range Estimation Based on Over-the-Air Measured System EVM
Performance," 2019 IEEE MTT-S International Microwave Symposium
(IMS), Boston, MA, USA, 2019, pp. 476-479
[LWZ20] Y. Li, X. Wang and A. Zhu, "Sampling Rate Reduction for Digital Predistortion
of Broadband RF Power Amplifiers," in IEEE Transactions on Microwave
Theory and Techniques, vol. 68, no. 3, pp. 1054-1064, March 2020, doi:
10.1109/TMTT.2019.2944813.
[LXA+21] V. Lampu et al., "Energy-Efficient Array Transmitters Through Outphasing
and Over-the-Air Combining," 2021 55th Asilomar Conference on Signals,
Systems, and Computers, 2021, pp. 608-615.
[MMC+20] C. Madapatha, B. Makki, C. Fang, et al., “On integrated access and backhaul
networks: Current status and potentials,” IEEE open j. Commun. Soc., vol. 1,
pp. 1374–1389, Sep. 2020.
[MMH+22] C. Madapatha, B. Makki, H. Guo, and T. Svensson., Constrained deployment
optimization in integrated access and backhaul networks. arXiv preprint
arXiv:2210.05253.