Cell Selection Technique for Millimeter-Wave Cellular Systems with Hybrid Beamforming

In this paper, a cell selection technique for millimeter-wave (mm-wave) cellular systems with hybrid beamforming is proposed. To select a serving cell, taking into account hybrid beamforming structures in a mm-wave cellular system, the angles of arrival and departure for all candidate cells need to be estimated in the initialization stage, requiring a long processing time. To enable simultaneous multi-beam transmissions in a multi-cell environment, a cell and beam synchronization signal (CBSS) is proposed to carry beam IDs in conjunction with cell IDs. A serving cell maximizing the channel capacity of the hybrid beamformer is selected with the estimated channel information and the optimum precoder. The performance of the proposed technique is evaluated by a computer simulation with a spatial channel model in a simple model of a mm-wave cellular system. It is shown by simulation that the proposed technique with the CBSS can significantly reduce the processing time for channel estimation and cell selection, and can achieve additional gains in channel capacity, or in bit error rate, compared to that obtained by conventional techniques.


Introduction
Owing to the wide use of smartphones and sensor nodes for various applications and services, mobile data traffic is growing rapidly. In particular, sensor nodes such as CMOS sensors and charge-coupled device (CCD) cameras require a high data rate transmission. Very large arrays and/or higher frequency spectrum are recognized to accommodate the increase in mobile traffic. Massive multiple-input multiple-output (MIMO) is expected to bring large gains in terms of spectral and energy efficiency and robustness to hardware failures and impairments [1][2][3]. Further, millimeter wave (mm-wave) communication is a promising technology that allows the use of a wide bandwidth for the high data rates required for next-generation outdoor cellular systems [4,5]. However, mm-wave signals experience severe large-scale path losses compared to lower-frequency band signals [5,6]. To compensate for the large path losses, and to increase the transmission range in the mm-wave frequency band, highly directional beamforming antennas are necessary at both the base station (BS) and the mobile station (MS). Fortunately, a decrease in wavelength enables a large number of antenna arrays to be packed into small form factors, making it feasible to realize the large arrays needed for high beamforming gains [7,8].
Traditionally, MIMO processing is often performed digitally in baseband to support multi-stream service in multipath environments. To perform digital beamforming at mm-wave frequencies, each antenna element requires dedicated baseband and radio frequency (RF) hardware, resulting in high costs and a substantial amount of power consumption. This has constrained mm-wave systems to rely heavily on analog or RF processing [8,9]. Analog beamforming has hitherto been used in the mm-wave systems; however, it is suboptimal compared to digital beamforming. As analog beamforming relies entirely on RF domain processing, the performance is limited by quantization of phase angles and the support of only single stream transmissions [10,11]. In order to overcome the issues faced by digital and analog beamforming, an alternative beamforming architecture, known as hybrid beamforming, has been introduced [12,13]. In hybrid beamforming, beamforming is partitioned into analog and digital domains. Because of the additional digital processing, more degrees of freedom are available in designing the hybrid precoder compared to analog-only beamformers. In addition, hybrid beamforming allows for the supporting of multi-stream and multi-user transmissions.
However, in mm-wave systems, analog beamforming has been used to date for beam training and cell selection in the initialization stage. In the standards of wireless LAN and wireless PAN (IEEE 802.11ad and IEEE 802.15c) widely used for wireless sensor networks (WSNs), analog beamforming is used for both beam training and data transmission using the 60 GHz unlicensed spectrum [14][15][16][17]. In the pre-5G specification for 5G mm-wave cellular systems, called KT Pyeongchang 5G Technical Specification (5G-SIG), analog beamforming is used for beam training and cell selection in the initialization stage. However, both analog beamforming and hybrid beamforming are allowed in the data transmission stage [18,19]. The 5G-SIG, approved by leading world companies in wireless communication, was released in 2016 to demonstrate 5G pilot services for the Pyeongchang Winter Olympic Games in February 2018. In 5G-SIG directional beamforming, antennas are used to maintain the link budget during the initialization stage. The serving cell is found by time-division beam switching (TDBS) in which transmit (Tx) beams of candidate BSs, and receive (Rx) beams of MS, are swept to measure signal-to-noise ratios (SNRs) for all possible links. The serving cell and Tx-Rx beam pair with the highest SNR in the initialization stage are selected. However, the conventional cell selection technique, in which the serving cell is selected based on the measurements of Tx-Rx beam pairs for all possible links (candidate cells), may not be optimal in multipath environments. For example, in a multipath channel where one line-of-sight (LoS) and multiple non-line-of-sight (NLoS) paths are available in candidate cells, only a single path providing the highest SNR (possibly LoS) will be used to select the serving cell in the conventional technique. However, the selected cell in the initialization stage may not be optimal in the data transmission stage where hybrid beamforming is employed. So far, to the best of our knowledge, there has been no attempt to select a serving cell in the initialization stage by considering the hybrid beamforming structure to be used in the data transmission stage.
In this paper, a cell selection technique for mm-wave cellular systems in the initialization stage is proposed for hybrid beamforming structures in data transmission stage. To select a serving cell, taking into consideration the hybrid beamforming structure in the initialization stage, we need to estimate angles of arrival (AoA) and angles of departure (AoD) for all candidate cells. The processing time required for channel estimation increases proportionally to the product of the number of Tx beams, the number of Rx beams, and the number of candidate cells [20]. This long processing time will create significant overheads for a moving MS because the channel estimation should be performed periodically for a possible handover. A new channel estimation technique is proposed by transmitting multiple beams simultaneously, where each beam is mapped with a unique proposed cell and beam synchronization signal (CBSS), all in the physical layer. In the proposed CBSS, the beam ID (BID) is designed in conjunction with cell ID (CID), because the BID must be detected in a multicell environment. It is shown that multiple beams with different BIDs can be transmitted simultaneously with minimal inter-beam interference in a multicell environment, resulting in a significant reduction in processing time for channel estimation and cell selection. The hierarchical structure of CBSS gives the advantage of reusing the same BIDs in adjacent cells, allowing us to generate a large set of sequences. In the proposed cell selection technique, the serving cell is selected based on the maximum channel capacity of the hybrid beamformer with the estimated channel information and optimum precoder. Moreover, the serving cell and the corresponding AoDs and AoAs are determined by taking into consideration the number of RF chains that will be used in the data transmission stage.
The remainder of the paper is organized as follows. Section 2 summarizes the conventional cell selection technique for mm-wave cellular systems and the system model that will be used in this paper. In Section 3, the cell selection technique for mm-wave cellular systems with hybrid beamforming is proposed, after describing the properties of CBSS and channel estimation technique using the CBSS. In Section 4, the performance of the proposed technique is verified by computer simulation using a simple model of a mm-wave cellular system. Finally, conclusions are drawn in Section 5. The following notations are used throughout the paper: A is a matrix, a is a vector, and A is set.
[·] * , (·) , · , · , and < > are conjugate, conjugate transpose, floor operation, ceiling operation, and Jacobi symbol, respectively. N TX and N B denote the number of antennas and beams at the BS while N RX and N I denote the number of antennas and beams at the MS, respectively. N S is the number of data streams and N RF is the number of RF transmit chains. c is the cell ID, c is reference cell ID, b is the beam ID, and b is reference beam ID. I is the identity matrix and 1 N is the N-length all-ones vector. A ⊗ B is the Kronecker product of A and B.

Conventional Cell Selection Technique
A cell search procedure that will be used for an MS in conventional LTE systems can be summarized as follows. The MS first needs to scan by measuring the received signal strength indicators (RSSIs) of all supported carriers for entire bandwidth, each with 20 MHz. The measured RSSIs include the signal strength of the desired cell, interference cells, and noise. Then, RSSIs lower than the predefined threshold (or signal strength sensitivity) are filtered, resulting in a list of surviving candidates. Each surviving candidate is processed by cell identification block until the MS finds a cell on which to camp. Cell identification involves primary synchronization signal (PSS) and secondary synchronization signal (SSS) detection. The SSS carries the information of cell ID group, while the PSS carries the information of cell ID within the cell ID group. A physical cell ID is constructed by combining the cell ID group and cell ID within the cell ID group after PSS and SSS detection.
However, in mm-wave cellular systems, the selection of the best beam pair, as well as the cell selection, is important because a maximum array gain is obtained in these systems when Tx and Rx beams are perfectly aligned [14,15]. Especially for systems with narrow beams, a small mismatch or misalignment between Tx and Rx beams may result in a substantial loss in the received power. In mm-wave cellular systems, the TDBS scheme has been used to find the serving cell with the best beam pair in the initialization and handover stages [18,19]. In the TDBS scheme, individual Tx beams are transmitted from the BS until all the Tx beams are transmitted. The Rx beam sweep is performed at the MS for each Tx beam to measure the SNR for each Tx-Rx beam pair. The measurement of the SNRs for all possible Tx-Rx beam pairs must be performed for all candidate cells to select a serving BS with the best beam pair. Therefore, the serving cell and best beam pair are selected based on the measurement results obtained by all possible Tx-Rx beam pairs of candidate cells. A physical cell ID of the serving cell is obtained by detecting the PSS and SSS transmitted from the serving cell. A beam pair with maximum SNR is selected as the best beam pair and the corresponding cell with the physical cell ID as the serving cell. However, the selected cell and beam pair may not be optimal for hybrid beamforming system because the effect of the multipath channel is not considered in the initialization stage [10,11,14]. Figure 1 shows an example of a mm-wave cellular system with two BSs with different CIDs. Here, it is assumed that a multipath (NLoS) component exists in cell 0 and an MS is located at the cell boundary. It is also assumed that the LoS path is stronger than the NLoS path in cell 0. Here, switched beamforming is used at both the BSs and the MS. In this example, a serving cell will be selected by measuring the signal strengths of LoS paths from two cells. The cell with the higher signal strength of the LoS path will be selected [18,19]. Here, the multipath (NLoS) component is not used for cell selection. However, in hybrid beamforming systems, the multipath channel is used to increase the SNR in power-limited conditions (e.g., cell edge), or to increase user data rates, as long as the multipath channel has sufficiently different propagation paths. Therefore, the serving cell selected in the initialization stage may not be optimal for hybrid beamforming system in the data transmission stage.

System Model
Consider a mm-wave cellular system in which a transmitter with N TX antennas transmitting N S data streams to a receiver with N RX antennas. It is assumed that, in the data transmission stage, hybrid beamforming is used at the BS, and analog beamforming with a set of predefined angles is used at the MS. The BS is equipped with N RF transmit chains to enable N S data streams transmission, such that N S ≤ N RF ≤ N TX [12,21]. This architecture enables the transmitter to apply an N RF × N S baseband precoder by using N RF transmit chains, followed by an N TX × N RF RF precoder using analog circuitry. In the proposed technique, it is assumed that all Tx beams are simultaneously transmitted from all candidate BSs, and Rx beams are swept over in this period, to reduce the processing time required for channel estimation. All Tx beams are transmitted repeatedly until one round of Rx beam sweep is completed. In the proposed technique, CBSSs are mapped to the corresponding Tx beams, simultaneously transmitted from candidate BSs. Because the BID needs to be detected in a multicell environment with multiple beams, the CBSS is designed such that the BID is mapped on a synchronization sequence in association with its CID. Both the BID and CID are transmitted on a single CBSS to facilitate joint detection of the BID and CID. This hierarchical structure enables the MS to identify the CID of the beam from the received CBSS. In a mm-wave cellular system based on orthogonal frequency division multiplexing (OFDM), the signals received by the i-th Rx beam at the MS in the initialization stage can be expressed in the frequency domain as follows: where q i ∈ C N RX ×1 and p c,b ∈ C N TX ×1 denote an analog combiner for the i-th Rx beam at the MS, and an analog precoder for the b-th beam in the c-th cell, respectively. Here, N I , N B , and N C denote the number of Rx beams at the MS, the number of Tx beams at the BS, and the number of candidate cells, respectively. H c (k) is an N RX × N TX matrix at subcarrier k that represents the mm-wave channel between the N B beams of the c-th cell and the N I Rx beams of the MS. X c,b (k) and W(k) denote the CBSS transmitted from the b-th beam of the c-th cell and additive white Gaussian noise (AWGN), respectively. If a total number of Tx beams at the BS is N B and the number of available RF transmit chains is N RF (≤ N B ), then only N RF beams can be simultaneously transmitted, and N I Rx beams are swept over in this period. The same procedure needs to be repeated N B /N RF times by switching the set of N RF beams until all the Tx beams are transmitted. For example, if the values of N B , N RF , and N I are 16, 1, and 8, respectively, then only one beam will be transmitted until one round of Rx beam sweep is completed. This procedure will be repeated 16 times. If the value of N RF is 8, then 8 beams will be simultaneously transmitted until one round of Rx beam sweep is completed. This procedure will be repeated only 2 times. Thus, the processing time required for channel estimation increases as the number of RF transmit chains decreases. In this paper, it is assumed that the number of RF transmit chains is the same as the number of Tx beams for notational convenience.
In Equation (1), the matrix H c (k) represents the channel at subcarrier k, which is derived from the delay-τ channel matrix H c (τ). Here, the delay-τ channel matrix can be considered as a geometric wideband channel model, which can be viewed as the sum of N M scattering clusters each with time delay τ m . Each cluster consists of N L propagation rays with relative time delay τ l . The delay-τ channel matrix H c (τ) can be written as [22]: and the channel at subcarrier k, H c (k), can be expressed as [23]: where N L and N M denote the number of propagation rays and the number of clusters, respectively. p rc (τ) is the pulse shaping function for T S -spaced signaling evaluated at τ second, β c m,l is the complex gain of the l-th ray in the m-th cluster, whereas a RX α c ml , ϕ c ml and a TX φ c ml , θ c ml represent the Rx and Tx array response vectors at an azimuth (elevation) angle of α c ml ϕ c ml and φ c ml θ c ml , respectively. The array response vector for the N-element uniform linear array (ULA) on the y-axis and uniform planar array (UPA) in the xy-plane with N e and N f elements on the x and y axes, respectively, is given by [24]: where d is the element spacing and λ is the wavelength. Also, 0 ≤ e < N e and 0 ≤ f < N f are antenna element indices for x and y, respectively.

Cell and Beam Synchronization Signal (CBSS)
In order to select a serving cell, we need to estimate channels (AoAs and AoDs) for all candidate cells in the initialization stage. In the proposed technique, CBSSs are simultaneously transmitted from all Tx beams of candidate BSs to reduce the processing time required for channel estimation. Here, it is assumed that the synchronization and CID detection are performed by the conventional primary synchronization signal (PSS) and secondary synchronization signal (SSS), used in LTE systems [25]. In addition, it is assumed that one additional OFDM symbol carries the proposed CBSS, next to the SSS. The CBSS needs to carry the CID, as well as the BID, in order to be detected in the multi-cell environment. If the CBSS carries only BID information, it will be impossible to determine which cell had transmitted the BID. Therefore, a large number of different sequences (CID × BID) need to be generated. In order to transmit the information in one OFDM symbol, the BID is mapped in association with the CID in the CBSS. The proposed CBSS can be viewed as a combination of a Chu sequence and a polyphase sequence [26]. In the proposed CBSS, a new sequence is generated by mapping the CID to the root index of the Chu sequence and the BID to the index of the polyphase sequence. Figure 2 shows the generation process of the CBSS. The CBSS is defined by the product of a prime-length Chu sequence X c (k) and a polyphase sequence P b (k) in the frequency domain as follows: where, c ∈ N C , b ∈ N B , r c ∈ N R , N C = 0 1 . . . where Nz r c χγ 2 , χ = N z +1 2 , γ = N z −1 2 . The autocorrelation function of the CBSS in the frequency domain is given by: Nz (r c kl) = N z e jπ Nz (r c l(l+1)−2Lbl) δ(mod(l, N z )) (8) where δ(·), (r c ) −1 , and χr c /N z denote the Kronecker delta function, the multiplicative inverse of r c , and the Legendre symbol with a value of −1, 0, or 1, respectively. As seen in Equation (7), the time domain version of the CBSS can be represented as a Chu sequence with a cyclic shift and its amplitude is given by x c,b (n) = 1/ √ N z for all n. From Equations (7) and (8) The intra-cell correlation (correlation between two beams b and b in the same cell) and inter-cell correlation (correlation between a beam in one cell and a beam in the adjacent cell) are analyzed in an ideal condition, ignoring the channel and noise effects, as follows: and where c,c Equation (10) can be derived using the Gauss sum property [27]. From Equations (9) and (10), it can be observed that the intra-cell correlation has a peak amplitude N z only when b = b for b, b ∈ N B , and the amplitude of the inter-cell correlation is given by √ N z . As can be seen in Equations (9) and (10), CBSSs with the same Chu's root indices are perfectly orthogonal, and CBSSs with different Chu's root indices are not orthogonal. However, the interference term caused by the beam transmitted from neighboring BSs is relatively small compared to the signal from the reference beam. For example, the ratio between √ N z and N z is 3% when N z = 1021. When an STO (δ t ) is present, the intra-cell correlation is given by: The correlation value in Equation (12) becomes high when the following condition is satisfied: where ∆b is defined as b − b. To avoid the condition in Equation (13), the parameter L needs to be adjusted and can be set to the following value: where δ t max and L max denote the maximum STO and maximum value of parameter L, respectively. The value of L must increase proportionally with respect to the STO, because of the time delay between the BS and MS. This can result in the reduction of the number of BIDs to N z /L .

Channel Estimation and Cell Selection
Using the CBSSs, channels of candidate cells are estimated in the initialization stage. The channel estimation begins with the correlation between the received signal in Equation (1) and the known CBSSs. The correlation value between the signal received from the i-th Rx beam at the MS and reference CBSS is given by: where X c,b (k) denote the reference CBSS with the reference CID c and BID b. The correlation in Equation (15) should be performed for N C cells and N B Tx beams. Here, the vector q i is constructed to form an i-th Rx beam in the direction of m-th cluster, thereby receiving the signals only from the m-th cluster. If the N L propagation rays contributed by the mth cluster are received within the sampling interval (T S ), the channel H c (k) with Tx and Rx beamforming can be written as H c . For example, when only two cells (cell 0 and cell 1) are present, the correlation in Equation (15) is given by: where the reference CID, c, is set to 0 and the noise effect is ignored. From the first two terms in Equations (15) and (16), it can be observed that there is no interference from other beams in the same cell. That is, the correlation value is high (N z ) when BIDs are the same, and the correlation value is zero when other BIDs are used. The third term in Equations (15) and (16) represents the interference term ( √ N z ) from beams in the other cells. In a vector form, q i H c p c,b can also be written as follows: where β c is unknown path gains of N M N L propagation paths. The array response matrices (unitary discrete Fourier transform [DFT] matrices) at the BS and MS are defined for ULAs and UPAs, as follows: where A TXD ∈ C N TX ×G TX and A RXD ∈ C N RX ×G RX are defined for ULA. Here, G TX and G RX denote uniform fine grid points in [−π/2, π/2) at the BS and MS, respectively. Also, and A RXD ∈ C N RX ×(G φ RX ×G θ RX ) are defined for the UPA. Here, G φ TX G θ TX and G α RX G ϕ RX denote uniform fine grid points in the azimuth (elevation) domain with [−π/2, π/2) ([0, π)) at the BS and MS, respectively. In Equation (17), A c is an N B N I × G TX G RX and matrix for ULA and UPA, respectively. In Equation (2), it is assumed that azimuth (elevation) AoA and AoD, α c ml ϕ c ml and φ c ml θ c ml , are taken from the above-mentioned uniform fine grids. Next, all the N B N I correlation values are stacked in a vector form as follows: Equations (17) and (19) are used to obtain AoAs and AoDs of channel paths. The orthogonal matching pursuit (OMP) algorithm is used to obtain these values because of its simplicity and fast implementation [23,28]. Stagewise OMP (StOMP) or generalized OMP (gOMP) algorithms can also be considered for estimation of AoAs and AoDs of channel paths [29,30]. Using the vector a c obtained by OMP, comprised of indices of AoAs and AoDs for all N M N L propagation paths, multipath gains β c can be estimated as follows: Then, the estimated channel matrix can be constructed for each candidate cell as follows: Finally, a serving cell should be selected in the initialization stage, taking into consideration the hybrid beamforming structure, because a hybrid beamforming system will be used in the data transmission stage. Here, we select the serving cell that maximizes channel capacity in the hybrid beamforming system. The unconstrained optimal precoder F c opt for hybrid beamformer is obtained by performing a singular value decomposition (SVD) of the estimated channel in Equation (22) as H c = U c Σ c V c * [13]. Here, U c is an N RX × N RX unitary matrix, Σ c is an N RX × N TX diagonal matrix, and V c is an N TX × N TX unitary matrix. The matrix V c can be expressed as: The optimal unconstrained precoder for H c is given by F c opt = V c 1 . In general, the channel capacity of a MIMO system is expressed in terms of the channel matrix H c and precoder F c opt , as follows [31]: When N RF Tx chains are available in the hybrid beamformer at the BS, Equation (24) can be rewritten as: By repeating the same procedure for all candidate cells, the cell providing the maximum channel capacity is selected as the optimal serving cell as follows: After the serving cell is selected, the estimated vectors of AoDs (φĉ) and AoAs (αĉ), corresponding to the indices of the array response vector, in the azimuth plane are given as follows: where 1 1×N RF is a 1 × N RF vector with all values of 1.

Simulation
In this section, the performance of the proposed cell selection technique is evaluated by computer simulation in a simple scenario shown in Figure 1. Here, an MS is located at the cell boundary between two cells (CID = 0 with root index = 54, and CID = 1 with root index = 59). The system bandwidth, carrier frequency, and fast Fourier transform (FFT) size are set to 500 MHz, 28 GHz, and 2048, respectively [9,18,32]. For ULA, the number of antenna elements is set to 16 and 8 in the BS and MS, respectively. For UPA, the number of antenna elements is set to 64 N e = N f = 8 and 16 N e = N f = 4 in the BS and MS, respectively. In addition, N B = 16 and N I = 8 for both ULA and UPA. The length of CBSS sequence, N z , is set to 1021. It is assumed that an LoS path exists in both cells. In addition to the LoS path, there exist NLoS paths in cell 0 and the number of NLoS paths (clusters) is set to 4 [33]. For simulation, the spatial channel model (SCM) is used with a k-factor of 6 dB [34]. The number of rays in the NLoS path in a Rician fading channel is set to 20 (N L = 20). In Figure 3, beampatterns of unitary DFT matrices with N B = 16, used in the simulation, are shown for ULA and UPA. Here, G TX = G φ TX = 180 and G θ TX = 1. The fine grid point is increased with 1 degree of resolution in azimuth angle. For UPA, the elevation angle is fixed to π/4 for convenience.
Before applying the proposed technique to the mm-wave cellular system, the properties of the proposed CBSS, discussed in Section 3, are verified. In Figure 4, the correlation property of the CBSS is shown when STO is not present. Here, the reference CID and BID are set to 0 r c = 54 and 7, respectively. The parameter L is set to 1. From Figure 4, it can be seen that the intra-beam correlation in the same cell (CID = 0) is zero except in the case of the same BID (BID = 7), as given by Equation (9). The normalized inter-cell correlation with the CID equal to 1 r c = 59 is equal to 0.03 1/ √ N z for all values of BID, as given by Equation (10). Figure 5 shows the correlation property of the CBSS in the time domain when STO exists. The reference CID and BID are set to 0 r c = 54 and 0, respectively. It can be seen that when the STO is zero, correlation values are zero except for the lag l = 0. When the values of STO are 1 and 2, the peak shifts by one and two lags, respectively. So, the shift increases proportionally to the STO value. If the parameter L is set to 1, there may be an error in BID detection. Therefore, the parameter L satisfying Equation (13) should not be selected when the CBSS is designed. To avoid misdetection in the channel environment (δ t max = 6), the parameter L is set to 7 in the simulation. Figure 5 shows the correlation results when the parameter is used. Peaks can be found at 7 and 14, when CBSSs with BIDs equal to 1 and 2 are used, respectively.    Figure 6 shows the achievable channel capacity when the conventional and proposed techniques are used. In the conventional technique, the serving cell is selected by measuring SNRs for all candidate cells with Tx/Rx beam sweeping, and choosing the cell with the highest SNR. Cell 1 is selected as the serving cell in the conventional technique because cell 1 has a stronger LoS path than cell 0. However, in the proposed technique (LoS + NLoS(1), LoS + NLoS(2), LoS + NLoS(4)), cell 0 will be selected as the serving cell. That is because more multipath channels are considered for cell selection in the proposed technique as the number of RF chains (N RF ) increases. When N RF is set to {2, 3, 5}, an additional gain of {0.64, 1.8, 3.23} and {1.24, 3.64, 6.6} bps/Hz at an SNR of 9 dB can be achieved for ULA and UPA, respectively. In the proposed technique, the serving cell is selected using Equations (25) and (26).  Figure 7 shows the achievable channel capacity when the proposed technique is used with varying k-factor and different antenna configurations. Figure 7a shows the results for ULA with 16 × 8 and 32 × 16 antenna configurations and Figure 7b shows the results for UPA with 64 × 16 and 16 × 8 antenna configurations. Here, the channel capacity is measured at SNR 9 dB when k-factor varies from 2 to 20 with a step size of 2. As can be seen in Figure 7, spectral efficiencies for all cases (LOS, LOS + NLOS(1), LOS + NLOS(2), LOS + NLOS(4)) converge towards LOS case, irrespective of ULA/UPA. That is because LOS path becomes stronger than other NLOS paths as k-factor increases. The spectral efficiency increases as k-factor decreases, because NLOS components become more significant. When the k-factor is equal to 6 dB, the ULA configuration with 32 × 16 obtains 2.0-2.5 bps/Hz gain in spectral efficiency over the ULA with 16 × 8, while the UPA configuration with 64 × 16 obtains 4.0-5.5 bps/Hz gain over the UPA with 16 × 8.  Figure 8 shows the bit error rate (BER) performance of the mm-wave system in the data transmission stage after a serving cell is selected by two different techniques. Here, QPSK is used for modulation. For reference, the analytic curve for a Rician fading channel with a k-factor equal to 6 dB is inserted. Note that the conventional technique (LoS) is better than the analytic result, because the SCM channel model can be approximated as an AWGN channel in the simulation when the beamforming is performed in the LoS direction. In the analytic result (Rician fading), only the array gain (16 × 8 in ULA, 64 × 16 in UPA) is taken into account. In the conventional technique, BER is measured in cell 1, because cell 1 is selected as the serving cell in the initialization stage. However, the BER is measured in cell 0 in the proposed technique (LoS + NLoS(1), LoS + NLoS(2), LoS + NLoS(4)). As can be seen in Figure 8, an additional maximum-ratio combining (MRC) gain of {1, 2, 3} dB can be achieved in the proposed technique when N RF is set to {2, 3, 5}, respectively.  Table 1 compares the processing time and computational complexity required for channel estimation and cell selection when the conventional and proposed techniques are used. Here, T SS and T CBSS represent the lengths of the synchronization signal (SS) and CBSS symbol, respectively. It is assumed that T SS is equal to the length of two OFDM symbols (PSS and SSS), as in LTE systems, and T CBSS is equal to one OFDM symbol length. In this example, the number of candidate cells, the number of Tx beams at BS, and the number of Rx beams at MS, are assumed to be N C = 3, N B = 16, and N I = 8, respectively. From Table 1, it can be seen that the processing time required for channel estimation and cell selection with the proposed technique is much smaller than with the conventional technique. In the proposed technique, the processing time decreases proportionally to the number of RF transmit chains. In the conventional technique, the processing time increases linearly with the increase in the number of candidate cells, whereas it does not change in the proposed technique. Next, the computational complexities required for channel estimation are compared. In the conventional technique, CIDs are carried by PSS and SSS, which are transmitted from each beam. Beam sweep is performed at both the BS and MS. However, in the proposed technique, CIDs are carried by PSS and SSS, and BIDs are carried by CBSS in a hierarchical manner. In addition to PSS and SSS, CBSS is also transmitted from each beam simultaneously. Here, beam sweep is performed only at the MS. In the proposed technique, the number of RF transmit chains at the BS is considered to compute computational complexity. Here, N represents a prime number less than minimum bandwidth in LTE and is set to 61. Table 1 shows that, as the number of RF transmit chains increases, the complexity required for channel estimation decreases because multiple beams can be transmitted simultaneously from the BS with minimal interbeam interference, and channel estimation is performed jointly at the MS.

Conclusions
In this paper, a cell selection technique for mm-wave cellular systems with hybrid beamforming is proposed using the CBSS in the initialization stage. It is shown by simulation that the CID and BID can be detected correctly in a multicell environment with STOs, by assigning an appropriate value to the parameter (L) in the CBSS. The proposed technique with the CBSS is shown to require only 6% of the processing time for channel estimation and cell selection, compared to the conventional technique, when the number of RF transmit chains is 8. It is also shown that the proposed cell section technique can achieve 3.64 bps/Hz (UPA) gain in channel capacity or 2 dB gain in BER, when N RF = 3. Although the cell selection technique is described in this paper for mm-wave cellular systems with hybrid beamforming in the initialization stage, it can also be applied to cell searching in the handover period.