Toraldo’s Composed Pupil: A Theoretical Analysis of the Near Field

: Over the years, there has been much speculation to understand whether (and how) it was possible to go below the diffraction limit. An advance in knowledge was achieved with the development of microwave techniques. In fact, more than fifty years after the publication of Toraldo’s article dealing with this topic, some experimental measurements in the range of microwaves confirmed the validity of his model. Since some measurements were performed in the region of near field, while Toraldo’s model refers to the far field, the need for a theoretical analysis in the framework of the Fresnel optics arose. The main goal of the present paper is to describe the problem of propagation in the near field (Fresnel optics) by using the same theoretical model already proposed by Toraldo. In order to test the validity of this new approach, the theoretical model has been compared with the FEKO simulation. The comparison of the theoretical model with the FEKO simulation in the far field for an open pupil (an open circular aperture) shows perfect agreement, as expected. We will demonstrate that there is also good agreement in the near field, although it is limited to the region around the main lobe, which is usually the region of main physical interest. Moving away from the main lobe, namely away from the optical axis, the agreement becomes less significant.


Introduction
The topic of super-resolution is undoubtedly one of the main topics studied in the field of physics because of its implications in theoretical physics (related to Heisenberg's uncertainty principle), as well as its possible practical applications.
After the pioneering paper of Schelkunoff [1] (in which it was demonstrated that the radiation pattern of a linear end-fire array of a given number of elements can be narrowed by a factor of more than two by carefully selecting the phase distribution of the currents relative to the several elements), in 1952, Giuliano Toraldo di Francia presented a theoretical study in which he demonstrated the possibility to increase the resolving power of an optical instrument with the use of special types of pupils (composed by concentric circular coronae), provided that specific requirements of the phase and amplitude of the single corona were satisfied [2].
Overcoming the diffraction limit is one of the main challenges that physicists have faced over the last hundred years.There has been much discussion as to whether this limit could be exceeded or not, and it was common opinion that, like the uncertainty principle in quantum mechanics, the diffraction limit seemed to be absolutely unsurmountable.
For this reason, Toraldo's model could not be ignored and found its first experimental confirmation in the range of microwaves [3,4] in 2003 and 2004.More than ten years later, some experimental measurements were performed at a frequency of 20 GHz [5], namely still in the range of microwaves but at a frequency different with respect to that of Refs.[3,4].The interest in the 20 GHz measurement stems from the fact that it belongs to one of the most common wavebands used in radio astronomy [5].As in the case of Refs.[3,4], measurements at 20 GHz also deal with the diffraction field of a composed pupil which a plane wave affects.Since those measures were performed in the near field and Toraldo's model refers to the far field (see Appendix A), the need to have a theoretical model able to interpret those results arose.
Many studies have been dedicated to the evaluation of the near field in different physical systems, from antenna array technology to sensing and communications technology in view of new sixth-generation (6G) networks [6,7], to the characterization of bright and dark foci in the diffraction field [8], and to topics more related to quantum mechanics such as, theoretical statistics for systems describing the near field for photon antibunching and entanglement [9].The literature relating to near-field propagation is extremely rich as regards super-resolution in the construction of images [10][11][12][13][14][15].However, as far as we know, very little can be found in the literature regarding the super-resolution effects of electromagnetic propagation in the range of microwaves.
In this paper, we present, in Section 2, a theoretical model in order to analyze the diffraction in the near field for an open pupil as well as for a composed one.Then, in Section 3, we compare the theoretical results with the FEKO simulation; in Section 4, we present our conclusions, and finally, in Appendix A, we report a brief description of Toraldo's model.

Theoretical Analysis
Let us consider a plane wave normally incident on an opaque diaphragm (theoretically of infinite dimensions) with a generic aperture C whose extension is much smaller than the size of the diaphragm (see Figure 1).
Optics 2024, 5, FOR PEER REVIEW 2 GHz [5], namely still in the range of microwaves but at a frequency different with respect to that of Refs.[3,4].The interest in the 20 GHz measurement stems from the fact that it belongs to one of the most common wavebands used in radio astronomy [5].As in the case of Refs.[3,4], measurements at 20 GHz also deal with the diffraction field of a composed pupil which a plane wave affects.Since those measures were performed in the near field and Toraldo's model refers to the far field (see Appendix A), the need to have a theoretical model able to interpret those results arose.
Many studies have been dedicated to the evaluation of the near field in different physical systems, from antenna array technology to sensing and communications technology in view of new sixth-generation (6G) networks [6,7], to the characterization of bright and dark foci in the diffraction field [8], and to topics more related to quantum mechanics such as, theoretical statistics for systems describing the near field for photon antibunching and entanglement [9].The literature relating to near-field propagation is extremely rich as regards super-resolution in the construction of images [10][11][12][13][14][15].However, as far as we know, very li le can be found in the literature regarding the superresolution effects of electromagnetic propagation in the range of microwaves.
In this paper, we present, in Section 2, a theoretical model in order to analyze the diffraction in the near field for an open pupil as well as for a composed one.Then, in Section 3, we compare the theoretical results with the FEKO simulation; in Section 4, we present our conclusions, and finally, in Appendix A, we report a brief description of Toraldo's model.

Theoretical Analysis
Let us consider a plane wave normally incident on an opaque diaphragm (theoretically of infinite dimensions) with a generic aperture  whose extension is much smaller than the size of the diaphragm (see Figure 1).Under the condition that the distance between the diaphragm and the screen (or, more generally, the observation region) is much greater with respect to the dimension of the aperture , the diffracted field at a generic point (, , ) in the near region, and in the scalar approximation, can be wri en as [16]  Under the condition that the distance between the diaphragm and the screen (or, more generally, the observation region) is much greater with respect to the dimension of the aperture C, the diffracted field at a generic point P(x, y, z) in the near region, and in the scalar approximation, can be written as [16] where λ is the wavelength, k = 2π/λ is the wave number, A is the amplitude of the incident wave, and x ′ and y ′ are the Cartesian coordinates of a point over the diaphragm (Equation (1) can be generalized also for a diaphragm with more than one aperture by simply adding other integrals for additional apertures).We wish to recall that in the case of a circular aperture, the near region can be defined as the region of radius R ≲ 2D 2 /λ , where D is the diameter of the source and λ is the wavelength.In the case in which the aperture is characterized by rotational symmetry, it is convenient to use polar coordinates.Let us denote with ρ ′ and θ ′ the coordinates of a point over the diaphragm, and with ρ and θ, the ones of a point on the screen, such as (2) therefore, we have the following: In Equation (1), we can change the integration variables from dx ′ , dy ′ to dρ ′ , dθ ′ by using the Jacobian determinant, and we obtain By substituting previous equations into Equation (1), the diffraction field for a circular aperture becomes with D being the diameter of the aperture.By recalling the integral representation of the Bessel function of first kind, that is and by considering that the cosine function is an even function, we obtain where R = z + ρ 2 /2z (for ρ ≪ z, R is approximately equal to the distance between the center of the diaphragm and a given point on the screen) and J 0 is the zero-order Bessel function of first kind.Equation ( 3) can be integrated by parts by using the following well-known formula: and by choosing the integration of Equation ( 3) gives Optics 2024, 5

322
By evaluating the term in the square brackets, we obtain We can now continue by solving the integral in Equation ( 5), still by using the integration by parts, whose solution gives rise to an integral which can be solved again by the integration by parts, and so on.Generally speaking, we can see that the solution of the integral like that of Equation ( 3) gives rise to a series of the type In the present case, a = k/2z, b = kρ/z , and ε = D/2 ; therefore, the final result is If we limit ourselves to the first term E 1 inside brackets in Equation (7), namely we obtain an amplitude |E 1 | which has the same behavior as the far field A f f (as expected), which can be written as (see Appendix A) where x ≡ ρ, and where we have made the substitution sin θ = x/ √ x 2 + z 2 ≃ x/z (z ≫ x).This means that the near field is described by an equation in which the first term is similar to that of the far field (the only difference lies in the numerical factor), which is dominant for z ≫ D, while the other terms constitute a correction that becomes more and more important when the value of z approaches D. The results are shown in Figure 2: in (a), we can see that the first term of Equation (7) (dotted line) is exactly equal to the value of the far field (solid line), while in (b), the total value of Equation ( 7) (dotted line) shows a clear discrepancy between the near and far field when z approaches D.
Once the diffracted field in the near field is known for a circular aperture, namely for an open pupil, it is easy to estimate the diffracted field for a circular corona which can be evaluated by the difference between the transmitted and the reflected fields.If we indicate with D 1 and D 2 the diameters of the internal and external circumferences, respectively (see Figure 3), then the diffracted amplitude A c for the circular corona is given by (still for an incident plane wave) where E is given by Equation (7).
The material surrounding the pupil is assumed to be a perfect electric conductor (PEC), thus totally reflecting the electromagnetic radiation (blue color in Figure 3).In the microwave range, the material through which the field is transmitted (white color in Figure 3) is constituted, for instance, by a dielectric whose thickness is selected to produce the dephasing required by the theoretical model in order to obtain super-resolution.7) till the 7th term of the series.As expected, the first term of the series of ( 7), which refers to the near field, is equal to that of the far field (when normalized), while the total near field is equal to the far field only for values of  sufficiently greater than the diameter of the aperture.When the value of  approaches , the discrepancy between the near and far field becomes more and more evident.
Once the diffracted field in the near field is known for a circular aperture, namely for an open pupil, it is easy to estimate the diffracted field for a circular corona which can be evaluated by the difference between the transmi ed and the reflected fields.If we indicate with  and  the diameters of the internal and external circumferences, respectively (see Figure 3), then the diffracted amplitude  for the circular corona is given by (still for an incident plane wave) where  is given by Equation ( 7).The material surrounding the pupil is assumed to be a perfect electric conductor (PEC), thus totally reflecting the electromagnetic radiation (blue color in Figure 3).In the microwave range, the material through which the field is transmi ed (white color in Figure 3) is constituted, for instance, by a dielectric whose thickness is selected to produce the dephasing required by the theoretical model in order to obtain super-resolution.7) till the 7th term of the series.As expected, the first term of the series of (7), which refers to the near field, is equal to that of the far field (when normalized), while the total near field is equal to the far field only for values of z sufficiently greater than the diameter of the aperture.When the value of z approaches D, the discrepancy between the near and far field becomes more and more evident.7) till the 7th term of the series.As expected, the first term of the series of (7), which refers to the near field, is equal to that of the far field (when normalized), while the total near field is equal to the far field only for values of  sufficiently greater than the diameter of the aperture.When the value of  approaches , the discrepancy between the near and far field becomes more and more evident.
Once the diffracted field in the near field is known for a circular aperture, namely for an open pupil, it is easy to estimate the diffracted field for a circular corona which can be evaluated by the difference between the transmi ed and the reflected fields.If we indicate with  and  the diameters of the internal and external circumferences, respectively (see Figure 3), then the diffracted amplitude  for the circular corona is given by (still for an incident plane wave) where  is given by Equation ( 7).The material surrounding the pupil is assumed to be a perfect electric conductor (PEC), thus totally reflecting the electromagnetic radiation (blue color in Figure 3).In the microwave range, the material through which the field is transmi ed (white color in Figure 3) is constituted, for instance, by a dielectric whose thickness is selected to produce the dephasing required by the theoretical model in order to obtain super-resolution.The result of Equation ( 10) can be generalized for any system constituted by N circular coronae plus the open pupil in the center.By denoting with E n the contribution of each corona, and by E 0 the contribution due to the central pupil (coincident with E(D 2 ) in Equation ( 10)), the total diffracted amplitude A t is given by where the plus sign refers to the transmitted field, and the negative one to the reflected field.

Theoretical Model and FEKO Simulation Results
In order to evaluate the reliability of the theoretical model, we can compare it with the results of the commercial software FEKO © [17], a comprehensive EM simulation software tool used for the electromagnetic field analysis of 3D structures.
Let us compare the normalized amplitude evaluated by the theoretical model, Equation ( 7), and by FEKO, for an open pupil of diameter D = 9 cm, at a frequency of 20 GHz, and for a distance z = 18 cm (Figure 4a).The choice of these values is due to the fact that some experimental measurements have been performed precisely for these values [5].
the results of the commercial software FEKO © [17], a comprehensive EM simulation software tool used for the electromagnetic field analysis of 3 structures.
Let us compare the normalized amplitude evaluated by the theoretical model, Equation (7), and by FEKO, for an open pupil of diameter  = 9 cm, at a frequency of 20 GHz, and for a distance  = 18 cm (Figure 4a).The choice of these values is due to the fact that some experimental measurements have been performed precisely for these values [5].As we can see in Figure 5, the results are in good agreement almost everywhere, showing only an opposition of the phase in the tail, around the value of 8 cm, i.e., very close to the edge of the pupil.This modest discrepancy may be due (apart from the different way to evaluate the field) to the fact that the value of  (=18 cm) cannot be considered much greater with respect to  (=9 cm), a condition required by the theoretical model described in the previous section.As we can see in Figure 5, the results are in good agreement almost everywhere, showing only an opposition of the phase in the tail, around the value of 8 cm, i.e., very close to the edge of the pupil.This modest discrepancy may be due (apart from the different way to evaluate the field) to the fact that the value of z (=18 cm) cannot be considered much greater with respect to D (=9 cm), a condition required by the theoretical model described in the previous section.7)), while the do ed line is the result obtained by the FEKO simulation.
As a further step, let us consider a composed pupil (hereafter named TP3) constituted by two coronae plus the central pupil, as shown in Figure 4b.In this case, if we denote with  and  the contribution of the transmi ed and reflected fields by the external and internal coronae, respectively (each one evaluated through Equation ( 10)), and with  the contribution due to the central pupil (Equation ( 7)), the total transmi ed amplitude  can be wri en as As required by Toraldo's model, in order to have super-resolution, each corona must satisfy specific conditions for the phase and amplitude (see Appendix A).However, the condition of the phase alone is a necessary and sufficient condition in order to have superresolution; thus, in the case of TP3, the previous condition is satisfied by adding a  factor to the phase of the innermost corona.
In Figure 6, we report the results of the theoretical model for a TP3 composed pupil together with the results of the FEKO simulation; we can see that both the TP3 theoretical model and the FEKO simulation clearly show a narrower main lobe compared to the open pupil.We can also see that there is rather good agreement between the theoretical model Other parameter values are the same as in Figure 2. The solid red line is the theoretical model (Equation ( 7)), while the dotted line is the result obtained by the FEKO simulation.
As a further step, let us consider a composed pupil (hereafter named TP3) constituted by two coronae plus the central pupil, as shown in Figure 4b.In this case, if we denote with E 1 and E 2 the contribution of the transmitted and reflected fields by the external and internal coronae, respectively (each one evaluated through Equation ( 10)), and with E 0 the contribution due to the central pupil (Equation ( 7)), the total transmitted amplitude A t can be written as As required by Toraldo's model, in order to have super-resolution, each corona must satisfy specific conditions for the phase and amplitude (see Appendix A).However, the condition of the phase alone is a necessary and sufficient condition in order to have superresolution; thus, in the case of TP3, the previous condition is satisfied by adding a π factor to the phase of the innermost corona.
In Figure 6, we report the results of the theoretical model for a TP3 composed pupil together with the results of the FEKO simulation; we can see that both the TP3 theoretical model and the FEKO simulation clearly show a narrower main lobe compared to the open pupil.We can also see that there is rather good agreement between the theoretical model and FEKO simulation, especially in the region of physical interest, namely the main lobe.The discrepancy in the area beyond 6 cm may be due to the fact that, as previously discussed, Equation ( 7) has been obtained under the condition that the distance z is much greater with respect to the dimension of the aperture.Actually, it is not easy to establish when this condition is fully satisfied, and this may be the reason why there is not perfect agreement between the two approaches.Moreover, we have to recall that Toraldo's model in the near field is described by an infinite series, which cannot be totally evaluated, and the FEKO simulator is not able to reproduce the exact theoretical model presented here (we will discuss this issue in our conclusions).In Figure 7, we show the diffraction amplitude for a TP4 pupil, i.e., composed of three coronae plus the central open aperture (Figure 7a), which again shows a narrower main lobe.The agreement between the theoretical model and the FEKO simulation is still good, except for x ≳ 6 cm, where, except for a slight difference in the amplitude, there is a shift in the position of the maxima and minima.As a final test, we have also performed a comparison between the analytical model and the FEKO simulations in the case of a five-coronae TP (or TP5) in the far field (Figure 8).We note again reasonable agreement in the region of the main lobe, while the sidelobes In Figure 7, we show the diffraction amplitude for a TP4 pupil, i.e., composed of three coronae plus the central open aperture (Figure 7a), which again shows a narrower main lobe.The agreement between the theoretical model and the FEKO simulation is still good, except for x ≳ 6 cm, where, except for a slight difference in the amplitude, there is a shift in the position of the maxima and minima.In Figure 7, we show the diffraction amplitude for a TP4 pupil, i.e., composed of three coronae plus the central open aperture (Figure 7a), which again shows a narrower main lobe.The agreement between the theoretical model and the FEKO simulation is still good, except for x ≳ 6 cm, where, except for a slight difference in the amplitude, there is a shift in the position of the maxima and minima.As a final test, we have also performed a comparison between the analytical model and the FEKO simulations in the case of a five-coronae TP (or TP5) in the far field (Figure 8).We note again reasonable agreement in the region of the main lobe, while the sidelobes appear quite different between the theoretical model and the numerical simulation.The As a final test, we have also performed a comparison between the analytical model and the FEKO simulations in the case of a five-coronae TP (or TP5) in the far field (Figure 8).We note again reasonable agreement in the region of the main lobe, while the sidelobes appear quite different between the theoretical model and the numerical simulation.The possible reasons for this disagreement will be discussed in our conclusions.As previously mentioned, simulations at 20 GHz were chosen because of its applications in radio astronomy.In fact, the work presented here is part of a larger project, named "Pupille Toraldo" (or PUTO for short), which is devoted to the designing, testing, and implementing of a super-resolution optical system on a radio telescope.As part of this project, and based on the theoretical model and our FEKO simulations, we conducted more comprehensive laboratory tests at the selected frequency of 20 GHz [5].These tests clearly showed that the super-resolution effect was achieved and also showed good agreement with what was expected.After confirming the presence of super-resolving effects, we designed an optical module based on a two-lens collimator placed after the Cassegrain focus of the telescope and before the receiver dewar.The first lens of the collimator generates an image of the primary reflector which is then brought into a subsequent focus by the second lens.Toraldo's pupil is placed at the image of the entrance pupil where it can modify the incident wavefront [18].
After further tests (inside an anechoic chamber) of the two-lens collimator, which showed results quite consistent with the FEKO simulation, the Toraldo's pupil optical module was successively mounted on a small Cassegrain antenna 125 cm in diameter [19] with the same focal ratio (≃3) as the Medicina 32 m radio telescope [20].Electromagnetic simulations and preliminary field tests of the antenna plus optical module have shown that super-resolution can be achieved with this optical arrangement.An improved Toraldo's pupil optical module is now being designed for applications with the Istituto Nazionale di Astrofisica (INAF) radio telescopes.

Conclusions
The concept of super-resolution refers to various methods for improving the angular resolution of an optical imaging system beyond the classical diffraction limit.Variable transmi ance pupils constitute an a ractive method to achieve super-resolution, and Toraldo's pupils represent the simplest example of such filters.
We have analyzed the field by Toraldo's pupils and found a formula for the complex amplitude in the near-field region.We have then compared our result with FEKO numerical simulations at 20 GHz for a few Toraldo's pupils.We found reasonable agreement between the FEKO results and our theoretical model, with some deviations outside the main lobe.As previously mentioned, simulations at 20 GHz were chosen because of its applications in radio astronomy.In fact, the work presented here is part of a larger project, named "Pupille Toraldo" (or PUTO for short), which is devoted the designing, testing, and implementing of a super-resolution optical system on a radio telescope.As part of this project, and based on the theoretical model and our FEKO simulations, we conducted more comprehensive laboratory tests at the selected frequency of 20 GHz [5].These tests clearly showed that the super-resolution effect was achieved and also showed good agreement with what was expected.After confirming the presence of super-resolving effects, we designed an optical module based on a two-lens collimator placed after the Cassegrain focus of the telescope and before the receiver dewar.The first lens of the collimator generates an image of the primary reflector which is then brought into a subsequent focus by the second lens.Toraldo's pupil is placed at the image of the entrance pupil where it can modify the incident wavefront [18].
After further tests (inside an anechoic chamber) of the two-lens collimator, which showed results quite consistent with the FEKO simulation, the Toraldo's pupil optical module was successively mounted on a small Cassegrain antenna 125 cm in diameter [19] with the same focal ratio (≃3) as the Medicina 32 m radio telescope [20].Electromagnetic simulations and preliminary field tests of the antenna plus optical module have shown that super-resolution can be achieved with this optical arrangement.An improved Toraldo's pupil optical module is now being designed for applications with the Istituto Nazionale di Astrofisica (INAF) radio telescopes.

Conclusions
The concept of super-resolution refers to various methods for improving the angular resolution of an optical imaging system beyond the classical diffraction limit.Variable transmittance pupils constitute an attractive method to achieve super-resolution, and Toraldo's pupils represent the simplest example of such filters.
We have analyzed the field diffracted by Toraldo's pupils and found a formula for the complex amplitude in the near-field region.We have then compared our result with FEKO numerical simulations at 20 GHz for a few Toraldo's pupils.We found reasonable agreement between the FEKO results and our theoretical model, with some deviations outside the main lobe.
On the one hand, the theoretical model suffers some approximations and therefore has some limitations in the near field (but in the far field, it proved to work well).On the other hand, the FEKO simulator is not totally reliable because of difficulties in modeling complex systems such as those represented by Toraldo's pupils.More specifically, the FEKO solver is based on the full-wave technique of the Method of Moment, which solves Maxwell's equations numerically.Therefore, FEKO provides an exact solution of the EM problem, but with a finite accuracy determined by the selected meshing of the set-up geometry.
In the case under consideration, the diaphragm consists of a finite screen illuminated by an incident plane wave, therefore producing significant diffraction from the edges of the screen.This deviates from the theoretical model where external diffraction does not exist thanks to the assumption of an infinite screen.Using a directive horn antenna instead of a plane wave has been tested to reduce the illumination of the edges of the plane, but it caused a non-uniform phase in the diaphragm plane.Another difference between the theoretical model and FEKO is in the dielectric thickness, which is infinitesimal and finite in the two approaches, respectively.Finally, in the FEKO model, the electric field transported by the EM wave features two orthogonal components, namely the main component (copolar) and the unwanted component (cross-polar), while in the theoretical model, the latter component is neglected.
Therefore, the theoretical model results and the FEKO results are expected to show some discrepancies, although both approaches are able to describe quite well the propagation in the case of Toraldo's pupils.For an ideal system (namely infinite screen, an ideal incident plane wave, zero thickness for materials composing pupils, etc.), the theoretical model is probably more reliable.On the contrary, for a real system, in which previous conditions are not satisfied, FEKO simulation represents the better tool to be utilized.
The convenience of the theoretical model is that it can be used to predict the profile of the near field generated by Toraldo's pupil in order to select the best geometry for a specific application.For optical systems more complex compared to the simple Toraldo's pupil, the theoretical model can be useful for speeding up the optimization procedure, which is not always easy to implement and often takes a long time with electromagnetic software tools.

Figure 1 .
Figure 1.A schematic illustration of an electromagnetic diffraction process for a generic aperture .The field is detected over a screen positioned at a given distance from the diaphragm.

Figure 1 .
Figure 1.A schematic illustration of an electromagnetic diffraction process for a generic aperture C. The field is detected over a screen positioned at a given distance from the diaphragm.

Figure 2 .
Figure 2. The normalized intensity at a point  on the screen as a function of the radial coordinate for diffraction by an open pupil (a circular aperture) and for different values of z, for  = 9 cm and  = 1.5 cm.The do ed lines refer to the far field, while the solid lines refer to the near field evaluated by (a) the first term of Equation (7); (b) all of Equation (7) till the 7th term of the series.As expected, the first term of the series of (7), which refers to the near field, is equal to that of the far field (when normalized), while the total near field is equal to the far field only for values of  sufficiently greater than the diameter of the aperture.When the value of  approaches , the discrepancy between the near and far field becomes more and more evident.

Figure 3 .
Figure 3.A drawing of a single corona whose external and internal diameters are  and  , respectively.

Figure 2 .
Figure 2. The normalized intensity at a point P on the screen as a function of the radial coordinate for diffraction by an open pupil (a circular aperture) and for different values of z, for D = 9 cm and λ = 1.5 cm.The dotted lines refer to the far field, while the solid lines refer to the near field evaluated by (a) the first term of Equation (7); (b) all of Equation (7) till the 7th term of the series.As expected, the first term of the series of(7), which refers to the near field, is equal to that of the far field (when normalized), while the total near field is equal to the far field only for values of z sufficiently greater than the diameter of the aperture.When the value of z approaches D, the discrepancy between the near and far field becomes more and more evident.

Figure 2 .
Figure 2. The normalized intensity at a point  on the screen as a function of the radial coordinate for diffraction by an open pupil (a circular aperture) and for different values of z, for  = 9 cm and  = 1.5 cm.The do ed lines refer to the far field, while the solid lines refer to the near field evaluated by (a) the first term of Equation (7); (b) all of Equation (7) till the 7th term of the series.As expected, the first term of the series of(7), which refers to the near field, is equal to that of the far field (when normalized), while the total near field is equal to the far field only for values of  sufficiently greater than the diameter of the aperture.When the value of  approaches , the discrepancy between the near and far field becomes more and more evident.

Figure 3 .
Figure 3.A drawing of a single corona whose external and internal diameters are  and  , respectively.

Figure 3 .
Figure 3.A drawing of a single corona whose external and internal diameters are D 2 and D 1 , respectively.

Figure 4 .
Figure 4.A drawing of an open pupil (a) and a system constituted by two coronae plus the central aperture (b).

Figure 4 .
Figure 4.A drawing of an open pupil (a) and a system constituted by two coronae plus the central aperture (b).

7 Figure 5 .
Figure 5.The normalized amplitude diffracted by an open pupil: the field is evaluated at  = 18 cm.Other parameter values are the same as in Figure 2. The solid red line is the theoretical model (Equation (7)), while the do ed line is the result obtained by the FEKO simulation.

Figure 5 .
Figure 5.The normalized amplitude diffracted by an open pupil: the field is evaluated at z = 18 cm.Other parameter values are the same as in Figure2.The solid red line is the theoretical model (Equation (7)), while the dotted line is the result obtained by the FEKO simulation.

8 Figure 6 .
Figure 6.The normalized amplitude, as a function of the transversal coordinate , evaluated by Equation (7) ( ≡ , solid black line) and by the FEKO simulation (do ed magenta line).The dimension of the composed corona is the same as in Figure 4b.For comparison, the amplitude of the open pupil is also shown (solid red line).The parameter values are the same as in Figure 2.

Figure 7 .
Figure 7. Composed pupil TP4 (a) and relative diffraction amplitude (b).Symbols and parameter values are same as in Figure 6.

Figure 6 .
Figure 6.The normalized amplitude, as a function of the transversal coordinate §, evaluated by Equation (7) (ρ ≡ x, solid black line) and by the FEKO simulation (dotted magenta line).The dimension of the composed corona is the same as in Figure 4b.For comparison, the amplitude of the open pupil is also shown (solid red line).The parameter values are the same as in Figure 2.

8 Figure 6 .
Figure 6.The normalized amplitude, as a function of the transversal coordinate , evaluated by Equation (7) ( ≡ , solid black line) and by the FEKO simulation (do ed magenta line).The dimension of the composed corona is the same as in Figure 4b.For comparison, the amplitude of the open pupil is also shown (solid red line).The parameter values are the same as in Figure 2.

Figure 7 .
Figure 7. Composed pupil TP4 (a) and relative diffraction amplitude (b).Symbols and parameter values are same as in Figure 6.

Figure 7 .
Figure 7. Composed pupil TP4 (a) and relative diffraction amplitude (b).Symbols and parameter values are same as in Figure 6.

Figure 8 .
Figure 8.(a) A composed TP5 pupil and (b) an evaluation of the relative amplitude in the far field (Fraunhofer region) as a function of the diffraction angle .The red do ed line represents the FEKO simulation, while the solid blue one is the result of the theoretical model evaluated by Equation (7) with  = 5000 cm.Starting from the central pupil, the diameters are 1.5, 3.6, 7.2, 9, and 16 cm.

Figure 8 .
Figure 8.(a) A composed TP5 pupil and (b) an evaluation of the relative amplitude in the far field (Fraunhofer region) as a function of the diffraction angle θ.The red dotted line represents the FEKO simulation, while the solid blue one is the result of the theoretical model evaluated by Equation (7) with z = 5000 cm.Starting from the central pupil, the diameters are 1.5, 3.6, 7.2, 9, and 16 cm.