# Directional amplitude backscatter modulation with suppressed Doppler based on rotating resonant loop

Ashkan Azarfar, Nicolas Barbot, Etienne Perret

► To cite this version:

Ashkan Azarfar, Nicolas Barbot, Etienne Perret. Directional amplitude backscatter modulation with suppressed Doppler based on rotating resonant loop. Scientific Reports, 2022, 12 (1), 10.1038/s41598-022-26609-w. hal-04775552

HAL Id: hal-04775552

https://hal.science/hal-04775552v1

Submitted on 10 Nov 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

# Directional amplitude backscatter modulation with suppressed Doppler based on rotating resonant loop

Ashkan Azarfar✉, Nicolas Barbot & Etienne Perret

The directional amplitude backscatter modulation with suppressed Doppler is demonstrated based on the scattering from a symmetrically rotating resonant loop. The concept is studied theoretically and experimentally with perfectly compatible results. The symmetrical rotation of the scatterer and the effect of radial resonance, as the two crucial points to realize the idea, are highlighted through the comparison between the symmetric and non-symmetric cases, and the results obtained for scatterers with and without radial resonance. The presented backscattering modulation technique provides an amplitude modulating waveform which is uniquely linked to the directional reradiation pattern of the rotating loop scatterer in a definite resonant mode. With the pure directional amplitude modulation (DAM) induced on the backscattered wave, the envelope waveform can be accurately retrieved from the received signal using the In-phase and Quadrature (IQ) representation. The contribution of the background in a real environment can be detected and removed to obtain the exact modulating waveform. This property of the proposed backscattering modulation method can be applied for sensing, localization, and identification purposes with high sensitivity, read range, and robustness.

Radio Frequency (RF) backscatter communication based on modulated reflectors has been firstly introduced in 1948¹. Stockman¹ has considered the concept with a quite general point of view in his work, and he has demonstrated that a modulated reflector can be realized based on any phenomenon which changes the scatterer properties so that, the backscattered wave is modified in terms of magnitude, phase, or polarization during the time. For a time-harmonic incident wave, any modification on these three backscattered wave parameters will be translated to an amplitude modulation, phase modulation, or a mixed amplitude-phase modulation on the received signal. From a structural view, the addressed modulation methods in Ref.¹ can be classified to two main types which are applied for non-loaded scatterers (i.e. any scattering structure in which there is no loading point or port)² and loaded scatterers³. The latter, called "variable-damping modulation" in¹, is mainly realized by load-modulated antennas⁴, and it has been widely employed in classical UHF Radio Frequency IDentification (RFID)⁵ and Ambient Backscattering Communication (AmBC)⁶ where a modulating electronic component modifies the magnitude and phase⁷, or polarization⁸,⁹ of the reflection from the antenna during time. Usually, in this type, the energy required for the modulation is harvested from the incident electromagnetic wave, or the built-in battery, to power up the electronic component. On the other hand, in the former type based on moving non-loaded scatterers, which is more compatible with chipless RFID technology¹⁰, the required energy for modulation is provided by the motion source¹,¹¹.

In contrast to the loaded scatterer approach in which the data is stored in the electronic chip, for the non-loaded scatterer method the data can be associated either to the physical properties of the movement source (e.g. frequency and amplitude of the sound wave or periodic motion) or to the electromagnetic scattering properties of the scatterer (e.g. resonance frequency and reradiation pattern of the scatterer defined by its geometry). However, there are exceptional works like¹²,¹³ which cannot be fitted in this classification, since a load-modulated antenna has been used in Ref.¹²,¹³ while the modulation source is motion.

Following the first test examples demonstrated in Ref.¹ with rotationally modulated scatterers (i.e. corner reflectors), rotational motion has been considered in wide variety of research works¹⁴⁻²⁶ to utilize rotation induced backscattering modulation for sensing, classification, and identification goal. Most of the time, a rotating scatterer modifies a combination of the phase, magnitude and polarization on the backscattered wave, which

Grenoble INP, LCIS, University of Grenoble Alpes, 26000 Valence, France. ✉email: ashkan.azarfar@lcis.grenoble-inp.fr
yields a mixed amplitude-phase modulation on the received signal. In particular, the phase modulation caused by the radial velocity of the rotating scatterer towards the wave incidence direction (Doppler modulation), which is usually inevitable, has been extensively studied as micro-Doppler effect²⁷. However, in most researches based on micro-Doppler effect¹⁵,¹⁹⁻²³, the contribution of the phase modulation has not been clearly distinguished from the variation of the other wave parameters like magnitude and polarization, while they may present, and they have been considered jointly with the phase (micro-Doppler) modulation. In fact, since the dimension of the rotating scatterers studied in Refs.¹⁹⁻²³ (respectively helicopter blades, wind turbine, corner reflector, star-shaped wheel, and Crookes radiometer) are much larger than the operating wavelength, they cannot be modeled as a point scatterer, and therefore they tends to induce a mixed amplitude-phase modulation on the received signal, even in a symmetrically rotating and non-depolarizing configuration. It should be mentioned that, this concept has been cleverly addressed far earlier by Ref.¹⁷ for scattering from jet engines in terms of received signal analysis, however, not in terms of scatterers behaviour. Taking into account the scatterer behaviour, the scatterer can be specially designed and configured to modulate only a specific wave parameter during the rotation to provide a pure amplitude or phase modulation. For example, rotating resonant dipole scatterers in Ref.²⁴ have been configured such that the backscattered wave is modulated only in phase based on the micro-Doppler effect, which provides a pure analog phase modulation employed for RFID application. Moreover, to just modify the polarization of the reflected wave, dipole scatterers rotating perpendicularly respect to the incidence direction, which yields zero phase (Doppler) modulation, have been used in Ref.²⁶ to achieve a pure amplitude modulation for sensing and identification.

According to the authors' best knowledge, there is a lack of study to consider some special rotating scatterers which modulate the backscattered field only in magnitude, without affecting the phase¹⁵,¹⁹⁻²⁴ or the polarization²⁵,²⁶ of the wave. "Directional scatterers" with definite reradiation pattern can be used to modify only the magnitude of the backscattered wave during the rotation, provided that the Doppler modulation is suppressed somehow. In other words, such kind of scatterer modulates the magnitude of the backscattered field during rotation with a waveform linked to the reradiation pattern of the scatterer. Worth mentioning that, although Stockman has addressed the DAM induced by "angular displacement" in his paper, neither him nor afterwards research works have discussed this concept analytically with a feasible example.

In this paper, we demonstrate that resonant scatterers, specifically the rectangular loop here, can provide a pure DAM when they are symmetrically rotated in a non-depolarizing configuration. As the key point, at resonant mode of the structure, the phase of the induced current on the scatterer is imposed by the mode, and the phase variation due to Doppler will be suppressed by the resonance effect, while the different points on the scatterer have a radial velocity towards the incidence direction. Moreover, the reradiation pattern of the scatterer is determined by the resonant current, and therefore, remains invariant during the rotation. This concept leads to a DAM on the backscattered wave, while the modulation waveform is uniquely defined by the reradiation pattern of the scatterer at the resonant mode.

# Results

**Directional amplitude backscatter modulation.** For a given scatterer impinged by a time-harmonic plane wave with defined polarization and incidence direction, the spatial distribution of the magnitude and phase of the scattered field components (the field components defines the polarization state of the scattered field) perfectly describes the modulation induced on the wave reflected from the rotating scatterer. Accordingly, for a definite rotation axis, to equivalently observe the variations during the rotation of the scatterer, the magnitude and phase variation of the scattered field components can be considered for different incidence angles corresponding to different rotation angles during the time.

The considered configuration in this paper is presented by Fig. 1 where a z-polarized plane wave impinges on the scatterers at different incidence angles $\phi = \phi_i$ (equivalent to rotation of the scatterers around z-axis), and the variation of the scattered field in terms of polarization, magnitude, and phase is investigated on xy-plane for observation angles $\phi = \phi_o$ ($0^\circ < \phi_o < 360^\circ$). To preserve the polarization of the reflected wave during the rotation same as that of the incident wave, the scatterers should be designed in a non-depolarizing configuration such that the current path on them is fixed for all incidence angles, and the reradiated field has no cross-polarized component with respect to the incident wave. For the shown configuration in Fig. 1, vertically (z-) aligned plates (with relatively large dimensions compared to the wavelength), dipoles, and rectangular loops can provide a non-depolarizing condition. Moreover, in terms of phase variation, the geometrical symmetry of the scatterer about the rotation axis is a necessary condition leading to minimize the phase variation of the induced current on the scatterer for different incidence angles (and equivalently during the rotation of the scatterer). However, as the key contribution of the work, it will be shown that the geometrical symmetry and non-depolarizing configuration of a rotating scatterer are not sufficient conditions to modulate the backscattered wave only in magnitude, and the resonance effect should be included to suppress the phase variation due to Doppler.

To demonstrate the concept, two symmetric and non-depolarizing scattering configurations are considered as shown in Fig. 1b,e. The first one is composed of two equal-length z-aligned dipoles resonating in z-direction which are positioned symmetrically about the origin (Fig. 1b). The second one is a symmetric z-aligned rectangular loop resonating in x-direction (Fig. 1e). The resonance of the loop scatterer along x-direction (as shown in Fig. 1e) or more generally along the radial direction when the loop is rotated around z-axis is therefore called radial resonance since its direction is aligned with radial propagation of the incident wave in the proposed scattering configuration. Both structures are impinged at two different incidence angles ($\phi = \phi_{i1}, \phi_{i2}$), namely, oblique and normal incidence. For the oblique and normal incidence, the phase and magnitude distribution of the scattered field from the two dipoles are respectively shown in Fig. 1a,c, and those of the scattered field from the rectangular loop is respectively shown in Fig. 1d,f, while the scattered wave is co-polarized with the incident
<figure><img src="image_3.png" /><figcaption>**Figure 1.** The azimuth distribution ($0° < \phi_o < 360°$) of the phase and normalized magnitude of the scattered electric field from two dipoles for (a) oblique incidence (c) normal incidence. The azimuth distribution ($0° < \phi_o < 360°$) of the phase and normalized magnitude of the scattered electric field from rectangular loop for (d) oblique incidence (f) normal incidence. Symmetrical and non-depolarizing scattering structure (b) without radial resonance [two equal-length coupled dipoles] (e) with radial resonance [rectangular loop]. All the data are obtained in simulation.</figcaption></figure>

wave in the *xy*-plane for both structures. The magnitude of the scattered field have been normalized for all the cases to provide a better comparison between them. As it is illustrated in Fig. 1b, the phase of the induced current on the two dipoles depends on the wave path distances associated to each dipole (linked to the Doppler effect), which cause two dipoles to be in-phase for normal incidence (not shown in Fig. 1b) and to be out of phase for oblique incidence (shown in Fig. 1b). Consequently, the distribution of the phase and magnitude of the scattered field from two dipoles drastically changes as the incident angle varies (Fig. 1a,c). Nevertheless, for the rectangular loop with the radial resonance (resonating in *x*-direction), the phase of the induced current is imposed by the resonant mode, and it does not change for different incidence directions (for both normal and oblique incidence are the same, and it is shown in Fig. 1e) which means the Doppler effect has been suppressed due to the resonance. This phenomenon causes the azimuth distribution of the magnitude and phase of the scattered field to not vary for different incidence angles, while in addition, the phase distribution is constant as it is shown in Fig. 1d,f. Accordingly, this property can be used to realize a pure directional amplitude backscattering modulation. It should be mentioned that although the concept has been described for the fundamental resonant mode of the loop structure (the one shown in Fig. 1b), it is valid for all the resonant modes which produce a radial resonance.

The modulation process is presented in Fig. 2 for the monostatic scattering where $\bar{k}_s = -\bar{k}_i$ (Fig. 2a), and the rectangular loop is rotating around its symmetry axis (z-axis). As it is shown in Fig. 2b, the magnitude of the backscattered wave will be modulated during the rotation, while the envelope waveform of the received amplitude-modulated signal is exactly proportional to the directional reradiation pattern of the loop at the resonant mode. Applying a coherent IQ demodulation on the received signal yields the IQ constellation diagram
<figure><img src="image_4.png" /><figcaption>**Figure 2.** (a) Rotating rectangular loop in the monostatic scattering configuration. (b) Directional amplitude modulation induced on the backscattered wave by rotating loop. (c) The IQ constellation diagram of the demodulated received signal. The linear modulation path crosses the origin, and is limited in between the minimum and maximum of the reradiation pattern corresponding to the positions of the loop shown in (b).</figcaption></figure>

presented in Fig. 2c. The modulation path is limited between two points corresponding to the maximum and minimum of the reradiation pattern, and coincided on a line which passes through the origin.

**Theoretical analysis.** To verify the introduced concept, scattering from a rotating resonant rectangular loop is studied theoretically. A circuit model for the rectangular loop structure is proposed, and it is applied to develop an analytical model for the rotating loop. The model will be derived for general case (non-symmetrically rotating loop), and the effect of symmetry is highlighted for the special case (symmetrically rotating loop).

**Circuit model.** Based on the induction theorem²⁸, for a scatterer made of good conductors which is impinged by an electromagnetic wave, the induced current on the scatterer can perfectly describes the scattered field. Applying a circuit model to calculate the induced current on the scatterer (usually wire scatterers) is an efficient way to obtain the scattered field, provided that the circuit model is achievable. For a rectangular loop of wire, the circuit model is presented for the specific incident wave configuration shown in Fig. 3a. The rectangular loop with the length of 2L (in x-direction), the height of 2h (in z-direction), and wire radius of a/2 is non-symmetrically positioned in xz-plane such that $d_L$ and $d_R$ are respectively the distance between the left and right sides of the loop and the coordinate origin ($d_L + d_R = 2L$) as it is depicted in Fig. 3c. The rectangular loop is impinged by a time-harmonic z-polarized plane wave at the frequency of $f_0$ with the wave vector of $k_i = k\hat{k}_i$ ($k_i$ is placed at xy-plane) where $k = 2\pi f_0/c_0$ and $c_0$ is the free-space light velocity (Fig. 3a). For $h \ll L$ and $kh \ll 1$, it is shown in Ref.²⁹ that the wave-impinged rectangular loop can be modeled by a section of two-wire transmission line with propagation constant of $\gamma \approx jk$ and characteristic impedance of $Z_c$, which is excited by two voltage generators at both left and right ends with respective complex ElectroMotive Force (EMF) of $V_L$ and $V_R$ (Fig. 3e). The voltage sources are followed with two identical resistors $R_f(\lambda)$ which models the radiation contribution of the two short end sides of the loop, while the negligible radiation of the transmission line section is ignored. In terms of reradiation, the two currents flowing on the left and right ends of the transmission line, namely $I_L$ and $I_R$, are the only current components which contribute in the scattered field. $I_L$ and $I_R$ can be calculated for given $V_L$ and $V_R$ based on the superposition principle and by using the even-odd mode decomposition technique (see Supplementary A). However, the EMF of the voltage sources ($V_L$ and $V_R$) should be obtained from the incident electric field, which leads to define a vector effective length associated to each resonant mode of the loop structure. Without lack of generality, and remarking that the model can be developed for all the even and odd resonant modes of the loop structure, at the rest of development we choose the dominant even mode of the rectangular loop as the case study (see Supplementary A).

The vector effective length of the loop at the fundamental resonant mode ($L = (c_0/f_0)/4 = \lambda/4$) is derived based on the reciprocity theorem²⁸ using the method proposed in Ref.³⁰ (see Supplementary B). As the rectangular loop has been modeled by two reradiating current components on the left and right sides ($I_L = I_R = I_0$), there are two vector effective lengths ($\vec{l}_e^L$ and $\vec{l}_e^R$) obtained as

$$ \vec{l}_e^L(\theta, \phi) = 2h \sin \theta \cos(kL \sin \theta \cos \phi) e^{-jkL \sin \theta \cos \phi \hat{\theta}}, \quad (1) $$

$$ \vec{l}_e^R(\theta, \phi) = 2h \sin \theta \cos(kL \sin \theta \cos \phi) e^{jkL \sin \theta \cos \phi \hat{\theta}}, \quad (2) $$

which are respectively associated to $I_L$ and $I_R$, and yield the induced EMF on each side in terms of the incident electric field as

$$ V_L = \tilde{E}_i(r, \theta, \phi) \cdot \vec{l}_e^L(\theta, \phi) \Big|_{\substack{\theta=\theta_i, \phi=\phi_i \\ x=d_L, y=0, z=0}} \quad (3) $$
<figure><img src="image_5.png" /><figcaption>**Figure 3.** (a) Non-symmetrically positioned rectangular loop impinged by vertically polarized plane wave. (b) Scattering from non-symmetrically rotating rectangular loop. (c) Geometrical dimensions of the rectangular loop. (d) Scattering from two symmetrically rotating coupled dipoles. The left and right side radiating elements are respectively remarked with red and green dots. (e) The proposed circuit model for the wave-impinged loop.</figcaption></figure>

$$V_R = \vec{E}_i(r, \theta, \phi) \cdot \vec{l}_e^R(\theta, \phi) \Big|_{\substack{\theta=\theta_i, \phi=\phi_i \\ x=-d_R, y=0, z=0}} \quad (4)$$

where $\theta_i$ and $\phi_i$ are determined by incidence direction (i.e. $\hat{k}_i$). As the induced EMFs are determined by (3) and (4) in terms of the incident field, the induced current on the loop and consequently the scattered field can be calculated using the proposed circuit model. The model will be applied to obtain the scattering from rotating loop in the next step.

*Scattering from rotating rectangular loop.* For a time-harmonic incident plane wave, scattering from a moving scatterer, when its velocity is much smaller than $c_0$, can be solved for time-dependent sequences of stationary scattered fields using the quasi-stationary approximation theory$^{31}$. Thus, for the rotating loop, obtaining the current induced on the loop at each time instant during the rotation, while the loop is assumed fixed at the corresponding coordinate, makes it possible to derive the quasi-stationary scattered field.

Based on the developed circuit model, the induced current on the rotating loop can be obtained during the rotation by calculating the induced EMFs ($V_L$ and $V_R$) based on quasi-stationary approximation. Consider that the rectangular loop is non-symmetrically rotating around z-axis with angular rotation frequency of $\omega_r$ as it is shown in Fig. 3b, where the xyz is the reference Coordinate System (CS) and $x'y'z'$ is the local rotating CS attached to the rotating loop. The two coordinate systems are related by the rotation transformation with time-varying azimuth angle of $\phi_s(t) = \omega_r t = 2\pi f_r t$ which can be expressed as

$$\begin{bmatrix} x \\ y \\ z \end{bmatrix} = \begin{bmatrix} \cos \phi_s(t) & -\sin \phi_s(t) & 0 \\ \sin \phi_s(t) & \cos \phi_s(t) & 0 \\ 0 & 0 & 1 \end{bmatrix} \begin{bmatrix} x' \\ y' \\ z' \end{bmatrix}. \quad (5)$$

The electric field incident at the direction of $[\theta_i = \pi/2, \phi_i]$ can be written as

$$\vec{E}_i(\vec{r}) = E_0 e^{ik(x \cos \phi_i + y \sin \phi_i)} \hat{z}, \quad (6)$$

and the scattered field is observed at $[\theta_o, \phi_o]$ direction. The induced EMFs during the rotation can be easily obtained by (3) and (4), provided that the incident field is expressed in the rotating CS. Using (5), the incident electric field can be written as
$$ \vec{E}_i(\vec{r}') = E_0 e^{jk[x' \cos(\phi_i - \phi_s) + y' \sin(\phi_i - \phi_s)]} \hat{z}, \quad (7) $$

in the rotating CS. The time-varying complex envelope of the EMFs on the rotating loop ($\tilde{V}_L(t)$ and $\tilde{V}_R(t)$) can be obtained using (1)–(4) and (7), while (1)–(4) should be rewritten in the rotating CS ($x'y'z'$) by replacing $r, \theta, \phi$ with $r', \theta', \phi'$. The resultant expressions for $\tilde{V}_L(t)$ and $\tilde{V}_R(t)$ are given by

$$ \tilde{V}_L(t) = \vec{E}_i(r', \theta', \phi') \cdot \vec{l}_e^L(\theta', \phi') \Big|_{\substack{x'=d_L, y'=0, z'=0 \\ \theta'=\theta_i=\pi/2, \phi'=\phi_i-\phi_s}} = 2hE_0 \cos[kL \cos(\phi_i - \phi_s)] e^{jk(\frac{d_L-d_R}{2})} \cos(\phi_i - \phi_s), \quad (8) $$

$$ \tilde{V}_R(t) = \vec{E}_i(r', \theta', \phi') \cdot \vec{l}_e^R(\theta', \phi') \Big|_{\substack{x'=-d_R, y'=0, z'=0 \\ \theta'=\theta_i=\pi/2, \phi'=\phi_i-\phi_s}} = 2hE_0 \cos[kL \cos(\phi_i - \phi_s)] e^{jk(\frac{d_L-d_R}{2})} \cos(\phi_i - \phi_s), \quad (9) $$

which are identical to each other, and yield in-phase and equal induced currents $\tilde{I}_L(t) = \tilde{V}_L(t)/R_r$ and $\tilde{I}_R(t) = \tilde{V}_R(t)/R_r$ on the rotating loop during the time. Therefore, the model is perfectly compatible with the proposed physical concept that the radial resonance of the rectangular loop will suppress the phase variation caused by the Doppler effect along the rotating loop.

The quasi-stationary scattered field form the rotating loop at the fundamental resonance can be calculated based on the currents ($\tilde{I}_L(t)$ and $\tilde{I}_R(t)$) and by using the radiation integrals$^{28}$. After some mathematical manipulations, the final expression for the time-varying complex envelope of the scattered field is obtained as

$$ \tilde{\vec{E}}_s(\vec{r}, t) = \frac{j\eta k E_0 (2h)^2 e^{-jkr}}{2\pi r R_r} \sin \theta_o \cos[kL \cos(\phi_i - \phi_s)] \cos[kL \sin \theta_o \cos(\phi_o - \phi_s)] \\ \times e^{jk(\frac{d_L-d_R}{2})[\cos(\phi_i - \phi_s) + \sin \theta_o \cos(\phi_o - \phi_s)]} \hat{\theta}, \quad (10) $$

where $\eta$ is the free-space intrinsic impedance and $\phi_s = \omega_r t$ is the time-varying rotation angle. The time-varying complex envelope can be rewritten as

$$ \tilde{\vec{E}}_s(\vec{r}, t) = \frac{j\eta k E_0 (2h)^2 e^{-jkr}}{2\pi r R_r} \sin \theta_o m(t) e^{j\psi(t)} \hat{\theta}, \quad (11) $$

while $m(t)$ and $\psi(t)$ are respectively associated to the DAM and the Doppler phase modulation, and can be described as

$$ m(t) = \cos[kL \cos(\phi_i - \phi_s)] \cos[kL \sin \theta_o \cos(\phi_o - \phi_s)], \quad (12) $$

$$ \psi(t) = k(\frac{d_L - d_R}{2})[\cos(\phi_i - \phi_s) + \sin \theta_o \cos(\phi_o - \phi_s)]. \quad (13) $$

Obviously, for the symmetrical rotation ($d_L = d_R$), the Doppler modulation term will be completely cancelled $\psi(t) = 0$ which provides a pure DAM with the waveform of $m(t)$ linked to the reradiation pattern of the rectangular loop in its fundamental resonant mode. Thus, the proposed concept is proven theoretically by this demonstration.

**Simulation.** The presented analytical model and the idea of the paper are verified by simulation. The two important effects to achieve a pure DAM, the symmetrical rotation and radial resonance of the scatterer, will be addressed well through the simulation results. For this, the two structures presented in Fig. 1 are considered. A rotating rectangular loop of wire with $L = 15$ mm, $h = 2$ mm, and $a = 100$ µm (Fig. 3b,d), and two rotating coupled dipoles with $H = 30$ mm, $D = 30$ mm, and $a = 100$ µm (Fig. 3e) are simulated in the monostatic scattering configuration with $\theta_i = \pi/2, \phi_i = \pi/2$ and $\theta_o = \pi/2, \phi_o = \pi/2$. The simulation process is described in "Methods".

**Symmetrical and non-symmetrical rotation.** The symmetry of the scatterer about the rotation axis is the necessary condition to realize a pure DAM. In other words, even for the rectangular loop scatterer with radial resonance, it will be shown that the symmetrical rotation is essential for pure amplitude modulation, as it is predicted by the analytic model.

The phase and normalized magnitude of the induced currents on the left and right side of the rotating loop are extracted by simulation during one period of rotation $0 < \phi_s < 2\pi$. The simulation results for symmetrically ($d_L = d_R$) and non-symmetrically ($d_L \neq d_R$, without lacking generality $d_L = 3d_R$) rotating resonant loop at $f_0 = 5$ GHz are respectively shown in Figs. 4 and 5, both in perfect agreement with the results obtained from the model. Due to the radial resonance of the loop, for both cases (symmetrical and non-symmetrical rotation), the variation of the magnitude (Figs. 4a, 5a) and phase (Figs. 4b, 5b) of the left and right side currents [$\tilde{I}_L(t)$ and $\tilde{I}_R(t)$] are identical to each other. However, the phase variation during the time is constant only for the symmetrical rotation (Fig. 4b), while it is varying sinusoidally for the non-symmetric rotation (Fig. 5b) resulting a phase modulation on the backscattered field. This fact has been clearly demonstrated in terms of IQ representation for the normalized backscattered electric field, and it is shown in Figs. 4c and 5c, where $I = real(\tilde{E}_s)$ and $Q = imag(\tilde{E}_s)$. As predicted by the model, the limiting points of the modulation path is determined by the maximum and minimum of the amplitude modulating waveform $m(t)$ which is associated to the reradiation
<figure><img src="image_6.png" /><figcaption>**Figure 4.** (a) Normalized magnitude, (b) phase of the induced currents on the left and right side of the symmetrically rotating rectangular loop during one period of rotation, obtained from the model and simulations. (c) IQ constellation of the time-varying complex envelope of the normalized backscattered field from symmetrically rotating loop in simulation and theory.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 5.** (a) Normalized magnitude, (b) phase of the induced currents on the left and right side of the non-symmetrically rotating rectangular loop during one period of rotation, obtained from the model and simulations. (c) IQ constellation of the time-varying complex envelope of the normalized backscattered field from non-symmetrically rotating loop in simulation and theory.</figcaption></figure>

pattern of the resonant mode in the loop. Finally, it is worth-mentioning that the modulation path for the symmetric case ($d_L = d_R$, $\psi(t) = 0$) is coincided on a line which pass through the coordinate origin (pure amplitude modulation), while that for the non-symmetric case ($d_L = 3d_R$, $\psi(t) \neq 0$) follows a drop shape path as it is shown in Fig. 5c (mixed amplitude-phase modulation).

**Symmetrical rotation: rectangular loop vs. two coupled dipoles.** The importance of the radial resonance to achieve directional amplitude modulation based on rotating scatterer is clarified by simulation for two symmetrically rotating scatterers, namely the rectangular loop and two equal length coupled dipoles as it is depicted in Fig. 3d.

For the rotating rectangular loop and two rotating dipoles, the simulation results are presented in Fig. 6a,b for the phase and normalized magnitude of the induced current on the left and right side elements during one period of rotation. Both simulations for the loop and the two dipoles have been done at $f_0 = 5$ GHz which is the fundamental resonance frequency of both scatterers. In contrast to the loop in which $\tilde{I}_L(t) = \tilde{I}_R(t)$, the phase and magnitude variation of the induced currents on the left and right side dipoles are completely different. This fact, as it was discussed in Fig. 1, proves that, even in a symmetrical rotation, without the radial resonance, it is impossible to induce a pure DAM on the backscattered wave using a rotating scatterer. The IQ diagram of the normalized backscattered field shown in Fig. 6c also confirms that, although the modulation path for the two rotating dipoles is a line, however, it does not pass through the origin, and consequently it represents a mixed amplitude-phase modulation.

**Experimental verification.** At the final step, the proposed concept is validated through the measurement results. In this regard, the rectangular loop has been realized with a 250 µm radius copper wire, while the length and height of the loop are respectively 2L = 29.8 mm and 2h = 3.9 mm. The fundamental resonance of the loop
<figure><img src="image_8.png" /><figcaption>Figure 6. (a) Normalized magnitude, (b) phase of left-side and right-side induced currents on the two symmetrically rotating dipoles and the rectangular loop. (c) IQ constellation of the time-varying complex envelope of the normalized backscattered field from symmetrically rotating loop and two symmetrically rotating dipole. All the data are obtained in simulation.</figcaption></figure>

has been measured at $f_0 = 4.98$ GHz, and the experimental study has been done using the monostatic setup which will be described in “Methods” section.

**Environment contribution.** Until this point, in the theoretical model and simulation, it has been assumed that the rotating loop is located in the free space. Nevertheless, to retrieve the same modulation waveform as the one obtained in theory from the measured backscattered signal, it is necessary to consider the contribution of the real environment in which the experiment has been done. Accordingly, it is supposed that in the measurement process the rotating loop is the only moving element in the environment, and all the other object in the environment are stationary. In addition, it is assumed that only the direct reflection of the modulated backscattered wave from the rotating loop is captured by the receiving antenna and the multi-reflection of the modulated wave with other present objects is much less significant compared to the direct reflection. With these assumptions, the backscattering contribution of the environment can be considered as a time-invariant term which is added to the complex envelope of the quasi-stationary backscattered field from the loop (10), and it is described as

$$ \vec{E}_s^{\text{meas}}(\vec{r}, t) = \vec{E}_s(\vec{r}, t) + \vec{E}_s^{\text{env}}(\vec{r}), \quad (14) $$

where $\vec{E}_s^{\text{env}}(\vec{r})$ will be observed as an added constant complex voltage $\rho_{\text{env}}$ in the demodulated received signal. Accordingly, in the IQ diagram representation, the contribution of the environment cause the modulation path to be translated along the constant vector of $\rho_{\text{env}}$ in the IQ plane. This is an important point for retrieving the directional amplitude modulation waveform associated to the reradiation pattern of the loop at its resonant modes. Actually, since for all the resonant modes of the loop (and of course for the fundamental mode) the reradiation pattern has some null points (zero magnitude as shown in Fig. 2b), the amplitude of the received signal will be more affected by noise at these points which makes the null points detectable on the modulation path in the IQ plane. Thus, the contribution of the environment ($\rho_{\text{env}}$) can be easily removed by putting the null point at the origin of the IQ plane, and then the modulation waveform is retrieved accurately. Note that this compensation cannot be done for the two symmetrically rotating dipoles since the modulation path does not cross the origin of the IQ plane. The measurement results in the next step will be processed based on this fact, and the DAM is validated experimentally.

**Measurement results.** The measurements have been done for symmetrically and non-symmetrically rotating loop in time and frequency domain using the bench described in “Methods” section. After coherent IQ demodulation, the captured raw IQ data are presented for symmetrical and non-symmetrical rotation respectively in Fig. 7a,d. To represent the modulation path and to obtain the modulation waveform in a convenient form, the noise fluctuation on the raw IQ data have been reduced by applying moving average, and the resultant averaged data are presented in Fig. 7b,e for both cases. Based on the averaged IQ data, the modulation path is depicted for symmetrical and non-symmetrical rotation respectively in Fig. 7c,f [black curves], which are compatible with the predicted path in theory for both cases (line path for the symmetric rotation and drop shape path for the non-symmetric rotation).

As it has been discussed previously, to retrieve the modulation waveform accurately, the contribution of the stationary environment $\rho_{\text{env}}$ should be removed from the data to achieve the original modulation path (Fig. 7c,f green curves). For each case, the original modulation path and the waveform are investigated on one period of rotation indicated in Fig. 7b,e by shadowed region ($27 < t < 52$ ms for symmetric rotation and $56 < t < 81$ ms for non-symmetric rotation), and are shown in Fig. 8. The modulation path can be divided to four quarter sections which correspond to the four quarter of the rotation with the respective orientation of the loop which is clearly illustrated by color markers in Fig. 8a,e for the symmetric case, and in Fig. 8c,f for the non-symmetric case. Moreover, these four sections are related to the modulation waveform in each case (Fig. 8b,d), while the
<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 7.** Obtained raw IQ data in measurement for scattering from (a) symmetrically, (d) non-symmetrically rotating rectangular loop, and the averaged IQ data for scattering from (b) symmetrically (e) non-symmetrically rotating rectangular loop. The measured IQ constellation diagram plotted based on the averaged data for (c) symmetrically (f) non-symmetrically rotating rectangular loop. The original measured modulation paths are shown in black color, and after the subtraction of the environment contribution (red arrow) they are shown in green color.</figcaption></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /></figure>

<figure><img src="image_18.png" /></figure>

<figure><img src="image_19.png" /></figure>

<figure><img src="image_20.png" /><figcaption>**Figure 8.** The modulation path during one period of rotation for (a) symmetrically, (c) non-symmetrically rotating loop. The magnitude and phase variation of the retrieved modulating waveform from measured signal for (b) symmetrically, (d) non-symmetrically rotating loop during one period of rotation. Each period of rotation is divided to four sections shown in green, red, blue and black color on the modulation path and modulating waveform. The respective orientation of the loop is given at each section for (e) symmetrical, (f) non-symmetrical rotation.</figcaption></figure>

<figure><img src="image_21.png" /></figure>

<figure><img src="image_22.png" /></figure>

<figure><img src="image_23.png" /></figure>

<figure><img src="image_24.png" /><figcaption>**Figure 9.** The magnitude and phase of the measured amplitude modulating waveform for (a) symmetrical, (d) non-symmetrical rotation of the loop compared with those predicted by the model. The measured PSD of the received signal backscattered from (b) symmetrically, (e) non-symmetrically rotating loop. The calculated Fourier coefficients of the time-varying complex envelope of the backscattered field from (c) symmetrically, (f) non-symmetrically rotating loop.</figcaption></figure>

magnitude variation is associated to the reradiation pattern of the fundamental resonance of the loop, and the phase variation is constant for the symmetric rotation (pure DAM with suppressed Doppler), and is sinusoidal for the non-symmetric rotation (mixed amplitude-phase modulation).

Finally, the normalized modulation waveform obtained from the measurement are compared with those predicted by the model in Fig. 9a,d, which shows a perfect agreement between the model and measurements. Furthermore, the Power Spectral Density (PSD) of the demodulated signal $S_r(f)$ has been measured for the
<figure><img src="image_25.png" /><figcaption>**Figure 10.** Backscatter communication based on rotating scatterers with (a) pure Doppler modulation, (b) pure polarization modulation, and (c) pure DAM. (d,e) backscatter communication with combined DAM and Doppler modulation. The phase change of the induced current on the scatterer during the rotation in (d,e) is shown by two red dots for normal incidence, and one blue and one green dot for oblique incidence.</figcaption></figure>

symmetric and non-symmetric cases, and are respectively shown in Fig. 9b,e. The measured $S_r(f)$ can be compared with the calculated Fourier coefficients [$C_n$] of the backscattered field expressed in (10), which are presented in Fig. 9c,f. The measured and calculated frequency representation of the backscattered field contain components at the harmonics of the rotation frequency ($nf_r$; $n = 0, \pm 1, \pm 2, ...$). For the symmetric case, both in measurement and model, only the even harmonics exist (i.e. $2nf_r$), while for the non-symmetric case all the harmonics are present. Therefore, the frequency domain measurement results are also in accord with the proposed model in terms of harmonics and their respective ratios.

## Discussion

The backscatter modulation based on rotating resonant loop scatterer was proposed to introduce the first feasible example for DAM. The effect of radial resonance was utilized to suppress the phase modulation induced by Doppler which leads to a pure DAM on the backscattered wave. The concept was studied with an analytical model, and the model was verified by numerical simulations and measurements.

Worth mentioning that, from the physical point of view, the proposed concept is quite specific since it clarifies that in a general backscatter modulation based on a moving scatterer, the modulation can be originated from three different physical phenomena as Doppler modulation, polarization modulation, and DAM, while the latter has not been clearly addressed before, and it has been usually merged with Doppler modulation. For rotating scatterers which have a radial velocity towards the incident wave, the Doppler modulation should be apparently involved. However, the suggestion of a backscatter modulation based on a rotating scatterer in which the modulation is done just proportionally to the directional property of the scatterer, perfectly clarifies the pure DAM and the Doppler suppression due to resonance effect, and separates its effect from the two other phenomena (Doppler and polarization modulation). To clarify this fact, Fig. 10a-c illustrates the three different rotation configurations which respectively provides a pure Doppler modulation²⁴, pure polarization modulation²⁶, and pure DAM (this work). Figure 10a presents the pure Doppler (phase) modulation in which the wire scatterer is aligned with the z-axis and is rotating around it. In this configuration, the backscattered field keeps the same polarization state as the incident field in all azimuth and elevation directions. Also, for when the radius of rotation is negligible compared to observation distance, since the reradiation pattern of the scatterer does not change during the rotation, there is no modulation in terms of magnitude (DAM) and polarization on the backscattered field. However, the phase (Doppler) modulation is shown in Fig. 10a by the phase change of the induced current during the rotation while the two different phases associated with two different radial distances are depicted by red and green color dots on the two endpoints of the scatterer. After, Fig. 10b shows the configuration for pure polarization modulation where the wire scatterer is lying in the xz-plane and is rotating around the y-axis and the incident wave impinges normally on the scatterer. In this case, the polarization state of the backscattered field is modulated during the rotation and the linear polarization state of the backscattered wave is rotating in its wave-front plane, while the phase and magnitude of the induced current on the scatterer and its reradiation
<figure><img src="image_26.png" /><figcaption>**Figure 11.** (a) The measurement bench in the anechoic chamber. (b) The circular foam support used for rotating the rectangular loop. (c) Fabricated rectangular loop.</figcaption></figure>

pattern remain unchanged during the rotation. Since the radial distance of the scatterer respect to the incidence direction is fixed, the invariant phase of the induced currents at the two endpoints are depicted with red dots during the rotation. Finally, the resonant symmetrically rotating loop (proposed by this work) is shown in Fig. 10c where the backscattered wave keeps the same polarization state as the incident wave. The Doppler modulation due to the different radial distances is suppressed by the resonance effect and consequently, the reflected wave is modulated only due to the variation of the reradiation pattern of the scatterer during the rotation. In addition, to have some general examples (with the same considered scatterers) in which a combined modulation is involved, Fig. 10d,e presents two configurations (respectively horizontally rotating non-resonant dipole and symmetrically rotating non-resonant loop) while the backscattered wave experiences both Doppler modulation and DAM. Note that both configurations shown in Fig. 10d,e will also provide a pure DAM if the scatterer is resonant so that the Doppler modulation is suppressed.

Last but not least, since the pure DAM induced on the backscattered wave is proportional to the reradiation pattern of the resonant scatterer which usually has some null points, using coherent demodulation, the contribution of the environment can be easily removed to accurately retrieve the associated reradiation pattern from the received signal. This fact raises some potential applications for this kind of pure amplitude backscatter modulation such as pattern identification (where the reradiation pattern of the scatterer at each resonant mode can be assigned to a specific ID or different scatterers with specific patterns can be used for each ID like²²), localization (where the general bistatic expression of the time-varying backscattered field (Eq. 10) is utilized to find the incidence $\phi_i$ and observation $\phi_o$ angles based on the retrieved amplitude profile $m(t)$ which can be applied to determine the location of the rotating scatterer with respect to the transmitter and receiver in a 2D localization scenario based on Angle-of-Arrival (AoA) algorithm.), and wireless sensors to detect the asymmetry of rotation in non-destructive monitoring of rotating mechanical systems (for example the loop can be attached to a rotating shaft and whenever the shaft is not rotating around its center [$d_L \neq d_R$ in Eq. (10)] the undesired phase modulation $\psi(t) \neq 0$ can be detected as a problem occurred to the rotation of the shaft). It is a key point that in all these proposed applications, the communication process is based on a pure DAM which occupies less bandwidth compared to the Doppler modulation for a definite radius of rotation (theoretically infinite number of micro-Doppler harmonics¹⁵,²⁷) and it can be demodulated even with simple envelope detectors. It should be mentioned that, although this kind of DAM for backscatter communication will be more suffered by the environment interference compared to phase (Doppler) modulation, this fact can be justified by the natural trade-off between the bandwidth and noise immunity which has always existed in telecommunication.

## Methods

**Frequency domain simulations.** As it was mentioned previously, the quasi-stationary approximation can be used to solve scattering from moving scatterers. This means that the problem can be analyzed by a sequence of frequency domain simulations which correspond to the sequence of time instants along the movement. Accordingly, in each simulation the location of the scatterer should be changed such that during the full sequence of simulations, the full trajectory of the scatterer is covered. For this purpose, since the considered scatterer in this work is a wire scatterer, Method of Moments (MOM) has been selected as the faster computational approach to run the sequence of simulations with enough accuracy. The NEC2 package based on MOM has been used for doing simulation in this work. To obtain the result for a full rotation of the structure, the sequence of the simulations have been done at the desired frequency for each 1° of rotation. The currents have been extracted at the middle point of the left and right radiating elements ($\tilde{I}_L$ and $\tilde{I}_R$) in each structure (rectangular loop and two coupled dipoles) as it is shown in Fig. 3 by the green and red dots.

**Measurement setup.** The measurement bench has been organized in an electromagnetic anechoic chamber as it is shown in Fig. 11a. Two horn antennas (A.H. Systems, INC. SAS-571) with vertical linearly polarized
radiation have been utilized as transmitting and receiving elements which are respectively connected to the signal generator (HP8720D in CW mode) and spectrum analyzer (Tektronix RSA3408A). The antennas are closely positioned to realize the monostatic scattering configuration. Moreover, the measurement bench and the instruments are located in the side lobe of the TX and RX antennas to reduce the undesired effect of the bench presence in the anechoic chamber. To preserve the phase variation of the backscattered signal accurately (coherent transmission and reception) the local oscillator of the signal generator and that of the spectrum analyzer have been synchronized by using a 10 MHz common reference signal. The rotation has been realized by a motor which rotates a circular foam support with a programmed speed (Fig. 11a,b). A rectangular loop has been fabricated using 250 µm copper wire (Fig. 11c), and it was inserted carefully into the slot on the foam support such that the orientation of the loop is fixed during the rotation, and it is exactly in accord with the model (Fig. 11b). The rotational frequency has been set at $f_r = 40$ Hz during the measurement, while it can be set at any other frequency without affecting the results, as long as the transmitter and the receiver are coherent. A Continuous Wave (CW) signal at the frequency of $f_0 = 4.98$ GHz (the measured fundamental resonance frequency of the fabricated rectangular loop) has been injected into the transmitting antenna. The backscattered wave from the rotating loop has been captured by the receiving antenna, while the distance between the antennas and the rotating support is 1 m to fulfill the far-field radiation condition for both antennas and the scatterer. The received signal has been acquired by the spectrum analyzer in time (for IQ data) and frequency domain (for PSD data). The baseband demodulated I and Q data have been recorded for 160 msec acquisition time with a sampling rate of 64 kS/s, and the PSD of the signal has been measured across a 800 Hz span with a resolution of 1 Hz.

## Data availability

The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.

Received: 17 August 2022; Accepted: 16 December 2022

Published online: 20 December 2022

## References

1. Stockman, H. Communication by means of reflected power. Proc. IRE **36**, 1196–1204 (1948).

2. Vena, A., Perret, E. & Tedjini, S. Chipless RFID Based on RF Encoding Particle: Realization, Coding and Reading System (Elsevier, Berlin, 2016).

3. Harrington, R. F. Theory of loaded scatterers. In Proc. Institution Electrical Engineers, Vol. 111, 617–623 (IET, 1964).

4. Harrington, R. Electromagnetic scattering by antennas. IEEE Trans. Antennas Propag. **11**, 595–596 (1963).

5. Nikitin, P. V. & Rao, K. S. Antennas and propagation in UHF RFID systems. In 2008 IEEE International Conference on RFID, 277–288 (IEEE, 2008).

6. Van Huynh, N. et al. Ambient backscatter communications: A contemporary survey. IEEE Commun. Surv. Tutor. **20**, 2889–2922 (2018).

7. Nikitin, P., Rao, K. & Martinez, R. Differential RCS of RFID tag. Electron. Lett. **43**, 431–432 (2007).

8. Lietzén, J., Liljemark, A., Duan, R., Jäntti, R. & Viikari, V. Polarization conversion-based ambient backscatter system. IEEE Access **8**, 216793–216804 (2020).

9. Fara, R. et al. Polarization-based reconfigurable tags for robust ambient backscatter communications. IEEE Open J. Commun. Soc. **1**, 1140–1152 (2020).

10. Tedjini, S. et al. Hold the chips: Chipless technology, an alternative technique for RFID. IEEE Microw. Mag. **14**, 56–65 (2013).

11. Newburgh, R. G. & Borgiotti, G. V. Backscattered Spectra from Rotating and Vibrating Short Wires and Their Relation to the Identification Problem. Tech. Rep. (Air Force Cambridge Research Labs, 1975).

12. Reynolds, M. S. A 500 °C tolerant ultra-high temperature 2.4 GHz 32 bit chipless RFID tag with a mechanical BPSK modulator. In 2017 IEEE International Conference on RFID (RFID), 144–148 (IEEE, 2017).

13. Nikitin, P. Leon theremin (lev termen). IEEE Antennas Propag. Mag. **54**, 252–257 (2012).

14. Chaudhari, K. & Pradeep, T. Spatiotemporal mapping of three dimensional rotational dynamics of single ultrasmall gold nanorods. Sci. Rep. **4**, 1–9 (2014).

15. Kozlov, V., Filonov, D., Yankelevich, Y. & Ginzburg, P. Micro-Doppler frequency comb generation by rotating wire scatterers. J. Quant. Spectrosc. Radiat. Transfer **190**, 7–12 (2017).

16. Kozlov, V., Vovchuk, D., Kosulnikov, S., Filonov, D. & Ginzburg, P. Micro-Doppler signatures of subwavelength nonrigid bodies in motion. Phys. Rev. B **104**, 054307 (2021).

17. Bell, M. R. & Grubbs, R. A. JEM modeling and measurement for radar target identification. IEEE Trans. Aerosp. Electron. Syst. **29**, 73–87 (1993).

18. Suresh, P., Thayaparan, T., Obulesu, T. & Venkataramaniah, K. Extracting micro-Doppler radar signatures from rotating targets using Fourier-Bessel transform and time-frequency analysis. IEEE Trans. Geosci. Remote Sens. **52**, 3204–3210 (2013).

19. Thayaparan, T. et al. Analysis of radar micro-Doppler signatures from experimental helicopter and human data. IET Radar Sonar Navig. **1**, 289–299 (2007).

20. Krasnov, O. A. & Yarovoy, A. G. Radar micro-Doppler of wind turbines: Simulation and analysis using rotating linear wire structures. Int. J. Microw. Wirel. Technol. **7**, 459–467 (2015).

21. Lazaro, A., Lazaro, M., Villarino, R. & De Paco, P. New radar micro-Doppler tag for road safety based on the signature of rotating backscatters. IEEE Sens. J. **21**, 8604–8612 (2021).

22. Verdú, J., Guerrero, E., Parrón, J., Lázaro, A. & de Paco, P. Star-shaped wheel for mechanical micro-Doppler modulation. IEEE Antennas Wirel. Propag. Lett. **20**, 2452–2456 (2021).

23. Henry, D., Pons, P. & Aubert, H. Range-Doppler imaging approach for the wireless reading of mechanical sensors. In 2021 15th European Conference on Antennas and Propagation (EuCAP), 1–5 (IEEE, 2021).

24. Azarfar, A., Barbot, N. & Perret, E. Chipless RFID based on micro-Doppler effect. IEEE Trans. Microw. Theory Tech. **70**, 766–778 (2021).

25. Wong, J., Reed, I. S. & Kaprielian, Z. A. A model for the radar echo from a random collection of rotating dipole scatterers. IEEE Trans. Aerospace Electron. Syst. **2**, 171–178 (1967).

26. Barbot, N. & Perret, E. linear time-variant chipless RFID sensor. IEEE J. Radio Freq. Identif. **6**, 104–111 (2021).

27. Chen, V. C., Li, F., Ho, S.-S. & Wechsler, H. Micro-Doppler effect in radar: Phenomenon, model, and simulation study. IEEE Trans. Aerosp. Electron. Syst. **42**, 2–21 (2006).
28. Harrington, R. F. *Time-Harmonic Electromagnetic Fields* (Wiley-IEEE Press, Berlin, 2001).

29. King, R. The rectangular loop antenna as a dipole. *IRE Trans. Antennas Propag.* **7**, 53–61 (1959).

30. Wunsch, A. The vector effective length of slot antennas. *IEEE Trans. Antennas Propag.* **39**, 705–709 (1991).

31. Cooper, J. & Wilcox, C. H. Scattering by moving bodies: The quasi stationary approximation. *Math. Methods Appl. Sci.* **2**, 131–148 (1980).

**Acknowledgements**

This work was supported by the European Research Council (ERC) through the European Union's Horizon 2020 Research and Innovation Program (ScattererID) under Grant No. 772539.

**Author contributions**

All authors conceived the idea. A.A. developed the analytical model, performed the experiments, and wrote the manuscript. All authors analyzed the results and contributed to the development of the idea. The manuscript was reviewed and proofread by all authors.

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-022-26609-w.

**Correspondence** and requests for materials should be addressed to A.A.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_27.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022