I Introduction
Devicetodevice (D2D) communication underlaying cellular networks, which allows a device to communicate with its proximity device over the licensed cellular bandwidth, is recognized as a promising technology in future networks due to its advantages such as high spectrum efficiency, high energy efficiency (EE) and low transmission delay [asadi2014survey]. Interference management is the most important challenge for underlaying D2D communication[tehrani2014device]. The D2D link and the cellular link operating in the same licensed band interfere with each other severely[DSABook2019Liang], and the interference needs to be carefully suppressed via efficient interference management[feng2013device] and resource allocation[OlavD2DTWCHC2011].
Recently, reconfigurable intelligent surfaces (RIS) have emerged as a new and revolutionary technology to achieve spectrum, energy and costefficient wireless networks[LiangLISA2019][huang2019holographic]. Specifically, RIS consist of a large number of passive lowcost reflecting elements, each of which can adjust the phase and amplitude of the incident electromagnetic wave and reflect it passively[8936989][YangLiangNOMAIRS19]. Thus, RIS are able to enhance desired signals and suppress interference by designing the reflecting coefficient (including phase and amplitude) of each reflecting element. For instance, the weightedsum rate of an RISaided multiuser multipleinput singleoutput downlink communication system was maximized in [8982186], by jointly optimizing the base station’s (BS’s) active beamforming and the RIS’s passive beamforming (i.e., reflecting coefficients). The EE of an RISempowered downlink multiuser communication system was maximized in [8741198], by jointly optimizing the BS’s transmit power and the RIS’s passive beamforming.
RIS can be explored to enhance the strengths of desired signals for both D2D links and cellular links, and suppress the severe interference between each paired D2D link and cellular link. This motivates us to study RISempowered D2D communication underlaying a cellular network in this paper. The main contributions are summarized as follows.

We formulate a problem to maximize the overall network sum rate, by jointly optimizing the resource reuse indicators (i.e., user pairing between D2D users and cellular users (CUs)), the transmit power and the RIS’s passive beamforming, subject to the signaltointerferenceplusnoise ratio (SINR) constraints for both D2D links and cellular links. However, the problem is challenging to be solved optimally, since the user pairing (involving integer variables) and the resource allocation are closely coupled.

To decouple the problem, we first propose an efficient relativechannelstrength based userpairing scheme with low complexity. Under the obtained userpairing design, an iterative algorithm based on alternating optimization is further proposed. The successive convex approximation technique is exploited to optimize the transmit power; while the Lagrangian dual transform and quadratic transform techniques are utilized to optimize the passive beamforming. The algorithm’s convergency is proved and its complexity is analyzed.

Numerical results show that the proposed design achieves significant sumrate enhancement compared to traditional underlaying D2D without RIS, and suffers from slight degradation compared to the bestachievable performance under ideal user pairing.
Ii System Model
As shown in Fig. 1, we consider an RISempowered cellular network with underlay D2D, which consists of an RIS, D2D transmitters (TXs) denoted as , D2D receivers (RXs) denoted as , active CUs (i.e., cellular users) denoted as , and a cellular BS. The RIS has reflecting elements, while each D2D TX, D2D RX, CU and the BS are equipped with a single antenna. A controller is attached to the RIS to control the reflecting coefficients and communicate with other network components through separate wireless links. We assume that the D2D links share the uplink (UL) spectrum of the cellular network, since the UL spectrum is typically underutilized compared to the downlink spectrum. To alleviate interference, we consider that a D2D link shares at most one CU’s spectrum resource, while the resource of a CU can be shared by at most one D2D link [feng2013device] [ramezani2017joint] [8299474].
All channels are assumed to experience quasistatic flatfading. The channels from to and RIS are denoted by and , respectively. For notational clarity, we represent each channel related to the cellular network with a tilde. The channels from to BS and RIS are denoted by and , respectively; the channels from RIS to RX l and BS are denoted by and , respectively; the interference channels from TX i to BS and from CU k to RX l are denoted by and , respectively.
The transmitted signals from TX i and CU k are denoted as and
, respectively, which follow independent circularly symmetric complex Gaussian (CSCG) distribution with zero mean and unit variance, i.e.,
, . Denote the index set of active D2D pairs as . The corresponding SINR for RX n decoding from is(1) 
where and are the transmit power of TX i and CU k, respectively; denotes the reflecting coefficient matrix, where and ; is the resource reuse indicator for cellular link and D2D link , when D2D link reuses the resource of CU , and otherwise; is the power of additive white Gaussian noise (AWGN) at RX n.
The SINR for BS decoding from CU k is
(2) 
where is the power of AWGN at BS.
Hence, the overall network’s sum rate in bps/Hz is
(3) 
where the length(
) resource reuse indicator vector
, and the length() power allocation vector .Iii problem formulation
This paper aims to maximize the sum rate in (3), by jointly optimizing the resource reuse indicator vector , the transmit power vector and the reflecting coefficients matrix . The optimization problem is formulated as
(4a)  
s.t.  (4b)  
(4c)  
(4d)  
(4e)  
(4f)  
(4g)  
(4h)  
(4i) 
where (4b) and (4c) indicate the required minimum SINRs (i.e., qualityofservice) and for the D2D links and cellular links, respectively; (4d) ensures that a D2D link shares at most one CU’s resource, while (4e) indicates that the resource of a CU can be shared by at most one D2D link; (4f) and (4g) are the maximum transmit power constraints on the TXs and CUs, respectively; (4h) and (4i) are the practical constraints on the reflecting coefficients.
Notice that (P1) is a nonconvex problem. First, (P1) involves integer variables and thus is NPhard. Moreover, the objective function and the constraint functions of (4b) and (4c) are nonconcave with respect to the variables , and , and these variables are all coupled. There is no standard method to solve such a nonconvex problem.
Iv Solution to (P1)
In order to solve (P1) effectively, we first propose an efficient userpairing scheme to determine integer variables , then optimize and in an iterative manner.
Iva RelativeChannelStrength based Pairing Scheme
Since the userpairing design involves integer programming which is hard to solve, we propose a relativechannelstrength (RCS) based lowcomplexity pairing scheme to design the resource reuse indictors .
Notice that there are different possible pairings denoted as a set . Each possible pairing can be viewed as an index mapping denoted as , for , i.e., the mapping maps each CU index to a D2Dlink index . The RCSbased pairing scheme determines the pairing by the following criterion
(5) 
This heuristic pairing scheme chooses the pairing mapping which maximizes the sum of the ratios of each pairedCUtoBS channel strength over the paired CUtoRX interference channel strength. Clearly, this heuristic pairing scheme that requires only simple comparison features low complexity, but fortunately its resultant design only suffers from slight performance degradation compared to the design with ideal pairing achieved by exhaustive search, as numerically shown in Section
V.IvB Algorithm for Solving (P1) with Given Pairing Design
After heuristic pairing, we apply the alternating optimization (AO) [hong2015unified] algorithm to decouple the variables and . For given , we optimize based on the successive convex approximation (SCA) technique [beck2010sequential]. For given , we optimize based on the Lagrangian dual transform and quadratic transform techniques [shen2018fractional].
IvB1 Optimizing Transmit Power Vector
In each iteration , for given reflecting coefficient matrix , the transmit power vector can be optimized by solving
(6a)  
s.t.  (6b) 
The objective function of (P1.1) is nonconvex due to its dependence. We exploit the SCA technique to solve (P1.1). Specifically, we need to find a concave lower bound to approximate the objective function. Since any convex function can be lower bounded by its firstorder Taylor expansion at any point, we obtain the following concave lower bound at the point
(7) 
where , , , , and .
With given and , (P1.1) is approximated as
(8a)  
s.t.  (8b) 
Problem (P1.2) is a convex problem which can be efficiently solved with standard methods, e.g., CVX[grant2008cvx].
IvB2 Optimizing Reflecting Coefficient Matrix
In each iteration , for given transmit power vector , the reflecting coefficient matrix can be optimized by solving
(9a)  
s.t.  (9b) 
We tackle the logarithm in the objective function via the Lagrangian dual transform technique. Introducing auxiliary variables and , the new objective function can be equivalently expressed as
(10) 
It is easy to validate that the optimal values of and are and , respectively. We define , , and . From (1) and (2), optimizing the reflecting coefficient can be equivalently transformed into optimizing in the following objective function
(11) 
where , , and .
Then, utilizing the quadratic transform method proposed in [shen2018fractional], we introduce an auxiliary variable , transforming (11) to
(12) 
We first optimize with fixed , then optimize with fixed . It can be easily confirmed that is a concave differentiable function over with fixed , so the optimal solution of can be obtained by setting . The optimal value of is given by
(13)  
(14) 
Then, we replace and with and , respectively. Denote , , , , , , and . Optimizing for given , the objective function is transformed as follows
(15) 
where is a constant, the matrix and vector are
(16)  
(17) 
with , , , and .
Similarly, leveraging the quadratic transform method, the constraints (4b) and (4c) can be transformed into
(18) 
(19) 
where , , , , and .
Hence, (P2.1) is transformed into the following problem
(20)  
s.t.  (21) 
The resulting (P2.2) is a quadratic constrained quadratic programming (QCQP) problem. Thus, (P2.2) can also be effectively solved by standard methods.
IvC Overall Algorithm
Theorem 1.
Algorithm 1 is guaranteed to converge.
Proof.
First, in Step 3, since the optimal solution is obtained for given , we have the following inequality on the sum rate
(22) 
where (a) and (c) hold since the Taylor expansion in (7) is tight at given local points and , respectively, and (b) holds since is the optimal solution to (P1.2).
Second, in Step 4, since is the optimal solution to (P2.2), we can obtain the following inequality
(23) 
Since the objective value is nondecreasing after each iteration and is upper bounded by some positive constant, the overall Algorithm 1 is guaranteed to converge. ∎
Problems (P1.2) and (P2.2) are alteratively solved in each outerlayer AO iteration. Specifically, (P1.2) can be solved in operations by the extended waterfilling algorithm[liu2013complexity], while (P2.2) is a convex QCQP which can be solved by using interior point methods with complexity [hassanien2008robust]. Hence, the complexity of Algorithm 1 is , where denotes the number of outerlayer AO iterations.
V Numerical Results
This section provides numerical results to validate the performance of the proposed RISempowered underlaying D2D network. We assume that , , and are independently Rayleigh fading distributed, while , , and follow independent Rician fading distribution, i.e.,
(25) 
where is the Rician factor of , is the line of sight (LoS) component, and is the nonLoS (NLOS) component each element of which follows distribution . Set . Similarly, , and are generated in the same way as with Rician factors .
We assume that the CUs are uniformly distributed in the cellular cell with radius
meters (m). We adopt the clustered distribution model in [feng2013device] for D2D users, i.e., the clusters are randomly located in the cell, and each D2D link is uniformly distributed in one cluster with radius m. We set and . The RIS is located between two D2D clusters. The locations of CUs and D2D users are generated by the above method. The largescale path losses from TXs and CUs to RXs are , from TXs and CU to BS are , the other passlosses of RISrelated channels are [feng2013device][8982186], where is the distance in meters. We set , dBm , bps/Hz, bps/Hz, and dBm[feng2013device]. The simulation results are based on 1000 channel realizations.For comparison, we consider the following two benchmarks, i.e., (1) underlaying D2D without RIS, (2) RISempowered D2D with ideal user pairing. For the first benchmark, we consider the traditional RISempowered D2D underlaying a cellular network without RIS. We maximize the sum rate by jointly optimizing the resource reuse indicator and the transmit power vector . We use the solving algorithm in [feng2013device] for this problem, and omit the details herein. For the second benchmark, we exhaustively search over possible userpairings, and jointly optimize as well as under each pairing. This benchmark gives an achievable upperbound sumrate performance of the RISempowered underlaying D2D network.
Fig. 2 plots the sum rate versus the maximum transmit power of the transmitters. The proposed design achieves significant sum rate gain compared to the first benchmark. For instance, the sum rate of the proposed RCSbased design is and higher than that of the first benchmark when is 10 and 20 dBm, respectively. In addition, compared to the second benchmark based on exhaustive search, the proposed design with low complexity has slight degradation of performance. Furthermore, the sum rate of the proposed design increases as increases, since more reflecting elements can further enhance equivalent channel strengthes and suppress the interlink interferences. Also, the finiteresolution phase shifters of reflecting elements usually degrade the sumrate performance. The sum rate increases with the increase of the phaseshift quantization bits . In particular, the 2bit phase shifter can obtain sufficiently high performance gain with a slight performance degradation compared to the ideal case of continuous phase shifters.
Fig. 3 plots the sum rate versus the CUs’ minimum rate requirement . The sum rate decreases as increases, which reveals the rate tradeoff between the D2D links and the cellular links. Compared to the first benchmark, the proposed design achieves significant sumrate gain by introducing the RIS. Compared to the second benchmark, the proposed design suffers from slight sumrate performance degradation, but it obviously outperforms this benchmark in terms of computational complexity. The proposed design solves the jointresourceallocation optimization problem only once, while this benchmark needs to solve such problem for times under all possible pairings, resulting into unaffordable complexity especially for large numbers of D2D and cellular links.
Vi Conclusion
This paper has studied an RISempowered underlaying D2D communication network. The overall network sum rate is maximized by jointly optimizing the resource reuse indicators, the transmit power and the passive beamforming. First, an efficient relativechannelstrength based userpairing scheme with low complexity is proposed to determine the resource reuse indicators. Then, the transmit power and the passive beamforming are optimized by utilizing the proposed alternatingoptimization based iterative algorithm. Numerical results show that the proposed design achieves significant performance enhancement compared to traditional underlaying D2D network without RIS, and suffers from slight performance degradation compared to RISempowered underlying D2D with ideal userpairing. This work can be extended to other scenarios such as multiantenna BS/users and multiple RISs.
Comments
There are no comments yet.