Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Optimal Diffractive Focusing of Quantum Waves

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 (𝐼𝑄𝑆𝑇𝐼𝑄𝑆𝑇\it IQSTitalic_IQST), 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 (𝐼𝑄𝑆𝑇𝐼𝑄𝑆𝑇\it IQSTitalic_IQST), Universität Ulm, 89081 Ulm, Germany Hagler Institute for Advanced Study at Texas A&\&&M University, Texas A&\&&M AgriLife Research, Institute for Quantum Science and Engineering (IQSE), and Department of Physics and Astronomy, Texas A&\&&M 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 ψ0subscript𝜓0\psi_{0}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT in two spatial dimensions that maximizes the intensity |ψ|2superscript𝜓2\left|\psi\right|^{2}| italic_ψ | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT 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 ψ0subscript𝜓0\psi_{0}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is radially symmetric, as it provides the best diffractive focusing [3], and write the solution as

ψ(ρ,τ)=2π0ρdρG(2)(ρ,τ|ρ,0)ψ0(ρ)𝜓𝜌𝜏2𝜋superscriptsubscript0superscript𝜌differential-dsuperscript𝜌superscriptG2𝜌conditional𝜏superscript𝜌0subscript𝜓0superscript𝜌\psi(\rho,\tau)=2\pi\int\limits_{0}^{\infty}\rho^{\prime}{\rm d}\rho^{\prime}% \,{\rm G}^{(2)}(\rho,\tau|\rho^{\prime},0)\,\psi_{0}(\rho^{\prime})italic_ψ ( italic_ρ , italic_τ ) = 2 italic_π ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_ρ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT roman_d italic_ρ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT roman_G start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT ( italic_ρ , italic_τ | italic_ρ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , 0 ) italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_ρ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) (1)

of the time-dependent two-dimensional Schödinger equation of a free particle in terms of the corresponding Green function

G(2)(ρ,τ|ρ,0)=12πiτexp(iρ2+ρ22τ)J0(ρρτ)superscriptG2𝜌conditional𝜏superscript𝜌012𝜋i𝜏isuperscript𝜌2superscript𝜌22𝜏subscript𝐽0𝜌superscript𝜌𝜏{\rm G}^{\rm(2)}(\rho,\tau|\rho^{\prime},0)=\frac{1}{2\pi{\rm i}\tau}\exp\left% ({\rm i}\frac{\rho^{2}+\rho^{\prime 2}}{2\tau}\right)J_{0}\left(\frac{\rho\rho% ^{\prime}}{\tau}\right)roman_G start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT ( italic_ρ , italic_τ | italic_ρ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , 0 ) = divide start_ARG 1 end_ARG start_ARG 2 italic_π roman_i italic_τ end_ARG roman_exp ( roman_i divide start_ARG italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_ρ start_POSTSUPERSCRIPT ′ 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_τ end_ARG ) italic_J start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( divide start_ARG italic_ρ italic_ρ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_τ end_ARG ) (2)

with the Bessel function J0subscript𝐽0J_{0}italic_J start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT of the first kind [7]. Here ρr/R𝜌𝑟𝑅\rho\equiv r/Ritalic_ρ ≡ italic_r / italic_R and τt/(MR2)𝜏Planck-constant-over-2-pi𝑡𝑀superscript𝑅2\tau\equiv\hbar t/(MR^{2})italic_τ ≡ roman_ℏ italic_t / ( italic_M italic_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) are the dimensionless radial coordinate and time, respectively, wherein M𝑀Mitalic_M and R𝑅Ritalic_R denote the mass of the particle and the radius of the circular aperture. In the case of the two-dimensional paraxial Helmholtz equation, τ𝜏\tauitalic_τ is equivalent to the longitudinal distance z(kR2)τ𝑧𝑘superscript𝑅2𝜏z\equiv(kR^{2})\tauitalic_z ≡ ( italic_k italic_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) italic_τ from the screen, where k𝑘kitalic_k denotes the wave number.

We consider only wave functions ψ0subscript𝜓0\psi_{0}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT that are truncated by the aperture ρ1𝜌1\rho\leq 1italic_ρ ≤ 1 and vanish elsewhere, ψ0(ρ1)=0subscript𝜓0𝜌10\psi_{0}(\rho\geq 1)=0italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_ρ ≥ 1 ) = 0. As a result, for a prescribed focusing time τfsubscript𝜏f\tau_{\rm f}italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT, or focal distance zfkR2τfsubscript𝑧f𝑘superscript𝑅2subscript𝜏fz_{\rm f}\equiv kR^{2}\tau_{\rm f}italic_z start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT ≡ italic_k italic_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT, the intensity I[ψ0]𝐼delimited-[]subscript𝜓0I[\psi_{0}]italic_I [ italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ] along the symmetry axis, ρ=0𝜌0\rho=0italic_ρ = 0, takes the form

I[ψ0]=1τf201udu01vdvcos(u2v22τf)ψ0(u)ψ0(v),𝐼delimited-[]subscript𝜓01superscriptsubscript𝜏f2superscriptsubscript01𝑢differential-d𝑢superscriptsubscript01𝑣differential-d𝑣superscript𝑢2superscript𝑣22subscript𝜏fsubscript𝜓0𝑢subscript𝜓0𝑣\displaystyle I[\psi_{0}]=\frac{1}{\tau_{\rm f}^{2}}\int\limits_{0}^{1}u{\rm d% }u\int\limits_{0}^{1}v{\rm d}v\cos\left(\frac{u^{2}-v^{2}}{2\tau_{\rm f}}% \right)\psi_{0}(u)\psi_{0}(v),\qquaditalic_I [ italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ] = divide start_ARG 1 end_ARG start_ARG italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT italic_u roman_d italic_u ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT italic_v roman_d italic_v roman_cos ( divide start_ARG italic_u start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_v start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT end_ARG ) italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_u ) italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_v ) , (3)

where we have used that ψ0subscript𝜓0\psi_{0}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is real.

In order to solve the optimization problem, we first construct the Lagrange function

[ψ0]I[ψ0]λ[2π01uduψ02(u)1],delimited-[]subscript𝜓0𝐼delimited-[]subscript𝜓0𝜆delimited-[]2𝜋superscriptsubscript01𝑢differential-d𝑢superscriptsubscript𝜓02𝑢1{\mathcal{L}}[\psi_{0}]\equiv I[\psi_{0}]-\lambda\left[2\pi\int\limits_{0}^{1}% u{\rm d}u\,\psi_{0}^{2}(u)-1\right],caligraphic_L [ italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ] ≡ italic_I [ italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ] - italic_λ [ 2 italic_π ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT italic_u roman_d italic_u italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_u ) - 1 ] , (4)

where the Lagrange multiplier λ𝜆\lambdaitalic_λ takes into account the normalization condition for ψ0subscript𝜓0\psi_{0}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, and then perform the variation of [ψ0]delimited-[]subscript𝜓0{\mathcal{L}}[\psi_{0}]caligraphic_L [ italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ] with respect to ψ0subscript𝜓0\psi_{0}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, to arrive at the eigenvalue problem

12πτf201vdvcos(u2v22τf)ψ0(v)=λψ0(u)12𝜋superscriptsubscript𝜏f2superscriptsubscript01𝑣differential-d𝑣superscript𝑢2superscript𝑣22subscript𝜏fsubscript𝜓0𝑣𝜆subscript𝜓0𝑢\frac{1}{2\pi\tau_{\rm f}^{2}}\int\limits_{0}^{1}v{\rm d}v\cos\left(\frac{u^{2% }-v^{2}}{2\tau_{\rm f}}\right)\psi_{0}(v)=\lambda\psi_{0}(u)divide start_ARG 1 end_ARG start_ARG 2 italic_π italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT italic_v roman_d italic_v roman_cos ( divide start_ARG italic_u start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_v start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT end_ARG ) italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_v ) = italic_λ italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_u ) (5)

for the optimal wave function ψ0subscript𝜓0\psi_{0}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT corresponding to the eigenvalue λ𝜆\lambdaitalic_λ.

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 τfsubscript𝜏f\tau_{\rm f}italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT, we obtain the maximum eigenvalue

λ+(τf)=18πτf2[1+2τf|sin(12τf)|]subscript𝜆subscript𝜏f18𝜋superscriptsubscript𝜏f2delimited-[]12subscript𝜏f12subscript𝜏f\lambda_{+}(\tau_{\rm f})=\frac{1}{8\pi\tau_{\rm f}^{2}}\left[1+2\tau_{\rm f}% \left|\sin{\left(\frac{1}{2\tau_{\rm f}}\right)}\right|\,\right]italic_λ start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ( italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT ) = divide start_ARG 1 end_ARG start_ARG 8 italic_π italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ 1 + 2 italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT | roman_sin ( divide start_ARG 1 end_ARG start_ARG 2 italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT end_ARG ) | ] (6)

and the normalized optimal initial wave function

ψ0(opt)(ρ)=N[1+acos(ρ22τf)+1asin(ρ22τf)].superscriptsubscript𝜓0opt𝜌𝑁delimited-[]1𝑎superscript𝜌22subscript𝜏f1𝑎superscript𝜌22subscript𝜏f\psi_{0}^{\rm(opt)}(\rho)=N\left[\sqrt{1+a}\cos\left(\frac{\rho^{2}}{2\tau_{% \rm f}}\right)+\sqrt{1-a}\sin\left(\frac{\rho^{2}}{2\tau_{\rm f}}\right)\right].italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT ( italic_ρ ) = italic_N [ square-root start_ARG 1 + italic_a end_ARG roman_cos ( divide start_ARG italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT end_ARG ) + square-root start_ARG 1 - italic_a end_ARG roman_sin ( divide start_ARG italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT end_ARG ) ] . (7)

Here, acos[1/(2τf)]sign{sin[1/(2τf)]}𝑎12subscript𝜏fsign12subscript𝜏fa\equiv\cos[1/(2\tau_{\rm f})]{\rm sign}\left\{\sin[1/(2\tau_{\rm f})]\right\}italic_a ≡ roman_cos [ 1 / ( 2 italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT ) ] roman_sign { roman_sin [ 1 / ( 2 italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT ) ] } and N1/8π2τf2λ+(τf)𝑁18superscript𝜋2superscriptsubscript𝜏f2subscript𝜆subscript𝜏fN\equiv 1/\sqrt{8\pi^{2}\tau_{\rm f}^{2}\lambda_{+}(\tau_{\rm f})}italic_N ≡ 1 / square-root start_ARG 8 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_λ start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ( italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT ) end_ARG are the amplitude parameter and the normalization constant, respectively, with sign(x)sign𝑥{\rm sign}(x)roman_sign ( italic_x ) being the sign function.

Substituting ψ0(opt)superscriptsubscript𝜓0opt\psi_{0}^{\rm(opt)}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT given by Eq. (7) into the expression, Eq. (3), for the intensity at ρ=0𝜌0\rho=0italic_ρ = 0, we prove that the intensity, indeed, achieves its maximum value Imax(opt)(τf)I[ψ0(opt)]=λ+(τf)superscriptsubscript𝐼maxoptsubscript𝜏f𝐼delimited-[]superscriptsubscript𝜓0optsubscript𝜆subscript𝜏fI_{\rm max}^{\rm(opt)}(\tau_{\rm f})\equiv I[\psi_{0}^{\rm(opt)}]=\lambda_{+}(% \tau_{\rm f})italic_I start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT ( italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT ) ≡ italic_I [ italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT ] = italic_λ start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ( italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT ) for any given focusing time τfsubscript𝜏f\tau_{\rm f}italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT, or the dimensionless distance zfsubscript𝑧fz_{\rm f}italic_z start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT from the screen (within the paraxial approximation). In particular, for

τn012πn0,subscript𝜏subscript𝑛012𝜋subscript𝑛0\tau_{n_{0}}\equiv\frac{1}{2\pi n_{0}},italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ≡ divide start_ARG 1 end_ARG start_ARG 2 italic_π italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG , (8)

where the integer n0subscript𝑛0n_{0}italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT counts the number of Fresnel zones that fit in the circular aperture 0ρ10𝜌10\leq\rho\leq 10 ≤ italic_ρ ≤ 1, Eq. (6) yields

Imax(opt)(τn0)=π2n02.superscriptsubscript𝐼maxoptsubscript𝜏subscript𝑛0𝜋2superscriptsubscript𝑛02I_{\rm max}^{\rm(opt)}(\tau_{n_{0}})=\frac{\pi}{2}n_{0}^{2}.italic_I start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT ( italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) = divide start_ARG italic_π end_ARG start_ARG 2 end_ARG italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT . (9)

Fresnel zones: Next we compare the maximum focusing intensity, Eq. (9), of the optimal state ψ0(opt)superscriptsubscript𝜓0opt\psi_{0}^{\rm(opt)}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT 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 (n=1,3,5,𝑛135n=1,3,5,\ldotsitalic_n = 1 , 3 , 5 , …) annular zones are transparent, whereas even (n=2,4,6,𝑛246n=2,4,6,\ldotsitalic_n = 2 , 4 , 6 , …) zones are opaque, that is absorbing the incoming waves, with n=1𝑛1n=1italic_n = 1 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 π𝜋\piitalic_π.

Refer to caption
Figure 1: Implementation of Optimal Diffractive Focusing. We show the theoretical amplitude (a) and intensity (b) of the n0=14subscript𝑛014n_{0}=14italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 14 optimal wave together with the measured optical transmission through the fabricated focusing element (c) for the focusing time τn0=1/(2πn0)subscript𝜏subscript𝑛012𝜋subscript𝑛0\tau_{n_{0}}=1/(2\pi n_{0})italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = 1 / ( 2 italic_π italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ). (d) Experimental apparatus used to generate the optimal state with a PBOE. A linearly polarized 633nm633nm633\,{\rm nm}633 roman_nm Gaussian beam passes a half-wave plate (λ/2𝜆2\lambda/2italic_λ / 2) which rotates it to the horizontal polarization. The beam is then expanded by a factor of five by two lenses (L) to obtain a relatively flat profile before it goes through the PBOE followed by a polarizer. The latter is then imaged by a 4f-system in order to examine its propagation dynamics by a CCD camera. Numerically simulated (e) and experimentally observed (f) intensity distributions for a cross section of the beam as it propagates from the plane of the device to the focus for the optimal state. Numerical simulation taking into account contributions (g) from both the horizontal (Cosine) and vertical (Sine) polarization components of the modulated beam. Experimental setup (h) for focusing with Fresnel zone patterns using the SLM. The 633 nm laser source is expanded to cover the SLM. The 4f-lens system then images the SLM onto the CCD camera with an iris (I) placed at the focus to select the first order of diffraction.

For the focusing time τn0subscript𝜏subscript𝑛0\tau_{n_{0}}italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT, we derive in the Methods section the maximal intensities

Imax(AFZ)(τn0)=2π{n0(n0+1),n0=1,3,5,n02,n0=2,4,6,superscriptsubscript𝐼maxAFZsubscript𝜏subscript𝑛02𝜋casessubscript𝑛0subscript𝑛01subscript𝑛0135superscriptsubscript𝑛02subscript𝑛0246I_{\rm max}^{\rm(AFZ)}(\tau_{n_{0}})=\frac{2}{\pi}\begin{cases}n_{0}(n_{0}+1),% &\quad n_{0}=1,3,5,\ldots\\ n_{0}^{2},&\quad n_{0}=2,4,6,\ldots\end{cases}italic_I start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_AFZ ) end_POSTSUPERSCRIPT ( italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) = divide start_ARG 2 end_ARG start_ARG italic_π end_ARG { start_ROW start_CELL italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + 1 ) , end_CELL start_CELL italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1 , 3 , 5 , … end_CELL end_ROW start_ROW start_CELL italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , end_CELL start_CELL italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 2 , 4 , 6 , … end_CELL end_ROW (10)

and

Imax(PFZ)(τn0)=4πn02superscriptsubscript𝐼maxPFZsubscript𝜏subscript𝑛04𝜋superscriptsubscript𝑛02I_{\rm max}^{\rm(PFZ)}(\tau_{n_{0}})=\frac{4}{\pi}n_{0}^{2}italic_I start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_PFZ ) end_POSTSUPERSCRIPT ( italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) = divide start_ARG 4 end_ARG start_ARG italic_π end_ARG italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT (11)

at the symmetry axis ρ=0𝜌0\rho=0italic_ρ = 0.

A comparison of Eqs. (9), (10) and (11) reveals that the optimal state ψ0(opt)superscriptsubscript𝜓0opt\psi_{0}^{\rm(opt)}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT defined by Eq. (7) gives rise to focusing improved by the factor π2/8superscript𝜋28\pi^{2}/8italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / 8 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 ψ0(opt)superscriptsubscript𝜓0opt\psi_{0}^{\rm(opt)}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT given by Eq. (7). The space-varying amplitude for the focusing time τn0subscript𝜏subscript𝑛0\tau_{n_{0}}italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT defined by Eq. (8) with n0=14subscript𝑛014n_{0}=14italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 14 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 1920×1080192010801920\times 10801920 × 1080 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 μm𝜇m\mu\text{m}italic_μ m steps for 25.0 mm.

Whereas Fig. 1(e) shows the exact evolution of the beam originating from the optimal state ψ0(opt)superscriptsubscript𝜓0opt\psi_{0}^{\rm(opt)}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT 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.

Refer to caption
Figure 2: Comparison of the three methods for two-dimensional diffractive focusing. For a fixed focusing time τn0=1/(2πn0)subscript𝜏subscript𝑛012𝜋subscript𝑛0\tau_{n_{0}}=1/(2\pi n_{0})italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = 1 / ( 2 italic_π italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ), with n0=2,3,,10subscript𝑛02310n_{0}=2,3,\ldots,10italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 2 , 3 , … , 10, we display the theoretical optimal maximum intensities, Imax(opt)(τn0)superscriptsubscript𝐼maxoptsubscript𝜏subscript𝑛0I_{\rm max}^{\rm(opt)}(\tau_{n_{0}})italic_I start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT ( italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) (blue line), given by Eq. (9), as well as Imax(AFZ)(τn0)superscriptsubscript𝐼maxAFZsubscript𝜏subscript𝑛0I_{\rm max}^{\rm(AFZ)}(\tau_{n_{0}})italic_I start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_AFZ ) end_POSTSUPERSCRIPT ( italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) (red and green), and Imax(PFZ)(τn0)superscriptsubscript𝐼maxPFZsubscript𝜏subscript𝑛0I_{\rm max}^{\rm(PFZ)}(\tau_{n_{0}})italic_I start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_PFZ ) end_POSTSUPERSCRIPT ( italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) (orange), associated with the Fresnel zone plates, Eqs. (10) and (11). The points represent the corresponding measurements.

To compare the propagation of ψ0(opt)superscriptsubscript𝜓0opt\psi_{0}^{\rm(opt)}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT 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 1272×1024127210241272\times 10241272 × 1024 resolution and a pixel size of 12.5μ12.5𝜇12.5\leavevmode\nobreak\ \mu12.5 italic_μ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 Imax(τn0)subscript𝐼maxsubscript𝜏subscript𝑛0I_{\text{max}}(\tau_{n_{0}})italic_I start_POSTSUBSCRIPT max end_POSTSUBSCRIPT ( italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) for nine different patterns corresponding to focusing times τn0subscript𝜏subscript𝑛0\tau_{n_{0}}italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT, as defined in Eq. (8), with n0=2,3,,10subscript𝑛02310n_{0}=2,3,\ldots,10italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 2 , 3 , … , 10. 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 τ=τf/2𝜏subscript𝜏f2\tau=\tau_{\rm f}/2italic_τ = italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT / 2 arising from the modulation by our PBOE. Indeed, an imperfect polarization alignment in our generation apparatus leads to contributions of the sin(2α)2𝛼\sin(2\alpha)roman_sin ( 2 italic_α )-term in Eq. (.4). The angular orientation of the liquid crystals, α=α(x,y)𝛼𝛼𝑥𝑦\alpha=\alpha(x,y)italic_α = italic_α ( italic_x , italic_y ), ranges from 0 to π/2𝜋2\pi/2italic_π / 2 such that the horizontally polarized component of the field oscillates from +11+1+ 1 to 11-1- 1 according to cos(2α)2𝛼\cos(2\alpha)roman_cos ( 2 italic_α ), while the term sin(2α)2𝛼\sin(2\alpha)roman_sin ( 2 italic_α ) oscillates from 0 to +1 and back to 0. As a consequence, the contribution from the sin(2α)2𝛼\sin(2\alpha)roman_sin ( 2 italic_α )-component, which does not have negative amplitudes, behaves like the Fresnel zones, giving rise to a focus at τf/2subscript𝜏f2\tau_{\rm f}/2italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT / 2. In Fig. 1(g), we show the expected propagation including the contribution from the sin(2α)2𝛼\sin(2\alpha)roman_sin ( 2 italic_α )-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 n0subscript𝑛0n_{0}italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT 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 n0subscript𝑛0n_{0}italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, corresponding to a tighter focusing time τn0subscript𝜏subscript𝑛0\tau_{n_{0}}italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT, 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 τn0subscript𝜏subscript𝑛0\tau_{n_{0}}italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT with n0=14subscript𝑛014n_{0}=14italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 14, 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 n0subscript𝑛0n_{0}italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is increased, the number of oscillations in ψ0(opt)superscriptsubscript𝜓0opt\psi_{0}^{\rm(opt)}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT from +11+1+ 1 to 11-1- 1 increases and in particular, the outer rings of the pattern become ever thinner.

Summary. We have derived the optimal real-valued matter wave ψ0(opt)superscriptsubscript𝜓0opt\psi_{0}^{\rm(opt)}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT 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 ψ0(opt)superscriptsubscript𝜓0opt\psi_{0}^{\rm(opt)}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT using liquid crystal devices, verifying the superior focusing properties of ψ0(opt)superscriptsubscript𝜓0opt\psi_{0}^{\rm(opt)}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT 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 ψ0(opt)superscriptsubscript𝜓0opt\psi_{0}^{\rm(opt)}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT 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

ψ0(u)=12πτf2λ[Acos(u22τf)+Bsin(u22τf)],subscript𝜓0𝑢12𝜋superscriptsubscript𝜏f2𝜆delimited-[]𝐴superscript𝑢22subscript𝜏f𝐵superscript𝑢22subscript𝜏f\psi_{0}(u)=\frac{1}{2\pi\tau_{\rm f}^{2}\lambda}\left[A\cos\left(\frac{u^{2}}% {2\tau_{\rm f}}\right)+B\sin\left(\frac{u^{2}}{2\tau_{\rm f}}\right)\right],italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_u ) = divide start_ARG 1 end_ARG start_ARG 2 italic_π italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_λ end_ARG [ italic_A roman_cos ( divide start_ARG italic_u start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT end_ARG ) + italic_B roman_sin ( divide start_ARG italic_u start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT end_ARG ) ] , (M1)

where

A01vdvcos(v22τf)ψ0(v)𝐴superscriptsubscript01𝑣differential-d𝑣superscript𝑣22subscript𝜏fsubscript𝜓0𝑣A\equiv\int\limits_{0}^{1}v\,{\rm d}v\cos\left(\frac{v^{2}}{2\tau_{\rm f}}% \right)\psi_{0}(v)italic_A ≡ ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT italic_v roman_d italic_v roman_cos ( divide start_ARG italic_v start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT end_ARG ) italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_v ) (M2)

and

B01vdvsin(v22τf)ψ0(v)𝐵superscriptsubscript01𝑣differential-d𝑣superscript𝑣22subscript𝜏fsubscript𝜓0𝑣B\equiv\int\limits_{0}^{1}v\,{\rm d}v\sin\left(\frac{v^{2}}{2\tau_{\rm f}}% \right)\psi_{0}(v)italic_B ≡ ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT italic_v roman_d italic_v roman_sin ( divide start_ARG italic_v start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT end_ARG ) italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_v ) (M3)

are functions solely of τfsubscript𝜏f\tau_{\rm f}italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT.

Next, we insert ψ0subscript𝜓0\psi_{0}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT given by Eq. (M1) into Eqs. (M2) and (M3), and obtain the system

[λ1+τfsin(1/τf)8πτf2]Asin2[1/(2τf)]4πτfB=delimited-[]𝜆1subscript𝜏f1subscript𝜏f8𝜋superscriptsubscript𝜏f2𝐴superscript212subscript𝜏f4𝜋subscript𝜏f𝐵absent\displaystyle\left[\lambda-\frac{1+\tau_{\rm f}\sin(1/\tau_{\rm f})}{8\pi\tau_% {\rm f}^{2}}\right]A-\frac{\sin^{2}[1/(2\tau_{\rm f})]}{4\pi\tau_{\rm f}}B=[ italic_λ - divide start_ARG 1 + italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT roman_sin ( 1 / italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT ) end_ARG start_ARG 8 italic_π italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ] italic_A - divide start_ARG roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [ 1 / ( 2 italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT ) ] end_ARG start_ARG 4 italic_π italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT end_ARG italic_B = 00\displaystyle 0 (M4)
sin2[1/(2τf)]4πτfA+[λ1τfsin(1/τf)8πτf2]B=superscript212subscript𝜏f4𝜋subscript𝜏f𝐴delimited-[]𝜆1subscript𝜏f1subscript𝜏f8𝜋superscriptsubscript𝜏f2𝐵absent\displaystyle-\frac{\sin^{2}[1/(2\tau_{\rm f})]}{4\pi\tau_{\rm f}}A+\left[% \lambda-\frac{1-\tau_{\rm f}\sin(1/\tau_{\rm f})}{8\pi\tau_{\rm f}^{2}}\right]B=- divide start_ARG roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [ 1 / ( 2 italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT ) ] end_ARG start_ARG 4 italic_π italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT end_ARG italic_A + [ italic_λ - divide start_ARG 1 - italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT roman_sin ( 1 / italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT ) end_ARG start_ARG 8 italic_π italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ] italic_B = 00\displaystyle 0 (M5)

of algebraic equations for A𝐴Aitalic_A and B𝐵Bitalic_B, which has non-trivial solutions, only when its determinant is zero, that is

(λ18πτf2)2(sin[1/(2τf)]4πτf)2=0.superscript𝜆18𝜋superscriptsubscript𝜏f22superscript12subscript𝜏f4𝜋subscript𝜏f20\left(\lambda-\frac{1}{8\pi\tau_{\rm f}^{2}}\right)^{2}-\left(\frac{\sin[1/(2% \tau_{\rm f})]}{4\pi\tau_{\rm f}}\right)^{2}=0.( italic_λ - divide start_ARG 1 end_ARG start_ARG 8 italic_π italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - ( divide start_ARG roman_sin [ 1 / ( 2 italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT ) ] end_ARG start_ARG 4 italic_π italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 0 . (M6)

This elementary quadratic equation has the two solutions

λ±=18πτf2[1±2τf|sin(12τf)|].subscript𝜆plus-or-minus18𝜋superscriptsubscript𝜏f2delimited-[]plus-or-minus12subscript𝜏f12subscript𝜏f\lambda_{\pm}=\frac{1}{8\pi\tau_{\rm f}^{2}}\left[1\pm 2\tau_{\rm f}\left|\sin% \left(\frac{1}{2\tau_{\rm f}}\right)\right|\,\right].italic_λ start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 8 italic_π italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ 1 ± 2 italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT | roman_sin ( divide start_ARG 1 end_ARG start_ARG 2 italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT end_ARG ) | ] . (M7)

By inserting the maximal eigenvalue λ+subscript𝜆\lambda_{+}italic_λ start_POSTSUBSCRIPT + end_POSTSUBSCRIPT into Eq. (M4), we find the relation between A𝐴Aitalic_A and B𝐵Bitalic_B, and thus the normalized optimal initial wave function ψ0(opt)superscriptsubscript𝜓0opt\psi_{0}^{\rm(opt)}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT given by Eq. (7).

.2 Amplitude and phase Fresnel zones:
Maximal intensity

For a given value of τfsubscript𝜏f\tau_{\rm f}italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT the radii

ρn2πτfnsubscript𝜌𝑛2𝜋subscript𝜏f𝑛\rho_{n}\equiv\sqrt{2\pi\tau_{\rm f}n}italic_ρ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ≡ square-root start_ARG 2 italic_π italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT italic_n end_ARG (M8)

of the Fresnel zones, with n=1,2,3,,n0𝑛123subscript𝑛0n=1,2,3,\ldots,n_{0}italic_n = 1 , 2 , 3 , … , italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, extend to the maximum number n0subscript𝑛0n_{0}italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT of zones fitting within the circular aperture 0ρ10𝜌10\leq\rho\leq 10 ≤ italic_ρ ≤ 1 [1]. Therefore, the initial wave functions ψ0(AFZ)superscriptsubscript𝜓0AFZ\psi_{0}^{\rm(AFZ)}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_AFZ ) end_POSTSUPERSCRIPT and ψ0(PFZ)superscriptsubscript𝜓0PFZ\psi_{0}^{\rm(PFZ)}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_PFZ ) end_POSTSUPERSCRIPT for the amplitude and phase Fresnel zone patterns read

ψ0(AFZ)(ρ)=N1[𝕌(ρ;0,ρ1)+n=1𝕌(ρ;ρ2n,ρ2n+1)]superscriptsubscript𝜓0AFZ𝜌subscript𝑁1delimited-[]𝕌𝜌0subscript𝜌1superscriptsubscript𝑛1𝕌𝜌subscript𝜌2𝑛subscript𝜌2𝑛1\psi_{0}^{\rm(AFZ)}(\rho)=N_{1}\left[{\mathbb{U}}(\rho;0,\rho_{1})+\sum_{n=1}^% {\infty}{\mathbb{U}}(\rho;\rho_{2n},\rho_{2n+1})\right]italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_AFZ ) end_POSTSUPERSCRIPT ( italic_ρ ) = italic_N start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ blackboard_U ( italic_ρ ; 0 , italic_ρ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) + ∑ start_POSTSUBSCRIPT italic_n = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT blackboard_U ( italic_ρ ; italic_ρ start_POSTSUBSCRIPT 2 italic_n end_POSTSUBSCRIPT , italic_ρ start_POSTSUBSCRIPT 2 italic_n + 1 end_POSTSUBSCRIPT ) ] (M9)

and

ψ0(PFZ)(ρ)=1π[𝕌(ρ;0,ρ1)+n=1(1)n𝕌(ρ;ρn,ρn+1)]superscriptsubscript𝜓0PFZ𝜌1𝜋delimited-[]𝕌𝜌0subscript𝜌1superscriptsubscript𝑛1superscript1𝑛𝕌𝜌subscript𝜌𝑛subscript𝜌𝑛1\psi_{0}^{\rm(PFZ)}(\rho)=\frac{1}{\sqrt{\pi}}\left[{\mathbb{U}}(\rho;0,\rho_{% 1})+\sum_{n=1}^{\infty}(-1)^{n}{\mathbb{U}}(\rho;\rho_{n},\rho_{n+1})\right]italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_PFZ ) end_POSTSUPERSCRIPT ( italic_ρ ) = divide start_ARG 1 end_ARG start_ARG square-root start_ARG italic_π end_ARG end_ARG [ blackboard_U ( italic_ρ ; 0 , italic_ρ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) + ∑ start_POSTSUBSCRIPT italic_n = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT ( - 1 ) start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT blackboard_U ( italic_ρ ; italic_ρ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT , italic_ρ start_POSTSUBSCRIPT italic_n + 1 end_POSTSUBSCRIPT ) ] (M10)

with 𝕌(ρ;a,b)Θ(bρ)Θ(aρ)𝕌𝜌𝑎𝑏Θ𝑏𝜌Θ𝑎𝜌{\mathbb{U}}(\rho;a,b)\equiv\Theta(b-\rho)-\Theta(a-\rho)blackboard_U ( italic_ρ ; italic_a , italic_b ) ≡ roman_Θ ( italic_b - italic_ρ ) - roman_Θ ( italic_a - italic_ρ ) and a<b𝑎𝑏a<bitalic_a < italic_b. Here, Θ(ρ)Θ𝜌\Theta(\rho)roman_Θ ( italic_ρ ) denotes the Heaviside function and N1subscript𝑁1N_{1}italic_N start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT is a normalization constant depending on τfsubscript𝜏f\tau_{\rm f}italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT.

To derive an analytical formula for the maximal intensity, we choose the focusing times τfτn01/(2πn0)subscript𝜏fsubscript𝜏subscript𝑛012𝜋subscript𝑛0\tau_{\rm f}\equiv\tau_{n_{0}}\equiv 1/(2\pi n_{0})italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT ≡ italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ≡ 1 / ( 2 italic_π italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ). In this case, Eq. (M8) reduces to ρn=n/n0subscript𝜌𝑛𝑛subscript𝑛0\rho_{n}=\sqrt{n/n_{0}}italic_ρ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT = square-root start_ARG italic_n / italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG, and the normalization condition

πN12[ρ12+(ρ32ρ22)+(ρ52ρ42)+]=1𝜋superscriptsubscript𝑁12delimited-[]superscriptsubscript𝜌12superscriptsubscript𝜌32superscriptsubscript𝜌22superscriptsubscript𝜌52superscriptsubscript𝜌421\pi N_{1}^{2}\left[\rho_{1}^{2}+\left(\rho_{3}^{2}-\rho_{2}^{2}\right)+\left(% \rho_{5}^{2}-\rho_{4}^{2}\right)+\ldots\right]=1italic_π italic_N start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [ italic_ρ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_ρ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ρ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) + ( italic_ρ start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ρ start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) + … ] = 1 (M11)

for ψ0(AFZ)superscriptsubscript𝜓0AFZ\psi_{0}^{\rm(AFZ)}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_AFZ ) end_POSTSUPERSCRIPT, Eq. (M9), defines the constant

N12π{n0n0+1,n0=1,3,5,1,n0=2,4,6,subscript𝑁12𝜋casessubscript𝑛0subscript𝑛01subscript𝑛01351subscript𝑛0246N_{1}\equiv\sqrt{\frac{2}{\pi}}\begin{cases}\sqrt{\frac{n_{0}}{n_{0}+1}},&% \quad n_{0}=1,3,5,\ldots\\ 1,&\quad n_{0}=2,4,6,\ldots\end{cases}italic_N start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ≡ square-root start_ARG divide start_ARG 2 end_ARG start_ARG italic_π end_ARG end_ARG { start_ROW start_CELL square-root start_ARG divide start_ARG italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + 1 end_ARG end_ARG , end_CELL start_CELL italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1 , 3 , 5 , … end_CELL end_ROW start_ROW start_CELL 1 , end_CELL start_CELL italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 2 , 4 , 6 , … end_CELL end_ROW (M12)

as a function of n0subscript𝑛0n_{0}italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT.

Next, we insert the initial profile ψ0(AFZ)superscriptsubscript𝜓0AFZ\psi_{0}^{\rm(AFZ)}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_AFZ ) end_POSTSUPERSCRIPT given by Eq. (M9) into Eq. (3), and obtain the expression

Imax(AFZ)(τn0)=1τn02|01uduexp(iu22τn0)ψ0(AFZ)(u)|2superscriptsubscript𝐼maxAFZsubscript𝜏subscript𝑛01superscriptsubscript𝜏subscript𝑛02superscriptsuperscriptsubscript01𝑢differential-d𝑢isuperscript𝑢22subscript𝜏subscript𝑛0superscriptsubscript𝜓0AFZ𝑢2I_{\rm max}^{\rm(AFZ)}(\tau_{n_{0}})=\frac{1}{\tau_{n_{0}}^{2}}\left|\int% \limits_{0}^{1}u{\rm d}u\,\exp\left({\rm i}\frac{u^{2}}{2\tau_{n_{0}}}\right)% \psi_{0}^{\rm(AFZ)}(u)\right|^{2}italic_I start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_AFZ ) end_POSTSUPERSCRIPT ( italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) = divide start_ARG 1 end_ARG start_ARG italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG | ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT italic_u roman_d italic_u roman_exp ( roman_i divide start_ARG italic_u start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT end_ARG ) italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_AFZ ) end_POSTSUPERSCRIPT ( italic_u ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT
=N12|exp(iρ122τn0)1+exp(iρ322τn0)exp(iρ222τn0)+|2,absentsuperscriptsubscript𝑁12superscriptisuperscriptsubscript𝜌122subscript𝜏subscript𝑛01isuperscriptsubscript𝜌322subscript𝜏subscript𝑛0isuperscriptsubscript𝜌222subscript𝜏subscript𝑛02=N_{1}^{2}\left|\exp\left({\rm i}\frac{\rho_{1}^{2}}{2\tau_{n_{0}}}\right)-1+% \exp\left({\rm i}\frac{\rho_{3}^{2}}{2\tau_{n_{0}}}\right)-\exp\left({\rm i}% \frac{\rho_{2}^{2}}{2\tau_{n_{0}}}\right)+\ldots\right|^{2},= italic_N start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT | roman_exp ( roman_i divide start_ARG italic_ρ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT end_ARG ) - 1 + roman_exp ( roman_i divide start_ARG italic_ρ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT end_ARG ) - roman_exp ( roman_i divide start_ARG italic_ρ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT end_ARG ) + … | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ,

that is

Imax(AFZ)(τn0)=N12{(n0+1)2,n0=1,3,5,n02,n0=2,4,6,,superscriptsubscript𝐼maxAFZsubscript𝜏subscript𝑛0superscriptsubscript𝑁12casessuperscriptsubscript𝑛012subscript𝑛0135superscriptsubscript𝑛02subscript𝑛0246I_{\rm max}^{\rm(AFZ)}(\tau_{n_{0}})=N_{1}^{2}\begin{cases}(n_{0}+1)^{2},&% \quad n_{0}=1,3,5,\ldots\\ n_{0}^{2},&\quad n_{0}=2,4,6,\ldots,\end{cases}italic_I start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_AFZ ) end_POSTSUPERSCRIPT ( italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) = italic_N start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT { start_ROW start_CELL ( italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , end_CELL start_CELL italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1 , 3 , 5 , … end_CELL end_ROW start_ROW start_CELL italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , end_CELL start_CELL italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 2 , 4 , 6 , … , end_CELL end_ROW (M13)

where we have used the fact that ρn2/(2τn0)=nπsuperscriptsubscript𝜌𝑛22subscript𝜏subscript𝑛0𝑛𝜋\rho_{n}^{2}/(2\tau_{n_{0}})=n\piitalic_ρ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / ( 2 italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) = italic_n italic_π.

As a result, Eq. (M13) combined with Eq. (M12) gives rise to the maximum intensity Imax(AFZ)(τn0)superscriptsubscript𝐼maxAFZsubscript𝜏subscript𝑛0I_{\rm max}^{\rm(AFZ)}(\tau_{n_{0}})italic_I start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_AFZ ) end_POSTSUPERSCRIPT ( italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ), Eq. (10), produced by the amplitude Fresnel zones. Analogously, we derive the corresponding maximum intensity Imax(PFZ)(τn0)superscriptsubscript𝐼maxPFZsubscript𝜏subscript𝑛0I_{\rm max}^{\rm(PFZ)}(\tau_{n_{0}})italic_I start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_PFZ ) end_POSTSUPERSCRIPT ( italic_τ start_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ), 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 φ0φ0(x)subscript𝜑0subscript𝜑0𝑥\varphi_{0}\equiv\varphi_{0}(x)italic_φ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ≡ italic_φ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) that maximizes the intensity |φ(0)|2superscript𝜑02\left|\varphi(0)\right|^{2}| italic_φ ( 0 ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT of the field at x=0𝑥0x=0italic_x = 0, at the focusing time tfsubscript𝑡ft_{\rm f}italic_t start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT.

In this case we use the one-dimensional Green function

G(1)(ξ,τ|ξ,0)=12πiτexp[i(ξξ)22τ]superscriptG1𝜉conditional𝜏superscript𝜉012𝜋i𝜏isuperscript𝜉superscript𝜉22𝜏{\rm G}^{\rm(1)}(\xi,\tau|\xi^{\prime},0)=\frac{1}{\sqrt{2\pi{\rm i}\tau}}\exp% \left[{\rm i}\frac{(\xi-\xi^{\prime})^{2}}{2\tau}\right]roman_G start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT ( italic_ξ , italic_τ | italic_ξ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , 0 ) = divide start_ARG 1 end_ARG start_ARG square-root start_ARG 2 italic_π roman_i italic_τ end_ARG end_ARG roman_exp [ roman_i divide start_ARG ( italic_ξ - italic_ξ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_τ end_ARG ] (M14)

for the time-dependent one-dimensional Schrödinger equation of a free particle. Here, ξx/L𝜉𝑥𝐿\xi\equiv x/Litalic_ξ ≡ italic_x / italic_L and τt/(ML2)𝜏Planck-constant-over-2-pi𝑡𝑀superscript𝐿2\tau\equiv\hbar t/(ML^{2})italic_τ ≡ roman_ℏ italic_t / ( italic_M italic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) are the dimensionless position and time, respectively, and M𝑀Mitalic_M and L𝐿Litalic_L 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

12πτf11dξcos(ξ2ξ22τf)φ0(ξ)=μφ0(ξ)12𝜋subscript𝜏fsuperscriptsubscript11differential-dsuperscript𝜉superscript𝜉2superscript𝜉22subscript𝜏fsubscript𝜑0superscript𝜉𝜇subscript𝜑0𝜉\frac{1}{2\pi\tau_{\rm f}}\int\limits_{-1}^{1}{\rm d}\xi^{\prime}\cos\left(% \frac{\xi^{2}-\xi^{\prime 2}}{2\tau_{\rm f}}\right)\varphi_{0}(\xi^{\prime})=% \mu\varphi_{0}(\xi)divide start_ARG 1 end_ARG start_ARG 2 italic_π italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT end_ARG ∫ start_POSTSUBSCRIPT - 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT roman_d italic_ξ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT roman_cos ( divide start_ARG italic_ξ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ξ start_POSTSUPERSCRIPT ′ 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT end_ARG ) italic_φ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_ξ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) = italic_μ italic_φ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_ξ ) (M15)

for the optimal initial wave function φ0subscript𝜑0\varphi_{0}italic_φ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT with the eigenvalue μ𝜇\muitalic_μ, that determines the maximum intensity achieved at τfsubscript𝜏f\tau_{\rm f}italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT. Here, we have assumed that φ0(ξ)=0subscript𝜑0𝜉0\varphi_{0}(\xi)=0italic_φ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_ξ ) = 0 for |ξ|>1𝜉1|\xi|>1| italic_ξ | > 1.

Since we are interested in the maximum of the intensity, we solve the integral equation (M15) only for the largest eigenvalue μ+(τf)subscript𝜇subscript𝜏f\mu_{+}(\tau_{\rm f})italic_μ start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ( italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT ). As a result, for a given τfsubscript𝜏f\tau_{\rm f}italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT, we find the optimal initial wave function

φ0(opt)(ξ)=1+r4πτfμ+cos(ξ22τf)+1r4πτfμ+sin(ξ22τf)superscriptsubscript𝜑0opt𝜉1𝑟4𝜋subscript𝜏fsubscript𝜇superscript𝜉22subscript𝜏f1𝑟4𝜋subscript𝜏fsubscript𝜇superscript𝜉22subscript𝜏f\varphi_{0}^{\rm(opt)}(\xi)=\sqrt{\frac{1+r}{4\pi\tau_{\rm f}\mu_{+}}}\cos% \left(\frac{\xi^{2}}{2\tau_{\rm f}}\right)+\sqrt{\frac{1-r}{4\pi\tau_{\rm f}% \mu_{+}}}\sin\left(\frac{\xi^{2}}{2\tau_{\rm f}}\right)italic_φ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT ( italic_ξ ) = square-root start_ARG divide start_ARG 1 + italic_r end_ARG start_ARG 4 italic_π italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT italic_μ start_POSTSUBSCRIPT + end_POSTSUBSCRIPT end_ARG end_ARG roman_cos ( divide start_ARG italic_ξ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT end_ARG ) + square-root start_ARG divide start_ARG 1 - italic_r end_ARG start_ARG 4 italic_π italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT italic_μ start_POSTSUBSCRIPT + end_POSTSUBSCRIPT end_ARG end_ARG roman_sin ( divide start_ARG italic_ξ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT end_ARG ) (M16)

for |ξ|1𝜉1|\xi|\leq 1| italic_ξ | ≤ 1, with φ0(opt)(ξ)=0superscriptsubscript𝜑0opt𝜉0\varphi_{0}^{\rm(opt)}(\xi)=0italic_φ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT ( italic_ξ ) = 0 for |ξ|>1𝜉1|\xi|>1| italic_ξ | > 1, and the corresponding maximal eigenvalue

μ+(τf)=Imax(1D)(τf)=[2/(πτf)].subscript𝜇subscript𝜏fsuperscriptsubscript𝐼max1Dsubscript𝜏fdelimited-[]2𝜋subscript𝜏f\mu_{+}(\tau_{\rm f})=I_{\rm max}^{\rm(1D)}(\tau_{\rm f})=\mathcal{I}\left[% \sqrt{2/(\pi\tau_{\rm f})}\right].italic_μ start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ( italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT ) = italic_I start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 roman_D ) end_POSTSUPERSCRIPT ( italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT ) = caligraphic_I [ square-root start_ARG 2 / ( italic_π italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT ) end_ARG ] . (M17)

Here we have expressed the intensity

(z)14{z2+z[C(z)]2+[S(z)]2}𝑧14superscript𝑧2𝑧superscriptdelimited-[]𝐶𝑧2superscriptdelimited-[]𝑆𝑧2\mathcal{I}(z)\equiv\frac{1}{4}\left\{z^{2}+z\sqrt{\left[C\left(z\right)\right% ]^{2}+\left[S\left(z\right)\right]^{2}}\right\}caligraphic_I ( italic_z ) ≡ divide start_ARG 1 end_ARG start_ARG 4 end_ARG { italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_z square-root start_ARG [ italic_C ( italic_z ) ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + [ italic_S ( italic_z ) ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG }

and the parameter

r(τf)C[2/(πτf)]{C[2/(πτf)]}2+{S[2/(πτf)]}2𝑟subscript𝜏f𝐶delimited-[]2𝜋subscript𝜏𝑓superscript𝐶delimited-[]2𝜋subscript𝜏𝑓2superscript𝑆delimited-[]2𝜋subscript𝜏𝑓2r(\tau_{\rm f})\equiv\frac{C\left[\sqrt{2/(\pi\tau_{f})}\,\right]}{\sqrt{\left% \{C\left[\sqrt{2/(\pi\tau_{f})}\,\right]\right\}^{2}+\left\{S\left[\sqrt{2/(% \pi\tau_{f})}\,\right]\right\}^{2}}}italic_r ( italic_τ start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT ) ≡ divide start_ARG italic_C [ square-root start_ARG 2 / ( italic_π italic_τ start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT ) end_ARG ] end_ARG start_ARG square-root start_ARG { italic_C [ square-root start_ARG 2 / ( italic_π italic_τ start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT ) end_ARG ] } start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + { italic_S [ square-root start_ARG 2 / ( italic_π italic_τ start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT ) end_ARG ] } start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG (M18)

in terms of the Fresnel integrals C𝐶Citalic_C and S𝑆Sitalic_S [7].

.4 Pancharatnam-Berry optical element

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

𝒰^qsubscript^𝒰𝑞\displaystyle\hat{\mathcal{U}}_{q}over^ start_ARG caligraphic_U end_ARG start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT \displaystyle\cdot (𝐞H𝐞V)=cos(δ2)(𝐞H𝐞V)matrixsubscript𝐞Hsubscript𝐞V𝛿2matrixsubscript𝐞Hsubscript𝐞V\displaystyle\begin{pmatrix}\mathbf{e}_{\rm H}\\ \mathbf{e}_{\rm V}\end{pmatrix}=\cos\left(\frac{\delta}{2}\right)\begin{% pmatrix}\mathbf{e}_{\rm H}\\ \mathbf{e}_{\rm V}\end{pmatrix}( start_ARG start_ROW start_CELL bold_e start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL bold_e start_POSTSUBSCRIPT roman_V end_POSTSUBSCRIPT end_CELL end_ROW end_ARG ) = roman_cos ( divide start_ARG italic_δ end_ARG start_ARG 2 end_ARG ) ( start_ARG start_ROW start_CELL bold_e start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL bold_e start_POSTSUBSCRIPT roman_V end_POSTSUBSCRIPT end_CELL end_ROW end_ARG )
+\displaystyle++ isin(δ2)(cos[2α(x,y)]sin[2α(x,y)]sin[2α(x,y)]cos[2α(x,y)])(𝐞H𝐞V)𝑖𝛿2matrix2𝛼𝑥𝑦2𝛼𝑥𝑦2𝛼𝑥𝑦2𝛼𝑥𝑦matrixsubscript𝐞Hsubscript𝐞V\displaystyle i\sin\left(\frac{\delta}{2}\right)\begin{pmatrix}\cos\left[2% \alpha(x,y)\right]\;\;\;\;\;\sin\left[2\alpha(x,y)\right]\\ \sin\left[2\alpha(x,y)\right]\;-\cos\left[2\alpha(x,y)\right]\end{pmatrix}% \begin{pmatrix}\mathbf{e}_{\rm H}\\ \mathbf{e}_{\rm V}\end{pmatrix}italic_i roman_sin ( divide start_ARG italic_δ end_ARG start_ARG 2 end_ARG ) ( start_ARG start_ROW start_CELL roman_cos [ 2 italic_α ( italic_x , italic_y ) ] roman_sin [ 2 italic_α ( italic_x , italic_y ) ] end_CELL end_ROW start_ROW start_CELL roman_sin [ 2 italic_α ( italic_x , italic_y ) ] - roman_cos [ 2 italic_α ( italic_x , italic_y ) ] end_CELL end_ROW end_ARG ) ( start_ARG start_ROW start_CELL bold_e start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL bold_e start_POSTSUBSCRIPT roman_V end_POSTSUBSCRIPT end_CELL end_ROW end_ARG )

on the horizontal 𝐞Hsubscript𝐞H{\bf e}_{\rm H}bold_e start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT and vertical 𝐞Vsubscript𝐞V{\bf e}_{\rm V}bold_e start_POSTSUBSCRIPT roman_V end_POSTSUBSCRIPT polarization components of an optical beam. Here δ𝛿\deltaitalic_δ is the optical retardation of the liquid crystal molecules and αα(x,y)𝛼𝛼𝑥𝑦\alpha\equiv\alpha(x,y)italic_α ≡ italic_α ( italic_x , italic_y ) is the device’s spatially dependent liquid crystal axis orientation expressed in terms of the transverse Cartesian coordinates x𝑥xitalic_x and y𝑦yitalic_y.

When the device is perfectly tuned, that is for δ=π𝛿𝜋\delta=\piitalic_δ = italic_π, 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 cos[2α(x,y)]2𝛼𝑥𝑦\cos\left[2\alpha(x,y)\right]roman_cos [ 2 italic_α ( italic_x , italic_y ) ]. This procedure was employed to generate our real-valued optimal state ψ0(opt)ψ0(opt)(x,y)superscriptsubscript𝜓0optsuperscriptsubscript𝜓0opt𝑥𝑦\psi_{0}^{\rm(opt)}\equiv\psi_{0}^{\rm(opt)}(x,y)italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT ≡ italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT ( italic_x , italic_y ) by means of a device defined by an optical axis of α(x,y)=(1/2)arccos[ψ0(opt)(x,y)/ψmax(opt)]𝛼𝑥𝑦12superscriptsubscript𝜓0opt𝑥𝑦superscriptsubscript𝜓maxopt\alpha(x,y)=(1/2)\arccos\left[\psi_{0}^{\rm(opt)}(x,y)/\psi_{\rm max}^{\rm(opt% )}\right]italic_α ( italic_x , italic_y ) = ( 1 / 2 ) roman_arccos [ italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT ( italic_x , italic_y ) / italic_ψ start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT ], where ψmax(opt)superscriptsubscript𝜓maxopt\psi_{\rm max}^{\rm(opt)}italic_ψ start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( roman_opt ) end_POSTSUPERSCRIPT 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. A 65, 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 s𝑠sitalic_s-wave. Phys. Rev. Lett. 89, 060404 (2002).
  • [4] Case, W. B., Sadurni, E. & Schleich, W. P. A diffractive mechanism of focusing. Opt. Express 20, 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 Optics 18, 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 Photonics 2, 501–505 (2008).
  • [12] Karimi, E., Piccirillo, B., Marrucci, L. & Santamato, E. Improved focusing with hypergeometric-gaussian type-ii optical modes. Opt. Express 16, 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 A&\&&M University for a Faculty Fellowship, and to Texas A&\&&M 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.