This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Polarization Aware Movable Antenna

Runxin Zhang, Yulin Shao, Yonina C. Eldar R. Zhang and Y. Shao are with the State Key Laboratory of Internet of Things for Smart City and the Department of Electrical and Computer Engineering, University of Macau, Macau S.A.R., 999078, China. Y. Shao is also with the Department of Electrical and Electronic Engineering, Imperial College London, London SW7 2AZ, U.K. Y. C. Eldar is with the Department of Mathematics and Computer Science, Weizmann Institute of Science, Rehovot 7610001, Israel.Correspondence to: ylshao@um.edu.mo.
Abstract

This paper presents a polarization-aware movable antenna (PAMA) framework that integrates polarization effects into the design and optimization of movable antennas (MAs). While MAs have proven effective at boosting wireless communication performance, existing studies primarily focus on phase variations caused by different propagation paths and leverage antenna movements to maximize channel gains. This narrow focus limits the full potential of MAs. In this work, we introduce a polarization-aware channel model rooted in electromagnetic theory, unveiling a defining advantage of MAs over other wireless technologies such as precoding: the ability to optimize polarization matching. This new understanding enables PAMA to extend the applicability of MAs beyond radio-frequency, multipath-rich scenarios to higher-frequency bands, such as mmWave, even with a single line-of-sight (LOS) path. Our findings demonstrate that incorporating polarization considerations into MAs significantly enhances efficiency, link reliability, and data throughput, paving the way for more robust and efficient future wireless networks.

Index Terms:
Movable antenna, fluid antenna system, polarization, mmWave, MU-MISO.

I Introduction

The relentless pursuit of high-speed and reliable wireless communication has become a cornerstone of modern technological advancement. Emerging applications such as autonomous vehicles, the artificial intelligence of things, and augmented reality demand seamless connectivity with unprecedented data rates and minimal latency[1, 2]. These evolving requirements are pushing the boundaries of existing wireless communication systems, necessitating innovative solutions to address challenges like spectrum scarcity, signal degradation, and dynamic environmental conditions.

One promising technique that has emerged to meet these challenges is the use of movable antennas (MAs) [3, 4, 5, 6], also known as fluid antenna systems (FAS) [7, 8, 9, 10, 11]. By dynamically adjusting the position and orientation of antennas, it is possible to achieve enhanced signal-to-noise ratio (SNR), wider coverage, and adaptability to changing environments and user locations, without solely relying on traditional methods like increasing transmission power, beamforming, or deploying additional infrastructure[12, 13, 14, 15, 16, 17].

Refer to caption
Figure 1: Movable antennas with the capability for 3D translation and rotation.

Previous research has explored the applications of MAs and FAS in radio frequency bands characterized by rich multipath environments [3, 7]. In these settings, antenna movement is leveraged to exploit spatial diversity arising from amplitude fluctuations caused by multipath propagation.

In [7], the authors propose FAS that allows the receiving antenna to move to positions to capture the strongest signal. Analysis of outage probability demonstrates that, within a movement range of 0.2λ\lambda, FAS can achieve performance comparable to a maximal ratio combining system with five antennas spanning 2λ\lambda. This flexibility in antenna positioning offers significant advantages, especially in scenarios requiring massive connectivity. It allows each user to find an optimal position where instantaneous interference is minimized, effectively creating deep nulls. This spatial multiplexing gain facilitates multi-user access, a concept introduced as Fluid Antenna Multiple Access (FAMA) in [8]. Building upon this concept, [9] extends FAS to a two-dimensional (2D) implementation and discusses reconfigurable pixel technologies as a practical approach for realization. Additionally, [18] introduces a continuous FAS (CFAS) architecture, demonstrating that CFAS outperforms discrete FAS with limited ports, particularly at medium to high threshold levels.

The MA system presented in [3] exhibits additional phase variations due to antenna translation, significantly surpassing fixed-position antennas in both deterministic and stochastic channels. The system’s performance further improves with increased positioning accuracy, larger movement ranges, and enhanced multipath effects. Building on 2D MA system, advancements in six-dimensional (6D) systems, such as those discussed in [4], enable simultaneous three-dimensional (3D) translation and rotation, enhancing the system’s directional capabilities. By imposing spatial constraints on surfaces equipped with 6D movable antennas (6DMA), user interference is reduced, improving the signal-to-interference-plus-noise ratio (SINR). In [5], the 6DMA model is simplified by integrating conventional fixed-position antenna arrays with 6DMA surfaces capable of rotating along a circular path. By limiting the degrees of freedom for surface movement to the azimuth angle along the circular track, computational complexity is significantly reduced while enhancing the system’s maximum capacity.

Despite these advancements, the impact of polarization changes due to antenna movement and rotation has been largely overlooked. Polarization is a fundamental characteristic of electromagnetic waves, describing the orientation of the electric field vector as the wave propagates through space [19, 20, 21, 22, 23]. In wireless communication systems, proper alignment of polarization between transmitting and receiving antennas is crucial for efficient signal transmission and reception. When polarizations are properly aligned, the receiving antenna can maximally capture the transmitted energy, resulting in stronger signal strength, improved SNR, and enhanced overall system performance. Conversely, polarization mismatch occurs when there is a misalignment between the polarization orientations of the transmitting and receiving antennas, leading to significant signal attenuation as the misaligned component of the electromagnetic wave is effectively lost. Neglecting polarization effects can lead to incomplete system models and suboptimal performance assessments, hindering a full understanding of the impact and advantages of MAs.

This oversight is particularly critical as communication systems move towards higher frequency bands, such as the millimeter-wave (mmWave) spectrum [24, 25]. At these higher frequencies, signal propagation becomes more directional with reduced diffraction around obstacles, often resulting in a single line-of-sight (LOS) link between transmitter and receiver [26]. Previously, it was believed that MAs offer no gains in such LOS scenarios due to the lack of multipath diversity [7, 3]. However, our formulation incorporating polarization effects indicates otherwise – MAs can achieve significant performance gains even in LOS scenarios by facilitating polarization matching. As such, there is a clear gap in current research: the need to incorporate polarization effects into the design and optimization of MA systems, especially as the wireless industry advances towards higher frequency bands to achieve greater bandwidths and data rates.

To fill the gap, this paper makes the following main contributions:

  • We introduce a polarization-aware movable antenna (PAMA) framework, which incorporates polarization effects into the design and optimization of antenna movement and rotation. PAMA uncovers a previously overlooked, but defining benefit of antenna movement – facilitating optimal polarization matching between the transmitter and receiver. By accounting for the impact of polarization changes due to movement and rotation, PAMA provides a more refined system model and deepens our understanding of the dynamic capabilities of MAs.

  • Building upon the electromagnetic theory, we develop a polarization-aware channel model that describes the impact of 3D movement and rotation of antennas on the channel gains. Grounded in Maxwell’s equations, this model allows for precise characterization of how antenna positioning and orientation affect polarization alignment and, consequently, signal reception quality. By integrating these factors into the system design, we enable the optimization of antenna movements to maintain polarization alignment, thereby maximizing the sum rate of the communication system.

  • Contrary to prior beliefs that MAs are more effective in low-frequency, multipath-rich environments, our PAMA framework demonstrates that MAs are highly effective in higher frequency bands such as the mmWave band, even with a single LOS link. The smaller size of antennas at higher frequencies makes them easier to move and rotate, facilitating practical implementation of PAMA systems in these bands. This new understanding extends the potential and application of MA to next-generation wireless communication systems that demand higher bandwidths.

Organization: The remainder of this paper is structured as follows. Section II introduces the system model and essential geometric concepts used throughout the work. Section III presents our PAMA framework and the polarization-aware channel model. Section IV explores the impacts of 3D antenna translation and rotation on the PAMA channel gains. Section V demonstrates the potential of PAMA through a sum-rate optimization problem. Finally, Section VI concludes the paper.

Notations: We use boldface lowercase letters to denote column vectors (e.g., 𝒉\bm{h}, 𝒘\bm{w}) and boldface uppercase letters to denote matrices (e.g., 𝑯\bm{H}, 𝑷\bm{P}). For a vector or matrix, ()(\cdot)^{\top} denotes the transpose, ()(\cdot)^{*} denotes the complex conjugate, and ()H(\cdot)^{H} denotes the conjugate transpose. \mathbb{C} stand for the sets of complex values. The imaginary unit is represented by jj. The Euclidean norm of a vector 𝒑\bm{p} is denoted by |𝒑||\bm{p}|. For two vectors 𝒂\bm{a} and 𝒃\bm{b}, 𝒂,𝒃\langle\bm{a},\bm{b}\rangle denotes the angle between them, and 𝒂𝒃\bm{a}\cdot\bm{b} denotes their inner product.

II System Model

Refer to caption
Figure 2: The two coordinate systems and their unit vectors.

We consider a mobile edge network within a 3D environment. The network comprises a base station (BS) and KK mobile users distributed within the 3D space. The BS is equipped with LL half-wave dipole antennas, and each mobile user is equipped with a single half-wave dipole antenna. We focus on the downlink transmission from the BS to the users, operating in the mmWave frequency band. Both the antennas at the BS and mobile users are movable.

Definition 1 (Movable Antennas).

A movable antenna is capable of performing 3D translations and 3D rotations, allowing it to navigate through and orient within 3D space, as depicted in Fig. 1.

Translation typically requires larger operational spaces, often spanning tens of wavelengths, whereas rotation can be accommodated within a more compact area, approximately the size of a half-wavelength sphere, as shown in Fig. 1. In this work, we assume that the antennas on mobile users are restricted to 3D rotation, due to the spatial constraint. In contrast, the LL antennas at the BS can perform both 3D rotation and translation independently within a designated region. The extension of this model to include translatable antennas for mobile users is straightforward.

In this work, we consider the far-field scenario, which implies the following:

  • Plane wave. The waves received by the users can be treated as plane waves because the propagation distance is much greater than the wavelength and the dimensions of the movable regions.

  • Point representation of antennas. The physical sizes of the transmitting and receiving antennas are negligible compared to the propagation distance, allowing us to represent their positions as points specified by their center coordinates.

To mathematically represent the network and the positions of the antennas, we establish two coordinate systems: the Cartesian Coordinate System (CCS) and the Spherical Coordinate System (SCS), as depicted in Fig. 2.

Definition 2 (Coordinate Systems).

We define a CCS to describe the positions of the transmitting and receiving antennas, with three orthogonal unit vectors denoted as 𝐢\bm{i}, 𝐣\bm{j}, and 𝐤\bm{k}.

Correspondingly, an SCS is established with its origin at point 𝐨=(0,0,0)\bm{o}=(0,0,0) in the CCS, and the direction of 𝐤\bm{k} serving as the reference axis. The SCS has three unit vectors: the radial unit vector 𝛎\bm{\nu}, the polar angle unit vector ϑ\bm{\vartheta}, and the azimuthal angle unit vector 𝛗\bm{\varphi}.

TABLE I: Conversion formulas between the unit vectors of CCS and SCS.
Dot products 𝒊\bm{i} 𝒋\bm{j} 𝒌\bm{k}
𝝂\bm{\nu} sinθcosϕ\sin{\theta}\cos{\phi} sinθsinϕ\sin{\theta}\sin{\phi} cosθ\cos{\theta}
ϑ\bm{\vartheta} cosθcosϕ\cos{\theta}\cos{\phi} cosθsinϕ\cos{\theta}\sin{\phi} sinθ-\sin{\theta}
𝝋\bm{\varphi} sinθ-\sin{\theta} cosϕ\cos{\phi} 0

Table I provides the conversion formulas between the unit vectors of the two coordinate systems. For instance, the radial unit vector 𝝂\bm{\nu} in the SCS can be expressed in terms of the CCS unit vectors as 𝝂=sin(θ)cos(ϕ)𝒊+sin(θ)sin(ϕ)𝒋+cos(θ)𝒌\bm{\nu}=\sin(\theta)\cos(\phi)\bm{i}+\sin(\theta)\sin(\phi)\bm{j}+\cos(\theta)\bm{k}.

Given the two coordinate systems, we introduce several definitions regarding the antennas’ positions and orientations.

Definition 3 (Positions of the Antennas).

In the far-field scenario, the position of an antenna can be represented by its center coordinates. We denote the center coordinate of the \ell-th transmitting antennas at the BS by 𝐩t,=(xt,,yt,,zt,)\bm{p}_{t,\ell}=\left(x_{t,\ell},y_{t,\ell},z_{t,\ell}\right), =1,2,,L\ell=1,2,\ldots,L. Similarly, the center coordinate of the receiving antenna at the κ\kappa-th user is denoted by 𝐩r,κ=(xr,κ,yr,κ,zr,κ)\bm{p}_{r,\kappa}=\left(x_{r,\kappa},y_{r,\kappa},z_{r,\kappa}\right) for κ=1,2,,K\kappa=1,2,\ldots,K.

Definition 4 (Directions of Antennas).

The direction of an antenna, characterized by a unit vector 𝐧\bm{n}, is the orientation in which the antenna is pointed. We denote the direction for the \ell-th transmitting antenna at the BS by 𝐧t,=(xt,(n),yt,(n),zt,(n))\bm{n}_{t,\ell}=\big{(}x_{t,\ell}^{(n)},y_{t,\ell}^{(n)},z_{t,\ell}^{(n)}\big{)}, and for the receiving antenna at the κ\kappa-th user by 𝐧r,κ=(xr,κ(n),yr,κ(n),zr,κ(n))\bm{n}_{r,\kappa}=\big{(}x_{r,\kappa}^{(n)},y_{r,\kappa}^{(n)},z_{r,\kappa}^{(n)}\big{)}.

Definition 5 (Polar and Azimuthal Angles of Antennas).

The polar angle of an antenna, denoted by θ(n)[0,π]\theta^{(n)}\in[0,\pi], is the angle between 𝐤\bm{k} and the direction of the antenna. The azimuthal angle of an antenna, denoted by ϕ(n)[0,2π)\phi^{(n)}\in[0,2\pi), is the angle measured counterclockwise from 𝐢\bm{i} to the projection of the antenna onto the x𝐨yx\bm{o}y plane. For the \ell-th transmitting antenna, these angles are denoted by θt,(n)\theta_{t,\ell}^{(n)} and ϕt,(n)\phi_{t,\ell}^{(n)}, respectively. For the receiving antenna at the κ\kappa-th user, they are denoted by θr,κ(n)\theta_{r,\kappa}^{(n)} and ϕr,κ(n)\phi_{r,\kappa}^{(n)}, respectively.

Refer to caption
(a) \ell-th transmitting antenna
Refer to caption
(b) κ\kappa-th receiving antenna
Figure 3: Illustration of the polar and azimuthal angles, along with the geometric representation of the propagation link between the \ell-th transmitting antenna and the κ\kappa-th receiving antenna.

The polar and azimuthal angles are illustrated in Fig. 3. A summary of the definitions and their geometric relationships is given in Tab. II.

Lemma 1.

The direction of an antenna in the CCS can be represented by its polar and azimuthal angles. For the \ell-th transmitting antenna, we have

𝒏t,=(xt,(n),yt,(n),zt,(n))\displaystyle\bm{n}_{t,\ell}=\left(x_{t,\ell}^{(n)},y_{t,\ell}^{(n)},z_{t,\ell}^{(n)}\right)
=(sinθt,(n)cosϕt,(n),sinθt,(n)sinϕt,(n),cosθt,(n)).\displaystyle=\left(\sin{\theta_{t,\ell}^{(n)}}\cos{\phi_{t,\ell}^{(n)}},\sin{\theta_{t,\ell}^{(n)}}\sin{\phi_{t,\ell}^{(n)}},\cos{\theta_{t,\ell}^{(n)}}\right). (1)

The CCS representation facilitates mathematical operations such as dot products, providing computational simplicity. In contrast, the SCS representation offers a clearer depiction of the relationship between channel gains and antenna orientation.

III The PAMA Framework

Building upon the system model, this section presents our PAMA framework. At the core of wireless communication is the interaction between transmitting and receiving antennas through electromagnetic fields. An alternating current in the transmitting antenna generates electromagnetic waves, creating an electromagnetic field that propagates through space to the receiver. The receiving antenna interacts with this field, specifically sensing the electric field component, which induces a current within the antenna used to recover the transmitted signal. This fundamental interaction is governed by Maxwell’s equations, which describe how electric and magnetic fields are generated by charges and currents, and how they propagate in space.

In this section, we begin by modeling the radiation of the electric field from movable transmitting antennas, characterizing how movements – 3D translations and rotations – affect the spatial distribution and polarization of the emitted electromagnetic waves. We then explore how movable receiving antennas detect and interact with the electric field. Ultimately, the analysis leads to a polarization-aware channel representation for movable antennas.

III-A Electric field radiation

Given that the electric field at the receiver is a superposition of waves from all transmitting antennas, we can focus on the radiation analysis of a single antenna. To facilitate the exposition, we introduce the emission angle of a transmitting antenna.

Definition 6 (Emission Angle).

The emission angle θ(e)[0,π]\theta^{(e)}\in[0,\pi] is the angle between the direction of the transmitting antenna and the propagation direction of the electromagnetic wave from the antenna to an observation point 𝐩\bm{p}. The emission angle of the \ell-th transmitting antenna can be written as

θ,𝒑(e)𝒑𝒑t,,𝒏t,,\theta_{\ell,\bm{p}}^{(e)}\triangleq\left\langle\bm{p}-\bm{p}_{t,\ell},\bm{n}_{t,\ell}\right\rangle, (2)

where cosθ,𝐩(e)=(𝐩𝐩t,)𝐧t,|(𝐩𝐩t,)||𝐧t,|\cos{\theta^{(e)}_{\ell,\bm{p}}}=\frac{\left(\bm{p}-\bm{p}_{t,\ell}\right)\cdot\bm{n}_{t,\ell}}{\left|\left(\bm{p}-\bm{p}_{t,\ell}\right)\right||\bm{n}_{t,\ell}|}.

Theorem 2 below summarizes the main result of this section.

Theorem 2 (Electric Field Radiation).

Consider the \ell-th transmitting antenna located at 𝐩t,\bm{p}_{t,\ell} and oriented in the direction 𝐧t,\bm{n}_{t,\ell}. Suppose the transmitted signal at \ell-th antenna is ss_{\ell}. At any observation point 𝐩\bm{p} in space, the electric field induced by this antenna is given by

𝓔,𝒑,𝒑𝒏𝓔,𝒑,\displaystyle\bm{\mathcal{E}}_{\ell,\bm{p}}\triangleq\mathcal{E}_{\ell,\bm{p}}\bm{n}_{\bm{\mathcal{E}}_{\ell,\bm{p}}}, (3)
,𝒑=2jcμ1sej2πλ|𝒑|4π|𝒑|cos(π2cosθ,𝒑(e))sinθ,𝒑(e)ej2πλ𝒑𝒑t,|𝒑|,\displaystyle\mathcal{E}_{\ell,\bm{p}}=2jc\mu_{1}s_{\ell}\frac{e^{-j\frac{2\pi}{\lambda}|\bm{p}|}}{4\pi|\bm{p}|}\frac{\cos{\left(\frac{\pi}{2}\cos{\theta^{(e)}_{\ell,\bm{p}}}\right)}}{\sin{\theta^{(e)}_{\ell,\bm{p}}}}e^{j\frac{2\pi}{\lambda}\frac{\bm{p}\cdot\bm{p}_{t,\ell}}{|\bm{p}|}},
𝒏𝓔,𝒑=𝒏t,(𝒏t,𝒑|𝒑|)𝒑|𝒑||𝒏t,(𝒏t,𝒑|𝒑|)𝒑|𝒑||,\displaystyle\bm{n}_{\bm{\mathcal{E}}_{\ell,\bm{p}}}=\frac{\bm{n}_{t,\ell}-\left(\bm{n}_{t,\ell}\cdot\frac{\bm{p}}{|\bm{p}|}\right)\frac{\bm{p}}{|\bm{p}|}}{\left|\bm{n}_{t,\ell}-\left(\bm{n}_{t,\ell}\cdot\frac{\bm{p}}{|\bm{p}|}\right)\frac{\bm{p}}{|\bm{p}|}\right|},

where ,𝐩\mathcal{E}_{\ell,\bm{p}} denotes the complex amplitude of the induced electric field, incorporating both magnitude and phase; 𝐧𝓔,𝐩\bm{n}_{\bm{\mathcal{E}}_{\ell,\bm{p}}} denotes the polarization direction of the electromagnetic wave at point 𝐩\bm{p}; cc is the speed of light; μ1\mu_{1} is the magnetic permeability of air, λ\lambda is the wavelength.

TABLE II: A summary of the definitions and their geometric relationships.
For \ell-th transmitting antenna For κ\kappa-th receiving antenna
Position 𝒑t,=(xt,,yt,,zt,,)\bm{p}_{t,\ell}=\left(x_{t,\ell},y_{t,\ell},z_{t,\ell},\right) Position 𝒑r,κ=(xr,κ,yr,κ,zr,κ)\bm{p}_{r,\kappa}=\left(x_{r,\kappa},y_{r,\kappa},z_{r,\kappa}\right)
Direction 𝒏t,=(xt,(n),yt,(n),zt,(n),)\bm{n}_{t,\ell}=\left(x_{t,\ell}^{(n)},y_{t,\ell}^{(n)},z_{t,\ell}^{(n)},\right) Polar Angle: θt,(n)\theta_{t,\ell}^{(n)}, Azimuthal Angle: ϕt,(n)\phi_{t,\ell}^{(n)} Direction 𝒏r,κ=(xr,κ(n),yr,κ(n),zr,κ(n),)\bm{n}_{r,\kappa}=\left(x_{r,\kappa}^{(n)},y_{r,\kappa}^{(n)},z_{r,\kappa}^{(n)},\right) Polar Angle: θr,κ(n)\theta_{r,\kappa}^{(n)}, Azimuthal Angle: ϕr,κ(n)\phi_{r,\kappa}^{(n)}
Emission Angle θ,κ(e)\theta_{\ell,\kappa}^{(e)}, with cosθ,κ(e)=𝒑r,κ𝒏t,|𝒑r,κ||𝒏t,|\cos{\theta_{\ell,\kappa}^{(e)}}=\frac{\bm{p}_{r,\kappa}\cdot\bm{n}_{t,\ell}}{|\bm{p}_{r,\kappa}||\bm{n}_{t,\ell}|} Incident Angle θκ(i)\theta_{\kappa}^{(i)}, with sinθκ(i)=𝒑r,κ𝒏r,κ|𝒑r,κ||𝒏r,κ|\sin{\theta_{\kappa}^{(i)}}=\frac{\bm{p}_{r,\kappa}\cdot\bm{n}_{r,\kappa}}{|\bm{p}_{r,\kappa}||\bm{n}_{r,\kappa}|}
Transmitted Angle θκ(t)\theta_{\kappa}^{(t)}, with sinθκ(t)sinθκ(i)=ε1μ1ε2μ2\frac{\sin{\theta_{\kappa}^{(t)}}}{\sin{\theta_{\kappa}^{(i)}}}=\sqrt{\frac{\varepsilon_{1}\mu_{1}}{\varepsilon_{2}\mu_{2}}}
Proof.

We start by considering a reference scenario where the dipole antenna is stationary, located at the origin 𝒑t,=𝒐\bm{p}_{t,\ell}=\bm{o}, and oriented along the direction 𝒏t,=𝒌\bm{n}_{t,\ell}=\bm{k}. When transmitting signal ss_{\ell}, the current density 𝑱\bm{J} applied to the dipole can be written as [19]

𝑱=I0ssin(2πλ(λ4|z|))𝒌,|z|λ4.\displaystyle\bm{J}=I_{0}s_{\ell}\sin{\left(\frac{2\pi}{\lambda}\left(\frac{\lambda}{4}-|z|\right)\right)}\bm{k},\quad|z|\leq\frac{\lambda}{4}. (4)

where I0I_{0} is the maximum current. Without loss of generality, we set I0=1I_{0}=1.

In the mmWave band, signal diffraction and scattering are minimal, making non-line-of-sight (NLOS) paths significantly weaker than the LOS path.111As a concrete example, LOS signals can be over 5050 dB stronger than the NLOS signal at a distance of 100 meters, according to the mmWave channel specified in 3GPP TR 38.901 [26]. Therefore, we only consider the LOS component.

The electric field radiated by the stationary transmitting antenna at any point 𝒑=(xp,yp,zp)\bm{p}=(x_{p},y_{p},z_{p}) in space is given by

𝓔𝒐,𝒑=2jcμ1sej2πλ|𝒑|4π|𝒑|cos(π2cosθ𝒐,𝒑(e))sinθ𝒐,𝒑(e)ϑ,\bm{\mathcal{E}}_{\bm{o},\bm{p}}=2jc\mu_{1}s_{\ell}\frac{e^{-j\frac{2\pi}{\lambda}|\bm{p}|}}{4\pi|\bm{p}|}\frac{\cos{\left(\frac{\pi}{2}\cos{\theta_{\bm{o},\bm{p}}^{(e)}}\right)}}{\sin{\theta_{\bm{o},\bm{p}}^{(e)}}}\bm{\vartheta}, (5)

where θ𝒐,𝒑(e)\theta_{\bm{o},\bm{p}}^{(e)} is the emission angle between 𝒌\bm{k} and 𝒑\bm{p}, and cosθ𝒐,𝒑(e)=𝒑𝒌|𝒑||𝒌|=zp|𝒑|\cos{\theta_{\bm{o},\bm{p}}^{(e)}}=\frac{\bm{p}\cdot\bm{k}}{|\bm{p}||\bm{k}|}=\frac{z_{p}}{|\bm{p}|}. The unit vector ϑ\bm{\vartheta} is the polarization direction at 𝒑\bm{p}. Eq. (5) is a classical result in electromagnetic theory [19, 27] and provides the foundation for our analysis. In Appendix A we provide a detailed derivation of this equation from Maxwell’s equations.

We next consider the antenna movement. After translation and rotation, the \ell-th antenna moves to 𝒑t,\bm{p}_{t,\ell} with a direction 𝒏t,\bm{n}_{t,\ell}. The complex amplitude of the induced electric field follows from (5) and can be written as

,𝒑=2jcμ1sej2πλ|𝒑𝒑t,|4π|𝒑𝒑t,|cos(π2cosθ,𝒑(e))sinθ,𝒑(e),\mathcal{E}_{\ell,\bm{p}}=2jc\mu_{1}s_{\ell}\frac{e^{-j\frac{2\pi}{\lambda}\left|\bm{p}-\bm{p}_{t,\ell}\right|}}{4\pi\left|\bm{p}-\bm{p}_{t,\ell}\right|}\frac{\cos{\left(\frac{\pi}{2}\cos{\theta^{(e)}_{\ell,\bm{p}}}\right)}}{\sin{\theta^{(e)}_{\ell,\bm{p}}}}, (6)

where θ,𝒑(e)\theta^{(e)}_{\ell,\bm{p}} is the emission angle from 𝒑t,\bm{p}_{t,\ell} to 𝒑\bm{p}. Under the far-field condition, 𝒑\bm{p} and 𝒑𝒑t,\bm{p}-\bm{p}_{t,\ell} can be treated parallel, hence we have

|𝒑𝒑t,|\displaystyle\left|\bm{p}-\bm{p}_{t,\ell}\right| \displaystyle\approx |𝒑||𝒑t,|cos𝒑,𝒑t,=|𝒑|1|𝒑|𝒑𝒑t,,\displaystyle|\bm{p}|-\left|\bm{p}_{t,\ell}\right|\cos{\left\langle\bm{p},\bm{p}_{t,\ell}\right\rangle}=|\bm{p}|-\frac{1}{|\bm{p}|}\bm{p}\cdot\bm{p}_{t,\ell},
cosθ,𝒑(e)\displaystyle\cos{\theta^{(e)}_{\ell,\bm{p}}} =\displaystyle= (𝒑𝒑t,)𝒏t,|(𝒑𝒑t,)||𝒏t,|1|𝒑|𝒑𝒏t,.\displaystyle\frac{\left(\bm{p}-\bm{p}_{t,\ell}\right)\cdot\bm{n}_{t,\ell}}{\left|\left(\bm{p}-\bm{p}_{t,\ell}\right)\right||\bm{n}_{t,\ell}|}\approx\frac{1}{|\bm{p}|}\bm{p}\cdot\bm{n}_{t,\ell}.

Therefore, (6) can be refined as

,𝒑2jcμ1sej2πλ|𝒑|4π|𝒑|cos(π2cosθ,𝒑(e))sinθ,𝒑(e)ej2πλ𝒑𝒑t,|𝒑|,\mathcal{E}_{\ell,\bm{p}}\approx 2jc\mu_{1}s_{\ell}\frac{e^{-j\frac{2\pi}{\lambda}|\bm{p}|}}{4\pi|\bm{p}|}\frac{\cos{\left(\frac{\pi}{2}\cos{\theta^{(e)}_{\ell,\bm{p}}}\right)}}{\sin{\theta^{(e)}_{\ell,\bm{p}}}}e^{j\frac{2\pi}{\lambda}\frac{\bm{p}\cdot\bm{p}_{t,\ell}}{|\bm{p}|}}, (7)

where the approximation follows from lim|𝒑t,|/|𝒑|0(|𝒑|𝒑𝒑t,|𝒑|)=|𝒑|\lim\limits_{|\bm{p}_{t,\ell}|/|\bm{p}|\rightarrow 0}\big{(}|\bm{p}|-\frac{\bm{p}\cdot\bm{p}_{t,\ell}}{|\bm{p}|}\big{)}=|\bm{p}|. Note that the term 2πλ\frac{2\pi}{\lambda} in the exponent indicates that the phase varies with distance on the order of the wavelength, hence the term 𝒑𝒑t,|𝒑|\frac{\bm{p}\cdot\bm{p}_{t,\ell}}{|\bm{p}|} in the exponent cannot be neglected.

At point 𝒑\bm{p}, the polarization direction of the electric field is perpendicular to the propagation direction, and lies in the plane formed by the direction of the transmitting antenna 𝒏t,\bm{n}_{t,\ell} and the propagation direction. Thus, the polarization direction 𝒏𝓔,𝒑\bm{n}_{\bm{\mathcal{E}}_{\ell,\bm{p}}} can be obtained by subtracting the component along the propagation direction from the direction of the transmitting antenna, yielding

𝒏𝓔,𝒑=𝒏t,(𝒏t,𝒑|𝒑|)𝒑|𝒑||𝒏t,(𝒏t,𝒑|𝒑|)𝒑|𝒑||.\bm{n}_{\bm{\mathcal{E}}_{\ell,\bm{p}}}=\frac{\bm{n}_{t,\ell}-\left(\bm{n}_{t,\ell}\cdot\frac{\bm{p}}{|\bm{p}|}\right)\frac{\bm{p}}{|\bm{p}|}}{\left|\bm{n}_{t,\ell}-\left(\bm{n}_{t,\ell}\cdot\frac{\bm{p}}{|\bm{p}|}\right)\frac{\bm{p}}{|\bm{p}|}\right|}. (8)

Combining (6) and (8) gives us the induced electric field at observation point 𝒑\bm{p} in (3). ∎

Eq. (3) reveals that, at a given observation point 𝒑\bm{p}, both the complex amplitude and the polarization of the electric field are affected by the movement of the transmitting antenna.

Corollary 3.

The magnitude of the electric field

|,𝒑|2cμ1s4π|𝒑|cos(π2cosθ,𝒑(e))sinθ,𝒑(e)|\mathcal{E}_{\ell,\bm{p}}|\triangleq\frac{2c\mu_{1}s_{\ell}}{4\pi|\bm{p}|}\frac{\cos{\left(\frac{\pi}{2}\cos{\theta^{(e)}_{\ell,\bm{p}}}\right)}}{\sin{\theta^{(e)}_{\ell,\bm{p}}}} (9)

first increases and then decreases as a function of θ,𝐩(e)\theta^{(e)}_{\ell,\bm{p}}, reaching its maximum value when θ,𝐩(e)=π2\theta^{(e)}_{\ell,\bm{p}}=\frac{\pi}{2}. This variation depends solely on the emission angle, which results from the 3D orientation of the transmitting antenna.

Proof.

See Appendix B. ∎

This paper focuses on the half-wave dipole, the most fundamental type of antenna, to develop the PAMA framework. The dipole antenna is directional in the 3D space: when positioned vertically, it radiates uniformly in the horizontal plane but has directional properties in the vertical plane. Our PAMA framework is sufficiently versatile to extend to more complex directional antennas, provided their radiation patterns are known.

III-B Electric Field Sensing

Theorem 2 establishes the electric field induced by a movable transmitting antenna. In this section, we explore how a movable receiving antenna senses the electric field.

The process of capturing radio waves with a receiving antenna involves the transmission of electromagnetic waves from air into the antenna. Dipoles emit linearly polarized radio waves, and when such a plane wave reaches the boundary between the air and the antenna, part of the wave is reflected back into the air, while the remainder is transmitted into the antenna.

Definition 7 (Incident Angle).

The incident angle θ(i)[0,π2]\theta^{(i)}\in[0,\frac{\pi}{2}] is the angle between the direction of the incoming electromagnetic wave and the normal to the surface of the antenna at the point of impact. The direction of the normal 𝐧r,κ\bm{n}_{r,\kappa}^{\perp} can be obtained by subtracting the component that is collinear with the receiving antenna 𝐧r,κ\bm{n}_{r,\kappa} from the direction of the incident electromagnetic wave 𝐩r,κ\bm{p}_{r,\kappa}:

𝒏r,κ=𝒑r,κ(𝒑r,κ𝒏r,κ)𝒏r,κ|𝒑r,κ(𝒑r,κ𝒏r,κ)𝒏r,κ|.\bm{n}_{r,\kappa}^{\perp}=\frac{\bm{p}_{r,\kappa}-(\bm{p}_{r,\kappa}\cdot\bm{n}_{r,\kappa})\bm{n}_{r,\kappa}}{\left|\bm{p}_{r,\kappa}-\left(\bm{p}_{r,\kappa}\cdot\bm{n}_{r,\kappa}\right)\bm{n}_{r,\kappa}\right|}. (10)

We denote by θ,κ(i)\theta_{\ell,\kappa}^{(i)} the incident angle of the wave from the \ell-th transmitting antenna to the κ\kappa-th receiving antenna.

Definition 8 (Transmitted Angle).

The transmitted angle θ(t)[0,π2]\theta^{(t)}\in[0,\frac{\pi}{2}] is the angle between the normal of the receiving antenna and the new propagation direction of the electromagnetic wave inside the antenna after it has entered.

Refer to caption
Figure 4: Illustration of the incident angle θ(i)\theta^{(i)} and transmitted angle θ(t)\theta^{(t)}.

An illustration of the incident angle and transmitted angle is shown in Fig. 4

Definition 9 (Permittivity and Permeability).

Permittivity ε\varepsilon quantifies the ability of a material to store electrical energy in an electric field. Permeability μ\mu is a measure of how easy or difficult it is for a magnetic field to pass through a material [27]. We define ε1\varepsilon_{1} and ε2\varepsilon_{2} as the permittivity of air and antenna, respectively, and μ1\mu_{1} and μ2\mu_{2} as the permeability of air and antenna, respectively. We further define a relative permittivity

εrε2ε1>1.\varepsilon_{r}\triangleq\frac{\varepsilon_{2}}{\varepsilon_{1}}>1. (11)
Lemma 4 (Snell’s law [28]).

Let ε1\varepsilon_{1} and ε2\varepsilon_{2} be the permittivity of air and antenna, respectively, and μ1\mu_{1} and μ2\mu_{2} be the permeability of air and antenna, respectively. We have

sinθ(t)sinθ(i)=ε1μ1ε2μ2.\frac{\sin{\theta^{(t)}}}{\sin{\theta^{(i)}}}=\sqrt{\frac{\varepsilon_{1}\mu_{1}}{\varepsilon_{2}\mu_{2}}}. (12)

A critical aspect of electric field sensing is polarization matching, which refers to the alignment of the transmitting and receiving antennas’ polarization states. The efficiency of polarization matching determines how much energy is transferred into the antenna. Therefore, effective polarization matching is essential for maximizing the power transfer and minimizing signal degradation, especially in systems with movable antennas where orientation can dynamically change.

Definition 10 (Polarization Matching Angle).

The polarization matching angle α[0,π]\alpha\in[0,\pi] is the angle between the polarization direction of the incident wave and the direction of the receiving antenna. We denote by α,κ=𝐧𝓔,𝐩r,κ,𝐧r,κ\alpha_{\ell,\kappa}=\langle\bm{n}_{\bm{\mathcal{E}}_{\ell,\bm{p}_{r,\kappa}}},\bm{n}_{r,\kappa}\rangle the polarization matching angle of the wave emitted from the \ell-th transmitting antenna to the κ\kappa-th receiving antenna.

Theorem 5 below summarizes the main result of this section.

Theorem 5 (Electric Field Sensing).

Consider the signal transmitted from the \ell-th transmitting antenna. The received signal at the κ\kappa-th user’ antenna can be written as

r,κ=,𝒑r,κAF(θ,κ(i),α,κ),\displaystyle r_{\ell,\kappa}=\frac{\mathcal{E}_{\ell,\bm{p}_{r,\kappa}}}{A_{F}}\mathcal{M}\left(\theta_{\ell,\kappa}^{(i)},\alpha_{\ell,\kappa}\right), (13)
(θ,κ(i),α,κ)=1|Γ,κ|2cos2α,κ|Γ,κ|2sin2α,κ,\displaystyle\mathcal{M}\left(\theta_{\ell,\kappa}^{(i)},\alpha_{\ell,\kappa}\right)\!=\!\sqrt{1\!-\!\left|\Gamma_{\parallel,\kappa}\right|^{2}\cos^{2}{\alpha_{\ell,\kappa}}\!-\!\left|\Gamma_{\perp,\kappa}\right|^{2}\sin^{2}{\alpha_{\ell,\kappa}}},

where ,𝐩r,κ\mathcal{E}_{\ell,\bm{p}_{r,\kappa}} is the complex amplitude of the induced electric field by the \ell-th antenna at 𝐩r,κ\bm{p}_{r,\kappa}; the constant AFA_{F} is an antenna factor for converting field amplitude to voltage; (θ,κ(i),α,κ)\mathcal{M}\left(\theta_{\ell,\kappa}^{(i)},\alpha_{\ell,\kappa}\right) denotes the polarization matching efficiency, considering both the polarization matching angle and energy reflection; Γ,κ\Gamma_{\parallel,\kappa} and Γ,κ\Gamma_{\perp,\kappa} are the energy reflection coefficients for the horizontal and vertical polarization:

Γ,κ=εr1+cos2θ,κ(i)εrcosθ,κ(i)εr1+cos2θ,κ(i)+εrcosθ,κ(i),\displaystyle\Gamma_{\parallel,\kappa}=\frac{\sqrt{\varepsilon_{r}-1+\cos^{2}{\theta_{\ell,\kappa}^{(i)}}}-\varepsilon_{r}\cos{\theta_{\ell,\kappa}^{(i)}}}{\sqrt{\varepsilon_{r}-1+\cos^{2}{\theta_{\ell,\kappa}^{(i)}}}+\varepsilon_{r}\cos{\theta_{\ell,\kappa}^{(i)}}}, (14)
Γ,κ=εr1+cos2θ,κ(i)cosθ,κ(i)εr1+cos2θ,κ(i)+cosθ,κ(i).\displaystyle\Gamma_{\perp,\kappa}=\frac{\sqrt{\varepsilon_{r}-1+\cos^{2}{\theta_{\ell,\kappa}^{(i)}}}-\cos{\theta_{\ell,\kappa}^{(i)}}}{\sqrt{\varepsilon_{r}-1+\cos^{2}{\theta_{\ell,\kappa}^{(i)}}}+\cos{\theta_{\ell,\kappa}^{(i)}}}. (15)
Proof.

The electric field induced by the \ell-th transmitting antenna possesses a polarization direction 𝒏𝓔,𝒑r,κ\bm{n}_{\bm{\mathcal{E}}_{\ell,\bm{p}_{r,\kappa}}}. The κ\kappa-th receiving antenna, oriented at 𝒏r,κ\bm{n}_{r,\kappa}, has a polarization matching angle α,κ\alpha_{\ell,\kappa} w.r.t. the incident wave. We can partition the incident polarization into horizontal and vertical components.

Due to the electromagnetic boundary conditions at the interface of different media (air to antenna), part of the energy of each polarization component is reflected. The reflection coefficients Γ,κ\Gamma_{\parallel,\kappa} and Γ,κ\Gamma_{\perp,\kappa}, corresponding to horizontal and vertical polarizations respectively, determine how much of the incident wave’s energy is transmitted into the receiver.

The reflection coefficient Γ\Gamma for any polarization can generally be expressed in terms of the characteristic impedances of the interacting media [27, 21]:

Γ=Z2Z1Z2+Z1,\Gamma=\frac{Z_{2}-Z_{1}}{Z_{2}+Z_{1}}, (16)

where

Z1={120πμ1ε1cosθ,κ(i),for horizontal polarization,120πμ1ε1/cosθ,κ(i),for vertical polarization,\displaystyle Z_{1}=\begin{cases}120\pi\sqrt{\frac{\mu_{1}}{\varepsilon_{1}}}\cos{\theta^{(i)}_{\ell,\kappa}},&\text{for horizontal polarization},\\ 120\pi\sqrt{\frac{\mu_{1}}{\varepsilon_{1}}}/\cos{\theta^{(i)}_{\ell,\kappa}},&\text{for vertical polarization},\end{cases}
Z2={120πμ2ε2cosθ,κ(t),for horizontal polarization,120πμ2ε2/cosθ,κ(t),for vertical polarization.\displaystyle Z_{2}=\begin{cases}120\pi\sqrt{\frac{\mu_{2}}{\varepsilon_{2}}}\cos{\theta^{(t)}_{\ell,\kappa}},&\text{for horizontal polarization},\\ 120\pi\sqrt{\frac{\mu_{2}}{\varepsilon_{2}}}/\cos{\theta^{(t)}_{\ell,\kappa}},&\text{for vertical polarization}.\end{cases}

Given (16), we compute the reflection coefficients for the horizontal and vertical components:

Γ=120πμ2ε21sin2θ,κ(t)120πμ1ε1cosθ,κ(i)120πμ2ε21sin2θ,κ(t)+120πμ1ε1cosθ,κ(i)\displaystyle\Gamma_{\parallel}=\frac{120\pi\sqrt{\frac{\mu_{2}}{\varepsilon_{2}}}\sqrt{1-\sin^{2}{\theta^{(t)}_{\ell,\kappa}}}-120\pi\sqrt{\frac{\mu_{1}}{\varepsilon_{1}}}\cos{\theta^{(i)}_{\ell,\kappa}}}{120\pi\sqrt{\frac{\mu_{2}}{\varepsilon_{2}}}\sqrt{1-\sin^{2}{\theta^{(t)}_{\ell,\kappa}}}+120\pi\sqrt{\frac{\mu_{1}}{\varepsilon_{1}}}\cos{\theta^{(i)}_{\ell,\kappa}}}
=(a)εrsin2θ,κ(i)εr1sin2θ,κ(i)εrsin2θ,κ(i)+εr1sin2θ,κ(i),\displaystyle\overset{(a)}{=}\frac{\sqrt{\varepsilon_{r}-\sin^{2}{\theta^{(i)}_{\ell,\kappa}}}-\varepsilon_{r}\sqrt{1-\sin^{2}{\theta^{(i)}_{\ell,\kappa}}}}{\sqrt{\varepsilon_{r}-\sin^{2}{\theta^{(i)}_{\ell,\kappa}}}+\varepsilon_{r}\sqrt{1-\sin^{2}{\theta^{(i)}_{\ell,\kappa}}}}, (17)
Γ=εrsin2θ,κ(i)1sin2θ,κ(i)εrsin2θ,κ(i)+1sin2θ,κ(i),\displaystyle\Gamma_{\perp}=\frac{\sqrt{\varepsilon_{r}-\sin^{2}{\theta^{(i)}_{\ell,\kappa}}}-\sqrt{1-\sin^{2}{\theta^{(i)}_{\ell,\kappa}}}}{\sqrt{\varepsilon_{r}-\sin^{2}{\theta^{(i)}_{\ell,\kappa}}}+\sqrt{1-\sin^{2}{\theta^{(i)}_{\ell,\kappa}}}}, (18)

where (a) follows because the permeability of most materials remains constant and is unaffected by changes in frequency or temperature [21]. Therefore, we have μ1=μ2\mu_{1}=\mu_{2}.

Refer to caption
Figure 5: Relationships between the PAMA system, the electric field, and the PAMA channel gain under the far-field condition with a single LOS path. (a) the impact of rotating the transmitting antenna on the radiation pattern relative to the CCS; (b) variation in received signal strength as a function of the emission angle; (c), (d): changes in polarization matching efficiency based on the incident angle θ(i)\theta^{(i)} and the polarization matching angle α\alpha.

Without loss of generality, this paper represents the received signal by the voltage induced at the receiver. In practice, the difference between the amplitude of the received signal in terms of voltage and current is simply a scaling factor: the resistance, which can be considered part of the antenna factor. Therefore, the received signal, i.e., the induced voltage, at the κ\kappa-th user from the \ell-th transmitting antenna can be written as

r,κ=V0=,𝒑r,κAF(θ,κ(i),α,κ),\displaystyle r_{\ell,\kappa}=V_{0}=\frac{\mathcal{E}_{\ell,\bm{p}_{r,\kappa}}}{A_{F}}\mathcal{M}\left(\theta_{\ell,\kappa}^{(i)},\alpha_{\ell,\kappa}\right),

where the constant AFA_{F} is the antenna factor for converting field amplitude to voltage, as shown in (13). ∎

III-C Polarization-aware channel gain

Building upon Theorems 2 and 5, this section characterizes the polarization-aware channel gain for movable transmitting and receiving antennas.

Theorem 6 (PAMA Channel Gain).

Consider the transmission from the \ell-th transmitting antenna to the κ\kappa-th receiving antenna. The polarization-aware channel gain is given by

hκ,(𝒑t,,𝒏t,,𝒑r,κ,𝒏r,)=r,κs\displaystyle h_{\kappa,\ell}\left(\bm{p}_{t,\ell},\bm{n}_{t,\ell},\bm{p}_{r,\kappa},\bm{n}_{r,\ell}\right)=\frac{r_{\ell,\kappa}}{s_{\ell}} (19)
=Gκcos(π2cosθ,κ(e))sinθ,κ(e)(θk(i),α,κ)ej2πλ𝒑r,κ𝒑t,|𝒑r,κ|,\displaystyle=G_{\kappa}\frac{\cos{\left(\frac{\pi}{2}\cos{\theta^{(e)}_{\ell,\kappa}}\right)}}{\sin{\theta^{(e)}_{\ell,\kappa}}}\mathcal{M}\left(\theta_{k}^{(i)},\alpha_{\ell,\kappa}\right)e^{j\frac{2\pi}{\lambda}\frac{\bm{p}_{r,\kappa}\cdot\bm{p}_{t,\ell}}{|\bm{p}_{r,\kappa}|}},

where Gκ2jcμ1AFej2πλ|𝐩r,κ|4π|𝐩r,κ|G_{\kappa}\triangleq\frac{2jc\mu_{1}}{A_{F}}\frac{e^{-j\frac{2\pi}{\lambda}|\bm{p}_{r,\kappa}|}}{4\pi|\bm{p}_{r,\kappa}|}.

The PAMA channel gain model captures the changes induced by the radiation pattern, the polarization matching efficiency, and phase shifts as the antenna translates and rotates. Compared to existing channel models for MAs that primarily focus on the phase changes, our PAMA model offers a more comprehensive and rigorous framework to explore the potential of MAs, enabling deeper analysis and performance assessment.

IV Impacts of Antenna Movements

This section analyzes the PAMA channel model and examines how each variable influences channel quality, providing insights into the key characteristics of PAMA and identifying the sources of its performance gains. Fig. 5 summarizes the relationships between antenna movement and their corresponding impacts on the channel. It is essential to highlight that these properties are derived under the conditions of a far-field scenario and a dominant LOS path. For clarity, the indices \ell and κ\kappa in the subscripts are omitted in this section.

Proposition 1.

Antenna translation only contributes a phase term to the PAMA channel gain. This implies that the effect of antenna translation can be fully realized by precoding.

In the far-field LOS scenario, the electromagnetic waves can be approximated as plane waves, and the spatial fading characteristics are negligible. The relative orientation of the electric field vectors between the transmitter and receiver remains unchanged despite antenna translation. This means that antenna translation does not affect the amplitude or polarization alignment of the received signal – only the phase changes due to the difference in path length, leading to Proposition 1.

In practical systems, digital precoding is often easier to implement than mechanical antenna movements. Therefore, when studying MAs, it is essential to consider the unique advantages that MAs introduce beyond what precoding can achieve. Proposition 1 indicates that the benefits of antenna translation are more significant in near-field or rich multi-path environments, in which case antenna translation can impact the amplitude and spatial characteristics of the channel gain, offering advantages that precoding alone cannot replicate.

In the context of this paper, which focuses on far-field and mmWave dominant LOS paths, the unique gains that MAs can offer over precoding are primarily obtained through antenna rotation rather than translation.

Proposition 2.

The rotation of the transmitting antenna alters the amplitude of the channel gain by modifying both the radiation pattern of the induced electric field and the polarization matching efficiency at the receiver.

When the transmitting antenna rotates, it changes the emission angle θ(e)\theta^{(e)} of the electromagnetic wave. This alteration affects the antenna’s radiation pattern – the distribution of the electric field w.r.t. the CCS – which in turn influences the electric field intensity received. As illustrated in Fig. 5 (a) and (b), different emission angles lead to variations in the radiation pattern, resulting in changes to the received signal’s strength. According to Corollary 3, the magnitude or intensity of the electric field reaches its maximum when the emission angle is θ(e)=π2\theta^{(e)}=\frac{\pi}{2}. To achieve this, the transmitting antenna should be positioned within the plane that is perpendicular to the direction of wave propagation.

Furthermore, the rotation of the transmitting antenna also modifies its polarization orientation, thereby changing the polarization matching angle α\alpha between the transmitter and receiver, and hence the matching efficiency (θ(i),α)\mathcal{M}\left(\theta^{(i)},\alpha\right) in the PAMA channel gain.

Proposition 3.

The rotation of the receiving antennas alters the amplitude of the channel gain by modifying the polarization matching efficiency.

Rotating the receiving antennas affects both the incident angle θ(i)\theta^{(i)} of the incoming electromagnetic wave and the polarization matching angle α\alpha, both of which impacts the polarization matching efficiency.

  1. 1.

    The change of the incident angle θ(i)\theta^{(i)} alters the energy reflection coefficients Γ\Gamma_{\parallel} and Γ\Gamma_{\perp}, as defined in (14) and (15). These reflection coefficients characterize the proportion of electromagnetic energy that enters the antenna. A smaller reflection coefficient indicates that more energy is transmitted into the antenna, enhancing the received signal strength.

  2. 2.

    The orientation of the receiving antenna relative to the electric field vector of the incoming wave determines the polarization matching angle α\alpha, which describes the ratio of the horizontal to vertical components of the electromagnetic wave’s polarization relative to the antenna’s orientation. A better alignment (i.e., a smaller mismatch angle) leads to higher polarization matching efficiency.

Figs. 5 (c) and (d) illustrate how the polarization matching efficiency varies with the incident angle θ(i)\theta^{(i)} and the polarization matching angle α\alpha for antennas with different relative dielectric constants εr\varepsilon_{r}. As can be seen, without proper optimization of the receiving antenna’s rotation, the matching efficiency can decrease significantly, potentially dropping to zero. In such cases, the antenna fails to effectively receive the incoming signal, resulting in a loss of information. Therefore, careful control of the receiving antenna’s orientation is essential for maximizing the channel gain and ensuring reliable communication in the PAMA framework.

Refer to caption
(a) |h|2|h|^{2}
Refer to caption
(b) Geometry relationship.
Figure 6: Channel amplitude variations and the orientations of the transmitting antenna in Scenario 1, where 𝒑t,1=[0,0,0]\bm{p}_{t,1}=[0,0,0], 𝒑r,1=[75,40,50]\bm{p}_{r,1}=[75,-40,50], and 𝒏r,1=[0,0,1]\bm{n}_{r,1}=[0,0,1].
Refer to caption
(a) |h|2|h|^{2}
Refer to caption
(b) Geometry relationship.
Figure 7: Channel amplitude variations and the orientations of the receiving antenna in Scenario 2, where 𝒑t,1=[0,0,0]\bm{p}_{t,1}=[0,0,0], 𝒏t,1=[0,0,1]\bm{n}_{t,1}=[0,0,1], and 𝒑r,1=[75,40,50]\bm{p}_{r,1}=[75,-40,50].

Propositions 2 and 3 indicate that the amplitude of the channel gain is influenced by the orientations of both the transmitting and receiving antennas. To illustrate the extent to which antenna rotations affect the channel amplitude, we simulate two scenarios:

  • Scenario 1: Rotatable transmitting antenna with fixed receiving antenna;

  • Scenario 2: Fixed transmitting antenna with rotatable receiving antenna.

In the simulations, the positions of the antennas are set as 𝒑t=[0,0,0]\bm{p}_{t}=[0,0,0], 𝒑r=[75,40,50]\bm{p}_{r}=[75,-40,50]. The fixed antennas are oriented vertically upward, with a direction vector of [0,0,1][0,0,1]. For visualization purposes, we vary the polar and azimuthal angles of the rotatable antennas over two periods, even though their actual values lie within a single period. The channel amplitude variations for the two scenarios are presented in Figs. 6(a) and 7(a), respectively.

For scenario 1 with a fixed receiving antenna, Fig. 6(a) shows that a randomly oriented transmitting antenna has only a 67.5%67.5\% probability of delivering at least half of the maximum possible energy to the user. Each period can yield two sets of maximum values due to the symmetry of the dipole antennas, where the collinear reverse sets of polar and azimuthal angles provide equivalent performance.

For scenario 2 with a fixed transmitting antenna, Fig. 7(a) shows that a randomly oriented receiving antenna has a 99.0%99.0\% probability of capturing at least half of the maximum energy. However, the remaining 1.0%1.0\% of orientations present a significant risk of receiving virtually no energy. This outcome highlights the critical importance of antenna orientation in ensuring reliable communication.

Figs. 6(b) and 7(b) depict the positions and orientations of the transmitting and receiving antennas for both scenarios. In these illustrations, fixed Antennas are shown with vertical black arrows; the optimal orientations of the rotatable antennas are shown with solid red arrows, indicating the directions that maximize the channel gain; the worst-case orientations of the rotatable antennas are shown with dashed green arrows, representing the directions that minimize the channel gain.

Proposition 4.

In a single-input single-output (SISO) PAMA system:

  1. 1.

    Optimal orientation: The best channel gain is achieved when the transmitting and receiving antennas are oriented such that they, along with the propagation path, lie within the same plane.

  2. 2.

    Worst orientation: The poorest channel gain occurs when one antenna points directly at the other.

When the electromagnetic field reaches the receiving antenna, it can be decomposed into three orthogonal components along: the direction of the receiving antenna 𝒏r\bm{n}_{r}; the normal direction of the incident signal 𝒏r\bm{n}_{r}^{\perp}; a direction perpendicular to both 𝒏r\bm{n}_{r} and 𝒏r\bm{n}_{r}^{\perp}, denoted by 𝒏3\bm{n}_{3}. This decomposition is illustrated in Fig. 8. The component along 𝒏r\bm{n}_{r} represents the horizontal polarization component, while the components along 𝒏r\bm{n}_{r}^{\perp} and 𝒏3\bm{n}_{3} represent the vertical polarization components.

Refer to caption
Figure 8: The electromagnetic field at the receiving antenna can be decomposed into three orthogonal components: along the direction of the receiving antenna 𝒏r\bm{n}_{r}, the normal direction of the incident signal 𝒏r\bm{n}_{r}^{\perp}, and a third direction 𝒏3\bm{n}_{3} that is perpendicular to both 𝒏r\bm{n}_{r} and 𝒏r\bm{n}_{r}^{\perp}.

According to (14) and (15), we have ΓΓ\Gamma_{\parallel}\geq\Gamma_{\perp}, meaning that the horizontal polarization component has a smaller reflection coefficient, and thus a better matching efficiency, than the vertical component. Therefore, increasing the component along 𝒏r\bm{n}_{r} by adjusting antenna orientations enhances the channel gain. However, achieving complete horizontal polarization requires the polarization matching angle α\alpha to be zero, which also corresponds to an incident angle θ(i)\theta^{(i)} of zero. Under this condition, the reflection coefficient Γ\Gamma_{\parallel} increases, reducing the matching efficiency. Therefore, simply minimizing α\alpha is suboptimal. Instead, we need to balance the horizontal and vertical polarization components to maximize the overall matching efficiency.

Since the polarization direction 𝒏𝓔\bm{n}_{\bm{\mathcal{E}}} is perpendicular to the propagation direction, and 𝒑r\bm{p}_{r} is perpendicular to the antenna direction, the polarization matching angle satisfies α=𝒏𝓔,𝒏rθ(i)=𝒏r,𝒑r\alpha\allowbreak=\langle\bm{n}_{\bm{\mathcal{E}}},\allowbreak\bm{n}_{r}\rangle\allowbreak\geq\theta^{(i)}=\allowbreak\langle\bm{n}_{r}^{\perp},\allowbreak\bm{p}_{r}\rangle. To increase the incident angle and reduce the reflection coefficient without introducing additional vertical polarization components along 𝒏3\bm{n}_{3}, it is optimal to maintain α=θ(i)\alpha=\theta^{(i)}. For linearly polarized antennas, this configuration implies that the transmitting antenna, the receiving antenna, and the propagation path are all aligned within the same plane, confirming Proposition 4. It is important to note that lying in the same plane does not necessarily mean that the transmitting and receiving antennas are parallel, as depicted in Figs. 6(a) and 7(a).

Figs. 6(b) and 7(b) demonstrate that the worst-case orientation consistently occurs when one antenna is directed toward the other. This phenomenon is explained by:

  • For the transmitting antenna in Figs. 6(b), when the emission angle θ(e)=0\theta^{(e)}=0, meaning the transmitting antenna points directly at the receiver, the electric field intensity in that direction is effectively zero due to the dipole antenna’s radiation pattern.

  • For the receiving antenna in Figs. 7(b), when the incident angle θ(i)=π2\theta^{(i)}=\frac{\pi}{2}, the reflection coefficients Γ,1=Γ,1=1\Gamma_{\parallel,1}=\Gamma_{\perp,1}=1, indicating total reflection with no energy entering the antenna. Additionally, the polarization matching angle α=π2\alpha=\frac{\pi}{2}, indicating complete polarization mismatch and resulting in zero capturable energy.

V Optimization with PAMA

By incorporating polarization, we have redefined movable antennas and established more refined, polarization-aware channel models. This enhancement serves as a foundation for further optimizations, specifically tailored to harness the dynamic adjustments of antenna orientations and positions to maximize system performance. In this paper, we focus on the fundamental sum-rate optimization [13] to reveal the potential of PAMA.

V-A Sum-rate optimization

We consider a multi-user multiple input single output (MU-MISO) configuration where the BS serves KK users simultaneously, KLK\leq L, using transmit beamforming. Let the source data symbols intended for the users be 𝒙=[x1,x2,,xK]\bm{x}=[x_{1},x_{2},\dots,x_{K}]^{\top}, where the symbol power is normalized such that |xκ|2=1|x_{\kappa}|^{2}=1, κ=1,2,,K\kappa=1,2,\dots,K. The beamforming vector for the κ\kappa-th user is denoted by 𝒘κL×1\bm{w}_{\kappa}\in\mathbb{C}^{L\times 1}, with individual elements w,κw_{\ell,\kappa}. A power scaling matrix, defined as 𝑷=diag{P1,P2,,PK}\bm{P}=\text{diag}\{P_{1},P_{2},\dots,P_{K}\}, is employed to control the power distribution among the users. After beamforming, the signal transmitted from the BS antennas, denoted by 𝒔[s1,s2,,sL]\bm{s}\triangleq[s_{1},s_{2},\dots,s_{L}]^{\top}, is written as

s=κ=1KPκw,κxκ,\displaystyle s_{\ell}=\sum_{\kappa=1}^{K}\sqrt{P_{\kappa}}w_{\ell,\kappa}x_{\kappa}, (20)

subject to the constraint κ=1KPκ=P\sum_{\kappa=1}^{K}P_{\kappa}=P, where PP denotes the total power constraint.

The symbols received by the users, denoted by 𝒚[y1,y2,,yK]\bm{y}\triangleq[y_{1},y_{2},\dots,y_{K}]^{\top}, are given by

𝒚=𝑯𝑾𝑷𝒙+𝒛,\displaystyle\bm{y}=\bm{H}\bm{W}\bm{P}\bm{x}+\bm{z}, (21)

where 𝑯=[𝒉1,𝒉2,,𝒉K]\bm{H}=[\bm{h}_{1}^{\top},\bm{h}_{2}^{\top},\dots,\bm{h}_{K}^{\top}]^{\top} is the channel matrix with the elements defined by (19); 𝑾[𝒘1,𝒘2,,𝒘K]\bm{W}\triangleq[\bm{w}_{1},\bm{w}_{2},...,\bm{w}_{K}] is the precoding matrix; 𝒛=[z1,z2,,zK]\bm{z}=[z_{1},z_{2},\dots,z_{K}]^{\top} is the additive white Gaussian noise (AWGN) vector, where zk𝒞𝒩(0,σκ2)z_{k}\sim\mathcal{CN}(0,\sigma^{2}_{\kappa}).

At each user κ\kappa, the SINR is defined as

γκ=Pκ|𝒉κ𝒘κ|2σκ2+j=1,jκK(Pj|𝒉κ𝒘j|2),\displaystyle\gamma_{\kappa}=\frac{P_{\kappa}\left|\bm{h}_{\kappa}\bm{w}_{\kappa}\right|^{2}}{\sigma_{\kappa}^{2}+\sum_{j=1,j\neq\kappa}^{K}{\left(P_{j}\big{|}\bm{h}_{\kappa}\bm{w}_{j}\big{|}^{2}\right)}}, (22)

and the achievable rate is

Rκ=12log2(1+γκ).\displaystyle R_{\kappa}=\frac{1}{2}\log_{2}{\left(1+\gamma_{\kappa}\right)}. (23)

The average achievable rate across all users is

R=1Kk=1KRκ12log2(1+γtotal),R=\frac{1}{K}\sum\limits_{k=1}^{K}{R_{\kappa}}\triangleq\frac{1}{2}\log_{2}{\left(1+\gamma_{\text{total}}\right)}, (24)

with the equivalent total SINR for the system:

γtotal=[κ=1K(1+γκ)]1K1.\gamma_{\text{total}}=\left[\prod_{\kappa=1}^{K}\left(1+\gamma_{\kappa}\right)\right]^{\frac{1}{K}}-1. (25)

Our goal is to maximize the equivalent total SINR of the system by jointly optimizing the positions 𝒑t,\bm{p}_{t,\ell} and orientations 𝒏t,\bm{n}_{t,\ell} of the transmitting antennas, the orientations 𝒏r,κ\bm{n}_{r,\kappa} of the receiving antennas, beamforming matrix 𝑾\bm{W}, and the power scaling matrix 𝑷\bm{P}. The optimization problem is formulated as

(P1): max𝑾,𝑷,{𝒏t,},{𝒑t,},{𝒏r,κ}γtotal,\displaystyle\max_{\bm{W},\bm{P},\{\bm{n}_{t,\ell}\},\{\bm{p}_{t,\ell}\},\{\bm{n}_{r,\kappa}\}}\gamma_{\text{total}}, (26a)
s.t. κ=1KPκ=P\displaystyle\sum_{\kappa=1}^{K}P_{\kappa}=P (26b)
|𝒑t,𝒑t,i|λ2,1iL,i,\displaystyle\left|\bm{p}_{t,\ell}-\bm{p}_{t,i}\right|\leq\frac{\lambda}{2},\quad\forall 1\leq i\leq L,i\neq\ell, (26c)
𝒑t,𝒞t.\displaystyle\bm{p}_{t,\ell}\in\mathcal{C}_{t}. (26d)

In constraint (26c), a minimum separation of λ2\frac{\lambda}{2} between any two antennas is imposed to prevent collisions. The term 𝒞t\mathcal{C}_{t} in (26d) represents the movable region of the transmitting antennas.

To solve (P1), we employ zero-forcing (ZF) beamforming, which can fully eliminate the inter-channel interferences [29]. The ZF beamformer is given by

𝑾[𝒘1,𝒘2,,𝒘K]=𝑯H(𝑯𝑯H)1,\bm{W}^{\prime}\triangleq[\bm{w}_{1}^{\prime},\bm{w}_{2}^{\prime},\dots,\bm{w}_{K}^{\prime}]=\bm{H}^{H}\left(\bm{H}\bm{H}^{H}\right)^{-1},
𝑾=[𝒘1|𝒘1|,𝒘2|𝒘2|,,𝒘K|𝒘K|].\bm{W}=\left[\frac{\bm{w}_{1}^{\prime}}{|\bm{w}_{1}^{\prime}|},\frac{\bm{w}_{2}^{\prime}}{|\bm{w}_{2}^{\prime}|},\dots,\frac{\bm{w}_{K}^{\prime}}{|\bm{w}_{K}^{\prime}|}\right]. (27)

Since ZF beamforming completely removes interferences, the power allocation matrix can be easily obtained by water-filling [30], yielding

Pκ=max{γth1γκ,0},\displaystyle P_{\kappa}=\max{\left\{\gamma_{th}-\frac{1}{\gamma_{\kappa}},~{}0\right\}}, (28)

where γth\gamma_{th} is the water level that can be determined iteratively to satisfy the total power constraint.

Next, we transform the non-convex feasible region of 𝒑t,\bm{p}_{t,\ell} into a convex set by linearizing the non-convex constraint in (26c) [4]. During the tt-th iteration of the gradient descent, considering the position 𝒑t,(t1)\bm{p}_{t,\ell}^{(t-1)} from the previous iteration and the position 𝒑t,i\bm{p}_{t,i} (1iL1\leq i\leq L, ii\neq\ell), we linearize by determining the intersection of the displacement vector 𝒑t,(t1)𝒑t,i\bm{p}_{t,\ell}^{(t-1)}-\bm{p}_{t,i} with the spherical boundary |𝒑t,ib𝒑t,i|=λ2\left|\bm{p}_{t,i}^{\text{b}}-\bm{p}_{t,i}\right|=\frac{\lambda}{2}. The boundary point is given by

𝒑ib=𝒑t,i+λ2𝒑t,(t1)𝒑t,i|𝒑t,(t1)𝒑t,i|.\displaystyle\bm{p}_{\ell i}^{\text{b}}=\bm{p}_{t,i}+\frac{\lambda}{2}\frac{\bm{p}_{t,\ell}^{(t-1)}-\bm{p}_{t,i}}{\left|\bm{p}_{t,\ell}^{(t-1)}-\bm{p}_{t,i}\right|}.

A hyperplane is then defined w.r.t. 𝒑ib\bm{p}_{\ell i}^{\text{b}} and the vector 𝒑t,(t1)𝒑t,i\bm{p}_{t,\ell}^{(t-1)}-\bm{p}_{t,i} as

{𝒘|(𝒑ib𝒑t,i)𝒘=c},\displaystyle\left\{\bm{w}\left|\left(\bm{p}_{\ell i}^{\text{b}}-\bm{p}_{t,i}\right)^{\top}\bm{w}=c\right.\right\},

and the corresponding closed halfspace divided by this hyperplane is given by

{𝒘|(𝒑t,(t1)𝒑t,i)(𝒘𝒑ib)0},\displaystyle\left\{\bm{w}\left|\left(\bm{p}_{t,\ell}^{(t-1)}-\bm{p}_{t,i}\right)^{\top}\left(\bm{w}-\bm{p}_{\ell i}^{\text{b}}\right)\geq 0\right.\right\},

which includes 𝒑ib\bm{p}_{\ell i}^{\text{b}} and any vector that forms an acute angle with the vector 𝒑t,(t1)𝒑t,i\bm{p}_{t,\ell}^{(t-1)}-\bm{p}_{t,i}.

By setting 𝒘=𝒑t,\bm{w}=\bm{p}_{t,\ell} in (V-A), the non-convex constraint (26c) can be approximated as the following linear inequality:

(𝒑t,(t1)𝒑t,i)(𝒑t,𝒑ib)0,1iL,i.\displaystyle\left(\bm{p}_{t,\ell}^{(t-1)}-\bm{p}_{t,i}\right)^{\top}\left(\bm{p}_{t,\ell}-\bm{p}_{\ell i}^{\text{b}}\right)\geq 0,~{}\forall 1\leq i\leq L,i\neq\ell. (29)

Thus, problem (P1) is redefined to a more tractable form

(P2): max{𝒏t,},{𝒑t,},{𝒏r,κ}γtotal,\displaystyle\max_{\{\bm{n}_{t,\ell}\},\{\bm{p}_{t,\ell}\},\{\bm{n}_{r,\kappa}\}}\gamma_{\text{total}}, (30a)
s.t. (27),(28),(26b),(29),(26d).\displaystyle\eqref{e:w},~{}\eqref{e:P},~{}\eqref{e:P1power},~{}\eqref{e:P1ref2},~{}\eqref{e:P1region}. (30b)

Problem (P2) can be solved efficiently using alternating optimization and projected gradient descent.

Remark 1 (Translation versus beamforming).

With ZF beamforming, the phase shifts introduced by antenna translation are entirely compensated. The equivalent channel matrix 𝐇𝐖\bm{HW} reduces to a diagonal form, with each diagonal element given by 1|𝐰κ|\frac{1}{|\bm{w}^{\prime}_{\kappa}|}.

In this MU-MISO scenario, the effect of antenna translation can be fully replicated through precoding. Therefore, the optimization problem (P1) can be solved without considering the optimization of the antenna positions {𝐩t,}\{\bm{p}_{t,\ell}\}. The constraints (26c) and (26d) can also be omitted. This conclusion is supported by our simulations depicted in Fig. 10, which show that antenna translation brings no benefits during the optimization of (P2).

Remark 2 (Channel estimation in PAMA).

Channel estimation poses a significant challenge in MA systems because determining the optimal antenna position requires knowledge of the channel characteristics at every possible position and orientation within the movement region 𝒞t\mathcal{C}_{t}. This task is particularly complex in radio frequency bands with rich multipath environments.

However, within our PAMA framework operating in the mmWave band, channel estimation is considerably simplified due to the dominance of the LOS path. PAMA only requires information about user positions and antenna orientations. With this information, we can accurately infer the channel gains at any position and orientation using (19).

Remark 3 (The K>LK>L case).

ZF beamforming effectively eliminates inter-user interference only when the number of users KK does not exceed the number of transmit antennas LL. When K>LK>L, we can partition the users into groups of size up to LL and assign orthogonal resources to each group. Alternatively, exploring different optimization criteria, such as minimizing interference leakage [31], can lead to alternative beamforming designs.

V-B Numerical simulations

TABLE III: Simulation parameter settings.
Symbol Description Value
ff Carrier frequency 30 GHz
λ\lambda Wave length 0.01 m
σ2\sigma^{2} Noise power -20 dBm
PP Signal power 0.5 W
εr\varepsilon_{r} Relative permittivity 2
μ\mu Permeability 4π×1074\pi\times 10^{-7}
AFA_{F} Antenna factor 1
𝒞t\mathcal{C}_{t} Transmitter region xr,κ,yr,κ,zr,κ[100λ,100λ]x_{r,\kappa},y_{r,\kappa},z_{r,\kappa}\in[-100\lambda,100\lambda]
LL #transmitting antennas 8
KK #users {1,2,4,8}\{1,2,4,8\}
Refer to caption
(a) K=1K=1
Refer to caption
(b) K=2K=2
Refer to caption
(c) K=4K=4
Refer to caption
(d) K=8K=8
Figure 9: Optimal positions and orientations of the transmitting antennas, optimal orientations of the receiving antennas, and the corresponding electric field isosurfaces in the MU-MISO system for different user counts.

This section presents numerical simulations to evaluate the performance of our PAMA framework in the context of sum-rate optimization. In the simulations, we consider users uniformly distributed within a 3D coverage area – a cube with dimensions of 200200 meters by 200200 meters by 200200 meters – centered at the origin 𝒐\bm{o} of the CCS. Unless specified otherwise, the simulation parameters used are summarized in Tab. III. The BS is equipped with L=8L=8 antennas, and the number of users KK varies among {1,2,4,8}\{1,2,4,8\}.

V-B1 Impact of polarization matching

Our PAMA framework integrates polarization effects into MAs and reveals a significant advantage of MAs: the ability to achieve optimal polarization alignment between the transmitter and receiver. To validate this capability, we consider the following antenna configurations within the PAMA framework:

  • Configuration 1: Transmitting antennas are randomly positioned and oriented; receiving antennas have random orientations.

  • Configuration 2: Transmitting antennas have optimized positions but random orientations; receiving antennas have random orientations.

  • Configuration 3: Both the positions and orientations of the transmitting antennas are optimized; receiving antennas have random orientations.

  • Configuration 4: Transmitting antennas are randomly positioned and oriented; the orientations of the receiving antennas are optimized.

  • Configuration 5: Both the positions and orientations of the transmitting antennas are optimized; the orientations of the receiving antennas are also optimized.

In Configurations 1 and 2, there is no optimization of polarization matching in the system. In Configuration 3, polarization matching is optimized by adjusting the orientations of the transmitting antennas. In Configuration 4, polarization matching is optimized by adjusting the orientations of the receiving antennas. In Configuration 5, polarization matching is optimized by adjusting the orientations of both transmitting and receiving antennas.

Refer to caption
Figure 10: The equivalent total SNR γtotal\gamma_{\text{total}} versus the number of users KK for various polarization matching configurations.

The simulation results for the configurations described above are illustrated in Fig. 10, where we evaluate the equivalent total SNR γtotal\gamma_{\text{total}} against the number of users, with the number of transmitting antennas set to L=8L=8. The following observations can be made:

  • As the number of users increases, the equivalent total system SNR decreases. With more users sharing the same transmission resources, each user receives less focused beamforming gain.

  • The performance curves of configurations 1 and 2 overlap, indicating that antenna translation does not provide any additional benefits. This occurs because the effect of antenna translation can be fully replicated through precoding, as explained in Remark 1.

  • When comparing configuration 3 to configuration 1, we find that optimizing polarization matching by rotating the transmitting antennas effectively enhances the equivalent SNR. The SNR gain increases with the number of users and averages around 33 dB. If we aim for a target SNR of 4040dB, the system can support one additional user by employing polarization matching through transmitting antenna rotation.

  • Comparing configuration 4 with configuration 1 reveals that optimizing polarization matching by rotating the receiving antennas provides minimal gain. This observation aligns with the analysis in Fig. 7, which shows that randomly oriented receiving antennas have approximately a 99%99\% probability of capturing at least half of the maximum signal energy. However, optimizing the orientation of the receiving antennas remains important for reducing outage probability, as there is still a small chance – less than 1%1\% – of complete signal loss without such optimization.

  • As expected, configuration 5 achieves the highest system gain and effectively mitigates the decline in equivalent SNR as the number of users increases. This outcome confirms the effectiveness of the PAMA framework in optimizing polarization matching between the transmitter and receiver, thereby enhancing overall system performance.

To provide a clearer understanding of the optimized PAMA in configuration 5, we present in Fig. 9 the optimal orientations of both the transmitting and receiving antennas as the number of users increases, along with the corresponding electric field isosurfaces.

  • In the single-user scenario illustrated in Fig. 9(a), all transmitting antennas are oriented to focus their maximum radiation patterns toward the target. The combined electric field distribution from the eight antennas resembles that of a single half-wavelength dipole antenna, forming a characteristic doughnut-shaped pattern.

  • As the number of users increases, the electric field distribution becomes more isotropic to effectively accommodate all users’ requirements.

Refer to caption
Figure 11: The equivalent total SNR γtotal\gamma_{\text{total}} versus the antenna rotation granularity for varying numbers of users KK.
Refer to caption
Figure 12: The equivalent total SNR γtotal\gamma_{\text{total}} versus the total transmission power PP for varying numbers of users KK.

V-B2 Impact of antenna rotation granularity

Implementing PAMA in practice faces the challenge of precisely controlling the rotations of both transmitting and receiving antennas. While the simulation results in Fig. 10 assume infinite control precision, real-world systems are constrained by finite rotation granularity. To assess how rotation granularity affects system performance, we quantize the polar and azimuthal angles of the transmitting and receiving antennas, i.e., θt,(n)\theta_{t,\ell}^{(n)}, ϕt,(n)\phi_{t,\ell}^{(n)}, θr,κ(n)\theta_{r,\kappa}^{(n)}, and ϕr,κ(n)\phi_{r,\kappa}^{(n)}. We then evaluate the decline in the equivalent total SNR, γtotal\gamma_{\text{total}}, across different quantization levels. The simulation results are depicted in Fig. 11, where a resolution level of 0 corresponds to the optimal orientation with infinite precision.

The findings show that the negative impact of reduced quantization precision on γtotal\gamma_{\text{total}} becomes more pronounced as the number of users increases. For example, with a rotation precision limited to 8080^{\circ}, a system serving K=8K=8 users experiences a 33 dB decrease in γtotal\gamma_{\text{total}}. When the rotation resolution is improved to 3030^{\circ}, the achievable γtotal\gamma_{\text{total}} stays within 0.50.5 dB of the maximum possible gain. In this scenario, each angular dimension – both the polar and azimuthal angles – requires adjustments in only 360/30=12360^{\circ}/30^{\circ}=12 discrete directions. This significantly reduces the implementation complexity, as controlling antenna rotations becomes much more practical with a finite set of orientations. These results demonstrate that PAMA can perform effectively even with finite control precision, confirming that high-resolution rotational control is not necessary to achieve near-optimal performance.

V-B3 Transmit power reduction and convergence behavior

Next, we examine how the equivalent total SNR γtotal\gamma_{\text{total}} varies with the total transmit power PP. The simulation results are shown in Fig. 12 and explained as follows.

  • For a fixed number of users, increasing PP leads to a higher γtotal\gamma_{\text{total}}, whereas increasing the number of users results in a lower γtotal\gamma_{\text{total}}.

  • With polarization matching through rotations of both transmitting and receiving antennas, the PAMA framework demonstrates greater power efficiency. To achieve the same total SNR, PAMA significantly reduces the required transmit power. This outcome aligns with expectations, as the transmitter can adjust its orientation to direct higher field intensity toward each user, while the receiver can optimize matching efficiency by rotating its antenna for optimal alignment.

Refer to caption
Figure 13: Comparison of relative gains between configuration 5 and configuration 1 with the increase in KK, and the convergence behavior of the optimization algorithm.

Finally, Fig. 13 illustrates the convergence behavior of our optimization algorithm by showing the relative gain in total SNR. The results reveal that:

  • The total SNR increases with each iteration across various values of KK. For K<5K<5, the algorithm typically converges within about 1010 iterations, achieving a gain of less than 33 dB. However, when K5K\geq 5, more iterations are required for convergence.

  • Despite the increased number of iterations needed for larger KK, the additional receiving antennas introduce greater degrees of freedom, enabling higher gains within a reasonable number of iterations. For instance, with K=8K=8, the gain can reach up to 77 dB within 20 iterations.

VI Conclusion

This paper presented the PAMA framework, which integrates polarization effects into the translation and rotation of antennas in wireless communication systems. By factoring in polarization, PAMA uncovers a defining advantage of MA: the ability to optimize polarization matching between the transmitter and receiver. This enhancement debunks the belief that MA is only effective in low-frequency multipath environments and extends the applicability of MA to higher frequency bands like the mmWave spectrum.

At mmWave frequencies, antennas are physically smaller, making them easier to move and rotate, which facilitates the practical implementation of PAMA systems. The dominance of LOS paths in the mmWave band simplifies channel estimation compared to traditional MA systems operating in environments with rich multipath. Additionally, PAMA achieves near-optimal performance with low rotation granularity. These factors collectively enhance the feasibility of deploying PAMA in real-world scenarios.

The insights from this work underscore the critical role of polarization in maximizing the performance of MA. Future research may focus on integrating PAMA with other advanced communication technologies, evaluating its performance in various high-frequency settings, and addressing practical implementation challenges such as antenna control mechanisms and mobility constraints. By advancing the understanding of how polarization can be harnessed in movable antenna systems, this work lays the foundation for developing more efficient and reliable wireless communication networks.

Appendix A Electric Field Radiation

As a preliminary, the appendix explains the derivation of (5), the radiation patterns of a half-wave dipole.

Maxwell’s equations can be written as

×𝓔=jωμ𝓗,\displaystyle\nabla\times\bm{\mathcal{E}}=-j\omega\mu\bm{\mathcal{H}}, (31)
×𝓗=jωε𝓔+𝑱,\displaystyle\nabla\times\bm{\mathcal{H}}=j\omega\varepsilon\bm{\mathcal{E}}+\bm{J}, (32)
𝓔=ρε,\displaystyle\nabla\cdot\bm{\mathcal{E}}=\frac{\rho}{\varepsilon}, (33)
𝓗=0,\displaystyle\nabla\cdot\bm{\mathcal{H}}=0, (34)

where

  • 𝓔\bm{\mathcal{E}} and 𝓗\bm{\mathcal{H}} refer to the electric field and magnetic field.

  • μ\mu is the magnetic permeability of air, describing its ability to conduct magnetic field lines. ε\varepsilon is the permittivity, characterizing a material’s capacity to store electric charge.

  • 𝑱\bm{J} and ρ\rho represent the source current density and source charge density, respectively.

  • ω\omega is the angular frequency of the wave. The expression 𝒂×𝒃\bm{a}\times\bm{b} denotes the cross product; ×𝒃\nabla\times\bm{b} denotes the curl of 𝒃\bm{b}; 𝒂𝒃\bm{a}\cdot\bm{b} denotes the dot product; 𝒃\nabla\cdot\bm{b} denotes the divergence of 𝒃\bm{b}, and \nabla is the Hamiltonian operator.

Since the divergence of 𝓗\bm{\mathcal{H}} in (34) is zero, the vector field 𝓗\bm{\mathcal{H}} has no sources or sinks; it only circulates. This allows us to define a magnetic vector potential 𝑨\bm{A} so that 𝓗\bm{\mathcal{H}} can be described as the curl of 𝑨\bm{A}:

𝓗=1μ×𝑨.\displaystyle\bm{\mathcal{H}}=\frac{1}{\mu}\nabla\times\bm{A}. (35)

Substituting (35) into (31) yields ×(𝓔+jω𝑨)=0\nabla\times(\bm{\mathcal{E}}+j\omega\bm{A})=0. In vector calculus, if the curl of a vector field is zero, that field can be expressed as the gradient of a scalar potential. Therefore, we can define an electric scalar potential Φ\Phi such that

𝓔+jω𝑨=Φ.\displaystyle\bm{\mathcal{E}}+j\omega\bm{A}=-\nabla\Phi. (36)

It is important to note that ×Φ\nabla\times\nabla\Phi is always zero for any scalar potential Φ\Phi.

Equations (35) and (36) show how 𝓗\bm{\mathcal{H}} and 𝓔\bm{\mathcal{E}} can be represented in terms of their potential functions. This means that once we find these potential functions, we can determine the corresponding field expressions. The next step is to solve for these potential functions.

First, substituting (35) into (32) gives

1μ××𝑨=jωε𝓔+𝑱.\displaystyle\frac{1}{\mu}\nabla\times\nabla\times\bm{A}=j\omega\varepsilon\bm{\mathcal{E}}+\bm{J}. (37)

Next, we use the vector identity ××𝑨(𝑨)2𝑨\nabla\times\nabla\times\bm{A}\equiv\nabla(\nabla\cdot\bm{A})-\nabla^{2}\bm{A}, where 2=\nabla^{2}=\nabla\cdot\nabla represents the Laplace operator, which calculates the divergence of the gradient. By applying (36), we can rewrite (37) as

(𝑨)2𝑨=jωμε(jω𝑨Φ)+μ𝑱.\displaystyle\nabla(\nabla\cdot\bm{A})-\nabla^{2}\bm{A}=j\omega\mu\varepsilon\left(-j\omega\bm{A}-\nabla\Phi\right)+\mu\bm{J}. (38)

After a simple rearrangement and using the Lorentz condition 𝑨=jωμεΦ\nabla\cdot\bm{A}=-j\omega\mu\varepsilon\Phi[21], we get

2𝑨+ω2με𝑨=μ𝑱.\displaystyle\nabla^{2}\bm{A}+\omega^{2}\mu\varepsilon\bm{A}=-\mu\bm{J}. (39)

This equation is known as the vector wave equation. With this differential equation, if the current 𝑱\bm{J} is known, we can solve for the vector potential 𝑨\bm{A}. Once 𝑨\bm{A} is determined, the magnetic field 𝓗\bm{\mathcal{H}} can be found using (35). Then, by using the Lorentz condition along with (36), the electric field 𝓔\bm{\mathcal{E}} can also be solved.

To solve for 𝑨\bm{A}, we break down the vector wave equation in (39) into three separate scalar equations. Substituting these into (39), we obtain

2Ax+β2Ax=μJx,\displaystyle\nabla^{2}A_{x}+\beta^{2}A_{x}=-\mu J_{x},
2Ay+β2Ay=μJy,\displaystyle\nabla^{2}A_{y}+\beta^{2}A_{y}=-\mu J_{y},
2Az+β2Az=μJz,\displaystyle\nabla^{2}A_{z}+\beta^{2}A_{z}=-\mu J_{z}, (40)

where β=ω2με=2πλ\beta=\omega^{2}\mu\varepsilon=\frac{2\pi}{\lambda} is the phase constant of the transmitted wave, representing the phase change per unit distance.

The three equations in (A) are structurally the same. Once we solve one of them, the solutions for the other two follow easily. To start, we find the solution for a point source. This solution, known as the unit impulse response, can be used to build a general solution by considering an arbitrary source as a combination of many point sources. The differential equation for a point source is

2ψ+β2ψ=δ(x)δ(y)δ(z),\displaystyle\nabla^{2}\psi+\beta^{2}\psi=-\delta(x)\delta(y)\delta(z), (41)

where ψ\psi represents the response to a point source located at the origin, and δ()\delta(\cdot) is the unit impulse function. Although a point source is infinitesimally small, the current it represents has a specific direction. In practical problems, the point source models a small segment of current with a defined direction. If we consider the point source current to be directed along the zz-axis, then ψ=Az\psi=A_{z}.

Since the point source only exists at the origin and is zero everywhere else, (41) can be simplified to

2ψ+β2ψ=0.\displaystyle\nabla^{2}\psi+\beta^{2}\psi=0. (42)

The two solutions to this equation are e±jβ|𝒑||𝒑|\frac{e^{\pm j\beta|\bm{p}|}}{|\bm{p}|}, where |𝒑||\bm{p}| is the distance from the point source to the observation point 𝒑\bm{p}. These solutions represent waves moving outward and inward from the source, respectively. The physically relevant solution is the one that describes waves propagating outward from the point source. By determining the constant of proportionality, the solution for the point source becomes

ψ=ejβ|𝒑|4π|𝒑|.\displaystyle\psi=\frac{e^{-j\beta|\bm{p}|}}{4\pi|\bm{p}|}. (43)

For a current density directed along the zz-axis, the resulting vector potential will also be zz-directed. If we think of the source as a collection of point sources with a distribution given by JzJ_{z}, the response AzA_{z} can be found by summing the responses of each point source using (43). This summation is represented by an integral over the source volume vv:

Az=vμJzejβ|𝒑|4π|𝒑|𝑑v.\displaystyle A_{z}=\iiint\limits_{v}{\mu J_{z}\frac{e^{-j\beta|\bm{p}|}}{4\pi|\bm{p}|}dv}. (44)

Similar integrals apply for the xx- and yy-components. The overall solution is the sum of all components, expressed as

𝑨=vμ𝑱ejβ|𝑹|4π|𝑹|𝑑v,\displaystyle\bm{A}=\iiint\limits_{v}{\mu\bm{J}\frac{e^{-j\beta|\bm{R}|}}{4\pi|\bm{R}|}dv}, (45)

where 𝑹\bm{R} is the vector from the location of the point source 𝒓\bm{r} to the observation point 𝒑\bm{p}, hence is given by 𝑹=𝒑𝒓\bm{R}=\bm{p}-\bm{r}. Therefore, we arrive at the solution to the vector wave equation (39) as shown in (45).

In summary, when the current density 𝑱\bm{J} is known, the steps to find the induced field are:

  1. 1.

    Use (45) to find the magnetic vector potential 𝑨\bm{A};

  2. 2.

    Determine the magnetic field 𝓗\bm{\mathcal{H}} using (35);

  3. 3.

    Apply (32) and consider that 𝑱=0\bm{J}=0 at the observation point 𝒑\bm{p} to derive the electric field as 𝓔=1jωε×𝓗\bm{\mathcal{E}}=\frac{1}{j\omega\varepsilon}\nabla\times\bm{\mathcal{H}}.

For a half-wave dipole located at the origin and placed along the zz-axis, the current density 𝑱\bm{J} distribution is given by

J(z)=sin[β(λ4|z|)],\displaystyle J(z)=\sin{\left[\beta\left(\frac{\lambda}{4}-|z|\right)\right]}, (46)

where |z|λ4|z|\leq\frac{\lambda}{4}.

The magnetic vector potential 𝑨\bm{A} can be derived from (45) as follows:

𝑨=vμ𝑱ejβ|𝑹|4π|𝑹|𝑑v\displaystyle\bm{A}=\iiint\limits_{v}{\mu\bm{J}\frac{e^{-j\beta|\bm{R}|}}{4\pi|\bm{R}|}dv} (47)
=(a)μejβ|𝒑|4π|𝒑|v𝑱ejβ𝒑|𝒑|𝒓𝑑v\displaystyle\overset{(a)}{=}\mu\frac{e^{-j\beta|\bm{p}|}}{4\pi|\bm{p}|}\iiint\limits_{v}{\bm{J}e^{-j\beta\frac{\bm{p}}{|\bm{p}|}\cdot\bm{r}}dv}
=(b)𝒌μejβ|𝒑|4π|𝒑|J(z)ejβzcosθ𝑑z\displaystyle\overset{(b)}{=}\bm{k}\mu\frac{e^{-j\beta|\bm{p}|}}{4\pi|\bm{p}|}\int{J(z)e^{-j\beta z\cos{\theta}}dz}
=(c)𝒌μejβ|𝒑|4π|𝒑|sin[β(λ4|z|)]ejβzcosθ𝑑z\displaystyle\overset{(c)}{=}\bm{k}\mu\frac{e^{-j\beta|\bm{p}|}}{4\pi|\bm{p}|}\int{\sin{\left[\beta\left(\frac{\lambda}{4}-|z|\right)\right]}e^{-j\beta z\cos{\theta}}dz}
=(d)2𝒌μejβ|𝒑|4π|𝒑|jcos(π2cosθ)+sin(π2cosθ)cosθjβsin2θ,\displaystyle\overset{(d)}{=}2\bm{k}\mu\frac{e^{-j\beta|\bm{p}|}}{4\pi|\bm{p}|}\frac{j\cos{\left(\frac{\pi}{2}\cos{\theta}\right)}+\sin{\left(\frac{\pi}{2}\cos{\theta}\right)}-\cos{\theta}}{j\beta\sin^{2}{\theta}},

where (a) follows from |𝑹|=|𝒑|𝒓𝒑|𝒑||\bm{R}|=|\bm{p}|-\bm{r}\cdot\frac{\bm{p}}{|\bm{p}|} in the far field; (b) follows because the polar angle of the observation point 𝒑\bm{p} is the emission angel of the transmitting antenna when the dipole is oriented along the zz-axis; (c) follows from (46); (d) follows from β=2πλ\beta=\frac{2\pi}{\lambda}, cosx=ejx+ejx2\cos{x}=\frac{e^{jx}+e^{-jx}}{2}, and ejπ2+ejπ2=0e^{j\frac{\pi}{2}}+e^{-j\frac{\pi}{2}}=0, ejπ2ejπ2=2je^{j\frac{\pi}{2}}-e^{-j\frac{\pi}{2}}=2j.

Finally, the electric field induced at 𝒑\bm{p} can be derived as

𝓔=jωRe{Aθ}ϑ=jωsinθRe{Az}ϑ\displaystyle\bm{\mathcal{E}}=-j\omega\text{Re}\{A_{\theta}\}\bm{\vartheta}=j\omega\sin{\theta}\text{Re}\{A_{z}\}\bm{\vartheta} (48)
=jωsinθ2μejβ|𝒑|4π|𝒑|cos(π2cosθ)βsin2θϑ\displaystyle=j\omega\sin{\theta}2\mu\frac{e^{-j\beta|\bm{p}|}}{4\pi|\bm{p}|}\frac{\cos{\left(\frac{\pi}{2}\cos{\theta}\right)}}{\beta\sin^{2}{\theta}}\bm{\vartheta}
=jωμ2βejβ|𝒑|4π|𝒑|cos(π2cosθ)sinθϑ.\displaystyle=j\omega\mu\frac{2}{\beta}\frac{e^{-j\beta|\bm{p}|}}{4\pi|\bm{p}|}\frac{\cos{\left(\frac{\pi}{2}\cos{\theta}\right)}}{\sin{\theta}}\bm{\vartheta}.

Appendix B Proof of Corollary 3

Since the first multiplicative factor in (9) is a constant, the behavior of the magnitude is determined by the second term: Acos(π2cosθ,𝒑(e))sinθ,𝒑(e)A\triangleq\frac{\cos{\left(\frac{\pi}{2}\cos{\theta^{(e)}_{\ell,\bm{p}}}\right)}}{\sin{\theta^{(e)}_{\ell,\bm{p}}}}. We analyze the derivative of AA with respect to (w.r.t) θ,𝒑(e)\theta^{(e)}_{\ell,\bm{p}} to examine its variation.

The derivative can be expressed as

Aθ,𝒑(e)=1sin2θ,𝒑(e)[(π412)cos(π2cosθ,𝒑(e)+θ,𝒑(e))\displaystyle\frac{\partial A}{\partial\theta^{(e)}_{\ell,\bm{p}}}=\frac{1}{\sin^{2}{\theta^{(e)}_{\ell,\bm{p}}}}\Bigg{[}\left(-\frac{\pi}{4}-\frac{1}{2}\right)\cos{\left(\frac{\pi}{2}\cos{\theta^{(e)}_{\ell,\bm{p}}}+\theta^{(e)}_{\ell,\bm{p}}\right)}
+(π412)cos(π2cosθ,𝒑(e)θ,𝒑(e))],\displaystyle\hskip 28.45274pt+\left(\frac{\pi}{4}-\frac{1}{2}\right)\cos{\left(\frac{\pi}{2}\cos{\theta^{(e)}_{\ell,\bm{p}}}-\theta^{(e)}_{\ell,\bm{p}}\right)}\Bigg{]}, (49)

where θ,𝒑(e)[0,π]\theta^{(e)}_{\ell,\bm{p}}\in[0,\pi]. To simplify the notations, we define A1π2cosθ,𝒑(e)+θ,𝒑(e)A_{1}\triangleq\frac{\pi}{2}\cos{\theta^{(e)}_{\ell,\bm{p}}}+\theta^{(e)}_{\ell,\bm{p}} and A2π2cosθ,𝒑(e)θ,𝒑(e)A_{2}\triangleq\frac{\pi}{2}\cos{\theta^{(e)}_{\ell,\bm{p}}}-\theta^{(e)}_{\ell,\bm{p}}.

When θ,𝒑(e)(0,π2)\theta^{(e)}_{\ell,\bm{p}}\in(0,\frac{\pi}{2}), we have

A1θ,𝒑(e)=π2sinθ,𝒑(e)+1.\frac{\partial A_{1}}{\partial\theta^{(e)}_{\ell,\bm{p}}}=-\frac{\pi}{2}\sin{\theta^{(e)}_{\ell,\bm{p}}}+1.

Thus, A1θ,𝒑(e)\frac{\partial A_{1}}{\partial\theta^{(e)}_{\ell,\bm{p}}} is a monotonically decreasing function of θ,𝒑(e)\theta^{(e)}_{\ell,\bm{p}} from 1 to π2+1<0-\frac{\pi}{2}+1<0. The minimum value of A1=π2A_{1}=\frac{\pi}{2} occurs at θ,𝒑(e)=0\theta^{(e)}_{\ell,\bm{p}}=0 or θ,𝒑(e)=π2\theta^{(e)}_{\ell,\bm{p}}=\frac{\pi}{2}. Since

A1π2+θ,𝒑(e)<π,A_{1}\leq\frac{\pi}{2}+\theta^{(e)}_{\ell,\bm{p}}<\pi,

we have π2<A1<π\frac{\pi}{2}<A_{1}<\pi for θ,𝒑(e)(0,π2)\theta^{(e)}_{\ell,\bm{p}}\in(0,\frac{\pi}{2}). Therefore, cos(A1)<0\cos{\left(A_{1}\right)}<0, and

(π412)cos(A1)>0.\left(-\frac{\pi}{4}-\frac{1}{2}\right)\cos{\left(A_{1}\right)}>0.

On the other hand, we have

A2θ,𝒑(e)=π2sinθ,𝒑(e)1<0.\frac{\partial A_{2}}{\partial\theta^{(e)}_{\ell,\bm{p}}}=-\frac{\pi}{2}\sin{\theta^{(e)}_{\ell,\bm{p}}}-1<0.

A2A_{2} is a monotonically decreasing function when θ,𝒑(e)(0,π2)\theta^{(e)}_{\ell,\bm{p}}\in(0,\frac{\pi}{2}), with A2(0)=π2A_{2}(0)=\frac{\pi}{2} and A2(π2)=π2A_{2}(\frac{\pi}{2})=-\frac{\pi}{2}. Therefore, |A2|<π2|A_{2}|<\frac{\pi}{2}, implying that cos(π2cosθθ)>0\cos\left(\frac{\pi}{2}\cos\theta-\theta\right)>0, and thus

(π412)cos(A2)>0.\left(\frac{\pi}{4}-\frac{1}{2}\right)\cos{\left(A_{2}\right)}>0.

Additionally, sin2θ,𝒑(e)>0\sin^{2}{\theta^{(e)}_{\ell,\bm{p}}}>0 when θ,𝒑(e)(0,π2)\theta^{(e)}_{\ell,\bm{p}}\in(0,\frac{\pi}{2}). Overall, we have Aθ,𝒑(e)>0\frac{\partial A}{\partial\theta^{(e)}_{\ell,\bm{p}}}>0 when θ,𝒑(e)(0,π2)\theta^{(e)}_{\ell,\bm{p}}\in(0,\frac{\pi}{2}).

Similarly, it can be shown that Aθ,𝒑(e)<0\frac{\partial A}{\partial\theta^{(e)}_{\ell,\bm{p}}}<0 when θ(π2,π)\theta\in(\frac{\pi}{2},\pi). Since Aθ,𝒑(e)|θ=π2=0\frac{\partial A}{\partial\theta^{(e)}_{\ell,\bm{p}}}\big{|}_{\theta=\frac{\pi}{2}}=0, both AA and the magnitude |,𝒑||\mathcal{E}_{\ell,\bm{p}}| reach their maximum when θ,𝒑(e)=π2\theta^{(e)}_{\ell,\bm{p}}=\frac{\pi}{2}, proving the corollary.

References

  • [1] S. Dang, O. Amin, B. Shihada, and M.-S. Alouini, “What should 6G be?” Nature Electronics, vol. 3, no. 1, pp. 20–29, 2020.
  • [2] C.-X. Wang, X. You, X. Gao, X. Zhu, Z. Li, C. Zhang, H. Wang, Y. Huang, Y. Chen, H. Haas et al., “On the road to 6G: Visions, requirements, key technologies, and testbeds,” IEEE Communications Surveys & Tutorials, vol. 25, no. 2, pp. 905–974, 2023.
  • [3] L. Zhu, W. Ma, and R. Zhang, “Modeling and performance analysis for movable antenna enabled wireless communications,” IEEE Transactions on Wireless Communications, 2023.
  • [4] X. Shao, Q. Jiang, and R. Zhang, “6D movable antenna based on user distribution: Modeling and optimization,” arXiv preprint arXiv:2403.08123, 2024.
  • [5] X. Shi, X. Shao, and R. Zhang, “Capacity maximization for base station with hybrid fixed and movable antennas,” arXiv preprint arXiv:2405.07176, 2024.
  • [6] B. Ning, S. Yang, Y. Wu, P. Wang, W. Mei, C. Yuen, and E. Björnson, “Movable antenna-enhanced wireless communications: General architectures and implementation methods,” arXiv preprint arXiv:2407.15448, 2024.
  • [7] K.-K. Wong, A. Shojaeifard, K.-F. Tong, and Y. Zhang, “Fluid antenna systems,” IEEE Transactions on Wireless Communications, vol. 20, no. 3, pp. 1950–1962, 2020.
  • [8] K.-K. Wong and K.-F. Tong, “Fluid antenna multiple access,” IEEE Transactions on Wireless Communications, vol. 21, no. 7, pp. 4801–4815, 2021.
  • [9] K.-K. Wong, K.-F. Tong, Y. Shen, Y. Chen, and Y. Zhang, “Bruce Lee-inspired fluid antenna system: Six research topics and the potentials for 6G,” Frontiers in Communications and Networks, vol. 3, pp. 853–416, 2022.
  • [10] K. K. Wong, A. Shojaeifard, K.-F. Tong, and Y. Zhang, “Performance limits of fluid antenna systems,” IEEE Communications Letters, vol. 24, no. 11, pp. 2469–2472, 2020.
  • [11] M. Khammassi, A. Kammoun, and M.-S. Alouini, “A new analytical approximation of the fluid antenna system channel,” IEEE Transactions on Wireless Communications, vol. 22, no. 12, pp. 8843–8858, 2023.
  • [12] T. Yoo and A. Goldsmith, “On the optimality of multiantenna broadcast scheduling using zero-forcing beamforming,” IEEE Journal on Selected Areas in Communications, vol. 24, no. 3, pp. 528–541, 2006.
  • [13] D. Tse and P. Viswanath, Fundamentals of wireless communication.   Cambridge university press, 2005.
  • [14] W. Tang, M. Z. Chen, X. Chen, J. Y. Dai, Y. Han, M. Di Renzo, Y. Zeng, S. Jin, Q. Cheng, and T. J. Cui, “Wireless communications with reconfigurable intelligent surface: Path loss modeling and experimental measurement,” IEEE Transactions on Wireless Communications, vol. 20, no. 1, pp. 421–439, 2020.
  • [15] Q. Wu and R. Zhang, “Intelligent reflecting surface enhanced wireless network via joint active and passive beamforming,” IEEE Transactions on Wireless Communications, vol. 18, no. 11, pp. 5394–5409, 2019.
  • [16] L. Zhu, W. Ma, B. Ning, and R. Zhang, “Movable-antenna enhanced multiuser communication via antenna position optimization,” IEEE Transactions on Wireless Communications, 2023.
  • [17] W. Ma, L. Zhu, and R. Zhang, “MIMO capacity characterization for movable antenna systems,” IEEE Transactions on Wireless Communications, 2023.
  • [18] C. Psomas, P. J. Smith, H. A. Suraweera, and I. Krikidis, “Continuous fluid antenna systems: Modeling and analysis,” IEEE Communications Letters, 2023.
  • [19] W. L. Stutzman and G. A. Thiele, Antenna theory and design.   John Wiley & Sons, 2012.
  • [20] M. Costa, A. Richter, and V. Koivunen, “DoA and polarization estimation for arbitrary array configurations,” IEEE Transactions on Signal Processing, vol. 60, no. 5, pp. 2330–2343, 2012.
  • [21] Y. Huang, Antennas: from theory to practice.   John Wiley & Sons, 2021.
  • [22] R. U. Nabar, H. Bolcskei, V. Erceg, D. Gesbert, and A. J. Paulraj, “Performance of multiantenna signaling techniques in the presence of polarization diversity,” IEEE Transactions on Signal Processing, vol. 50, no. 10, pp. 2553–2562, 2002.
  • [23] M. Kurum, M. Deshpande, A. T. Joseph, P. E. O’Neill, R. H. Lang, and O. Eroglu, “SCoBi-Veg: A generalized bistatic scattering model of reflectometry from vegetation for signals of opportunity applications,” IEEE Transactions on Geoscience and Remote Sensing, vol. 57, no. 2, pp. 1049–1068, 2018.
  • [24] Y. Niu, Y. Li, D. Jin, L. Su, and A. V. Vasilakos, “A survey of millimeter wave communications (mmWave) for 5G: opportunities and challenges,” Wireless networks, vol. 21, pp. 2657–2676, 2015.
  • [25] R. W. Heath, N. Gonzalez-Prelcic, S. Rangan, W. Roh, and A. M. Sayeed, “An overview of signal processing techniques for millimeter wave MIMO systems,” IEEE Journal of Selected Topics in Signal Processing, vol. 10, no. 3, pp. 436–453, 2016.
  • [26] Study on channel model for frequencies from 0.5 to 100 GHz (3GPP TR 38.901 version 18.0.0 Release 18), 3GPP Std. 3GPP TR 38.901, 2024.
  • [27] J. A. Stratton, Electromagnetic theory.   John Wiley & Sons, 2007, vol. 33.
  • [28] C. Parazzoli, R. Greegor, K. Li, B. Koltenbah, and M. Tanielian, “Experimental verification and simulation of negative index of refraction using snell’s law,” Physical Review Letters, vol. 90, no. 10, p. 107401, 2003.
  • [29] O. Tervo, L.-N. Tran, and M. Juntti, “Optimal energy-efficient transmit beamforming for multi-user MISO downlink,” IEEE Transactions on Signal Processing, vol. 63, no. 20, pp. 5574–5588, 2015.
  • [30] W. Yu, W. Rhee, S. Boyd, and J. M. Cioffi, “Iterative water-filling for gaussian vector multiple-access channels,” IEEE Transactions on Information Theory, vol. 50, no. 1, pp. 145–152, 2004.
  • [31] M. Sadek, A. Tarighat, and A. H. Sayed, “A leakage-based precoding scheme for downlink multi-user MIMO channels,” IEEE Transactions on Wireless Communications, vol. 6, no. 5, pp. 1711–1721, 2007.