FOCUSING OF LIGHT PULSES BY DIFFERENT TYPES OF WAVY PARABOLIC MIRRORS

: The propagation of an ultra-short light pulse is studied in the framework of scalar diffraction theory. Light pulses are focused by different types of wavy parabolic surfaces. The temporal-spatial behavior of the two-dimensional wave field is computed in the vicinity of the focal plane. It is shown that the slightly perturbation from the perfect parabolic shape leads a space-time dispersion of the pulse in the neighborhood of the focus.


Introduction
Since the advent of the laser over sixty years ago, there has been intense interest in the quest to generate ultra-short laser pulses in the picosecond and femtosecond time scale. Nowadays great efforts have been devoted to generate and study pulses in the attosecond regime, based on highly nonlinear frequency conversion of femtosecond sources into the soft X-ray spectral region. As the period of one optical cycle in the near-infrared is just 2.7 femtoseconds, lasers with pulse widths in this range up to about 100 femtoseconds are generally termed, 'ultra-fast' lasers. As a pulse shorter than the period of the carrier frequency is basically a DC signal, optical pulses are now about as they are going to get. On the one hand a shorter optical pulse affords a higher time resolution in pump-probe measurements, allowing a wealth of physical processes to be studied as ultra-fast processes in chemical and biological materials, optical 50 G. NYITRAY, P. ELTER Pollack Periodica 15, 2020, 3 communications, medical imaging, ultra-fast time-resolved spectroscopy, and radio frequency signal processing. On the other hand a short optical pulse can produce extremely high peak powers with only moderate energies. A typical amplified Ti:Al 2 O 2 ultra-fast laser system can produce 1 mJ pulses of 100 fs duration, which corresponds to roughly 10 GW of power. When these pulses are focused, it is simple to obtain intensities greater than 10 15 W/cm 2 , where the electric field of the optical pulses is comparable to the binding field of electrons to ions. Many new nonlinear optical phenomena manifest themselves in this regime. The study of most nonlinear optical processes becomes routine, and even more exotic effects become accessible: highharmonic generation, attosecond X-ray generation and relativistic nonlinear optics.
Over the past two decades, the propagation characteristics of ultra-short, and hence ultra-wide bandwidth pulses, including spectral changes, the distortion of the pulse time, shape etc. have been investigated. The temporal and spatial behavior of a focused short pulse was studied in [1]- [7]. The spatial profile of a femtosecond pulse front can suffer a considerable distortion upon propagation through a lens. The delay between the pulse and phase fronts is proportional to the difference between the group and phase velocity of the light in the material of the lens. It is also proportional to the path length of the ray in the material of the lens, therefore the delay is different for various regions of the lens cross section and it can be described by a parabolic function of the distance from the optical axis. According to the numerical calculations based on geometric optics the pulse front even may form a loop, which means that the marginal rays have already passed the focus, while axially propagating part has not reached it yet [2]. Calculations based on scalar diffraction theory supported the prediction of geometric optics and show that the intensity distribution in the focal plane is different from the Airy pattern. The diffraction pattern consists of rings, with continuously increasing radius. Perhaps, the most unexpected result of that investigation is that a spike appears on the optical axis in front of the horseshoe-shaped pulse front. The properties of the boundary wave pulse are studied in [2]. A direct explanation of the formation of the boundary wave pulse is given in [8].
The current work aims to visualize the light field in the focal plane given by the imperfect (wavy) parabolic mirror. The transverse intensity distribution of the focused pulse is controlled by the undulation of the mirror surface.

Focusing of continuous waves
The paraboloid mirror is an ideal focusing device. An optical instrument that collects the parallel incident light into a single focal point, within the accuracy of the Gaussian optics, is called an ideal focusing device. Paraboloid mirrors have neither chromatic nor spherical aberration, which is inevitably property of a common lens. The nearly full paraboloid capable of focusing within 4 solid angle as compared with lens that can focus light only within 2 solid angle. A rigorous diffraction theory of focused light from parabolic mirrors has been developed over the course of nearly a century, beginning with a treatment in 1920 by Ignatovsky [9]. Another theoretical approach was provided by Richards and Wolf [10], [11], where strongly focused beams were accurately described.

LIGHT PULSES BY DIFFERENT TYPES OF WAVY PARABOLIC MIRRORS 51
Pollack Periodica 15, 2020, 3 For further discussion about focusing and diffraction the reader may refer to books e.g. [12], [13]. The paraboloid and spherical mirrors as dispersion-less optical instruments play important role in the field of generation short pulses of a few optical cycles. As it is shown in Fig. 1, the mirror is assumed to lie in the ( , ) plane and the observation plane ( ′, ′) is shifted by from it.
where 1, 2 / is the wavenumber, stands for wavelength and . The diameter of the mirror is 2 . The value of parameter is characterizes the paraboloid mirror shape.
The perfect paraboloid mirror converts the monochromatic plane waves into spherical-like waves having radius of curvature of converging to the focal point. If the mirror is unperturbed the Fourier transformation of the pupil function can be observed in the focal plane (see Fig. 1). The surface imperfection (waviness) of the mirror modifies the energy distribution of the focal plane as it can be seen in Fig. 1. The input is assumed to be uniformly illuminated by a normally incident, monochromatic plane wave of amplitude , in which case the disturbance incident on the mirror is where , is the pupil function yielding an effective illuminated region in the input plane. Usually the pupil function has a circular symmetry hence , circ / , where is a radial coordinate in the plane of the aperture. By definition if the G. NYITRAY, P. ELTER

Pollack Periodica 15, 2020, 3
circle function takes unit value. In case > the value of the circle function drops to zero and if its value is 1/2. The pupil function expresses the finite extent of the mirror in the transverse direction (limits the region of integration). One of the most remarkable and useful properties of a paraboloid is its inherent ability to perform two dimensional Fourier transforms in the focal plane. Let the radius of the mirror be . The Fourier-Bessel transform of a cylindrical pupil function is where , is a radial coordinate in the plane of the observation and -. is the first order Bessel function of the first kind. The amplitude distribution in the focal plane is seen to be The | / | intensity distribution in the focal plane is referred to as the Airy pattern. The | / | field distribution is circularly symmetric and consists of a central lobe and a series of concentric rings of diminishing amplitude. In case when the pupil function is one dimensional it can be expressed by a rectangle function 0#A /2 . The field distribution in the focal plane takes the following form: where 2V 2 is the ideal parabolic phase distribution, and WX> characterizes the imperfection (waviness) of the mirror. The quantity ℎ stands for the spatial modulation depth. It is assumed that the maximum value of the path deviation from the ideal mirror is small (10%). The function WX> can be called phase perturbation function, which can be periodic, quasi-periodic and aperiodic either.
Diffraction is a cornerstone of optical physics and has implications for the design of all optical system. According to the Huygens-Fresnel principle each point on the wavefront of a disturbance can be considered to be a new source of a secondary spherical disturbance 0 23> / , where is the distance between points Y and Y / . For a monochromatic field . Y, A . Y, Z 0 [2\] of frequency Z at a point Y, the diffraction field at a point Y / is given by the Fresnel diffraction integral as follows: Pollack Periodica 15, 2020, 3 where ab is an infinitesimal patch of area, assuming the incident field plane and the observation plane to be , and / , / coordinate planes, respectively, . Y, Z is the complex amplitude of the incident field, and the distance is given exactly by =cd / / . d is the distance of the observation plane from the aperture. To reduce the Huygens-Fresnel principle to a more simple and usable expression, one can introduce the Fresnel approximation. The approximation is based on the binominal expansion of the square root in the expression of . Retaining only the first two terms of the expansion takes the following form: This approximation means that the spherical wavelets of the Huygens-Fresnel principle have been replaced by wavelets with parabolic wavefronts. It is worth noting that in the Fresnel approximation Eq. (7) is seen to be a convolution < . ⨂i, where a ⨂ symbol between any two functions indicates that those functions are to be convolved and i 0 23> / is referred as the point-spread function (impulse-response) [12]. A wide-frequency-bandwidth . Y, A field can be represented in the form of a Fourier integral: where . Y, Z is the complex amplitude of monochromatic component with frequency Z and is given by the Fourier inversion formula for the monochromatic field component . Y, Z , and from the Fresnel diffraction integral the diffraction field < Y / , Z can be derived. The spectral intensity of the diffraction field is l < Y′, Z m < Y′, Z m . Let us consider the following incident Gaussian-shaped pulse with amplitude n : where Z n is the pulse central frequency and s is its time duration. The spectral distribution corresponding to the incident Gaussian-shaped pulse is G. NYITRAY, P. ELTER

Pollack Periodica 15, 2020, 3
The paraboloid mirror converts the monochromatic plane waves into spherical-like waves having radius of curvature of converging to the focal point, where is the wavenumber 2 / . For the sake of simplicity only two dimensional mirrors are studied (y dimension vanishes) hence the diffraction integral takes a simpler form The total phase delay takes the following form: 2 . The amplitude of the field of the pulse around the focus of the mirror can be calculated as the superposition of the fields of individual spectral components as where Z n is the central frequency, s stands for the pulse duration, A / expresses the time retardation A / A /# N/# and is the amplitude distribution of the incident wave and N 2V / . Taking into account wavy mirrors WX> in Eq. (6), it takes the following form: where ƒ n is the spatial frequency and the Taylor series of cosine function is needed for numerical calculation. Aperiodic phase distribution can be achieved by changing the local spatial frequency of the cosine function along the surface of the mirror, Fig. 2 where the parameter OE is responsible for the rate of instantaneous frequency changing. It should be noted that in terms of computer implementation [14]- [16] the usage of the cubic spline interpolation of the chirped cosine function is appropriate for numerical calculation.
The perturbation can be emphasized in the center (or in the edges) of the mirror by applying the following functions: . ℎ• .
WX> and ℎ• WX> , where    3e shows that an unperturbed mirror (Fig. 3a) creates the Fourier transform (sinc function) of the rectangular pupil function as expected. Numerical calculations show that the slightly periodic deviation from the perfect mirror shape leads a spatial dispersion of the pulse in the neighborhood of the focus as it can be seen in Fig. 3f -Fig. 3h. The phase disturbance can be characterized by a harmonic function. Since the value of modulation depth is ℎ 0.1 the maximum deviation from the perfect parabolic phase delay is 10%. It turns out that the effects of the mirror waviness are similar to the effects of sinusoidal phase grating. A thin sinusoidal phase grating can be characterized by the amplitude transmittance function as follows: where for simplicity it is assumed that the grating structure is bounded by an aperture of width 2a. The parameter • represents the peak to peak change of phase transmittance across the screen, and ƒ ¬W is the spatial frequency of the grating. If the grating is illuminated by a unit-amplitude, normally incident plane wave, then the field distribution immediately behind the screen is given precisely by Eq. (19). As it can be seen in Fig. 3f the introduction of the harmonic phase perturbation has thus deflected energy out of the zero order into multitude of higher orders. Note that the strengths of the various orders are symmetric about the zero order. Introduction of the harmonic phase disturbance leads the redistribution of energy between the side lobes. Due to the field modulation the energy is scattered to a larger volume surrounding the focal point compared to the unperturbed case. As mentioned earlier there are many other possibilities to introduce phase-disturbance deviating from the perfect mirror shape, for example, the envelope of the harmonic phase perturbation function can be modulated by a Gaussian type bell-curve (see Fig. 3c). Aperiodic waviness can be realized by introducing spatially chirped cosine function instead of periodic one as it can be seen in Fig. 4a. In order to emphasize the perturbation in the center of the mirror let the envelope of the chirp function is a Gaussian type bell-curve (see Fig. 4b). Due to the aperiodic phase perturbation with the effect of amplitude diminishing towards the edges the energy may leave the central region (zero order) and forms a double peak in the vicinity of the optical axis as it depicted by Fig. 4e and Fig. 5f. Creating dark central region on the optical axis can be useful because a beam of this kind has applications in optical tweezers for trapping particles with refractive index lower than the surrounding medium, where the gradient force will be directed toward the region of lowest intensity. Fig. 5 shows the complete temporalspatial distributions of all pulses, which were shown only in the spatial regime in Fig. 3 and in Fig. 4.  If the amplitude of the chirp function increases towards the edges of the mirror two almost merged peak can be observed in the central region (see Fig. 4f and Fig. 5g). During the simulation significant temporal broadening was not observed, however in case of extremely short pulses of a few optical cycles incident on a perfectly conducting spherical mirror under the focusing process significant temporal spreading was predicted, which avoids a fully compensated pulse [17]. In the current study the number of complete optical cycles is exactly six, which cannot be considered as a few optical cycles. The results of the current study may contribute to the field of generation and propagation of attosecond pulses [18]- [21].

Conclusion
The propagation of a focused ultra-short light pulse was studied in the framework of scalar diffraction theory. The light pulse was focused by an imperfect (wavy) parabolic mirror. Due to the harmonic and aperiodic phase perturbation the temporal and spatial behavior of the wave field is distorted in the vicinity of the focal plane. It turned out that the effects of the mirror waviness are similar to the effects of sinusoidal phase grating. The phase perturbation leads the redistribution of energy between the side lobes. The harmonic phase perturbation has thus deflected energy out of the zero order into multitude of higher orders. It was shown that in case of aperiodic phase perturbation a double peak is formed close to the optical axis, creating dark central region on the optical axis can be useful because a pulse of this kind has applications in optical tweezers for trapping particles. Other possible areas of application are the surface processing and micro lithography.