Downlink multi-user algorithms for millimeter-wave wideband linear arrays on PD-NOMA-based squint steering beams

This article investigates using a phased linear antenna array instead of the planar array to circumvent the problem that two frequency squint steering main beams cannot cover any two beam directions simultaneously. First, we approximate the donut-shaped main beam of the linear array by means of multiple pencil-shaped main beams of a virtual planar array for matching the steering main beam of the linear array with the multi-path sparse scattering channel model mathematically and give a method for calculating the number of antenna elements of the virtual array. Second, we cope with possible inter-user interference on a single squint main beam of the linear array in some scenarios by means of the power-domain non-orthogonal multiple access (PD-NOMA) technique, making it possible to support communication with two users on a single squint main beam at the base station (BS) side. The feasible domain of PD-NOMA is given when a single antenna is used for both the BS and the user end, assuming a two-user successive interference cancellation (SIC) decoding power ratio limit. Third, three algorithms are given for serving multi-user at the BS via squint beams of the linear array. Finally, numerical results show that the second proposed algorithm supporting PD-NOMA pairing within a single donut-shaped squint main beam significantly increases the number of simultaneous users served within a single cellular system.


Introduction
Wideband millimeter-wave (mmWave) communications have received increasing attention in recent years, and the beam squint problem caused by phase differences between frequencies in the radio frequency (RF) analog domain induced by the utilization of nonideal wideband phase shifters is one of the important problems encountered [1][2][3]. While the use of frequency squint multi-beam to serve multi-user within a cellular cell has been studied [3], the restricted horizontal azimuth of the squint steering beams introduced by employing the planar antenna array limits its application in cellular communications. In contrast, the steering beam of a linear array arranged in the z-axis can cover all horizontal azimuths in the x-y plane while changing the elevation angle, thus avoiding the problem of restricted horizontal azimuth of the squint beam of a planar array. However, the hollow conical main beam of the linear array raises two new issues: (i) the matching of the donut-shaped steering main beam and the sparse scattering channel model and (ii) the handling of the additional interference from such main beams over the corresponding pencil-shaped beams of the planar array.
As for the former issue, the mmWave sparse scattering Saleh-Valenzuela channel model represents one channel path according to the form of the product of beam steering vectors of antenna arrays at both the transmitting and receiving ends [4,5]. For linear arrays, the array factor and channel path are generally described by one angular parameter in a 2-D scenario [6], which is more convenient and concise to describe the nature of the array, but the main direction of the steer beam of a linear array is donut-shaped in 3-D space, so there may exist scattering paths at different horizontal azimuths, where the angle-of-departure (AoD) of a particular channel path is within the angular range of the transmitting donut-shaped main beam. Therefore, it is not complete to describe the channel information of a linear array in the 3-D scenario with only one angular information, which raises the problem of matching the steer beam and the channel path. To the extent of the authors' knowledge, this issue has not received sufficient attention so far.
On the latter problem, we draw on the joint spatial division and multiplexing (JSDM) grouping approach [7,8] and the power-domain non-orthogonal multiple access (PD-NOMA) scheme [9][10][11] to deal with the multi-user interference between steering beams and on a donut-shaped main beam of the linear array at the base station (BS) side, respectively. On the one hand, the work in [7] assumes that the BS is high placed, leading to scatterers mainly at the user end, and therefore, the majority path of a user has a relatively close AoD, so the one-ring model is adopted to model the channel for each user. Afterwards, the users are divided into several orthogonal or near-orthogonal groups based on the channel state information, and then the interference between users is processed within the groups. Whereas the work in [8] assumes that multi-user can share a common scatterer in the mmWave cellular communication scenario, the user grouping is then performed by exploiting the channel covariance information, which can also be interpreted as a grouping by the AoD of the channel paths at the BS side.
On the other hand, the downlink PD-NOMA scheme multiplexes time-frequency blocks and mitigates interference by means of signal superposition at the BS end and successive interference cancellation (SIC) at some receivers [9]. In the mmWave cellular scenario, it is currently evolving mainly from NOMA pairing between near and far users within a single main beam to the more novel NOMA pairing among multi-beam [12]. The work in [13] investigates the beam design and power allocation problems in the two-user paired PD-NOMA scheme, without considering the effect of the donut-shaped main beam of the linear array. As the AoD of the scattering path of multi-user is potentially on a single donut-shaped main beam of the linear array at the transmitter end, in this paper, we use the PD-NOMA to enable a single squint beam of a wideband phased array at the BS end to communicate with two users simultaneously in some scenarios, further increasing the number of users served concurrently in the system. To briefly illustrate the frequency squint property of the wideband array, we only employ the steering beam and do not consider the possibility of more complex frequency squint precoding schemes, which also avoids the effect of additional non-convexity introduced into the precoding design optimization problem by the constant modulus constraint of phase shifters in the array. In this case, the problem of finding suitable users for pairing and allocating the transmit power for the two paired users needs to be addressed in order to implement the PD-NOMA opportunity communication.
The remainder of this paper is organized as follows. In Section 2, we address the problem of matching the donut-shaped main beam of the linear array with the sparse scattering channel model in a 3-D scenario, as well as providing both the channel model and the system model. In Section 3, we first discuss NOMA user pairing strategies for a steering beam, after which three beam squint algorithms for mmWave multi-user communications are presented. In Section 4, simulation results show the performance of the proposed algorithms, before concluding the paper in Section 5.
Notation: Boldface uppercase letters, boldface lowercase letters, and lowercase letters are used to denote matrices, vectors, and scalars, respectively. The superscripts (·) T and (·) H denote the transpose and conjugate transpose, respectively; · F represents the Frobenius norm of a matrix; · denotes the rounding function; | · | is the amplitude of a complex value.

Matching of the scattering path and the steering main beam
It is assumed that the linear array is placed in the z-axis, thus having a unique elevation angle in the main direction of its steering beam, which facilitates the formulation of channel modeling and algorithm design when using the steering beam. As shown in Fig. 1, the main direction of the steering beam of the linear array is donut-shaped and can be divided into multiple pieces by vertical slicing, each piece being approximated by one steering main beam of the same virtual planar array in the x-y plane. Moreover, due to Fig. 1 The linear array is arranged on the z-axis, and its donut-shaped steered main beam is approximated by several steering beams of a virtual planar array on the x-y plane the limitations of scatterer size in the channel and beamwidth at the transmitter side, the scattering path of the mmWave channel in practical cellular communication scenarios is more directional and can generally be described by the array response of a planar array. Accordingly, the channel path can be specified by the array response of the virtual planar array corresponding to the linear array at the BS side, which is thus mathematically coherent from the transmitter, channel, and receiver side. In addition, it is important to ensure that the power in one direction of the donut-shaped main beam is not excessive by the angular separation of the horizontal azimuth φ of the two adjacent slices. Furthermore, it is not necessary for this virtual ring main beam to be complete. If there is no scattering path in some horizontal azimuth range, the slices in the corresponding range can be removed to give a clearer picture of the connection between the transmitter beam and the channel scattering path. If there are multiple AoDs within the half-power beamwidth (HPBW) of the virtual planar array corresponding to the same donut-shaped main beam, the relationship to these channel paths can be represented by a slice in the main beam instead.

Number of elements in the virtual planar array
As shown in Fig. 1, the normalized array factor for a linear array of N uniformly spaced array elements placed on the z-axis and with the origin set as the reference point can be written as where where k is the wave number, d is the array element spacing, θ is the elevation angle, and β is the phase difference of the emitted signals from adjacent array elements. After adjusting the reference point to the array midpoint [6], the array factor can be approximated as In (2), by setting ψ = 0, we have where θ m is the elevation angle corresponding to any horizontal azimuth angle in the donut-shaped steering main beam direction of the linear array, which can be kept to a uniform value of θ ring by selecting a suitable coordinate system. Since it is possible for the donut-shaped main beam to be irradiating several scatterers in the physical channel, the connection between the transmitting main beam and the channel path can be described by a set of steering main beams of a planar array with disjoint HPBW ranges when assuming that the energy radiated by the side lobes in the transmitting beam at the receiving end is highly attenuated and negligible due to atmospheric absorption of mmWave. Thus, by approximating the part of the array factor of the linear array in the first quadrant through the combination of the normalized array factors of a virtual planar array in the x-y plane in Fig. 1, and more precisely, by approximating the main lobe of the steering beam of the linear array by the main lobes of several steering beams of the planar array and neglecting the effect of the side lobes, we have The half-power direction θ h of the main beam of the linear array can be determined by Therefore, the HPBW of the linear array can be approximated as On the other hand, the HPBW of a planar array in the x-y plane can be expressed as [6] where x0 and y0 are the HPBWs of linear arrays equally spaced on the x-and y-axes with the number of array elements M x and M y , respectively, with the beam direction toward the z-axis, with in which c is the velocity of electromagnetic wave. Similarly, replacing M x with M y in (10) gives y0 . When a square planar array is employed, such that M = M x = M y , (8) and (9) can be further simplified to To approximate the steering main beam of the linear array by steering main beams of the planar array, let (7) and (11) have equal HPBWs, which gives Since the HPBW is not perfectly symmetrical on both sides in the main beam direction, the relationship in (13) can be written exactly as Substituting (4) where (a) is based on the assumption that the two terms in square brackets in the denominator are the half-power direction of the main beam, (b) is based on the definition of the HPBW, and (c) is due to substituting (11). If (10) is substituted into (c) again, we get M, i.e., deriving (15) from the top down and from the bottom up, for the assumption to hold sufficient necessary for the half-power direction θ ± hi and satisfies θ − hi > θ + hi . Therefore, the half-power direction of the frequency squint main beam in elevation can also be noted as Alternatively, the approximate expression of the above equation can also be obtained by substituting β = kd cos θ m into (13) and by a simple derivation, namely Since the number of elements on each side of a planar antenna array is an integer, simply rounding to the nearest whole number, we have . Furthermore, since the definition of x0 in (c) of (15) is the HPBW of the linear array in the vertical direction, which is independent of the actual beam direction θ mi , it is related to the number of antenna elements. Therefore, it can be assumed that changing the beam direction will cause a change in the number of elements of the virtual planar array. Due to the small angular spread of the mmWave band and the narrow HPBW due to the high number of elements in the array, we use the same virtual planar array to approximate paths in the HPBW range of the same scattering cluster for ease of handling.
Next, we study the relationship between the direction of the frequency squint main beam θ mi and the parameter M of the virtual planar antenna array once the main beam direction θ m at frequency f c has been determined. Since a linear array on the z-axis is considered, there is As the main direction of the beam at frequency f c is θ m , substituting into (23) gives β = π cos θ m ; again substituting into (23) gives By substituting (24) into (17), we get By fixing the array spacing of the virtual planar array to half the wavelength corresponding to the central frequency f c and rewriting (25), we get

System model
Following the multi-path scattering Saleh-Valenzuela channel model [5], the channel of a particular user can generally be represented as where N c is the number of all clusters from the BS to that user, N pc is the number of paths within the c th cluster, and a r (φ r cl , θ r cl ) and a t (φ t cl , θ t cl ) are the steering beam vectors represented by array responses at the receiving and transmitting ends for the l th path in the c th channel cluster, respectively.
For convenience, due to the small angular spread of the beam in the actual mmWave sparse scattering environment, it is assumed that the paths within each cluster can be represented by the array response of the same virtual planar array; if the user has multiple scattering clusters on the same donut-shaped main beam, we still use the array response of the same virtual planar array to approximate the paths in multiple scattering clusters; when the difference in elevation angle between clusters is large, corresponding to several steering donut-shaped beams at the transmitting side, the difference in size of the virtual planar array between clusters is more pronounced, and we assume that the donut-shaped main beam affects a range within the HPBW of the transmitting elevation angle θ t ring in the real channel environment. According to the approximation of the array factor in (5), the actual channel caused by a particular steering beam of the linear array can then be equivalently expressed as whereN c andN pc are the number of clusters and paths within the HPBW range of θ t ring , respectively; the elevation angleθ t cl in the AoD of each path in the channel needs to be within the HPBW range of θ t ring at the transmitter side, i.e., satisfying |θ t ring − θ t cl | ≤ h f i /2,θ t cl ∈ {θ t cl };ã t is the vector form of the virtual planar array response corresponding to a certain steering beam of the linear array.
When the user at the receiving end uses a square antenna array with M r elements on each side and placed in the x-y plane, its array response at frequency f s with the main beam orientated to (φ r s , θ r s ) is And the array response of the virtual planar array at the transmitter side can be written asã where M t is the number of elements per side of the virtual array obtained from (18), f s is the operating frequency, and (φ t s , θ t s ) is the virtual main beam direction. Consequently, the equivalent baseband signal at the user receiver side can be represented as where P is the transmitted power of the individual steering beams. Hence, the rate for this user can be expressed as where (d) is obtained by representing the channel according to (28) as the actual channel generated by a single donut-shaped transmitting main beam in the scattering environment and assuming that the path l with the highest received energy in the actual channel is taken through the steering beam at the receiver end; (e) is assumed that the quantization resolution of the phase shifter of the wideband analog array at the receiver side is high enough to satisfy φ r i = φ r l and θ r i = θ r l . Figure 4 shows the diagram of a two-user paired NOMA scheme in a downlink multi-user scenario when the BS side is equipped with a linear array and each user end with a planar array, where the linear array is arranged in the z-axis and one of its transmitting steering beam has two scatterers in its donut-shaped main beam direction, and the corresponding two channel clusters are scattered and directed toward the far user 1 and the near user 2, respectively. The average received energy in the corresponding scattering cluster is smaller for user 1 due to the longer path experienced by the scattering path and the fact that the energy of the longer path is more absorbed in the atmosphere. Due to the use of the phased linear array at the transmitter side, there is no particularly well-handled solution for inter-user interference on the donut-shaped main beam at the transmitter side [14][15][16]. On the one hand, we assume that a grouping strategy has been used to serve these users with serious interference using orthogonal resources. On the other hand, it is also possible to pair users with similar departure angles on some of the beams to serve more users through PD-NOMA in the face of unavoidable interference. We classify NOMA pairing strategies into following three categories and discuss their feasibility. Among them, the first strategy is employed in both Algorithms 2 and 3, the second strategy is only used in Algorithm 3, and the last strategy is not included in our algorithms and we have explained its reason. Furthermore, we provide Algorithm 1 for steering the squint beam without any NOMA scheme as a benchmark for comparison.

User pairing strategies and squint multi-beam algorithms
• NOMA pairing within the direction of a donut-shaped main beam This type of NOMA pairing corresponds to Algorithm 2, which selects two users whose channels contain paths with similar departure elevation angles, and the energy difference between these two paths is large enough to consider the PD-NOMA scheme to serve user 2 as well while guaranteeing the minimum required rate for user 1.
• NOMA pairing in the direction of the main beam and its first side lobe We attempted to utilize the side lobe direction of a steering beam for pairing to break the distance limit between the two paired users and the BS side in the conventional PD-NOMA scheme. The amplitude of the first side lobe in its dominant direction is −13.465 dB, with a normalized amplitude of about 0.2 (0.212 is more accurate according to equation (6)(7)(8)(9)(10)(11)(12)(13)(14)(15)(16)(17) in [6]). It is possible to consider user 2 as the main beam direction user, while there is a scattering path of user 1 in the elevation direction of the donut-shaped first side lobe of the steering beam. In this case, the energy of the corresponding path of user 1 is not necessarily weaker, but since it is the energy of the equivalent channel in the direction of the first side lobe that is significantly weaker than that of user 2, it is also possible to serve simultaneously these two users with significantly different departure elevation angles of the channel path through the PD-NOMA scheme. The corresponding power allocation coefficients β 1 and β 2 are then obtained by applying 0.04|h 1 | 2 to replace the original |h 1 | 2 as the channel energy for user 1, and the corresponding SIC decoding feasible domain conditions should reformulate (66) and (72) in Appendix 2 as It is found by simulation that the channel coefficients for both user 1 and user 2 are too large in this case, and according to the parameter settings in Fig. 8, approximately  The BS rejects the enhanced beam service for this SU, reports and feeds back the request result 13: end for |h 1 | 2 > 4 needs to be satisfied. This magnitude requirement is too large for a single-antenna scenario. For scenarios using array antennas, there is some feasibility due to the array gain at the transceiver end, which will be demonstrated in Algorithm 3 and verified by simulation experiments. • NOMA pairing in the direction of the main beam and its grating lobe Pairing between the main beam and grating lobe directions is also possible; however, the presence of a grating lobe would reduce the transmit power in the main beam direction while maintaining a uniform transmit power for each beam. This pairing is not adopted for the time being, since in this paper the squint beam is applied and there is no need to provide multi-beam with grating lobes. It is possible to circumvent this by limiting the maximum elevation angle of the main beam and the maximum working frequency of the squint beam, according to the conditions given in the forepart of Appendix 1 for a wideband array not to create a grating lobe.

Benchmark: downlink frequency squint multi-beam algorithm without PD-NOMA
Algorithm 1 draws on the idea of grouping in [8] and does not use the scheme of the PD-NOMA. For convenience, the user operating at the central frequency f c is referred to as the primary user (PU) and the users operating at other frequencies are referred to as secondary users (SUs).
In step 1, the BS updates the array response of its steering beam at frequency f c according to the elevation angle θ t new in the AoD of the strongest path of the PU. According to Pan  (30) and (31), the steering beam of the linear array is approximated by the steering beams of a virtual planar array. In step 2, the PU calculates the receiving array response from (29) based on the new arrival angle (θ r new , φ r new ) and operating at frequency f c to generate the received beam. In step 5, the elevation angle space at the transmitter end is divided into angular regions by the HPBW h at the central frequency of the PU, with each region supporting at most one donut-shaped squint main beam within the angular range, and each steering beam serving at most one user at a time. It is obtained from (6) and (7)  IsGroupActivated(·) indicates whether a particular transmit elevation interval group is occupied by true or false, and the default value is false. In step 6, since the BS side uses a wideband phased array and the phase shifters therein are unquantized ideal ones, the operating frequency of the desired squint beam f s can be calculated from (54).
In step 7, the BS reuses the PU's array response at frequency f s , marks this transmit elevation interval group as occupied, and feeds back f s to that user. IsUserServed(·) indicates whether a particular SU is served by true or false, and the default value is false.
In step 8, the SU calculates the array response based on the new arriving angle (θ r s , φ r s ) and f s from (29) and works at frequency f s to produce the reception beam.

The enhanced beam squint algorithm supporting PD-NOMA in the donut-shaped main beam
Algorithm 2 improves on Algorithm 1 by enabling each donut-shaped squint main beam at the transmitting side to serve two users simultaneously via the PD-NOMA.
The work in [17] gives the achievable rate of NOMA users when a single antenna is employed at both the transmitter and receiver, while in this paper, since the BS uses a linear array and each user uses a planar array, the achievable rate of users in the two-user NOMA pairing scheme can be written as The BS rejects the enhanced beam service for this SU, reports and feeds back the request result 19: end for Note that NOMA pairing requires that user 2 be able to decode user 1's signal through the SIC, so the maximum transmit power coefficient of user 2 is limited to β max 2 in step 1, and after selecting user 2, user 1 is chosen exhaustively until a suitable size of β 2 is picked.
In step 1, the BS divides the departure elevation angle into intervals according to θ t s ∈ θ t new − i h /2, θ t new + i h /2 , and in each interval, if there is a user path, then choose Note that since the interval is divided only by the HPBW of the elevation angle corresponding to the frequency f c , it is possible that when f s > f c could result in the departure angle of the strongest path of the selected paired user being outside the HPBW range where the actual squint main beam is located, thus reducing the rate performance of the user in the interval. This problem is potentially inherent in both Algorithms 2 and 3, and a more reasonable way of grouping by angular space division is subject to further investigation. However, it is ignored for the moment in this paper as it is more concerned with increasing the number of simultaneous users served through PD-NOMA groups.
In step 4, the selection of the group G i with i = 0 is due to the assumption that the PD-NOMA scheme is used only in squint beams.
In step 5, the departure angles of the strongest paths of users 1 and 2 in the NOMA pair are denoted as θ t s1 , φ t s1 and θ t s2 , φ t s2 , respectively, and the arrival angles of the corresponding paths of the two users are denoted as θ r s1 , φ r s1 and θ r s2 , φ r s2 , respectively. In step 6, the BS takes θ t s1 + θ t s2 /2 as the transmit elevation angle θ t s of the steering beam in this interval and calculates its operating frequency from f s = f c cos θ t new / cos θ t s . In step 7, the BS allocates transmitted power coefficients β 1 and β 2 within the beam to the users in this NOMA pair according to (38) and (39).
In step 9, the BS reuses the array response of the PU at frequency f s and feedback f s to the users within this NOMA pair. IsUserPaired(·) indicates whether a particular SU is served in a NOMA group by true or false, and the default value is false.
In step 10, users 1 and 2 within the NOMA pair calculate the array response based on the new arrival angles θ r s1 , φ r s1 and θ r s2 , φ r s2 as well as the new squint frequency f s from (29) and work on the frequency f s to generate the receiving beam, respectively. After that, it jumps out of the inner loop at the BS side.

The enhanced beam squint algorithm supporting NOMA pairing of the main beam direction with the first side lobe direction
Algorithm 3 improves on Algorithm 2 by supporting the PD-NOMA opportunity communication for pairing between a user in the direction of the squint main beam and a user in the direction of its first side lobe. The general idea is to select a group with only one active user and find a suitable unserved user in the first side lobe direction of the steering beam in the group for pairing as NOMA user 1.
To avoid the concept of defining the HPBW in the direction of the first side lobe, we assume that the direction of the first side lobe points exactly to the strongest path of NOMA user 1. By exhaustively enumerating the unserved users as NOMA user 1, the operating frequency f i and the main direction θ mi of the squint beam corresponding to the direction of its strongest path can be obtained from the direction θ s1 of NOMA user 1, assuming that NOMA user 1 has been selected and from equation (6-16a) in [6], the first side lobe direction of this steering beam is given by , feeds back f i to users 1 and 2 in the NOMA pair 8: The BS rejects the enhanced beam service for this SU, reports and feeds back the request result 16: end for Substituting (54) into (40) gives the relationship between the direction of the first side lobe and the operating frequency, as In step 2, the elevation angle in the AoD of this SU's strongest path as the first side lobe direction θ ± s1 , the corresponding squint frequency f i of the new squint main beam is obtained from (41). To guarantee the validity of the HPBW range and power distribution coefficients in the simulation, set f i = cos θ m + 3 N f c / cos θ + s1 . After the operating frequency f i has been determined, the half-power beam angle range θ + hi , θ − hi for the main direction of this squint beam can be obtained, and within this angle range, the user who occupies a steering beam alone is matched as user 2 in the NOMA pair, and then the array response of the PU's transmitted beam is multiplexed at frequency f i to serve this NOMA pair, where there is According to (21) the half-power direction of this main beam is alternatively it can be written as In step 5, treat the selected SU as user 2 in the first side lobe NOMA pair and calculate the power allocation coefficients β 1 and β 2 from (45) and (46). When calculating the power allocation coefficients, the equivalent channel energy of user 1 should be multiplied by 0.04, since it is in the direction of the first side lobe of the transmitting beam. Therefore, the power allocation coefficients in the NOMA pair should vary from (38) and (39) to Since the user in the first side lobe direction is selected as user 1 in the NOMA pair before selecting user 2 in Algorithm 3, there is no need to adjust the frequency of the squint beam in the process of selecting NOMA user 2, thus avoiding the situation where the operating frequency determined by user 1 in the NOMA pair causes the strongest path of user 2 to be outside the HBPW of the squint main beam, resulting from the selection of NOMA user 2 first.
In step 8, users 1 and 2 within the NOMA pair calculate the array response based on the new arrival angles θ r s1 , φ r s1 and θ r s2 , φ r s2 as well as the new squint frequency f i from (29) and work on the frequency f i to generate the receiving beam, respectively.
In step 9, the BS marks users 1 and 2 as users in a NOMA pair and jumps out of the inner loop.

Results and discussion
The simulations in this section are based on a uni-cell scenario with one BS and multiuser, where the BS side consists of a linear array of N t antenna elements, the total number of users is N u and each user end uses a square planar array with M r elements on each side, and all transceivers are wideband phased antenna arrays with central frequency f c of 45 GHz and element spacing of λ c /2. It is assumed that the linear array of the transmitting BS is placed in the z-axis, while the planar array of the receiving end is placed in the x-y plane of the user's respective coordinate system and that both the elevation angle θ mi and the azimuth angle φ m (f i ) of the main beam at both the transceiver ends are within [ 0 • , 90 • ]. The elevation angle θ m of a PU transmitted beam is limited to [ 30 • , 60 • ] for more efficient use of the squint beam, and the operating frequency interval of the squint beam is limited to [ f c /2, 3f c /2] for avoiding the appearance of grating lobes and an excessive number of elements in the virtual planar array. Similar to [3], we use the extended Saleh-Valenzuela channel model [4] for the mmWave band to model the user's channel, with each user operating at its own frequency and its channel consisting of several clusters with a number of paths in each cluster, each of which is described by multiplying the product of the array response of the transmitting and receiving steering beams by the channel coefficient. Specifically, the difference in the scattering environment between users is represented by the number of clusters per user channel which is randomly chosen between 1 and 8, and the number of paths per cluster which is randomly picked between 1 and 10. The mean cluster angle of the scattering clusters of each user follows a uniform distribution, and the arrival and departure angles of the paths within the clusters each follows a Laplacian distribution with a standard deviation of this cluster angle spread, which is held for all clusters at 7.5 • . The signal and noise both independently obey a Gaussian distribution with mean 0 and standard deviation of 1, i.e. , CN (0, 1), and assume that the scenario is noise limited such that P Note that the maximum number of steering beams at the BS side is constrained by its power constraint in practice, which is not limited in the simulation; the channel of the i th cluster of a user satisfies E H i 2 F = N t N r , where N t is the actual number of antennas in the transmitting linear array, N r = M 2 r , α il is the complex gain coefficient of the l th path within its i th cluster, and its value obeys the i.i.d. complex Gaussian distribution is the average power of the internal path of this cluster. And to simplify the design, only the direction of the maximum energy path of each user is actually selected by the transceiver for the direction of the steering beam, which avoids the problem of dimension mismatch of the transmit steering beam between multi-cluster caused by inconsistent virtual planar array sizes by assuming that only one cluster in the channel needs to be handled. Finally, all results are reported for averaged over 5000 random channel realizations. Unless otherwise mentioned, these parameters are used for the experiments. Figure 5 compares the number of simultaneous users served when the BS is equipped with a linear array and a planar array, respectively, using benchmark Algorithm 1. For fairness, when the number of antennas of the linear array is determined, the number of antennas of the planar array is selected in such a way that the steering beams of both are pointed at the same path and the HPBWs in the vertical direction of the two squint main beams are basically the same. The number of users served of the linear array is significantly better than that of the planar array due to the constraint of the horizontal main direction of the squint steering beam of the planar array, and the more users in the system, the more significant it is. Figures 6 and 7 compare the number of simultaneous users served and the sum rate performance of the system for the three proposed algorithms in two scenarios with the number of linear array antennas at the BS side of 32 and 64 while keeping N r = 64 for all users and with t 1 taken at 2 bits/s/Hz. As the transmit elevation angle is divided into groups by the HPBW of the steering beam of the array at the BS side, the HPBW is larger when the number of antennas is smaller, and the angle interval per group is larger and the number of groups divided is smaller. Specifically, the number of groups varies depending on the direction of the main user's transmit beam, with about 15 intervals at N t = 32 and more than 40 intervals at N t = 64. Therefore, when the number of antennas at the BS side is smaller, it is easier to have serviceable users within a group for the same number of users in the system. Consequently, with the limitation on the number of users that can be served simultaneously in each group, it is easier to reach saturation with active users in each group as the number of users to be served in the system increases.
The sum rate performance of the three proposed algorithms is compared in Fig. 6. It can be seen that the sum rates of all three algorithms are significantly higher in the configuration for N t = 64 than for N t = 32; the sum rates of Algorithm 1 are higher than those of Algorithms 2 and 3 in each configuration, while the sum rates of Algorithm 2 are slightly higher than those of Algorithm 3, which are closer to each other. When the number of antennas at the BS side is N t = 32, the increasing trend of the sum rate in the system smoothes out significantly as the number of users in the system increases, which is due to the fact that each group is basically served by users as the number of users increases. For N t = 64, there is only a moderating trend in the sum rate as the number of users increases. This is due to the fact that more antennas at the transmitter end divide more intervals, and the number of users in the graph is not yet such that there is a greater probability that each interval will have a serving user.
The performance of the proposed three algorithms in terms of the number of simultaneous users served is compared in Fig. 7. For the total number of users to be served simultaneously in the system, the performance of all three algorithms in the N t = 64 configuration is higher than when N t = 32; as the number of users to be served in the system increases, the curves of the three algorithms in the N t = 32 configuration flatten out significantly; thanks to the larger number of groups and that each group can serve more users at the same time through the PD-NOMA, at N t = 64, the moderating trend of the boost in the curves of the number of users served simultaneously by Algorithms 2 and 3 is not as pronounced as in Algorithm 1. In addition, there is an intersection point in the curves of the number of NOMA served simultaneously in the system for both antenna configurations. At the beginning, when the number of users to be served in the system is low, the number of NOMA users served is smaller for N t = 64 than for N t = 32. This is due to the fact that with more antennas at the transmitter there are more grouping intervals and the users served are more likely to be in separate groups and therefore have less opportunity for NOMA pairing. As the number of NOMA users served rises, starting at about 15 users in the figure, the number of NOMA users in the two algorithms for N t = 32 is smaller than for N t = 64, due to the larger range of angles per grouping interval at N t = 32, and as the number of active grouping pairs in the system becomes larger, additional users to be served with the strongest path of elevation angles in the AoD have a significantly higher probability of appearing in the paired groupings. In contrast, with N t = 64, there is still a higher probability that a newly added user to be served will appear in an unpaired grouping when there are the same number of grouping pairs in the system because of the greater number of group intervals divided. Comparing Figs. 6 and 7, compared to Algorithm 1, Algorithms 2 and 3 increase the number of simultaneous served users but decrease the sum rate in the system, due to the fact that when NOMA pairing is successful, part of the transmit power is allocated to NOMA user 1, which has a poor channel condition, reducing the rate of NOMA user 2, which has a better channel condition. Compared to Algorithm 2, Algorithm 3 slightly increases the number of NOMA users served simultaneously in the system; however, it also further reduces the sum rate of the system due to the increased power allocated to NOMA user 1 in the side lobe NOMA pair. Since the number of simultaneous NOMA users served simultaneously in Algorithms 2 and 3 is similar for the same antenna configuration, it can be seen that the number of successful side lobe NOMA pairs in Algorithm 3 is extremely small, and as Algorithm 3 is more complex than Algorithm 2, it is generally better to employ only Algorithm 2 to increase the number of simultaneous users served in the system.

Conclusion
In this article, we studied the downlink cellular communication scheme for serving multiuser through a wideband phased linear array at the BS side with frequency squint beams. First, a method using steering main beams of a virtual planar array to simulate the donutshaped main beam of the linear array was proposed to solve the matching problem of the steering main beam and the channel model. Then, an algorithm drawing on the concept of JSDM grouping to divide the transmitter elevation angle interval into multiple small intervals and serve multi-user by means of squint beams from the linear array was proposed, where at most one user can be served by a single squint beam in each small interval. Subsequently, the algorithm supporting PD-NOMA for two users whose transmit channel path directions are both within the half-power beam range of a single quint donut-shaped main beam and the algorithm which supports PD-NOMA pairing of the main direction of a single squint beam and its first side lobe direction were raised for serving more users, for which the feasible domain of a two-user PD-NOMA channel for a given SIC decoding power constraint in the single-antenna scenario was briefly discussed and placed in Appendix 2 for consistency and conciseness in the bulk of the paper. The feasibility of the proposed scheme and the performance of the proposed algorithms were finally verified by simulation.
where {1, 2, · · · , N} is the set of users in the NOMA pair, R i is the rate of the i th user, and β i is the proportion of transmit power allocated to that user's signal at the BS side. Firstly, it is necessary to select the appropriate NOMA pair such that the channel quality differs significantly across users. Secondly, the power allocation strategy of the NOMA pair is to maximize the rate of the users with good channel quality while guaranteeing the minimum QoS requirements of the users with poor channel quality. If it is not possible to simultaneously serve the users in the pair to guarantee their QoS, then the lower priority users in the NOMA group are dropped or the users with good channel conditions are served. In a two-user NOMA pairing scheme in which the BS side and the subscriber end are equipped with a single-antenna scenario, it is assumed without loss of generality that the channel quality of user 2 is better than that of user 1, i.e., |h 1 | 2 < |h 2 | 2 , and this optimization problem is equivalent to satisfying the minimum rate requirement t 1 for the far user with poor channel quality to maximize the rate of the near user with a better channel condition, that is max β i i∈{1,2} R 2 , s.t. β 1 + β 2 ≤ 1, This problem has been solved, and according to the work in [17], it follows that Substituting β 2 = 1 − β 1 into (62), we get The power allocation coefficient for user 2 in the pair needs to satisfy where the left inequality is due to the fact that the effective power distribution needs to be positive and the right inequality is a restriction on the proportion of power that can be successfully decoded by the SIC; substituting these two constraints into (64) respectively gives 1 − β max 2 2 t 1 |h 1 | 2 ≤ 2 t 1 − 1 σ 2 n /P.
• When 1 − β max 2 2 t 1 ≤ 0, since the left side of the inequality in (67) is non-positive, (67) holds constant, it follows that It means that the constraint imposed by the right-hand inequality in (65) vanishes as long as (68) is satisfied. • When 1 − β max 2 2 t 1 > 0, it follows from (67) that If this NOMA pair matches, then there is As can be seen when the NOMA pair is matched, the worse the channel quality of user 1, the smaller the rate of user 2; the better the channel quality of user 2, the larger the rate of user 2. Let the minimum rate requirement for user 2 be t 2 , i. e., R 2 ≥ t 2 needs to be satisfied; substituting (71), we have (72) Figure 8 shows the feasible domain of the channel state for the two-user PD-NOMA pair in (72), illustrated by t 1 = t 2 = 0.2 and P/σ 2 n = 1, where the curved part is one of the hyperbolas determined in (72), the vertical dashed line on the left is determined by the left inequality in (65), and the vertical dashed line on the right is set according to the SIC decoding requirements and determined by the right inequality in (65). When (68) holds, the entire shaded region is the NOMA feasible domain, whereas when (68) is not satisfied, the shaded region on the left is the feasible domain of the actual matched pair.