MIMO with $1$-bit Pre/Post-Coding Resolution: A Quantum Annealing Approach

Afundamental enabling technology for future 6G wireless communication systems is massive/large multiple-input multiple-output (MIMO) systems [1]. By exploiting the spatial degrees of freedom associated with the antenna arrays and through appropriate signal processing, MIMO systems are able to support the extremely high requirements of capacity, throughput and reliability. However, as the number of antennas increases, the implementation complexity/cost of conventional digital MIMO in terms of radio frequency (RF) chains, digital-to-analogue converters (DACs) and/or analogue-to-digital converters (ADCs), base-band signal processing etc. becomes a bottleneck for practical applications [2]. To address this challenge, analogue/digital hybrid MIMO architectures which split the processing into digital and analogue domain to balance the trade-off between performance and implementation complexity/cost have emerged as a promising solution [2]. Other relevant techniques consider analogue-only signal processing where low-resolutions phase shifters are used to implement analogue pre/post-coding (mainly suitable for high frequency bands e.g., millimeter wave), or the employment of low-resolution DACs/ADCs [3, 4].

In contrast to existing approaches, in this work, we study a basic (digital) point-to-point MIMO system with low-resolution pre/post-coding processing; specifically, we assume that both the transmitter and receiver employ $1$-bit signal processing resolution ($1$-bit per complex dimension) in order to extremely reduce complexity and power consumption¹¹1It is worth noting that although $1$-bit processing at the transmitter is similar to $1$-bit DAC, the assumption of $1$-bit processing at the receiver is fundamentally different to conventional $1$-bit ADC; this assumption, makes the considered architecture different from the existing $1$-bit DAC/ADC MIMO architectures [4].. This new MIMO architecture is promising for future low-power/low-computation devices in the era of the Internet of Things. Due to the complex-valued binary resolution of the new MIMO architecture, the design of pre/post-coding vectors becomes an NP-hard problem; it can be solved optimally through exhaustive search (ES) over all the possible transmit/receive configurations corresponding to an exponential complexity. To the best of the author’s knowledge, the considered MIMO design problem is reported for the first time in this paper.

Adiabatic quantum computing is a promising tool to solve NP-hard problems by using the principles of the Adiabatic Theorem [5]. Specifically, if a quantum system is initially in the ground state (quantum state with the lowest energy) of an initial Hamiltonian and the system evolves/changes slowly (gradually), it will converge to the ground state of the final Hamiltonian. By encoding the final Hamiltonian to the desired problem, the Adiabatic evolution can be used to solve complex combinatorial optimization problems which are represented as instances of the Ising model (equivalent to quadratic unconstrained binary optimization (QUBO) problem) [6]. D-WAVE is a commercial analogue quantum device that implements a noisy approximation of the quantum Adiabatic algorithm called quantum annealing (QA) [5, 7]. This quantum device has received considerable interest lately due to the high number of available qubits (more than $5,000$ qubits in latest hardware architectures) and its friendly interface for remote access. It is worth noting that in QA, instead of analyzing the computation time/complexity of a given algorithm, we mainly study the trade-off between time and the probability that the output is correct.

D-WAVE QA is introduced as a promising computation tool to solve complex (NP-hard) combinatorial problems in wireless communication system [8]. The associated literature mainly focuses on the solution of the maximum-likelihood detection problem in large multi-user MIMO systems [9, 10]; recent research studies employ D-WAVE QA in other wireless communication problems such as channel decoding (polar codes) [11], passive beamforming design in intelligent reflecting surfaces [12], configuration selection in reconfigurable/fluid-antenna systems [13], beam assignment in satellite systems [14] etc.. The integration of D-WAVE QA in the design of wireless communication systems is a new research area of paramount importance.

In this paper, we extend our work in [15] towards the design of $1$-bit digital (complex-valued) pre/post coding vectors. Specifically, we solve the problem of pre/post-coding design in a MIMO system with $1$-bit resolution per complex dimension. We focus on the maximization of the received signal-to-noise ratio (SNR) which due to the complex binary assumption relies to the solution of an NP-hard complex-valued combinatorial problem. By using the principles of alternating optimization and the computation capabilities of the D-WAVE QA, an iterative algorithm is proposed which solves appropriate QUBO instances at each iteration to gradually compute the pre/post-coding vectors. An appropriate mathematical framework that converts the original complex-valued combinatorial problem to real-valued QUBO instances is investigated. Experimental results in a state-of-the-art D-WAVE QA device show that the proposed iterative algorithm achieves near-optimal performance (similar to ES) while ensuring polynomial complexity.

To further improve the quality of the returned distinct D-WAVE QA solutions, a new pre-processing technique that preserves the quadratic QUBO matrix from the detrimental effects of the Hamiltonian noise is proposed. The proposed pre-processing technique is motivated by the non-uniform compression in pulse code modulation (PCM) communication systems. Specifically, motivated by the experimental observation that small magnitude elements ($\approx 0$) of the quadratic matrix are more sensitive to the Hamiltonian noise [16] (similar to quantization noise in PCM [17]), a non-linear companding ($\mu$-law) pre-processing is proposed to boost small-magnitude elements without qualitatively affecting the structure of the original QUBO problem. Experimental results show that the proposed pre-processing technique improves the quality of the D-WAVE solutions (we have less distinct solutions of higher quality) and increases the probability of the optimal solution.

The remainder of this paper is organized as follows. The system model and the problem formulation are described in Section II. Section III presents the proposed QA-based iterative algorithm for the design of the $1$-bit pre/post-coding vectors. The QUBO pre-processing technique that improves the efficiency of the D-WAVE QA is presented in Section IV. Pertinent experimental examples are presented in Section V. Finally, concluding remarks are discussed in the last section.

Notation: Boldface letters refer to vectors (lower case) or matrices (upper case); $(\cdot)^{T}$, $(\cdot)^{H}$, $\|\cdot\|$, and $\|\cdot\|_{\max}$ denote transpose, complex conjugate transpose, Euclidean norm, and max-norm, respectively; $\textrm{diag}(\boldsymbol{x})$ denotes a diagonal matrix whose main diagonal is $\boldsymbol{x}$; $j=\sqrt{-1}$ is the imaginary unit; $\mathcal{R}(\cdot)$ and $\mathcal{I}(\cdot)$ return the real and imaginary part of their complex arguments, respectively; $\boldsymbol{x}(k:n)$ refers to the vector with elements $x_{k},x_{k+1},\ldots,x_{n}$; the space of $x\times y$ complex matrices is denoted by $\mathbb{C}^{x\times y}$, $\mathbbm{1}$ denotes an all-ones column vector of appropriate dimension; $\textrm{sgn}(\cdot)$ denotes the sign function; $\mathcal{CN}(\mu,\sigma^{2})$ denotes complex Gaussian distribution with mean $\mu$ and variance $\sigma^{2}$.

We consider a point-to-point MIMO system consisting of $N_{T}$ and $N_{R}$ transmit and receive antennas, respectively [18]. Let $\boldsymbol{H}\in\mathbbm{C}^{N_{R}\times N_{T}}$ denote the MIMO channel matrix, where entries correspond to the channels between the transmit and and receive antennas i.e., $h_{i,j}$ is the channel coefficient between the $j$-th transmit antenna and the $i$-th receive antenna. Without loss of generality, we assume normalized Rayleigh block fading channels i.e., $h_{i,j}\sim\mathcal{CN}(0,1)$.

The MIMO system operates in the beamforming mode [18] (single data flow) to maximize the received SNR and thus digital pre/post- processing vectors are adjusted accordingly. To reduce complexity and power consumption, both the transmitter and the receiver employ vector processing of $1$-bit resolution (i.e., $1$-bit for the real part and $1$-bit for the imaginary part). Fig. 1 schematically presents the system model.

Let $\boldsymbol{f}\in\mathcal{S}^{N_{T}\times 1}$ and $\boldsymbol{g}\in\mathcal{S}^{N_{R}\times 1}$ denote the (unnormalized) complex $1$-bit pre-coding and post-coding vectors, respectively, taking values in the quantized complex set $\mathcal{S}\triangleq\{\pm 1\pm 1j\}$; the normalized pre/post-coding vectors are given by $\boldsymbol{f}/\|\boldsymbol{f}\|=\boldsymbol{f}/\sqrt{2N_{T}}$ and $\boldsymbol{g}/\|\boldsymbol{g}\|=\boldsymbol{g}/\sqrt{2N_{R}}$, respectively. We assume a perfect channel state information at both the transmitter and the receiver. The received SNR is given by

where $\sigma^{2}$ is the variance of the additive white Gaussian noise, $P$ is the transmit power, while the extra terms in the denominator are due to the power normalization of the pre/post-coding vectors. Since the objective of the MIMO system is to maximize the received SNR, we introduce the following design problem

Due to the binary nature of the complex pre/post-coding vectors, the above combinatorial problem is NP-hard; the optimal solution requires ES over all the possible pre/post-coding vectors.

We propose a new algorithm that incorporates a QA solver in the design of the $1$-bit complex pre/post-coding vectors. By taking into account the binary combinatorial structure of the problem, we introduce an alternating optimization algorithm that solves a QUBO problem at each iteration. The QUBO formulations are solved through a state-of-the-art D-WAVE QA device [7]. More specifically, based on the equivalent expression in (2b), we firstly fix the vector $\boldsymbol{g}$ (by using a random solution in the initial step) and we study the following combinatorial problem with respect to the pre-coding vector $\boldsymbol{f}$ i.e.,

Given that the involved matrices/vectors in (3) are complex valued (which it is not compatible with current QA solvers), we transform the problem in (3) to an equivalent real-valued form

where $\mathcal{P}\triangleq\{+1,-1\}$ denotes the spin set. Then, we convert the pre-coding vector $\boldsymbol{f}_{r}$ (containing spin variables $\{+1,-1\}$) to the binary vector $\boldsymbol{b}_{f}$ (with entries in the set $\{0,1\}$) by using the transformation [9] $\boldsymbol{b}_{f}=\frac{1}{2}(\boldsymbol{f}_{r}+\mathbbm{1})$ and thus the problem in (4a) is converted to a standard QUBO formulation

The last mathematical step transforms the above (symmetric) quadratic matrix $\boldsymbol{V}_{0}$ to a form that is compatible with the D-WAVE QA solver [10]; we normalize the matrix such as all its entries take values in the range $[-1,+1]$. The calibrated matrix is written by $\boldsymbol{V}_{n}=\frac{\boldsymbol{V}_{0}}{\|\boldsymbol{V}_{0}\|_{\max}}$ and thus the D-WAVE compatible QUBO problem is written as

The returned binary vector $\boldsymbol{b}_{f}$ is converted back to the spin vector $\boldsymbol{f}$. By using the principles of alternating optimization, we then fix $\boldsymbol{f}$ (by using the previous solution), and we solve the equivalent problem in (2c) i.e.,

which is transformed to an equivalent real-valued spin form

The above spin-variable formulation is converted to a QUBO formulation by using similar analytical steps like before

where $\boldsymbol{b}_{g}$ is the binary representation of the vector $\boldsymbol{g}$, while the normalized matrix $\boldsymbol{R}_{n}$ is defined similar to $\boldsymbol{V}_{n}$.

The above alternating optimization process is repeated until convergence or a maximum number of iterations ($K$ iterations) is achieved. Due to the discrete/binary nature of the problem, the final solution is sensitive to the initial conditions (i.e., initial vector $\boldsymbol{g}^{(0)})$ which could result in convergence to a local maximum (SNR objective function). To overcome this limitation and facilitate the iterative algorithm to escape local maxima and get closer to the optimal solution, we introduce $L$ independent execution of the algorithm by using different initial conditions ($\boldsymbol{g}_{l}^{(0)}$ with $l=1,\ldots,L$); from the $L$ returned solutions, we keep the one corresponding to the maximum SNR objective function. The pseudocode of the QA-based iterative algorithm is given in Algorithm 1.

The time complexity of the Algorithm 1 mainly refers to the number of anneals, which is a system parameter that is tuned empirically. Obviously, as the size of the problem increases, the number of anneals should also increase (in a polynomial way) in order to successfully determine the optimal solution. Specifically, the time complexity of the proposed QA scheme can be written as $K\times L\times(2\times N_{a})\times T_{a}$ (time-units) where $N_{a}$ is the number of anneals, $T_{a}$ is the anneal time, while the factor $2$ is due to the fact that we solve two QUBO problems at each basic algorithmic iteration. It is obvious that the time complexity of the proposed scheme is polynomial with respect the main system parameters and thus much lower than the exponential time complexity of the classical ES scheme.

In this section, we study a new pre-processing QA technique that converts the quadratic QUBO matrix in a form that is more robust to Hamiltonian noise [19]. Due to the interaction of the D-WAVE QA device with the external environment as well as due to hardware limitations, the QA process suffers from noise effects that decrease the quality of the solution. These noise effects are collectively represented by the term integrated control errors (ICE) which captures all the sources of solution infidelity (e.g., cross-talk with adjacent qubits, temperature fluctuations, flux noise, quantization/precision, bandwidth of the analogue circuits etc.). An informative overview of the ICE sources is given in [16, Ch. 3]. Due to the ICE effects, the QA device solves a slightly altered QUBO problem given by

where $\boldsymbol{Q}$ is the original symmetric quadratic QUBO matrix, while the random matrix $\boldsymbol{E}$ represents the Hamiltonian/ICE noise in the bias and couplers coefficients (the elements of the matrix $\boldsymbol{E}$ are random variables e.g., coloured Gaussian distribution [16]). Therefore, it is obvious that the QA device solves a modified/noisy problem whose ground state may be different that the one of the original QUBO problem. The associated literature mainly focuses on the characterization of the noise by using statistical signal processing tools [19] or the investigation of heuristic post-processing techniques that improve the quality of the returned solution [20]. In practice, to combat the noise/ICE effects and overcome potential local minima, the same QUBO instance is solved multiple times (each execution is called anneal) in a row and the best solution among all the anneals is adopted²²2It is worth noting that this concept is very similar to retransmission techniques (e.g., repetition code) in wireless communications systems, where the same information signal is transmitted multiple times through independent noisy channels to achieve diversity and increase reliability [21, Ch. 3]..

From the modified QUBO problem in (12), it is obvious that the elements of the original matrix $\boldsymbol{Q}$ with small values/magnitudes ($\approx 0$) are more sensitive to noise effects. Our proposed pre-processing QUBO technique converts the $\boldsymbol{Q}$ matrix elements with small magnitude to higher values without affecting the qualitative interconnection between the qubits (i.e., the sign and the order in the strength of the bias/couplers coefficients remain the same). Specifically, inspired by the non-uniform quantization process in PCM systems which boosts the lower amplitudes and reduces the dynamic range of the input speech signal [17, Ch. 7], we introduce a nonlinear compression block (named compander by using the terminology from PCM systems) to pre-process $\boldsymbol{Q}$ elements before the D-WAVE QA operation. Without loss of generality, we consider the $\mu$-low (logarithmic) compander which is used in PCM systems in the United States and Canada, defined as follows [17, Ch. 7]

The positive parameter $\mu$ controls the amount of compression. Fig. 2 plots the characteristic function of the $\mu$-law for various values of $\mu$; it can be seen that as the $\mu$ parameter increases, the low magnitude elements are mapped to much higher values. In our experimental studies (without loss of generality), we consider $\mu=255$ which is the value that is used in commercial PCM systems. After QUBO pre-processing, the D-WAVE QA device solves the following noise-robust QUBO problem

Fig. 3 schematically presents the proposed QA pre-processing technique to combat ICE effects; it is worth noting, that the proposed approach is generic and be used with other physics/quantum-inspired heuristic solvers (e.g. coherent Ising machines).

Experimental results in a state-of-the-art D-WAVE QA device are carried-out in order to validate the efficiency of the proposed techniques.

In this paper, we studied a new low-complexity digital MIMO architecture with $1$-bit (processing) resolution per complex dimension. The design of pre/post-coding that maximizes the received SNR relies to an NP-hard combinatorial problem with exponential complexity for the optimal ES scheme. An iterative D-WAVE QA-based algorithm that solves appropriate real-valued QUBO instances at each iteration is proposed; the iterative scheme gradually achieves near-optimal (similar to ES) performance while ensuring polynomial complexity. To further boost D-WAVE efficiency, a pre-processing non-linear companding technique that transforms the quadratic matrix to a form that is more robust to Hamiltonian/ICE noise is also investigated. Experimental results show that the new pre-processing technique improves the quality of the D-WAVE solutions and could decrease the required number of the total anneals.