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,
[
1], which gives the second-order correlations of the field at two different points
and
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,
A normalized version of the CSD function is the complex degree of coherence, which is defined as
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
where
n and
m are indices denoting the eigenfunctions
of the homogeneous Fredholm integral equation of the second kind whose kernel is
, and
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
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
and, consequently,
2.1. Shape Invariant Partially Coherent Fields: LG Modes
A particular set of orthonormal mode is formed by the LG functions, defined as
where
,
,
, and
are the associated Laguerre polynomials [
23,
24]. The parameter
(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
Assuming paraxial approximation, the free-space propagated LG modes can be calculated as [
23]
where the parameters
,
and
are
being the Rayleigh distance. From Equation (
5) the resulting propagated CSD turns out to be
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
constant, provided that
is replaced by
. This is briefly expressed by saying that the partially coherent fields radiated by our sources are shape invariant. In particular, on letting
in Equation (
11), we obtain the intensity
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])
with
and
positive parameters, into Equation (
3), with the modes given in Equation (
7), the following CSD is obtained:
where
A is a proportionality factor having dimensions of an intensity, while
are two parameters related to the widths of the intensity profile and of the degree of coherence, respectively, which are
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
) 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
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
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
and this makes such sources of the shape-invariant type. The relations between (
a,
b) and (
w,
) are still given by Equation (
15), while
which implies
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
, with equal eigenvalues was chosen. In particular, if
with
the Kronecker symbol, the following CSD is obtained:
where
is the incomplete Gamma function [
24]. The corresponding intensity reads
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.,
with
and
.
The resulting CSD is of the form
with
and
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
with
.
Using the expression of the modes in Equation (
7), the following expression for the resulting CSD turns out:
where
Note that the case 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
on the (finite or infinite) interval
, which are orthonormal with respect to the non-negative weight function
. The theorem states that, if
is the coefficient of
in the polynomial and
then
Using the above result, Equation (
30) can be written as
if
. In the limit of
, the resulting indeterminate can be evaluated by means of the de l’Hôpital theorem [
24], which gives
A source characterized by the CSD in Equation (
29), with the function
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
, and (
34) we obtain the following expression for the intensity:
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,
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,
where
is the Bessel function of the first kind and order
m [
24], the following asymptotic expression is obtained for the source intensity:
The usefulness of introducing the Bessel functions is apparent if we consider their property according to which
is negligibly small for
[
24,
49]. This allows us to roughly estimate the radius of the central dark region, say
. In fact, if we neglect the difference between
m and
, as well as that between
N and
, we can state that all the Bessel functions appearing in Equation (
37) vanishes if
so that
Radial profiles and 2D plots of the intensity profile are shown in
Figure 3 for
and
, 50, and 100, corresponding to
values of
w,
w, and
w, respectively. It is worth noting that the intensity drops practically to zero outside the donut region. For example, the intensity
is below
of its maximum for
or
w. The same happens for
if
w or
w, and for
if
w or
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
It is seen that the absolute value of
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.,
for any choice of
m,
and
. 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
are shown in
Figure 4. It can be observed that setting
(
) close to zero and varying
(
), it passes through
N zeros. The same happens for large values of
(
), i.e., at points outside the intensity donut. The positions of such minima along
(
) remain almost the same on varying
(
) and are located within the intensity donut (for
).
Maps of the degree of coherence across the source plane can be seen in
Figure 5, which shows the absolute value of
as a function of
for three different values of
. The three positions
are denoted by small dots. Chosen parameters are
and
, 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
fixed and increasing
we obtain
In particular, if
is close to zero,
The absolute value of the degree of coherence along a radius is also represented in
Figure 6 as a function of
for three different values of
. The order
N is fixed (equal to 5) and several values of
m are considered (
). The three values of
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 (
) with topological charge
. 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
the degree of coherence is over 0.75 for any pair of points located within a circle of radius
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
and
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
and
, with the same weight, i.e.,
the intensity profile does not change with respect to the previous case, but an angular modulation is introduced into the degree of coherence, being
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
. In fact, if
the degree of coherence is zero for points that satisfy
, while is unitary for points with
. 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
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 for any pairs of diametrically opposite points on a circle constant.
Sources with the above characteristics could have potential applications in several fields, such as particle trapping, optical communications, imaging, and more.