# Optimal Multiuser Loading in Quantized Massive MIMO under Spatially Correlated Channels

Low-resolution digital-to-analog converter (DAC) has shown great potential in facilitating cost- and power-efficient implementation of massive multiple-input multiple-output (MIMO) systems. We investigate the performance of a massive MIMO downlink network with low-resolution DACs using regularized zero-forcing (RZF) precoding. It serves multiple receivers equipped with finite-resolution analog-to-digital converters (ADCs). By taking the quantization errors at both the transmitter and receivers into account under spatially correlated channels, the regularization parameter for RZF is optimized with a closed-form solution by applying the asymptotic random matrix theory. The optimal regularization parameter increases linearly with respect to the user loading ratio while independent of the ADC quantization resolution and the channel correlation. Furthermore, asymptotic sum rate performance is characterized and a closed-form expression for the optimal user loading ratio is obtained at low signal-to-noise ratio. The optimal ratio increases with the DAC resolution while it decreases with the ADC resolution. Numerical simulations verify our observations.

## Authors

• 9 publications
• 136 publications
• 3 publications
• 12 publications
• 35 publications
• ### Secure Massive MIMO Communication with Low-resolution DACs

In this paper, we investigate secure transmission in a massive multiple-...
01/28/2019 ∙ by Jindan Xu, et al. ∙ 0

• ### Analog Versus Hybrid Precoding for Multiuser Massive MIMO with Quantized CSI Feedback

In this letter, we study the performance of a downlink multiuser massive...
06/01/2020 ∙ by Yaqiong Zhao, et al. ∙ 0

• ### SLNR Based Precoding for One-Bit Quantized Massive MIMO in mmWave Communications

Massive multiple-input multiple-output (MIMO) is a key technology for 5G...
02/12/2019 ∙ by Yavuz Yapici, et al. ∙ 0

• ### Spatially Correlated Massive MIMO Relay Systems with Low-Resolution ADCs

In this paper, we investigate the massive MIMO relay system, where the r...
04/16/2020 ∙ by Peihao Dong, et al. ∙ 0

• ### Massive MIMO with Dense Arrays and 1-bit Data Converters

We consider wireless communication systems with compact planar arrays ha...
12/04/2020 ∙ by Amine Mezghani, et al. ∙ 0

• ### Low Resolution Digital-to-Analog Converter with Digital Dithering for MIMO Transmitter

Based on an equivalent model for quantizers with noisy inputs recently p...
04/19/2019 ∙ by Arkady Molev-Shteiman, et al. ∙ 0

• ### On the Performance of Massive MIMO Systems With Low-Resolution ADCs Over Rician Fading Channels

This paper considers uplink massive multiple-input multiple-output (MIMO...
06/24/2019 ∙ by Tianle Liu, et al. ∙ 0

##### This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

## I Introduction

Massive multiple-input multiple-output (MIMO) has gained significant attention as a candidate technique for the next generation wireless system [2]-[4]. In massive MIMO, a large amount of antennas equipped at base station (BS) can provide high spectral and energy efficiencies [5]. Despite these merits of massive MIMO, it suffers from a challenging issue of high cost and power consumption for applications even at the BS. This is due to the fact that each antenna has to be driven by a separate radio-frequency (RF) chain.

Considering that the power consumption of each RF chain can decrease dramatically by reducing the resolutions of digital-to-analog converters (DACs), one of the potential solutions is to employ low-resolution DACs for downlink transmissions [6]-[11]. More specifically in [6]-[8], nonlinear precoding schemes were proposed for massive MIMO downlink transmission with low-resolution DACs. In [6], a novel precoding technique using 1-bit DACs was presented to mitigate multiuser interference and quantization distortions. In [7], a computationally-efficient 1-bit beamforming algorithm, named POKEMON (short for PrOjected downlinK bEaMfOrmiNg), was proposed. The authors of [8]

studied perturbation methods which minimize the probability of errors at receivers in a massive MIMO downlink with 1-bit DAC. Alternatively, linear precoding techniques were studied for the downlink transmissions using low-precision DACs

[9]-[11]. A linear precoding approach was studied in [9], where the output data of conventional linear precoders were directly quantized by low-resolution DACs. Considering 1-bit DACs and zero-forcing (ZF) precoder, [10] analyzed the system performance for massive MIMO downlink transmission. While in [11], the system performance was investigated for the downlink using multi-bit DACs with common precoding schemes, including regularized zero-forcing (RZF) and maximal-ratio-combining (MRC) precoders.

As for uplink channels, low-resolution analog-to-digital converters (ADCs) were adopted to reduce the hardware and power cost [12]-[15]

. The impact of 1-bit ADC on channel estimation was considered in

[12] and satisfactory system performance was observed in terms of both symbol error rate (SER) and mutual information. The performance analysis was then extended in [13] for a multiuser relay network using mixed ADCs. Spectral efficiency of massive MIMO with low-resolution ADCs was studied in [14] while an extension to mixed-ADC architecture was then conducted in [15] under energy constraint.

Moreover, spatial correlation at the BS side generally exists in massive MIMO since the large number of antennas could not be set far away enough from each other due to the restriction of device size. It was revealed in [20] that the capacity decreased about 20 percent with a separation of four wavelengths between adjacent antennas, compared to an uncorrelated MIMO channel assumed in most of the above literature. Separable correlation model for a single-user communication system was studied in [21]. A more general correlation model, i.e., Unitary-Independent-Unitary (UIU) model, was introduced in [22]

. By studying the distribution of the eigenvalues of a correlated MIMO channel matrix, closed-form approximations of asymptotic channel capacity were derived in

[23],[24]. Considering the transmit-side channel correlation, large system analysis for a downlink multiuser broadcast channel was conducted in [25]. Using low-resolution ADCs at the BS, uplink performance of a massive MIMO network with spatial correlations was studied in [26].

Most of the existing studies focused on low-resolution converters at BS, i.e., few-bit DACs for downlink or ADCs for uplink. In this paper, we investigate the joint effect of low-resolution DACs at BS along with finite-resolution ADCs at users, extending the previous conference paper in [1]. We adopt a separable spatial correlation model, which is well-known as the Kronecker model [21], to characterize the spatial correlation at the BS. At the user sides, the correlation is not considered under the assumption that single-antenna users lie far enough from each other in a rich scattering environment. Perfect channel state information (CSI) is assumed known at the BS for analytical simplicity, considering there are already a number of studies, i.e., [27]-[29], specifically dealing with the problem of channel estimation. For performance enhancement, we first optimize the regularization parameter for RZF precoding design using the asymptotic random matrix theory. Given the optimal RZF precoder, we then analyze the asymptotic sum rate performance of the massive MIMO system. By maximizing the sum rate per antenna, the optimal user loading ratio is further derived. Especially for low SNRs, deemed as an important scenario in massive MIMO, we obtain a closed-form solution of the optimal user loading ratio. Engineering insights for the optimal ratio are accordingly concluded.

The rest of this paper is structured as follows. Quantization and channel models are introduced in Section II. In Section III, we derive the optimal regularization parameter, asymptotic achievable rate, and the optimal user loading ratio under spatially correlated channels. Section IV analyzes a special case of uncorrelated channels and gives a closed-form solution of the optimal user loading for low SNR. Section V presents simulation results. Conclusions are drawn in Section VI.

Notations: , , and represent the transpose, conjugate, and conjugate transpose of , respectively. denotes the trace of and keeps only the diagonal entries of . is the expectation operator. denotes the almost sure convergence.

## Ii System Model

We consider a multiuser massive MIMO downlink network as illustrated in Fig. 1. The BS equips antennas and it simultaneously serves users. In order to reduce power consumption for driving the massive antenna array, each antenna at the BS connects to a pair of low-resolution DACs to separately process real and imaginary parts of complex signals. Each user equips a single antenna connected to a pair of finite-resolution ADCs. Although it is in general difficult to exactly characterize the nonlinear operation of both DACs and ADCs, the well-known Bussgang theorem [30],[31] is widely used to approximately decompose the quantized data into two uncorrelated parts. One part represents the linearly distorted signal and the other is quantization noise. Let be the Gaussian source data with and denote the precoding matrix with where amounts to the transmit power budget. By applying the Bussgang theorem, the transmit signal after DAC conversion can be modeled as [32]

 x=√1−ρDAPs+nDA, (1)

where denotes the distortion factor which reflects the number of DAC quantization bits , and represents the quantization noise which satisfies

As in most existing works, e.g., [13], [31]-[33], we consider the optimal non-uniform quantization model since it provides a tractable and effective way of well characterizing the quantization performance. For uniform quantizers, the quantization model is more involved [11]. For performance analysis, it has been shown that the resulting performance difference is marginal, by applying the two models especially for popular quantization levels.

At user sides, the stacked received vector,

, of all users equals

 y=Hx+n, (3)

where represents the additive white Gaussian noise (AWGN), and denotes the downlink channel with spatial correlation among transmit antennas, which can be further expressed as

 H=D~HR12, (4)

where is a diagonal matrix and the diagonal entries denote large-scale fadings, the entries of

are assumed independent and identically distributed complex Gaussian random variables with zero mean and unit variance, and

denotes a transmit correlation matrix [34]. Assume that power control based on statistical channel-inverse is conducted to compensate large-scale fadings of different users, as done in [35]. We consider a normalized channel model with . Typically, is Hermitian and positive definite with . We assume that each user experiences the same transmit correlation since it heavily relies on the antenna array and scattering distribution at BS [25],[26]. Similar to DACs, the output data vector after the ADC quantization can be decomposed according to the Bussang Theorem, as follows

where is the distortion factor of ADC quantization, and is the quantization error, which can be expressed as

 (6)

Compared to DACs, an additional scalar factor, , is multiplied for ADCs and we have . Power normalization is not a necessity here at the receiver because it does not affect the equivalent signal-to-interference-quantization-and-noise ratio (SIQNR) [32]. The minimum distortion parameters for the optimal quantization has been studied in [36] and the typical values of and corresponding to various quantization resolutions are exemplified in Table I. In particular for the special resolution of 1-bit quantization, the author in [30] has derived that , which is exactly the same as in Table I for this special case. In the condition of moderate to high-resolution quantizations, we have [31]. The calculated values for are verified, approximately the same as these in Table I.

Now substituting (1) and (3) into (5), the received data vector after ADC quantization equals

Let and denote the th row of and the th column of , respectively. Then, the received signal of the th user can be extracted from (7) as

where , , and are, respectively, the th elements of , , and . In (II), the first summation term is the desired signal of user while the second term represents multiuser interference. The third and forth terms come from the DAC and ADC quantization distortions, respectively. Then, the equivalent SIQNR of the th user, , can be readily characterized as

 (9)

## Iii Asymptotic Performance Analysis

In this section, we analyze the performance of the massive MIMO systems using the asymptotic random matrix theory [37], assuming both and simultaneously grow large while the user loading ratio, , remains invariant. Transmit-side channel correlation is considered and the impacts of both DAC and ADC quantizations are analyzed.

### Iii-a Optimal Regularization Parameter

Without loss of generality, we first characterize the asymptotic behavior of SIQNR at user . Consider a RZF precoder as

 P=c(HHH+αIN)−1HH, (10)

where is the regularization parameter and is a constant guaranteeing the power constraint as

 c= ⎷PTr{H(HHH+αIN)−2HH}. (11)

Note that a typical value of is in conventional MIMO systems [17],[38]. For the massive MIMO with low-resolution DACs at BS and finite-resolution ADCs at users, however, can be further optimized as will be discussed later. Before that, we present the asymptotic behavior of SIQNR with a Toeplitz correlation matrix as given in the following theorem.

###### Theorem 1.

The asymptotic SIQNR, , can be expressed as

where denotes the average system SNR, is defined as the unique solution of the following function

 ξ =Eλ{λ(1+ξ)ρ(1+ξ)+βλ}, (13)

where is defined as the normalized regularization parameter, denotes the eigenvalue of , and is defined as

 Eij≜Eλ{λi[ρ(1+ξ)+βλ]j}. (14)
###### Proof.

See Appendix C. ∎

From (12), is affected by distortion parameters and due to quantizations. In particular, decreases with both and because and . Since and are distortion factors decreasing with quantization resolution, we say that higher-resolution DAC or ADC (smaller or ) achieves larger , which corresponds to better performance. The impact of spatial correlation lies in the values of , , and , which depends on the eigenvalues of .

Note that the normalized regularization parameter, , which is usually set as for MIMO systems in existing literature, plays an important role in performance improvement with RZF. Now, we are ready to optimize for rate maximization considering the impacts of channel correlation, low-resolution DACs, and finite-resolution ADCs. The optimization problem follows

 maxρ   γ. (15)

This optimization problem is equivalent to maximize the sum rate of all users, since the instantaneous SIQNR of each user converges to the same in (12). The closed-form solution to (15) is given in Lemma 1.

###### Lemma 1.

The optimal regularization parameter, , is obtained as

 ρ∗=(ρDAγ0+1)β(1−ρDA)γ0. (16)
###### Proof.

In order to solve the problem in (15), we use (12) and set , which yields

 ~E[(ρDAγ0+1)∂∂ρ[(1+ξ)2E12]+(1−ρDA)γ0∂E22∂ρ]−∂~E∂ρ[(ρDAγ0+1)(1+ξ)2E12+(1−ρDA)E22γ0]=0, (17)

where we define for notational brevity and it follows

 ∂~E∂ρ= ∂ξ∂ρ[E22+ρβ(1+ξ)2E12]+ξ[∂E22∂ρ+(1+ξ)2E12β+ρβ∂∂ρ[(1+ξ)2E12]] = ξ∂E22∂ρ+ξρβ∂∂ρ[(1+ξ)2E12], (18)

where (47) and (58) in Appendix C are utilized. Then, by substituting (III-A) into (17) and after some basic manipulations, Lemma 1 is obtained. ∎

Comparing (16) with the conventional typical value of , we can rewrite , where can be regarded as an equivalent system SNR affected by low-resolution DACs. Analogous to the typical , the derived replaces the original SNR with the equivalent one, i.e., .

Remark 1. It is interesting that in (16) is independent of correlation matrix , which implies that the optimal RZF precoding is not affected by the spatial correlation at the BS. Besides, is independent of which implies that the resolution of receiving ADCs does not affect the optimal RZF precoding. This looks promising since it allows the BS to determine the optimal RZF parameter without knowing the ADC setups at user sides. From (16), increases with while it decreases with .

In general, serving more users directly increases instantaneous sum rate while leading to heavier pilot overhead for channel estimation. It is interesting to optimize the user loading ratio, , by characterizing the rate performance within a period no longer than coherence interval.

Substituting (16) into (12), the maximized SIQNR of each user can be expressed as

where is computed from (13) by substituting in (16). On one hand, the effect of ADC quantization is evident. The more quantization bits ADCs provide, the smaller is and the larger will be. One the other hand, the impact of low-resolution DAC lies in as indicated in (16). Besides, the spatial correlation affects the asymptotic performance through .

By applying the assumption of the worst-case Gaussian interference, the asymptotic achievable rate of each user equals

in (20) represents the achievable rate of each user per channel use without considering the pilot overhead. Assume that each user needs an overhead of to transmit precoded pilot signal for user detection within a coherence interval . The achievable sum rate per antenna including the pilot overhead can then be evaluated from (20) as

where . Then, we can acquire the optimal user loading ratio, , by maximizing , which yields

 β∗=arg maxβ ¯¯¯¯R. (22)

Since in (21) is a continuous function w.r.t , is obtained when the derivative of w.r.t. equals zero, or tends to the critical values that or . Consider the scenario of low SNR with . Substituting to (21), it is obvious that . For and , we have by using (16). Substituting this into (13) yields , which implies . As a result, we can conclude that the maximum of can only be obtained by solving the equation

 ∂¯¯¯¯R∂β=0. (23)

Although it is hard to find a simple closed-form solution of (23) in general cases, we can resort to offline numerical methods, e.g., the bisection method, as in Section V. Since in (22) depends only on system parameters, lookup tables of can be constructed.

Note that in maximizing , we obtain for any value of . Then, we optimize by substituting the obtained in (16), which can be considered as a function w.r.t . Therefore, it is safe to conclude that both the achieved and are globally optimal.

## Iv Performance under Uncorrelated Channels

In this section, we investigate a special case of uncorrelated channels which allows a closed-form solution to the optimal user loading ratio at low SNR.

### Iv-a Asymptotic Achievable Rate

Under the special case of uncorrelated channels, i.e., in (4), the value of in (13) can be further simplified. It reduces to the unique positive solution of the following function

 ξ=(1+ξ)ρ(1+ξ)+β. (24)

Note that the solution to is a function of and . By solving the quadratic equality in (24), the positive one of its solutions, referred to as , equals

 ξUC(β,ρ)=12[√(1−β)2ρ2+2(1+β)ρ+1+1−βρ−1]. (25)

From the definition of in (14), we have with for an uncorrelated massive MIMO channel. Thus, the asymptotic SIQNR per user in (12) can be simplified as

From the conclusions in Remark 1, it indicates that the optimal regularization parameter is independent of the spatial correlation. Thus, remains the same as expressed in (16), which maximizes in (26). One typical scenario for massive MIMO is widely acknowledged as the low SNR case. We focus on this special case with . Substituting in (16) into in (25), it yields

 ξ∗UC(β,ρ∗) =12  ⎷(1−β)2[1(1−ρDA)γ0+ρDA1−ρDA]2β2+2(1+β)[1(1−ρDA)γ0+ρDA1−ρDA]β+1+1−β2[1(1−ρDA)γ0+ρDA1−ρDA]β−12 ≈(1−ρDA)γ0β+o(γ0), (27)

where (27) applies the Taylor’s expansion for small , i.e., , and uses and for small . is an ignorable higher-order term w.r.t . Replacing in (20) with and substituting (27), the achievable rate of each user at low SNR can be rewritten as:

We observe that a smaller value of , or provides higher . It implies that high-resolution ADCs, DACs, or a small user loading ratio provides high achievable rate for each user.

Then, we characterize the rate loss caused by the finite-resolution ADCs for a fixed quantization resolution of DACs.

###### Lemma 2.

Under the uncorrelated massive MIMO channel, the normalized rate loss per energy caused by finite-resolution ADCs with fixed-bit DACs can be characterized at low SNR as

###### Proof.

Given the ideal infinite-resolution ADCs, the benchmark asymptotic rate can be evaluated from (28) by letting . It gives

 ~RUC=log2β+(1−ρDA)γ0β. (30)

Then, the normalized rate loss per energy can be characterized by subtracting (28) from (30) as follows

where we use the fact that for . ∎

Remark 2. Given , is approximately linearly affected by and at low SNR. Particularly, for as set in the following simulations, we have with and according to Table I. This implies that, given 1-bit DACs at the BS, the normalized rate loss in (29) caused by 3-bit ADCs at the user side is relatively negligible.

Under the assumption of low SNR, we are able to obtain the optimal user loading ratio in closed form by maximizing the sum rate per antenna, given in Lemma 3.

###### Lemma 3.

For an uncorrelated massive MIMO system, the problem in (22) has a closed-form solution for asymptotically low SNR as

###### Proof.

Utilizing (21), we first rewrite the sum rate per antenna without spatial correlation as

where is expressed in (27). It directly yields that

where (34) comes from according to (27). Letting , we get

where we use the equality for small . Substituting (27) into (35) and dropping since , we arrive at a quadratic equation w.r.t. , i.e, , which has only one positive solution given in (32).

Remark 3. On one hand, is an increasing function w.r.t . More users should be loaded at higher SNR in order to achieve the maximum sum rate per antenna, while it increases the total pilot overhead. On the other hand, increases with while it decreases with . It implies that the optimal user loading ratio increases as the ADC resolution decreases at the users or the DAC resolution increases at the BS.

## V Simulation Results

In this section, we evaluate the asymptotic SIQNR, , the asymptotic achievable rate of each user, , and the asymptotic sum rate per antenna, , via simulations. Optimal regularization parameter is verified by numerical methods. Moreover, the closed-form solution of in (32) is also tested at low SNRs. Consider a well-known exponential Toeplitz correlation model as follows

 rij=ν|i−j|, (36)

where denotes the correlation coefficient and represents the entry of at the th row and th column. We use for simulation in the following. Generally, small correlation coefficient is representative in a rich scattering environment, e.g., for indoor scenarios. While large values of, e.g., represents the scenario where the BS locates in an unobstructed environment.

We calculate in (12) and test our derived in (16) using the values of , , and calculated in Appendix D. Fig. 2 displays the asymptotic SIQNR with various correlation coefficients. Obviously, first increases but then decreases with increasing . The pentagrams mark the maximum with , which verifies obtained from (16). It can be observed that is independent of as analyzed in Remark 1, which implies that the optimal RZF precoding design is independent of the antenna spatial correlation. Fig. 3 shows the asymptotic SIQNR with various ADC resolutions, i.e., bits. is marked by pentagrams and is independent of the ADC resolution as indicated by (16). Fig. 4 displays the asymptotic SIQNR with varying DAC resolutions. According to (16), marked by pentagrams are , and for , and , respectively. It can be observed that increases with decreasing as discussed in Remark 1. Fig. 5 and Fig. 6 show the asymptotic SIQNR with various and . In Fig. 5, , and for , and , respectively. We verify that proportionally increases with . In Fig. 6, , and for , and dB, which verifies that decreases with .

Fig. 7 compares the achievable rate per user with the optimal and other three typical values, i.e., leading to a conventional RZF precoder, leading to a ZF precoder, and leading to an MRC precoder. The performance gap between the optimal RZF and the conventional RZF is generally not that significant, but it increases with SNR. This is because the effect of quantization noise becomes to dominate the performance at high SNR. Thus, the regularization parameter is better to be optimized according to quantization distortion factors. In addition, the optimal RZF precoding provides better rate performance than ZF precoding with ranging from to dB. As for MRC precoding, it achieves almost the same rate as the optimal RZF at low SNRs since the channel noise significantly overwhelms the inter-user interference under this condition. However, the multiuser interference becomes dominating as SNR increases and thus the achievable rate with MRC precoder becomes much lower than the optimal RZF at high SNRs.

Fig. 8 compares the bit error rate (BER) of various precoding schemes using QPSK modulation. Error floor exists due to the use of 1-bit DACs and 3-bit ADCs. Compared to conventional RZF and ZF precoders, the optimal RZF precoder improves the BER performance more significantly with larger . It implies that the optimal RZF precoding is more effective for spatially correlated massive MIMO channels using low-resolution DACs and finite-resolution ADCs. Note that the BER values are high because we consider a heavily loaded scenario with and no channel coding is employed. This allows us to better focus on evaluating the performance of the proposed optimal RZF precoder.

Fig. 9 shows the achievable rate of each user affected by finite-resolution ADCs. Solid lines correspond to our derived asymptotic expression in (20) and it matches well with the simulation results. Note that converges with increasing SNR since only one-bit information can be transmitted by each antenna with 1-bit DACs. Besides, higher ADC resolution provides larger and 5-bit ADCs achieve almost the same rate as infinite-resolution ones. Note that at low SNRs, the rate loss caused by ADC quantization is marginal.

Fig. 10 compares the sum rate per antenna of the optimal and conventional RZF, ZF, and MRC precodings. The SNR is  dB. The optimal RZF precoding provides higher than the other three precoders with ranging from 0 to 1, which implies that it is necessary to optimize the regularization paramater . With increasing from 0 to 1, first increases since more users are served but then decreases due to higher pilot overhead. Obviously, an optimal exists by maximizing . Compared to the conventional RZF precoder, the optimal RZF precoder significantly improves the sum rate for . This is because the optimal differs more with the conventional for larger . While for even larger , we have for and thus the performance gap tends small.

Although the closed-form in (32) is obtained for a special case without spatial correlation, Fig. 11 tests it for a massive MIMO channel with an exponential Toeplitz correlation matrix as in (36). Comparing the rates with the derived and exact , we observe that the derived is still precise at low SNR with less significant correlations. It implies that (32) can be applied to massive MIMO systems where the correlation is generally not severe.

We show the asymptotic rate performance without spatial correlation in the following. Fig. 12 shows the sum rate per antenna under an uncorrelated channel. The pentagram marks the maximum with . Specifically, equals , and , for , and , respectively. We find that increases slightly when ADC resolution decreases as discussed in Remark 3. Fig. 13 compares the performance of derived from (32) with exact value by numerical methods. It can be observed that with our derived is extremely close to that with exact one at low SNR.

## Vi Conclusion

In this paper, we investigate a downlink multiuser massive MIMO network and study the performance of low-resolution DACs equipped at BS and finite-resolution ADCs at use sides. RZF precoding is used and transmit-side spatial correlation is considered. By maximizing the asymptotic SIQNR, we derive the optimal regularization parameter which is found independent of the channel correlation and ADC quantization. It implies that the optimal RZF precoding can be conducted at the BS without any knowledge of the ADC resolution at the user sides. Moreover, by maximizing the sum rate per antenna, a closed-form solution of the optimal user loading ratio is obtained at low SNR. Although the obtained is derived under a special case without correlation, it is verified also valid for slightly correlated channels. Some extensions of this study, e.g., considering imperfect CSI and multi-antenna users with low-resolution ADCs, are of interest in future work.

## Appendix A

The following lemma is recalled from [25, Lemma 1] for completeness.

###### Lemma 4.

In the large system limit, the following convergence of terms in (9) can be stated as follows:

1) Let denote the channel matrix removing the th row and we have

 hTk(HHkHk+αIN)−1h∗ka.s.−−→ξ, (37)

where follows the same definition in (13).

2) Moreover,

 hTk(HHkHk+αIN)−1HHkHk (38)

where follows the same definition in (14).

3) As for constant , it converges to

 c2a.s.−−→−P(1+ξ)2β∂ξ∂ρ. (39)

## Appendix B

###### Lemma 5.

Considering a typical positive definite Hermitian correlation matrix with a Toeplitz structure, converges to

 diag(PPH)a.s.−−→PNIN. (40)
###### Proof.

Using (10), we have

 diag(PPH) =c2diag((HHH+αIN)−1HHH(HHH+αIN)−1) =c2diag((HHH+αIN)−1(HHH+αIN−αIN)(HHH+αIN)−1) =c2diag((HHH+αIN)−1−α(HHH+αIN)−2) a.s.−−→c2diag((MR+αIN)−1−α(MR+αIN)−2) (41) ≜dIN, (42)

where (41) uses (4) and the fact that

due to the Central Limit Theorem. According to the Continuous Mapping Theorem

[39], the convergence preserves for continuous matrix functions like, and , of nonsingular matrices with positive definite Hermitian matrix and . The matrix within the function in (41) is still a Toeplitz matrix since the inverse of a positive definite Toeplitz matrix is asymptotically still Toeplitz [39, Theorem 4.3]

. Thus, it becomes a scaled identity matrix after taking the diagonal entries as in (

42), where is defined as a constant (absorbing ) which can thus be determined by checking the power constraint of the term within the brackets, i.e., , yielding in (40). ∎

## Appendix C Proof of Theorem 1

The proof of Theorem 1 applies Lemmas 4 and 5. The value of (9) is calculated term-by-term. First concerning the numerator in (9), the energy of the desired signal follows

 ∣∣hTkpk∣∣2 =c2∣∣∣hTk(HHH+αIN)−1h∗k∣∣∣2 (43) (44) a.s.−−→−Pξ2β∂ξ∂ρ, (45)

where (43) is obtained by substituting (10), (44) uses the matrix inversion lemma that

 (HHH+αIN)−1h∗k= (HHkHk+αIN+h∗khTk)−1h∗k = = (HHkHk+αIN)−1h∗k1+hTk(HHkHk+αIN)−1h∗k, (46)

and (45) utilizes 1) and 3) in Lemma 4, applying that convergence preserves for the continuous function according to the Continuous Mapping Theorem [39]. According to (13) and (14), we have

 ∂ξ∂ρ=−(1+ξ)2E121−βE22. (47)

Substituting (47) into (45), the desired signal power converges to

 ∣∣hTkpk∣∣2a.s.−−→Pξ2(1−βE22)(1+ξ)2βE12. (48)

Then, we consider the multiuser interference in the denominator of (9). It follows that

 ∑j≠k∣∣hTkpj∣∣2