The increasing usage of electromagnetic polarization in telecommunication systems is well known and documented, ranging from applications such as diversity and frequency reuse to electromagnetic compatibility and future 4 k/8 k resolution Digital Television Broadcasting  -  . The polarization of a wave is defined by the path described by the electric field vector at a fixed point in space as time varies, and the polarization of an antenna is defined as the polarization of the wave it radiates when operating in transmitting mode  . The possibility of generating 3D polarizations may lead to enhancements in the existing telecommunication systems or even the development of new configurations and applications.
This work presents a way of generating 3D polarizations and field distributions, where the electric field intensity vector can be oriented in any direction in space. Therefore, it is possible to generate 3D polarizations and/or fields with different shapes, as time varies, such as cylinders, spheres, and others. In addition, the orthogonality between 3D polarizations is also discussed. Antennas naturally generate a polarization that is orthogonal to the desired one, normally referred to as cross polarization, but the idea here is to generate the orthogonal polariza- tions in a controlled manner, similarly to the use of orthogonal polarizations in frequency reuse systems. Cross polarization is normally undesired, with the exception of the application introduced in  .
It was previously demonstrated analytically and numerically by the authors in  that 3D polarizations can be generated at an observation point inside a 3D antenna array, where the elements are facing inwards. Volumetric control of the intensity of electric fields and sound has been previously published     , but as far as the authors know, generation and application of 3D polarizations and fields have not yet been investigated in the way herein presented.
It was demonstrated in  that ideal antennas can generate 3D polarizations, and here the concept is extended for pyramidal horn antennas, as well as spatially distributed electromagnetic fields. A new method introduced by the authors for optimum design of these antennas was also employed in order to ensure greater gain stability with respect to manufacturing variations  .
A brief summary of the concept discussed in  is presented for completeness, also illustrating different scenarios than the ones covered in  , and it is shown that antenna arrays can be used for spatial multiplexing in wireless communi- cation systems by controlling the fields in 3D space, in order to reduce interference among users. This might be critical in the actual trending of the evolution of cellular systems towards nano and phantom cells  . Another possible application for the field and polarization control herein discussed is the generation of torus-knotted polarizations and fields, where the latter can be encountered in nature around planets, stars and other celestial objects, as indicated historically in Radioastronomy as well as possibly to be further studied with state-of-the-art Radiotelescopes  , and may be of importance for plasma confinement and nuclear magnetic resonance devices   .
2. Three-Dimensional Polarization with Ideal Dipoles (Near-Field) and Orthogonality
The antenna array shown in Figure 1 is used to generate 3D polarizations at the observation point O, and was derived from the antenna array configurations introduced by one the authors in   . All the computer simulations employed the computer package MathCad  . Note that with this configuration, each
Figure 1. Array of 3 ideal dipole antennas in free space. Although not discussed in this work, further consideration can also be given to array configurations of moving elements, including ones in known paths.
antenna produces fields only in the radial direction. This was done to simplify the equations discussed in  , even though it is possible to change the configuration for any desired direction. If the amplitude and phase of the antennas are constant, the resulting fields obtained by superposition (it is assumed a linear, homogenous, time invariant medium) normally lead to an elliptical polarization. However, by controlling the variations of the amplitude and/or phase of the elements it is possible to rotate this ellipse with time, and hence yielding a 3D polarization.
The electric field generated by one antenna is frequency dependent, however the figure formed by the fields generated in the array of ideal dipole are frequency independent if the same frequency is used for the all antennas. This occurs because the proportionality between the electric fields generated by individual antennas is kept constant when the frequency is changed. Thus, for illustration purposes, without loss of generality, we use f = 100 MHz, I = 1 mA, h = λ/100 (the length of all 3 antennas placed at a distance equal to λ from point O). It is important to note that for the scenarios considered in this study there is no propagation, so the 3D polarization needs to be visualized as the polarization of the resulting electric field, not the wave. Alternate configurations that allow the 3D polarization to move within a desired path are currently being investigated (by determining different sets of feed currents for each desired location along the trajectory, for example), but the theory does not lose its importance nor generality as herein presented. By varying continuously the phase of the ideal dipole P2 in Figure 1, i.e. ψ2 in the interval 0 to 2π, with ψ1 = ψ3 = 0˚ respectively for the other two dipoles, the resulting polarization is shown in Figure 2. In general, 3 components (basis) of electric fields are necessary in order to generate a 3D polarization. In the case of the components being orthogonal, they form an orthogonal base; the reader is referred to  for further details.
3. Three-Dimensional Polarization with Pyramidal Horn Antennas (Far-Field)
In order to show that it is also possible to generate 3D polarizations with antennas other than the ideal dipoles, we employed the array of pyramidal horns shown in Figure 3. The geometry of each horn is depicted in Figure 4.
Figure 2. Cylindrical 3D polarization obtained with the configuration of Figure 1 with ψ1 = ψ3 = 0˚ and ψ2 varying continuously in the interval 0 to 2π.
Figure 3. Array of pyramidal horn antennas.
Figure 4. Geometry of a pyramidal horn antenna.
The electric field of each horn can be approximated by (far-field) 
where E0 is the absolute value of the intensity of the electric field at the aperture (assuming the mode TE01). The electric field at the aperture is given by
where. The functions I1 and I2 are given by
where C(x) and S(x) are the Fresnel functions, and are given in  . Although the equations are not exact, the generation of 3D polarization from them can be demonstrated without loss of generality since a more exact approach would lead to different values in the radius of the sphere, which can always be compensated by adjusting properly the amplitude and/or phases of the element excitations (which in this example are the electric fields in the aperture of each antenna). Using electric fields distributions in the aperture of each horn as
where ζ is a variable controlled by the amplitude variation of the electric field in the aperture of each antenna, it can be shown that the total electrical field (far-field analysis) at the observation point O is given by
which in the time domain is
where C0 and φc are respectively the absolute value and phase of (valid for O at the origin)
The method introduced by the authors in  was employed to design each pyramidal horn (all equal to each other), and is herein omitted for simplicity. The optimum design ensures greater gain stability with respect to variations of the aperture dimensions  . The design gain was 18 dBi, frequency of 6 GHz and the waveguide WR137 (R1 = 19.94 cm, R2 = 18.2 cm, A = 18.46 cm, B = 14.19 cm). The antennas are away 103λ from the origin (far-field), and E0 = 1 V/m. The spherical polarization shown in Figure 5 is obtained by varying ζ continuously.
4. Spatial Multiplexing by Electric Field Control
In the previous sections the field control was done in only one point in the space; the reader is referred to  for a sensitivity analysis of degradations around that point. For spatial multiplexing, the fields should be controlled in many points in the space. This section shows the development of a method for controlling the electric field using three-dimensional array of antennas in multiple points in space, and how it can be applied in spatial multiplexing during the transmission of information in wireless digital communication systems.
The number of users in a given region can be increased with the new method, which is based on space division multiplexing, using the same frequency and time resources, with the least possible interference between users. Thus, the channel capacity can be increased. The interference can be further reduced between nearby users, increasing the signal to noise ratio (S/N) with the control of orthogonal polarizations.
4.1. Method for Electric Field Control
When m antennas are used for the control of n points in space, individual electric fields generated by each antenna are added to obtain the resulting field. The
Figure 5. Spherical polarization with pyramidal horns (far-field).
most direct method for this is to decompose the spherical components ar and aθ directions in cartesian coordinates ax, ay and az:
where arx, ary and arz respectively represent the decomposition of the components of the unit vector ar in the directions ax, ay and az. Also, aθx, aθy and aθz respectively represent the decomposition of the components of the unit vector aθ in the directions ax, ay e az.
Considering an antenna located at point P1, and an observation point in O1, the equations below are obtained in cartesian coordinates:
where P1 is the position vector of antenna 1, S1 is the vector direction of antenna 1, and O1 is the position vector of the observation point 1. Superscript T in the equations was used for simplification, and indicates the transpose of the row vector. All these vectors are referenced to the origin of the cartesian system. P1 and O1 start from the origin and go to the specified points, and the distances represent the amplitudes. S1 starts from point P1 and is defined here as a unit vector, and its amplitude equals 1. The distance vector D11 between O1 and P1 is. The angle formed between the vector D11 and S1 is given by
The distance vector D11 has the same direction of the vector ar11. In general, establishing the variable u as being equal to x, y or z, aru11 can be determined as:
With arx11, ary11, and arz11 values obtained through (17), aθx11, aθy11, and aθz11 can be determined as:
Using, the contribution of the excitation source can be separated of the contribution responsible for the position and direction of the antenna:
where φEu11 is the argument of Eu11. The coefficient Cu11 has the unit Ω/m. Cu11 and φEu11 are calculated as
The phase and amplitude can be grouped in a single complex variable σu11 as shown below
Once the values of the length of the antenna and the frequency of the excitation source are determined in the vacuum, the coefficient σu11 will vary only with changes in the position and direction of the antenna, as well as the distance of this to the observation point. With m antennas and n observation points (m = 3n), Equations (25) and (26) can be generalized to
where i and l are integers number, with,. Also, u = x, y or z. Since n can take any integer value, the electric field control in space can exist in the infinite points. In full matrix notation, this can be written as:
Calling s the matrix of coefficients, I the current vector and E the vector fields, the values of the excitation sources of the antennas can be determined by:
where is the inverse matrix of s. In practice, numerical methods such as Gaussian elimination or LU decomposition are used because they are more efficient and accurate than the numerical computation of the inverse matrix  . Thus, the values of the amplitudes and phases of the current phasor are obtained with the desired values of E at each point, as well with the calculation of the complex coefficient matrix s.
4.2. Digital Signal Transmission and Spatial Multiplexing Using the Method for Electric Field Control
In this section the volumetric antennas array and fields control can be employed for the transmission of digital signals as well as to minimize the interference by spatial multiplexing. For the simulation, the antenna array presented in Figure 6 was used.
For the observation point O1, the main antennas are 1, 2 and 3 because they contribute more significantly to the resulting electric field at this point, since the value of the coefficient Cu has a higher value. Similarly, for point O2, the main antennas are 4, 5 and 6. In the configuration of Figure 6, in each observation point the antennas are oriented such that the θ angle calculated by (26) is 90˚.
Assuming, at an operating frequency of 100 MHz, the observation points are defined as and. The vectors positions and directions of the antennas are:,
, , ,
, , , and
Note that the mutual couplings between the transmitting antenna and the reception points are negligible, since the distances between the points of observations and the antennas are in far-fields. The amplitudes of the mutual impedances of the antennas in transmission and reception are considerably reduced if the distances between them are large  . Within this context, we obtain a matrix of complex coefficients s:
Figure 6. Array with 6 ideal antennas, working together for the electric field control of two observation points.
It was considered for the simulation digital transmissions ASK with binary signals with ma < 1. Although other modulations can be employed, we adopted this modulation for being simpler than the others, given that the focus of the simulation is to illustrate the transmission of signals to two different receiving points, sustaining isolation between the points at a same frequency and polarization.
For the simulation, it was also considered as bit 0 an amplitude of electric field equals to 2 × 10−3 V/m, and bit 1 an amplitude equals to 4 × 10−3 V/m. In the absence of transmission, the amplitude must be zero. A circular polarization was
used in the plane, noting that all antennas need to be excited
to generate the desired electric field.
The simplest condition happens when there is no transmission at both observation points. In this configuration, the desired fields E0 and E1 are zero, and current feeds of the antennas are zero. For other conditions, Equation (30) was used to obtain the amplitudes and phases of the currents. Table 1 presents all results.
In the transition times of the desired signals, the currents must be changed in these moments. According to the table, all antennas must be fed, even when there is transmission at only one point. This is due to the fact that all antennas should work together to nullify the field in one of the reception points. It is also noted that the electric current amplitudes of main antennas for the point where
Table 1. Results of the antennas feed current for ASK transmission, achieving isolation between the two observations points.
there is no transmission is small compared with the current amplitudes of the main antennas for the point where there is transmission.
Thus, all antennas work together and are responsible for the transmission of all signals. As each antenna carries a small portion of the information, the fading can be mitigated, due to the path diversity.
For illustration, Figure 7 shows the electric fields obtained when the desired signal is bit 1 at point 1 and null at point 2.
Figure 7 shows that the electric fields were generated according to the desired
signals, that is, a circular polarization in the plane, with an
amplitude equals to 4 × 10−3 V/m for the point of observation O1 (a), and an amplitude equals zero for the point of observation O2 (b).
In this way, any electric field is generated at the observation point, regardless of the desired field at another point, resulting in a transmission isolated by spatial separation or spatial multiplexing.
To observe how the signal is degraded when users move away from the observation point used in the design, consider other example using the same conditions, except for, which is approximately 3 km for 100 MHz. In order to obtain the same desired signals at O1 e O2, with the new current values obtained by (30), electric fields were simulated for different points
The time-average amplitude of the resulting field at the observation point i is calculated as follows:
In possession of the value, the error of time-average amplitude of the resulting field regarding the value of the desired field is calculated as:
Figure 7. Instantaneous electric fields resulting in points of observations O1 (a) and O2 (b), given in V/m, when the desired signals are bit 1 at O1 and null at O2, and the antenna feed current shown in Table 1.
Instead of calculating the difference of electric field amplitude at each instant of time, (32) was used in order not to consider the propagation delays between observations points. As the polarizations are in different planes, and the observa- tions points are separated by multiple wavelengths, the Equation (32) becomes more appropriate to the intended results. Table 2 shows these results.
In Table 2, the positive or negative sign of the value of the error indicates that the time-average amplitude of the field is increased or reduced. It is evident that, when the desired signals are not zero at the points of observations, the errors are minor compared to the lack of transmission. In both cases, the errors are small, though the observations points used in the design are about 300 m away from the receptions points. Figure 8 shows the resulting electric fields for one case.
Despite that the O1n and O2n points are approximately 300 m from the points O1
Table 2. Errors of time-average amplitude of the resulting electric field for the O1n and O2n points.
Figure 8. Instantaneous electric fields resulting at points of observations O1n (a) and O2n (b), given in V/m, when the desired signals are null at O1 and bit 0 at O2.
and O2, the amplitudes of the electric fields generated and the polarization planes almost did not suffer degradation.
Considering the small decoupling between the receiving antennas due to the difference between the planes of polarization, the received power will be slightly smaller. Consequently, the largest error obtained in the example will have a small reduction.
In fact, spatial regions whose deviations are small are formed around the observation point Oi, and the isolation is preserved. Thus, the reception point could move within that region without suffer interference. The region depends on the threshold detection to identify the signals.
Although the distance between the new observations points and those used in the design are large, it represents only a percentage of the distance from the observation points O1 and O2 to the origin, which equals 3.000 m. Thus, as these points move away from their main antennas, the greater will be the region of reception, where users can move without any significant errors. This is because the amplitude of the electric fields suffers smooth transitions between the two observation points.
Ignoring the effects of fading, the theoretical perfect isolation can be maintained if the system knows the locations of all the reception points. So with the new position of the observation point, the matrix of complex coefficients s could be updated, if necessary and new values of currents are obtained using (30). The sensitivity of the signals with respect to various distances from the target points, as well as case examples with pyramidal horn antennas, is discussed in depth in  . Furthermore, this approach can be of substantial interest when considering the fact that cellular systems are moving towards micro and phantom cells  .
5. Torus-Knots Polarizations and Fields
The theory and possible applications related to torus-knotted electromagnetic fields is discussed in   , which can be encountered in nature around celestial objects and are used for plasma confinement and with nuclear magnetic resonance devices. Here we only comment for completeness that antenna array configurations similar to the ones presented in the previous sections can also be employed to generate these distributions in a controlled and consistent fashion. Furthermore, we show that not only the torus-knotted fields but polarizations as well can be obtained with the procedure  . Figure 9 and Figure 10 show an example of torus-knots for an electromagnetic three-dimensional polarization distribution (as time varies). The synthesis procedure is herein omitted for simplicity, as well as the torus-knotted field distributions (at a given instant of time), and the reader is referred to  for more information.
This work extended substantially a procedure introduced by the authors in  to generate three-dimensional (3D) polarizations, where the electric field intensity
Figure 9. An illustrative example of torus-knotted polarization in the xy-plane, produced by an array of 3 ideal dipole antennas. Although not herein shown similar distributions for the electromagnetic fields at a given instant of time can also be obtained  .
Figure 10. An illustrative example of torus-knotted polarization (three-dimensional view of Figure 9). Although not herein shown similar distributions for the electromagnetic fields at a given instant of time can also be obtained  .
vector can be oriented in any direction in space. These polarizations were referred to as “3D” due to their distribution in space, although time is also inherent. It was shown analytically and numerically that it is possible to generate 3D polarizations with different shapes, as time varies, such as cylinders, spheres, and others. In contrast to  where only ideal dipoles were employed for generating 3D polarizations, here pyramidal horns were also used, not only for polarizations but also to generate controlled electromagnetic field distributions. In addition, the orthogonality between 3D polarizations was also mentioned. It is worth noting that for the scenarios considered in this study, there is no propagation, so the 3D polarization needs to be visualized as the polarization of the resulting electric field, not the wave. Alternate configurations that allow the 3D polarization to move within a desired path are currently under investigation (by determining different sets of feed currents for each desired location, for example), but the basic theory as presented here does not lose its generality. In fact, more accurate results can be obtained with the numerical modeling using commercial and/or in-house antenna codes, which will also allow determining the correct excitations to compensate for effects not accounted here, such as mutual coupling in smaller cells, near-field effects for non-ideal antennas, reflections, refractions, etc…
Although the main purpose was to introduce a mathematically consistent synthesis procedure for controlling the three-dimensionally distributed electromagnetic fields and polarizations, with applications yet to be adapted or discovered, an effort was made to illustrate that the technique can be used for digital transmission of information in wireless communication systems with spatial multiplexing. For each observation point, three antennas will be necessary. If more antennas are used, the fading will be smaller and the results will be closer to the desired. Volumetric regions are created with this array, and isolation could be preserved if the reception points are within these regions. The interference for spatial multiplexing can be further reduced if the method of electric field control is combined with the use of orthogonal 3D polarizations between close users, as in nano and phantom cellular systems  . Finally, another possible use of the theory was briefly mentioned for torus-knotted electromagnetic fields, as well as illustrated for three-dimensional polarization distributions. Torus-knotted fields can be found naturally around various celestial objects, as historically shown in Radioastronomy as well as possibly to be further investigated with advanced Radiotelescopes  , and can be of importance for plasma confinement and nuclear magnetic resonance devices   . Measurement procedures for three-dimensional electromagnetic fields in space (at a given instant of time) and polarizations (as time varies) are subject of further discussion, as for mathematical multi-dimensional results.
 Forenza, A. and Heath, R.W. (2006) Benefit of Pattern Diversity via Two-Element Array of Circular Patch Antennas in Indoor Clustered MIMO Channels. IEEE Transactions on Communications, 54, 943-954.
 Terada, M.A.B. (2003) Novel Beam Scanning Technique for Circularly Polarized Reflector Antennas. Proceedings of the IEEE International Symposium on Antennas and Propagation, Columbus, 22-27 June 2003, 770-773.
 Terada, M.A.B. (2004) Distributed Three-Dimensional Array Configuration for Volumetric-Controlled Field Coverage: A New Concept for High-Efficiency, Distance-Selective Wireless Local Area Network. URSI-B Symposium on Electromagnetic Theory (URSI_EMTS), Pisa, 23-27 May 2004, 975-977.
 Terada, M.A.B. (2004) Distance-Selective Wireless Network with Enhanced Isolation Characteristics: A New Concept Based on Three-Dimensional Wave Synthesis for Volumetric-Controlled Field Coverage. Proceedings of the 2004 IEEE International Symposium on Antennas and Propagation, Monterey, 20-25 June 2004, 2, 1736-1739.
 Pereira L.P.S. and Terada, M.A.B. (2011) New Method for Optimum Design of Pyramidal Horn Antennas. Journal of Microwaves, Optoelectronics and Electromagnetic Applications, 10, 266-277.
 Wanderley, P. and Terada, M. (2012) Assessment of the Applicability of the Ikegami Propagation Model in Modern Wireless Communication Scenarios. Journal of Electromagnetic Waves and Applications, 26, 1483-1491.
 Young, André, Terada, Marco A.B., de Villiers, Dirk I.L. and Davidson, David B. (2012) Assessment of the Sensitivity of the South African KAT-7 and MeerKAT/SKA Radio Telescope Reflector Antennas. IEEE 2012 International Conference on Electromagnetics in Advanced Applications, 1, 486-489.
 Arrayás, M. and Trueba, J.L. (2015) A Class of Non-Null Toroidal Electromagnetic Fields and Its Relation to the Model of Electromagnetic Knots. Journal of Physics A: Mathematical and Theoretical, 48, 1-14.
 Mathcad (2017)