Partially coherent beam generation with metasurfaces

Roman Calpe; Atri Halder; Meilan Luo; Matias Koivurova; Jari Turunen

doi:doi:10.1364/PRJ.492233

Photonics Research, 2023, 11 (9): 1535, Published Online: Aug. 28, 2023

Partially coherent beam generation with metasurfaces

论文大纲

Roman Calpe ^1,*Atri Halder ¹Meilan Luo ^1,2Matias Koivurova ^3,4Jari Turunen ¹

Author Affiliations

¹ Center for Photonics Sciences, University of Eastern Finlandhttps://ror.org/00cyydd11, FI-80101 Joensuu, Finland

² Department of Physics and Synergetic Innovation Center for Quantum Effects and Applications, Hunan Normal University, Changsha 410081, China

³ Tampere Institute for Advanced Study, Tampere University, 33100 Tampere, Finland

⁴ Faculty of Engineering and Natural Sciences, Tampere University, FI-33720 Tampere, Finland

Abstract

An optical system for the generation of partially coherent beams with genuine cross-spectral density functions from spatially modulated globally incoherent sources is presented. The spatial intensity modulation of the incoherent source is achieved by quasi-planar metasurfaces based on spatial-frequency modulation of binary Bragg surface-relief diffraction gratings. Two types of beams are demonstrated experimentally: (i) azimuthally periodic, radially quasi-periodic beams and (ii) rotationally symmetric Bessel-correlated beams with annular far-zone radiation patterns.

1. INTRODUCTION

We explore implementations of an experimental scheme, proposed by Gori and Santarsiero [1,2], for the generation of a wide range of spatially partially coherent optical beams. The scheme is based on a globally incoherent field with a specific spatially varying intensity distribution $p (v)$ at the input plane of an optical system that performs a linear transformation described by an arbitrary spatial kernel function $H (ρ, v)$ . Partially coherent fields generated in this way are guaranteed to have genuine (non-negative definite) cross-spectral density functions (CSDs) $W (ρ_{1}, ρ_{2})$ .

In our implementation the required spatial variation of the intensity distribution is created by first reducing the spatial coherence of a laser beam with a random rotating diffuser [3 –5] and then imaging the scattered field onto an intensity-modulating metasurface. The kernel function of the optical system is taken to be of Fourier form, which is sufficient here since the fields to be demonstrated follow the Schell-model [6]. The generation scheme and its practical limitations are discussed in Section 2.

Among numerous possible options, we concentrate on a recently introduced class of azimuthally periodic and radially quasi-periodic partially coherent beams [7]. In a special case these beams have annular far-zone distributions, and in a certain limit they reduce to fundamental Laguerre–Gaussian correlated beams [8,9]. To make the present paper self-contained, we cover the main properties of these fields in Section 3.

The metasurfaces are realized as binary diffractive structures based on Bragg diffraction by surface-relief gratings as information carriers [10]. The spatial intensity information for coherence control is included by spatial-frequency modulation with a continuous variation of the local fill factor of the carrier grating [11,12]. The operating principle, design, and fabrication of the metasurfaces are the subject of Section 4.

The experimental results are presented in Section 5, along with simulations based on the theory in Section 3. Results of beam propagation over several distances $Δ z$ in the half-space $z > z_{0}$ are presented, and the complex degrees of coherence at the plane $z = 0$ , which give rise to the propagation properties, are measured with a mirror-based wavefront folding interferometer (WFI) [13].

While only particular examples of partially coherent beams are demonstrated here, the generation scheme is applicable to any beam covered by the model in Refs. [1,2]. It can also be extended to vector beams. These issues are discussed in Section 6 before conclusions.

2. DEVISING GENUINE PARTIALLY COHERENT BEAMS

Gori and Santarsiero [1] introduced a class of genuine spatially partially coherent fields, for which the CSD has a non-negative definite form $W (ρ_{1}, ρ_{2}) = \int_{- \infty}^{\infty} p (v) H^{*} (ρ_{1}, v) H (ρ_{2}, v) d^{2} v .$ (1)Here, $ρ_{1}$ and $ρ_{2}$ are transverse position coordinates in a plane $z = constant$ , $p (v)$ is a non-negative function of a 2D variable $v$ , and $H (ρ, v)$ is a kernel of an arbitrary linear transformation. The frequency dependence of the CSD is implicit. A wide range of physically relevant CSDs is obtained if the kernel function is taken to be of the form $H (ρ, v) = H (ρ) \exp (- i 2 π ρ \cdot v),$ (2)which can be realized with a 2D Fourier-transforming optical system. Upon inserting Eq. (2) into Eq. (1) we get $W (ρ_{1}, ρ_{2}) = H^{*} (ρ_{1}) H (ρ_{2}) \int_{- \infty}^{\infty} p (v) \exp (- i 2 π Δ ρ \cdot v) d^{2} v,$ (3)where $Δ ρ = ρ_{2} - ρ_{1}$ . This implies that the complex degree of spatial coherence depends only on the coordinate difference, resulting in fields which are of the Schell-model form. Fields that depart from this model require more sophisticated kernels, as discussed in Section 6.

The schematic setup for generating spatially partially coherent beams with genuine CSDs is illustrated in Fig. 1(a). One starts from a globally incoherent field at the input plane $A$ of the linear system $T$ , where the intensity distribution is proportional to the kernel function $H (ρ, v)$ . The weight function $p (v)$ is embedded in the metasurface at $A$ to modulate the intensity profile; then the CSD given by Eq. (1) is produced at the output plane $B$ . If desired, the output field can be modulated by adding a complex-amplitude transmission filter with transmission function $t (ρ)$ at plane $B$ , which transforms the CSD as $W^{'} (ρ_{1}, ρ_{2}) = t^{*} (ρ_{1}) t (ρ_{2}) W (ρ_{1}, ρ_{2})$ (4)and therefore adds freedom to the choice of realizable fields. In the simplest case the filter is just a hard or smoothly apodized aperture.

(a) Experimental setup for transforming a globally incoherent field at plane A into a custom-designed partially coherent beam at plane B using a linear optical transformation system T. (b) Generation of the globally incoherent field at plane A: D is a dynamic diffuser to reduce spatial coherence of the incident focused laser beam, L1 is a collimating lens, L2 and L3 form an afocal imaging system of magnification m′=1/5, G is a Gaussian filter, and the metasurface is placed at plane A. The beam radii at planes D, G, and A are wD, wG, and wA, and the coherence radii at G and B are σG and σB.

Fig. 1. (a) Experimental setup for transforming a globally incoherent field at plane $A$ into a custom-designed partially coherent beam at plane $B$ using a linear optical transformation system $T$ . (b) Generation of the globally incoherent field at plane $A$ : $D$ is a dynamic diffuser to reduce spatial coherence of the incident focused laser beam, $L_{1}$ is a collimating lens, $L_{2}$ and $L_{3}$ form an afocal imaging system of magnification $m^{'} = 1 / 5$ , $G$ is a Gaussian filter, and the metasurface is placed at plane $A$ . The beam radii at planes $D$ , $G$ , and $A$ are $w_{D}$ , $w_{G}$ , and $w_{A}$ , and the coherence radii at G and B are $σ_{G}$ and $σ_{B}$ .

下载图片查看所有图片

The system we use to generate the incoherent field at plane $A$ is illustrated schematically in Fig. 1(b). In our implementation a linearly polarized HeNe laser beam (wavelength $λ_{0} = 632.8 nm$ ) is focused to a width $w_{D}$ on a standard rotating diffuser to reduce its spatial coherence area. The scattered field is collimated by lens $L_{1}$ and passes through a Gaussian apodization filter $G$ to produce a Gaussian Schell-model (GSM) beam [5]. The transmitted beam is then demagnified by an afocal system consisting of lenses $L_{2}$ and $L_{3}$ to illuminate the metasurface with a GSM beam of radius $w_{A}$ .

Admittedly, the use of a rotating diffuser is a bulky and somewhat cumbersome approach to reduce spatial coherence. For future needs we therefore put forward an alternative technical solution: moving the diffuser randomly in the $x y$ plane with a piezoelectric device. Random movement over a sufficient area is possible since piezoelectric devices allow fast 2D scanning over distances larger than the required coherence area. Strictly speaking, the field on $A$ is not completely incoherent, but it is quasi-homogeneous. Therefore the conditions for the validity of Eq. (1) do not hold exactly, but non-negative definiteness is nevertheless ensured since the field is physically realizable.

3. PARTIALLY COHERENT MODEL BEAMS

We demonstrate the scheme in Fig. 1 by applying it to a recently introduced class of radially quasi-periodic and azimuthally periodic (RAP) beams [7]. These are generally Bessel-correlated beams with a Gaussian intensity profile at the source plane $z = z_{0}$ , which, in a special case, can have an annular distribution of radiant intensity. The main characteristics of RAP beams are briefly covered to make the paper self-contained and to summarize the formulas that enable numerical simulations of experimental results.

3.2 A. Radially Quasi-Periodic, Azimuthally Periodic Beams

We introduce RAP beams by assuming that the function $H (ρ)$ in Eq. (2) has a Gaussian form $H (ρ) = \exp (- ρ^{2} / w_{0}^{2})$ . The weight function $p (v)$ is written in polar coordinates $v = (v, ϕ)$ as $p (v, ϕ) = p_{m} \exp (- 2 π^{2} σ_{0}^{2} v^{2}) [1 + C_{m} J_{m} (2 π ρ_{0} v) \cos (m ϕ)]$ (5)with $p_{m} = 2 π σ_{0}^{2} {\begin{cases} {1 + C_{m} \exp [- ρ_{0}^{2} / (2 σ_{0}^{2})]}^{- 1} & when m = 0, \\ 1 & when m \neq 0 . \end{cases}$ (6)Here, $σ_{0}$ and $ρ_{0}$ are positive constants, $J_{m} (z)$ denotes the Bessel function of the first kind and order $m$ , $C_{0} = \pm 1$ , and $C_{m} = 1$ for $m > 1$ . Equations (1) and (2) result in CSDs of the form $W (ρ_{1}, ρ_{2}) = \sqrt{S (ρ_{1}) S (ρ_{2})} μ (Δ ρ)$ , where the spectral density at $z = z_{0}$ , using Eqs. (5) and (6), has a Gaussian profile $S (ρ) = \exp (- 2 ρ^{2} / w_{0}^{2}) .$ (7)

The complex degree of coherence (DOC) for solutions $m \geq 1$ is $μ (Δ ρ) = \exp (- Δ ρ^{2} / 2 σ_{0}^{2}) \times [1 + C_{m} {(- i)}^{m} \cos (m Δ φ) \exp (- \frac{ρ_{s}^{2}}{2}) I_{m} (\frac{ρ_{s} Δ ρ}{σ_{0}})],$ (8)where $Δ ρ = (Δ ρ, Δ φ)$ , $ρ_{s} = ρ_{0} / σ_{0}$ , and $I_{m} (z)$ denotes the modified Bessel function of the first kind and order $m$ .

A closed-form expression for the propagated CSD in the paraxial domain can be found with the standard Fresnel formula [7]. In particular, the spectral density, which we will measure in Section 5, evolves as $S (ρ, z) = \frac{w_{0}^{2}}{w^{2} (Δ z)} \exp [- \frac{2 ρ^{2}}{w^{2} (Δ z)}] {1 + C_{m} {(- 1)}^{m} \cos (m φ) \times \exp [- \frac{R^{2}}{2 q^{2}} c (Δ z)] J_{m} [\frac{2 R}{q \sqrt{1 + q^{2}}} \frac{Δ z / z_{R}}{1 + {(Δ z / z_{R})}^{2}} \frac{ρ}{w_{0}}]},$ (9)where $ρ = (ρ, φ)$ , $Δ z = z - z_{0}$ is the propagation distance from plane $B$ , $q = σ_{0} / w_{0}$ , and $R = ρ_{0} / w_{0}$ . Further, $w (Δ z) = w_{0} {[1 + {(Δ z / z_{R})}^{2}]}^{1 / 2}$ (10)is the beam width at $z > 0$ , and the quantity $c (Δ z)$ is $c (Δ z) = \frac{1 + {(Δ z / z_{G})}^{2}}{1 + {(Δ z / z_{R})}^{2}},$ (11)where $z_{G} = k w_{0}^{2} / 2$ is the Rayleigh range of a fully coherent Gaussian beam, and $z_{R} = z_{G} {(1 + q^{- 2})}^{- 1 / 2}$ (12)is the Rayleigh range of a partially coherent GSM beam [14] characterized by coherence parameter $q$ .

3.3 B. Beams with Annular Far-Zone Radiation Pattern

Let us consider RAP beams with $m = 0$ and $C_{0} = \pm 1$ in Eq. (5), which gives a weight function $p (v, ϕ) = p_{0} \exp (- 2 π^{2} σ_{0}^{2} v^{2}) [1 \pm J_{0} (2 π ρ_{0} v)],$ (13)where $p_{0} = 2 π σ_{0}^{2} B_{0}$ and $B_{0} = {[1 \pm \exp (- ρ_{s}^{2} / 2)]}^{- 1} .$ (14)

The spectral density at the source plane is still given by Eq. (7), but the complex DOC takes the form $μ (Δ ρ) = B_{0} \exp (- Δ ρ^{2} / 2 σ_{0}^{2}) \times [1 \pm \exp (- \frac{ρ_{s}^{2}}{2}) I_{0} (\frac{ρ_{s} Δ ρ}{σ_{0}})] .$ (15)

The spectral density at $z > 0$ , $S (ρ, z) = B_{0} \frac{w_{0}^{2}}{w^{2} (Δ z)} \exp [- \frac{2 ρ^{2}}{w^{2} (Δ z)}] \times {1 \pm \exp [- \frac{R^{2}}{2 q^{2}} c (Δ z)] J_{0} [\frac{2 R}{q \sqrt{1 + q^{2}}} \frac{Δ z / z_{R}}{1 + {(Δ z / z_{R})}^{2}} \frac{ρ}{w_{0}}]},$ (16)takes an annular (donut-shaped) form upon propagation to the far field. It therefore shares some properties of other ‘hollow’ beams, including fully coherent higher-order Bessel beams [15,16], orbital angular momentum, azimuthally polarized beams [17,18], generalized Bessel–Gauss beams [19], partially coherent dark propagation-invariant beams [20,21], and various Bessel-correlated beams that may have donut shapes either initially or in the far-zone [22]. Donut-shaped beams have proven important in particle trapping [23] and stimulated emission depletion (STED) microscopy [24].

In the limiting case $ρ_{s} \to 0$ , and when $C_{0} = - 1$ , the complex DOC given in Eq. (15) reduces to $μ (Δ ρ) = \exp (- \frac{Δ ρ^{2}}{2 σ_{0}^{2}}) (1 - \frac{Δ ρ^{2}}{2 σ_{0}^{2}}),$ (17)which is the DOC of a Laguerre–Gauss Schell-model (LGSM) source of order $n = 1$ [8]. The spectral density of the LGSM beams takes the form $S (ρ, z) = \frac{w_{0}^{2}}{w^{2} (Δ z)} \exp [- \frac{2 ρ^{2}}{w^{2} (Δ z)}] \times [c (Δ z) + \frac{2 Δ z^{2} z_{R}^{2} ρ^{2}}{(w_{0}^{2} + σ_{0}^{2}) {(Δ z^{2} + z_{R}^{2})}^{2}}] .$ (18)

If $C_{0} = 1$ and we again let $ρ_{s} \to 0$ , Eq. (15) becomes the DOC of a GSM beam. Figure 2 illustrates some properties of the $m = 0$ solutions not discussed in Ref. [7]. Cross sections of the DOC plotted in Fig. 2(a) show central peaks and side lobes of expanding widths when the value of $ρ_{s}$ in the $I_{0}$ term in Eq. (15) increases. The far-field intensity in Fig. 2(b) is also affected: the central minimum is close to zero at small values of $ρ_{s}$ , which also provides the widest central dark region.

Properties of RAP beams of order m=0 with λ=633 nm, w0=300 μm, and q=1/30. (a) Cross sections of the DOC and (b) far-field (at z=10zR) spectral densities when C0=−1, ρs→0 (black), ρs=1 (red), ρs=2 (green), and ρs=4 (blue).

Fig. 2. Properties of RAP beams of order $m = 0$ with $λ = 633 nm$ , $w_{0} = 300 μm$ , and $q = 1 / 30$ . (a) Cross sections of the DOC and (b) far-field (at $z = 10 z_{R}$ ) spectral densities when $C_{0} = - 1$ , $ρ_{s} \to 0$ (black), $ρ_{s} = 1$ (red), $ρ_{s} = 2$ (green), and $ρ_{s} = 4$ (blue).

下载图片查看所有图片

4. METASURFACE DESIGN AND REALIZATION

In Ref. [7] we presented an on-axis encoding technique of the metasurface at plane $A$ , which can be realized in the axially symmetric geometry of Fig. 1. However, our experiments showed that this scheme has axial noise, which proved difficult to avoid. We therefore opted for an alternative, off-axis coding scheme. Our method is effectively similar to that used in Ref. [25], with the difference that we modulate the field amplitude with a more compact metasurface, instead of filtering the diffracted light from a spatial light modulator. In addition, we can obtain extremely high diffraction efficiency with the off-axis encoding. The drawback of our technique is the somewhat trickier alignment of the optical components compared to the on-axis method.

4.2 A. Operating Principle and Carrier Design

The metasurface is realized using a high-frequency carrier grating as introduced in Ref. [10]. This is a binary surface-relief grating with depth $h$ and fill factor $f = c / d$ as design parameters, as defined in Fig. 3(a). The grating operates at large angles of incidence $θ \sim 30 °$ . Hence the geometry of Fig. 1 is folded at $A$ by twice this angle. The design angle of incidence is the Bragg angle $θ = θ_{B}$ given by $\sin θ_{B} = \frac{λ_{0}}{2 nd} .$ (19)

Assuming that the incident field is linearly polarized in the direction of the grating grooves, we design the carrier grating for maximum efficiency of order $l = - 1$ in TE polarization for $λ_{0} = 633 nm$ , $n = 1.49$ , and $d = λ_{0}$ .

$(a) Operating geometry and design parameters of the Bragg carrier grating. θ is the first Bragg angle of incidence from air to the substrate of refractive index ns and θ′ from the substrate to the grating layer with refractive index n. h is the groove depth, c is the ridge width, and d is the period. (b) Efficiency η−1 of order l=−1 as a function of the fill factor in TE polarization. (c) Efficiency η−1 as a function of angle of incidence for different fill factors.$

Fig. 3. (a) Operating geometry and design parameters of the Bragg carrier grating. $θ$ is the first Bragg angle of incidence from air to the substrate of refractive index $n_{s}$ and $θ^{'}$ from the substrate to the grating layer with refractive index $n$ . $h$ is the groove depth, $c$ is the ridge width, and $d$ is the period. (b) Efficiency $η_{- 1}$ of order $l = - 1$ as a function of the fill factor in TE polarization. (c) Efficiency $η_{- 1}$ as a function of angle of incidence for different fill factors.

下载图片查看所有图片

The efficiency of the carrier grating as a function of fill factor $f = c / d$ was calculated using the Fourier modal method (FMM) [26]. The results are shown in Fig. 3(b). The efficiency has two maxima reaching almost unity, at $f \sim 0.16$ and $f \sim 0.62$ . We employ the second since the efficiency drops smoothly to zero when $f$ approaches unity. Since the incident field is not fully collimated, we plot in Fig. 3(c) the efficiency as a function of the angle of incidence $θ$ (Bragg selectivity of the grating) for several fill factors over the range $0 ° < θ < 40 °$ . It turns out that the angular selectivity is not a serious limiting factor for operation of the grating.

The intensity coding is performed as local variations of the fill factor $f$ as in Refs. [10,11]. Fill factors between 0.62 and 0.9, marked with dashed vertical lines in Fig. 3(b), are used for efficiency encoding since they are practically feasible. The target grating depth was chosen to be $h = 1350 nm$ , giving $η_{- 1} \approx 0.95$ at the second maximum $f = 0.62$ .

4.3 B. Fabrication

The sample was fabricated by a standard electron beam lithography procedure. A 1-inch-diameter and 0.5-mm-thick ${SiO}_{2}$ substrate was cleaned and spin-coated with a 1.4-μm-thick positive-tone resist layer (PMMA), which was coated with a conductive 30-nm-thick copper layer via thermal evaporation. After creating the grating patterns with electron beam lithography, the copper layer was removed in ${HNO}_{3}$ , followed by resist development in MIBK:IPA (1:3 solution) for 75 s and 30 s rinsing with IPA.

Scanning electron microscope images of one of the samples are shown in the Fig. 4. In Fig. 4(a) we show a cross section around a location where the fill factor changes abruptly from 0.8 to 0.62, and in Fig. 4(b) a top view of the region is shown where $f$ is fairly constant. The fill factor is fairly well reproduced, but in view of Fig. 4(a) the groove depth is somewhat reduced especially at large $f$ . This could be improved by increasing the electron beam dose. However, its main effect is to reduce $η_{- 1}$ , which should be low at large $f$ according to Fig. 3. In fact the modulation of $h$ helps us to achieve the low values of $η_{- 1}$ required beyond the second vertical line in Fig. 3(b).

Fig. 4. SEM images of a fabricated grating. (a) Side view. (b) Top view.

下载图片查看所有图片

5. EXPERIMENTAL RESULTS

The focal lengths of the lenses in Fig. 1(b) before plane A are flexible, but we chose $L_{1} = 150 mm$ , $L_{2} = 200 mm$ , and $L_{3} = 40 mm$ . We measured $w_{D} = 31 μm$ and $w_{A} \approx 880 μm$ . The GSM was used as in Ref. [5] to estimate the coherence radius at the filter $G$ (Thorlabs NDYR20A) to be $σ_{G} \approx 210 μm$ , which means that the field is quasi-homogenous with $σ_{1} / w_{G} \approx 1 / 21$ . Measurement with a WFI gives the results in Fig. 5, where we show cross sections of the spectral density and the absolute value of the spatial DOC at plane $A$ , from which $w_{A} \approx 880 μm$ . While $| μ (Δ x) |$ is not completely Gaussian but has some side lobes, we can estimate $σ_{A} \approx 42 μm$ . The field at plane $A$ is definitely quasi-homogeneous and, since the coherence area is larger than $d$ by a factor of 60, the metasurface works well as a Bragg deflector.

Fig. 5. Measured absolute value of the spatial DOC $μ (Δ x, z_{0})$ (blue) and the spectral density $S (x, z_{0})$ (orange) at plane $A$ .

下载图片查看所有图片

The Gaussian spot at plane $A$ , with $w_{A} \approx 880 μm$ , could be used directly to approximate the Gaussian factor in $p (v)$ . Hence, to encode $p (v)$ , only the residual RAP factor needed to be realized with the metasurface. Explicitly, this factor is $M (v, ϕ) = M_{0} [1 + J_{m} (2 π ρ_{0} v) \cos (m ϕ)],$ (20)for beams with $m > 0$ and $M (v, ϕ) = M_{0} [1 - J_{0} (2 π ρ_{0} v)],$ (21)for $m = 0$ . Here, $M_{0}$ is chosen such that the maximum does not exceed unity. Sampled $100 \times 100$ versions of the metasurfaces with spatially varying transmission efficiencies are illustrated in Fig. 6. In particular, we show the $m = 3$ , $ρ_{0} = 3 mm$ solution in Eq. (20) and the $m = 0$ , $ρ_{0} = 2 mm$ solution in Eq. (21), which are demonstrated experimentally. The modulation contrast of $p (v)$ is much weaker in the former case than in the latter. The linear system $T$ , which does the Fourier transform from plane $A$ to plane $B$ , is effectively a positive lens of focal length $f_{H} = 50 mm$ . Therefore, the spatial frequency $v$ at plane $A$ is connected to the real-space coordinate $ρ$ via $v = ρ / λ_{0} f_{H}$ .

Fig. 6. Spatially varying transmission efficiency $η_{- 1}$ of the designed metasurface for the RAP beams of order (a) $m = 3$ and (b) $m = 0$ .

下载图片查看所有图片

5.5 A. Source-Plane Complex Degree of Coherence

The measured amplitude and phase of the DOC at plane $z = z_{0}$ are shown in Fig. 7 for the two beams considered. The distributions of $μ (Δ x, Δ y)$ are evaluated using a mirror-based WFI as in Ref. [27]. We show the measured phases only across a circular area of radius $Δ ρ < 75 μm$ for $m = 3$ and $Δ ρ < 42 μm$ for $m = 0$ , where the absolute value of the DOC is significant and reliable phase measurements could be made.

Measured and simulated complex degrees of source-plane spatial coherence of the RAP beam with m=3 (top row) and m=0 (bottom row). (a)–(d) Absolute values and (e)–(h) phases of μ(Δx,Δy). Here, (a), (b), (e), and (f) are the measured results, while (c), (d), (g), and (h) are the simulated results.

Fig. 7. Measured and simulated complex degrees of source-plane spatial coherence of the RAP beam with $m = 3$ (top row) and $m = 0$ (bottom row). (a)–(d) Absolute values and (e)–(h) phases of $μ (Δ x, Δ y)$ . Here, (a), (b), (e), and (f) are the measured results, while (c), (d), (g), and (h) are the simulated results.

下载图片查看所有图片

The DOC of the $m = 3$ solution has an azimuthal periodicity of 60° in the side lobe pattern, whereas the phase in the central region $Δ ρ < 35 μm$ is constant. These features, as well as the sharp variations of $| μ (Δ x, Δ y) |$ , are evident in the measurements even though there are some distortions due to experimental imperfections particularly in the phase measurements within the central region of Fig. 7(e). The abrupt $2 π$ phase transitions in the central region of Fig. 7(f) are due to wrapping into the $[- π, π)$ interval.

5.6 B. Propagation of Beam Profiles

The propagation of beams that the DOCs in Fig. 7 generate is illustrated in Figs. 8 and 9, where the experimental and simulated results are shown in top and bottom rows, respectively. The experimental results show some speckle-like granularity that could be contributed by the finite coherence area at plane $A$ , where the field is not strictly speaking globally incoherent as assumed in the model. In particular, the discrepancies caused by the granularity and additional side lobes in the coherence function of the $m = 3$ RAP beam in the Fig. 7 upper panel are visible in the measured intensity profiles in Fig. 8 at $z = 90 mm$ and $z = 140 mm$ .

Illustration of propagation characteristics of RAP-beam with m=3. The measured intensity profiles at four different propagation distances are in the top row, and the simulation is in the bottom row.

Fig. 8. Illustration of propagation characteristics of RAP-beam with $m = 3$ . The measured intensity profiles at four different propagation distances are in the top row, and the simulation is in the bottom row.

下载图片查看所有图片

Fig. 9. Same as Fig. 8 but for the beam with $m = 0$ . The parameters used for the calculations were $w_{0} = 210 μm$ , $R = 1 / 7.5$ , and $q = 1 / 30$ .

下载图片查看所有图片

The evolution of the RAP beam with $m = 3$ is shown in Fig. 8 at the source plane and three different propagation distances. We assume parameters $R = 1 / 4.5$ and $q = 1 / 40$ , which gives $z_{G} = 567 mm$ and $z_{R} = 14.2 mm$ . At $z = 0$ the beam shape is Gaussian with a beam waist $w_{0} = 338 μm$ . As the propagation distance grows the beam acquires its RAP features. The annular $m = 0$ beam in Fig. 9, with $w_{0} = 210 μm$ , $R = 1 / 7.5$ , $q = 1 / 30$ , $z_{G} = 218 mm$ , and $z_{R} = 7.3 mm$ , begins to show its dark center at distances $z \sim 10 mm$ . This property is preserved throughout further propagation.

6. DISCUSSION

The RAP beams demonstrated here are not the only possible classes of partially coherent fields that we can generate with the experimental technique discussed in this paper. While Fourier-type kernels allow the demonstration of any Schell-model beam, the use of more general kernels opens the way to produce more general fields. For instance, metasurfaces acting as vortex phase filters with transmittance $t (ρ) = t (ρ, ϕ) = \exp (i l ϕ)$ , where $l$ is an integer, can be placed at plane $B$ to produce orbital angular momentum. Likewise, polarizing optical components can contribute to the linear transformations, allowing the realization of vector beams. Space-variant imaging systems with beam shaping or polarization-control metasurfaces (such as polarization gratings [28,29]) in the intermediate Fourier plane of a $4 F$ imaging system are options, as are transformations involving interferometers. For example, a WFI transforms any Schell-model beam into a corresponding specular beam [30,31]. We have recently demonstrated such specular transformations using a mirror-based WFI [27].

7. CONCLUSIONS

We have presented an implementation of an experimental system for the generation of spatially partially coherent beams with genuine correlation functions, using quasi-planar metasurfaces. The method was applied to particular radially quasi-periodic, azimuthally periodic beams. However, it is applicable to the generation of any physically realizable beams, whether scalar or vectorial.

References

[1] F. Gori, M. Santarsiero. Devicing genuine spatial correlation functions. Opt. Lett., 2007, 32: 3531-3533.

[2] F. Gori, M. Santarsiero. Devicing genuine twisted cross-spectral densities. Opt. Lett., 2018, 43: 595-598.

[3] W. Martiensen, E. Spiller. Coherence and fluctuations in light beams. Am. J. Phys., 1964, 32: 919-926.

[4] J. D. Farina, L. M. Narducci, E. Collett. Generation of highly directional beams from a globally incoherent source. Opt. Commun., 1980, 32: 203-208.

[5] Q. He, J. Turunen, A. T. Friberg. Propagation and imaging experiments with Gaussian Schell-model beams. Opt. Commun., 1988, 67: 245-250.

[6] A. C. Schell. A technique for the determination of the radiation pattern from a partially coherent aperture. IEEE Trans. Ant. Propag., 1967, 15: 187-188.

[7] M. Luo, M. Koivurova, J. Turunen. Azimuthally periodic and radially quasi-periodic Bessel-correlated fields. Opt. Express, 2022, 30: 11754-11766.

[8] Z. Mei, O. Korotkova. Random sources generating ring-shaped beams. Opt. Lett., 2013, 38: 91-93.

[9] M. Luo, D. Zhao. Elliptical Laguerre Gaussian Schell-model beams with a twist in random media. Opt. Express, 2019, 27: 30044-30054.

[10] J. Turunen, P. Blair, J. M. Miller, M. R. Taghizadeh, E. Noponen. Bragg holograms with binary synthetic surface-relief profile. Opt. Lett., 1993, 18: 1022-1024.

[11] E. Noponen, J. Turunen. Binary high-frequency-carrier diffractive optical elements: electromagnetic theory. J. Opt. Soc. Am. A, 1994, 11: 1097-1109.

[12] E. Tervonen, J. Turunen, J. Pekola. Pulse-frequency-modulated high-frequency-carrier diffractive elements for pattern projection. Opt. Eng., 1994, 33: 2579-2587.

[13] A. Halder, H. Partanen, A. Leinonen, M. Koivurova, T. K. Hakala, T. Setälä, J. Turunen, A. T. Friberg. Mirror-based scanning wavefront-folding interferometer for coherence measurements. Opt. Lett., 2020, 45: 4260-4263.

[14] A. T. Friberg, R. Sudol. Propagation parameters of Gaussian Schell model beams. Opt. Commun., 1982, 41: 383-387.

[15] J. Durnin. Exact solutions for nondiffracting beams: I. The scalar theory. J. Opt. Soc. Am. A, 1987, 4: 651-654.

[16] J. Durnin, J. J. Miceli, J. H. Eberly. Diffraction-free beams. Phys. Rev. Lett., 1987, 58: 1499-1501.

[17] R. H. Jordan, D. H. Hall. Free-space azimuthal paraxial wave equation: the azimuthal Bessel–Gauss beam solution. Opt. Lett., 1994, 19: 427-429.

[18] D. H. Hall. Vector–beam solutions of Maxwell’s wave equation. Opt. Lett., 1996, 21: 9-11.

[19] V. Bagini, F. Frezza, M. Santarsiero, G. Schettini, G. S. Spagnolo. Generalized Bessel–Gauss beams. J. Mod. Opt., 1996, 43: 1155-1166.

[20] J. Turunen, A. Vasara, A. T. Friberg. Propagation-invariance and self-imaging in variable-coherence optics. J. Opt. Soc. Am. A, 1991, 8: 282-289.

[21] S. A. Ponomarenko, W. Huang, M. Cada. Dark and antidark diffraction-free beams. Opt. Lett., 2007, 32: 2508-2510.

[22] A. Halder, M. Koivurova, H. Partanen, J. Turunen. Paraxial propagation of a class of Bessel-correlated fields. Opt. Express, 2018, 26: 11055-11067.

[23] M. Woerdemann, C. Alpmann, M. Esseling, C. Denz. Advanced optical trapping by complex beam shaping. Laser Photon. Rev., 2013, 7: 839-854.

[24] S. W. Hell, J. Wichmann. Breaking the diffraction resolution limit by stimulated emission: stimulated-emission-depletion fluorescence microscopy. Opt. Lett., 1994, 19: 780-782.

[25] Y. Cai, Y. Chen, F. Wang. Generation and propagation of partially coherent beams with nonconventional correlation functions: a review [Invited]. J. Opt. Soc. Am. A, 2014, 31: 2083-2096.

[26] KimH.ParkJ.LeeB., Fourier Modal Method and Its Applications in Computational Nanophotonics (CRC Press, 2017).

[27] D. Das, A. Halder, H. Partanen, M. Koivurova, J. Turunen. Propagation of Bessel-correlated specular and antispecular beams. Opt. Express, 2022, 30: 5709-5721.

[28] F. Gori. Measuring Stokes parameters by means of a polarization grating. Opt. Lett., 1999, 24: 584-586.

[29] J. Tervo, J. Turunen. Paraxial-domain diffractive elements with 100% efficiency based on polarization gratings. Opt. Lett., 2000, 25: 785-786.

[30] F. Gori, G. Guattari, C. Palma, C. Padovani. Specular cross-spectral density functions. Opt. Commun., 1988, 68: 239-243.

[31] H. Partanen, N. Sharmin, J. Tervo, J. Turunen. Specular and antispecular light beams. Opt. Express, 2015, 23: 28718-28727.

1. INTRODUCTION

2. DEVISING GENUINE PARTIALLY COHERENT BEAMS

3. PARTIALLY COHERENT MODEL BEAMS

3.2 A. Radially Quasi-Periodic, Azimuthally Periodic Beams

3.3 B. Beams with Annular Far-Zone Radiation Pattern

4. METASURFACE DESIGN AND REALIZATION

4.2 A. Operating Principle and Carrier Design

4.3 B. Fabrication

5. EXPERIMENTAL RESULTS

5.5 A. Source-Plane Complex Degree of Coherence

5.6 B. Propagation of Beam Profiles

6. DISCUSSION

7. CONCLUSIONS

Roman Calpe, Atri Halder, Meilan Luo, Matias Koivurova, Jari Turunen. Partially coherent beam generation with metasurfaces[J]. Photonics Research, 2023, 11(9): 1535.

Partially coherent beam generation with metasurfaces

1. INTRODUCTION

2. DEVISING GENUINE PARTIALLY COHERENT BEAMS

3. PARTIALLY COHERENT MODEL BEAMS

3.2 A. Radially Quasi-Periodic, Azimuthally Periodic Beams

3.3 B. Beams with Annular Far-Zone Radiation Pattern

4. METASURFACE DESIGN AND REALIZATION

4.2 A. Operating Principle and Carrier Design

4.3 B. Fabrication

Fig. 4. SEM images of a fabricated grating. (a) Side view. (b) Top view.

5. EXPERIMENTAL RESULTS

Fig. 5. Measured absolute value of the spatial DOC $μ (Δ x, z_{0})$ (blue) and the spectral density $S (x, z_{0})$ (orange) at plane $A$ .

Fig. 6. Spatially varying transmission efficiency $η_{- 1}$ of the designed metasurface for the RAP beams of order (a) $m = 3$ and (b) $m = 0$ .

5.5 A. Source-Plane Complex Degree of Coherence

5.6 B. Propagation of Beam Profiles

Fig. 8. Illustration of propagation characteristics of RAP-beam with $m = 3$ . The measured intensity profiles at four different propagation distances are in the top row, and the simulation is in the bottom row.

Fig. 9. Same as Fig. 8 but for the beam with $m = 0$ . The parameters used for the calculations were $w_{0} = 210 μm$ , $R = 1 / 7.5$ , and $q = 1 / 30$ .

6. DISCUSSION

7. CONCLUSIONS

Article Outline

关于本站 Cookie 的使用提示

全站搜索

Partially coherent beam generation with metasurfaces

1. INTRODUCTION

2. DEVISING GENUINE PARTIALLY COHERENT BEAMS

3. PARTIALLY COHERENT MODEL BEAMS

3.2 A. Radially Quasi-Periodic, Azimuthally Periodic Beams

3.3 B. Beams with Annular Far-Zone Radiation Pattern

Fig. 2. Properties of RAP beams of order m=0 with λ=633 nm, w0=300 μm, and q=1/30. (a) Cross sections of the DOC and (b) far-field (at z=10zR) spectral densities when C0=−1, ρs→0 (black), ρs=1 (red), ρs=2 (green), and ρs=4 (blue).

4. METASURFACE DESIGN AND REALIZATION

4.2 A. Operating Principle and Carrier Design

4.3 B. Fabrication

Fig. 4. SEM images of a fabricated grating. (a) Side view. (b) Top view.

5. EXPERIMENTAL RESULTS

Fig. 5. Measured absolute value of the spatial DOC μ(Δx,z0) (blue) and the spectral density S(x,z0) (orange) at plane A.

Fig. 6. Spatially varying transmission efficiency η−1 of the designed metasurface for the RAP beams of order (a) m=3 and (b) m=0.

5.5 A. Source-Plane Complex Degree of Coherence

5.6 B. Propagation of Beam Profiles

Fig. 8. Illustration of propagation characteristics of RAP-beam with m=3. The measured intensity profiles at four different propagation distances are in the top row, and the simulation is in the bottom row.

Fig. 9. Same as Fig. 8 but for the beam with m=0. The parameters used for the calculations were w0=210 μm, R=1/7.5, and q=1/30.

6. DISCUSSION

7. CONCLUSIONS

Article Outline

相关论文

相关资讯

关于本站 Cookie 的使用提示

全站搜索

Fig. 5. Measured absolute value of the spatial DOC $μ (Δ x, z_{0})$ (blue) and the spectral density $S (x, z_{0})$ (orange) at plane $A$ .

Fig. 6. Spatially varying transmission efficiency $η_{- 1}$ of the designed metasurface for the RAP beams of order (a) $m = 3$ and (b) $m = 0$ .

Fig. 8. Illustration of propagation characteristics of RAP-beam with $m = 3$ . The measured intensity profiles at four different propagation distances are in the top row, and the simulation is in the bottom row.

Fig. 9. Same as Fig. 8 but for the beam with $m = 0$ . The parameters used for the calculations were $w_{0} = 210 μm$ , $R = 1 / 7.5$ , and $q = 1 / 30$ .