Physical Modelling and Cancellation of External Passive Intermodulation in FDD MIMO ^††thanks: The research was carried out at Skoltech and supported by the Russian Science Foundation (project no. 24-29-00189).

Despite the wide variety of PIM compensation methods [4, 16, 17, 18, 19], modern approaches do not allow one to simulate external PIM in the MIMO system physically. All known papers use a compensation model approach. In this case, the compensation model design follows the physical mechanism of external PIM generation with significant simplifications. Therefore, these research results are limited by real data measurements. Such kinds of measurement are not available for many researchers. Also, none of these methods provides a comprehensive process by which artificial interference caused by an external PIM source could be simulated directly, especially in an arbitrary MIMO system. The authors of [20, 11, 21, 6, 22] also mention that they are unaware of any works explicitly aiming to model and cancel air-induced PIM in FDD MIMO scenarios.

It is also worth noting that in almost all recent articles, external PIM compensation methods intended for large-scale MIMO systems are tested only on base stations with 2 transceiver paths. Thus, works [6, 22] limit their research by a dual TX/RX chain MIMO system and 3 external sources of PIM. This is not enough for a comprehensive MIMO system test. Some works [17] consider more complex systems with more transceiver chains. Unfortunately, the results of such works are difficult to reproduce since no data and simulation codes are available in open sources.

Consequently, developing a unified physical model that allows a realistic simulation of the external PIM phenomenon in FDD MIMO systems is an important direction from both theoretical and practical points of view.

All the external PIM cancellation algorithms mentioned in the literature require real hardware-measured testing data. Accessing the real BS hardware or measured data is impossible in most cases. Additionally, setting up an experiment and environment with an artificial PIM source is a rather expensive procedure requiring many specialists and additional preparations.

This paper presents a new physical model of artificial external PIM generation based on electromagnetic theory. The model has the following features:

1) The ability to generate any number of external point PIM sources based on arbitrary Uniform Rectangular Array (URA) structure;

2) The ability of near- or far-field zone effects and polarization effects to be considered due to the near-field dipole antenna model;

3) The ability to test both UL and DL PIM compensation methods;

4) A relatively simple model that does not require large computational resources and provides acceptable modelling accuracy;

This work is based on a finite-length dipole antenna model of a zero-radius wire [24]. This allows for avoiding significant computational resource usage and provides good numerical agreement with a real antenna. Effects associated with the influence of magnetic fields are assumed to be negligible and are not considered in this article. For a dipole antenna oriented along the $z$-axis, the electric field at each frequency is given in cylindrical coordinates (1). A sinusoidal current distribution is assumed.

where $x,y,z$ are the observation point coordinates in the cartesian coordinate system, $E_{\rho},E_{\phi},E_{z}$ are electric field components in the cylindrical coordinate system, $r$ is observation point distance, $\eta$ is free space impedance, $I$ is the complex amplitude of the current feeding the antenna and $l$ is antenna length. The average current feeding $n$-th antenna can be calculated using the Ohm law:

where $P_{n}$ is antenna element radiating power, and $Z_{n}$ is antenna impedance. All antenna impedances are assumed to be the same, and antennas are considered not to affect each other. All antennas have the same radiating power.

Thus, the instantaneous value, if the input signal can be calculated as follows:

It is worth noting that the normalization to maintain power ratios is done as $\mathbb{E}[|\bm{U}_{n}(t)|^{2}]=1$ (for each TX chain separately).

The following approach with coordinate transformation has been applied to model a dipole oriented along other axes (crossed dipoles in the x-y plane). $x,y,z$ are observation point coordinates associated with $\vec{\mathbf{r}}$ vector, $x_{0}^{i},y_{0}^{i},z_{0}^{i}$ are $n$-th antenna coordinates from the array associated with the $\vec{\mathbf{r}}_{n}$ vector. The initial dipole antenna is aligned along the $z$-axis. The electric field calculated in cylindrical coordinate system at the observation point $\vec{\mathbf{r}}_{n}(x-x_{0}^{i},y-y_{0}^{i},z-z_{0}^{i})$ is $\vec{\mathbf{E}}(\vec{\mathbf{r}}-\vec{\mathbf{r}}_{n})^{cyl}=\vec{\mathbf{E}}(\Delta\vec{\mathbf{r}}_{n})^{cyl}$. Suppose that the dipole has been rotated in a local coordinate system under a linear transformation effect: $\mathbf{\Phi}_{coord}=\mathbf{\Phi}_{coord}(\Delta\vec{\mathbf{r}}_{n})$ without changes in the observation point. In this case, the observation point field will be calculated as $\vec{\mathbf{E}}(\Delta\vec{\mathbf{r}}_{n}^{*})^{cyl}=\vec{\mathbf{E}}(\mathbf{\Phi}_{coord}(\Delta\vec{\mathbf{r}}_{n}))^{cyl}$. At the same time, this field can be transformed into Cartesian coordinate system using linear transformation for vectors $\mathbf{\Phi}_{vect}=\mathbf{\Phi}_{vect}(\vec{\mathbf{r}}_{n}^{*})$: $\vec{\mathbf{E}}(\vec{\mathbf{r}}_{n}^{*})^{cart}=\mathbf{\Phi}_{vect}\vec{\mathbf{E}}(\Delta\vec{\mathbf{r}}_{n}^{*})^{cyl}$. Inverse coordinate transformation can be used to obtain the field in the original coordinate system: $\vec{\mathbf{E}}(\vec{\mathbf{r}}_{n})^{cart}=\mathbf{\Phi}_{coord}^{-1}\vec{\mathbf{E}}(\Delta\vec{\mathbf{r}}_{n}^{*})^{cart}$.

A typical MIMO antenna array consists of rows and columns of dual-polarized antenna elements radiating from the same point (crossed dipoles to consider polarisation properties). The whole array can be divided into subarrays. Each subarray is connected to two radio chains, normally one per polarization. Antennas within the subarray utilize tunable phase shifts to provide electric antenna directivity pattern tilt.

In practice, even if the antenna array consists of identical antenna elements, their current distributions may differ from the law of current distribution in a separate antenna due to the mutual coupling effect. In the current work, the mutual coupling effect is assumed to be negligible. Assume the current density distribution on the $n$-th antenna element as

where $I_{n}$ is the complex amplitude of the current applied to the $n$-th antenna element, $\vec{\bm{J}}$ is current density distribution in a local coordinate system of the antenna, $\vec{\bm{r}}$ and $\vec{\bm{r_{n}}}$ are observation point radius vectors and $n$-th antenna position in a global coordinate system, respectively. This current density distribution produces an electric field $\vec{\bm{E}}_{n}^{f}[\vec{\bm{J}_{n}}]$ at the observation point at frequency $f$, which can be calculated using (1).

Thus, the total electric field $\vec{\bm{E}^{f}}$ produced by $N_{TX}$ antenna elements at the observation point $\vec{\bm{r}}$ can be represented as a superposition of the fields of individual antenna elements, taking into account the exciting currents:

Assume that the electric field at the PIM source location (observation point) is $\vec{\bm{E}}^{f}(\vec{\bm{r}})$ is calculated using (5). The total field exciting the PIM source can be represented as a product of PIM source orientation vector $\vec{\bm{p}}=\vec{\bm{p}}(p_{x},p_{y},p_{z}),|\vec{\bm{p}}|=1$ and electric field vector at the given point:

The induced voltage at the observation point is proportional to the electric field:

The voltage after nonlinear distortion is proportional to a non-linear function of $U_{eff}^{t}=\operatorname{IDFT}(\sum_{f}U_{eff}^{f})$ taken in the time domain:

where $\operatorname{DFT}$ and $\operatorname{IDFT}$ are direct and inverse Discrete Fourier Transforms respectively. The radiation field in the reverse direction can be represented similarly to the forward path. In this case, the signal at $n$-th receiving antenna terminal will be defined as:

Assume the level of PIM generated by a point source is specified by normalization relative to the level of thermal noise. This does not affect the amplitude and phase relationship between the generating signals and the generated PIM. This is because, physically, a point source emits a tiny amount of energy, which cannot be considered using the point source model.

Nonlinear elements, depending on their nature, may also have memory effects [9]. In this case, a nonlinear element producing IMD harmonics can be modelled in various ways. One common approach is Generalized Memory Polynomial (GMP) [25]. Such nonlinear function $\mathcal{F}_{NL}=\mathcal{F}_{NL}[\mathcal{M},\mathcal{K},\mathcal{P}]$ has the following form of (10).

where $\mathcal{M}$ is the number of global precursor ($\mathcal{M}<0$) and postcursor ($\mathcal{M}>0$) taps, respectively, $\mathcal{K}$ - the number of envelope lead ($\mathcal{K}<0$) and lag ($\mathcal{K}>0$) samples; $\mathcal{P}$ is polynomial order, $g$ is amplitude coefficient, $u$ is input signal, $t$ is time index.

The modelling was carried out for 5G NR band n3 shown in Fig. 1. The component CCs of 5MHz bandwidths are located at the 1819 MHz and 1866.5 MHz frequencies respectively and their intermodulation products lie at the 1771.5 MHz central frequency (-71.25 MHz w.r.t. TX zero frequency). In baseband digital processing, where the PIM cancellation occurs, the TX and RX signals are oversampled by a factor of 16 for the BF generation. Artificially generated data is divided into 2 sets: a training set of 131072 samples and a testing set of 65536 samples. Signal simulation parameters are summarized in Table I.

In this work, the nonlinearity for PIM generation follows the structure of (10) with the following parameters for simplicity: $\mathcal{P}=\{3\},\mathcal{M}=\{0\},\mathcal{K}=\{0\}$. Thus, a simple, memoryless, nonlinear model of IMD-3 is assumed. The nonlinear model of the compensation algorithm follows the same structure as well. However, more complex parameters are chosen under the assumption of unknown nonlinearity behaviour at the external PIM source: $\mathcal{P}=\{3,5\},\mathcal{M}=\{-2,-1,0,1,2\},\mathcal{K}=\{-1,0,1\}$. The orientation of the PIM source $\vec{\bm{p}}$ is random for each scenario. The total PIM level is artificially scaled to a reasonable value.

Each antenna in the simulations is assumed to be ideally matched to the transceiver path, has a 50 Ohm impedance, and has an average radiating power of 37 dBm. The conversion of energy into radiation occurs completely.

Two scenarios of external PIM source location are considered:

1) Scenario 1: A single external PIM source is located in front of the centre of the antenna array at a distance of 2.5 m relative to it. PIM cartesian coordinates $\vec{\bm{r}}$ are: {0, 0, 2.5}.

2) Scenario 2: Three external PIM sources are located in front of the antenna with a 1 m separation along the $x$-axis and a 0.5 m separation along the $y$-axis. PIM cartesian coordinates $\vec{\bm{r}}$ are: {-1, -0.5, 3}, {0, 0, 3}, {1, 0.5, 3}.

The result for a large-scale MIMO system with 16 transmitting and receiving antennas is shown. Each transceiver chain in the system is represented by a set of half-wavelength dipoles grouped into unit cells of 4 antennas. The signal from each channel is divided equally among all antennas within a unit cell. The distance between the antennas is equal to half of the wavelength (w.r.t. the central frequency of the TX signal band). Even channels have vertical polarization; odd channels have horizontal polarization (crossed dipoles). The antenna layout is shown in Fig. 3 where numbering means antenna number which is not related to the number of TX chains directly. A dashed line limits unit cells.

The location of the PIM source significantly affects the received power distribution between the channels. This is especially noticeable in the near-field zone of the antenna array since the wave here is spherical and the electric field varies a lot. As the source moves away from the antenna array, the wave becomes flat, and such a power variation is no longer noticeable (Fig. 7). Polarization effects can also be noticed since the signal power at even antennas differs from the signal of odd antennas.