Next Article in Journal
Atmospheric Pre-Compensation of Ground-to-Space Communications with Adaptive Optics: Past, Present and Future—A Field Review
Next Article in Special Issue
Simple Method of Light Field Calculation for Shaping of 3D Light Curves
Previous Article in Journal
Vortex Beam in a Turbulent Kerr Medium for Atmospheric Communication
Previous Article in Special Issue
Multi-User Nonlinear Optical Cryptosystem Based on Polar Decomposition and Fractional Vortex Speckle Patterns
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Nonparaxial Focusing of Partially Coherent Gaussian Schell-Model and Bessel-Correlated Beams in Free Space

by
Nikolai I. Petrov
Scientific and Technological Centre of Unique Instrumentation, Russian Academy of Sciences, Moscow 117342, Russia
Photonics 2023, 10(7), 857; https://doi.org/10.3390/photonics10070857
Submission received: 6 July 2023 / Revised: 19 July 2023 / Accepted: 20 July 2023 / Published: 24 July 2023
(This article belongs to the Special Issue Light Focusing and Optical Vortices)

Abstract

:
The nonparaxial focusing of partially coherent beams in free space has been studied using the coherent-state and coherent-mode decomposition methods. Analytical expressions for the width and angular divergence of partially coherent Gaussian Schell-model (GSM) beams have been obtained using the coherent-state method. It has been shown that the focusing plane is shifted in the opposite axial direction compared to the geometric focusing plane. The influence of the nonparaxiality and spatial coherence of Bessel-correlated vortex beams on the intensity distribution and displacement of the focus plane has been analyzed. It has been shown that the shift of the focus plane increases with a decrease in the coherence radius of the source. A smaller diffraction spread has been shown for partially coherent Bessel-correlated beams compared to GSM beams.

1. Introduction

The theory of the propagation of partially coherent waves in free space, as well as in inhomogeneous media, has now been developed quite fully [1,2,3,4,5]. Partially coherent sources can be used to overcome speckles in imaging systems and increase the signal-to-noise ratio in free-space optical communication, remote sensing, etc. Gaussian Schell-model (GSM) beams are familiar partially coherent sources which represent natural sources [1]. The properties of these beams were studied extensively until recently [6,7,8,9,10,11,12,13]. Partially coherent sources with specific propagation properties that lead to the formation of highly directional light beams were considered in the last two decades [14,15,16,17,18,19,20,21,22,23].
Tightly focused beams find important applications in optical data storage, microscopy, lithography, particle trapping, etc. The focusing of partially coherent GSM and Laguerre–Gaussian-correlated Schell-model beams was studied theoretically and experimentally [24,25,26,27,28,29,30]. The focusing properties of partially coherent beams with specific correlation functions were also investigated. In [31], the focusing of partially coherent Bessel–Gaussian beams was considered.
The focusing of a Gaussian beam in free space has been long and well studied using various methods, including analytical, asymptotic, and numerical approaches to integration. Usually, the problem is solved by calculation of diffraction integrals. However, the numerical simulation of the diffraction integrals is often time-consuming, and asymptotic methods have been developed for the calculations. A hybrid integration method combining numerical integration with asymptotic methods was proposed in [32]. In [33], the accuracy and computational savings of this hybrid technology were examined. Because of the existence of high-frequency oscillatory terms in the Rayleigh–Sommerfeld integral, the calculations are usually performed numerically on a submicron scale. Therefore, other time-saving modeling approaches are needed to facilitate the solution of the problem. Coherent-mode decomposition is a powerful method to simulate the propagation of partially coherent beams [34]. Another effective analytical method for studying the propagation and focusing of light beams is the quantum theoretical method of coherent states [35,36,37]. Density matrix formalism is useful for describing partially coherent sources in inhomogeneous media [38]. Most publications are devoted to the analysis of beam characteristics in paraxial approximation. However, it is well known that nonparaxial effects significantly influence the characteristics of a tightly focused beam [23,39,40,41,42,43,44].
While the focusing of Gaussian Schell-model beams is well studied [27,28,29,30], the nonparaxial focusing of Bessel-correlated beams has not been considered yet.
In this paper, the nonparaxial focusing and spreading of partially coherent Gaussian Schell-model and Bessel-correlated vortex beams in free space are investigated using the coherent-state and -mode decomposition methods. The influence of the effects of nonparaxiality on the field intensity distributions in the axial and radial directions near the focusing plane is analyzed.

2. Focusing of Gaussian Schell-Model Beam

The cross-spectral density (CSD) function of the Gaussian Schell-model beam is given by [1]
Γ ( x 1 , x 2 , 0 ) = I 0 exp x 1 2 + x 2 2 a 2 x 1 x 2 2 r 0 2 i α 0 2 x 2 2 x 1 2
where a is the width of the beam, r 0 is the coherence length, α 0 = k / R f , k = 2 π / λ is the wavenumber, and R f is the wavefront curvature radius.
Below we consider the coherence function (1) in the coherent-state (CS) representation [35,36,37,38,39,40]:
Γ α , β , 0 = α Γ ^ β = α x x Γ ^ x x β d x d x
The coherent state |α⟩ can be defined as an eigenfunction of the annihilation operator a ̑ :
a ̑ α = α α ,
where a ̑ = x ^ w 0 + i k w 0 2 p ^ x , p ̑ x = i k x .
The explicit form of CS is given by a Gaussian beam function
α = 2 π w 0 2 1 4 exp x 2 w 0 2 + 2 x w 0 α α 2 2 α 2 2 ,
where the complex eigenvalues α = x 0 w 0 + i k w 0 2 p x 0 = α exp ( i ϑ ) determine the initial coordinates x0 of the center of the elementary beam and the angle of its inclination px0 = n0sinθ0 to the z-axis, n0 is the refractive index of the medium, and w0 is the elementary beam width.
The coherent states are the generating functions for the modes:
Ψ α ( x ) = α = exp α 2 2 m α m m ! 1 / 2 m
where | m = Ψ m x = k ω π 1 / 4 1 2 n m ! e x p k ω 2 x 2 H m ( k ω x ) is the Gauss–Hermite function.
Note that the Gauss–Hermite functions are the modal solutions of the Helmholtz equation describing the propagation of light in a planar graded-index waveguide with the refractive index in the transverse direction n 2 x = n 0 2 ω 2 x 2 , where n 0 is the refractive index on the axis and ω is the gradient parameter.
The average of any beam value M ^ can be obtained using the following expression:
M = T r Γ M T r Γ = α Γ ^ β β M ^ α d 2 α d 2 β α Γ ^ β d 2 α d 2 β
where d 2 α = d R e α d I m α , and Tr denotes the trace.

2.1. Beam Width

The average beam width is determined by
Δ x 2 = x ^ 2 x ^ 2 = d 2 α d 2 β α Γ ^ β β x ^ 2 α / d 2 α d 2 β α Γ ^ β β α
where x ^ = 1 2 k ω a ^ + a ^ + ; x ^ 2 = 1 2 k ω 1 + 2 a ^ + a ^ + a ^ 2 + a ^ + 2 , a ^ a ^ + = a ^ + a ^ + 1 ;
β x ^ 2 ( z ) α = 1 2 k ω 1 + e 2 i Ω z e x p e 2 χ 1 β * α e χ α 2 + e 2 i Ω z e x p e 2 χ 1 β * α e χ β * 2 + β * α β α ,
χ = i ω 2 k n 0 3 z ,   β α = e x p α 2 2 β 2 2 + α β * ,
α Γ ^ β = I 0 k ω π 1 / 2 e x p α * 2 + α 2 + β 2 + β 2 2 π p ~ e x p q ˇ 2 4 p ~ π p e x p k ω α * 2 2 p ,
where p = k ω 2 + 1 a 2 + 1 r 0 2 i α 0 2 , p ~ = k ω 2 + 1 a 2 + 1 r 0 2 + i α 0 2 1 p r 0 4 ,
q ~ = 2 k ω β 2 k ω α * p r 0 2 ,   w 0 2 = 2 k ω .
Calculating the integral (6), we obtain
Δ x 2 = 1 2 k ω 1 + 2 γ γ 2 + 4 ξ η 1 γ 2 4 ξ η 2 ξ 1 γ 2 4 ξ η 1 / 2 exp i 2 ω n 0 + 3 ω 2 k n 0 3 z 1 4 ξ η exp ( 2 χ ) + γ 2 exp ( 2 χ ) 2 γ exp ( χ ) 3 / 2 2 η 1 γ 2 4 ξ η 1 / 2 exp i 2 ω n 0 + 3 ω 2 k n 0 3 z 1 4 ξ η exp ( 2 χ ) + γ 2 exp ( 2 χ ) 2 γ exp ( χ ) 3 / 2 ,
where η = 1 2 k ω 2 κ κ 2 1 / r 0 4 , ξ = η * , κ = k ω 2 + 1 a 2 + 1 r 0 2 i α 0 2 , γ = k ω r 0 2 κ 2 1 / r 0 4 ,
χ = 2 i ω 2 k n 0 3 z ,   w 0 2 = 2 k ω
Note that the third and fourth terms on the right side of Equation (9) are complex values. However, these two values are complex conjugate quantities, so the sum of these two values is a real quantity. If we put n 0 = 1 and ω 0 , then the expressions obtained will describe the focusing of the light beam in free space.
Figure 1 shows the change in beam width during focusing for different coherence lengths of the incident beam at a wavefront curvature radius R f = 100 μm. Here and below on the graphs, the wavelength λ = 0.63 μm.
It follows that the beam width in the plane of focus increases with decreasing coherence length. The beam width is 0.72 μm for r 0 = 100   μ m , 1.547 μm for r 0 = 10   μ m , and 2.704 μm for r 0 = 5   μ m . In addition, the focus plane shifts δ z = 0.5   μ m at r 0 = 100   μ m , δ z = 4   μ m at r 0 = 10   μ m , and δ z = 14   μ m at r 0 = 5   μ m , which occurs due to the effects of nonparaxiality and coherence.
In Figure 2, the spreading of a partially coherent focused beam with R f = (plane wavefront) is shown for different initial coherence lengths. In contrast to Figure 1, where beam spreading occurs due to the wavefront curvature and diffraction, here only diffraction affects the beam width.
The beam width increases faster with distance with a shorter coherence length of the incident beam. In addition, the refractive index of the medium affects the spreading of the beam (Figure 2b). The spreading of the beam weakens with an increase in the refractive index of the medium.

2.2. Angular Divergence

The angular divergence of the partially coherent Gaussian Schell-model beam (1) in the initial plane z = 0 can be represented in the form
Δ θ 2 = Δ p 2 = 1 k 2 1 a 0 2 + 2 r 0 2 + k 2 a 0 2 4 R f 2 ,
where a 0 is the beam radius, r 0 is the coherence radius, and R f is the wavefront curvature radius.
In Figure 3, the angular divergence dependence on the beam waist width is presented for different coherence lengths.
It can be seen that the angular divergence of the GSM beam is higher for a source with less coherence.
The change in the angular divergence with the distance is determined by the formula
Δ p 2 = p ^ 2 p ^ 2 = T r Γ 0 p 2 z T r Γ 0 T r 2 Γ 0 p z T r Γ 0 ,
where p ^ = i ω 2 k a ^ a ^ + ; p ^ 2 = ω 2 k 1 + 2 a ^ + a ^ a ^ 2 a ^ + 2 .
Calculating the traces in (11), we obtain
Δ p 2 = ω 2 k 1 + 2 γ γ 2 + 4 ξ η 1 γ 2 4 ξ η + 2 ξ 1 γ 2 4 ξ η 1 / 2 exp i 2 ω n 0 + 3 ω 2 k n 0 3 z 1 4 ξ η exp ( 2 χ ) + γ 2 exp ( 2 χ ) 2 γ exp ( χ ) 3 / 2 + 2 η 1 γ 2 4 ξ η 1 / 2 exp i 2 ω n 0 + 3 ω 2 k n 0 3 z 1 4 ξ η exp ( 2 χ ) + γ 2 exp ( 2 χ ) 2 γ exp ( χ ) 3 / 2

2.3. Uncertainty Relationship

The threshold for the maximum focusing of the light is determined by the following uncertainty relationship:
Δ x 2   1 4 k 2 Δ p 2 ,
where the angular divergence Δ p is given by (12).
Relationship (13) is valid for a beam with a flat wavefront. In the case of beams with a spherical wavefront, the following uncertainty relationship can be obtained [39,40]:
R f 2   Δ x 4 Δ x 2 Δ p 2 1 / 4 k 2
and because Δ p 1 , R f Δ x , i.e., the minimum value of the wavefront curvature radius is diffraction-limited.
This indicates that the numerical aperture (convergence angle) is also diffraction-limited.
Relationship (14) can be obtained from the following Schrodinger–Robertson uncertainty relation [42]:
Δ x 2 Δ p 2 σ x p 2 σ x 2 σ p 2 1 r 2 1 / 4 k 2 ,
where r = σ x p σ x σ p , R f = σ x 2 σ x p , σ x p = 1 2 x p + p x x p ,
σ x p = 1 2 x p + p x = i 2 k a ^ 2 a ^ + 2 ,   x p = 0 ,
σ x p = i 2 k 2 ξ 1 γ 2 4 ξ η 1 / 2 exp i 2 ω n 0 + 3 ω 2 k n 0 3 z 1 4 ξ η exp ( 2 χ ) + γ 2 exp ( 2 χ ) 2 γ exp ( χ ) 3 / 2 2 η 1 γ 2 4 ξ η 1 / 2 exp i 2 ω n 0 + 3 ω 2 k n 0 3 z 1 4 ξ η exp ( 2 χ ) + γ 2 exp ( 2 χ ) 2 γ exp ( χ ) 3 / 2
In Figure 4, the dependence of the uncertainty values on the propagation distance is shown.
The uncertainty changes with distance, taking a minimum value in the plane of focus.

3. Focusing of Bessel-Correlated Beams

Consider now a partially coherent Bessel-correlated vortex beam with the wavefront curvature radius Rf in the plane z = 0 [23]:
Γ r 1 , r 2 = A ξ l / 2 ( 1 ξ ) e x p 1 + ξ 1 ξ r 1 2 + r 2 2 2 w 2 i k 2 R f ( r 2 2 r 1 2 ) I l 4 ξ 1 ξ r 1 r 2 w 2 e x p i l ( φ 1 φ 2 ) ,
where ξ = r 0 4 w 4 1 + w 4 r 0 4 1 / 2 1 2 , w is the spot size at the waist of the beam, r 0 is the coherence length, R f is the wavefront curvature, and l is the topological charge.
Note that when r 0 (fully coherent case) ξ 0 , and when r 0 0 (completely incoherent case) ξ 1 .
Unlike the CSD proposed in [15], here we have introduced a new parameter—the wavefront curvature radius Rf. In contrast to the conventional LG modes, here we consider the generalized Laguerre–Gauss modes with spherical wavefronts as the eigenfunctions
Ψ n l r , φ = 1 w 2 n ! π ( n + l ) ! 1 / 2 2 r w l e r 2 / w 2 i k r 2 / 2 R f L n l 2 r 2 w 2 e i l φ .
The CSD (17) can be represented in the form of coherent mode decomposition as [23]
Γ ( r 1 , r 2 , 0 ) = n λ n Ψ n * r 1 Ψ n r 2 = n λ n p p a p * a p ψ p * ψ p
where λ n = n ! ( n + l ) ! ξ n ; ξ = r 0 4 w 4 1 + w 4 r 0 4 1 / 2 1 2 ; ψ p l r , φ = J l ( μ p l r R 0 ) exp ( i l φ ) π R 0 J l + 1 ( μ p l ) ;
a p = B 0 ( 1 ) n μ p l 2 R 0 l q * n q n + l + 1 e x p μ p l 2 4 R 0 2 q L n l μ p l 2 2 w 2 R 0 2 q 2 ;   B 0 = 2 w 2 ( l + 1 ) / 2 n ! ( n + 1 ) ! 1 R 0 J l + 1 ( μ p l ) ;
q = 1 w 2 + i k 2 R f ; q * = 1 w 2 i k 2 R f ; and μ p l are the positive zeros of the Bessel function J l ( z ) .
The evolution of the modified CSD (17) can be presented as
Γ r 1 , r 2 , z = n l λ n l p p c p l c p l J l ( μ p r 1 / R 0 ) J l ( μ p r 2 / R 0 ) e i ( β p l β p l ) z ,
where c p l = a p π R 0 J l + 1 ( μ p l ) , β p l = k 1 μ p l k R 0 2 1 / 2 .

3.1. Average Beam Radius

The average value of the beam radius can be expressed as the mean square deviation Δ r :
Δ r 2 = r ^ 2 r ^ 2 = T r Γ 0 r 2 z T r Γ 0 T r 2 Γ 0 r z T r Γ 0
Due to the axial symmetry of the incident beam, the average value r ^ 2 = 0 .
The average value r ^ 2 is determined by the expression
r 2 ( z ) = T r Γ r 2 T r Γ = p p ψ p Γ ^ ψ p ψ p r 2 ψ p ,
where Γ ^ = n λ n p a p | J l ( μ p r / R 0 ) e i β p z p a p * J l ( μ p r / R 0 ) | e i β p z ,
ψ p r 2 ψ p = 0 R 0 r 2 J 0 ( μ p r / R 0 ) J 0 ( μ p r / R 0 ) r d r = 4 μ p μ p Δ 2 J 1 ( μ p ) J 1 ( μ p ) ,   if   l = 0 , ψ p r 2 ψ p = 0 R 0 r 2 J 1 ( μ p r / R 0 ) J 1 ( μ p r / R 0 ) r d r = 2 R 0 2 Δ 2 σ J 1 ( μ p ) J 1 ( μ p ) ,   if   l = 1 ,
Δ = μ p R 0 2 μ p R 0 2 , σ = μ p R 0 2 + μ p R 0 2 . Analytical expressions can also be obtained for l > 1 . Note that expression (22) describes the nonparaxial radius of the beam when the propagation constants are given by β p l = k 1 μ p l k R 0 2 1 / 2 . Expression (22) describes the paraxial radius of the beam when the propagation constants are given by β p l = k 1 1 2 μ p l k R 0 2 .
In Figure 5, the dependences of the nonparaxial (1) and paraxial (2) beam radii w = Δ r on the propagation distance are presented.
It can be seen that the focusing plane corresponding to the minimum value of the beam radius is shifted in the opposite axial direction compared to the geometric focusing plane z f g = R f , i.e., the distance to the focusing plane decreases with decreasing coherence radius. The focus plane shift and the average beam radius increase with a decrease in the radius of coherence. Note that the focusing plane in the paraxial approximation coincides with the geometric focusing plane.

3.2. Diffraction Spreading

In Figure 6, the diffraction spreading values of the beam radii of the focused plane wavefront GSM and the Bessel-correlated beams with distance are shown.
It can be seen that a partially coherent Bessel-correlated beam has less beam spread than the Gaussian Schell-model beam. This advantage of Bessel-correlated beams can be useful for free-space optical communication.

3.3. Beam Intensity Profile

The intensity distribution is determined with the following expression (20):
I r , z = Γ r , r , z = n l λ n l Ψ n l * ( r , z ) Ψ n l ( r , z ) = n λ n p p a p * a p ψ p * ψ p e i ( β p l β p l ) z
This expression describes the effects of nonparaxiality when the propagation constants are given by β p l = k 1 μ p l k R 0 2 1 / 2 and it does not take into account nonparaxial effects when the propagation constants are given by β p l = k 1 1 2 μ p l k R 0 2 .
In Figure 7, the intensity distributions along the propagation axis I 0 , z = Γ ( 0,0 , z ) are presented for nonparaxial and paraxial beams with different coherent radii and R f = 100   μ m , w = 30   μ m . The simulations show that intensity oscillations in front of the focus plane occur for coherent radiation (Figure 7a). However, these oscillations disappear if the coherence radius becomes smaller than the incident beam radius (Figure 7b). Note that there are no oscillations in the paraxial approximation. This indicates that nonparaxiality causes oscillations in the intensity distribution. However, the oscillations disappear for sources with low coherence (Figure 7b, curve 1). The focus plane shift increases, and the intensity decreases with a decrease in the radius of coherence.
In Figure 8, the intensity distributions in the initial plane and the focus plane are presented for a Bessel-correlated beam with different coherence radii and topological charge l = 0.
In the initial plane z = 0 , the intensity distribution becomes Gaussian when the radius of coherence r 0 , i.e., the source is completely coherent. The intensity distribution at r 0 0 (a completely incoherent source) becomes asymmetric with a long tail whose length increases with decreasing coherence radius.
In Figure 9, the intensity distributions in the focal plane are presented for a Bessel-correlated vortex beam with l = 1.
An increase in the coherence radius leads to an annular distribution with sharp boundaries (Figure 9a,b). The thickness of the ring is equal to FWHM = 0.86 μm in the focus plane (Figure 9b). Long tails in the intensity distribution appear with a decrease in the coherence radius (Figure 9c,d). In this case, the value FWHM = 0.78 μm in the focus plane. Although the bright central part of the ring is smaller in the case of a low coherence beam, a significant amount of the power of the focused beam is outside this bright ring.

4. Conclusions

Thus, the nonparaxial focusing of partially coherent Gaussian Schell-model and Bessel-correlated vortex beams in free space has been investigated using the coherent-state and the mode decomposition methods.
It is found that the average radii in the focusing plane increase with decreasing coherence radius. However, calculations of the intensity profiles show that bright spots and rings of smaller sizes can be obtained using the sources with low-coherence Bessel-correlated sources. The beam width estimated on the basis of the intensity distribution, i.e., the FWHM (full width at half maximum) value, shows a behavior different from that based on the calculation of the standard deviation. Indeed, the FWHM value of the focusing spot of a low-coherent light beam determined by the intensity distribution can be smaller than that of a coherent light beam. It is shown that the effect of shifting the focusing plane in the axial direction opposite to the geometric focusing plane is enhanced with a decrease in coherence.
It follows from the simulation results that the partially coherent Bessel-correlated beams have less spread than the Gaussian Schell-model beams. This preferred property can be useful in specific fields, especially in optical communication in free space.
Future research can be related to the consideration of focusing partially coherent and partially polarized vector beams [45]. Of practical interest is the consideration of structured vortex beams with spin angular momentum (SAM) and orbital angular momentum (OAM) [46,47,48], partially coherent accelerating beams [49], beams possessing self-healing abilities in a turbulent media [50], and spatiotemporal optical vortex beams [51].
In summary, the nonparaxial focusing of the GSM beam has been studied using the CS decomposition method. The nonparaxial focusing and spreading of a Bessel-correlated beam in free space were analyzed using the coherent-mode decomposition method. It is shown that partially coherent Bessel-correlated beams have a smaller spatial spread than GSM beams. The average spot size in the focusing plane increases with a decrease in the degree of coherence. However, the FWHM value of the focusing spot of a low-coherence light beam, determined by the intensity distribution, can be less than that of a coherent light beam. It is shown that the shift of the focusing plane increases with decreasing coherence.
The results obtained can be useful for microparticle trapping, transmitting information, optical visualization, and optical communication in free space.

Funding

This research was funded by the Ministry of Science and Higher Education of the Russian Federation under the State contract No. FFNS-2022-0009.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Wolf, E. Introduction to the Theory of Coherence and Polarization of Light; Cambridge University Press: Cambridge, UK, 2007. [Google Scholar]
  2. Mandel, L.; Wolf, E. Optical Coherence and Quantum Optics; Cambridge University Press: Cambridge, UK, 1995. [Google Scholar]
  3. Wolf, E. Unified theory of coherence and polarization of random electromagnetic beams. Phys. Lett. A 2003, 312, 263–267. [Google Scholar] [CrossRef]
  4. Tervo, J.; Setala, T.; Friberg, A.T. Theory of partially coherent electromagnetic fields in the space-frequency domain. J. Opt. Soc. Am. A 2004, 21, 2205–2215. [Google Scholar] [CrossRef]
  5. Korotkova, O.; Hoover, B.G.; Gamiz, V.L.; Wolf, E. Coherence and polarization properties of far fields generated by quasi-homogeneous planar electromagnetic sources. J. Opt. Soc. Am. A 2005, 22, 2547–2556. [Google Scholar] [CrossRef]
  6. Liu, Z.; Zeng, J.; Cai, Y. Review on vortex beams with low spatial coherence. Adv. Phys. 2019, 4, 1626766. [Google Scholar] [CrossRef] [Green Version]
  7. 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]
  8. Dong, M.; Zhao, C.; Cai, Y.; Yang, Y. Partially coherent vortex beams: Fundamentals and applications. Sci. China Phys. Mech. Astron. 2021, 64, 224201. [Google Scholar] [CrossRef]
  9. Zhang, H.; Lu, X.; Wang, Z.; Konijnenberg, A.P.; Wang, H.; Zhao, C.; Cai, Y. Generation and Propagation of Partially Coherent Power-Exponent-Phase Vortex Beam. Front. Phys. 2021, 9, 781688. [Google Scholar] [CrossRef]
  10. Chen, Y.; Wang, F.; Cai, Y. Partially coherent light beam shaping via complex spatial coherence structure engineering. Adv. Phys. X 2022, 7, 2009742. [Google Scholar] [CrossRef]
  11. Krivoshlykov, S.G.; Petrov, N.I.; Sisakyan, I.N. Spatial coherence of optical fields in longitudinally inhomogeneous media with the square-law index profile. Sov. J. Quantum Electron. 1985, 15, 330–338. [Google Scholar] [CrossRef]
  12. Petrov, N.I.; Petrova, G.N. Diffraction of partially coherent light beams by micro-lens arrays. Opt. Express 2017, 25, 22545–22564. [Google Scholar] [CrossRef] [PubMed]
  13. Petrov, N.I. Holographic diffuser with controlled scattering indicatrix. Comp. Opt. 2017, 41, 831–836. [Google Scholar] [CrossRef]
  14. Gori, F.; Palma, C. Partially coherent sources which give rise to highly directional light beams. Opt. Commun. 1978, 27, 185–188. [Google Scholar] [CrossRef]
  15. Ponomarenko, S.A. A class of partially coherent beams carrying optical vortices. J. Opt. Soc. Am. A 2001, 18, 150–156. [Google Scholar] [CrossRef]
  16. Marasinghe, M.L.; Premaratne, M.; Paganin, D.M.; Alonso, M.A. Coherence vortices in Mie scattered nonparaxial partially coherent beams. Opt. Expr. 2012, 20, 2858–2875. [Google Scholar] [CrossRef]
  17. Alieva, T.; Rodrigo, J.A.; Camara, A.; Abramochkin, E. Partially coherent stable and spiral beams. J. Opt. Soc. Am. A 2013, 30, 2237–2243. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Ostrovsky, A.S.; García-García, J.; Rickenstorff-Parrao, C.; Olvera-Santamaria, M.A. Partially coherent diffraction-free vortex beams with a Bessel-mode structure. Opt. Lett. 2017, 42, 5182–5185. [Google Scholar] [CrossRef]
  19. Santamaría, M.A.O.; Garciapiña, J.L.R.; García, J.G.; Ostrovsky, A.S. Generation of a partially coherent secondary source with Bessel-mode vortex structure by liquid crystal spatial light modulator. Opt. Commun. 2019, 439, 312–316. [Google Scholar] [CrossRef]
  20. 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]
  21. Xu, Z.; Li, X.; Liu, X.; Ponomarenko, S.A.; Cai, Y.; Liang, C. Vortex preserving statistical optical beams. Opt. Exp. 2020, 28, 8475–8483. [Google Scholar] [CrossRef]
  22. Ferlic, N.A.; van Iersel, M.; Davis, C.C. Simulation of Im-Bessel beam propagation through time-correlated atmospheric turbulence. Proc. SPIE 2021, 11834, 118340L. [Google Scholar]
  23. Petrov, N.I. Nonparaxial Propagation of Bessel Correlated Vortex Beams in Free Space. Micromachines 2023, 14, 38. [Google Scholar] [CrossRef]
  24. Fischer, D.G.; Visser, T.D. Spatial correlation properties of focused partially coherent light. J. Opt. Soc. Am. A 2004, 21, 2097–2102. [Google Scholar] [CrossRef]
  25. Foreman, M.R.; Torok, P. Focusing of spatially inhomogeneous partially coherent, partially polarized electromagnetic fields. J. Opt. Soc. Am. A 2009, 26, 2470–2479. [Google Scholar] [CrossRef] [Green Version]
  26. Chen, B.; Zhang, Z.; Pu, J. Tight focusing of partially coherent and circularly polarized vortex beams. J. Opt. Soc. Am A 2009, 26, 862–869. [Google Scholar] [CrossRef] [PubMed]
  27. Wang, F.; Zhu, S.; Cai, Y. Experimental study of the focusing properties of a Gaussian Schell-model vortex beam. Opt. Lett. 2011, 36, 3281–3283. [Google Scholar] [CrossRef] [PubMed]
  28. Xu, H.F.; Zhang, Z.; Qu, J.; Huang, W. The tight focusing properties of Laguerre–Gaussian-correlated Schell-model beams. J. Mod. Opt. 2016, 63, 1429–1437. [Google Scholar] [CrossRef]
  29. Ping, C.; Liang, C.H.; Wang, F.; Cai, Y.J. Radially polarized multi-Gaussian Schell-model beam and its tight focusing properties. Opt. Express 2017, 25, 32475–32490. [Google Scholar] [CrossRef]
  30. Lin, H.; Zhou, X.; Chen, Z.; Sasaki, O.; Li, Y.; Pu, J. Tight focusing properties of a circular partially coherent Gaussian beam. J. Opt. Soc. Am. A 2018, 35, 1974–1980. [Google Scholar] [CrossRef]
  31. Zhang, Z.; Pu, J.; Wang, X. Focusing of partially coherent Bessel–Gaussian beams through a high numerical-aperture objective. Opt. Lett. 2008, 33, 49–51. [Google Scholar] [CrossRef]
  32. Stamnes, J.J. Hybrid integration technique for efficient and accurate computation of diffraction integrals. J. Opt. Soc. Am. A 1989, 6, 1330–1342. [Google Scholar] [CrossRef]
  33. Eide, H.A.; Stamnes, J.J. Exact and approximate solutions for focusing of two-dimensional waves. J. Opt. Soc. Am. A 1998, 15, 1285–1319. [Google Scholar]
  34. Gori, F.; Santarsiero, M.; Simon, R.; Piquero, G.; Borghi, R.; Guattari, G. Coherent-mode decomposition of partially polarized, partially coherent sources. J. Opt. Soc. Am. A 2003, 20, 78–84. [Google Scholar] [CrossRef]
  35. Krivoshlykov, S.G.; Sisakyan, I.N. Coherent States and Light Propagation in Inhomogeneous Media. Sov. J. Quantum Electron. 1980, 10, 312–318. [Google Scholar] [CrossRef]
  36. Krivoshlykov, S.G.; Petrov, N.I.; Sisakyan, I.N. Correlated coherent states and propagation of arbitrary Gaussian beams in longitudinally homogeneous quadratic media exhibiting absorption or amplification. Sov. J. Quantum Electron. 1986, 16, 933–941. [Google Scholar] [CrossRef]
  37. Krivoshlykov, S.G.; Petrov, N.I.; Sisakyan, I.N. Mode excitation in waveguides with a quadratic refractive-index profile by partially coherent radiation sources. Sov. Tech. Phys. 1985, 30, 1031–1036. [Google Scholar]
  38. Krivoshlykov, S.G.; Petrov, N.I.; Sisakyan, I.N. Density-matrix formalism for partially coherent optical fields propagating in slightly inhomogeneous media. Opt. Quant. Electr. 1986, 18, 253–264. [Google Scholar] [CrossRef]
  39. Petrov, N.I. Nonparaxial focusing of wave beams in a graded-index medium. Rus. J. Quant. Electron. 1999, 29, 249–255. [Google Scholar] [CrossRef]
  40. Petrov, N.I. Focusing of beams into subwavelength area in an inhomogeneous medium. Opt. Exp. 2001, 9, 658–673. [Google Scholar] [CrossRef]
  41. Yuan, Y.; Du, S.; Dong, Y.; Wang, F.; Zhao, C.; Cai, Y. Nonparaxial propagation properties of a vector partially coherent dark hollow beam. J. Opt. Soc. Am. A 2013, 30, 1358–1372. [Google Scholar] [CrossRef]
  42. Petrov, N.I. Macroscopic quantum effects for classical light. Phys. Rev. A 2014, 90, 043814. [Google Scholar] [CrossRef] [Green Version]
  43. Petrov, N.I. Evanescent and propagating fields of a strongly focused beam. J. Opt. Soc. Am. A 2003, 20, 2385–2389. [Google Scholar] [CrossRef]
  44. Petrov, N.I. Reflection and transmission of strongly focused light beams at a dielectric interface. J. Mod. Opt. 2005, 52, 1545–1556. [Google Scholar] [CrossRef]
  45. Zhang, Z.; Fan, H.; Xu, H.F.; Qu, J.; Huang, W. Three-dimensional focus shaping of partially coherent circularly polarized vortex beams using a binary optic. J. Opt. 2015, 17, 065611. [Google Scholar] [CrossRef] [Green Version]
  46. Petrov, N.I. Depolarization of Light in Optical Fibers: Effects of Diffraction and Spin-Orbit Interaction. Fibers 2021, 9, 34. [Google Scholar] [CrossRef]
  47. Petrov, N.I. Depolarization of vector light beams on propagation in free space. Photonics 2022, 9, 162. [Google Scholar] [CrossRef]
  48. Li, H.; Ma, C.; Wang, J.; Tang, M.; Li, X. Spin-orbit Hall effect in the tight focusing of a radially polarized vortex beam. Opt. Exp. 2021, 29, 39419–39427. [Google Scholar] [CrossRef]
  49. Dev, V.; Reddy, A.N.K.; Ustinov, A.V.; Khonina, S.N.; Pal, V. Autofocusing and Self-Healing Properties of Aberration Laser Beams in a Turbulent Media. Phys. Rev. Appl. 2021, 16, 014061. [Google Scholar] [CrossRef]
  50. Zeng, Z.; Zhao, D. Partially coherent dual nonparaxial accelerating beams. Ann. Phys. 2020, 533, 2000438. [Google Scholar] [CrossRef]
  51. Hyde, M.W., IV; Korotkova, O.; Spencer, M.F. Partially coherent sources whose coherent modes are spatiotemporal optical vortex beams. J. Opt. 2023, 25, 035606. [Google Scholar] [CrossRef]
Figure 1. Beam width change with propagation distance for different coherence lengths of the incident beam.
Figure 1. Beam width change with propagation distance for different coherence lengths of the incident beam.
Photonics 10 00857 g001
Figure 2. Beam width change with distance for different coherence lengths of the incident beam: 1— r 0 = 5   μ m ; 2— r 0 = 10   μ m ; 3— r 0 = 100   μ m . (a) n 0 = 1.0 ; (b) n 0 = 1.5 .
Figure 2. Beam width change with distance for different coherence lengths of the incident beam: 1— r 0 = 5   μ m ; 2— r 0 = 10   μ m ; 3— r 0 = 100   μ m . (a) n 0 = 1.0 ; (b) n 0 = 1.5 .
Photonics 10 00857 g002
Figure 3. Angular divergence as a function of beam waist width: 1— r 0 = 10   μ m ; 2— r 0 = 50   μ m ; 3— r 0 = 1000   μ m . R f = .
Figure 3. Angular divergence as a function of beam waist width: 1— r 0 = 10   μ m ; 2— r 0 = 50   μ m ; 3— r 0 = 1000   μ m . R f = .
Photonics 10 00857 g003
Figure 4. The uncertainty value as function of distance for different coherence radii of the incident beam: 1— r 0 = 5   μ m ; 2— r 0 = 10   μ m ; 3— r 0 = 100   μ m . R f = 100   μ m .
Figure 4. The uncertainty value as function of distance for different coherence radii of the incident beam: 1— r 0 = 5   μ m ; 2— r 0 = 10   μ m ; 3— r 0 = 100   μ m . R f = 100   μ m .
Photonics 10 00857 g004
Figure 5. Nonparaxial (1) and paraxial (2) change in the beam radius with distance for different coherence radii of the incident beam: R f = 15 μm (a,b); R f = 100 μm (c,d).
Figure 5. Nonparaxial (1) and paraxial (2) change in the beam radius with distance for different coherence radii of the incident beam: R f = 15 μm (a,b); R f = 100 μm (c,d).
Photonics 10 00857 g005
Figure 6. Beam radius change with distance for GSM and Bessel-correlated beams: r 0 = 5   μ m . 1—GSM beam; 2—Bessel-correlated beam.
Figure 6. Beam radius change with distance for GSM and Bessel-correlated beams: r 0 = 5   μ m . 1—GSM beam; 2—Bessel-correlated beam.
Photonics 10 00857 g006
Figure 7. Intensity distributions along the propagation axis: r 0 = 100   μ m (a); r 0 = 10   μ m (b). 1—nonparaxial; 2—paraxial.
Figure 7. Intensity distributions along the propagation axis: r 0 = 100   μ m (a); r 0 = 10   μ m (b). 1—nonparaxial; 2—paraxial.
Photonics 10 00857 g007
Figure 8. Intensity distributions of Bessel-correlated beams with w = 30 μm, R f = 100   μ m , and l = 0 in initial (a,c) and focus (b,d) planes; r 0 = 100   μ m (a,b); r 0 = 10   μ m (c,d); z = 95.35 μm (b); z = 87 μm (d). The numbers below the images of the transverse intensity distribution indicate the position of the vertical red lines.
Figure 8. Intensity distributions of Bessel-correlated beams with w = 30 μm, R f = 100   μ m , and l = 0 in initial (a,c) and focus (b,d) planes; r 0 = 100   μ m (a,b); r 0 = 10   μ m (c,d); z = 95.35 μm (b); z = 87 μm (d). The numbers below the images of the transverse intensity distribution indicate the position of the vertical red lines.
Photonics 10 00857 g008aPhotonics 10 00857 g008b
Figure 9. Intensity distributions of Bessel-correlated vortex beams with w = 20 μm, R f = 100   μ m , and l = 1 in initial (a,c) and focus (b,d) planes; r 0 = 100   μ m (a,b); r 0 = 10   μ m (c,d); z = 95.35 μm (b); z = 91.8 μm (d). The numbers below the images of the transverse intensity distribution indicate the position of the vertical red lines.
Figure 9. Intensity distributions of Bessel-correlated vortex beams with w = 20 μm, R f = 100   μ m , and l = 1 in initial (a,c) and focus (b,d) planes; r 0 = 100   μ m (a,b); r 0 = 10   μ m (c,d); z = 95.35 μm (b); z = 91.8 μm (d). The numbers below the images of the transverse intensity distribution indicate the position of the vertical red lines.
Photonics 10 00857 g009aPhotonics 10 00857 g009b
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Petrov, N.I. Nonparaxial Focusing of Partially Coherent Gaussian Schell-Model and Bessel-Correlated Beams in Free Space. Photonics 2023, 10, 857. https://doi.org/10.3390/photonics10070857

AMA Style

Petrov NI. Nonparaxial Focusing of Partially Coherent Gaussian Schell-Model and Bessel-Correlated Beams in Free Space. Photonics. 2023; 10(7):857. https://doi.org/10.3390/photonics10070857

Chicago/Turabian Style

Petrov, Nikolai I. 2023. "Nonparaxial Focusing of Partially Coherent Gaussian Schell-Model and Bessel-Correlated Beams in Free Space" Photonics 10, no. 7: 857. https://doi.org/10.3390/photonics10070857

APA Style

Petrov, N. I. (2023). Nonparaxial Focusing of Partially Coherent Gaussian Schell-Model and Bessel-Correlated Beams in Free Space. Photonics, 10(7), 857. https://doi.org/10.3390/photonics10070857

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