X-ray vortex beams : A theoretical analysis

The recent demonstration that an optical vortex could be generated at x-ray wavelengths brings this interesting topological phenomenon into an entirely new regime with several possible applications. We examine the analytic propagation of an optical vortex generated in a synchrotron x-ray beam line. We compare the results obtained with the existing experimental data and further consider the generation and interpretation of mixed vortex-edge discontinuities which might be considered as non-integer charge vortices. ©2003 Optical Society of America OCIS codes: (340.6720) Synchrotron radiation; (350.7420) Waves References and Links 1. J. F. Nye and M. V. Berry, “Dislocations in wave trains,” Proc. R. Soc. London, Ser. A 336, 165 – 190 (1974). 2. N. B. Baranova, B. Ya, Zel’dovich, A. V. Mamayev, N. F. Pilipetskii, and V. V. Shkukov, “Dislocations of the wavefront of a speckle-inhomogeneous field (theory and experiment),” Pis’ma Zh. Eks. Teor. Fiz. 33, 206 – 210 (1981) [JETP Lett. 33, 195 – 199 (1981)]. 3. Z. S. Sacks, D. Rozas, and G. A. Swartzlander, Jr., “Holographic formation of optical-vortex filaments,” J. Opt. Soc. Am. B 15, 2226 – 2234 (1998). 4. M. S. Naschie, ed., “Special issue on nonlinear optical structures, patterns, chaos,” Chaos Solitons Fractals 4(8/9) (1994.). 5. V. Yu. Bazhenov, M. S. Soskin and M. V. Vasnetsov, “Screw dislocations in light wavefronts,” J. Mod. Opt. 39, 985 990 (1992). 6. I. Freund, “Critical point explosions in two-dimensional wave fields,” Opt. Commun. 159, 99 117(1999). 7. G. S. Agarwal and J. Banerji, “Spatial coherence and information entropy in optical vortex fields,” Opt. Lett. 27, 800 – 802 (2002). 8. N. B. Simpson, L. Allen and M. J. Padgett, “Optical tweezers and optical spanners with Laguerre-Gaussian modes,” J. Mod. Opt. 43, 2485 – 2491 (1996). 9. E. Wolf, Progress in optics 42, (Elsevier, 2001), 10. M. Harris, “Light-field fluctuations in space and time,” Contemporary Phys. 36, 215 – 233 (1995). 11. M. V. Vasnetsov, I. G. Marienko and M. S. Soskin, “Self-reconstruction of an optical vortex” JETP Lett. 71, 130 – 133 (2000). 12. Z. Bouchal, “Resistance of nondiffracting vortex beam against amplitude and phase pertubations,” Opt. Commun. 210, 155 – 164 (2002). 13. D. Rozas, Z. S. Sacks, and G. A. Swartzlander, Jr., “Experimental observation of fluidlike motion of optical vortices,” Phys. Rev. Lett. 79, 3399 – 3402 (1997). 14. D. Rozas and G. A. Swartzlander, Jr., “Observed rotational enhancement of nonlinear optical vorticies,” Opt. Lett. 25, 126 – 128 (2000). 15. M. Brambilla, F. Battipede, L. A. Lugiato, V. Penna, F. Prati, C. Tamm, C. O. Weiss, “Transverse laser patterns.m I. Phase singularity crystals,” Phys. Rev. A 43, 5090 – 5117 (1991) 16. W. Whewell, “ Essay towards a first approximation to a map of cotidal lines,” Phil. Trans. R. Soc. Lond. 123, 147 – 236 (1833). 17. M. W. Beijersbergen, L. Allen, H. E. L. O. van der Veen, and J. P. Woerdman, “Astigmatic laser mode converters and transfer of orbital angular momentum,“ Opt. Commun. 96, 123 – 132 (1993). 18. N. R. Heckenberg, R. McDuff, C. P. Smith, and A. G. White, “Generation of optical phase singularities by computer-generated holograms,” Opt. Lett. 17, 221 – 223 (1992). 19. M. W. Beijersbergen, R. P. C. Coerwinkel, M. Kristensen, and J. P. Woerdman, „Helical-wavefront laser beams produced with a spiral phaseplate,“ Opt. Commun. 112, 321 – 327 (1994). 20. G. A. Turnball, D. A. Robertson, G. M. Smith, L. Allen, and M. J. Padgett, “The generation of free-space Laguerre-Gaussian modes at millimetre-wave frequencies by use of a spiral phaseplate,” Opt. Commun. 127, 183 – 188 (1996). 21. G. H. Kim, J. H. Jeon, K. H. Ko, H. J. Moon, J. H. Lee, and J. S. Chang, “Optical vortices produced with a nonspiral phase plate,” Appl. Opt. 36, 8614 – 8621 (1997). (C) 2003 OSA 22 September 2003 / Vol. 11, No. 19 / OPTICS EXPRESS 2315 #2793 $15.00 US Received July 28, 2003; Revised August 25, 2003 22. A. G. Peele, P. J. McMahon, D. Paterson, C. Q. Tran, A. P. Mancuso, K. A. Nugent, J. P. Hayes, E. C. Harvey, B. Lai, I. McNulty, “Observation of an x-ray vortex,” Opt. Lett. 27, 1752 – 1754 (2002) 23. G. A. Swartzlander, Jr., “Peering into darkness with a vortex spatial filter,” Opt. Lett. 26, 497 – 499 (2001). 24. M. D. Levenson, G. Dai, T. Ebihara, "Vortex Mask: Making 80nm contacts with a twist!" Proc. SPIE 4889, 1293-1303 (2002). 25. K. A. Nugent, T. E. Gureyev, D. F. Cookson, D. Paganin, and Z. Barnea, “Quantitative phase imaging using hard x rays,” Phys. Rev. Lett. 77, 2961 – 2964 (1996). 26. L. J. Allen, H. M. L. Faulkner, K. A. Nugent, M. P. Oxley, and D. Paganin, “Phase retrieval from images in the presence of first-order vortices,” Phys. Rev. E 63, 037602 (2001). 27. D. Dragoman, “Unambiguous coherence retrieval from intensity measurements,” J. Opt. Soc. Am. A 20, 290 – 295 (2003). 28. D. Paterson, B. E. Allman, P. J. McMahon, J. Lin, N. Moldovan, K. A. Nugent, I. McNulty, C. T. Chantler, C. C. Retsch, T. H. K. Irving, D. C. Mancini, “Spatial coherence measurement of X-ray undulator radiation,” Opt. Commun. 195, 79 – 84 (2001). 29. G. Indebetouw, “Optical vortices and their propagation,” J. Mod. Opt. 40, 73 – 87 (1993). 30. I. V. Basistiy, V. Yu Bazhenov, M. S. Soskin, and M. V. Vasnetov, “Optics of light beams with screw dislocations,” Opt. Commun. 103, 422 – 428 (1993). 31. D. Rozas, C. T. Law, and G. A. Swartzlander, Jr., “Propagation dynamics of optical vortices,” J. Opt. Soc. Am. B 14, 3054 – 3065 (1997). 32. U. T. Schwarz, S. Sogomonian, M. Maier, “Propagation dynamics of phase dislocations embedded in a Bessel light beam,” Opt. Commun. 208, 255 – 262 (2002). 33. S. Orlov, K.Regelskis, V. Smilgevičius, and A. Stabinis, “Propagation of Bessel beams carrying optical vortices,” Opt. Commun. 209, 155 – 165 (2002). 34. V. Pyragaite and A. Stabinis, “Free-space propagation of overlapping light vortex beams,” Opt. Commun. 213, 187 – 191 (2002). 35. R. P. Singh and S. R. Chowdhury, “Trajectory of an optical vortex: canonical vs. non-canonical,” Opt. Commun. 215, 231 – 237 (2003). 36. J. Masajada, “Half-plane diffraction in the case of Gaussian beams containing an optical vortex,” Opt. Commun. 175, 289 – 294 (2000). 37. I. V. Basistiy, M. S. Soskin, and M. V. Vasnetsov, “Optical wavefront dislocations and their properties,” Opt. Commun. 119, 604 – 612 (1995). 38. D. V. Petrov, “Vortex-edge dislocation interaction in a linear medium,” Opt. Commun. 188, 307 – 312 (2001). 39. A. E. Siegman, An Introduction to Lasers and Masers, (McGraw-Hill, 1971). 40. I. S. Gradshteyn and I. M. Ryzhik, Table of integrals, series, and products, 5 ed., (Academic Press,1994).


Introduction
The ubiquity of wavefront dislocations, including point singularities in the transverse phase commonly referred to as screw-type dislocations or optical vortices, is, following the analysis of Nye and Berry [1], now generally recognized [2,3].Likewise, optical vortices have exercised the imagination of the optics community to a considerable degree over the past decade [4].Potential applications include wave-guides for atom optics [5], optical computing [6], optical communication [7], optical tweezers, optical spanners [8], high-resolution imaging [9], and probing scattering media [10].Optical vortices also exhibit the ability to regenerate their initial phase distribution after propagation through an obstacle [11,12] which offers promise in micromanipulation and in lithography.Several studies [13,14] have also been undertaken because of insights that may be gained into vortex systems in media other than light such as fluids, superfluids and superconductors [15].Indeed, one of the earliest descriptions of a vortex system appears in a study on amphidronic tidal points [16].Optical vortices will be generated when coherent light is scattered by a rough medium [10], hence their ubiquity.While, for the purposes of investigation, controlled production via laser cavity modes [15], cylindrical lens mode conversion of laser modes [17], computer generated holograms -including spiral zone plates [18] and forked diffraction gratings [5] -and phase plates, both spiral [19,20] and non-spiral [21] have all been demonstrated.
A recent publication made the perhaps surprising demonstration that production and detection of optical vortices was possible at x-ray wavelengths [22].This opens up a whole new parameter space of possible application for optical vortices.Possible x-ray applications include using the dark vortex core to mask a bright source, as in an astronomical x-ray source [23], or in a synchrotron-based small angle scattering experiment or, as a component in an x-ray phase-shifting mask, where the central zero in the vortex intensity pattern defines a small dot in the negative photoresist that creates a contact that can be used to connect layers [24].Another idea might be to take advantage of elemental characteristic absorption energies to allow the selective transfer of orbital angular momentum thus allowing the manipulation of components of a structure that contain elements matching the tuned beam energy -for instance microfabricated objects carrying certain metal compounds may be moved while leaving their surroundings unaffected.
The stability of optical vortices on propagation makes them particularly problematic in the context of non-interferometric phase measurements as these fail where there are zeros in the intensity or where there is rotational symmetry in the phase [25,26].It appears that singularities also play a significant role in certain methods to recover correlations in the wavefield.Techniques that use only measurements of intensity will fail when the intensity distribution is unchanged after propagation through a rotationally symmetric optical system, which is the case for a vortex beam [27].Both these types of measurements are of particular interest in the x-ray regime [28] and in the context of synchrotron experiments where experimental difficulties are greatly alleviated by minimal instrumentation.
The propagation of a vortex beam has been considered analytically elsewhere [5,29 -35].However, while some of the results of these analyses scale readily to x-ray wavelengths, there are some features, particularly in the context of a synchrotron experiment, which are not dealt with.For instance, all analyses that consider propagation of an optical vortex implanted in a Gaussian beam take the vortex as implanted at the beam waist, allowing the initial phase in the beam to be zero.A feature of the beam from a third-generation synchrotron source is high coherence and low divergence.We wish to retain this feature in a vortex experiment and, if we model the flux distribution as approximately Gaussian, then the beam waist occurs at the source position, which is at the electrons in the insertion device, and is therefore inaccessible.In practice the vortex is implanted by insertion of a phase plate some tens of metres [22] away from the source, consequently the phase terms in the beam at the phase plate must be retained.We consider this case in Section 2 below.In addition, there has been little analysis of the interaction of a propagated vortex beam with an object.The case of a Gaussian beam vortex generated by a diffraction grating, imaged by a lens and diffracted by a half plane or slit has been presented [36].Also vortex-edge interactions where both the vortex and the edge dislocation are implanted simultaneously in the beam have been examined [37,38].In Section 3 we analyse the Fresnel diffraction of a propagated x-ray vortex beam diffracted from a wire and compare with numerical and experimental results.In Section 4 we use the results of the previous sections to follow the evolution in the phase structure of a vortex beam as a function of energy.

Vortex beam propagation
Based on typical experimentally observed distributions (for instance as seen by the authors in the work reported in [2225,28]), we model the synchrotron source as a Gaussian beam with constant phase at the source and amplitude given by: where r and θ are cylindrical coordinates and ω 0 is the 1/e half-waist at z = 0.In the following derivations we explicitly retain the physical parameters so that their effect may be readily discerned.Typically, for a synchrotron ω 0 is of the order of 10 microns in the vertical direction and order 100's of microns in the horizontal direction.The Fresnel condition is satisfied if: where r and r′ are the detector and source positions respectively.For the synchrotron experiment reported in [22] where r is of the order of millimeters, r′ of the order of 10's or 100's of microns, the wavelength was 0.13 nm and z was approximately 40 m, this condition is satisfied.The Fresnel diffraction integral of Eq. ( 1) evaluated at z = z 1 gives [39]: where: ρ and φ are the cylindrical coordinates in the z 1 plane, k = 2π/λ is the wavenumber; Our technique for generating a vortex beam is to insert a phase plate that, ideally, is nonabsorbing and imparts a helical phase to the beam.It should be noted that material properties at x-ray wavelengths are such that the real part of the refractive index, 1 -δ, is less than unity, with δ < 10 -3 .Consequently, for a phase plate that consists of a spiral ramp cut into a substrate, the phase of the beam will be advanced through the thick parts of the ramp compared to the thin parts.This is opposite to the case in visible wavelength optics and results in a vortex with the opposite helicity to that which might be expected.After insertion of the phase plate so that the centre of the plate is coaxial with the beam the amplitude is: where m is an integer referred to as the charge of the vortex: We can then calculate the Fresnel propagation of this vortex beam by evaluating the Fresnel integral at some detector plane Z = Z′ -z 1 , where Z′ is the distance of the detector plane from the synchrotron source: where R and Θ are the cylindrical coordinates in the Z plane.Recalling the definition of a Bessel function of the first kind of order m: and inserting Eq. (4a) into Eq.( 5) we can write: This can be evaluated [40] to give: where: I n is a Modified Bessel function of the first kind of order n; Typically, at visible wavelengths, Eq. ( 1) is used as the input to Eq. ( 5).However, as discussed above, in the synchrotron context we use Eq. ( 3).The retention of the phase terms in the input beam before implanting the vortex via a phase plate greatly complicates Eq. ( 8) compared to the corresponding result for visible wavelength.
The detected intensity is given by: A plot of Eq. ( 9) is shown in Fig. 1 for the synchrotron parameters described above and with Z = 5.8 m, thus satisfying Eq. ( 2) above.We assume ω 0 is the same both vertically and horizontally.We also calculate the detector plane intensity using a standard numerical propagator to calculate the detector plane complex amplitude: Here F and F -1 are the Fourier transform of the argument, k x and k y are the Fourier transform spatial frequencies, a is the complex amplitude at coordinates x,y,z, and a(x,y,0) is given by Eq. (3a).The resultant intensity plotted as for Fig. 1 is indistinguishable from Fig. 1.

Vortex 'interferograms' produced by division of wavefront
In practice the existence of a vortex is not demonstrated merely by producing the type of intensity distribution shown in Fig. 1.Ideally some demonstration of the phase distribution should be made.This can be done interferometrically.In the visible regime division of amplitude interferometric techniques are commonly used, however at x-ray wavelengths this type of experiment is not straightforward.Accordingly, in the x-ray context, we use a division of wavefront technique such as the introduction of a wire to act as a secondary source.The wire is placed close to the phase plate and the Fresnel integral can be calculated for the resultant field.The result in [36] for diffraction of a vortex beam through an aperture with appropriate modification for the limits of the integral could be used to derive a complicated expression for the diffracted field.However, the parameters imposed by the typical synchrotron experiment allow us to analyse the experiment as if it were the interference of the vortex beam in Eq. ( 8) and a cylindrical wave generated at the wire which has a phase gradient along it due to the incident vortex phase.Strictly, the wave generated by the wire will not be a cylindrical wave due to the curvature in the beam and the vortex phase imparted by the phase plate, however, the approximation is shown to be valid under numerical simulation for the synchrotron parameters used above.The resulting intensity is: where A v is the real amplitude of the vortex beam given by Eq. ( 8), A cyl is the real amplitude of the cylindrical wave and there is a modulation to the intensity distribution given by the cosine of the interference term: where x offs and z offs are the offset coordinates of the wire with respect to the phase plate and where the first atan expression evaluates the phase of its complex argument.If we approximate the vortex beam as a plane wave with a vortex phase then Eq. ( 12), after switching to Cartesian coordinates, becomes: Accordingly, we see that the analytic form for the interferogram for x-rays at a synchrotron becomes similar to the form used in calculating a hologram for visible wavelength vortices [3,5].
Figure 2(a) shows the modulation of the intensity distribution produced by the cosine of θ interf given by Eq. ( 12), while Fig. 2(b) shows the result for Eq. ( 13). Figure 2(c) shows the result for a numerical simulation based on Eq. (10). Figure 2d shows an experimental result as reported in [22].In all cases it can be seen that the phase discontinuity introduces a fork in one of the fringes.The main difference between the first two and the last two images is that the intensity distribution due to the amplitude terms in Eq. ( 11) are not included in the first two images and that the finite extent of the wire, absorption through the wire and the effects of partial coherence are included in the latter pair.

Vortex and an edge discontinuity: non-integer charge
Consider now what happens if the spiral phase plate is not matched to the energy of the incident x-ray beam so that the phase step across the step on the ramp is not an integer multiple of 2π.In general the 'charge' of the spiral will be some non-integer, ν.In the derivation above Eq.( 6) is modified by replacing m with ν and the Bessel function of the first kind of order m becomes a sum of an Anger and a Weber function of the first kind of order ν.Note that the function is not replaced by a Bessel function of order ν as this has a quite different integral representation [40].We could proceed via the modified form of Eq. ( 7), however the success of Eq. ( 13) in predicting the qualitative form of the downstream 'interferogram' suggests that replacing m with ν in that equation will suffice.The approach here becomes similar to that taken elsewhere in calculating the hologram required to generate a vortex [3,5,37].
In the synchrotron context it is relatively straightforward to continually vary the incident energy over a broad range.The charge of the phaseplate will vary as a function of the refractive index decrement times the energy.Away from absorption edges the refractive index decrement is roughly inversely proportional to the square of the energy, so the charge will be roughly inversely proportional with energy in this regime.Figure 4 is a link to a movie that tracks the evolution in the 'interferogram' produced by Eq. ( 14), with α = 0, as a function of equally spaced energy steps.

Conclusions
We have developed a useful range of results dealing with vortex phenomena at x-ray energies and for experimental parameters appropriate to a synchrotron source.Analytic expressions that differ from those obtained with visible wavelength assumptions are obtained and care has been taken to retain physical parameters explicitly.It is shown that the analytic forms can usefully predict intensity and phase distributions obtained numerically and experimentally.Finally, we have generated a prediction of vortex phase structure as a function of energy for a given phaseplate.We will compare this with experimental results in future work.