Tri-Domain Multiuser MIMO Precoding Optimization and Channel Estimation with Spatial-EM Reconfigurable Antenna

For EMRA architectures, one appealing benefit is the ability to adapt the radiation field [Hasan2018TWC, Bahceci2017TWC, Wang2023TWC]. By electronically tuning loads or switch states, the current distribution of each antenna element can be changed, thereby reshaping its radiation pattern. At the system level, this effect can be modeled through antenna-mode selection or basis-function expansions whose coefficients serve as electromagnetic (EM)-domain precoders [Zhao2021TWCMode, Wang2023TWC, Ying2025TCOM]. In this setting, transmitter design depends not only on conventional spatial-domain channel state information (sCSI), but also on EM-domain channel state information (eCSI), which links the propagation geometry to the radiation-pattern basis and enables joint EM-domain and digital-domain precoding for substantial SE improvement over TFA systems [Wang2023TWC, Ying2025TCOM].

One core challenge for EMRA systems is channel estimation (CE). Recent studies have also examined channel state information (CSI) extrapolation across radiation states [Liang2024TVTExtrap]. In wideband parametric CE pipelines, the channel is first characterized by multipath delays, angle-of-departure (AoD), angle-of-arrival (AoA), and equivalent path gains, and the eCSI required for EM-domain design is then reconstructed from these geometric parameters [Bahceci2017TWC, Gao2014LCOMM, Liao2019TCOM, Ying2025TCOM]. Although this representation reduces the estimation dimension, its accuracy still hinges on spatial sampling. In orthogonal frequency division multiplexing (OFDM) systems, frequency diversity mainly assists delay estimation, whereas reliable AoD recovery still requires sufficiently rich spatial observations and is further affected by spatial-wideband and frequency-wideband effects [Wang2018TSPWideband]. For fixed-geometry arrays, simply enlarging the inter-element spacing is not always effective. Although a larger aperture can improve angular resolution, excessive spacing may also introduce spatial aliasing and grating-lobe ambiguity [Gao2014LCOMM, Liao2019TCOM]. Since AoD errors directly perturb the steering and basis matrices used to assemble eCSI, they propagate to EM-domain precoding and degrade the subsequent transmission design [Costa2010TSPSpherical, Nadeem2015TSP3D, Ying2025TCOM].

Work on SMA architectures now spans channel modeling, antenna-position optimization, and multiuser precoding for MA and FAS systems [Zhu2024TWC, Yang2025TSPFWMMSE, Xiao2024TWCJointPositioning, Zhu2023LCOMMNullSteering], and recent surveys provide broader accounts of these developments [Zhu2025COMSTTutorialMA, New2025COMST]. Antenna translation can also be combined with rotation to provide position and orientation DoFs [Ma2026COMSTSurveyRMA]. Linear minimum mean square error and Bayesian estimators have been studied for CE [Skouroumounis2023TCOM, Zhang2025TWC]. Compressed sensing (CS) methods recover AoDs, AoAs, and path gains from measurements at selected MA positions [Ma2023LCOMMCSMA, Xiao2024TWCCSMA]; with an EM radiation basis, these parameters can be used to assemble eCSI. For these methods, the angular estimates depend on the grid resolution and sensing geometry. In the proposed design, coordinated motion forms a uniform virtual array for two-dimensional estimation of signal parameters via rotational invariance techniques (ESPRIT), and the recovered parameters are used to assemble eCSI at the transmission positions.

Recent studies combine spatial and EM reconfiguration through electronically synthesized radiation locations [Ma2026COMSTSurveyRMA]. Pixel-based reconfigurable antennas emulate FAS ports on a fixed aperture [Zhang2025OJAPPRAFAS], and the reconfigurable pixel antenna (RPA)-based electronic movable-antenna array (REMAA) selects from a finite set of candidate radiation positions for multiuser transmission [Chen2025TCOMREMAA]. Radiation-center reconfigurable arrays similarly combine discrete radiation-center selection with analog and digital beamforming [Li2026TWCRCRAA]. These architectures provide electronic spatial flexibility through effective positions or radiation centers on a fixed aperture.

Physical repositioning directly controls the sampling coordinates and supplies a regular observation geometry for virtual-array CE. The EM coefficients remain available for radiation-pattern adaptation during data transmission, so sensing and transmission retain separate control variables [Zhang2025OJAPPRAFAS, Chen2025TCOMREMAA, Li2026TWCRCRAA].

The spatial-EM reconfigurable antenna (SEMRA) architecture considered here coordinates antenna repositioning across pilot slots to form a denser virtual array under the same per-user pilot overhead as EMRA. The resulting parametric eCSI model supports tri-domain precoding over antenna positions, EM-domain pattern weights, and digital precoders.

This paper presents a SEMRA framework that combines movement-aided parametric CE with tri-domain precoding. Antenna repositioning supplies the spatial samples for virtual-array eCSI reconstruction and the position variables for transmission optimization, while EM-domain pattern weights and digital precoders are optimized for data transmission. The contributions are summarized as follows.

• •

  Within the proposed SEMRA framework, we develop a concrete low-overhead movement-aided parametric CE procedure. By coordinating local repositioning of the $M_{y}\times M_{z}$ antennas over a prescribed multi-block pilot schedule, a denser virtual array is synthesized while preserving the same per-user pilot overhead as the EMRA baseline. The resulting virtual array supplies additional spatial samples for the parametric estimation pipeline, improves AoD estimation, and yields a position-independent parametric channel model from which the eCSI can be assembled at desired antenna positions without additional pilots.

• •

  We formulate the SEMRA SE-maximization problem and develop a tri-domain alternating optimization framework based on block coordinate descent (BCD) over antenna positions, EM-domain pattern weights, and digital precoders. A zero-forcing (ZF)-based joint design is first constructed as a low-complexity baseline to isolate the gain of spatial reconfiguration beyond existing EM-digital designs for EMRA systems. We then develop a weighted minimum mean square error (WMMSE)-based design to further improve the system SE.

• •

  Numerical results show that the proposed CE scheme improves eCSI accuracy over the EMRA baseline under the same per-user pilot overhead. Within the tested range, both SEMRA-ZF and SEMRA-WMMSE achieve higher SE than the EMRA baseline under both estimated and perfect eCSI scenarios. For the tested multiuser loads, SEMRA-WMMSE further improves over SEMRA-ZF. The SEMRA-over-EMRA advantage is especially evident at larger inter-element spacings, where the EMRA baseline is more sensitive to spatial aliasing.

The remainder of this paper is organized as follows. Section LABEL:S2 presents the system model and formulates the tri-domain SE-maximization problem. Sections LABEL:S3 and LABEL:S4 develop the SEMRA-ZF and SEMRA-WMMSE designs, respectively. Section LABEL:S5 introduces the movement-aided parametric CE procedure. Section LABEL:S6 reports the simulation results, and Section LABEL:S7 concludes the paper.

Notation: Matrices and column vectors are denoted by uppercase and lowercase boldface letters, respectively. $(\cdot)^{T}$, $(\cdot)^{*}$, $(\cdot)^{H}$, $(\cdot)^{-1}$, $(\cdot)^{\dagger}$, and $\mathbb{E}\{\cdot\}$ denote the transpose, conjugate, Hermitian transpose, inversion, Moore-Penrose pseudoinverse, and expectation. $[\bm{x}]_{k}$ and $[\bm{A}]_{i,j}$ denote the $k$-th entry of vector $\bm{x}$ and the $(i,j)$-th entry of matrix $\bm{A}$. $\mathrm{diag}(\bm{a})$ forms a diagonal matrix from vector $\bm{a}$, and $\mathrm{Blkdiag}\{\bm{a}_{1},\ldots,\bm{a}_{M}\}$ constructs a block diagonal matrix. $\odot$ denotes the Hadamard product. $\|\cdot\|_{2}$ and $\|\cdot\|_{\mathrm{F}}$ denote the $\ell_{2}$ norm and Frobenius norm, respectively. $\Re\{\cdot\}$ and $\Im\{\cdot\}$ denote the real and imaginary parts, respectively. $\mathrm{j}$ is the imaginary unit. $\bm{I}_{N}$ denotes the $N\times N$ identity matrix. $\mathcal{CN}(\mu,\sigma^{2})$ denotes a complex Gaussian distribution.

wbsphack @@writeaux“newlabelS2.1wcurrentlabel1wesphack

We consider a time-division duplex (TDD) multiuser MIMO-OFDM downlink system implemented with the SEMRA architecture illustrated in Fig. LABEL:fig:system_model. The base station (BS) employs an $M=M_{y}M_{z}$ uniform planar array (UPA) of SEMRA elements to serve $U$ single-antenna user equipments (UEs), each with a fixed position and radiation pattern. The inter-element spacings of the UPA along the y- and z-axes are $d_{y}$ and $d_{z}$, respectively. The number of OFDM subcarriers is $G$ with subcarrier spacing $\Delta f\triangleq B_{w}/G$, and $f_{g}\triangleq(g-1)\Delta f$ ($1\leq g\leq G$) denotes the frequency of the $g$-th subcarrier. The received signal at the $u$-th UE ($1\leq u\leq U$) and the $g$-th subcarrier ($1\leq g\leq G$) is given by

$$wbsphack@@writeaux{}{\newlabel{equ.rx}{{wcurrentlabel}{1}}}wesphacky_{u,g}=\bm{h}_{u,g}^{H}\bm{W}_{g}\bm{s}_{g}+n_{u,g},$$

(1)

where $\bm{h}_{u,g}=[h_{u,1,g},\ldots,h_{u,M,g}]^{T}\in\mathbb{C}^{M}$ is the downlink channel vector between the BS and the $u$-th UE, $\bm{W}_{g}=[\bm{w}_{1,g},\ldots,\bm{w}_{U,g}]\in\mathbb{C}^{M\times U}$ is the digital precoding matrix with $\bm{w}_{u,g}\in\mathbb{C}^{M\times 1}$ denoting the digital precoder vector for the $u$-th UE on the $g$-th subcarrier, $\bm{s}_{g}\in\mathbb{C}^{U}$ is the data symbol vector satisfying $\mathbb{E}\{\bm{s}_{g}\bm{s}_{g}^{H}\}=\bm{I}_{U}$, and $n_{u,g}\sim\mathcal{CN}(0,\sigma_{n}^{2})$ is additive white Gaussian noise (AWGN).

wbsphack @@writeaux“newlabelS2.2wcurrentlabel1wesphack The channel between the $m$-th BS antenna and the $u$-th UE on the $g$-th subcarrier is modeled as [38.901]

$\displaystyle\begin{split}wbsphack@@writeaux{}{\newlabel{equ.ch}{{wcurrentlabel}{1}}}wesphackh_{u,m,g}&=\sum_{i=1}^{L_{u}}\tilde{x}_{i,u}f_{\mathrm{rx},u}(\vartheta_{i,u},\varphi_{i,u})f_{\mathrm{tx},m}(\theta_{i,u},\phi_{i,u})\\ &\times e^{-\mathrm{j}\frac{2\pi}{\lambda}(\bm{k}_{\mathrm{tx},i,u}^{T}\bm{p}_{m}+\bm{k}_{\mathrm{rx},i,u}^{T}\bm{r}_{u})}\cdot e^{-\mathrm{j}2\pi\tau_{i,u}f_{g}},\end{split}$

(2)

where $L_{u}$ is the number of channel paths, $\tilde{x}_{i,u}$ and $\tau_{i,u}$ denote the complex gain and delay of the $i$-th path, respectively, $(\theta_{i,u},\phi_{i,u})$ and $(\vartheta_{i,u},\varphi_{i,u})$ are respectively the AoD and AoA, $f_{\mathrm{tx},m}(\cdot)$ and $f_{\mathrm{rx},u}(\cdot)$ denote the real-valued scalar radiation-pattern gains of the BS and UE antennas, respectively, and $\bm{p}_{m}$, $\bm{r}_{u}$ are the position vectors of BS antenna $m$ and UE $u$, respectively. Throughout this paper, the BS-side departure azimuth is restricted to the front half-space, i.e., $\phi_{i,u}\in[-\pi/2,\pi/2]$. The unit direction vectors $\bm{k}_{\mathrm{tx},i,u},\bm{k}_{\mathrm{rx},i,u}\in\mathbb{R}^{3}$ are defined as $\bm{k}_{\mathrm{tx},i,u}=[\sin\theta_{i,u}\cos\phi_{i,u},\sin\theta_{i,u}\sin\phi_{i,u},\cos\theta_{i,u}]^{T}$ and $\bm{k}_{\mathrm{rx},i,u}=[\sin\vartheta_{i,u}\cos\varphi_{i,u},\sin\vartheta_{i,u}\sin\varphi_{i,u},\cos\vartheta_{i,u}]^{T}$.

Incorporating the delay into the path gain as $x_{i,u,g}=\tilde{x}_{i,u}e^{-\mathrm{j}2\pi\tau_{i,u}f_{g}}$, the channel can be written in the compact matrix form as

$$wbsphack@@writeaux{}{\newlabel{equ.ch_{c}ompact}{{wcurrentlabel}{1}}}wesphackh_{u,m,g}=\bm{f}_{\mathrm{rx},u}^{T}\bm{A}_{u}\bm{\Sigma}_{u,g}\bm{B}_{u,m}\bm{f}_{\mathrm{tx},u,m},$$

(3)

where $\bm{f}_{\mathrm{rx},u}\in\mathbb{R}^{L_{u}}$ and $\bm{f}_{\mathrm{tx},u,m}\in\mathbb{R}^{L_{u}}$ collect the receive and transmit pattern gains along all $L_{u}$ paths, with the $i$-th elements being $f_{\mathrm{rx},u}(\vartheta_{i,u},\varphi_{i,u})$ and $f_{\mathrm{tx},m}(\theta_{i,u},\phi_{i,u})$, respectively, $\bm{A}_{u}=\mathrm{diag}(\bm{a}_{u})\in\mathbb{C}^{L_{u}\times L_{u}}$ and $\bm{B}_{u,m}=\mathrm{diag}(\bm{b}_{u,m})\in\mathbb{C}^{L_{u}\times L_{u}}$ capture the array steering phases, with $[\bm{a}_{u}]_{i}=e^{-\mathrm{j}\frac{2\pi}{\lambda}\bm{k}_{\mathrm{rx},i,u}^{T}\bm{r}_{u}}$ and $[\bm{b}_{u,m}]_{i}=e^{-\mathrm{j}\frac{2\pi}{\lambda}\bm{k}_{\mathrm{tx},i,u}^{T}\bm{p}_{m}}$, and $\bm{\Sigma}_{u,g}=\mathrm{diag}([x_{1,u,g},\ldots,x_{L_{u},u,g}])\in\mathbb{C}^{L_{u}\times L_{u}}$ collects the per-path gains.

wbsphack @@writeaux“newlabelS2.3wcurrentlabel1wesphack To enable EM-domain radiation pattern optimization, we represent the transmit pattern as a weighted combination of $K$ orthonormal basis functions

$$wbsphack@@writeaux{}{\newlabel{equ.sh}{{wcurrentlabel}{1}}}wesphackf_{\mathrm{tx},m}(\theta,\phi)=\sum_{k=1}^{K}\alpha_{k,m}\omega_{k}(\theta,\phi),$$

(4)

where $\{\omega_{k}(\theta,\phi)\}_{k=1}^{K}$ are real-valued orthonormal basis functions over the unit sphere, constructed from real spherical harmonics as detailed in [Ying2025TCOM]. The orthonormality condition reads $\int_{0}^{2\pi}\!\int_{0}^{\pi}\omega_{k}(\theta,\phi)\omega_{k^{\prime}}(\theta,\phi)\sin\theta\,\mathrm{d}\theta\,\mathrm{d}\phi=\delta_{k,k^{\prime}}$, and $\bm{\alpha}_{m}=[\alpha_{1,m},\ldots,\alpha_{K,m}]^{T}\in\mathbb{R}^{K}$ is the weight coefficient vector for the $m$-th antenna. Under the adopted decomposition, $f_{\mathrm{tx},m}(\theta,\phi)$ models the real-valued directional gain, and we further assume that $\|\bm{\alpha}_{m}\|_{2}=1$. Applying (LABEL:equ.sh) to the transmit pattern gain vector yields

$$wbsphack@@writeaux{}{\newlabel{equ.ftx}{{wcurrentlabel}{1}}}wesphack\bm{f}_{\mathrm{tx},u,m}=\bm{\Omega}_{u}\bm{\alpha}_{m}\in{\mathbb{R}^{L_{u}}},\quad\forall u,m,$$

(5)

where $\bm{\alpha}_{m}$ parameterizes real gain-pattern coefficients, and the $(i,k)$-th element of $\bm{\Omega}_{u}\in\mathbb{R}^{L_{u}\times K}$ is $[\bm{\Omega}_{u}]_{i,k}=\omega_{k}(\theta_{i,u},\phi_{i,u})$. Substituting (LABEL:equ.ftx) into (LABEL:equ.ch_compact), the channel can be expressed as

$$wbsphack@@writeaux{}{\newlabel{equ.ecsi}{{wcurrentlabel}{1}}}wesphackh_{u,m,g}=\bm{f}_{\mathrm{rx},u}^{T}\bm{A}_{u}\bm{\Sigma}_{u,g}\bm{B}_{u,m}\bm{\Omega}_{u}\bm{\alpha}_{m}=\bm{q}_{u,m,g}^{H}\bm{\alpha}_{m}.$$

(6)

We define the eCSI vector as

$$wbsphack@@writeaux{}{\newlabel{equ.q_{d}ef}{{wcurrentlabel}{1}}}wesphack\bm{q}_{u,m,g}^{H}\triangleq\bm{f}_{\mathrm{rx},u}^{T}\bm{A}_{u}\bm{\Sigma}_{u,g}\bm{B}_{u,m}\bm{\Omega}_{u}\in\mathbb{C}^{1\times K},$$

(7)

and define the corresponding position-independent equivalent channel gain vector as

$$wbsphack@@writeaux{}{\newlabel{equ.v_{d}ef}{{wcurrentlabel}{1}}}wesphack\bm{\chi}_{u,g}\triangleq\left(\bm{f}_{\mathrm{rx},u}^{T}\bm{A}_{u}\bm{\Sigma}_{u,g}\right)^{H}\in\mathbb{C}^{L_{u}}.$$

(8)

It collects the position-independent UE-side and path-gain terms. The corresponding column-vector form of the eCSI is $\bm{q}_{u,m,g}=\bm{\Omega}_{u}^{T}\bm{B}_{u,m}^{H}\bm{\chi}_{u,g}\in\mathbb{C}^{K\times 1}.$

By aggregating the eCSI across all $M$ antennas, we define the stacked eCSI vector $\bm{q}_{u,g}\triangleq[\bm{q}_{u,1,g}^{T},\ldots,\bm{q}_{u,M,g}^{T}]^{T}\in\mathbb{C}^{MK\times 1}$ and its conjugated form $\bar{\bm{q}}_{u,g}\triangleq\bm{q}_{u,g}^{*}$. We also define the block-diagonal EM precoding matrix $\bm{\Lambda}\triangleq\mathrm{Blkdiag}\{\bm{\alpha}_{1},\ldots,\bm{\alpha}_{M}\}\in\mathbb{R}^{MK\times M}$. The effective channel then satisfies

$$wbsphack@@writeaux{}{\newlabel{equ.h_{d}ef}{{wcurrentlabel}{1}}}wesphack\bm{h}_{u,g}^{H}=\bar{\bm{q}}_{u,g}^{H}\bm{\Lambda}.$$

(9)

Substituting (LABEL:equ.h_def) into (LABEL:equ.rx) yields $y_{u,g}=\bar{\bm{q}}_{u,g}^{H}\bm{\Lambda}\bm{W}_{g}\bm{s}_{g}+n_{u,g}$, so the signal-to-interference-plus-noise ratio (SINR) for the $u$-th UE on the $g$-th subcarrier is

$$wbsphack@@writeaux{}{\newlabel{equ.sinr}{{wcurrentlabel}{1}}}wesphack\mathrm{SINR}_{u,g}=\frac{|\bar{\bm{q}}_{u,g}^{H}\bm{\Lambda}\bm{w}_{u,g}|^{2}}{\sum_{\begin{subarray}{c}\ell=1\\ \ell\neq u\end{subarray}}^{U}|\bar{\bm{q}}_{u,g}^{H}\bm{\Lambda}\bm{w}_{\ell,g}|^{2}+\sigma_{n}^{2}}.$$

(10)

Accordingly, the sum SE is given by

$$wbsphack@@writeaux{}{\newlabel{equ.se}{{wcurrentlabel}{1}}}wesphackR=\sum_{g=1}^{G}\sum_{u=1}^{U}\log_{2}(1+\mathrm{SINR}_{u,g}).$$

(11)

wbsphack @@writeaux“newlabelS2.4wcurrentlabel1wesphack Unlike the EMRA with TFA architecture, the proposed SEMRA architecture jointly optimizes the per-antenna EM-domain pattern weights $\bm{\alpha}_{m}$ and antenna positions $\bm{p}_{m}$. Let $\bar{\bm{p}}_{m}$ denote the reference position of the $m$-th antenna in the reference UPA. In SEMRA, each transmit antenna is permitted to move within a local feasible region $\mathcal{S}_{m}$ centered at $\bar{\bm{p}}_{m}$ (depicted as the local movement region in Fig. LABEL:fig:system_model), which is defined as

$$wbsphack@@writeaux{}{\newlabel{equ.constraint}{{wcurrentlabel}{1}}}wesphack\bm{p}_{m}\in\mathcal{S}_{m}\triangleq\left\{\bm{x}\in\mathbb{R}^{3}\mid\|\bm{x}-\bar{\bm{p}}_{m}\|_{\infty}\leq D_{\max},\ [\bm{x}]_{1}=0\right\},$$

(12)

where $D_{\max}$ is the maximum allowable displacement for each antenna. To avoid antenna collisions, we require a minimum inter-antenna separation $D_{\mathrm{sep}}>0$ over all feasible positions, which is ensured for the reference UPA if $d_{y}\geq 2D_{\max}+D_{\mathrm{sep}}$ and $d_{z}\geq 2D_{\max}+D_{\mathrm{sep}}$. The position $\bm{p}_{m}$ enters the channel through the BS-side steering matrix $\bm{B}_{u,m}$ in (LABEL:equ.ch_compact), introducing an additional spatial degree of freedom for performance optimization.

According to (LABEL:equ.h_def), the effective channels are jointly determined by the antenna positions $\bm{p}\triangleq[\bm{p}_{1}^{T},\ldots,\bm{p}_{M}^{T}]^{T}\in\mathbb{R}^{3M}$ and the per-antenna EM-domain weights $\bm{\alpha}_{m}$, whose effects are captured by $\bar{\bm{q}}_{u,g}$ and $\bm{\Lambda}$, respectively. Let $\bm{W}\triangleq\{\bm{W}_{g}\}_{g=1}^{G}$ denote the collection of digital precoders across all subcarriers. We aim to maximize the sum SE of all UEs in (LABEL:equ.se) by jointly optimizing digital-domain precoder $\bm{W}$, EM-domain radiation-pattern weights $\bm{\alpha}_{m}$, and spatial-domain antenna position $\bm{p}$, leading to

$\displaystyle wbsphack@@writeaux{}{\newlabel{equ.problem}{{wcurrentlabel}{1}}}wesphack\begin{aligned} \mathcal{P}_{\mathrm{SE}}:\quad\max_{\{\bm{\alpha}_{m}\},\bm{W},\bm{p}}\quad&R\\ \text{s.t.}\quad&\sum_{g=1}^{G}\|\bm{W}_{g}\|_{\mathrm{F}}^{2}\leq P_{T},\\ &\|\bm{\alpha}_{m}\|_{2}=1,\ \bm{p}_{m}\in\mathcal{S}_{m},\ \forall m.\end{aligned}$

(13)

where $R$ is defined in (LABEL:equ.se) and depends on $\bm{W}$, $\bm{\alpha}_{m}$, and $\bm{p}$ through the effective channel vectors $\bm{h}_{u,g}$ in (LABEL:equ.h_def). The resulting problem is highly nonconvex due to the coupled variables and the unit-norm constraints, making it challenging to solve optimally.

wbsphack @@writeaux“newlabelS3.1zfwcurrentlabel1wesphack

Given the current effective channels $\bm{h}_{u,g}$ in (LABEL:equ.h_def), we adopt the ZF precoder. Define the normalized channel direction and the stacked matrix as in (LABEL:equ.O_g_def).

$$wbsphack@@writeaux{}{\newlabel{equ.O_{g}{}_{d}ef}{{wcurrentlabel}{1}}}wesphack\tilde{\bm{h}}_{u,g}\triangleq\frac{\bm{h}_{u,g}}{\|\bm{h}_{u,g}\|_{2}},\qquad\tilde{\bm{H}}_{g}\triangleq[\tilde{\bm{h}}_{1,g},\ldots,\tilde{\bm{h}}_{U,g}]\in\mathbb{C}^{M\times U}.$$

(14)

For $M\geq U$, the normalized-channel ZF precoder is computed as

$$wbsphack@@writeaux{}{\newlabel{equ.zf_{n}ormalize}{{wcurrentlabel}{1}}}wesphack\bm{W}_{g}=\sqrt{\frac{P_{T}}{G}}\frac{\tilde{\bm{W}}_{g}}{\|\tilde{\bm{W}}_{g}\|_{\mathrm{F}}},$$

(15)

where $\tilde{\bm{W}}_{g}\triangleq\tilde{\bm{H}}_{g}\big(\tilde{\bm{H}}_{g}^{H}\tilde{\bm{H}}_{g}\big)^{-1}\in\mathbb{C}^{M\times U}$ is the unnormalized ZF precoder built from the normalized channel directions in (LABEL:equ.O_g_def).

With $\bm{p}$ and the digital precoders fixed, we update the EM-domain precoder by optimizing the block-diagonal matrix $\bm{\Lambda}=\mathrm{Blkdiag}\{\bm{\alpha}_{1},\ldots,\bm{\alpha}_{M}\}$ under per-antenna unit-norm constraints. This is formulated on a masked oblique manifold. We define the oblique manifold

$$wbsphack@@writeaux{}{\newlabel{equ.oblique_{d}ef}{{wcurrentlabel}{1}}}wesphack\mathcal{OB}\triangleq\left\{\bm{\Lambda}\in\mathbb{R}^{KM\times M}:\left[\bm{\Lambda}^{T}\bm{\Lambda}\right]_{m,m}=1,\;\forall m\right\}.$$

(16)

Let $\bm{M}_{0}\triangleq\mathrm{Blkdiag}\{\bm{1}_{K},\ldots,\bm{1}_{K}\}\in\mathbb{R}^{KM\times M}$ and impose the structural constraint $\bm{\Lambda}=\bm{\Lambda}\odot\bm{M}_{0}$. Equivalently, the feasible EM set is $\{\bm{\Lambda}\in\mathcal{OB}:\bm{\Lambda}=\bm{\Lambda}\odot\bm{M}_{0}\}$. We define the cost function $\mathfrak{f}(\bm{\Lambda})\triangleq-R$. The Riemannian gradient at $\bm{\Lambda}$ is the orthogonal projection of the Euclidean gradient $\nabla\mathfrak{f}(\bm{\Lambda})$ onto the tangent space, given by

$$wbsphack@@writeaux{}{\newlabel{equ.oblique_{g}rad}{{wcurrentlabel}{1}}}wesphack\mathrm{grad}\mathfrak{f}(\bm{\Lambda})=\nabla\mathfrak{f}(\bm{\Lambda})-\bm{\Lambda}\,\mathrm{ddiag}\!\left(\bm{\Lambda}^{T}\nabla\mathfrak{f}(\bm{\Lambda})\right),$$

(17)

where $\mathrm{ddiag}(\cdot)$ sets all off-diagonal entries of its matrix argument to zero. The retraction is column-wise normalization $\mathrm{Retr}_{\bm{\Lambda}}(\gamma\bm{D})\triangleq\mathrm{normalize}(\bm{\Lambda}+\gamma\bm{D})$.

We adopt the Riemannian steepest-descent direction on this masked-oblique set, namely

$$wbsphack@@writeaux{}{\newlabel{equ.oblique_{c}g_{d}ir}{{wcurrentlabel}{1}}}wesphack\bm{D}^{(n)}=-\mathrm{grad}\mathfrak{f}\!\left(\bm{\Lambda}^{(n)}\right),$$

(18)

and choose a suitable step size $\gamma^{(n)}$ for the current EM subproblem with the digital precoders fixed. The EM precoder update is then

$$wbsphack@@writeaux{}{\newlabel{equ.lambda_{u}pdate_{z}f}{{wcurrentlabel}{1}}}wesphack\begin{split}\bm{\Lambda}^{(n+1)}&=\mathrm{Retr}_{\bm{\Lambda}^{(n)}}\!\left(\gamma^{(n)}\bm{D}^{(n)}\right)\\ &=\mathrm{normalize}\!\left(\bm{\Lambda}^{(n)}+\gamma^{(n)}\bm{D}^{(n)}\right).\end{split}$$

(19)

To keep the per-outer-iteration cost moderate, the EM block in both tri-domain algorithms below uses a single manifold-gradient step in each outer iteration rather than an inner loop to full convergence. Since the digital and spatial blocks are refreshed between successive EM updates, the EM descent direction is recomputed from the current gradient in every outer iteration. For compactness, variables without outer-iteration superscripts in the algorithm below denote the latest available block values within the current outer iteration. To express the Euclidean gradient compactly, we use (LABEL:equ.h_def), i.e., $\bm{h}_{u,g}^{H}=\bar{\bm{q}}_{u,g}^{H}\bm{\Lambda}$. Let $\zeta\triangleq 1/\sigma_{n}^{2}$, and $\bm{W}_{\bar{u},g}$ denote the digital precoder matrix obtained from $\bm{W}_{g}$ by removing its $u$-th column. The Euclidean gradient $\nabla\mathfrak{f}(\bm{\Lambda})$ is given by [Ying2025TCOM]

$$wbsphack@@writeaux{}{\newlabel{equ.grad_{e}uc_{z}f_{l}ambda}{{wcurrentlabel}{1}}}wesphack\nabla\mathfrak{f}\!\left(\bm{\Lambda}\right)=\left(-\sum_{g=1}^{G}\sum_{u=1}^{U}\frac{\bm{\Gamma}_{u,g}^{(1)}-\bm{\Gamma}_{u,g}^{(2)}}{\ln 2}\right)\odot\bm{M}_{0}.$$

(20)

The Hadamard masking enforces the block-diagonal sparsity of $\bm{\Lambda}$, such that each column update only acts on its corresponding $K$-dimensional pattern-weight block. Starting from a feasible block-diagonal initialization, the masked gradient in (LABEL:equ.grad_euc_zf_lambda) and the retraction in (LABEL:equ.lambda_update_zf) preserve this structure throughout the ZF-EM updates. Moreover,

	$\displaystyle wbsphack@@writeaux{}{\newlabel{equ.Gamma_{z}f_{d}ef}{{wcurrentlabel}{1}}}wesphack\bm{\Gamma}_{u,g}^{(1)}$	$\displaystyle\triangleq\frac{2\zeta\Re\!\left\{\bar{\bm{q}}_{u,g}\bar{\bm{q}}_{u,g}^{H}\bm{\Lambda}\bm{W}_{g}\bm{W}_{g}^{H}\right\}}{1+\zeta\|\bar{\bm{q}}_{u,g}^{H}\bm{\Lambda}\bm{W}_{g}\|_{2}^{2}},$		(21)
	$\displaystyle\bm{\Gamma}_{u,g}^{(2)}$	$\displaystyle\triangleq\frac{2\zeta\Re\!\left\{\bar{\bm{q}}_{u,g}\bar{\bm{q}}_{u,g}^{H}\bm{\Lambda}\bm{W}_{\bar{u},g}\bm{W}_{\bar{u},g}^{H}\right\}}{1+\zeta\|\bar{\bm{q}}_{u,g}^{H}\bm{\Lambda}\bm{W}_{\bar{u},g}\|_{2}^{2}}.$		(22)

wbsphack @@writeaux“newlabelS3.2zfwcurrentlabel1wesphack

Using the available parametric eCSI model, we derive the spatial gradient of the SE objective $R$ with respect to the conjugate channel coefficient $h_{u,m,g}^{*}$.