Next Article in Journal
Enhanced Yb:YAG Active Mirrors for High Power Laser Amplifiers
Next Article in Special Issue
Robust Topological Edge States in C6 Photonic Crystals
Previous Article in Journal
Photonic-Assisted Microwave Frequency Measurement Using High Q-Factor Microdisk with High Accuracy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Focusing of Radially Polarized Electromagnetic Waves by a Parabolic Mirror

by
Zerihun Tadele Godana
1,2,*,
János Hebling
1,2 and
László Pálfalvi
1,3
1
Institute of Physics, University of Pécs, H-7624 Pécs, Hungary
2
Szentágothai Research Centre, University of Pécs, H-7624 Pécs, Hungary
3
ELKH—PTE High-Field Terahertz Research Group, H-7624 Pécs, Hungary
*
Author to whom correspondence should be addressed.
Photonics 2023, 10(7), 848; https://doi.org/10.3390/photonics10070848
Submission received: 29 June 2023 / Revised: 17 July 2023 / Accepted: 18 July 2023 / Published: 21 July 2023
(This article belongs to the Special Issue Photonics: Theory, Technologies and Applications)

Abstract

:
It is well-known that a strong longitudinal electric field and a small spot size are observed when radially polarized beams are tightly focused using a high numerical aperture parabolic mirror. The longitudinal electric field component can accelerate electrons along the propagation axis at high intensities in the focal region, which opens an application in particle acceleration. In this paper, we present a rigorous derivation of the electric field obtained when a radially polarized, monochromatic, flat-top beam is focused by a parabolic mirror. The formulae were deduced from the Stratton–Chu integral known from vector diffraction theory. We examined the influence of the focusing parameters on the distribution of both the longitudinal and radial electric field components. In the small numerical aperture and short wavelength regimes, excellent agreement was found with the results obtained from the Rayleigh–Sommerfeld formula. The calculation method can be adapted for various beam types and for electromagnetic pulses as well.

1. Introduction

The paraboloid mirror can focus light nearly within a 4π solid angle [1,2]. Because of this special characteristic, researchers have been strongly interested in parabolic mirrors over the past several decades. The vector field nature of the light becomes essential in accurate description of a nonparaxial beam, when a beam is tightly focused by a parabolic mirror. The high-intensity laser community, which makes major efforts to achieve the highest laser intensities, is strongly interested in examinations of the vector field focusing characteristics [3]. By tightly focusing the beam in a vacuum using an off-axis parabolic (OAP) mirror, it is possible to attain extremely high intensity, which makes this of interest for laser-based particle acceleration [4]. For such a high-intensity, tightly focused electromagnetic field, a detailed description of the focused field is necessary in order to precisely identify the motion of charged particles [5]. Using a high numerical aperture parabolic mirror with a radially polarized beam is ideal for achieving a small focal spot size and strong longitudinal electric field, which opens an application possibility in particle acceleration [6,7].
Beginning with a study by Ignatovsky in 1920, a detailed diffraction theory of focused light from parabolic mirrors has been developed over the period of nearly a century. Ignatowsky transformed Maxwell’s equations to the parabolic coordinate system, set the boundary values, and then used these boundary values to solve Maxwell’s equations [8]. Richards and Wolf provided a different theoretical approach in which strongly focused beams were precisely defined in terms of the field distribution of the collimated input beam at the entry pupil of the focusing apparatus [9]. There are different approaches for evaluating tightly focused beams based on the Stratton–Chu formulation of Green’s theorem, which incorporates both the electric and magnetic fields after substituting values from Maxwell’s equations [10].
In 2000, Varga and Török derived a solution to the problem of a linearly polarized electromagnetic wave focused by a parabolic mirror using the Stratton–Chu integral by solving a boundary-value problem. They demonstrated that a segment of the paraboloid produces an intensity distribution that is similar to that obtained from a far-field-type approximation in the region of the focus [11]. In 2019, Xiahui Zegn and Xiya Chen proposed a precise analytical technique for the description of the vector electromagnetic fields created by an on- and off-axis parabolic mirror for circular and square incident beams based on the Stratton–Chu integral [3]. In 2021, the same group demonstrated an optimization process of OAP geometry to maximize the focused peak intensity based on precise knowledge of the tightly focusing properties of OAPs by employing the Stratton–Chu vector diffraction integrals and physical optics approximations and obtained the optimum configuration scale rule, allowing for the maximum peak intensity [12]. Endale et al. in 2022 examined focusing of radially polarized light using a Gaussian laser beam near its focal plane based on the Richards–Wolf diffraction method at various numerical apertures and Gaussian beam radii. They demonstrated that the longitudinal component becomes predominant at a high numerical aperture and high radius [13]. Recently, a more detailed study of the application of the vector field focusing properties of electromagnetic fields by a parabolic mirror based on the Stratton–Chu integral formalism was reported [14,15,16].
To the best of our knowledge, there is no literature that provides a detailed theoretical and analytical study of the vector field focusing properties of radially polarized beams by on- and off-axis parabolic mirrors based on the Stratton–Chu integral representation. In this paper, starting from the Stratton–Chu integral, we first derived general formulae to be used when a radially polarized electromagnetic plane wave is focused by a parabolic mirror. After validation, using these formulae we determined the electric field for various focusing conditions.

2. General Formulae of the Parabolic Mirror

Consider a parabolic mirror having its focus F at the origin O of a Cartesian coordinate system and a rotation axis that coincides with the z-axis, as shown in Figure 1. The paraboloid equation is given by:
z s = x s 2 + y s 2 4 f f ,
where x s , y s and z s are the coordinates of an arbitrary point S of the paraboloid and f is the focal length of the parabolic mirror. Let ρ s be the distance between point S and the focus F. Equation (1) becomes
ρ s = z s + 2 f
The unit vector pointing inward normal to the surface of the paraboloid is written as:
n x s , y s = 1 2 f x s e x + y s e y + e z 1 + x s 2   +   y s 2 4 f 2 .
Using the spherical coordinates, Equation (2) becomes
ρ s = 2 f 1 cos θ s ,
and the components of the unit vector n x s , y s are:
n x = sin θ s cos s 2 1 cos θ s 1 / 2 , n y = sin θ s sin s 2 1 cos θ s 1 / 2 , n z = 1 cos θ s 2 1 / 2 ,
where θ s is the polar, and s is the azimuthal angle as shown in Figure 1.

3. Rigorous Diffraction Theory of Radially Polarized Waves Focused by a Parabolic Mirror

First, we consider how the electric and magnetic fields are reflected from the surface of the parabolic mirror having perfect reflectance. Let us use E i and H i to represent the incident electric and magnetic fields, respectively. Based on the electromagnetic boundary conditions, upon reflection the tangential component of the electric field, E r , t , and the normal component of the magnetic field, H r , n , change sign, but the normal components of the electric field, E r , n , and the tangential components of the magnetic field, H r , t , remain unchanged [11]. Hence, the reflected fields are written as:
E r , n = E i , n = n E i · n , H r , t = H i , t = H i H i , n , E r , t = E i , t = E i E i , n , H r , n = H i , n = n H i · n ,
where E i , n and H i , n are the normal, and E i , t and H i , t are the tangential components of the incident electric and magnetic fields, respectively. Therefore, the total reflected fields are expressed as:
E r = E r , n + E r , t = 2 n E i · n E i , H r = H r , n + H r , t = H i 2 n H i · n .
The total field is given as the sum of the incident and reflected fields:
E = E i + E r = 2 n E i · n , H = H i + H r = 2 H i 2 n H i · n .
To derive an expression for the electric field focused by a segment of a parabolic mirror, we suppose that the segment of the paraboloid is limited by the apex V and a plane perpendicular to the optical axis positioned at z ≤ 0, as shown in Figure 2. The focusing angle δ, for which π/2 ≤ δ < π, determines the segment of the paraboloid [11].
We use the Stratton–Chu integral formula to derive an expression for the electric field upon focusing by a segment of a parabolic mirror [10].
E P = 1 4 π S i k n × H G + n × E × G + n · E G d A + 1 4 π i k C G H · d s ,
with
G u = e x p ( i k u ) u ,
where k = 2π/λ is wave number, λ is the wavelength of the incident beam and u is the distance between a point S x s , y s , z s on the surface of the paraboloid and the observation point P x , y , z . Hence,
u = x s x 2 + y s y 2 + z s z 2 1 / 2 = x 2 + y 2 + z 2 1 / 2 .
The first integral in Equation (9) represents the surface term, whereas the second integral represents the contour term, respectively.
G should be calculated at the points S x s , y s , z s of the paraboloid [11]. It can be expressed as:
G u = i k 1 1 i k u G ( u ) u x e x + y e y + z e z .
The surface element of the paraboloid dA is given by:
d A = 1 + z S x S 2 + z S y S 2 1 / 2 d x S d y S = ρ s 2 n z sin θ s d θ s d s .
To determine the electric field created upon focusing by a segment of the parabolic mirror, we assume that the radially polarized incident monochromatic plane wave propagates in the negative z direction. Considering the cylindrical coordinate system (ρ, , z), the radially polarized wave has an electric field in the radial direction with respect to the propagation axis while the magnetic field is aligned in the azimuthal orientation. There are no longitudinal components of the electric and magnetic fields. The incident electromagnetic field can be written as:
E i = a e x p i k z cos e x + sin e y , H i = a e x p i k z sin e x cos e y .
First, we discuss the surface integral in Equation (9). Using Equations (8) and (14), we can write quantities in Equation (9) as:
n × H = 2 a n z cos s e x + sin s e y n x n z cos s + n y n z sin s e z e x p i k z s , n × E = 0 , n · E = 2 a n x cos s + n y sin s e x p ( i k z s ) .
Substituting Equations (2), (10), (12), (13) and (15) into (9), we obtain the surface integral part in Equation (9):
E S P = E S ρ , z = a e x p 2 i k f 2 π i k δ π 0 2 π d θ s d s e x p i k u ρ s u ρ s 2 sin θ s × cos s cot θ s 2 1 1 i k u x u e ρ + cot θ s 2 1 1 1 i k u z u e z ,
where x = ρ s sin θ s cos s ρ , z = ρ s cos θ s z and u = x 2 + z 2 1 / 2 . P was supposed to lie on the x axis (see Figure 3). Due to the axial symmetry, E S P has no azimuthal component.
Now, we discuss the contour integral in Equation (9). According to Equation (4) and using the geometry shown in Figure 1 and Figure 2, the small vector line element on the circumference of the paraboloid is given by:
d s = 2 f cot δ 2 sin s e x + cos s e y d s .
From Equations (5), (8) and (14), the Cartesian coordinate components H x , H y ,   H z of the total magnetic field can be written as:
H x = 2 a sin s e x p i k z δ , H y = 2 a cos s e x p i k z δ , H z = 0 .
Hence:
H · d s = 4 a f cot δ 2 e x p i k z δ d s .
Substituting Equations (10), (12) and (19) into (9), we obtain the contour integral part in Equation (9):
E C P = E C ρ , z = a f e x p 2 i k f π cot δ 2 0 2 π e x p i k u ρ δ u × 1 1 i k u x u e ρ + z u e z d s ,
With x = ρ s ( δ ) sin δ cos s ρ , z = ρ s ( δ ) cos δ z and u = x 2 + z 2 1 / 2 . Due to the axial symmetry, E C P has no azimuthal component.
The total complex electric field can be given as:
E ρ , z = E S ρ , z + E C ρ , z .
We would like to call attention to the fact that in the following section, E ρ and E z represent the real parts of the complex field components.

4. Applications of the Theoretical Results for Different Focusing Conditions

To provide useful information for specialists who would like to apply the tightly focused field, we present and discuss some results of calculations using the previously introduced vector diffraction theory for various focusing geometries. Parabolic mirrors are standard devices in optical setups and are especially used at wavelength ranges (e.g., terahertz = THz) where conventional lenses and spherical mirrors fail. For the sake of generality, instead of specifying the focal length (f), the wavelength (λ) and r0 (that can be interpreted as the incident beam radius, see Figure 2) as absolute parameters, we introduce the λ/f and r0/f as relative parameters. This latter relates to the δ focusing angle (see Figure 2) as
δ = a r c c o s 1 4 r 0 f 2 1 1 4 r 0 f 2 + 1 ,
and
δ ¯ = π   -   δ
During the analysis, λ/f was kept below 0.1, since the λ/f > 0.1 case has less practical relevance. Supposing a typical value of f = 50 mm, the λ/f < 0.1 condition holds not only for the visible and (near-, mid-) infrared, but also for the THz frequency range (0.1–10 THz). Even in the case of λ/f = 0.1, the corresponding frequency is only 0.06 THz, so for λ/f > 0.1 the whole THz range is covered. The importance of the THz fields is outstanding due to their applicability for particle acceleration [17,18,19] because of their advantageous wavelength and because of the availability of pulses with extremely high pulse energies and electric field strengths owing to the tilted-pulse-front technique [20,21].
Before applying our above introduced theory based on the Stratton–Chu formulae, we must confirm its reliability. Therefore, as a validation, we compare its results with those of the commonly known scalar diffraction methods by treating separately two perpendicular (x and z) polarization components. The E ρ and E z field components determined with scalar diffraction methods (red, dashed curves in Figure 4) were deduced from the Rayleigh–Sommerfeld diffraction formula [22] adapted for the case shown in Figure 3 supposing a ‘thin’ focusing element. The amplitudes of the aperture functions belonging to the transversal and longitudinal field components just leaving the focusing element are
E F ρ = E 0 cos ε cos   F = E 0 f ρ F 2 + f 2 cos   F , E F z = E 0 sin ε = E 0 ρ F ρ F 2 + f 2
respectively, where E 0 is the electric field just before reaching the focusing element. The radially dependent part of the phase factor resulting from the focusing element is
e x p i φ F = e x p i k ρ F 2 2 f .
Hence
E ρ ρ 0 r 0 0 2 π E F ρ e x p i φ F exp i k u u 2 ρ F d ρ F d F , E z ρ 0 r 0 0 2 π E F z e x p i φ F exp i k u u 2 ρ F d ρ F d F ,
with
u = ρ F 2 + f 2 + ρ 2 2 ρ ρ F cos   F 1 2 .
We also determine the electric field distributions by using the formula of Endale et al. [13] for the case of a flat-top beam. This model is based on the diffraction theory of Richards–Wolf [9]. The curves computed this way are plotted by blue, dotted lines in Figure 4.
For parameters of λ/f = 10−3 and r0/f = 0.1 (with corresponding δ ¯ = 5.72 o a n g l e ), the contour term in our model is negligible beside the surface term. Furthermore, the parabolic mirror can be considered to be thin. In this paraxial regime, a particularly good agreement can be found between the results of our (Stratton–Chu) and the other two (Rayleigh–Sommerfeld, Richards–Wolf) models, as seen in Figure 4a,b.
For λ/f = 10−3 and for large numerical aperture (r0/f = 1.8 with corresponding δ ¯ = 83.9 o a n g l e ), the Rayleigh–Sommerfeld integral leads to a misleading result. Even the model of Endale et al. [13] (which is not restricted to small numerical apertures) works well only at the vicinity of the z axis, for low ρ/λ values.
These conclusions confirm the necessity of the use of our model in large numerical aperture regimes, when the spatial extension of the parabolic mirror in the z (longitudinal) direction becomes comparable with its transversal extension.
The contribution of the contour term is negligible in all cases examined in the following. However, for example in the case of r0/f = 0.2 and λ/f = 0.1, the surface and the contour terms have the same orders of magnitude, as illustrated in Figure 5. Here, the black dashed line belongs to the surface and the red dotted line to the contour term, while the blue solid line represents their sum, the total E z field component. Under the focusing conditions of Figure 5 for the case of the E ρ field component, the surface term dominates over the contour term.
In Figure 6, the normalized E ρ and E z distributions are plotted versus the radial coordinate normalized by the wavelength. This coordinate normalization makes it comfortable to plot the curves with different λ/f values on the same scale. Figure 6a,b belong to r0/f = 0.2 ( δ ¯ = 11.4°), Figure 6c,d to r0/f = 0.6 ( δ ¯ = 33.4°) and Figure 6d,e to r0/f = 1.8 ( δ ¯ = 83.9°). All curves of Figure 6 belong to the focal plane (z = 0). All curves are normalized to 1. However, in the graphs in brackets, one can find the field amplitude enhancement factors, h , relative to the amplitude of the incoming field, a , in order to provide information on how the magnitude of the field components scale with the r0/f and λ/f parameters, making it possible to compare the E ρ and E z amplitudes in case of a given focusing geometry.
All the curves of Figure 6 oscillate with decaying amplitude. At the focal point, E z has maxima, while E ρ is zero in each case, as expected. It was found that for the fixed r0/f value, the curves of different wavelengths could not be practically distinguished from each other below a threshold λ/f. These identical curves are indicated by the “≤” sign (see the labels). The threshold λ/f value increases with increasing r0/f.
From analyzing the enhancement factors (shown in the brackets), it is obvious that at constant r0/f the field strength increases with the decreasing wavelength. According to detailed calculations, below the threshold it scales with λ−1. If the wavelength is taken to be fixed, the field strength strongly increases with increasing r0/f, as expected. Comparing the amplitudes of curves with the same λ/f in Figure pairs Figure 6a,b, Figure 6c,d, Figure 6e,f, it is obvious that for the low value of r0/f (Figure 6a,b), the E ρ field component is larger than Ez, while for large r0/f (Figure 6e,f) the Ez is dominant, as expected intuitively. This analysis is especially interesting for those who plan particle acceleration applications. It is also seen in Figure 6 that the first zero crossing value decreases with increasing r0/f for both field components.
To more easily follow some characteristics, we plotted in the insets the absolute value of squared field components. The first maxima for the radial component and the first minima for both components shift to the left with increasing r0/f, although the degree of this shift decreases for larger r0/f.
Results like the ones shown in Figure 6 are needed when a waveguide-based electron accelerator is designed. If the aim is the efficient in-coupling into the waveguide by focusing, one should match the width (in ρ/λ) of the Ez distribution curve to the characteristic r1/λ value of the waveguide, where r1 is the core radius. For the optimized parameters (r1 = 380 µm core radius, d = 32 µm dielectric thickness, 0.6 THz frequency) of a dielectric coated metallic waveguide, r1/λ ≈ 0.7 [23]. Among the curves of Figure 6, the best agreement can be found in Figure 6d. This means that r0/f ≈ 0.6 can be an appropriate focusing geometry for in-coupling into the waveguide.
It is important to obtain information on the longitudinal distributions as well. As an example, the distribution of Ez along the optical axis was computed and plotted in Figure 7 for λ/f = 0.1 and r0/f = 0.6. It is clear from the graph that the peak of the curve is shifted from the focus towards the apex of the mirror, as already observed for the case of linear polarization [11]. The λ/f dependence of this shifting effect can be studied more clearly on the |Ez|2 curves (inset). Below λ/f = 10−2 this shift is negligible, while for λ/f = 10−1 the shift of the peak of the |Ez|2 curve is ~1.2 λ. The FWHM of the |Ez|2 curve is ~4.5 λ for λ/f = 10−1 and ~6 λ for λ/f ≤ 10−2. Examinations showed that for fixed λ/f, the shift is larger for lower r0/f values.

5. Conclusions

In this paper, based on the Stratton–Chu integral we derived general formulae to be used when a radially polarized electromagnetic wave is focused by a parabolic mirror with perfect reflectance. The theory was successfully validated within the limit of a small numerical aperture and negligible contour term when the Rayleigh–Sommerfeld diffraction formula served as a reliable alternative. We determined and analyzed the transversal and the longitudinal distribution of the radial and longitudinal electric field components for various focusing conditions. Our results can be interesting in the THz techniques, like THz imaging, linear and nonlinear THz spectroscopy, and especially for particle acceleration applications, where an intense longitudinal electric field is a requirement. Our work focused on monochromatic plane waves and simply connected paraboloid segments in on-axis geometry. But it can be further developed for the case of electromagnetic pulses, for the case of various input beams (e.g., vector Gaussian beam) and for various types of parabola segments (e.g., ring-like slice or off-axis geometry).

Author Contributions

Conceptualization, J.H. and L.P.; methodology, Z.T.G. and L.P.; software, Z.T.G. and L.P.; validation, Z.T.G. and L.P.; formal analysis, Z.T.G. and L.P.; investigation, Z.T.G. and L.P.; writing—original draft preparation, Z.T.G. and L.P.; writing—review and editing, Z.T.G., L.P. and J.H.; visualization, Z.T.G. and L.P.; supervision, L.P.; project administration, Z.T.G.; funding acquisition, J.H. and L.P. All authors have read and agreed to the published version of the manuscript.

Funding

The project has been supported by the Development and Innovation Fund of Hungary, financed under the TKP2021-EGA-17 funding scheme, and by the National Research, Development and Innovation Office (2018-1.2.1-NKP-2018-00010). This project has received funding from the European Union’s Horizon Europe research and innovation program under grant agreement No. 01046504.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the OPAL’ 2023 conference organizer and reviewer for their helpful remarks and recommendations.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bokor, N.; Davidson, N. 4π focusing with single paraboloid mirror. Opt. Commun. 2008, 28, 5499–5503. [Google Scholar] [CrossRef]
  2. April, A.; Piché, M. 4π Focusing of TM 01 beams under nonparaxial conditions. Opt. Express 2010, 18, 22128–22140. [Google Scholar] [CrossRef] [PubMed]
  3. Zeng, X.; Chen, X. Characterization of tightly focused vector fields formed by off-axis parabolic mirror. Opt. Express 2019, 27, 1179–1198. [Google Scholar] [CrossRef]
  4. Popov, K.; Bychenkov, V.Y.; Rozmus, W.; Sydora, R. Electron vacuum acceleration by a tightly focused laser pulse. Phys. Plasmas 2008, 15, 013108. [Google Scholar] [CrossRef]
  5. Jeong, T.M.; Weber, S.; Le Garrec, B.; Margarone, D.; Mocek, T.; Korn, G. Spatio-temporal modification of femtosecond focal spot under tight focusing condition. Opt. Express 2015, 23, 11641–11656. [Google Scholar] [CrossRef] [PubMed]
  6. Dorn, R.; Quabis, S.; Leuchs, G. Sharper focus for a radially polarized light beam. Phys. Rev. Lett. 2003, 91, 233901. [Google Scholar] [CrossRef]
  7. Marceau, V.; Varin, C.; Piché, M. Validity of the paraxial approximation for electron acceleration with radially polarized laser beams. Opt. Lett. 2013, 38, 821–823. [Google Scholar] [CrossRef] [Green Version]
  8. Ignatovsky, V.S. Diffraction by a parabolic mirror having arbitrary opening. Trans. Opt. Inst. Petrograd 1920, 1, 1–30. [Google Scholar]
  9. Richards, B.; Wolf, E. Electromagnetic diffraction in optical systems. II. Structure of the image field in an aplanatic system. Proc. R. Soc. Lond. Ser. A Math. Phys. Sci. 1959, 253, 358–379. [Google Scholar] [CrossRef]
  10. Stratton, J.A.; Chu, L. Diffraction theory of electromagnetic waves. Phys. Rev. 1939, 56, 99–107. [Google Scholar] [CrossRef]
  11. Varga, P.; Török, P. Focusing of electromagnetic waves by paraboloid mirrors. I. Theory. J. Opt. Soc. Am. A 2000, 17, 2081–2089. [Google Scholar] [CrossRef]
  12. Zeng, X. Configuration optimization of off-axis parabolic mirror for enhancing the focusability of a laser beam. Chin. Opt. Lett. 2021, 19, 032601. Available online: https://opg.optica.org/col/abstract.cfm?URI=col-19-3-032601 (accessed on 14 June 2023). [CrossRef]
  13. Endale, G.; Mohan, D.; Yadav, S. Focusing of Radially polarized Light Using Gaussian laser beam Near its focal length. IOP Conf. Ser. Mater. Sci. Eng. 2022, 1225, 012027. [Google Scholar] [CrossRef]
  14. Phan, T.; Sell, D.; Wang, E.W.; Doshay, S.; Edee, K.; Yang, J.; Fan, J.A. High-efficiency, large-area, topology-optimized metasurfaces. Light Sci. Appl. 2019, 8, 48. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Grand, J.; Le Ru, E.C. Practical implementation of accurate finite-element calculations for electromagnetic scattering by nanoparticles. Plasmonics 2020, 15, 109–121. [Google Scholar] [CrossRef]
  16. Guerboukha, H.; Shrestha, R.; Neronha, J.; Ryan, O.; Hornbuckle, M.; Fang, Z.; Mittlemana, D.M. Efficient leaky-wave antennas at terahertz frequencies generating highly directional beams. Appl. Phys. Lett. 2020, 117, 261103. [Google Scholar] [CrossRef]
  17. Zhang, D.; Fallahi, A.; Hemmer, M.; Wu, X.; Fakhari, M.; Hua, Y.; Cankaya, H.; Calendron, A.L.; Zapata, L.E.; Matlis, N.M.; et al. Segmented terahertz electron accelerator and manipulator (STEAM). Nat. Photonics 2018, 12, 336–342. [Google Scholar] [CrossRef]
  18. Nanni, E.A.; Huang, W.R.; Hong, K.H.; Ravi, K.; Fallahi, A.; Moriena, G.; Dwayne Miller, R.J.; Kärtner, F.X. Terahertz-driven linear electron acceleration. Nat. Commun. 2015, 6, 8486. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Hibberd, M.T.; Healy, A.L.; Lake, D.S.; Georgiadis, V.; Smith, E.J.; Finlay, O.J.; Pacey, T.H.; Jones, J.K.; Saveliev, Y.; Walsh, D.A. Acceleration of relativistic beams using laser-generated terahertz pulses. Nat. Photonics 2020, 14, 755–759. [Google Scholar] [CrossRef]
  20. Zhang, B.; Ma, Z.; Ma, J.; Wu, X.; Ouyang, C.; Kong, D.; Hong, T.; Wang, X.; Yang, P.; Chen, L. 1.4-mJ high energy terahertz radiation from lithium niobates. Laser Photonics Rev. 2021, 15, 2000295. [Google Scholar] [CrossRef]
  21. Fülöp, J.A.; Ollmann, Z.; Lombosi, C.; Skrobol, C.; Klingebiel, S.; Pálfalvi, L.; Krausz, F.; Karsch, S.; Hebling, J. Efficient generation of THz pulses with 0.4 mJ energy. Opt. Express 2014, 22, 20155–20163. [Google Scholar] [CrossRef] [PubMed]
  22. Goodman, J.W. Introduction to Fourier Optics, 3rd ed.; Roberts and Company Publisher: Englewood, CO, USA, 2005; pp. 46–49. [Google Scholar]
  23. Wong, L.J.; Fallahi, A.; Kärtner, F.X. Compact electron acceleration and bunch compression in THz waveguides. Opt. Express 2013, 21, 9792–9806. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Schematic diagram of the parabolic mirror with notations.
Figure 1. Schematic diagram of the parabolic mirror with notations.
Photonics 10 00848 g001
Figure 2. Schematic diagram with notations for the study of focusing radially polarized waves by a segment of a parabolic mirror.
Figure 2. Schematic diagram with notations for the study of focusing radially polarized waves by a segment of a parabolic mirror.
Photonics 10 00848 g002
Figure 3. Schematic diagram with notations for the study of focusing radially polarized wave with a focusing element (the blue circle in the x-y plane) having small numerical aperture. The observation point P is located at a distance ρ from the z axis. The dashed line is a ‘ray’ passing through the focus F.
Figure 3. Schematic diagram with notations for the study of focusing radially polarized wave with a focusing element (the blue circle in the x-y plane) having small numerical aperture. The observation point P is located at a distance ρ from the z axis. The dashed line is a ‘ray’ passing through the focus F.
Photonics 10 00848 g003
Figure 4. The transversal distribution of the radial (a,c) and longitudinal (b,d) component of the electric field at z = 0 for λ/f = 10−3 and r0/f = 0.1 (a,b), r0/f = 1.8 (c,d) parameter values. The black solid curves were computed using the Stratton–Chu integral for parabolic mirror, the red dashed ones using the Rayleigh–Sommerfeld integral and the blue dotted ones using theory of Endale et al. [13] (Richards–Wolf integral).
Figure 4. The transversal distribution of the radial (a,c) and longitudinal (b,d) component of the electric field at z = 0 for λ/f = 10−3 and r0/f = 0.1 (a,b), r0/f = 1.8 (c,d) parameter values. The black solid curves were computed using the Stratton–Chu integral for parabolic mirror, the red dashed ones using the Rayleigh–Sommerfeld integral and the blue dotted ones using theory of Endale et al. [13] (Richards–Wolf integral).
Photonics 10 00848 g004
Figure 5. The transversal distribution of the longitudinal component of the electric field at z = 0 for λ/f = 0.1 and r0/f = 0.2 parameter values. The total field (blue solid line) consists of a surface (black dashed line) and a contour (red dotted line) term.
Figure 5. The transversal distribution of the longitudinal component of the electric field at z = 0 for λ/f = 0.1 and r0/f = 0.2 parameter values. The total field (blue solid line) consists of a surface (black dashed line) and a contour (red dotted line) term.
Photonics 10 00848 g005
Figure 6. The transversal distribution of the radial (a,c,e) and longitudinal (b,d,f) components of the electric field at z = 0 for various λ/f and r0/f parameter values. The insets show the corresponding |E|2 distribution. Please note that the vertical scale range for (a) differs from the ones of (bf).
Figure 6. The transversal distribution of the radial (a,c,e) and longitudinal (b,d,f) components of the electric field at z = 0 for various λ/f and r0/f parameter values. The insets show the corresponding |E|2 distribution. Please note that the vertical scale range for (a) differs from the ones of (bf).
Photonics 10 00848 g006
Figure 7. The longitudinal distribution of the longitudinal component of the electric field at ρ = 0 for r0/f = 0.6 and λ/f = 0.1 parameter values. The inset shows the corresponding |Ez|2 distribution for λ/f ≤ 10−3, λ/f = 10−2 and 10−1 parameter values.
Figure 7. The longitudinal distribution of the longitudinal component of the electric field at ρ = 0 for r0/f = 0.6 and λ/f = 0.1 parameter values. The inset shows the corresponding |Ez|2 distribution for λ/f ≤ 10−3, λ/f = 10−2 and 10−1 parameter values.
Photonics 10 00848 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Godana, Z.T.; Hebling, J.; Pálfalvi, L. Focusing of Radially Polarized Electromagnetic Waves by a Parabolic Mirror. Photonics 2023, 10, 848. https://doi.org/10.3390/photonics10070848

AMA Style

Godana ZT, Hebling J, Pálfalvi L. Focusing of Radially Polarized Electromagnetic Waves by a Parabolic Mirror. Photonics. 2023; 10(7):848. https://doi.org/10.3390/photonics10070848

Chicago/Turabian Style

Godana, Zerihun Tadele, János Hebling, and László Pálfalvi. 2023. "Focusing of Radially Polarized Electromagnetic Waves by a Parabolic Mirror" Photonics 10, no. 7: 848. https://doi.org/10.3390/photonics10070848

APA Style

Godana, Z. T., Hebling, J., & Pálfalvi, L. (2023). Focusing of Radially Polarized Electromagnetic Waves by a Parabolic Mirror. Photonics, 10(7), 848. https://doi.org/10.3390/photonics10070848

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop