Energy Efficiency Optimization in Massive MIMO Secure Multicast Transmission

Abstract

Herein, we focus on energy efficiency optimization for massive multiple-input multiple-output (MIMO) downlink secure multicast transmission exploiting statistical channel state information (CSI). Privacy engineering in the field of communication is a hot issue under study. The common signal transmitted by the base station is multicast transmitted to multiple legitimate user terminals in our system, but an eavesdropper might eavesdrop this signal. To achieve the energy efficiency utility–privacy trade-off of multicast transmission, we set up the problem of maximizing the energy efficiency which is defined as the ratio of the secure transmit rate to the power consumption. To simplify the formulated nonconvex problem, we use a lower bound of the secure multicast rate as the molecule of the design objective. We then obtain the eigenvector of the optimal transmit covariance matrix into a closed-form, simplifying the matrix-valued multicast transmission strategy problem into a power allocation problem in the beam domain. By utilizing the Minorize-Maximize method, an iterative algorithm is proposed to decompose the secure energy efficiency optimization problem into a sequence of iterative fractional programming subproblems. By using Dinkelbach’s transform, each subproblem becomes an iterative problem with the concave objective function, and it can be solved by classical convex optimization. We guarantee the convergence of the two-level iterative algorithm that we propose. Besides, we reduce the computational complexity of the algorithm by substituting the design objective with its deterministic equivalent. The numerical results show that the approach we propose performs well compared with the conventional methods.

1. Introduction

With the great development of network information technology and the continuous improvement of wireless communication technology, modern society puts forward higher requirements for the current communication transmission. With the extensive amount of research, there are a number of breakthrough technologies, such as the multiple-input multiple-output (MIMO) technology [1]. This technology has great potential to improve both the spectrum efficiency and the energy efficiency of the system, which is the development direction of future 5G/B5G wireless communication systems [2,3]. In the long-term evolution (LTE) standard, the wireless multicast transmission is also widely studied. In wireless multicast transmission, one base station (BS) sends public information to multiple user terminals (UTs). That is to say, that the BS can effectively transmit group-oriented signals by utilizing an evolved multimedia broadcast multicast service (eMBMS) [4,5,6,7,8]. As a mature technology, wireless multicast technology has certain advantages in improving the network traffic and reducing the communication cost. The combination of wireless multicast and massive MIMO is widely concerned with improving the quality of service (QoS), owing to its performance of effectively shaping the multicast transmission signals [9,10,11,12].

Due to the openness of the wireless propagation media, security issues are critical indicators in wireless transmissions of all fields. To ensure the security of transmission, many works traditionally use key-based enciphering. However, this method meets its challenge for secure wireless communications with additional control signaling and feedback channel management. Therefore, the security and privacy of physical layer security in wireless channels have a wide range of research [13,14,15], thanks to its great ability to work in a keyless scenarios. In previous works, the physical layer security of massive MIMO unicast transmission has been studied [16,17,18,19]. Meanwhile, in some existing works, the physical layer security of multicast transmission is also considered. For example, secure multicast precoding has been investigated [20,21,22]. Besides, some works have studied the security by introducing artificial noise into the multicast transmission, e.g., in [23,24].

Energy efficiency is another important indicator of performance in designing wireless transmission, which has been widely concerned by both academia and industry [25,26]. About 5% of the energy consumption of the whole world is produced by information and communication technologies, owing to its connected devices and infrastructure constantly being consumed [27]. In the future of 5G/B5G development, energy demand will soon become unmanageable. Therefore, a promising answer to this issue is to find considerable metrics which can limit energy consumptions but meet the needs of secure and reliable communication. However, most of the works on secure multicast transmission are limited to optimizing the secure rate, and almost no energy efficiency issues are considered. Some existing works have investigated the energy efficiency of unicast transmission instead of multicast transmission [28,29,30]. In some previous works, the energy efficiency multicast transmission was investigated [31,32], but the issues of security were not solved well.

In some works, it is assumed that the legitimate UTs’ instantaneous channel state information (CSI) is known [17,21,22,23,33]. However, due to CSI’s capture/feedback overhead in frequency division duplex MIMO systems, there are some difficulties in obtaining instantaneous CSI at the BS. It is the reason why most of the previous work has assumed this. On the other hand, statistical CSI has enormous advantages over instantaneous CSI for BS to obtain. Because of its slower change, BS can get statistical CSI accurately [34]. Therefore, statistical CSI is more suitable for wireless transmission design in high mobility scenarios. Besides, as a vast antenna array with high spatial resolution, massive MIMO channel typically exhibits new statistical characteristics [35,36].

We study the energy efficiency of massive MIMO secure multicast transmission with the assumption that only the statistical CSI of legal UTs and the eavesdropper can be obtained at the BS. We summarize the contributions of the current paper as follows.

• We first considered a secure multicast transmission scenario, in which the BS can only obtain the statistical CSI of the legitimate UTs and eavesdropper, and formulate our energy efficiency maximization problem. A tight lower bound rate of secure multicast is presented to simplify the optimization objective. Therefore, our problem is transformed into finding the optimal solution of this equivalent problem.
• Via making the transmit signaling directions into a closed-form, we simplify the energy efficiency optimization of matrix-valued strategy design to a beam domain power allocation problem, reducing the optimization variables significantly.
• By invoking the Minorize–Maximize (MM) method and Dinkelbach’s transform, we propose a two-level iterative algorithm to solve the equivalent problem. In the outer iteration, we decompose the original problem into several fractional programming subproblems, where the numerator is concave and the denominator is linear. In the inner iteration, we introduce related variables to remove the denominator and transform each subproblem into a common convex optimization problem to obtain the optimal solution of the equivalent problem. Our simulation results reveal that this solution is almost equal to the original problem.
• According to the large-dimensional random matrix theory, we use the deterministic equivalent of the secure transmit rate instead of massive sample averaging to reduce the computational complexity, which is proved to be almost the same as the original result.

The rest of this paper is arranged as follows. Section 2 presents the massive MIMO channel model. In Section 3, we divide the algorithm into three parts. The energy efficiency of secure multicast transmission is investigated. The optimal energy efficiency multicast transmit covariance matrix is decomposed into eigenvectors and eigenvalues as pretreatment. We propose an iterative algorithm for some fractional programming subproblems. Each subproblem is also solved by the inner iterative algorithm of sequential beam domain power allocation. The last part is to simplify the complexity of the problem. In Section 4, numerical simulation results are given. Finally, we summarize the full work in Section 5.

In our work, the use of symbols is also an important part. The main adopted symbols are as follows.

• ${\mathbb{C}}^{M\times N}$ denotes $M\times N$ dimensional complex-valued vector space and ${\mathbb{R}}^{M\times N}$ denotes $M\times N$ dimensional real-valued vector space.
• ${\left[\mathbf{A}\right]}_{m,n}$ denotes the $(m,n)$ th element of matrix $\mathbf{A}$.
• det $\{·\}$ denotes the matrix determinant operation, tr $\{·\}$ denotes the matrix trace operation, and diag $\{·\}$ denotes the diagonalization of the matrix, respectively.
• $\mathsf{E}\left\{·\right\}$ denotes the expectation operation and ${\left[x\right]}^{+}$ denotes max$\left\{x,0\right\}$.
• ${(·)}^T$, ${(·)}^{*}$, ${(·)}^H$ and ${(·)}^{-1}$ denote the matrix transpose, conjugate, conjugate-transpose and inverse, respectively.
• ⊙ denotes the Hadamard product.
• $\mathbf{X}⪰\mathbf{0}$ denotes that $\mathbf{X}$ is a positive semidefinite matrix.

2. Channel Model

We consider a massive MIMO secure multicast downlink transmission system, in which public messages sent by the BS are transmitted to the legitimate UTs. There is one BS with M antenna transmits and K UTs, each with $N_r$ antennas. The public messages are secret to the eavesdropper with $N_e$ antennas.

$\mathbf{x}\in {\mathbb{C}}^{M\times 1}$ denotes the signal that is sent from the BS in multicast transmission. $\mathsf{E}\left\{\mathbf{x}\right\}=\mathbf{0}$ and $\mathsf{E}\left\{\mathbf{x}{\mathbf{x}}^H\right\}=\mathbf{Q}\in {\mathbb{C}}^{M\times M}$, where $\mathbf{Q}$ denotes the transmit covariance matrix. We model the signals received by legitimate UT k and the eavesdropper as

$$\begin{array}{cc}\hfill {\mathbf{y}}_k& ={\mathbf{H}}_k\mathbf{x}+{\mathbf{n}}_k\in {\mathbb{C}}^{N_{\mathrm{r}}\times 1},\hfill \end{array}$$

(1)

$$\begin{array}{cc}\hfill {\mathbf{y}}_{\mathrm{ev}}& ={\mathbf{H}}_{\mathrm{ev}}\mathbf{x}+{\mathbf{n}}_{\mathrm{ev}}\in {\mathbb{C}}^{N_{\mathrm{e}}\times 1},\hfill \end{array}$$

(2)

respectively, where ${\mathbf{H}}_k$ and ${\mathbf{H}}_{\mathrm{ev}}$ are defined as the channel matrices of BS-UT k and BS-eavesdropper, respectively.

We consider our MIMO channel model to follow jointly spatially correlated Rayleigh fading, which means there are joint correlation characteristics between the transmitter and the receivers [37]. Thus, the channel matrices ${\mathbf{H}}_k$ and ${\mathbf{H}}_{\mathrm{ev}}$ in (1) and (2) can be expressed as

$$\begin{array}{cc}\hfill {\mathbf{H}}_k& ={\mathbf{U}}_k{\mathbf{G}}_k{\mathbf{V}}_k^H\in {\mathbb{C}}^{N_{\mathrm{r}}\times M},\hfill \end{array}$$

(3)

$$\begin{array}{cc}\hfill {\mathbf{H}}_{\mathrm{ev}}& ={\mathbf{U}}_{\mathrm{ev}}{\mathbf{G}}_{\mathrm{ev}}{\mathbf{V}}_{\mathrm{ev}}^H\in {\mathbb{C}}^{N_{\mathrm{e}}\times M},\hfill \end{array}$$

(4)

respectively, where ${\mathbf{G}}_k\in {\mathbb{C}}^{N_{\mathrm{r}}\times M}$ and ${\mathbf{G}}_{\mathrm{ev}}\in {\mathbb{C}}^{N_{\mathrm{e}}\times M}$ are random matrices with zero mean and independent elements, which are defined as the channel matrices of BS-UT k and BS-eavesdropper in the beam domain, respectively [34,38]. ${\mathbf{U}}_k\in {\mathbb{C}}^{N_{\mathrm{r}}\times N_{\mathrm{r}}}$, ${\mathbf{V}}_k\in {\mathbb{C}}^{M\times M}$, ${\mathbf{U}}_{\mathrm{ev}}\in {\mathbb{C}}^{N_{\mathrm{e}}\times N_{\mathrm{e}}}$, and ${\mathbf{V}}_{\mathrm{ev}}\in {\mathbb{C}}^{M\times M}$ are unitary matrices of transmitter and receivers, respectively. The statistics of the beam domain channels can be described as

$$\begin{array}{cc}\hfill {\mathsf{\Omega }}_k& =\mathsf{E}\left\{{\mathbf{G}}_k\odot {\mathbf{G}}_k^{*}\right\}\in {\mathbb{R}}^{N_{\mathrm{r}}\times M},\hfill \end{array}$$

(5)

$$\begin{array}{cc}\hfill {\mathsf{\Omega }}_{\mathrm{ev}}& =\mathsf{E}\left\{{\mathbf{G}}_{\mathrm{ev}}\odot {\mathbf{G}}_{\mathrm{ev}}^{*}\right\}\in {\mathbb{R}}^{N_{\mathrm{e}}\times M},\hfill \end{array}$$

(6)

respectively. It can be proved that as the number of antennas in the BS, M, increases to infinity, the eigenvectors of the correlation matrices in (5) and (6) tend to be equal, which can be denoted by a deterministic unitary matrix $\mathbf{V}$ [38]. $\mathbf{V}$ is only effected by the topology of transmit array at the BS, which means the channel matrices can be rewritten as

$$\begin{array}{cc}\hfill {\mathbf{H}}_k& \stackrel{M\to \infty }{=}{\mathbf{U}}_k{\mathbf{G}}_k{\mathbf{V}}^H,\hfill \end{array}$$

(7)

$$\begin{array}{cc}\hfill {\mathbf{H}}_{\mathrm{ev}}& \stackrel{M\to \infty }{=}{\mathbf{U}}_{\mathrm{ev}}{\mathbf{G}}_{\mathrm{ev}}{\mathbf{V}}^H,\hfill \end{array}$$

(8)

respectively. It has been widely used in previous works and shows considerable accuracy for a finite number of BS antennas [35,36]. Thus, we adopt the channel models in (7) and (8) in this paper.

3. Energy Efficiency Optimization of Secure Multicast Transmission

We study the energy efficiency multicast optimization of secure massive MIMO in this section. We assume that the BS can only capture the statistical CSI of all legitimate UTs and the eavesdropper [34]. Besides, legitimate UTs and the eavesdropper can obtain the instantaneous CSI. Before constructing the energy efficiency model, we first analyze the secure transmit rate. The received signal transmission rate needs to be higher for any user terminal than that received by the eavesdropper. Focusing on delay-tolerant scenarios [39], the ergodic secure multicast rate can be expressed as

$$\begin{array}{c}\hfill R_{\mathrm{se}}\triangleq {\left[R_{\mathrm{mc}}-R_{\mathrm{ev}}\right]}^{+},\end{array}$$

(9)

where the ergodic multicast rate $R_{\mathrm{mc}}$ is expressed as

$$\begin{array}{cc}\hfill R_{\mathrm{mc}}& =\underset{k}{min}\mathsf{E}\left\{logdet\left\{{\mathbf{I}}_{N_{\mathrm{r}}}+{\mathbf{H}}_k\mathbf{Q}{\mathbf{H}}_k^H\right\}\right\}\hfill \\ \hfill & =\underset{k}{min}\mathsf{E}\left\{logdet\left\{{\mathbf{I}}_{N_{\mathrm{r}}}+{\mathbf{G}}_k{\mathbf{V}}^H\mathbf{Q}\mathbf{V}{\mathbf{G}}_k^H\right\}\right\},\hfill \end{array}$$

(10)

and the ergodic rate for eavesdropper $R_{\mathrm{ev}}$ is expressed as

$$\begin{array}{cc}\hfill R_{\mathrm{ev}}& =\mathsf{E}\left\{logdet\left\{{\mathbf{I}}_{N_{\mathrm{e}}}+{\mathbf{H}}_{\mathrm{ev}}\mathbf{Q}{\mathbf{H}}_{\mathrm{ev}}^H\right\}\right\}\hfill \\ \hfill & =\mathsf{E}\left\{logdet\left\{{\mathbf{I}}_{N_{\mathrm{e}}}+{\mathbf{G}}_{\mathrm{ev}}{\mathbf{V}}^H\mathbf{Q}\mathbf{V}{\mathbf{G}}_{\mathrm{ev}}^H\right\}\right\}.\hfill \end{array}$$

(11)

In (10) and (11), we exploit the identity $det\left\{\mathbf{I}+\mathbf{A}\mathbf{B}\right\}=det\left\{\mathbf{I}+\mathbf{B}\mathbf{A}\right\}$.

Meanwhile, we adopt a model of power consumption as follows. For large-scale MIMO systems, most of the power is consumed by the BS, we then focus mainly on the power consumption of the BS. Besides, the static and dynamic power consumption of the radio frequency chains is considered. Specifically, our power consumption model is adopted as

$$\begin{array}{cc}\hfill P\left(\mathbf{Q}\right)& \triangleq \zeta \mathsf{tr}\left\{\mathbf{Q}\right\}+M{P}_{\mathrm{c}}+P_{\mathrm{s}},\hfill \end{array}$$

(12)

where $\mathsf{tr}\left\{\mathbf{Q}\right\}$ denotes the transmit power of multicast transmission. Note that the coefficient $\zeta \ge 1$ is the reciprocal of the transmit amplifier drain efficiency, $P_{\mathrm{s}}$ and $P_{\mathrm{c}}$ are the static power consumption and the hardware dissipated power of the BS, respectively.

We denote W by the bandwidth of secure transmission, the energy efficiency of massive MIMO secure multicast transmission is formulated by the ratio of secure transmit rate to power consumption as follows,

$$\begin{array}{c}\hfill {\mathrm{EE}}_{\mathrm{se}}=\frac{W{R}_{\mathrm{se}}}{P}.\end{array}$$

(13)

Our aim is optimizing the transmit covariance matrix $\mathbf{Q}$, leading to maximizing the ${\mathrm{EE}}_{\mathrm{se}}$ of downlink secure multicast in (13). The problem is addressed as

$$\begin{array}{cc}\hfill \underset{\mathbf{Q}}{argmax}& \frac{W{R}_{\mathrm{se}}}{P},\hfill \\ \hfill \mathrm{s}.\mathrm{t}.& \mathsf{tr}\left\{\mathbf{Q}\right\}\le P_{\max },\mathbf{Q}⪰\mathbf{0},\hfill \end{array}$$

(14)

where $P_{\max }$ denotes the maximum multicast power budget at the transmitter.

$R_{\mathrm{ev}}$ in (11) can get its upper bound via using Jensen’s inequality, which is expressed as follows,

$$\begin{array}{cc}\hfill R_{\mathrm{ev}}& \le R_{\mathrm{ev},\mathrm{ub}}\triangleq logdet\left\{{\mathbf{I}}_{N_{\mathrm{e}}}+\underset{\triangleq {\mathbf{A}}_{\mathrm{ev}}\left({\mathbf{V}}^H\mathbf{Q}\mathbf{V}\right)}{\underbrace{\mathsf{E}\left\{{\mathbf{G}}_{\mathrm{ev}}{\mathbf{V}}^H\mathbf{Q}\mathbf{V}{\mathbf{G}}_{\mathrm{ev}}^H\right\}}}\right\},\hfill \end{array}$$

(15)

where ${\mathbf{A}}_{\mathrm{ev}}\left(\mathbf{X}\right)\triangleq \mathsf{E}\left\{{\mathbf{G}}_{\mathrm{ev}}\mathbf{X}{\mathbf{G}}_{\mathrm{ev}}^H\right\}\in {\mathbb{C}}^{N_{\mathrm{e}}\times N_{\mathrm{e}}}$ denotes a function of matrix-valued form, which can output a diagonal matrix. Its ith diagonal element can be expressed as

$$\begin{array}{c}\hfill {\left[{\mathbf{A}}_{\mathrm{ev}}\left(\mathbf{X}\right)\right]}_{i,i}=\mathsf{tr}\left\{\mathsf{diag}\left\{{\left({\left[{\mathsf{\Omega }}_{\mathrm{ev}}\right]}_{i,:}\right)}^T\right\}\mathbf{X}\right\}.\end{array}$$

(16)

By employing the function in (15), the secure multicast rate in (9) can obtain its lower bound as

$$\begin{array}{c}\hfill R_{\mathrm{se}}\ge R_{\mathrm{se},\mathrm{lb}}\triangleq {\left[R_{\mathrm{mc}}-R_{\mathrm{ev},\mathrm{ub}}\right]}^{+}.\end{array}$$

(17)

For a vast region of signal-to-noise-ratio (SNR), the lower bound of secure multicast rate in (17) is very tight, and has little influence on the energy efficiency of the original problem, which can be seen in Section 4. Therefore, we can reduce the computational complexity while obtaining very approximate solutions of the original problem. The equivalent problem is approximately substituted for the original problem as follows,

$$\begin{array}{cc}\hfill \underset{\mathbf{Q}}{argmax}& \frac{W{\left[R_{\mathrm{mc}}-R_{\mathrm{ev},\mathrm{ub}}\right]}^{+}}{P},\hfill \end{array}$$

(18)

$$\begin{array}{cc}\hfill \mathrm{s}.\mathrm{t}.& \mathsf{tr}\left\{\mathbf{Q}\right\}\le P_{\max },\mathbf{Q}⪰\mathbf{0}.\hfill \end{array}$$

(19)

The feasible solution $\mathbf{Q}=\mathbf{0}$ leads the value of $R_{\mathrm{se},\mathrm{lb}}$ to zero. Therefore, all other feasible solutions leading to negative secure multicast rate are not the optimal solutions of the problem in (18). Besides, SNR is defined as $P_{\max }$ in this paper. As a result, we omit the operator ${\left[·\right]}^{+}$ without loss of any optimality, leading the problem equivalent to follows,

$$\begin{array}{cc}\hfill \underset{\mathbf{Q}}{argmax}& \frac{W(R_{\mathrm{mc}}-R_{\mathrm{ev},\mathrm{ub}})}{P},\hfill \end{array}$$

(20)

$$\begin{array}{cc}\hfill \mathrm{s}.\mathrm{t}.& \mathsf{tr}\left\{\mathbf{Q}\right\}\le P_{\max },\mathbf{Q}⪰\mathbf{0}.\hfill \end{array}$$

(21)

We decompose the covariance matrix as $\mathbf{Q}=\mathsf{\Phi }\mathsf{\Lambda }{\mathsf{\Phi }}^H$. The columns of $\mathsf{\Phi }$ and diagonal elements of $\mathsf{\Lambda }$, respectively, denote the eigenvectors and the eigenvalues of $\mathbf{Q}$, which indicates the multicast signaling directions and the powers allocation. We investigate the eigenvectors of the covariance matrix in Theorem 1.

Theorem 1.

Columns of $\mathbf{V}$ are expressed as the eigenvectors of the transmit covariance matrix ${\mathbf{Q}}^{\mathrm{opt}}$ in (20) as follows,

$$\begin{array}{c}\hfill {\mathbf{Q}}^{\mathrm{opt}}=\mathbf{V}\mathsf{\Lambda }{\mathbf{V}}^H.\end{array}$$

(22)

Proof. 

Please refer to the Appendix A. □

According to Theorem 1, we can indicate that the energy efficiency oriented secure multicast transmit signal directions are consistent with the eigenvector of the BS transmission correlation matrix. Therefore, the optimal problem for maximizing the system energy efficiency of secure downlink multicast transmission can be simplified into a beam domain power allocation problem. In particular, the precoding optimization problem of energy efficiency multicast transmission in (20) can be simplified to the following question without loss of optimality,

$$\begin{array}{cc}\hfill \underset{\mathsf{\Lambda }}{argmax}& W\frac{(R_{\mathrm{mc}}\left(\mathsf{\Lambda }\right)-R_{\mathrm{ev},\mathrm{ub}}\left(\mathsf{\Lambda }\right))}{P\left(\mathsf{\Lambda }\right)},\hfill \\ \hfill \mathrm{s}.\mathrm{t}.& \mathsf{tr}\left\{\mathsf{\Lambda }\right\}\le P_{\max },\mathsf{\Lambda }⪰\mathbf{0},\mathsf{\Lambda }\mathrm{diagonal},\hfill \end{array}$$

(23)

where

$$\begin{array}{cc}\hfill R_{\mathrm{mc}}\left(\mathsf{\Lambda }\right)& \triangleq \underset{k}{min}{R}_k\left(\mathsf{\Lambda }\right),\hfill \end{array}$$

(24)

$$\begin{array}{cc}\hfill R_k\left(\mathsf{\Lambda }\right)& \triangleq \mathsf{E}\left\{logdet\left\{{\mathbf{I}}_{N_{\mathrm{r}}}+{\mathbf{G}}_k\mathsf{\Lambda }{\mathbf{G}}_k^H\right\}\right\},\hfill \end{array}$$

(25)

$$\begin{array}{cc}\hfill R_{\mathrm{ev},\mathrm{ub}}\left(\mathsf{\Lambda }\right)& \triangleq logdet\left\{{\mathbf{I}}_{N_{\mathrm{e}}}+{\mathbf{A}}_{\mathrm{ev}}\left(\mathsf{\Lambda }\right)\right\},\hfill \end{array}$$

(26)

$$\begin{array}{cc}\hfill P\left(\mathsf{\Lambda }\right)& \triangleq \zeta \mathsf{tr}\left\{\mathsf{\Lambda }\right\}+M{P}_{\mathrm{c}}+P_{\mathrm{s}}.\hfill \end{array}$$

(27)

In (24), we adopt the minimum transmission rate received by legitimate users into our model to meet the requirement of confidentiality.

3.1. Minorize–Maximize Method

We note that function $R_{\mathrm{mc}}\left(\mathsf{\Lambda }\right)-R_{\mathrm{ev},\mathrm{ub}}\left(\mathsf{\Lambda }\right)$ in (23) is the subtraction of two concave functions. The function is defined as the difference of convex functions (d.c.), which can not be solved by convex optimization directly. Via using the MM method, we adopt an iterative algorithm to deal with it [40,41]. The basic idea is to approximate $R_{\mathrm{ev},\mathrm{ub}}\left(\mathsf{\Lambda }\right)$ by the first-order Taylor series expansion and transform it into a linear function. We can solve the subproblem with further algorithm and get the next iteration. The iterative subproblem is as follows,

$$\begin{array}{cc}\hfill {\mathsf{\Lambda }}^{\left(\ell +1\right)}=\underset{\mathsf{\Lambda }}{argmax}& W\frac{R_{\mathrm{mc}}\left(\mathsf{\Lambda }\right)-\mathsf{tr}\left\{{\left(\frac{\partial }{\partial \mathsf{\Lambda }}{R}_{\mathrm{ev},\mathrm{ub}}\left({\mathsf{\Lambda }}^{\left(\ell \right)}\right)\right)}^T\mathsf{\Lambda }\right\}}{P\left(\mathsf{\Lambda }\right)},\hfill \\ \hfill \mathrm{s}.\mathrm{t}.& \mathsf{tr}\left\{\mathsf{\Lambda }\right\}\le P_{\max },\mathsf{\Lambda }⪰\mathbf{0},\mathsf{\Lambda }\mathrm{diagonal},\hfill \end{array}$$

(28)

where ℓ denotes the iteration index. Note that the gradient of $R_{\mathrm{ev},\mathrm{ub}}\left(\mathsf{\Lambda }\right)$ over $\mathsf{\Lambda }$ can be proved to be a diagonal matrix. We can get its kth diagonal element as follows,   

$$\begin{array}{c}\hfill {\left[\frac{\partial }{\partial \mathsf{\Lambda }}{R}_{\mathrm{ev},\mathrm{ub}}\left({\mathsf{\Lambda }}^{\left(\ell \right)}\right)\right]}_{k,k}=\sum _{i=1}^{N_{\mathrm{e}}}\frac{{\left[{\mathsf{\Omega }}_{\mathrm{ev}}\right]}_{i,k}}{1+{\sum }_{j=1}^M{\left[{\mathsf{\Omega }}_{\mathrm{ev}}\right]}_{i,j}{\left[{\mathsf{\Lambda }}^{\left(\ell \right)}\right]}_{j,j}},\end{array}$$

(29)

where ${\left\{{\mathsf{\Lambda }}^{\left(\ell \right)}\right\}}_{\ell =0}^{\infty }$ is by the MM method in (28), which can converge to a locally optimal solution of the original problem in (20) [42]. Therefore, our next step is to solve the subproblem of each iteration.

3.2. Dinkelbach’s Transform

At the ℓ th iteration of subproblem in (28), the fractional function of the objective is still difficult to solve with convex optimition directly. We have $R_{\mathrm{ee}}\left(\mathsf{\Lambda }\right)=R_{\mathrm{mc}}\left(\mathsf{\Lambda }\right)-\mathsf{tr}\left\{{\left(\frac{\partial }{\partial \mathsf{\Lambda }}{R}_{\mathrm{ev},\mathrm{ub}}\left({\mathsf{\Lambda }}^{\left(\ell \right)}\right)\right)}^T\mathsf{\Lambda }\right\}$, which is a concave-minus-linear function. Thus, $R_{\mathrm{ee}}\left(\mathsf{\Lambda }\right)$ is a concave functions with respect to $\mathsf{\Lambda }$. Meanwhile, $P\left(\mathsf{\Lambda }\right)$ is a linear function. Dinkelbach’s transform is usually used to approach such concave-linear fractional program, whose advantage is that it does not have to include additional constraints. It has been proved that this algorithm possessed superlinear convergence [43,44]. Besides, Dinkelbach’s transform We express the powers allocation $\mathsf{\Lambda }$ solved in (28) as $\mathsf{\Gamma }$. Via invoking Dinkelbach’s transform, we can solve the subproblem of each energy efficiency secure multicast optimization in (28) by iteratively solving the sequence of convex optimization problem as follows,

$$\begin{array}{cc}\hfill {\mathsf{\Gamma }}^{\left(t+1\right)}=\underset{\mathsf{\Gamma }}{argmax}& W{R}_{\mathrm{mc}}\left(\mathsf{\Gamma }\right)-W\mathsf{tr}\left\{{\left(\frac{\partial }{\partial \mathsf{\Lambda }}{R}_{\mathrm{ev},\mathrm{ub}}\left({\mathsf{\Lambda }}^{\left(\ell \right)}\right)\right)}^T\mathsf{\Gamma }\right\}-{\eta }^{\left(t\right)}P\left(\mathsf{\Gamma }\right),\hfill \\ \hfill \mathrm{s}.\mathrm{t}.& \mathsf{tr}\left\{\mathsf{\Gamma }\right\}\le P_{\max },\mathsf{\Gamma }⪰\mathbf{0},\mathsf{\Gamma }\mathrm{diagonal},\hfill \end{array}$$

(30)

where t denotes the iteration index. ${\mathsf{\Lambda }}^{(\ell )}$ is the solution of the ℓ th outer iteration of the MM method. Note that it is treated as a constant and will not update in the inner iteration in (30). ${\eta }^{\left(t\right)}$ denotes an auxiliary variable introduced by Dinkelbach’s transform, which can be updated as follows,

$$\begin{array}{c}\hfill {\eta }^{\left(t\right)}=W\frac{R_{\mathrm{mc}}\left({\mathsf{\Gamma }}^{\left(t\right)}\right)-\mathsf{tr}\left\{{\left(\frac{\partial }{\partial \mathsf{\Lambda }}{R}_{\mathrm{ev},\mathrm{ub}}\left({\mathsf{\Lambda }}^{\left(\ell \right)}\right)\right)}^T{\mathsf{\Gamma }}^{\left(t\right)}\right\}}{P\left({\mathsf{\Gamma }}^{\left(t\right)}\right)}.\end{array}$$

(31)

By optimizing $\mathsf{\Gamma }$ and $\eta$ alternately, the generated solution sequence ${\left\{{\mathsf{\Gamma }}^{\left(t\right)}\right\}}_{t=0}^{\infty }$ in (30) is also proved to converge and get the optimal solution of the subproblem in (28).

3.3. Deterministic Equivalent

Due to a large number of samples need to be averaged during the expectation operation in (30), the optimization complexity is still high, although each subproblem in (28) is a concave program. Thus, without using Monte Carlo averaging channel realizations, we calculate the deterministic equivalent of the secure transmit rate by adopting the large-dimensional random matrix theory, which can reduce the implementation complexity [45,46]. Note that the method of deterministic equivalent can keep the convexity of the transmission rate [47]. Therefore, the concave subproblem in (28) is given as

$$\begin{array}{cc}\hfill {\mathsf{\Gamma }}^{\left(t+1\right)}=\underset{\mathsf{\Gamma }}{argmax}& \underset{k}{min}W{\overline{R}}_k\left(\mathsf{\Gamma }\right)-W\mathsf{tr}\left\{{\left(\frac{\partial }{\partial \mathsf{\Lambda }}{R}_{\mathrm{ev},\mathrm{ub}}\left({\mathsf{\Lambda }}^{\left(\ell \right)}\right)\right)}^T\mathsf{\Gamma }\right\}-{\eta }^{\left(t\right)}P\left(\mathsf{\Gamma }\right),\hfill \\ \hfill \mathrm{s}.\mathrm{t}.& \mathsf{tr}\left\{\mathsf{\Gamma }\right\}\le P_{\max },\mathsf{\Gamma }⪰\mathbf{0},\mathsf{\Gamma }\mathrm{diagonal}.\hfill \end{array}$$

(32)

In (32), ${\overline{R}}_k\left(\mathsf{\Gamma }\right)$ denotes the deterministic equivalent expression of $R_k\left(\mathsf{\Gamma }\right)$, which is expressed as

$$\begin{array}{cc}\hfill {\overline{R}}_k\left(\mathsf{\Gamma }\right)& =logdet\left\{{\mathbf{I}}_M+{\mathsf{\Xi }}_k\mathsf{\Gamma }\right\}+logdet\left\{{\tilde{\mathsf{\Phi }}}_k\right\}-\mathsf{tr}\left\{{\tilde{\mathsf{\Xi }}}_k{\tilde{\mathsf{\Phi }}}_k^{-1}\right\},\hfill \end{array}$$

(33)

where ${\mathsf{\Xi }}_k\in {\mathbb{C}}^{M\times M}$, ${\tilde{\mathsf{\Xi }}}_k\in {\mathbb{C}}^{N_{\mathrm{r}}\times N_{\mathrm{r}}}$, and ${\tilde{\mathsf{\Phi }}}_k\in {\mathbb{C}}^{N_{\mathrm{r}}\times N_{\mathrm{r}}}$ can be obtained by solving the following three fixed-point equations,

$$\begin{array}{cc}\hfill {\mathsf{\Xi }}_k& ={\mathbf{B}}_k\left({\tilde{\mathsf{\Phi }}}_k^{-1}\right),\hfill \end{array}$$

(34a)

$$\begin{array}{cc}\hfill {\tilde{\mathsf{\Xi }}}_k& ={\mathbf{C}}_k\left(\mathsf{\Gamma }{\left({\mathbf{I}}_M+\mathsf{\Gamma }{\mathsf{\Xi }}_k\right)}^{-1}\right),\hfill \end{array}$$

(34b)

$$\begin{array}{cc}\hfill {\tilde{\mathsf{\Phi }}}_k& ={\mathbf{I}}_{N_{\mathrm{r}}}+{\tilde{\mathsf{\Xi }}}_k,\hfill \end{array}$$

(34c)

where ${\mathbf{B}}_k\left(\mathbf{X}\right)\triangleq \mathsf{E}\left\{{\mathbf{G}}_k^H\mathbf{X}{\mathbf{G}}_k\right\}\in {\mathbb{C}}^{M\times M}$ and ${\mathbf{C}}_k\left(\mathbf{X}\right)\triangleq \mathsf{E}\left\{{\mathbf{G}}_k\mathbf{X}{\mathbf{G}}_k^H\right\}\in {\mathbb{C}}^{N_{\mathrm{r}}\times N_{\mathrm{r}}}$ are matrix-valued functions, which can output matrices in diagonal forms, the corresponding ith diagonal elements are respectively obtained by

$$\begin{array}{cc}\hfill {\left[{\mathbf{B}}_k\left(\mathbf{X}\right)\right]}_{i,i}& =\mathsf{tr}\left\{\mathsf{diag}\left\{{\left[{\mathsf{\Omega }}_k\right]}_{:,i}\right\}\mathbf{X}\right\},\hfill \end{array}$$

(35)

$$\begin{array}{cc}\hfill {\left[{\mathbf{C}}_k\left(\mathbf{X}\right)\right]}_{i,i}& =\mathsf{tr}\left\{\mathsf{diag}\left\{{\left({\left[{\mathsf{\Omega }}_k\right]}_{i,:}\right)}^T\right\}\mathbf{X}\right\}.\hfill \end{array}$$

(36)

Via the deterministic equivalent approach instead of exhaustive averaging, the expression ${\overline{R}}_k\left(\mathsf{\Gamma }\right)$ is efficiently calculated by adopting the channel statistics ${\mathsf{\Omega }}_k$ in a small number of fixed-point iterations. As a result, the computational complexity to optimize the convex program in (28) can be reduced significantly.

3.4. Discussion

In the case of small numbers of antennas, the deterministic equivalent expression ${\overline{R}}_k\left(\mathsf{\Gamma }\right)$ and $R_k\left(\mathsf{\Gamma }\right)$ are approximately equal [45,46]. Besides, with respect to $\mathsf{\Gamma }$, ${\overline{R}}_k\left(\mathsf{\Gamma }\right)$ is still observed to be concave from (33). Thus, each subproblem in (32) does not change the convexity and still guarantees the convergence. Finally, we obtain our algorithm by two-level iteration. The optimal power allocation is obtained by the inner iteration and assigned to the outer iteration. Then, the initial power allocation of the inner layer is updated after the outer iteration. The proposed energy efficiency optimization algorithm of massive MIMO secure multicast transmission is described in detail in Algorithm 1.

Algorithm 1 Energy Efficiency Optimization Algorithm for Secure Multicast Transmission

Input: The beam domain channel transmit statistics CSI ${\mathsf{\Omega }}_k$ and ${\mathsf{\Omega }}_{\mathrm{ev}}$, initial power allocation ${\mathsf{\Lambda }}^{\left(0\right)}$, the threshold $ϵ$

Output: Optimal power allocation matrix $\mathsf{\Lambda }$

1: Initialization: $\ell =0$

2: Calculate ${\mathsf{\Lambda }}^{\left(0\right)}$ using (26), (27), and (33)

3: repeat

4:    Update $\ell \leftarrow \ell +1$

5:    Initialization: $t=0$, ${\mathsf{\Gamma }}^{\left(0\right)}={\mathsf{\Lambda }}^{(\ell -1)}$, ${\eta }^{\left(0\right)}$ using (31)

6:    repeat

7:     Update $t\leftarrow t+1$

8:     Calculate ${\mathsf{\Gamma }}^{\left(t\right)}$ by using (30) with ${\eta }^{\left(t-1\right)}$

9:     Update ${\eta }^{\left(t\right)}$ using (31)

10:    until $\parallel {\mathsf{\Gamma }}^{\left(t\right)}-{\mathsf{\Gamma }}^{(t-1)}\parallel \le ϵ$

11:    Return ${\mathsf{\Lambda }}^{\left(\ell \right)}={\mathsf{\Gamma }}^{\left(t\right)}$

12: until$\parallel {\mathsf{\Lambda }}^{(\ell )}-{\mathsf{\Lambda }}^{(\ell -1)}\parallel \le ϵ$

13: Return $\mathsf{\Lambda }={\mathsf{\Lambda }}^{\left(\ell \right)}$

4. Simulation Results

To illustrate the performance of the energy efficiency optimization for massive MIMO secure multicast transmission with known statistical CSI. A series of assumptions about the simulation environment, transmission channel and antenna configuration are made in our work [48]. The major parameters related to simulation, such as the power consumption, are listed in

Table 1. Major simulation parameters of massive multiple-input multiple-output (MIMO) secure multicast transmission.

Parameter	Value
Channel model	3GPP SCM
Array topology	ULA with half wavelength antenna spacing
Propagation scenario	Suburban macro-cell
The numbers of BS antennas	M = 64
The numbers of UTs	K = 8
The numbers of UTs antennas	$N_{\mathrm{r}}$ = 4 ($\forall k$)
The numbers of eavesdropper antennas	$N_{\mathrm{e}}$ = 4
Bandwidth	W = 10 MHz
Amplifier inefficiency factor	$\xi$ = 5
Static power consumption	$P_{\mathrm{s}}$ = 40 dBm
Hardware dissipated power per antenna	$P_{\mathrm{c}}$ = 30 dBm

.

First of all, we evaluate our proposed algorithm convergence performance. In order to observe the convergence of the whole algorithm, we choose the outer iteration to draw Figure 1, which plots energy efficiency versus the number of iterations under different SNR. From Figure 1, it can be observed that our algorithm exhibits fast convergence performance. Besides, the lower the SNR regime is, the quicker the iteration converges. Thus, we can indicate that power allocation might perform better under the high SNR regime on optimizing secure multicast energy efficiency. In the regime of Figure 1, with the SNR increasing, the gain of energy efficiency increases faster, which shows that the secure multicast rate is growing rapidly.

Next, we evaluate the energy efficiency performance of the proposed algorithm with a lower bound of secure multicast rate and the original optimization problem with the secure rate obtained by the Monte Carlo Method. It is shown in Figure 2 that the solutions of the two are quite tight. Therefore, we can handle the approximation problem to obtain the desired energy efficiency but significantly reduce the algorithm complexity. Then, we compare the energy efficiency performance optimized by the proposed algorithm with the traditional rate maximization one. In addition, the deterministic equivalent results are compared with those of the Monte Carlo method. According to Figure 2, it is not difficult to find that the gap between the deterministic equivalent and the Monte Carlo methods is almost zero. Besides, under the low transmit power regime, our proposed energy efficiency maximization and the traditional rate maximization algorithms have similar energy efficiency performance. It is also proved that the inference in Figure 1 is correct. However, under the high transmit power regime, the energy efficiency maximization algorithm is obviously superior to the traditional rate maximization one. It also shows that we can find a compromise power allocation to meet the two requirements at the same time in the goal of maximum rate and maximum energy efficiency for massive MIMO secure multicast transmission.

Finally, we study the energy efficiency performance of our algorithm versus the numbers of antennas at the BS and show the results in Figure 3. With the number of BS antennas increasing, the change of energy efficiency performance tends to be flat, which shows that the power consumption dominates the performance of energy efficiency when the number of BS antennas is large. Besides, it can be indicated that we can find a suitable power allocation to meet our requirement of high energy efficiency. At the same time, we can also find that the energy efficiency does not increase with the number of antennas. This is because under large antennas and low power constraint, the secure rate changes faster than power consumption at the BS with the increase of the maximum power constraint, while the dynamic power consumption of the BS antennas increases rapidly under the high power constraint, leading the situation on the contrary. Finally, when the number of antennas increases, the system needs more power to reach better energy efficiency. Besides, the energy efficiency performance of the Monte Carlo method is quite close to the proposed one, which confirms the accuracy of the analysis.

5. Conclusions

We studied the energy efficiency optimization problem in downlink secure multicast transmission for massive MIMO with known statistical CSI in this paper. The energy efficiency maximization problem was formulated as maximizing the ratio of the secure multicast rate to the power consumption. We then obtained the secure multicast rate by introducing its lower bound to simplify the objective problem. Via making the transmit signaling directions into a closed-form, we simplified the energy efficiency optimization of matrix-valued strategy design to a beam domain power allocation problem. Next, we designed an algorithm including two-level iteration in the beam domain, and guaranteed the algorithm to converge to a locally optimal point based on the MM method and Dinkelbach’s transform. In addition, via substituting the multicast rate for its deterministic equivalent, the complexity of computational was significantly reduced. We made extensive simulations to verify the effectiveness of our algorithm. A large number of simulation results showed that the algorithm has fast convergence and low complexity. Compared with the rate maximization method, the proposed algorithm had more efficient and stable performance under the high transmission power constraint.