1 Introduction

Almost all of the current wireless communication takes place in a relatively small region of the electromagnetic spectrum below 6  GHz. This region has been allocated by government agencies around the world for commercial, civilian, military, public safety, and experimental use. However, the proliferation of devices and services that use or depend on wireless technologies has caused an ever-increasing discrepancy between the demand and the available bandwidth or degrees of freedom (DoF) This discrepancy termed spectrum crunch, if not addressed, will lead to lower data rates and reduced quality of service. Spectrum crunch will become even more acute when Internet of Things39 and device-to-device communication traffic are added to the already overloaded networks.

Increasing the DoFs is the only option for the next-generation (5G) wireless systems. The use of millimeter-wave (mmWave) enables an increase in DoFs by adding more bandwidth due to the availability of large (in the order of a few GHz) unlicensed spectrum between 30 and 300 GHz. However, as explained by Friis’ law41, signals transmitted in mmWaves have high isotropic pathloss, i.e., they decay at a much higher rate with the traveled distance. This leads to a reduced communication range compared to sub-6 GHz systems. Furthermore, mmWaves exhibit characteristics resembling the visible light. For example, they have high penetration loss through most material and hence are easily blocked by the surrounding objects. mmWave systems can overcome these shortcomings by employing beamforming (BF) That is, arrays of multiple antenna elements are used to extend the communication range and avoid obstacles in the environment by directing the signal energy in an intended direction.

However, the reliance of mmWave communication on directional links through beamforming poses new challenges that do not exist in wireless systems over the microwave bands. The transmitter (Tx) and the receiver (Rx) must first discover each other directions before they can start the data communication. This process, known in cellular systems as initial access, is generally performed omni-directionally (or using very wide beams) in the sub-6 GHz bands. But, due to the high path loss, if mmWave systems were to follow the same paradigm, the range of mutual discovery would be much smaller than the range where directional high-rate communication would be possible. Therefore, mutual discovery must be performed in a directional manner.

The directional discovery phase can last for a long time when the base station and the user employ arrays with many antenna elements forming narrow beams. While searching for the base station, the battery-limited user is always “on” burning energy. We show here that this energy consumption can be reduced significantly by employing a low-resolution fully digital front-end on the user side. The reason is, beamforming through a digital front-end reduces the discovery latency (or delay) by an order of magnitude compared to an analog front-end. Hence, the user is “on” for a shorter span of time, leading to considerable energy savings.

While the focus of this paper is on directional discovery in initial access, directional discovery is expected to be triggered in other phases of mmWave communication too. For example, in recovery from link failures, which will be frequent due to the sensitivity of the mmWave links to obstacles and changes in the environment. Handovers to a new base station will also be frequent since mmWave cells will be smaller in size compared to the current ones, and a mobile user may go through more cells for the same distance traveled. More importantly, since the battery-dependent devices operate at such high frequencies and bandwidth, more aggressive use of sleep/idle mode (discontinuous reception) will be critical from an energy consumption standpoint. For each of these operations, the user device must often (re-) discover the direction to the connected and neighboring base station(s). Therefore, it is extremely crucial that directional discovery and beam alignment is fast and energy-efficient.

1.1 Contributions

In this paper, we look into the problem of latency and energy consumption in directional discovery in mmWave systems during initial access. Our focus is on 5G cellular systems where the issue of communication range is more important and challenging than short-range mmWave communication, e.g., 802.11ad WiFi. There are two key take-away points in our work:

  1. 1.

    Digital beamforming results in both low latency and low energy consumption during initial discovery compared to analog beamforming.

  2. 2.

    Employing low-resolution analog to digital converters (ADCs) in fully digital front-ends can achieve low latency and even lower power consumption for both control signaling and data transmissions.

Mutual directional discovery requires the transmitter and the receiver search in their surrounding angular space at a granularity inversely proportional to the size of their antenna arrays: larger antenna arrays result in narrower beams potentially leading to higher discovery delays. In a cellular environment, this latency will also affect user handovers between different base stations, alternative link discovery in case of blockage, and the “idle mode” to “on mode” circles.

We show that due to the large latency in directional discovery incurred by analog beamforming, its energy consumption is greater not only than low-resolution digital but also than high-resolution digital. This difference in energy consumption increases with the size of the antenna arrays. When a 4-by-4 antenna array is used, analog beamforming burns as much as six times more energy than digital. This is due to the fact that analog beamforming needs more time to sample an angular domain that increases in size with the number of antennas.

Leveraging our previous work11,12, we establish a relationship between the beamforming architecture and the mutual discovery delay within the context of 3GPP initial access for mmWave networks5. Specifically, we show that between analog and digital beamforming, digital outperforms analog by a large margin—in the order of 300–900 ms. Similar to20,21, we detail and compare the power consumption of four beamforming architectures, namely, analog, digital, low-resolution digital, and hybrid, by assessing the components and devices they are comprised of. It is known that by reducing the resolution of analog-to-digital converters (ADC) we can significantly reduce the power consumption of a fully digital beamforming circuit (henceforth referred to as low-resolution fully digital). However, we show that while this reduction comes at a penalty of less than \(1 ~ \rm{dB}\) SNR in the low SNR regime (up to \(10 ~ \rm{dB}\))20, the discovery delay is kept at the same low levels as digital beamforming with high resolution—20–80 ms. Thus, low-resolution fully digital beamformers outperform analog as they can be power efficient during both data transmissions as well as signaling control messages.

Interestingly, in most studies related to mmWave systems, analog beamforming has been preferred over digital due to its low device power consumption. However, as we show, when discovery delay in initial access is taken into account, analog beamforming can burn multiple times more energy than any other alternative.

1.2 Related Work

Due to the reliance on highly directional links at mmWave frequencies, efficient beam management is key to establishing and maintaining a reliable link. A critical component of beam management is the beam discovery procedure for initial access. Current technical literature assumes an analog or hybrid front-end which limits the number of usable spatial streams. For instance,18,37 present a heuristic method to design beam patterns for initial beam discovery for hybrid beamformers. Raghavan et al. 40 propose to vary the beam widths depending on the user link quality. They show that users with good link (high SNR) can be detected with wide beams leading to a decrease in detection delays. In another direction of research,29,31,53 aim at designing optimal beam code books, i.e., sets of directions for connection establishment. Another important area of research has been the use of side information; for instance location information, channel quality measurements at microwave bands, etc., have been studied in6,14,15,19,23,25,34,38,51. Moreover, the work in44 proposes the use of online machine learning algorithms for beam detection for vehicular communications at mmWave. The works in33,42,54 use game- theoretic methods, whereas in27,48, the authors use genetic algorithms for initial beam discovery.

In9, the authors present a theoretical analysis of the trade-off between spending resources for accurate beam alignment on the one hand and using them for actual data communication on the other. More interestingly, the analysis in9 shows that with large coherence block lengths, exhaustive beam search outperforms hierarchical search. This understanding is reflected in the current 3GPP specifications2,3,4,5 on initial access where initial beam discovery and alignment is achieved through exhaustive search. We will discuss the 3GPP new radio beam search procedure briefly in Sect. 2.2. A detailed overview of the beam management procedure for 5G systems can be found in26.

Critical to beam discovery and beam alignment is the efficient signaling of pilot or synchronization signals and channel estimation. To this end, the work in22 leverages the sparsity of the mmWave channel. The authors use a compressed sensing framework for estimating the number of measurements necessary for estimating the channel covariance matrix for beam/angle detection with high confidence. Similarly, in47 compressed sensing is used for fast angle of arrival/departure and second-order channel statistics estimation. In52, a compressed sensing-based algorithm robust to frequency offsets and phase noise is presented. The sparsity of the mmWave channel is exploited in13, where a novel algorithm based on multiple-armed bandit beam selection for both initial beam alignment and beam tracking is proposed. In28, Hashemi et al. exploit the channel correlation to reduce the searching space and subsequently, the delay of beam discovery.

To better understand the interplay between the hardware (and power) constraints at mmWave and the beam discovery delay, in our previous work11,12, we presented a comparative analysis of analog and digital beamforming in terms of synchronization signal detectability and delay. We show that digital beamforming, even with low-resolution quantizers, performs dramatically better compared to analog or hybrid. Furthermore, our recent work,20,21, presents a thorough study of various beamforming architectures at both the transmitter and the receiver sides. There, we show that employing low-resolution digital beamforming introduces negligible loss in the system rate under realistic mmWave cellular assumptions. Based on this, we argue that as low-resolution fully digital front-ends can have low control delays while having the same power efficiency as analog and hybrid beamformers for the data plane. To the best of our knowledge, this is the first work where energy consumption during initial beam discovery has been studied. In this work, we show that a low-resolution fully digital beamformer is energy-efficient in both control plane in general, initial access in particular, and data plane.

The rest of the paper is organized as follows. In Sect. 2 we give an overview of the beam discovery problem and present the 3GPP new radio (NR) discovery procedure. In Sect. 3, we present the system model. We derive a correlation-based detector for beam discovery, present four mmWave beamformers and comment on the power consumption of each one of them. In Sect. 4, we evaluate through simulations the performance of the analog and high/low-resolution beamformers in terms of discovery delay and energy consumption.

Notation We use small boldface letters \(({\mathbf {h}})\) to denote vectors and boldface capital letters \(({\mathbf {H}})\) for matrices. Conjugate and conjugate transpose are denoted with \({\mathbf {h}}^*\) and the \(\hbox {L}{-}2\) norm with \(\Vert {\mathbf {h}} \Vert\) and \(\Vert {\mathbf {H}} \Vert\) for vectors and matrices, respectively.

Figure 1:
figure 1

Beam discovery in mmWave bands. Panels read from left to right and top to bottom. a Search over narrow sectors: the receiver, larger shape, and the transmitter, smaller shape, sequentially visit each sector combination until mutual discovery. b Hierarchical search. The receiver gradually narrows its beams in a series of steps.

2 Beam Discovery Through Beamsweep

Consider a transmitter–receiver pair operating in mmWave bands. Since they communicate through directional links, they both must discover each other in the spatial domain before the data communication can begin. An intuitive approach to this problem is to divide the spatial domain around the transmitter (or receiver) into multiple non-overlapping sectors, as presented in Fig. 1a. Each of these sectors corresponds to a transmitter or receiver beam. All possible combinations of the pair of transmitter–receiver beams are termed as the beamspace. Increasing the number of antennas makes the beams narrower, which increases the beamforming gain but also makes the size of the beamspace larger. This poses a fundamental problem for mmWave systems. On the one hand, narrow beams are necessary for usable link budgets. On the other, beam discovery, i.e., finding the two matching beam pairs at the transmitter and the receiver, a difficult problem.

In this work, we assume a stand-alone communication model where the transmitter and the receiver are not assisted by out-of-band information regarding timing or position. The non-stand-alone model has been investigated in6,14,15,19,23,25,34,38,51 and is not discussed in this work. There are two main approaches to beam discovery/selection under this assumption. One, called beamsweep, requires both the transmitter and the receiver to exhaustively search over the entire beamspace by measuring the received power for every possible transmitter–receiver beam pair. In another, the receiver side starts listening on the channel with the widest possible beam and step by step converges to the narrowest one. This is called hierarchical search, see Fig. 1b.

Both these techniques assume that a known signal called the synchronization signal in cellular is transmitted periodically, and the receiver will have to determine the direction in the beamspace where the incoming signal is stronger. While hierarchical search in principle is superior to beamsweep in terms of search delay17, mmWave standards for both WiFi and cellular have adopted a beamsweep-based procedure due to its simplicity2,3,4,5,30. In this work, we will also follow the beamsweep paradigm and our analysis and results are derived based on this assumption.

2.1 Effect of Beamforming on the Beamspace

Under the beamsweep assumption, the effective size of the beamspace is a function of the beamforming scheme employed at the transmitter and the receiver. That is, the beamforming scheme dictates how many directions the receiver can inspect in a single observation. There are three beamforming schemes:

  • Digital beamforming In this architecture, each antenna element is connected to a radio frequency (RF) chain and a pair of data converters, analog-to-digital or digital-to-analog converters (ADC or DAC). The beamforming (or spatial filtering) is performed by the digital baseband processor. For wideband systems with a large number of antennas, this architecture can have high power consumption when high precision DACs and ADCs are used. One way to use digital beamformers at mmWave is to use DACs and ADCs with few bits of quantizer resolution. This is attractive for beamformed systems as having digital samples from every antenna element allows the receiver to inspect as many directions as the number of antenna elements with full directional gain using only one observation of the channel. This reduction of the number of needed observations becomes significant as the beamspace grows and can potentially reduce the total power consumed during the search procedure.

  • Analog beamforming To avoid the use of a large number of DACs and ADCs, analog beamformers perform beamforming (or spatial filtering) on the analog (in RF or intermediate frequency) signal using RF phase shifters and power combiners (or splitters). The use of just a pair of ADCs considerably reduces the power consumed by these front-ends and hence, these are considered a prime candidate for initial mmWave-based cellular devices. However, analog beamformers can point their beams only in one direction at a given time. This leads to potentially high delays when the beamspace is large.

  • Hybrid beamforming This scheme is a combination of the digital and analog beamforming. A part of the beamforming is performed by M analog RF beamformers. These beamformed signals are digitized and combined (or precoded) by the digital baseband circuitry. This allows the receiver to inspect M directions at each time instance. Now, at an extreme \(M=1\), where we have analog beamforming. At the other extreme, M equals the number of antenna elements, where the scheme is effectively digital beamforming. Choosing M trades off spatial multiplexing capabilities on the one hand and energy consumption on the other.

Consider a transmitter equipped with an antenna array comprised of \(N_{\rm{Tx}}\) antenna elements and a receiver equipped with \(N_{\rm{Rx}}\) elements. This means that the size of the beamspace is equal to the product of \(N_{\rm{Tx}}\) by \(N_{\rm{Rx}}\). If they both use analog beamforming, then they must visit all these directions at \(N_{\rm{Tx}} \times N_{\rm{Rx}}\) separate channel inspections. Hence, the effective size of the beamspace is,

$$\begin{aligned}&\text {analog:} \qquad L_{\rm{an}} = N_{\rm{Tx}} \times N_{\rm{Rx}}. \end{aligned}$$
(1)

Now, suppose that the transmitter still uses analog beamforming but the receiver uses digital. Then, the receiver can “listen” on all the \(N_{\rm{Tx}}\) directions simultaneously and hence in this effective case the size of the beamspace becomes

$$\begin{aligned}&\text {digital:} \qquad L_{\rm{di}} = N_{\rm{Tx}}. \quad \quad \quad \end{aligned}$$
(2)

Applying the same logic, for hybrid beamforming, it is easy to see that the effective size of the beamspace is

$$\begin{aligned}&\text {hybrid:} \qquad L_{\rm{hy}} = \frac{N_{\rm{Tx}} \times N_{\rm{Rx}}}{M} \end{aligned}$$
(3)

Other combinations of beamforming schemes on either the transmitter or the receiver yield an effective beamspace of various sizes. Adopting one beamforming scheme versus the other has a fundamental impact on the effective beamspace. This, in turn, affects the time needed for the receiver to determine the best direction of the incoming signal. Furthermore, during beamsweeping the receiver front-end is always on. A front-end architecture that can listen on one or a few directions at a time will hence need to be powered on for a longer period of time to measure all the possible beam pairs. This can mean a considerable increase in the effective power consumed by analog and hybrid beamformers. In the next sections, we quantify the impact of the size of the beamspace and the choice of the beamformers on the energy consumed by the beam discovery procedure.

Figure 2:
figure 2

Frame structure of 3GPP NR. A frame is divided into ten subframes, a subframe into eight slots, and a slot into 14 OFDM symbols. A single SS block takes up 4 OFDM symbols

2.2 The 3GPP NR Paradigm

As an illustrative example of a system using beamsweep, we will discuss the 3GPP NR physical layer standard for initial access (IA) in mmWave bands, 3GPP frequency range 2 (FR2). Our analysis and results in the upcoming sections will all have an NR system as an underlying assumption. We chose NR for two reasons. One, it is the standard defined for 5G cellular systems, so it is expected to be adopted by millions of devices in the coming years. Two, the NR standard offers a well-defined set of assumptions regarding the beam discovery process and time and frequency numerology. This allows us to evaluate the different beamforming schemes with respect to energy consumption within a widely accepted context.

We next present the beam discovery process which is one of the basic steps of NR initial access, i.e., the procedure of establishing a link-layer connection between a base station, referred to as NR 5G nodeB (gNB) in the NR standard, and a mobile device or user equipment (UE).

2.2.1 NR Beam Discovery

Initial access involves a set of message exchanges between the gNB and the user equipment, whereby the user equipment identifies a serving gNB and synchronizes to it. The user equipment learns the physical cell identity of the gNB, sends back its own ID, and finally attaches to the cell. In this work, we will focus only on the first part of NR initial access: beam and cell discovery on the user equipment side, since this is the most critical in terms of detection delay and energy consumption.

On the gNB side, all the \(N_{\rm{Tx}}\) directions are swept with a periodicity of \(T = 20~\text {ms}\). More specifically, every T and for an interval of duration \(T_{\rm{ssb}} = 5~\text {ms}\), the gNB transmits B blocks of four OFDM symbols in \(N_{\rm{TX}}\) directions. These transmissions during \(T_{\rm{ssb}}\) are called a synchronization signal (SS) burst and each block of the four OFDM symbols, an SS block. Each SS block is comprised of the primary synchronization signal (PSS), the secondary synchronization signal (SSS) and the physical broadcast channel (PBCH). The PSS and the SSS together make up the physical ID of the cell. The first one takes one of three possible values among \(\{ 0,1,2\}\) while the latter, one of \(\{ 0,1, \ldots , 335\}\), resulting in a total of 1008 unique cell IDs. Each of these signals takes up 127 sub-carriers in frequency, while an entire SS block, including PBCH, includes 240 sub-carriers. In Fig. 2, we depict the time–frequency resources occupied by the SS block within an NR frame, subframe and slot.

Figure 3 shows an example of several SS blocks within a single SS burst. For the mmWave bands, the NR standard provisions a total of \(B=32\) SS blocks during each SS burst. It is assumed that each SS block is transmitted simultaneously in two directions using hybrid beamforming. Hence, a mmWave gNB can support up to 64 non-overlapping directions. As an example, in Fig. 4, we depict a scenario of a gNB beamsweeping \(N_{\rm{Tx}} = 8\) directions using \(B = 4\) SS blocks per SS bursts.

Now, the user equipment is also sweeping \(N_{\rm{Rx}}\) directions in searching for SS burst, see Fig. 4. These signals are known to the user equipment. It has to determine which one of the possible three PSS sequences and 336 SSS sequences were sent. Since the structure of an SS block is also known, once the PSS is detected and the optimal direction is found, the user equipment will move on to detecting the SSS within the same SS block, which is the third OFDM symbol, as shown in Fig. 2. Thus, the most critical part of beamsweeping is the detection of the PSS which unlocks all the remaining steps of initial access.

The user equipment is assumed to use analog beamforming. Hence, unlike the gNB it can probe only one direction at a time. Since the gNB sends the SS blocks in two directions simultaneously, according to Sect. 2.1 the effective size of the beamspace around the gNB and the user equipment is \(L = \left( {N_{\rm{TX}} \times N_{\rm{RX}}}\right) /{2}\).

Figure 3:
figure 3

SS burst. Each burst lasts for \(5 \,\rm{ms}\) and is repeated every \(20\, \rm{ms}\). Each SS block within an SS burst is shown in a different color.

Figure 4:
figure 4

Example of beam discovery within four SS bursts. During each SS burst, the gNB covers \(N_{\rm{Tx}} = 8\) directions in 4 SS blocks. The user equipment, on the other hand, scans all the \(N_{\rm{Rx}} = 4\) directions in four SS bursts. Thus, the user equipment and the gNB together need four SS bursts to cover the beamspace of size \(L = 16\).

3 Signal and System Model

We will leverage the NR beam discovery process described earlier for our analysis and modeling, with a few simplifying assumptions. We consider a single cell of radius \(d_{r}\) with the gNB situated at the center transmitting the PSS signals periodically. The user equipment, through analog beamforming will sweep \(N_{\rm{Rx}}\) directions sequentially in search of the PSS signal.

Our first simplification is assuming analog beamforming at the gNB. Therefore, the user equipment and the gNB together will have to sample a beamspace of size \(L_{\rm{ang}} = N_{\rm{tx}} \times N_{\rm{rx}}\). Note that for user equipment at low SNR, i.e., those at the edge of the cell, it may be needed to cycle through the beamspace more than once. We will denote each cycle with \(k = 1, 2, \ldots ,K\).

Our next assumption is that similar to12, the dominant path between the gNB and the user equipment is a line of sight (LOS) path aligned with exactly one of the transmitter–RX directions/sectors \(\ell ^{*}\) in the beamspace. Although this assumption may seem unrealistic, since real channels are seldom comprised of a single path, it is only used to derive our detectors below. In our evaluation and simulations, we will test our detectors using a channel model derived from real measurements8.

With these assumptions in mind, we model the received post-analog–BF received complex signal \({\mathbf {y}}_{k\ell }\) at the user equipment during the \(\ell\)th sampling of the beamspace and kth sweeping cycle as:

$$\begin{aligned} {\mathbf {y}}_{k \ell } = {\mathbf {u}}_{k \ell }^{*}{\mathbf {H}}_{k \ell } {\mathbf {v}}_{k \ell }^{*} {\mathbf {x}}+ {\mathbf {n}}_{k \ell }, \quad {\mathbf {n}}_{k\ell } \sim {\mathcal {C}}N (0,\sigma _{k\ell }^2{\mathbf {I}}_D). \end{aligned}$$
(4)

The transmitted PSS signal \({\mathbf {x}}\in {{\mathbb {C}}}^D\) is in a signal space with \(D \approx T_{\rm{PSS}} \times W_{\rm{pss}}\) orthogonal degrees of freedom, where \(T_{\rm{PSS}}\) and \(W_{\rm{pss}}\) are the PSS duration in time and bandwidth occupied by the PSS signal, respectively. We assume the PSS signal to be unit norm, \(\Vert {\mathbf {x}}\Vert = 1\). The vectors \({\mathbf {u}}_{k \ell } \in {{\mathbb {C}}}^{N_{\rm{RX}}}\) and \({\mathbf {v}}_{k \ell } \in {{\mathbb {C}}}^{N_{\rm{TX}}}\) are, respectively, the user equipment and gNB side beamforming vectors along a transmitter–receiver direction \(\ell\) in the sectorized beamspace at the kth sweep cycle. They are assumed to be of fixed norm.

The MIMO channel \({\mathbf {H}}\) is assumed to be flat fading within the PSS bandwidth and constant over the PSS transmission time \(T_{\rm{PSS}}\). It is defined as

$$\begin{aligned} {\mathbf {H}}= h_{k \ell }{\mathbf {u}}_{\ell ^{*}} {\mathbf {v}}_{\ell ^{*}}^*, \end{aligned}$$
(5)

where \(h_{k \ell }\) is a small-scale fading coefficient. The vectors \({\mathbf {u}}_{\ell ^{*}}\) and \({\mathbf {v}}_{\ell ^{*}}\) are the spatial signatures of the user equipment and the gNB antenna arrays describing a single LOS path between them aligned with only one of the transmitter–RX directions in the beamspace. They are given as

$$\begin{aligned} {\mathbf {u}}_{\ell ^{*}}&= [1, e^{-j 2 \pi {\varDelta }\cos (\phi _{\ell ^{*}}^{\rm{Rx}} )}, \ldots , e^{-j 2 \pi {\varDelta }(N_{\rm{RX}} -1)\cos (\phi _{\ell ^{*}}^{\rm{Rx}} )} ]^{\rm{T}} \\ {\mathbf {v}}_{\ell ^{*}}&= [1, e^{-j 2 \pi {\varDelta }\cos (\phi _{\ell ^{*}}^{\rm{Tx}} )}, \ldots , e^{-j 2 \pi {\varDelta }(N_{\rm{TX}} -1)\cos (\phi _{\ell ^{*}}^{\rm{Tx}} )} ]^{\rm{T}}, \end{aligned}$$
(6)

where \(\phi _{\ell ^*}^{\rm{Rx}}\) and \(\phi _{\ell ^*}^{\rm{Tx}}\) are the angles of arrival and departure in the directions of \(\ell ^{*}\) and \({\varDelta }\) is the distance between elements of the antenna arrays measured in wavelengths. Finally, \({\mathbf {n}}\) in (4) is the i.i.d. complex additive white Gaussian noise vector with covariance \(\sigma _{k\ell }^2{\mathbf {I}}_D\) and \({\mathbf {I}}_D\) is the \(D\times D\) identity matrix.

The aim of the user equipment and the gNB is to mutually steer their beams along their spatial signatures defined in (6). That is, to apply beamforming vectors, \({\mathbf {u}}_{\ell }\) and \({\mathbf {v}}_{\ell }\), as close to \({\mathbf {u}}_{\ell ^{*}}\) and \({\mathbf {v}}_{\ell ^{*}}\) as possible. This, based on our models in Eqs. (4) and (5) and our fixed norm assumption over the beamforming vectors, leads to the maximization of the energy of the received signal vector \({\mathbf {y}}_{k \ell }\).

3.1 Signal Detection

Given the signal and transmission model above we now derive the PSS detector through hypothesis testing. We will follow the same line of analysis given in12. Since we have assumed that the dominant path from the gNB to the user equipment is along a single transmitter–receiver beamspace direction \(\ell ^{*}\) (here the asterisk denotes the true or best direction), for each probed direction in each sweeping cycle our null hypothesis, \(H_0\), is that the signal is not present. Conversely, our alternative hypothesis, \(H_1\), is that the signal is present, i.e., the user equipment is probing the direction in which the gNB is transmitting. Hence, Eq. (4) under \(H_0\) and \(H_1\) becomes

$$\begin{aligned} H_0: ~ {\mathbf {y}}_{k \ell }&= {\mathbf {n}}_{k \ell }&\text{(signal } \text{ not } \text{ present), } \\ H_1: ~ {\mathbf {y}}_{k \ell }&= \alpha _{k \ell ^{*}} {\mathbf {x}}+ {\mathbf {n}}_{k \ell ^{*}}&\text{(signal } \text{ present } \text{ at } \text{ correct } \text{ beamspace } \text{ direction), } \end{aligned}$$
(7)

where the scalar channel coefficient \(\alpha _{k \ell ^{*}}\) is the result of applying the beamforming vectors along the spacial signatures of the user equipment and the gNB and is given by

$$\begin{aligned} \alpha _{k \ell ^{*}} = {\mathbf {u}}_{\ell ^{*}}^* {\mathbf {H}}_{k \ell ^{*}} {\mathbf {v}}_{\ell ^{*}} = h_{k \ell ^{*}} {\mathbf {u}}_{\ell ^{*}}^* {\mathbf {u}}_{\ell ^{*}} {\mathbf {v}}_{\ell ^{*}}^* {\mathbf {v}}_{\ell ^{*}}. \end{aligned}$$
(8)

We define the probability density of the received signal \({\mathbf {y}}\) under the two hypotheses as \(p({\mathbf {y}}| H_i, {\varvec{\alpha }}, {\sigma ^2}, \ell ^{*} )\) for \(i = \{0,1\}\), where the set \({\mathbf {y}}\) contains all the observed signals \({\mathbf {y}}_{k \ell }\) in K beam sweeps. The model contains several unknown parameters, namely, \({\varvec{\alpha }}\) the set of all the channel coefficients \(\alpha _{k \ell ^{*}}\), and \({\sigma }^2\) the set of all noise power levels.

Due to these unknown parameters, we will rely on the widely used generalized likelihood ratio test (GLRT) method49. The GLRT takes the likelihood distribution under each hypothesis maximized with respect to the unknown parameters and compares their ratio to a threshold. Alternatively, we can use the minimum log-likelihoods of the signal distribution under each hypothesis. These are given as

$${\varLambda }_0 := \min _{{\varvec{\sigma }}^2} -\ln p({\mathbf {y}}|H_0, {\varvec{\sigma }}^2)$$
(9a)
$${\varLambda }_1 := \min _{{\varvec{\sigma }}^2,{\varvec{\alpha }}} - \ln p({\mathbf {y}}|H_1, {\varvec{\sigma }}^2,{\varvec{\alpha }}).$$
(9b)

We then use the test

$${\varLambda }:= {\varLambda }_0 - {\varLambda }_1 \underset{H_0}{\overset{H_1}{\gtrless }} \; \theta$$
(10)

for a threshold \(\theta\). We observe that under the assumption that the noise vectors are independent in different measurements,

$$\ln p({\mathbf {y}}|H_0, {\varvec{\sigma }}^2,\ell ^{*}) = \sum _{k=1}^K \sum _{\ell =1}^L \ln p({\mathbf {y}}_{k\ell }|\sigma _{k \ell }^2,\ell ^{*})$$
(11a)
$$\ln p({\mathbf {y}}|H_1, {\varvec{\sigma }}^2,{\varvec{\alpha }},\ell ^{*}) = \sum _{k=1}^K \sum _{\ell =1}^L \ln p({\mathbf {y}}_{k\ell }|\sigma _{k \ell }^2,\alpha _{k \ell ^{*}},\ell ^{*}).$$
(11b)

Therefore, the negative log-likelihoods in Eq. (9) can be re-written as

$${\varLambda }_0 = \sum _{k=1}^K \sum _{\ell =1}^L {\varLambda }^{k\ell }_0(\ell ^{*}).$$
(12a)
$${\varLambda }_1 = \sum _{k=1}^K \sum _{\ell =1}^L {\varLambda }^{k\ell }_1(\ell ^{*}),$$
(12b)

where \({\varLambda }^{k\ell }_0(\ell ^{*})\) and \({\varLambda }^{k\ell }_1(\ell ^{*})\) are the minimum negative log-likelihoods in each measurement and given as:

$${\varLambda }^{k\ell }_0(\ell ^{*}) = \min _{\sigma ^2} -\ln p({\mathbf {y}}_{k\ell }|\sigma _{k \ell }^2,\ell ^{*}).$$
(13a)
$${\varLambda }^{k\ell }_1(\ell ^{*}) = \min _{\sigma ^2,\alpha _{k \ell }} -\ln p({\mathbf {y}}_{k\ell }|\sigma _{k \ell }^2,\alpha _{k \ell }, \ell ^{*}).$$
(13b)

Since the received signal in our model, given by Eq. (7), is Gaussian conditional on the parameters, hence

$${ -\ln p({\mathbf {y}}_{k \ell }|\sigma _{k \ell }^2,\ell ^{*}) }= \frac{1}{\sigma _{k \ell }^2} \Vert {\mathbf {y}}_{k \ell }\Vert ^2 + D\ln (\pi \sigma _{k\ell }^2)$$
(14)

under the \(H_0\) hypothesis, and

$$\begin{aligned} { -\ln p({\mathbf {y}}_{k \ell }|\sigma _{k \ell }^2,\alpha _{k \ell ^{*}},\ell ^{*}) }&= \frac{1}{\sigma _{k \ell }^2} \Vert {\mathbf {y}}_{k \ell } - \alpha _{k \ell ^{*}} {\mathbf {x}}\Vert ^2 + D\ln (\pi \sigma _{k\ell }^2) \end{aligned}$$
(15)

under the \(H_1\) hypothesis. Now, we use the above Eqs. (14) and (15) to obtain estimates of the unknown variables which are their minimizers.

We first take \({\varLambda }^{k\ell }_1(\ell ^{*})\). The channel coefficient \(\alpha _{k \ell ^{*}}\) that minimizes the log negative likelihood when the signal is present is given as

$$\begin{aligned} {\widehat{\alpha }}_{k \ell ^{*}} = \frac{{\mathbf {x}}^*{\mathbf {y}}_{k\ell }}{\Vert {\mathbf {x}}\Vert ^2}. \end{aligned}$$
(16)

Next we move on to minimize over the variance \(\sigma ^2\), which occurs at

$$\begin{aligned} {\widehat{\sigma }}_{k \ell }^2 = \frac{1}{D} \left( \Vert {\mathbf {y}}_{k\ell }\Vert ^2 - \frac{|{\mathbf {x}}^*{\mathbf {y}}_{k\ell }|^2}{\Vert {\mathbf {x}}\Vert ^2} \right) . \end{aligned}$$
(17)

When the signal is not present, i.e., \(\ell \ne \ell ^{*}\), we minimize \({\varLambda }^{k\ell }_0(\ell ^{*})\) over the variance \(\sigma ^2\). The minimum is obtained at

$${\widehat{\sigma }}_{k\ell }^2 = \frac{\Vert {\mathbf {y}}_{k\ell }\Vert ^2}{D}.$$
(18)

Substituting Eqs. (16), (17) and (18), in Eq. (13) we obtain

$${\varLambda }^{k\ell }_0(\ell ^{*}) = D \ln \left( \frac{ \pi e}{D} \Vert {\mathbf {y}}_{k\ell }\Vert ^2\right) ,$$
(19a)
$${\varLambda }^{k\ell }_1(\ell ^{*}) = D \ln \left( \frac{ \pi e}{D} \left( \Vert {\mathbf {y}}_{k \ell } \Vert ^2 - \frac{| {\mathbf {x}}^* {\mathbf {y}}_{k \ell }|^2}{\Vert {\mathbf {x}}\Vert ^2}\right) \right) .$$
(19b)

Combining (19) and from the test (10), and Eq. (12) we have the following log-likelihood difference

$$\begin{aligned} {\varLambda }= \sum _k^K \sum _{\ell }^L - D \ln \left( 1 - \rho _{k \ell } \right) , \end{aligned}$$
(20)

where the \(\rho _{k \ell }\) is the normalized energy of the correlation of the known signal with the received signal, obtained as

$$\begin{aligned} \rho _{k \ell } = \frac{|{\mathbf {x}}^* {\mathbf {y}}_{k \ell }|^2}{\Vert {\mathbf {x}}\Vert ^2 \Vert {\mathbf {y}}_{k \ell }\Vert ^2}. \end{aligned}$$

This, not surprisingly, means that while the user equipment probes each direction in K beamsweeps, it should correlate the incoming signal with the local replica of the known signal. Then, since Eq.  (20) is an increasing function of \(\rho _{k \ell }\), and due to our assumption that the signal is present only in one direction within the beamspace, the best direction \({\widehat{\ell }}^{*}\) is the one in which the correlation was the highest, that is,

$$\begin{aligned} {\widehat{\ell }}^{*} = \mathop {\rm{arg\,min}}\limits _{\ell = 1, \ldots , L} \sum _k^K\sum _\ell ^L \ln (1 - \rho _{k \ell }). \end{aligned}$$
(21)

While the above detector was derived for analog beamforming under favorable beamspace assumptions, it can also be used for both hybrid and digital beamforming. The difference is that with these beamforming architectures, the receiver has access to \(N_{\rm{RF}} >1\) RF chains and can essentially operate as \(N_{\rm{RF}}\) parallel analog systems. In the digital beamforming case \(N_{\rm{RF}}\) is equal to \(N_{\rm{Rx}}\). Note that with a fully digital detector, the receiver can test all the angles, including those not perfectly aligned in the beamspace. However, for the sake of simplicity and better comparison with analog beamforming, we will not consider a more powerful detector but assume that in digital beamforming too the arrival angles are aligned with the beamspace directions.

3.2 RF Architectures and Power Consumption

While standards documents on the implementation of mmWave cellular systems assume user equipments with analog beamforming and only a single RF chain, we believe that there is much value in using fully digital beamforming. To make our case more concrete in terms of energy consumption, we give here an overview of the receiver’s mmWave front-end. Note that the part of the mmWave front-end before the variable gain amplifier (VGA) and the ADCs (see Figs. 567) is the radio frequency front-end (RFFE). We model the power consumption of analog, hybrid, and fully digital and fully digital with low-resolution quantization front tends.

3.2.1 mmWave Front-Ends

Analog front-end Consider the mmWave receiver’s analog front-end shown in Fig. 5. It is comprised of \(N_{\rm{RX}}\) low noise amplifiers (LNA)s, \(N_{\rm{RX}}\) phase shifter (PS), a mixer, a combiner and a pair of analog-to-digital converter (ADC). The D.C. power consumption of each LNA, \(P_{\rm{DC}}^{\rm{LNA}}\) is a function of its gain \(G_{\rm{LNA}}^{\rm{PS}}\), the noise figure \(N_{\rm{LNA}}\) and the figure of merit (FoM) and is given as7,21,46

$$\begin{aligned} P_{\rm{LNA}}^{\rm{DC}} = \frac{G_{\rm{LNA}}^{\rm{PS}}}{\text {FoM} (N_{\rm{LNA}} - 1)}, \end{aligned}$$
(22)

in linear scale.

Figure 5:
figure 5

Analog beamforming: receiver front-end architecture.

Figure 6:
figure 6

Digital beamforming: receiver front-end architecture.

The PS and mixer are considered to be passive elements. While they do not burn any power, they do, however, introduce insertion losses (IL), which need to be compensated for. For example, the LNA gain above needs to be high enough to offset the loss introduced by the phase shifter \(\rm{{IL}}_{\rm{PS}}\). That is, if an LNA in a circuit without any PS had a gain \(G_{\rm{LNA}}\), then \(G_{\rm{LNA}}^{\rm{PS}} = G_{\rm{LNA}} + \rm{{IL}}_{\rm{PS}}\).

For the ADCs, power consumption is a function of the sampling frequency \(F_{\rm{s}}\), ADC’s figure of merit and the resolution in bits:

$$\begin{aligned} P_{\rm{ADC}} = \text {FoM} \times F_{\rm{s}} \times 2^q, \end{aligned}$$
(23)

where q is the number of resolution bits. To keep the variable incoming signal within the ADCs dynamic range while keeping a constant baseband power \(P_{\rm{BB}}^{\rm{out}}\) it is necessary to apply gain control to the input of the ADC. This is performed by the variable gain amplifier (VGA) with a gain range between 0 and \(G_{\rm{VGA}}^{\rm{max}}\). For the analog beamforming case this maximum gain is calculated as:

$$\begin{aligned} G_{\rm{VGA}}^{\rm{max}} = P_{\rm{BB}}^{\rm{out}} - 10\log (N_{\rm{RX}}) + \rm{{IL}}_{\rm{mix}} - (G_{\rm{LNA}}^{\rm{PS}} - \rm{{IL}}_{\rm{PS}}) - P_{\rm{RX}}(d = d_{\rm{max}}), \end{aligned}$$
(24)

where \(\rm{{IL}}_{\rm{mix}}\) is the insertion loss introduced by the mixer and \(P_{\rm{RX}}(d = d_{\rm{max}})\) is the maximum distance between the transmitter and the receive, i.e., the cell edge.

Given the above gain, the D.C. power draw of the VGA is given as

$$\begin{aligned} P_{\rm{VGA}}^{\rm{DC}} = \frac{G_{\rm{VGA}}^{\rm{max}} F_{\rm{BW}}}{\text {FoM} A_{\rm{chip}}}, \end{aligned}$$
(25)

where FoM is the FoM of the VGA, \(F_{\rm{BW}}\) the bandwidth of the VGA in \(\rm{GHz}\) and \(A_{\rm{chip}}\) is the active area in \({\hbox {mm}}^2\). For our analysis we assume that \(F_{\rm{BW}} = F_{\rm{s}}\), i.e., the bandwidth of the VGA is equal to the sampling rate of the ADC.

Based on the above, the total power consumption of an analog mmWave front-end can be calculated as

$$\begin{aligned} P_{\rm{an}} = N_{\rm{RX}} P_{\rm{LNA}}^{\rm{DC}} + P_{\rm{LO}} + P_{\rm{ADC}} + P_{\rm{VGA}}^{\rm{DC}}, \end{aligned}$$
(26)

where \(P_{\rm{LO}}\) is the power draw of the local oscillator.

Digital front-end Consider the digital front-end in Fig. 6. Since we have introduced most of the components in the analog front-end case, computing the power for the digital case is much easier. Notice in the figure, that since beamforming is done in baseband there is no need for phase shifters. Neither do the LNA gains need to compensate for their insertion losses. Therefore, we will use \(G_{\rm{LNA}}\) instead of \(G_{\rm{LNA}}^{\rm{PS}}\). However, the number of the ADCs is increased to \(N_{\rm{RX}}\), so is the number of VGAs, mixers and local oscillators. Thus, the DC power drawn by the LNA and VGA is given as

$$\begin{aligned} P_{\rm{LNA}}^{\rm{DC}} = \frac{G_{\rm{LNA}}}{\text {FoM} (N_{\rm{LNA}} - 1)}, \end{aligned}$$
(27)

and

$$\begin{aligned} G_{{\rm{VGA}}_{\rm{dig}}}^{\rm{max}} = P_{\rm{BB}}^{\rm{out}} - 10\log (N_{\rm{RX}}) + \rm{{IL}}_{\rm{mix}} - (G_{\rm{LNA}}) - P_{\rm{RX}}(d = d_{\rm{max}}). \end{aligned}$$
(28)

With these, the total power consumed by a fully digital mmWave front-end is

$$\begin{aligned} P_{\rm{dig}} = N_{\rm{RX}} (P_{\rm{LNA}}^{\rm{DC}} +P_{\rm{LO}} + P_{\rm{ADC}} + P_{\rm{{VGA}}_{\rm{dig}}}^{\rm{DC}}). \end{aligned}$$
(29)

Digital front-end with low-resolution ADCs The front-end here remains the same as in Fig. 6. The only thing that changes is the resolution q of the ADCs. Notice in Eq. (23), that the power consumption is exponential in q. Since the power consumption of a fully digital mmWave front-end grows linearly in the number of antenna elements \(N_{\rm{RX}}\) as given in Eq. (29), reducing the resolution is the only meaningful way of bringing the power draw of a fully digital front-end. In our calculations below and in Sect. 4 we will show how reducing the quantizer resolution affects both the power drawn by the front-end and the energy consumption during the beam discovery process.

Figure 7:
figure 7

Sub-array hybrid beamforming: receiver front-end architecture.

Hybrid front-end For the sake of completeness, we also present the mmWave hybrid front-end in Fig. 7. There are a few options in designing a hybrid front-end. Here, we show the “industry-standard” sub-array architecture where the antenna array is divided by the number of supported digital streams. The depicted front-end can support \(M = 2\) digital streams each connected to sub-arrays of size \(N_{\rm{RX}}/2\). Notice that this circuit has elements from both the fully digital and analog front-ends: \(N_{\rm{RX}}\) phase shifters and \(M = 2\) pairs of ADCs. Thus, the total power consumption of the hybrid front-end is

$$\begin{aligned} P_{\rm{hy}} = N_{\rm{RX}} P_{\rm{LNA}}^{\rm{DC}} + M (P_{\rm{LO}} + P_{\rm{ADC}} + P_{\rm{VGA}}^{\rm{DC}}). \end{aligned}$$
(30)
Table 1: mmWave receiver’s power consumption (mW).

3.2.2 Front-End Power Consumption

To give a better picture of each front-end’s power consumption, we will now give numerical examples. We will use power draws and losses of individual components as reported in studies in the literature and give a total number for each front-end when put together.

For an LNA with gain \(G_{\rm{LNA}} = 10~\text {dB}\) with a noise figure of \(N_{\rm{LNA}} = 3~\text {dB}\) and \(\text {FoM} = 6.5\, {\text {mW}}^{-1}\) we have used the \(90~\text {nm}\) CMOS LNA reported in7,21,46. We have assumed a PS insertion loss \(\rm{{IL}}_{PS} = 10~\text {dB}\) which the analog front-end LNA gain \(G_{\rm{LNA}}^{\rm{PS}}\) needs to compensate for.

For the ADC, we have used a 4-bit Flash-based ADC reported in36 with an \(\text {FoM}\) of \(67.6 \text {fJ/conversion step}\). While the 4-bit ADC falls into the category of a low-resolution ADC, we use the same architecture and \(\text {FoM}\) for higher quantization of \(q=10\) bits for better comparison. For the VGA, we consider the \(90 \text {nm}\) CMOS reported in50. The VGA’s FoM is 5280 for an active area of \({\text {0.01mm}}^2\). Last, we have assumed an oscillator power draw of \(10 ~ \text {dBm}\), a mixer insertion loss of \(\rm{{IL}}_{\rm{mix}}\) of \(6 ~ \text {dBm}\) and baseband power \(P_{\rm{BB}}^{\rm{out}}\) of \(10 ~ \text {dBm}\).

Using the above device properties, we, furthermore, consider an operating bandwidth (BW) of \(400~\text {MHz}\), the maximum provisioned bandwidth by 3GPP NR and a maximum distance \(d_{\rm{max}}\) of 100 meters. We then use the channel model8 to calculate the pathloss needed for computing the received power \(P_{\rm{RX}}\) in Eqs. (24) and (28). Specifically, for an EIRP Tx power of \(43 ~ \text {dBm}\) the received power is \(-87 ~ \text {dBm}\) at \(d_{\rm{max}}\) where the channel is expected to be non-line of sight (NLOS). Finally, we assume a sampling rate, \(F_{\rm{s}} = 1~\hbox {GHz}\) and a receiver array size of \(N_{\rm{Rx}} = 16\). While we assume a system bandwidth of \(400~\rm{MHz}\), however, according to the NR specifications OFDM sampling rate is set at \(491~\rm{MHz}\). Furthermore, we consider an oversampling rate at \(2 \times\) the system bandwidth to avoid aliasing. Hence, the \(1~\rm{GHz}\) sampling rate assumption.

By plugging the above numbers into (22)–(30), we obtain the front-end power consumption presented in Table 1. With the term RFFE we denote the radio frequency front-end, which is the pre-VGA part of the mmWave front-end, i.e., the LNAs and the LO. We observe a couple of points:

  • The power consumption of a high-resolution digital front-end is almost six times the consumption of the analog front-end. However, it is below a factor of \(N_{\rm{RX}}\). This is because while there are \(N_{\rm{RX}}\) ADCs in a digital front-end, there are also \(N_{\rm{RX}}\) phase shifters in an analog one, the insertion loss of which must be compensated by the LNAs. We will see in the next section how the gap between analog and digital closes when we take into account the beam discovery delay.

  • Reducing the ADC resolution from \(q = 10\) to \(q = 4\) bits has a dramatic effect on the power consumption. It is, in fact, even below that of the analog front-end. However, this reduction in resolution adds a distortion to the signal which will reduce the effective SNR. This reduction has a small effect in the low to median SNR regimes20,21.

4 Evaluation and Simulation

We now evaluate the performance of our correlation-based detector for receivers performing analog, fully digital, and low-resolution fully digital beamforming. In a series of steps, we first define the simulation setup and parameters, assess the effect of lowering the quantization resolution in fully digital front-ends, and finally illustrate the interplay between the beamspace and the energy consumption of each the front-end architectures presented in Sect. 3.

Table 2: SNR distribution parameters.

4.1 Channel Model and User SNR Distribution

We consider a single cell of radius \(d_r = 100~\text {m}\), with the mmWave base station, or gNB, situated at the center. We then randomly drop user equipments in this cell and compute their pathloss to the gNB according to the widely used model in8. This urban multipath model is based on real measurements in \(28~\text {GHz}\) bands performed in New York City10,32,43,55. While our detector was derived for a single-path LOS channel perfectly aligned with one of the beamspace directions, we will use this multipath channel in our evaluation and simulations. According to this model, links between the user equipment and the gNB are in LOS or NLOS based on an exponential probability distribution parametrized by the distance separating the two, d. Close-by user equipments have a high probability of being in LOS while as the distance grows the probability of being in NLOS increases. The resulting omni-directional pathloss, i.e., before applying beamforming, in \(\text {dB}\) is computed as

$$\begin{aligned} P_{L} = \mu + 10 \nu \log _{10}(d) + \zeta [\text {dB}], \quad \zeta \sim {\mathcal {N}}(0, \xi ^2), \end{aligned}$$
(31)

where \(\mu\), \(\nu\) and \(\xi ^2\) are parameters defined by whether the link is LOS or NLOS—see Table 2.

We then use this pathloss to derive the omni-directional downlink (DL) SNR at the user equipment, \({\rm{SNR}}_{\rm{omni}}\), as

$$\begin{aligned} {\rm{SNR}}_{\rm{omni}} = \frac{P_{\text {RX}}}{N_0 W_{\rm{sys}}}, \end{aligned}$$
(32)

where \(P_{\rm{RX}}\) is the received omni-directional power resulted from subtracting, in dB scale, the pathloss in (31) from the transmitted power \(P_{\rm{TX}}\). Note that all the terms in the equation itself have been converted in a linear scale. The remaining parameters above are the noise power density plus noise figure, \(N_0\), and the system bandwidth \(W_{\rm{sys}}\).

Next, we derive the user equipments SNR distribution through 10,000 random drops in our cell and divide them into three regimes: cell edge user equipments, i.e., the first percentile of the CDF, median users and above, and those in between these regimes. However, in the sequel, we characterize the energy consumed by the analog and digital beamformers during beam discovery for the cell edge and the median users. Figure 8 shows the CDF of the downlink omni-directional SNR distribution, based on Eq. (32), converted in dB scale. A summary of the parameters and their values used to derive this distribution is presented in Table 2.

Figure 8:
figure 8

Downlink SNR distribution of 10,000 user equipment drops. Each point corresponds to a value calculated using Eq. (32) and then converted in dB scale. The solid red line and the dashed gray line denote the cell edge and the median user snr regimes, respectively.

4.2 Low-Resolution Fully Digital Beamforming

In Sect. 3.2, we presented the power consumption of four mmWave front-ends implementing analog, hybrid, and fully digital beamforming. We also showed that the power consumption of the fully digital front-end could be dramatically reduced by employing low-resolution ADCs.

However, reducing the resolution will degrade the quality of the received signal. We quantify this degradation using the method presented in11,16,24. According to24, the effect of low-resolution quantization, i.e., analog-to-digital conversion, is modeled as a reverse gain \((1-\gamma )\) multiplied with the received sample if it had gone through an infinite-resolution quantizer and an additive Gaussian noise term. That is if the high quantization sample is given as

$$\begin{aligned} y[n] = x[n] + w[n], \end{aligned}$$

where x[n] are the transmitted signal sample and w[n] AWG noise with some variance \({\mathbb {E}}|{w[n]|^2 = \sigma _n^2}\), then the imperfectly quantized received sample is given as

$$\begin{aligned} y_q[n] = (1 - \gamma )x[n] + (1 - \gamma ) w[n] + v[n]. \end{aligned}$$
(33)

The last term v[n] is a zero-mean Gausian random variable with the following variance11,24

$$\begin{aligned} \sigma _{\rm{q}}^2 ={\mathbb {E}}{|v [n]|^2} = \gamma (1-\gamma ) ( {\mathbb {E}}{|x[n]|^2} + \sigma _n^2). \end{aligned}$$
(34)

Now, using Eq. (33) the effective SNR \({\rm{SNR}}_{\rm{eff}}\) after quantization is obtained as

$$\begin{aligned} {\rm{SNR}}_{\rm{eff}} = \frac{(1 - \gamma ) \rm{SNR_{\infty }}}{ 1+ \gamma \rm{SNR_{\infty }}}, \end{aligned}$$
(35)

where \(\rm{SNR_{\infty }}\) is the SNR value of an infinite-resolution quantizer. The value of \(\gamma\) depends on the quantizer’s design, the input distribution, and the quantization resolution bits q. As an example, for a Gaussian input the value of \(\gamma\) is only 0.01.

Thus, we see that low-quantization fully digital beamforming comes with a penalty in the SNR, an observation corroborated by the theoretic work in35,45. However, this penalty can be negligible depending on the quantization bits q, and the SNR regime, as shown in20,21, where the same model was used for OFDM inputs. The effective SNR can be much lower, though, for very high values of \({{\rm{SNR}}_{\infty }}\) at very low quantization bits, say one or two. However, in cellular systems, the overwhelming majority of user equipments are at low to medium SNR regimes where the distance between the effective SNR and \({{\rm{SNR}}_{\infty }}\) is below \(3~\rm{dB}\). Furthermore, at a quantization of \(q = 4\) bits as we consider for our low-resolution fully digital front-end, even at \({{\rm{SNR}}_{\infty }}\) values as high as \(20~\rm{dB}\), this distance is still below \(3 ~\rm{dB}\)20,21.

Table 3: System parameters.

4.3 System Parameters

We will compare the mmWave front-ends presented in Sect. 3, within the context of 3GPP NR’s PSS discovery. We assume a PSS signal that is confined within one OFDM symbol and 127 sub-carriers. For a sub-carrier spacing of \({\varDelta }_f = 120 \rm{kHz}\), the PSS signal duration and bandwidth are \(T_{\rm{PSS}} = 8.91 \upmu \rm{s}\) and \(W_{\rm{PSS}} = 15.24 \, \rm{MHz}\), respectively. We will assume that the channel is almost flat within the duration and bandwidth of the signal.

In this work, we assume that during an SS burst, there is no multiplexing with data symbols. That is, all the power of the base station is put on the signals in the SS block. Thus, the omni-directional PSS signal SNR, \({\rm{SNR}}_{\rm{PSS}}\), is given as

$$\begin{aligned} {\rm{SNR}}_{\rm{PSS}} = \frac{P_{\rm{RX}} T_{\rm{PSS}}}{N_0}. \end{aligned}$$
(36)

The above equation is related to the omni-directional in Eq. (32) as11:

$$\begin{aligned} {\rm{SNR}}_{\rm{PSS}} = W_{\rm{sys}} \times T_{\rm{PSS}} \times {\rm{SNR}}_{\rm{omni}} . \end{aligned}$$

Both the user equipment and the gNB are equipped with uniform planar arrays that allow them to beamsteer in both azimuth (az) and elevation (el). Note that in the sequel, we denote the available directions at the gNB and the user equipment with \(N_D^{\rm{Tx}}\) and \(N_D^{\rm{Rx}}\), respectively. These are the products of \(N_{\rm{Tx}}\) and \(N_{\rm{Rx}}\) in each dimension of the 2D arrays. Hence, the beamspace size L is equal to \(N_D^{\rm{Tx}} \times N_D^{\rm{Rx}}\). The gNB is assumed to perform only hybrid beamforming with \(M = 2\) directions at a time, while the user equipment may employ one of analog, fully digital, and low-resolution fully digital beamforming. Note that we keep the transmission power of the gNB fixed at \(30~~ \rm{dBm}\) (Table 2). Therefore, transmitting the PSS in \(M = 2\) directions implies that the sum of the power transmitted in each direction must be \(30~ \rm{dBm}\). We will assume equal power transmission. Thus, the signal in each direction will carry half of the total \(30~ \rm{dBm}\) power. The resulting beamforming gains will be applied to the user equipments’ omni-directional PSS signal SNR given by Eq. (36) or to its corresponding effective SNR in case of the low-resolution digital front-end.

Also, we assume a beamspace discovery process, as described in Sect. 2.2.1, where the gNB transmits the PSS directionally in its \(N_D^{\rm{Tx}}\) direction within an SS burst of duration \(5 \,\rm{ms}\) every \(T = 20~\rm{ms}\). The user equipment also sweeps its \(N_D^{\rm{Rx}}\) directions in search of one of the \(N_{\rm{PSS}} = 3\) PSS sequences. Finally we compute the threshold \(\theta\) in test (10). We will do so through the false alarm probability \(P_{\rm{FA}}\). Suppose we would like to maintain a constant target false alarm rate of \(R_{\rm{FA}}\) during each searching period. Then, the false alarm probability and the false alarm rate are connected through

$$\begin{aligned} P_{\rm{FA}} = \frac{R_{\rm{FA}}}{N_{\rm{PSS}} N_{\rm{dly}} N_{\rm{FO}}}. \end{aligned}$$
(37)

Now, \(N_{\rm{dly}}\) is the number of delay hypotheses. That is at which delay tau the signal was present. Since the PSS transmission is periodic, the number of the delay hypotheses is bounded, and it is the number of the samples within each SS burst period. If the signal was present at a delay \(\tau _1\), then it is also present at delay \(T + \tau _1\), where T is the SS burst period.

To compensate for two other sources of uncertainty, we incorporate them into the \(P_{\rm{FA}}\) calculation. These are the number of possible PSS sequences \(N_{\rm{PSS}}\) and the number of frequency offsets hypotheses \(N_{\rm{FO}}\). There are \(N_{\rm{PSS}}\) sequences that the user equipment may mistake one with the other two and are factored in the false alarm calculation. Next, we approximate \(N_{\rm{FO}}\) as follows. First, we assume a local oscillator error of \(\pm 10\) parts per million (ppm) and a Doppler shift for a velocity up to \(30~\rm{km/h}\) at \(F_{\rm{c}} = 28~\rm{GHz}\). The LO error and the Doppler shift will define a range \(\pm {\varDelta }_{f\rm{{max}}}\). We then discretize this range so that in each interval the channel does not rotate more than \(\pi /4\) within the duration of one OFDM symbol. The number of frequency offset hypotheses is then the number of obtained intervals.

The threshold \(\theta\), therefore, is the value of the normalized correlation \(\rho _{k\ell }\) which corresponds to the target \(P_{\rm{FA}}\). The detector will consider everything above this threshold as received signal and anything below as noise.

All the parameters described above are tabulated in Table 3.

4.4 Comparison of Different Beamforming Architectures

Now that we have defined the channel model, system and signal parameters, and described the beamspace and beam discovery procedure, we can start the comparison between the mmWave front-ends presented in Sect. 3.2. We first look into the beam discovery delays with beamsweeping for the beamformers and then tie these delays to the power consumptions reported in Sect. 3.2.2.

4.4.1 Beamsweeping Delay

We estimate the initial discovery delay as the time needed to go through all the beams in the beamspace times the total number beamsweeps K necessary to determine the best direction for a target miss detection probability\(P_{\rm{MD}}^{\rm{tgt}}\). Hence, K is a function of \(P_{\rm{MD}}^{\rm{tgt}}\). Note that K is also a function of the SNR. A low SNR user at the cell edge will have to sweep the beamspace and visit the same correct direction multiple times to accumulate enough energy required for the signal detection.

In Fig. 9 we present the miss detection probability \(P_{\rm{MD}}\) vs \({\rm{SNR}}_{\rm{PSS}}\), i.e., the PSS signal SNR, . There, we show the number of beamsweeps necessary to discover the correct path to the gNB for a \(P_{\rm{MD}}^{\rm{tgt}} = 1\%\). A mobile user at the cell edge needs \(K = 3\) beamsweeps to get the right direction, while the rest of the users can detect the signal in one beamsweep. Take as an example the case of the analog receiver with \(N_D^{\rm{Rx}} = 16\) directions, while the transmitter has \(N_D^{\rm{Tx}} = 64\) directions. The size L of the beamspace is then \(16 \times 64 = 1024\). Then, the delay in initial discovery \(T_{\rm{delay}}(P_{\rm{MD}}^{\rm{tgt}})\), for such a user equipment to successfully discover the true direction is

$$\begin{aligned} \frac{(K(P_{\rm{MD}}^{\rm{tgt}}) - 1 )\times L \times T}{N_{\rm{D}}^{\rm{Tx}}} + T_{\rm{ssb}}\le & {} T_{\rm{delay}}(P_{\rm{MD}}^{\rm{tgt}}) \le \frac{K(P_{\rm{MD}}^{\rm{tgt}})\times L \times T}{N_{\rm{D}}^{\rm{Tx}}} \nonumber \\= & {} \frac{20480 K(P_{\rm{MD}}^{\rm{tgt}}) }{64} \rm{ms}, \end{aligned}$$
(38)

where \(T = 20~\rm{ms}\) is the SS burst period and \(T_{\rm{ssb}} = 5~\rm{ms}\) the duration of an SS burst. The reason we divided the quantity above with \(N_D^{\rm{Tx}}\) is that within each SS burst period T, the gNB sweeps all its directions in one SS burst. Also, notice that we did not divide L by the number of the RF chains \(M =2\) at the transmitter. Since the NR procedure defines a maximum of \(B= 32\) SS blocks within one SS burst and the number of directions the gNB scans is \(N_{\rm{D}}^{\rm{Tx}} = 64\), the reduction in the effective beamspace due to hybrid beamforming is accounted for. Otherwise, dividing Eq. (38) by M would imply that the user equipment probes \(M = 2\) directions each time, which is not correct since it employs an analog beamformer.

From Eq. (38), we can see that even for a user that is not at the cell edge, i.e., \(K (P_{\rm{MD}}^{\rm{tgt}}) = 1\), the initial discovery delay can go as high as \(320\, \rm{ms}\), when the entire control-plane delay of 4G LTE must be below \(50\,\rm{ms}\)1. For the cell edge user, on the other hand, this delay is between \(640 ~\rm{ms}\) and almost a second. This highlights the large impact of the beamspace size on the initial discovery delay. Now, during initial discovery, the user equipment must be constantly “on” searching for the correct beamspace. This amount of delay, as we will show next, has a large impact on the power consumption of the analog beamformer.

Figure 9:
figure 9

Analog receiver: miss detection probability \(P_{\rm{MD}}\) vs \({\rm{SNR}}_{\rm{PSS}}\) in dB, as given by Eq. (36). \(K = 1,2,3\) and target miss detection probability \(P_{\rm{MD}}^{\rm{tgt}} = 1\%\).

Table 4: Upper bounds of initial discovery delays for different mmWave front-ends.

In Table 4 following the same process and using Eq. (38), we give the upper bounds of the initial discovery delay for the other two kinds of front-ends, fully digital and low-resolution fully digital. All the numbers are for a user equipment array size of 16 elements. Notice that the SNR penalty due to low quantization is negligible in both the cell edge and the rest of the area. The large gain compared to analog was the result of shrinking the beamspace size L to only \(N_{\rm{D}}^{\rm{Tx}}\) which is all covered in one SS burst.

Figure 10:
figure 10

Energy consumption of the mmWave front-ends: analog, high-resolution digital and low-resolution digital. Number of antenna elements at the user equipment \(N_{\rm{Rx}} = 16\). Cell edge user equipment in blue, median user equipment in red.

4.4.2 Beamsweeping Energy Consumption

We have so far characterized the beam discovery delay of the analog, the high-resolution digital, and the low-resolution digital beamformers. The energy consumption of a front-end is tightly connected to the amount of time it needs to operate. Since the user equipment is assumed to be always “on” trying to attach to a gNB, it is necessary to go beyond the device power consumption and look at energy consumption during the time needed to establish a directional link between the user equipment and the gNB. In Fig. 10, we show the energy consumption of the three front-ends, where the user equipment array size is \(N_{\rm{RX}} = 16\). The difference between analog and the two versions of digital is astonishing. The energy consumption of analog exceeds even that of high-resolution digital. Remember, when we looked at the front-ends’ power consumption (Joules/s) in Sect. 3.2.2, the consumption of digital beamforming was six times that of analog beamforming. But, as Fig. 10 establishes, the analog beamforming front-end burns almost four times more energy than even the high-resolution digital front to compensate for the “spatial blindness” during the beam sweep procedure.

However, since after establishing the directional link, in the data communication phase, analog is expected to burn less power than fully digital, we believe that mmWave systems will need to use a low-resolution fully digital front-end; despite the fact that the difference in energy consumption between high- and low-resolution digital front-ends is quite small during beam discovery.

Figure 11:
figure 11

Energy consumption as a function of the user equipment array size. a Energy consumption at the cell edge. b Energy consumption of the median user equipment.

Next, we turn our attention to the number of antenna elements at the receiver. We would like to know how the energy consumption scales with the size of the antenna array in the initial beamsweeping phase. In other words, what is the impact of various beamspace sizes on energy consumption? Fig. 11 answers this question. We try three different configurations at the user equipment: a 2D array of size 4, a 2D array of size 8, and a 2D array of size 16. There are a few interesting points these two figures raise:

  • First, while fewer antennas do bring the energy consumption of the user equipment down, analog beamforming still burns considerably more energy than both the digital front-ends.

  • Second, the gap between analog and digital beamforming widens as more antennas are employed at the user equipment. The reason is the following: when the number of antennas is small, and the post-BF SNR is low, analog beamforming looks into fewer directions at each period. This means that the most crucial factor in the front-ends’ energy consumption is the size of the beamspace. This leads to higher beam discovery delays resulting in longer “on” time for the user equipment. This is more evident in Fig. 11b where all front-ends need no more than \(K= 1\) beamsweep to detect the true path to the gNB.

  • Related to the previous point, the value of digital beamforming is evident in every antenna configuration. The reduction in the size of the effective beamspace is the fundamental and necessary aspect of expediting beam discovery and, subsequently, lowering energy consumption. Beyond initial link establishment, digital beamforming will reduce the complexities involved in handovers and avoiding blockage in connected mode.

5 Conclusion and Future Directions

The realization of mmWave communication systems requires addressing two crucial challenges, namely beam discovery delay and energy consumption. In this work, we have revealed how closely these two are tied together. Specifically, we showed that employing analog beamforming can significantly increase both the discovery delay and, subsequently, energy consumption. Hence, analog beamforming does not buy us lower energy consumption.

Digital beamforming, on the other hand, achieves lower delays and lower energy consumption during beam discovery. In addition, depending on operating SNR, using digital beamforming can allow the use of multi-stream communication and allow joint, flexible scheduling of frequency resources and directional beams. This may be extremely useful in the case of small data packets. By expediting beam discovery, fully digital may enable more aggressive use of sleep mode since beam tracking and paging can be less frequent. In fact, studying the effect of beamforming architectures on higher layers, i.e., MAC and TCP/IP, in a multi-cell setup is an interesting direction for future research.

It is true that when the directional data link is established, analog burns less energy compared to fully digital. Reducing the quantization resolution of the ADC pairs, however, can bring the energy requirements in data communication at the level of analog beamforming at a penalty negligible in most SNR regimes. Thus, low-resolution digital beamforming enables the realization of mmWave cellular systems with low power, low control latency, and more flexible data scheduling.