Next Article in Journal
Metasurface Holography in the Microwave Regime
Next Article in Special Issue
Linear Combinations of the Complex Degrees of Coherence
Previous Article in Journal
Nonreciprocal and Topological Plasmonics
Previous Article in Special Issue
Evolution of Spatiotemporal Intensity of Partially Coherent Pulsed Beams with Spatial Cosine-Gaussian and Temporal Laguerre–Gaussian Correlations in Still, Pure Water
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Type of Shape-Invariant Beams with Structured Coherence: Laguerre-Christoffel-Darboux Beams

by
Rosario Martínez-Herrero
1,†,
Massimo Santarsiero
2,*,†,
Gemma Piquero
1,† and
Juan Carlos González de Sande
3,†
1
Departamento de Óptica, Universidad Complutense de Madrid, Ciudad Universitaria, 28040 Madrid, Spain
2
Dipartimento di Ingegneria, Università Roma Tre, Via V. Volterra 62, 00146 Rome, Italy
3
ETSIS de Telecomunicación, Campus Sur, Universidad Politécnica de Madrid, 28031 Madrid, Spain
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Photonics 2021, 8(4), 134; https://doi.org/10.3390/photonics8040134
Submission received: 1 April 2021 / Revised: 17 April 2021 / Accepted: 19 April 2021 / Published: 20 April 2021
(This article belongs to the Special Issue Structured Light Coherence)

Abstract

:
A new class of sources presenting structured coherence properties is introduced and analyzed. They are obtained as the incoherent superposition of coherent Laguerre-Gaussian modes with suitable coefficients. This ensures that the shape of the intensity profile and the spatial coherence features of the propagated beams are invariant during paraxial approximation. A simple analytical expression is obtained for the cross-spectral density of the sources of this class, regardless of the number of superposed modes. Properties of these sources are analyzed and described by several examples.

1. Introduction

The spatial coherence properties of a scalar planar source up to the second order are taken into account by its cross-spectral density function, or CSD, which gives the correlation between the field values at two distinct points of the plane at any temporal frequency. On the CSD of a source depend, in particular, also the intensity and coherence properties of the field it radiates [1].
The research on possible forms of CSD that give rise to fields with peculiar propagation properties has been going on continuously since the first seminal papers of the late 1970s [2,3,4,5]. One of the difficulties one encounters in devising new forms of CSDs is due to the fact that the latter are not generic functions of two points in space. In fact, necessary and sufficient condition for them to represent a valid CSD is that they are kernels of nonnegative definite integral hermitian operators [1]. When this happens, the CSD is said to be genuine, or bona-fide. If such requirement is not met, the function cannot represent the CSD of a possible physical source. In general, it is not easy to check the nonnegativity of an integral kernel, even though genuineness criteria have been introduced to this aim [6,7].
A different way, but equivalent to the above one, to envisage new forms of CSD consists in starting from the modal expansion of the CSD. In fact, on exploiting a fundamental result of coherence theory, any bona-fide CSD can be written as the (possibly infinite) sum of terms, each of them corresponding to the CSD of a perfectly coherent light beam (a mode of the source) [1,7,8,9]. Each term is weighted by a non-negative coefficient (the eigenvalue). In physical terms, this means that any partially coherent source can be thought of as the superposition of an ensemble of mutually uncorrelated, perfectly coherent fields. The eigenvalues then play the role of the powers of the various modes in the superposition.
On summing the CSDs of a number of coherent fields, of course, the resulting CSD is genuine but, on the down side, the obtained result is rarely in closed form, that is, it is generally not given in terms of simple mathematical functions. On the contrary, this is a fundamental requirement in order for the source features, as well as the ones of the propagated fields, to be deduced from the analytical properties of the involved functions. Moreover, when the number of the superimposed modes is very large, the availability of a closed formula for the CSD turns out to be useful even from a computational point of view.
As an example, the modes of the so-called Gaussian Schell-model (GSM) sources [10,11] that will be recalled later on, are Hermite–Gaussian (in rectangular coordinates) or Laguerre–Gaussian (LG for short, in 2D polar coordinates) functions, i.e., the ones that represent the steady states of a quantum harmonic oscillator [12]. The corresponding eigenvalues decrease geometrically with the mode order. The final result of the sum of infinite terms is a CSD whose intensity and coherence profiles have simple Gaussian shapes.
In the present paper we introduce a new class of partially coherent sources that are obtained as the superposition of a finite number of LG modes. Whatever the number of modes involved, the CSDs of this class can be always be written in a simple closed form, whose expression depends on the number of the modes, but also on a second parameter, namely, the angular index (or topological charge) of the LG functions. In such a way, the structuring characteristics of the source (concerning both its intensity and degree of coherence) can be varied to a great extent on changing the above parameters. Due to the mathematical expression giving rise to the resulting closed form, sources of this class will be called Laguerre-Christoffel-Darboux (LCD) sources. Analogous sources have been recently introduced for the onedimensional case [13]. It will be seen that in some cases the intensity of the resulting source has a donut-like profile and it is practically uniform and vanishing in the inner and outer regions of the donut.
A significant feature of LCD sources, due to the propagation properties of the underlying modes, is that the fields radiated from them keep their shape unchanged during paraxial propagation at any distance from the source, except for a transverse scaling factor and a spherical curvature term. Such invariance concerns the whole CSD of the field (and, therefore, both the intensity profile and the coherence properties) and we will say that their CSDs are shape-invariant.
Another reason of interest for sources having modes of the LG type comes from the fact that such functions present wave-front singularities, known as optical vortices, and carry optical angular momentum (OAM). The presence of structures of this kind in coherent and partially coherent light beams is nowadays the subject of many research works (see for example [14,15,16,17,18,19,20,21] and references therein).
The paper is structured as follows. In Section 2 the concepts that are at the basis of the modal theory of coherence are recalled, as well as some known examples of partially coherent source obtained as superpositions of LG modes. In Section 3 the LCD sources are introduced, and their main features are analyzed in details. Finally, the main results of this work are summarized in Section 4.

2. Preliminaries

Partially coherent planar sources can be appropriately described by their CSD, W ( r 1 , r 2 , 0 ) [1], which gives the second-order correlations of the field at two different points r 1 and r 2 of the source plane. Here and in the following we omit the explicit dependence of W on the temporal frequency. From such CSD function, the intensity can be derived by evaluating it at coincident points, that is,
I ( r , 0 ) = W ( r , r , 0 ) .
A normalized version of the CSD function is the complex degree of coherence, which is defined as
γ ( r 1 , r 2 , 0 ) = W ( r 1 , r 2 , 0 ) I ( r 1 , 0 ) I ( r 2 , 0 ) ,
whose absolute value ranges form 0 to 1, the two extremes denoting (spatial) complete incoherence and perfect coherence, respectively.
According to Mercer’s theorem [1,22], any valid CSD function can be expanded as a series of the following type
W ( r 1 , r 2 , 0 ) = n , m λ n m Φ n m * ( r 1 , 0 ) Φ n m ( r 2 , 0 ) ,
where n and m are indices denoting the eigenfunctions Φ n m ( r , 0 ) of the homogeneous Fredholm integral equation of the second kind whose kernel is W ( r 1 , r 2 , 0 ) , and λ n m the corresponding eigenvalues [1]. Due to the non-negative definiteness of W, all the eigenvalues are nonnegative. Using this modal expansion, the intensity distribution across the source can be expressed as
I ( r , 0 ) = n , m λ n m Φ n m ( r , 0 ) 2 .
The use of the modes allows one to evaluate the CSD of the field propagated from the source in an easy way, as the sum of the CSDs of the propagated modes. In fact we have
W ( r 1 , r 2 , z ) = n , m λ n m Φ n m * ( r 1 , z ) Φ n m ( r 2 , z )
and, consequently,
I ( r , z ) = n , m λ n m Φ n m ( r , z ) 2 .

2.1. Shape Invariant Partially Coherent Fields: LG Modes

A particular set of orthonormal mode is formed by the LG functions, defined as
Φ n m ( r , 0 ) = 1 w 0 2 n ! π ( n + | m | ) ! e r 2 / w 0 2 e i m φ 2 r w 0 | m | L n | m | 2 r 2 w 0 2 ,
where n = 0 , 1 , , , m = 0 , ± 1 , , ± , r = ( r , φ ) , and L n m are the associated Laguerre polynomials [23,24]. The parameter w 0 (the spot size) specifies the transverse width of the functions.
By substituting Equation (7) into the Mercer’s expansion of Equation (3), a bona-fide CSD is obtained as
W ( r 1 , r 2 , 0 ) = 2 π w 0 2 e ( r 1 2 + r 2 2 ) / w 0 2 × n , m λ n m n ! ( n + | m | ) ! 2 r 1 r 2 w 0 2 | m | L n | m | 2 r 1 2 w 0 2 L n | m | 2 r 2 2 w 0 2 e i m ( φ 1 φ 2 ) .
Assuming paraxial approximation, the free-space propagated LG modes can be calculated as [23]
Φ n m ( r , z ) = 1 w z 2 n ! π ( n + | m | ) ! e r 2 / w z 2 e i k r 2 / ( 2 R z ) e i m φ × e i [ k z ( 2 n + | m | + 1 ) α z ] 2 r w z | m | L n | m | 2 r 2 w z 2 ,
where the parameters w z , R z and α z are
w z = w 0 1 + z z R 2 ; R z = z 1 + z R z 2 ; α z = arctan z z R ,
z R = k w 0 2 / 2 being the Rayleigh distance. From Equation (5) the resulting propagated CSD turns out to be
W ( r 1 , r 2 , z ) = 2 π w z 2 e ( r 1 2 + r 2 2 ) / w z 2 e i k ( r 1 2 r 2 2 ) / ( 2 R z ) × n , m λ n m n ! ( n + | m | ) ! 2 r 1 r 2 w z 2 | m | L n | m | 2 r 1 2 w z 2 L n | m | 2 r 2 2 w z 2 e i m ( φ 1 φ 2 ) .
On comparing Equations (8) and (11) we see that, except for a spherical curvature term, the CSD maintains exactly the same form across any plane z = constant, provided that w 0 is replaced by w z . This is briefly expressed by saying that the partially coherent fields radiated by our sources are shape invariant. In particular, on letting r 1 = r 2 = r in Equation (11), we obtain the intensity
I ( r , z ) = 2 π w z 2 e 2 r 2 / w z 2 n , m λ n m n ! ( n + | m | ) ! 2 r 2 w z 2 | m | L n | m | 2 r 2 w z 2 2 ,
which, besides maintaining the same shape under propagation, presents rotational invariance.
By virtue of the above property we can refer most of our considerations to the plane z = 0, omit the explicit dependence on z of the functions and, for brevity, use the symbol w (without any subscript) for the spot size of the modes.

2.2. Some Closed-Form CSDs with LG Modes

The most celebrated example of shape-invariant partially coherent beams are, without any doubt, the Gaussian Schell-model (GSM) beams. They were introduced more that forty years ago [2,3,4,5] and since then they have represented in a countless number of cases the archetype of partially coherent beams. This is mainly due to the simple mathematical form of their CSD.
The shape invariance of such beams can be traced back, of course, to the mathematical form of their coherent modes, which are, in fact, Hermite–Gaussian [10,11] or, equivalently, LG functions [25]. On putting eigenvalues of the form (see, for example, [26])
λ n m = λ 0 β 2 n + | m | ( 0 β 1 ) ,
with λ 0 and β positive parameters, into Equation (3), with the modes given in Equation (7), the following CSD is obtained:
W ( r 1 , r 2 ) = A e a ( r 1 2 + r 2 2 ) e b ( r 1 r 2 ) 2 ,
where A is a proportionality factor having dimensions of an intensity, while
a = 1 w 2 1 β 1 + β ; b = 2 w 2 β 1 β 2
are two parameters related to the widths of the intensity profile and of the degree of coherence, respectively, which are
I ( r ) = A e 2 a r 2 ; γ ( r 1 , r 2 ) = e b ( r 1 r 2 ) 2 .
Therefore, both the intensity and the degree of coherence of such sources are described by Gaussian functions. Furthermore, the degree of coherence is shift-invariant across the transverse section of the beam, that is, it only depends on the difference of the position vectors of the two considered points. Such a characteristic identifies the so-called Schell-model sources. It is interesting to note that the coherence area across the source (proportional to 1 / b ) reduces on increasing β , i.e., on increasing the number of modes that significantly contribute to the source [1]. This is a rather general result of the modal theory of coherence, although not universal (see Section 3).
An important generalization of the class of GSM source, due to Simon and Mukunda [26,27,28], is represented by the Twisted Gaussian Schell-model (TGSM) sources, whose CSD reads
W ( r 1 , r 2 ) = A e a ( r 1 2 + r 2 2 ) e b ( r 1 r 2 ) 2 e i u ( r 1 × r 2 ) ,
where u is a real parameter called the twist parameter and the subscript ⊥ denotes the z-component of a vector. Physical interpretations and possible realization schemes of these twisted sources can be found in Refs. [29,30]. The corresponding intensity profile and degree of coherence take the forms
I ( r ) = A e 2 a r 2 ; γ ( r 1 , r 2 ) = e b ( r 1 r 2 ) 2 e i u ( r 1 × r 2 ) .
We can notice that the absolute value of γ is still shift-invariant, but γ itself is not.
Simon and Mukunda introduced such sources on the basis of symmetry arguments and also found their coherent modes and eigenvalues. Even in this case the modes are LG functions, with eigenvalues of the form
λ n m = λ 0 β 2 n + | m | μ | m | / 2 ( 0 β 1 ; 0 μ 1 ) ,
and this makes such sources of the shape-invariant type. The relations between (a, b) and (w, β ) are still given by Equation (15), while
u = 2 b μ 1 μ + 1 ,
which implies
| u | 2 b ,
that is, the maximum value of the twist parameter is bounded by the inverse of the (squared) width of the degree of coherence. Therefore, no twist can be present in a perfectly coherent source [26,27,28,29,30]. Research on partially coherent sources presenting a twist is still very active [31,32,33,34,35,36,37,38,39,40,41,42].
The next example we recall is that of the partially coherent Flattened Gaussian beams. They were introduced first in their coherent version [43], to provide a simple model to manage the paraxial propagation of light beams with a flattened profile. Their partially coherent counterpart was presented as an example of the superposition of mutually uncorrelated modes with helicoidal phases [44]. In that case, only a finite number of LG modes, all having n = 0 , with equal eigenvalues was chosen. In particular, if
λ n m = λ 0 δ n , 0 × 1 ( m = 0 , 1 , 2 , , M ) 0 ( m > M )
with δ n , n the Kronecker symbol, the following CSD is obtained:
W M ( r 1 , r 2 ) = A M ! e ( r 1 2 + r 2 2 ) / w 2 exp 2 r 1 r 2 e i ( φ 1 φ 2 ) w 2 Γ M + 1 2 r 1 r 2 e i ( φ 1 φ 2 ) w 2 ,
where Γ M + 1 ( x ) is the incomplete Gamma function [24]. The corresponding intensity reads
I M ( r ) = A M ! Γ M + 1 2 r 2 w 2 .
Plots of the intensity and degree of coherence for such kind of sources are shown in [44]. While the former always presents a flat-topped circularly symmetric profile, interesting spatial structures can be observed for the degree of coherence.
To conclude this roundup of examples of closed-form CSDs with LG modes, we recall the one introduced by Ponomarenko in Ref. [45], which has been also recently experimentally synthesized [46]. It is obtained as the incoherent superposition of an infinite number of LG modes of different radial order, but with the same azimuthal index m. More precisely, the coefficients decrease geometrically with n, i.e.,
λ n m = λ 0 δ m , m ξ n
with 0 < ξ < 1 and m = 0 , ± 1 , , ± .
The resulting CSD is of the form
W m ( r 1 , r 2 ) = A e a ( r 1 2 + r 2 2 ) e i m ( φ 1 φ 2 ) I m c r 1 r 2 ,
with
a = 1 w 2 1 + ξ 1 ξ ; c = 1 w 2 4 ξ 1 ξ
and I m is the modified Bessel function of the first kind and order m [24]. Note that in this case the modulus of the spectral degree of coherence is independent of the relative orientation of the points, i.e., it is circularly symmetric.

3. Laguerre-Christoffel-Darboux Sources

Here, we introduce a new type of shape-invariant CSDs whose expression can be given in closed form. We first limit ourselves to the case where the angular index (m) of the LG functions in the Mercer expansion in Equation (3) is kept fixed and take a finite number of modes with different n and equal eigenvalues. In particular we take
λ n m = λ 0 δ m , m × 1 ( n = 0 , 1 , 2 , , N ) ; 0 ( n > N ) .
with m = 0 , ± 1 , , ± .
Using the expression of the modes in Equation (7), the following expression for the resulting CSD turns out:
W N m ( r 1 , r 2 ) = A e ( r 1 2 + r 2 2 ) / w 2 e i m ( φ 1 φ 2 ) 2 r 1 r 2 w 2 | m | G N m ( r 1 , r 2 ) ,
where
G N m ( r 1 , r 2 ) = n = 0 N n ! ( n + | m | ) ! L n | m | 2 r 1 2 w 2 L n | m | 2 r 2 2 w 2 .
Note that the case N = 0 corresponds to a single Laguerre–Gaussian mode, so that the source is perfectly coherent.
The effect of taking the above coefficients is that the sum in Equation (30) can be calculated exactly on exploiting the Christoffel-Darboux formula [24,47]. The latter holds for any system of polynomials { f n ( x ) } on the (finite or infinite) interval ( a , b ) , which are orthonormal with respect to the non-negative weight function w ( x ) . The theorem states that, if k n is the coefficient of x n in the polynomial and
h n = a b f n 2 ( x ) w ( x ) d x ,
then
n = 0 N f n ( x ) f n ( y ) h n = k N k N + 1 h N f N + 1 ( x ) f N ( y ) f N ( x ) f N + 1 ( y ) x y .
Using the above result, Equation (30) can be written as
G N m ( r 1 , r 2 ) = ( N + 1 ) ! ( N + | m | ) ! L N | m | 2 r 1 2 w 2 L N + 1 | m | 2 r 2 2 w 2 L N + 1 | m | 2 r 1 2 w 2 L N | m | 2 r 2 2 w 2 2 w 2 ( r 1 2 r 2 2 )
if r 1 r 2 . In the limit of r 1 = r 2 = r , the resulting indeterminate can be evaluated by means of the de l’Hôpital theorem [24], which gives
G N m ( r , r ) = ( N + 1 ) ! ( N + | m | ) ! L N | m | 2 r 2 w 2 L N | m | + 1 2 r 2 w 2 L N + 1 | m | 2 r 2 w 2 L N 1 | m | + 1 2 r 2 w 2 .
A source characterized by the CSD in Equation (29), with the function G N m given in Equation (33) will be called a Laguerre-Christoffel-Darboux (LCD) source. Sources of this class could be experimentally generated following the procedures described in Refs. [41,46].
From Equations (29), with r 1 = r 2 = r , and (34) we obtain the following expression for the intensity:
I N m ( r ) = 2 π w 2 e 2 r 2 / w 2 2 r 2 w 2 | m | ( N + 1 ) ! ( N + | m | ) ! × L N | m | 2 r 2 w 2 L N | m | + 1 2 r 2 w 2 L N + 1 | m | 2 r 2 w 2 L N 1 | m | + 1 2 r 2 w 2 ,
whose behavior is shown in Figure 1 and Figure 2 for some values of N and m. Figure 1 shows the changes in the intensity profile when m is kept fixed and the number of modes is varied. It can be noticed that, due to the presence of the optical vortex, the intensity profile shows a central dark region where the intensity is practically zero. The radius of this circular dark region increases with m and decreases with N. We note that the presence of a central dark zone with a fast variation of the intensity at the border of this zone could be efficiently used in particle trapping. Furthermore, N + 1 relative maxima are observed in the axial profile, which become less and less evident on increasing the value of the indices. Some intensity profiles across the source plane are shown in Figure 2. Due to the shape invariance of the LCD sources, Figure 1 and Figure 2 are valid at any transverse plane, the only change being the scale factor w that varies with z [see Equation (10)]. The same holds for the forthcoming figures of intensity and degree of coherence.
Interesting results are obtained on considering the limit of the expression in Equation (35) for large values of N. Using Equation (8.978.2) of Ref. [48], that is,
lim N 1 N m L N m ξ N = ξ m / 2 J m ( 2 ξ ) ,
where J m is the Bessel function of the first kind and order m [24], the following asymptotic expression is obtained for the source intensity:
I N m ( r ) 2 π w e 2 r 2 / w 2 ( N + 1 ) ! ( N + | m | ) ! 1 r N | m | + 1 / 2 J | m | 2 r 2 N w J | m | + 1 2 r 2 N w N 1 ( N 2 1 ) | m | / 2 J | m | 2 r 2 ( N + 1 ) w J | m | + 1 2 r 2 ( N 1 ) w .
The usefulness of introducing the Bessel functions is apparent if we consider their property according to which J k ( x ) is negligibly small for x < k [24,49]. This allows us to roughly estimate the radius of the central dark region, say r 0 . In fact, if we neglect the difference between m and m + 1 , as well as that between N and N + 1 , we can state that all the Bessel functions appearing in Equation (37) vanishes if
2 r 2 N w | m | ,
so that
r 0 | m | 2 2 N w .
Radial profiles and 2D plots of the intensity profile are shown in Figure 3 for N = 100 and m = 25 , 50, and 100, corresponding to r 0 values of 0.9 w, 1.8 w, and 3.5 w, respectively. It is worth noting that the intensity drops practically to zero outside the donut region. For example, the intensity I 100 , 25 ( r ) is below 10 4 of its maximum for r < 0.6 w or r > 15.5 w. The same happens for I 100 , 50 ( r ) if r < 1.3 w or r > 16.3 w, and for I 100 , 100 ( r ) if r < 2.6 w or r > 17.5 w. The inset in Figure 3 shows the ripples of the intensity profile, that become smoother and smoother on increasing N.
As far as the complex degree of coherence is concerned, we have
γ N m ( r 1 , r 2 ) = G N m ( r 1 , r 2 ) e i m ( φ 2 φ 1 ) G N m ( r 1 , r 1 ) G N m ( r 2 , r 2 ) .
It is seen that the absolute value of γ N m depends only on the radial distances of the considered points from the source center, so that the source exhibits perfect coherence along any annulus that is concentric to the source center, i.e.,
γ N m ( r , φ 1 ; r , φ 2 ) = 1
for any choice of m, φ 1 and φ 2 . In this sense, the sources present circular coherence [49,50]. Conversely, coherence can be partial or even vanishing between two points at different distances from the center along a radius.
Plots of the absolute value of the degree of coherence across the plane ( r 1 , r 2 ) are shown in Figure 4. It can be observed that setting r 2 ( r 1 ) close to zero and varying r 1 ( r 2 ), it passes through N zeros. The same happens for large values of r 2 ( r 1 ), i.e., at points outside the intensity donut. The positions of such minima along r 1 ( r 2 ) remain almost the same on varying r 2 ( r 1 ) and are located within the intensity donut (for m 0 ).
Maps of the degree of coherence across the source plane can be seen in Figure 5, which shows the absolute value of γ N m ( r 1 , r 2 ) as a function of r 1 for three different values of r 2 . The three positions r 2 are denoted by small dots. Chosen parameters are N = 2 and m = 4 , for which the intensity map is the one shown in the center part of Figure 2.
Asymptotic expressions of the degree of coherence can be also deduced from Equations (33) and (40). In fact, keeping r 2 fixed and increasing r 1 we obtain
lim r 1 γ N m ( r 1 , φ 1 ; r 2 , φ 2 ) = 1 N + 1 L N | m | 2 r 2 2 w 2 L N | m | 2 r 2 2 w 2 L N | m | + 1 2 r 2 2 w 2 L N + 1 | m | 2 r 2 2 w 2 L N 1 | m | + 1 2 r 2 2 w 2 .
In particular, if r 2 is close to zero,
lim r 1 r 2 0 γ N m ( r 1 , φ 1 ; r 2 , φ 2 ) = lim r 2 r 1 0 γ N m ( r 1 , φ 1 ; r 2 , φ 2 ) = | m | + 1 N + | m | + 1
The absolute value of the degree of coherence along a radius is also represented in Figure 6 as a function of r 1 for three different values of r 2 . The order N is fixed (equal to 5) and several values of m are considered ( m = 5 , 10 , 20 ). The three values of r 2 are chosen as follows: the first one corresponds to a point close to the inner border of the intensity donut; the second one is in the central part of the donut; and the third one is near the outer border of the donut. For any of these sources, complete correlation is observed for any pair of points located on the same circle concentric with the source center, as it was expected from Equation (41). Furthermore, within the region where the intensity is not negligible, the absolute value of the degree of coherence shows N minima, where it vanishes, meaning that there could be complete incoherence between points belonging to different circles. Outside that region the absolute degree of coherence tends to a limit given by Equation (42).
Similar plots are shown in Figure 7 for a source where a larger number of modes are considered ( N = 20 ) with topological charge m = 0 , 1 , 2 . For these cases, it can be observed that the degree of coherence is considerably high for pairs of points in the region where the intensity is significant. For example, for m = 0 the degree of coherence is over 0.75 for any pair of points located within a circle of radius 0.2 w where the intensity is over 0.2 of its maximum value. That is, a relatively highly coherent source is obtained, despite we are considering the incoherent superposition of a high number of modes. A similar behavior is observed for topological charge m = 1 and m = 2 and a ring around the circle where the maximum intensity is reached.
A source with a more structured degree of coherence is obtained if further terms are included into the Mercer expansion. In particular, if we take the modes of orders both ( n , m ) and ( n , m ) , with the same weight, i.e.,
λ n m = λ 0 δ m , m + δ m , m 2 × 1 ( n = 0 , 1 , 2 , , N ) 0 ( n > N )
the intensity profile does not change with respect to the previous case, but an angular modulation is introduced into the degree of coherence, being
γ N m ( r 1 , r 2 ) = G N m ( r 1 , r 2 ) cos [ m ( φ 2 φ 1 ) ] G N m ( r 1 , r 1 ) G N m ( r 2 , r 2 ) .
Figure 8 and Figure 9 show the behavior of the degree of coherence for several sources of this type. In this case the absolute value of the degree of coherence does not only depend on the radial coordinates but also on the angle difference of the considered points, except when m = 0 . In fact, if m 0 the degree of coherence is zero for points that satisfy φ 2 φ 1 = π ( + 1 / 2 ) / m , while is unitary for points with φ 1 φ 2 = π / m . For other angles, the degree of coherence is modulated by a sinusoidal function.
Note that sources described by Equation (8) with the eigenvalues given by Equation (28) or Equation (44), with odd m, present a degree of coherence that reaches the value 1 for any pair of diametrically opposite points on the same concentric circle (see Figure 9). Partially coherent beams with analogous characteristics have been recently considered, both theoretically and experimentally, to overcome the classical Rayleigh diffraction limit in the imaging of two holes (see for example [51,52,53,54]).

4. Conclusions

In this paper we have introduced a new class of partially coherent fields whose cross-spectral density is shape invariant upon paraxial propagation. This means that both their transverse intensity profile and the structure of their degree of spectral coherence keep their form unchanged during propagation, except for a transverse scaling factor.
Fields of this class are obtained on superimposing a finite number of mutually uncorrelated Laguerre–Gaussian modes, having fixed topological charge and equal power, but the resulting CSD can be always expressed in a simple mathematical form regardless the number of involved modes. The closed-form expression of the CSD can be exploited for deducing properties of the beams analytically, or for making the problem of their propagation easier when a numerical approach has to be used (e.g., in non-ABCD or non-deterministic systems).
The intensities of the beams we have considered have rotational symmetry and (for non vanishing topological charge of the component modes) have a donut-like profile, whose size can be varied at will on changing the parameters of the modes. An interesting feature of these profiles is that the intensity is uniform and practically zero in the whole regions inside the donut, thus mimicking a circular-well potential in particle trapping experiments.
The degree of coherence presents a vortex-phase structure and its modulus only depends on the distances of the involved points from the propagation axis. It is interesting to remark that, contrary to what usually happens, for low values of the topological charge (including the case of vanishing charge), the coherence area across the transverse plane remains relatively large in whole region where the intensity is significant, even if a very large number of component modes is considered.
The vortex structure of the degree of coherence can be converted into an angularly periodic modulation on superimposing modes with opposite charge, leaving the radial dependence unchanged. For odd values of the charge, the degree of coherence equals 1 for any pairs of diametrically opposite points on a circle r = constant.
Sources with the above characteristics could have potential applications in several fields, such as particle trapping, optical communications, imaging, and more.

Author Contributions

All authors contributed equally. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Ministerio de Economía y Competitividad, grant number PID2019-104268GB-C21.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors wish to thank F. Gori for his useful suggestions.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
CSDCross-spectral density
GSMGaussian Schell-model
LGLaguerre–Gaussian
LCDLaguerre–Christoffel–Darboux
OAMOrbital angular momentum
TGSMTwisted Gaussian Schell-model

References

  1. Mandel, L.; Wolf, E. Optical Coherence and Quantum Optics; Cambridge University Press: Cambridge, UK, 1995. [Google Scholar] [CrossRef]
  2. Collett, E.; Wolf, E. Is complete spatial coherence necessary for the generation of highly directional light beams? Opt. Lett. 1978, 2, 27–29. [Google Scholar] [CrossRef]
  3. Wolf, E.; Collett, E. Partially coherent sources which produce the same far-field intensity distribution as a laser. Opt. Commun. 1978, 25, 293–296. [Google Scholar] [CrossRef]
  4. Foley, J.T.; Zubairy, M. The directionality of gaussian Schell-model beams. Opt. Commun. 1978, 26, 297–300. [Google Scholar] [CrossRef]
  5. Gori, F.; Palma, C. Partially coherent sources which give rise to highly directional light beams. Opt. Commun. 1978, 27, 185–188. [Google Scholar] [CrossRef]
  6. Gori, F.; Santarsiero, M. Devising genuine spatial correlation functions. Opt. Lett. 2007, 32, 3531–3533. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Martínez-Herrero, R.; Mejías, P.M.; Gori, F. Genuine cross-spectral densities and pseudo-modal expansions. Opt. Lett. 2009, 34, 1399–1401. [Google Scholar] [CrossRef] [PubMed]
  8. Wolf, E. New theory of partial coherence in the space–frequency domain. Part I: Spectra and cross spectra of steady-state sources. J. Opt. Soc. Am. 1982, 72, 343–351. [Google Scholar] [CrossRef]
  9. Martínez-Herrero, R. Expansion of complex degree of coherence. II Nuovo C. B 1979, 54, 205–210. [Google Scholar]
  10. Gori, F. Collett-Wolf sources and multimode lasers. Opt. Commun. 1980, 34, 301–305. [Google Scholar] [CrossRef]
  11. Starikov, A.; Wolf, E. Coherent-mode representation of Gaussian Schell-model sources and of their radiation fields. J. Opt. Soc. Am. 1982, 72, 923–928. [Google Scholar] [CrossRef]
  12. Cohen-Tannoudji, C.; Diu, B.; Laloë, F. Quantum Mechanics; Trans. of: Mécanique quantique. Paris: Hermann, 1973; Wiley: New York, NY, USA, 1977. [Google Scholar]
  13. Martínez-Herrero, R.; Gori, F. Christoffel–Darboux sources. Opt. Lett. 2021, 46, 973–976. [Google Scholar] [CrossRef]
  14. Soskin, M.S.; Gorshkov, V.N.; Vasnetsov, M.V.; Malos, J.T.; Heckenberg, N.R. Topological charge and angular momentum of light beams carrying optical vortices. Phys. Rev. A 1997, 56, 4064–4075. [Google Scholar] [CrossRef] [Green Version]
  15. Allen, L.; Padgett, M.; Babiker, M. IV The Orbital Angular Momentum of Light; Progress in Optics; Elsevier: Amsterdam, The Netherlands, 1999; Volume 39, pp. 291–372. [Google Scholar] [CrossRef]
  16. Allen, L.; Barnett, S.M.; Padgett, M.J. Optical Angular Momentum; CRC Press: Boca Raton, FL, USA, 2003. [Google Scholar] [CrossRef]
  17. Padgett, M.; Courtial, J.; Allen, L. Light’s Orbital Angular Momentum. Phys. Today 2004, 57, 35–40. [Google Scholar] [CrossRef]
  18. Zeng, J.; Lin, R.; Liu, X.; Zhao, C.; Cai, Y. Review on partially coherent vortex beams. Front. Optoelectron. 2019, 12, 229–248. [Google Scholar] [CrossRef]
  19. Shen, Y.; Wang, X.; Xie, Z.; Min, C.; Fu, X.; Liu, Q.; Gong, M.; Yuan, X. Optical vortices 30 years on: OAM manipulation from topological charge to multiple singularities. Light Sci. Appl. 2019, 8, 90. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Zhang, Y.; Cai, Y.; Gbur, G. Control of orbital angular momentum with partially coherent vortex beams. Opt. Lett. 2019, 44, 3617–3620. [Google Scholar] [CrossRef]
  21. Zhang, Y.; Korotkova, O.; Cai, Y.; Gbur, G. Correlation-induced orbital angular momentum changes. Phys. Rev. A 2020, 102, 063513. [Google Scholar] [CrossRef]
  22. Mercer, J. Functions of Positive and Negative Type, and their Connection with the Theory of Integral Equations. Philos. Trans. R. Soc. Lond. Ser. A 1909, 209, 415–446. [Google Scholar]
  23. Siegman, A.E. Lasers; University Science Books: Mill Valley, CA, USA, 1986. [Google Scholar]
  24. Abramowitz, M.; Stegun, I. (Eds.) Handbook of Mathematical Functions; Dover Publications Inc.: Mineola, NY, USA, 1972. [Google Scholar]
  25. Abramochkin, E.; Volostnikov, V. Beam transformations and nontransformed beams. Opt. Commun. 1991, 83, 123–135. [Google Scholar] [CrossRef]
  26. Simon, R.; Sundar, K.; Mukunda, N. Twisted Gaussian Schell-model beams. I. Symmetry structure and normal-mode spectrum. J. Opt. Soc. Am. A 1993, 10, 2008–2016. [Google Scholar] [CrossRef]
  27. Simon, R.; Mukunda, N. Twisted Gaussian Schell-model beams. J. Opt. Soc. Am. A 1993, 10, 95–109. [Google Scholar] [CrossRef]
  28. Sundar, K.; Simon, R.; Mukunda, N. Twisted Gaussian Schell-model beams. II. Spectrum analysis and propagation characteristics. J. Opt. Soc. Am. A 1993, 10, 2017–2023. [Google Scholar] [CrossRef]
  29. Ambrosini, D.; Bagini, V.; Gori, F.; Santarsiero, M. Twisted Gaussian Schell-model Beams: A Superposition Model. J. Mod. Opt. 1994, 41, 1391–1399. [Google Scholar] [CrossRef]
  30. Friberg, A.; Tervonen, E.; Turunen, J. Interpretation and experimental demonstration of twisted Gaussian Schell-model beams. J. Opt. Soc. Am. A 1994, 11, 1818–1826. [Google Scholar] [CrossRef]
  31. Mirhosseini, M.; Magaña-Loaiza, O.S.; O’Sullivan, M.N.; Rodenburg, B.; Malik, M.; Lavery, M.P.J.; Padgett, M.J.; Gauthier, D.J.; Boyd, R.W. High-dimensional quantum cryptography with twisted light. New J. Phys. 2015, 17, 033033. [Google Scholar] [CrossRef]
  32. Gori, F.; Santarsiero, M. Twisted Gaussian Schell-model beams as series of partially coherent modified Bessel-Gauss beams. Opt. Lett. 2015, 40, 1587–1590. [Google Scholar] [CrossRef]
  33. Borghi, R.; Gori, F.; Guattari, G.; Santarsiero, M. Twisted Schell-model beams with axial symmetry. Opt. Lett. 2015, 40, 4504–4507. [Google Scholar] [CrossRef]
  34. Mei, Z.; Korotkova, O. Random sources for rotating spectral densities. Opt. Lett. 2017, 42, 255–258. [Google Scholar] [CrossRef] [PubMed]
  35. Gori, F.; Santarsiero, M. Devising genuine twisted cross-spectral densities. Opt. Lett. 2018, 43, 595–598. [Google Scholar] [CrossRef] [PubMed]
  36. Borghi, R. Twisting partially coherent light. Opt. Lett. 2018, 43, 1627–1630. [Google Scholar] [CrossRef]
  37. Wang, H.; Peng, X.; Liu, L.; Wang, F.; Cai, Y.; Ponomarenko, S.A. Generating bona fide twisted Gaussian Schell-model beams. Opt. Lett. 2019, 44, 3709–3712. [Google Scholar] [CrossRef] [PubMed]
  38. Santarsiero, M.; Gori, F.; Alonzo, M. Higher-order twisted/astigmatic Gaussian Schell-model cross-spectral densities and their separability features. Opt. Express 2019, 27, 8554–8565. [Google Scholar] [CrossRef] [PubMed]
  39. Hyde, M.W. Twisted space-frequency and space-time partially coherent beams. Sci. Rep. 2020, 10, 12443. [Google Scholar] [CrossRef]
  40. Peng, X.; Wang, H.; Liu, L.; Wang, F.; Popov, S.; Cai, Y. Self-reconstruction of twisted Laguerre–Gaussian Schell-model beams partially blocked by an opaque obstacle. Opt. Express 2020, 28, 31510–31523. [Google Scholar] [CrossRef] [PubMed]
  41. Tian, C.; Zhu, S.; Huang, H.; Cai, Y.; Li, Z. Customizing twisted Schell-model beams. Opt. Lett. 2020, 45, 5880–5883. [Google Scholar] [CrossRef] [PubMed]
  42. Liu, Z.; Wan, L.; Zhou, Y.; Zhang, Y.; Zhao, D. Progress on Studies of Beams Carrying Twist. Photonics 2021, 8, 92. [Google Scholar] [CrossRef]
  43. Gori, F. Flattened Gaussian beams. Opt. Commun. 1994, 107, 335–341. [Google Scholar] [CrossRef]
  44. Gori, F.; Santarsiero, M.; Borghi, R.; Vicalvi, S. Partially coherent sources with helicoidal modes. J. Mod. Opt. 1998, 45, 539–554. [Google Scholar] [CrossRef]
  45. Ponomarenko, S.A. A class of partially coherent beams carrying optical vortices. J. Opt. Soc. Am. A 2001, 18, 150–156. [Google Scholar] [CrossRef]
  46. Chen, X.; Li, J.; Rafsanjani, S.M.H.; Korotkova, O. Synthesis of Im-Bessel correlated beams via coherent modes. Opt. Lett. 2018, 43, 3590–3593. [Google Scholar] [CrossRef]
  47. Szegö, G. Orthogonal Polynomials; American Mathematical Society: Providence, RI, USA, 1975. [Google Scholar]
  48. Gradshteyn, I.S.; Ryzhik, I.M. Table of Integrals, Series, and Products, 4th ed.; Academic Press: Cambridge, MA, USA, 1965. [Google Scholar]
  49. Santarsiero, M.; Martínez-Herrero, R.; Maluenda, D.; de Sande, J.C.G.; Piquero, G.; Gori, F. Synthesis of circularly coherent sources. Opt. Lett. 2017, 42, 4115–4118. [Google Scholar] [CrossRef] [PubMed]
  50. Santarsiero, M.; Martínez-Herrero, R.; Maluenda, D.; de Sande, J.C.G.; Piquero, G.; Gori, F. Partially coherent sources with circular coherence. Opt. Lett. 2017, 42, 1512–1515. [Google Scholar] [CrossRef] [Green Version]
  51. Tamburini, F.; Anzolin, G.; Umbriaco, G.; Bianchini, A.; Barbieri, C. Overcoming the Rayleigh Criterion Limit with Optical Vortices. Phys. Rev. Lett. 2006, 97, 163903. [Google Scholar] [CrossRef] [Green Version]
  52. Tong, Z.; Korotkova, O. Beyond the classical Rayleigh limit with twisted light. Opt. Lett. 2012, 37, 2595–2597. [Google Scholar] [CrossRef] [PubMed]
  53. Liang, C.; Wu, G.; Wang, F.; Li, W.; Cai, Y.; Ponomarenko, S.A. Overcoming the classical Rayleigh diffraction limit by controlling two-point correlations of partially coherent light sources. Opt. Express 2017, 25, 28352–28362. [Google Scholar] [CrossRef]
  54. Liang, C.; Monfared, Y.E.; Liu, X.; Qi, B.; Wang, F.; Korotkova, O.; Cai, Y. Optimizing illumination’s complex coherence state for overcoming Rayleigh’s resolution limit. Chin. Opt. Lett. 2021, 19, 052601. [Google Scholar]
Figure 1. Intensity profiles I N m given by Equation (35) for several values of N and m (as indicated in Figure labels).
Figure 1. Intensity profiles I N m given by Equation (35) for several values of N and m (as indicated in Figure labels).
Photonics 08 00134 g001
Figure 2. Maps of the intensity I N m (normalized to the maximum) across the source plane for three different pairs ( N , m ) (as indicated in the figure labels).
Figure 2. Maps of the intensity I N m (normalized to the maximum) across the source plane for three different pairs ( N , m ) (as indicated in the figure labels).
Photonics 08 00134 g002
Figure 3. Intensity profiles as functions of the radial coordinate and intensity maps (normalized to its maximum) for N = 100 and m = 25 , 50 , 100 .
Figure 3. Intensity profiles as functions of the radial coordinate and intensity maps (normalized to its maximum) for N = 100 and m = 25 , 50 , 100 .
Photonics 08 00134 g003
Figure 4. Absolute value of the complex degree of coherence and funtion in the r 1 r 2 plane for N = 1 , 2 , 5 , 10 (from left to right) and m = 0 (upper row) and m = 5 (lower row).
Figure 4. Absolute value of the complex degree of coherence and funtion in the r 1 r 2 plane for N = 1 , 2 , 5 , 10 (from left to right) and m = 0 (upper row) and m = 5 (lower row).
Photonics 08 00134 g004
Figure 5. Absolute degree of coherence (Equation (40)) relative to three different points, indicated by red dots, for N = 2 , m = 4 .
Figure 5. Absolute degree of coherence (Equation (40)) relative to three different points, indicated by red dots, for N = 2 , m = 4 .
Photonics 08 00134 g005
Figure 6. Upper row: intensity profile given by Equation (35) for N = 5 and several values of m. Second, third and fourth rows: absolute value of the degree of coherence for relative to a point located at a distance r 2 from the source center (indicated with an arrow on the intensity profile and a dashed line in the degree of coherence profile).
Figure 6. Upper row: intensity profile given by Equation (35) for N = 5 and several values of m. Second, third and fourth rows: absolute value of the degree of coherence for relative to a point located at a distance r 2 from the source center (indicated with an arrow on the intensity profile and a dashed line in the degree of coherence profile).
Photonics 08 00134 g006
Figure 7. Upper row: intensity profile given by Equation (35) for N = 20 and several values of m. Second, third and fourth rows: absolute value of the degree of coherence for relative to a point located at a distance r 2 from the source center (indicated with an arrow on the intensity profile and a dashed line in the degree of coherence profile).
Figure 7. Upper row: intensity profile given by Equation (35) for N = 20 and several values of m. Second, third and fourth rows: absolute value of the degree of coherence for relative to a point located at a distance r 2 from the source center (indicated with an arrow on the intensity profile and a dashed line in the degree of coherence profile).
Photonics 08 00134 g007
Figure 8. Degree of coherence for a source described by Equations (44) and (45), N = 2 , m = 4 , relative to three different points indicated by white dots and in the labels.
Figure 8. Degree of coherence for a source described by Equations (44) and (45), N = 2 , m = 4 , relative to three different points indicated by white dots and in the labels.
Photonics 08 00134 g008
Figure 9. Degree of coherence for a source described by Equations (44) and (45), N = 3 , m = 3 , relative to three different points indicated by white dots and in the labels.
Figure 9. Degree of coherence for a source described by Equations (44) and (45), N = 3 , m = 3 , relative to three different points indicated by white dots and in the labels.
Photonics 08 00134 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Martínez-Herrero, R.; Santarsiero, M.; Piquero, G.; González de Sande, J.C. A New Type of Shape-Invariant Beams with Structured Coherence: Laguerre-Christoffel-Darboux Beams. Photonics 2021, 8, 134. https://doi.org/10.3390/photonics8040134

AMA Style

Martínez-Herrero R, Santarsiero M, Piquero G, González de Sande JC. A New Type of Shape-Invariant Beams with Structured Coherence: Laguerre-Christoffel-Darboux Beams. Photonics. 2021; 8(4):134. https://doi.org/10.3390/photonics8040134

Chicago/Turabian Style

Martínez-Herrero, Rosario, Massimo Santarsiero, Gemma Piquero, and Juan Carlos González de Sande. 2021. "A New Type of Shape-Invariant Beams with Structured Coherence: Laguerre-Christoffel-Darboux Beams" Photonics 8, no. 4: 134. https://doi.org/10.3390/photonics8040134

APA Style

Martínez-Herrero, R., Santarsiero, M., Piquero, G., & González de Sande, J. C. (2021). A New Type of Shape-Invariant Beams with Structured Coherence: Laguerre-Christoffel-Darboux Beams. Photonics, 8(4), 134. https://doi.org/10.3390/photonics8040134

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop