Maxim A. Efremov
German Aerospace Center (DLR), Institute of Quantum Technologies, 89081 Ulm, Germany
Institut für Quantenphysik and Center for Integrated Quantum Science and Technology (), Universität Ulm, 89081 Ulm, Germany
Felix Hufnagel
Nexus for Quantum Technologies, University of Ottawa, K1N 5N6, Ottawa, ON, Canada
Hugo Larocque
Nexus for Quantum Technologies, University of Ottawa, K1N 5N6, Ottawa, ON, Canada
Research Laboratory of Electronics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
Wolfgang P. Schleich
Institut für Quantenphysik and Center for Integrated Quantum Science and Technology (), Universität Ulm, 89081 Ulm, Germany
Hagler Institute for Advanced Study at Texas AM University, Texas AM AgriLife Research, Institute for Quantum Science and Engineering (IQSE), and Department of Physics and Astronomy, Texas AM University, College Station, Texas 77843-4242, USA
Ebrahim Karimi
Nexus for Quantum Technologies, University of Ottawa, K1N 5N6, Ottawa, ON, Canada
(June 19, 2024)
Abstract
Following the familiar analogy between the optical paraxial wave equation and the Schrödinger equation, we derive the optimal, real-valued wave function for focusing in one and two space dimensions without the use of any phase component. We compare and contrast the focusing parameters of the optimal waves with those of other diffractive focusing approaches, such as Fresnel zones. Moreover, we experimentally demonstrate these focusing properties on optical beams using both reflective and transmissive liquid crystal devices. Our results provide an alternative direction for focusing waves where phase elements are challenging to implement, such as for X-rays, THz radiation, and electron beams.
Introduction
Fresnel zone plates [1] are optical elements that focus an incident beam due to binary variations in its amplitude and phase. They offer precise control over diffractive propagation and enable efficient beam focusing in systems, where traditional lensing elements are not immediately available. In this article, we address the fundamental question whether any other approach to wave shaping can surpass the limit set by a Fresnel zone plate. In particular, we show for the case of matter waves that the answer to this question is a clear ”Yes!”, by deriving analytical solutions of the corresponding variational problem. Moreover, we demonstrate the supremacy of our approach compared to Fresnel zone plates by an experiment with light.
A scalar wave, such as a matter wave or an unpolarized electromagnetic field, comprises two components: an amplitude and a phase. The common way to focus an electromagnetic wave is to modulate its phase using a lens by applying a parabolic phase variation in space. However, there are waves for which a phase-modulating lens does not exist due to technological limitations in implementing phase-altering components in such systems. For instance, implementing such components for X-rays and matter waves often requires subnanometer manufacturing.
More effective approaches to focus waves can be achieved via amplitude modulation in space. For example, blocking part of the wave by a circular aperture or annular rings known as Fresnel zones will focus it to the Arago-Poisson spot [1]. In these examples, the incoming waves are spatially selected without being modified by the materials. These diffractive focusing techniques are crucially determined by a non-Gaussian initial wave function, as well as by the underlying dimensionality of the problem [2, 3, 4], and have been employed for surface gravity water waves and plasmonic waves [5, 6].
We emphasize that while Fresnel zones provide one approach to focusing the waves by amplitude modulation, one may question whether other approaches, e.g. nonbinary amplitude modulations, provide even better focusing. In the present article, we obtain the optimal initial wave function for focusing a free particle, i.e. matter waves, in one and two dimensions, and compare and contrast the focusing parameters of the optimal two-dimensional wave function to those of the Fresnel zone approach. The analogy between the Schrödinger equation and the paraxial Helmholtz equation allows us to extend our results to electromagnetic waves. Finally, we experimentally verify the focusing properties of the two-dimensional pattern at optical wavelengths using a reflective spatial light modulator and a fabricated transmissive liquid crystal device.
Results Theory of optimal focusing:
Our goal is to determine the optimal initial real-valued, aperture-constrained, and normalized wave function in two spatial dimensions that maximizes the intensity of the field on the symmetry axis at a prescribed focusing time. Our choice of the number of dimensions results from the fact that in one dimension the focusing is weaker, as shown in the Methods section.
Hence, we assume that is radially symmetric, as it provides the best diffractive focusing [3], and write the solution as
(1)
of the time-dependent two-dimensional Schödinger equation of a free particle in terms of the corresponding Green function
(2)
with the Bessel function of the first kind [7]. Here and are the dimensionless radial coordinate and time, respectively, wherein and denote the mass of the particle and the radius of the circular aperture. In the case of the two-dimensional paraxial Helmholtz equation, is equivalent to the longitudinal distance from the screen, where denotes the wave number.
We consider only wave functions that are truncated by the aperture and vanish elsewhere, . As a result, for a prescribed focusing time , or focal distance , the intensity along the symmetry axis, , takes the form
(3)
where we have used that is real.
In order to solve the optimization problem, we first construct the Lagrange function
(4)
where the Lagrange multiplier takes into account the normalization condition for , and then perform the variation of with respect to , to arrive at the eigenvalue problem
(5)
for the optimal wave function corresponding to the eigenvalue .
Since Eq. (5) is a linear integral equation with a degenerate kernel, its solution can be found analytically, as shown in the Methods section. Indeed, for a fixed value of , we obtain the maximum eigenvalue
(6)
and the normalized optimal initial wave function
(7)
Here, and are the amplitude parameter and the normalization constant, respectively, with being the sign function.
Substituting given by Eq. (7) into the expression, Eq. (3), for the intensity at , we prove that the intensity, indeed, achieves its maximum value for any given focusing time , or the dimensionless distance from the screen (within the paraxial approximation). In particular, for
(8)
where the integer counts the number of Fresnel zones that fit in the circular aperture , Eq. (6) yields
(9)
Fresnel zones:
Next we compare the maximum focusing intensity, Eq. (9), of the optimal state with the Fresnel zones approach. For this purpose, we consider two different designs.
An amplitude Fresnel zone (AFZ) plate alters the amplitude, while the phase Fresnel zone (PFZ) plate modifies the phase of the incoming wave. In the AFZ, only odd () annular zones are transparent, whereas even () zones are opaque, that is absorbing the incoming waves, with being the innermost zone containing the origin. The AFZ plate is a nonunitary object, i.e. the input intensity is not conserved. In the PFZ, we keep even and odd zones fully transparent; however, the phase in the even zones is shifted by .
For the focusing time , we derive in the Methods section the maximal intensities
(10)
and
(11)
at the symmetry axis .
A comparison of Eqs. (9), (10) and (11) reveals that the optimal state defined by Eq. (7) gives rise to focusing improved by the factor compared to the best Fresnel method.
Experiments: Now we demonstrate experimentally optimal diffractive focusing for the two-dimensional case using optical light. For this purpose we have fabricated a transmissive, liquid crystal optical element, that is a Pancharatnam-Berry optical element (PBOE) [8], described in the Methods section, which can be operated at many different wavelengths. It generates the optimal state given by Eq. (7). The space-varying amplitude for the focusing time defined by Eq. (8) with is shown in Fig. 1(a) together with the expected and measured intensities, Figs. 1(b) and (c), respectively.
The complete experimental apparatus used to generate the optimal state is displayed in Fig. 1(d). The PBOE placed between a half-wave plate and a polarizer is illuminated by a 633 nm He-Ne laser with an expanded Gaussian profile. A 4-f lens system is used to image the device on a pixel CCD camera placed on a translation stage, which allows us to measure the intensity of the modulated beam along its propagation to the focus. We have obtained this intensity profile in 50 steps for 25.0 mm.
Whereas Fig. 1(e) shows the exact evolution of the beam originating from the optimal state given by Eq. (7), Figs. 1(f) and (g) display the experimentally measured and expected intensity along the propagation. As further elaborated in the Discussion section, we expect imperfections in our optical system to affect the propagation profile.
To compare the propagation of with the ones created by the Fresnel zone plates, Eqs. (M9) and (M10) in the Methods section, we replace our PBOE by a reflective spatial light modulator (SLM), as depicted in Fig. 1(h).
We used a Hamamatsu liquid crystal on silicon (LCOS) SLM with resolution and a pixel size of m. Moreover, we are able to encode both the intensity and the phase of the pattern on the incident beam using an amplitude masking technique [9]. A phase diffraction grating is added to the pattern on the SLM which produces the desired field in the first order of diffraction. We then select this first order with a 4f lens system and an iris, thereby allowing us to remove all other diffraction orders while imaging the SLM plane onto our moveable CCD camera. Although SLMs do not reach the spatial resolution of our PBOE, their programmability can more readily streamline experiments comparing various focusing approaches.
In Fig. 2, we display the maximal focusing intensity for nine different patterns corresponding to focusing times , as defined in Eq. (8), with . The results of these experiments involving the optimal wave function as well as both forms of the Fresnel zones are depicted in Fig. 2.
Discussion
The propagation of the beam, shown in Fig. 1(f), features an artificial peat at arising from the modulation by our PBOE. Indeed, an imperfect polarization alignment in our generation apparatus leads to contributions of the -term in Eq. (.4). The angular orientation of the liquid crystals, , ranges from 0 to such that the horizontally polarized component of the field oscillates from to according to , while the term oscillates from 0 to +1 and back to 0. As a consequence, the contribution from the -component, which does not have negative amplitudes, behaves like the Fresnel zones, giving rise to a focus at . In Fig. 1(g), we show the expected propagation including the contribution from the -term, which is in good agreement with our experimental results shown in Fig. 1(f).
The scalings of the peak intensities of the focusing methods, considered in our article, with are shown in Fig. 2 by solid lines together with the experimental results depicted by the differently coloured data points. The peak focal intensity increases with increasing , corresponding to a tighter focusing time , or equivalently, to a shorter focal length. Furthermore, the optimal wave function consistently outperforms both methods relying on Fresnel zone plates.
We conclude this discussion by emphasizing that with our PBOE we were able to achieve a better resolution in our pattern creation, allowing for a tighter focusing time with , than with the SLM. This advantage is primarily due to the fact that a diffraction grating is necessary when the SLM is used to form an arbitrary wave function, which thus limits the maximum spatial frequency of the phase oscillations corresponding to the desired pattern. In addition, the SLM has limited control over both phase and spatial modulation, as prescribed by its bit depth and pixel pitch, respectively. As is increased, the number of oscillations in from to increases and in particular, the outer rings of the pattern become ever thinner.
Summary.
We have derived the optimal real-valued matter wave for focusing in both one- and two- dimensions. The analogy between the Schrödinger equation and the paraxial wave equation allows us to transfer our treatment to light. In our optical experiment, we have realized the two-dimensional optimal wave function using liquid crystal devices, verifying the superior focusing properties of compared to diffractive focusing from Fresnel zone patterns.
The optimal diffractive patterns derived here may be of interest to many different communities where phase modulation, due to technological limitations, is not directly possible. We can also envision extending this technique to vector fields, such as spinors in both optical and matter waves, where combinations of amplitude masks and specially polarized vector modes bring highly structured variations in focused beams [10, 11, 12]. The application of to these tight focusing problems remains to be explored.
Methods
.1 Optimal state in two dimensions
To obtain the analytical solution of the integral equation (5), we cast it in the form
(M1)
where
(M2)
and
(M3)
are functions solely of .
Next, we insert given by Eq. (M1) into Eqs. (M2) and (M3), and obtain the system
(M4)
(M5)
of algebraic equations for and , which has non-trivial solutions, only when its determinant is zero, that is
(M6)
This elementary quadratic equation has the two solutions
(M7)
By inserting the maximal eigenvalue into Eq. (M4), we find the relation between and , and thus the normalized optimal initial wave function given by Eq. (7).
.2 Amplitude and phase Fresnel zones:
Maximal intensity
For a given value of the radii
(M8)
of the Fresnel zones, with , extend to the maximum number of zones fitting within the circular aperture [1]. Therefore, the initial wave functions and for the amplitude and phase Fresnel zone patterns read
(M9)
and
(M10)
with and . Here, denotes the Heaviside function and is a normalization constant depending on .
To derive an analytical formula for the maximal intensity, we choose the focusing times . In this case, Eq. (M8) reduces to , and the normalization condition
Next, we insert the initial profile given by Eq. (M9) into Eq. (3), and obtain the expression
that is
(M13)
where we have used the fact that .
As a result, Eq. (M13) combined with Eq. (M12) gives rise to the maximum intensity , Eq. (10), produced by the amplitude Fresnel zones. Analogously, we derive the corresponding maximum intensity , Eq. (11), for the Fresnel phase zones.
.3 Optimal state in one dimension
In this section we determine the optimal initial real-valued and normalized wave function that maximizes the intensity of the field at , at the focusing time .
In this case we use the one-dimensional Green function
(M14)
for the time-dependent one-dimensional Schrödinger equation of a free particle. Here, and are the dimensionless position and time, respectively, and and denote the mass of the particle and the slit width.
We again apply the method of the Lagrange multipliers and arrive at the eigenvalue problem
(M15)
for the optimal initial wave function with the eigenvalue , that determines the maximum intensity achieved at . Here, we have assumed that for .
Since we are interested in the maximum of the intensity, we solve the integral equation (M15) only for the largest eigenvalue . As a result, for a given , we find the optimal initial wave function
(M16)
for , with for , and the corresponding maximal eigenvalue
Our device consists of a patterned layer of birefringent nematic liquid crystals whose orientation locally determines that of the medium’s optical axis. This feature causes the element to have the action
on the horizontal and vertical polarization components of an optical beam. Here is the optical retardation of the liquid crystal molecules and is the device’s spatially dependent liquid crystal axis orientation expressed in terms of the transverse Cartesian coordinates and .
When the device is perfectly tuned, that is for , and followed by a horizontally oriented polarizer, it can effectively be used to mask the amplitude profile of incoming horizontally polarized light by a factor of . This procedure was employed to generate our real-valued optimal state by means of a device defined by an optical axis of , where is the maximum value of the optimal state.
References
[1]
Born, M. & Wolf, E.
Principles of Optics: Electromagnetic Theory of
Propagation, Interference and Diffraction of Light
(Elsevier, 2013).
[2]
Cirone, M. A., Rza¸żewski, K.,
Schleich, W. P., Straub, F. &
Wheeler, J. A.
Quantum anticentrifugal force.
Phys. Rev. A65,
022101 (2001).
[3]
Białynicki-Birula, I., Cirone, M. A.,
Dahl, J. P., Fedorov, M. &
Schleich, W. P.
In- and outbound spreading of a free-particle
-wave.
Phys. Rev. Lett.89, 060404
(2002).
[4]
Case, W. B., Sadurni, E. &
Schleich, W. P.
A diffractive mechanism of focusing.
Opt. Express20,
27253 (2012).
[5]
Weisman, D. et al.Diffractive focusing of waves in time and in space.
Phys. Rev. Lett.118, 154301
(2017).
[6]
Weisman, D. et al.Diffractive guiding of waves by a periodic array of
slits.
Phys. Rev. Lett.127, 014303
(2021).
[7]
Abramowitz, M. & Stegun, I. A.
Handbook of Mathematical Functions with
Formulas, Graphs, and Mathematical Tables, vol. 55
(US Government Printing Office, 1968).
[8]
Larocque, H. et al.Arbitrary optical wavefront shaping via spin-to-orbit
coupling.
Journal of Optics18, 124002
(2016).
[9]
Bolduc, E., Bent, N.,
Santamato, E., Karimi, E. &
Boyd, R. W.
Exact solution to simultaneous intensity and phase
encryption with a single phase-only hologram.
Opt. Lett.38,
3546–3549 (2013).
[10]
Dorn, R., Quabis, S. &
Leuchs, G.
Sharper focus for a radially polarized light beam.
Phys. Rev. Lett.91, 233901
(2003).
[11]
Wang, H., Shi, L.,
Lukyanchuk, B., Sheppard, C. &
Chong, C. T.
Creation of a needle of longitudinally polarized
light in vacuum using binary optics.
Nature Photonics2, 501–505
(2008).
[12]
Karimi, E., Piccirillo, B.,
Marrucci, L. & Santamato, E.
Improved focusing with hypergeometric-gaussian
type-ii optical modes.
Opt. Express16,
21069 (2008).
Data availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.
Code availability
The code used for the data analysis is available from the corresponding author upon reasonable request.
Ethics declarations The authors declare no competing interests.
Corresponding authors
Correspondence and requests for materials should be addressed to maxim.efremov@dlr.de.
Acknowledgement
Maxim A. Efremov and Felix Hufnagel contributed equally to this work. We thank P. Boegel for fruitful discussions. W.P.S. is grateful to the Hagler Institute for Advanced Study at Texas AM University for a Faculty Fellowship, and to Texas AM AgriLife Research for the support of this work. This project was conceived during a visit of E.K. to Ulm University made possible by IQST. F. H. and E.K. acknowledge the support of the Canada Research Chair (CRC) Program, NRC-uOttawa Joint Centre for Extreme Quantum Photonics (JCEP) and NSERC.