Next Article in Journal
Glycerol–Water Solution-Assisted Mach–Zehnder Temperature Sensor in Specialty Fiber with Two Cores and One Channel
Next Article in Special Issue
A New Type of Shape-Invariant Beams with Structured Coherence: Laguerre-Christoffel-Darboux Beams
Previous Article in Journal
SOA-Based Multilevel Polarization Shift On–Off Keying Transmission for Free-Space Optical Communication
Previous Article in Special Issue
Progress on Studies of Beams Carrying Twist
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evolution of Spatiotemporal Intensity of Partially Coherent Pulsed Beams with Spatial Cosine-Gaussian and Temporal Laguerre–Gaussian Correlations in Still, Pure Water

1
Henan Key Laboratory of Electromagnetic Transformation and Detection, College of Physics and Electronic Information, Luoyang Normal University, Luoyang 471934, China
2
Department of Physics, University of Miami, Coral Gables, FL 33146, USA
3
B. I. Stepanov Institute of Physics, National Academy of Sciences of Belarus (NASB), 220072 Minsk, Belarus
*
Author to whom correspondence should be addressed.
Photonics 2021, 8(4), 102; https://doi.org/10.3390/photonics8040102
Submission received: 17 February 2021 / Revised: 24 March 2021 / Accepted: 31 March 2021 / Published: 2 April 2021
(This article belongs to the Special Issue Structured Light Coherence)

Abstract

:
A new family of partially coherent pulsed beams with spatial cosine-Gaussian and temporal Laguerre–Gaussian correlations, named spatial cosine-Gaussian and temporal Laguerre–Gaussian correlated Schell-model (SCTLGSM) pulsed beams, is introduced. An analytic propagation formula is derived for the SCTLGSM pulsed beam through the spatiotemporal ABCD optical system characterizing a continuous dispersive medium. As an example, the evolution of spatiotemporal intensity of the SCTLGSM pulsed beam in a still, pure water column is then investigated. It is found that the SCTLGSM pulsed beams simultaneously exhibit spatiotemporal self-splitting and self-focusing phenomena, which can be attributed to the special spatial/temporal coherence structures and the presence of pulse chirper in the source plane. The physical interpretation of the obtained phenomena is given. The results obtained in this paper will be of interest in underwater optical technologies, e.g., directed energy and communications.

1. Introduction

In the past few years, a great deal of attention has been paid to partially coherent beams radiated by sources with non-Gaussian correlation functions [1,2,3,4,5,6,7,8,9]. This is due to a variety of extraordinary features that these beams exhibit on free-space propagation, such as self-splitting [10,11,12], self-focusing [13], self-steering [14], and useful shapes they acquire in the far field: flat-tops [15], rings [7,16], rectangles [17], grids [18], and dark-hollow intensity profiles [19]. These new properties have already stimulated new applications in optical particle trapping, free-space optical communications [20], classic imaging [21] and remote sensing [22]. It was also recently revealed that partially coherent sources with complex coherence states (non-trivial phase distributions) can be used for fine directional control of the produced far-field spectral density and polarization distributions [23,24].
All of the aforementioned investigations are confined to statistically stationary light sources/beams in the spatial domain. More recently, non-stationary optical pulses, i.e., those with partial temporal/spectral and spatial coherence states were introduced [25,26,27,28]. Such pulses are envisioned to benefit applications in pulse shaping, laser micromachining, medical diagnosis, ghost imaging, etc. Conventionally, in time domain, a partially coherent light pulse was characterized by Schell model with a temporal Gaussian correlation function. In 2013, Lajunen et al. introduced the non-uniformly partially coherent pulse sources, i.e., sources with a non-Gaussian correlation function [29]. The generated pulsed beams were shown to lead to self-focusing upon propagation in time domain. Subsequently, other models for partially coherent pulsed sources were developed [30,31,32,33,34] to illustrate temporal self-splitting, multispot self-focusing and formation of flat far-fields, among other effects. Partially coherent pulses with the structured phase of the correlation function were also recently studied to demonstrate the control in the average pulse arrival time [35]. However, investigations in [27,28,29,30,31,32,33,34,35,36] were solely confined to the time domain while spatial pulse profiles were assumed to be constant. The non-stationary pulses can be realized experimentally. A stationary, spatially partially coherent beam with a prescribed profile of the coherence function can be first realized with the help of a spatial light modulator or Fourier optics 2f system [12]. Then the spatial correlation can be converted to temporal correlation with the direct space-to-time pulse shaper [37].
Thus, an interesting question arises: what new phenomena can one observe on propagation of partially coherent pulsed beams radiated from sources having non-Gaussian, and in general different, profiles in temporal and spatial correlation functions at the same time? In fact, partially coherent pulses with both spatial and temporal/spectral structured correlations were already considered [38,39] but the analysis relied on Gaussian distributions alone. It would be interesting to consider non-Gaussian but analytically tractable models for the spatial and temporal partial coherence.
In this paper, we choose a combination of the cosine-Gaussian Schell model for the spatial correlations and the Laguerre–Gaussian correlated Schell-model for the temporal ones. Each of these two types of partial coherence was studied separately in the literature. The cosine-Gaussian Schell model for the spatial correlations was studied both theoretically [10] and experimentally [12] and it was demonstrated that the initially Gaussian intensity distribution changes to the dark-hollow distribution and remains shape-invariant in the far field. It was shown also that this model is much simpler for analytical treatment compared to other non-Gaussian models. The Laguerre–Gaussian correlated Schell-model for the temporal correlations was studied theoretically [31] and it was found that it exhibits temporal self-splitting into n + 1 lobes upon propagation in a dispersive medium, where n is the order of the Laguerre polynomial. A potentially high number of lobes distinguishes this model from the alternative ones, which may be attractive for some applications. Thus, we introduce a new family of spatially and temporally partially coherent pulsed beams with spatial cosine-Gaussian and temporal Laguerre–Gaussian correlations, named SCTLGSM pulses. The mathematical consistency and physical feasibility of such beams are briefly discussed in the Appendix A. We then explore the evolution of the spatiotemporal intensity of a typical SCTLGSM pulsed beam through an ABCD optical system describing a diffractive and a dispersive medium, generalizing the particular media considered in Refs. [10,12,31]. Using the second-order optical coherence theory of the non-stationary optical fields, we derive, in Section 2, the analytical expressions for the spatiotemporal intensity of the SCTLGSM pulsed beams passing through the ABCD optical systems. In Section 3, the detailed numerical results regarding the SCTLGSM pulse evolution are presented and the corresponding physical interpretations are given. In Section 4, the main results obtained in this paper are summarized. The self-splitting/self-focusing of a non-stationary pulse in the temporal and spatial domains can be potentially used in optical communication systems operating in dispersive environments, e.g., underwater, for information encoding. This phenomenon can also be applied in temporal ghost imaging [40,41] and optical coherence tomography [42] for obtaining the pulse replicas that can be used for further manipulation of imaged/sensed fields.

2. Theory and Method

In the space–time domain, the second-order correlation properties of a partially coherent pulsed beam are described by the mutual coherence function (MCF) [39,43,44]. Let for the SCTLGSM pulsed beam, the mutual coherence function at the source plane z = 0 be defined as
Γ 0 ( ρ 1 , τ 1 ; ρ 2 , τ 2 ) = S ( ρ 1 , ρ 2 ) T ( τ 1 , τ 2 )
where S and T characterize the spatial and the temporal parts, respectively:
S ( ρ 1 ,   ρ 2 ) = exp ( ρ 1 2 + ρ 2 2 4 w 0 2 ) exp [ ( x 2 x 1 ) 2 2 δ 2 ] cos [ m 2 π ( x 2 x 1 ) δ ] × exp [ ( y 2 y 1 ) 2 2 δ 2 ] cos [ m 2 π ( y 2 y ) δ ] ,
T ( τ 1 , τ 2 ) = exp ( τ 1 2 + τ 2 2 T 0 2 ) exp [ ( τ 2 τ 1 ) 2 2 T c 2 ] L n [ ( τ 2 τ 1 ) 2 T c 2 ] exp [ i ω 0 ( τ 2 τ 1 ) ]
Here ρ 1 = ( x 1 , y 1 ) and ρ 2 = ( x 2 , y 2 ) are the two-dimensional transverse position vectors, τ1 and τ2 are two time instants, cos(•) is the cosine function, and m is the beam order parameter. Further, Ln(•) denotes the Laguerre polynomial of mode order n, ω0 is the pulse’s carrier frequency, while w0 and δ denote the r.m.s. spatial beam width and the r.m.s. beam transverse coherence width, respectively. T0 and Tc represent the pulse duration and the temporal coherence length, respectively. The pulsed beam expressed by Equation (1) is therefore the spatial cosine-Gaussian and temporal Laguerre–Gaussian correlated Schell-model (SCTLGSM) pulsed beam. The SCTLGSM pulsed beam reduces to the conventional Gaussian Schell-model (GSM) pulsed beam when m = 0 and n = 0. For m ≠ 0 and n ≠ 0, the spatial and the temporal coherence parts of the mutual coherence function are modulated by the cosine function and the Laguerre function, respectively. The details of the model can be found in the Appendix A.
We will now consider propagation of the SCTLGSM pulsed beam in the spatiotemporal domain by means of the extended Collins formula, under the paraxial approximation. We consider an additional linear chirp s imposed on the pulse before it enters the dispersive medium [45]. This chirp is described by the factor exp ( i τ 2 s / 2 ) in the field amplitude, can be realized by an amplitude modulator, and results in pulse compression upon propagation in a dispersive medium in the case of Gaussian Schell model [26]. A similar effect is expected for the SCTLGSM pulsed beam. Adopting the ABCD characterization of the optical system embedded in a dispersive medium in the absence of spatiotemporal coupling, we employ the following integral formula [46,47]
Γ ( r 1 , r 2 , t 1 , t 2 , z ) = [ k ( ω 0 ) 2 π B S ] 2 ω 0 2 π B T Γ 0 ( ρ 1 , τ 1 ; ρ 2 , τ 2 ) × exp { i k ( ω 0 ) 2 B S [ A S ( ρ 1 2 ρ 2 2 ) 2 ( ρ 1 r 1 ρ 2 r 2 ) + D S ( r 1 2 r 2 2 ) ] } d 2 ρ 1 d 2 ρ 2 × exp { i ω 0 2 B T [ A T ( τ 1 2 τ 2 2 ) 2 ( τ 1 t 1 τ 2 t 2 ) + D T ( t 1 2 t 2 2 ) ] } d τ 1 d τ 2 .
here k(ω0) = n(ω0)ω0/c denotes the wave number, in which n(ω0) is the refractive index of the medium at carrier frequency ω0, c is the speed of light in vacuum, r1 = (x1, y1) and r2 = (x2, y2) are position vectors in any plane z > 0 and t1 and t2 are two time instants of the pulse profile. Further, AS, BS and DS are the transfer matrix elements of the optical system in the spatial domain [48], while AT, BT and DT are those in the temporal domain [45].
On substituting from Equations (1)–(3) into Equation (4), one can derive the following simplified expression
Γ ( r 1 , r 2 , t 1 , t 2 , z ) = H S ( r 1 , r 2 , z ) H T ( t 1 , t 2 , z ) ,
where
H S ( r 1 , r 2 , z ) = [ k ( ω 0 ) 2 π B S ] 2 S ( ρ 1 , ρ 2 ) × exp { i k ( ω 0 ) 2 B S [ A S ( ρ 1 2 ρ 2 2 ) 2 ( ρ 1 r 1 ρ 2 r 2 ) + D S ( r 1 2 r 2 2 ) ] } d 2 ρ 1 d 2 ρ 2 ,
H T ( t 1 , t 2 , z ) = ω 0 2 π B T T ( τ 1 , τ 2 ) exp { i ω 0 2 B T [ A T ( τ 1 2 τ 2 2 ) 2 ( τ 1 t 1 τ 2 t 2 ) + D T ( t 1 2 t 2 2 ) ] } d τ 1 d τ 2 .
After tedious calculation, one can obtain the following expression, for the spatial domain:
H S ( r 1 , r 2 , z ) = H S ( x 1 , x 2 , z ) H S ( y 1 , y 2 , z )
where
H S ( α 1 , α 2 , z ) = 1 Δ exp ( 2 π m 2 B S 2 2 k 2 δ 2 w 0 2 Δ 2 + i k α s α d R ) × exp [ 1 2 Δ 2 ( α s 2 w 0 2 + α d 2 Q 2 ) ] cosh [ 2 π m δ Δ 2 ( B S α s k w 0 2 i A s α d ) ] ,    
and
α s = ( α 1 + α 2 ) / 2 , α d = α 2 α 1 , ( α = x , y ) ,
Q = ( 1 4 w 0 2 + 1 δ 2 ) 1 2 ,
Δ = A S 2 + ( B S k w 0 Q ) 2 ,
R = B S ( B S k w 0 Q ) 2 + A S 2 B S D S ( B S k w 0 Q ) 2 + A S ( A S D S 1 ) ,
with the spatial ABCD matrix in the form
( A S B S C S D S ) = ( 1 z 0 1 ) .
On the other hand, in the temporal domain, one can obtain the expression
H T ( t 1 , t 2 , z ) = ω 0 T 0 2 2 B T exp ( i ω 0 D T B T t ¯ t d ) exp ( T 0 2 ω 0 2 8 B T 2 t d 2 ) exp ( b 2 4 a ) q = 0 n ( n q ) H 2 q ( b / 2 a ) q ! ( 2 T c ) 2 q a q + 1 / 2 ,
where H2q is Hermite polynomial of order 2q, and
t ¯ = ( t 1 + t 2 ) / 2 ,   t d = t 2 t 1 ,
a = 1 2 T 0 2 + 1 2 T c 2 + A T 2 T 0 2 ω 0 2 8 B T 2 ,
b = ω 0 B T t ¯ + i A T T 0 2 ω 0 2 4 B T 2 t d ,
with the temporal ABCD matrix given by expression
( A T B T C T D T ) = ( 1 + s β 2 z ω 0 β 2 z s / ω 0 1 ) .
Here we assume that the refractive index of the dispersive medium is given by expression n(ω) = naω + nb, where na = β2c and nb = c/vg − 2β2ω0c. Here β2 denotes the group velocity dispersion and vg is the group velocity of the pulse [49]. Furthermore, the time coordinate is a retarded time with respect to a frame moving with group velocity vg and s is the chirp coefficient of the pulse.
Let r1 = r2 = r and t1 = t2 = t in Equations (8) and (15), respectively. Then the spatiotemporal intensity of the SCTLGSM pulsed beam becomes
I ( r , t , z ) = Γ ( r , r , t , t , z ) = H S ( r , r , z ) H T ( t , t , z ) = 1 Δ 2 exp ( 2 π m 2 B S 2 k 2 δ 2 w 0 2 Δ 2 ) exp [ 1 2 Δ 2 ( x 2 + y 2 w 0 2 ) ] × cosh [ 2 π m δ Δ 2 ( B S x k w 0 2 ) ] cosh [ 2 π m δ Δ 2 ( B S y k w 0 2 ) ] × ω 0 T 0 2 2 B T exp ( ω 0 2 4 a B T 2 t 2 ) q = 0 n ( n q ) H 2 q ( ω 0 t / 2 a B T ) q ! ( 2 T c ) 2 q a q + 1 / 2 ,
where cosh(•) is the hyperbolic cosine function. Equation (20) is the main formula derived in this paper. It can be used to investigate the spatiotemporal intensity evolution of the SCTLGSM pulsed beams in any diffractive and dispersive medium.

3. Spatiotemporal Intensity Evolution of the SCTLGSM Pulsed Beams in Water

In this section, the interaction of the SCTLGSM pulsed beams in a column of still, pure water is examined by means of numerical examples. The group velocity refractive index of water at 20 °C and the standard atmospheric pressure is ng = c/vg = 1.3591 and its group velocity dispersion coefficient is β2 = 58.174 ps2km−1 [50]. In the following calculation, the pulses’ and medium’s parameters are chosen to be w0 = 2 mm, σ = 2 mm, T0 = 4 ps, Tc = 2 ps, ω0 = 3.667 rad/fs (λ0 = 514 nm [51]), unless different values are specified. Because of still water absorption [52] the propagation distances will be limited to 100 m throughout the text.
Figure 1a–d shows the evolution of the spatiotemporal intensity I(x, y, t, z) of the SCTLGSM pulsed beam as a function of propagation distance z and horizontal coordinate x with different beam orders, m = 0, 1, 2 and 3. The other values of parameters are y = 0, t = 0, s =0 and n = 0. Figure 2a,b shows the same as Figure 1 but at different time instants t = 2 ps, 5 ps, and for m = 1. As can be seen, beam order plays an important role in the spatial splitting of the SCTLGSM pulsed beam. Of course, for the GSM pulsed beam (m = 0), there is no beam splitting upon propagation. However, for m = 1, the beam splits into two sub-beams when the propagation distance z exceeds 30 m. For larger beam orders (see m = 2 and m = 3), the starting point of the split moves closer to the source plane z = 0. Note that time t has no impact on the beam splitting position. For larger time instants t, the spatiotemporal intensity is gradually reduced (see t = 2 ps and 5 ps in Figure 2). These results can be interpreted physically as follows. Equation (20) can be expressed as
I ( r , t , z ) = 1 4 Δ 2 { exp [ 1 2 Δ 2 w 0 2 ( x 2 π m B S δ k ) 2 ] + exp [ 1 2 Δ 2 w 0 2 ( x + 2 π m B S δ k ) 2 ] } × { exp [ 1 2 Δ 2 w 0 2 ( y 2 π m B S δ k ) 2 ] + exp [ 1 2 Δ 2 w 0 2 ( y + 2 π m B S δ k ) 2 ] } × ω 0 T 0 2 2 B T exp ( ω 0 2 4 a B T 2 t 2 ) q = 0 n ( n q ) H 2 q ( ω 0 t / 2 a B T ) q ! ( 2 T c ) 2 q a q + 1 / 2 .
Further, when y = 0, t = 0 and n = 0, Equation (21) can be expressed as
I ( r , t , z ) = 1 2 Δ 2 exp ( π m 2 z 2 k 2 δ 2 w 0 2 Δ 2 ) ω 0 T 0 2 2 a B T × { exp [ 1 2 Δ 2 w 0 2 ( x 2 π m z δ k ) 2 ] + exp [ 1 2 Δ 2 w 0 2 ( x + 2 π m z δ k ) 2 ] } .
It can be seen from Equation (22) that the spatiotemporal intensity I(x, y, t, z) is a superposition of two exponential functions symmetric with respect to the coordinate axis x = 0. When the propagation distance z is short, the two exponential functions are similar, hence the two beams cannot be distinguished. However, when z grows large enough, the two beams become visibly different. Further, when beam order m increases, the position where the beam starts to split moves closer to the source plane z = 0. Moreover, when t ≠ 0, the decrease of the spatiotemporal intensity is attributed to the last exponential function in Equation (21), which plays an important role in determining the value of spatiotemporal intensity.
Figure 3 and Figure 4 show the density plots for the normalized spatiotemporal intensity I(x, y, t, z) of the SCTLGSM pulsed beam in the xy plane at some propagation distances z = 0, 30 m, 40 m, 60 m and for m = 1 (Figure 3) and m = 2 (Figure 4), respectively. The other values of parameters are t = 0, s = 0 and n = 0. It is shown that the SCTLGSM pulsed beam exhibits self-splitting properties upon propagation in water, i.e., the initial beam spot evolves into four sub-beam spots in the far field depending on the value of beam order m. For larger m, the self-splitting time is smaller. These results can be explained from the Equation (22). When m increases, coefficient (2π)1/2m/δk of z in the exponential function increases as well. Hence, the rate of change in the spatiotemporal intensity increases with the increasing propagation distance.
Figure 5 illustrates the evolution of the spatiotemporal intensity I(x, y, t, z) of the SCTLGSM pulsed beam as a function of propagation distance z and time t, with different beam orders n = 0, 1, 2 and 3. The other calculation parameters are x = 0, y = 0, s = −0.4ps−2 and m = 0. One can see from Figure 5 that the beam exhibits self-splitting behavior when n > 0. Specifically, the beam splits into n + 1 sub-beams with increasing propagation distance z, being in good agreement with the results of the Laguerre–Gaussian Schell-model pulsed beam in time domain [31]. While for n = 0, i.e., for the conventional Gaussian-correlated Schell-model pulsed beam, there is no self-splitting upon propagation. These results can be interpreted as follows. For n = 0, the last term of Equation (21) can be re-expressed as
I n = 0 = ω 0 T 0 2 2 a B T exp ( ω 0 2 4 a B T 2 t 2 ) ,
while for n=1,
I n = 1 = ω 0 T 0 2 2 a B T exp ( ω 0 2 4 a B T 2 t 2 ) [ ω 0 2 4 a 2 B T 2 T c 2 t 2 + 1 1 2 a T c 2 ] ,
and for n = 2,
I n = 2 = ω 0 T 0 2 2 a B T exp ( ω 0 2 4 a B T 2 t 2 ) [ ω 0 4 32 a 4 B T 4 T c 4 t 4 + ( 1 3 2 a T c 2 ) ω 0 2 4 a 2 B T 2 T c 2 t 2 + 1 1 2 a T c 2 + 3 8 a 2 T c 4 ] ,
here Equations (24) and (25) can be re-written as:
I n = 1 = ε exp ( φ t 2 ) [ ξ 1 t 2 + ψ 1 ] ,
I n = 2 = ε exp ( φ t 2 ) [ ζ t 4 + ξ 2 t 2 + ψ 2 ] ,
ε = ω 0 T 0 2 2 a B T ,   φ = ω 0 2 4 a B T 2 ,   ξ 1 = ω 0 2 4 a 2 B T 2 T c 2 ,   ψ 1 = 1 1 2 a T c 2 ,
ζ = ω 0 4 32 a 4 B T 4 T c 4 ,   ξ 2 = ( 1 3 2 a T c 2 ) ω 0 2 4 a 2 B T 2 T c 2 ,   ψ 2 = 1 1 2 a T c 2 + 3 8 a 2 T c 4 .
From Equation (24) we see that In=1 is a product of functions exp ( φ t 2 ) and ξ 1 t 2 + ψ 1 . Thus, letting
d I n = 1 d t = 2 ε exp [ φ t 2 ] ( φ ξ 1 t 2 + ξ 1 φ ψ 1 ) = 0 ,
and solving Equation (30), In=1 will obtain two symmetric intensity maxima at t = ± ( ξ 1 φ ψ 1 ) / φ ξ 1 when ξ 1 φ ψ 1 . Therefore, the beam starts splitting into two sub-beams in Figure 5b. In addition, when z is very short, the influence of term ξ1 t2+Ψ1 can be neglected compared with that of exp(−φt2). This is why there is only one beam when z is much closer to the source plane. Similarly, for n = 2 and n = 3, we can obtain 3 intensity maxima and 4 intensity maxima in the far field, respectively.
Figure 6 shows the evolution of spatiotemporal intensity I(x, y, t, z) of the SCTLGSM pulsed beam as a function of propagation distance z and time t with different values of chirp coefficient s = 0, −0.4ps−2, −0.8ps−2, −1.2ps−2, and for n = 0, m = 0, x = 0, y = 0. It is shown that the self-focusing phenomenon takes place upon propagation when chirp coefficient s < 0, and it becomes more noticeable with decreasing chirp coefficient s, while the position of the focus shifts towards the source plane. This result can be explained as follows. When n = 0, m = 0, x = 0, y = 0, Equation (23) can be expressed as
I n = 0 = ω 0 T 0 2 2 a B T exp ( ω 0 2 4 a B T 2 t 2 ) = T 0 2 T ( z ) exp ( 1 T 2 ( z ) t 2 ) ,
where
T 2 ( z ) = ( 2 T 0 2 + 2 T c 2 ) β 2 2 z 2 + ( 1 + s β 2 z ) 2 T 0 2 2 = ( 2 T 0 2 + 2 T c 2 + T 0 2 s 2 2 ) β 2 2 z 2 + s T 0 2 β 2 z + T 0 2 2 .
One can see from Equation (32) that T 2 ( z ) is a quadratic function of z, and will reach minima when
z min = s T 0 2 β 2 T 0 2 s 2 + 4 / T 0 2 + 4 / T c 2 .
Hence, the self-focusing phenomenon appears upon propagation when chirp coefficient s < 0. Moreover, with increasing |s|, zmin will decrease, which can be seen from Equation (33), because for large |s|, 4 / T 0 2 + 4 / T c 2   can be omitted. Equation (33) can be expressed as zminβ2/|s| which is why the position of focus shifts towards the source plane with increasing |s|. Here the influence of s on propagation in the temporal ABCD matrix is equivalent to a lens. Hence, the self-focusing phenomenon appears.
Figure 7 shows the evolution of spatiotemporal intensity I(x, y, t, z) of SCTLGSM pulsed beam as a function of propagation distance z and time t with different chirp coefficients s = 0, −0.4ps−2, −0.8ps−2, −1.2ps−2, and for n = 1, m = 0, x = 0, y = 0. It is shown that the self-focusing and the self-splitting phenomena take place at the same time when chirp coefficient s < 0. Due to the shift of the position of self-focusing, there are three intensity peaks appearing for larger |s|.
Figure 8a–e shows the density plots for the normalized spatiotemporal intensity I(x, y, t, z) of the SCTLGSM pulsed beam in the xt plane at fixed propagation distances z = 0, 20 m, 30 m, 40 m, 50 m and for s = −0.4ps−2, n = 0, m = 0, y = 0. As can be seen, there is a circular intensity distribution in the xt plane for z = 0. With the increasing propagation distance z, the circular intensity distribution becomes elliptical (Figure 8c–e), where the self-focusing appears in the coordinate axis t. Figure 8f shows the normalized spatiotemporal intensity I(x, y, t, z) of the SCTLGSM pulsed beam as a function of propagation distance z for different values of the chirp coefficient s = 0, −0.4ps−2, −0.8ps−2, −1.2ps−2, where n = 0, m = 0, x = 0, y = 0. One can see that the self-focusing phenomenon becomes more and more noticeable with increasing chirp coefficient |s|. Moreover, for the case of m = n = 1, shown in Figure 9a–c, the beam splits first in the coordinate axis t, and then splits in the coordinate axis x (Figure 9d). In the far field, the beam splits into four sub-beam spots in the xt plane. Figure 9f shows the normalized spatiotemporal intensity I(x, y, t, z) of the SCTLGSM pulsed beam as a function of propagation distance z for different chirp coefficient values: s = 0, −0.4ps−2, −0.8ps−2, −1.2ps−2, where n = 1, m = 1, x = 0, y = 0, t = 0. One can find that the self-focusing phenomenon is more noticeable for s = −1.2ps−2. Comparison of Figure 8 with Figure 9 implies that the beam orders m and n play an important role in determining the spatiotemporal intensity distribution in xt plane.
Figure 10a–d includes the density plots for the normalized spatiotemporal intensity I(x, y, t, z) of the SCTLGSM pulsed beam in the xt plane for different beam orders m = n = 0, m = n = 1, m = n = 2, m = n = 3 at z = 50 m, and for s = −0.4 ps−2, y = 0. It is shown that no beam splitting occurs in the far field for m = n = 0. However, for m = n = 1, the initial beam starts to split into two sub-beams in coordinate axis t and coordinate axis x, respectively. With increasing beam order, such as m = n = 2, the initial beam starts to split into three sub-beams in coordinate axis t and two sub-beams in coordinate axis x, respectively. In addition, for m = n = 3, the initial beam starts to split into four sub-beams in coordinate axis t and two sub-beams coordinate axis x, respectively. These results are consistent with those of Figure 5.

4. Concluding Remarks

In this paper, a new class of partially coherent pulsed beams with spatial cosine-Gaussian and temporal Laguerre–Gaussian correlations is introduced. The analytical expressions of the spatiotemporal intensity of the SCTLGSM pulsed beam through the ABCD optical system on passing in a diffractive/dispersive medium are obtained and used to investigate the evolution of the spatiotemporal intensity of the SCTLGSM pulsed beam in still, pure water. It is found that the SCTLGSM pulsed beam exhibits the phenomena of simultaneous spatiotemporal self-splitting and self-focusing. When the chirp coefficient is trivial, s = 0, in the xz plane, the beam always splits into two sub-beams, when m > 0, while in the tz plane, the beam splits into n + 1 sub-beams. When s < 0, the spatiotemporal self-focusing phenomenon takes place. The influence of s on propagation in a dispersive medium is equivalent to the effect of a lens in a diffractive medium. Physically, the special structuring of the spatial coherence and temporal coherence states result in the spatiotemporal self-splitting phenomenon of the SCTLGSM pulsed beam. Furthermore, the pulse chirp leads to the spatiotemporal self-focusing phenomenon of the SCTLGSM pulsed beam. The results obtained in this paper can be readily extended to other dispersive media and/or other temporal ABCD optical systems.
More importantly, we have explored the possibility of shaping of optical pulsed signals by source correlations in one of the most frequently encountered natural environments—pure water. Such a method can be of interest for rapidly developing underwater optical technologies, including wireless communications and directed energy. In particular, one can obtain the desired beam profiles delivered at suitable time instants by modulating the spatial and temporal coherence structures of the pulsed beams and/or adding a linear chirp.
It should be stressed that our work was carried out under the assumption that the refractive index of the dispersive medium is linear. The self-focusing effect comes from the source-induced spatial and temporal interference effects (pulse chirp) and should not be associated with the medium’s non-linear effects. In cases when the dispersive medium is nonlinear, the four-wave-mixing [53] and formation of optical solitons must be additionally discussed [54]. We note that there is a report available on propagation of the spatially partially coherent beams with non-Gaussian correlation in oceanic turbulence [55] in which the correlation-induced phenomenon of self-focusing also appears, confirming the validity of our results.

Author Contributions

C.D. (Investigation, Methodology, Writing—Original Draft Preparation, Writing—Review & Editing); O.K. (Formal Analysis, Writing—Review & Editing); D.H. (Formal Analysis, Validation); Z.Z. (Formal Analysis, Validation); L.P. (Resources, Formal Analysis, Project Administration). All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded partially by the National Natural Science Foundation of China (61575091), China Scholarship Council (202008410578), 2020 Central Plains Talents Program of Henan, High-level talents international training program and Excellent Overseas Visiting Scholar Program of Henan.

Institutional Review Board Statement

This study did not involve humans or animals.

Informed Consent Statement

This study did not involve humans.

Data Availability Statement

The study did not report any data.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In this appendix, we will justify the validity of adopted mutual coherence function with spatial cosine-Gaussian and temporal Laguerre–Gaussian correlations. The mutual coherence function of the SCTLGSM pulsed beam at the source plane z = 0 is defined as
Γ 0 ( ρ 1 , τ 1 ; ρ 2 , τ 2 ) = S ( ρ 1 , ρ 2 ) T ( τ 1 , τ 2 ) ,
where the spatial and the temporal mutual coherence functions can be written in the following forms, respectively,
S ( ρ 1 , ρ 2 ) = p S ( ε ) H * ( ρ 1 , ε ) H ( ρ 2 , ε ) d 2 ε ,
T ( τ 1 , τ 2 ) = p T ( v ) h * ( τ 1 , v ) h ( τ 2 , v ) d v ,
where pS and pT are non-negative weight functions, and H(ρ, ε) and h(t, v) are arbitrary kernels in the space domain and the time domain, respectively. We assume that
H ( ρ , ε ) = i π δ σ exp ( ρ 2 4 w 0 2 ) exp [ i δ σ ( ε 2 2 ρ ε ) ] ,
h ( τ , v ) = exp ( τ 2 T 0 2 ) exp [ i 2 2 T c ( v 2 2 τ v ) ] exp ( i ω 0 τ ) ,
and choose the weighting functions pS and pT as follows:
p S ( ε ) = cosh ( 2 n 2 π ε x / σ ) cosh ( 2 n 2 π ε y / σ ) exp ( 2 ε 2 σ 2 ) ,
p T ( v ) = 1 2 n n ! π H n 2 ( v ) exp ( v 2 ) ,
where parameter σ is a factor relating to beam width [12], and Hn is a Hermite polynomial of order n [31]. Here, we choose the mutual coherence functions in the spatial domain and in time domain from Refs. [12,31], respectively.
The cross-spectral density function W 0 ( ρ 1 , ω 1 ; ρ 2 , ω 2 ) is obtained from the mutual coherence function using the generalized Wiener–Khintchine theorem
W 0 ( ρ 1 , ω 1 ; ρ 2 , ω 2 ) = 1 ( 2 π ) 2 Γ 0 ( ρ 1 , τ 1 ; ρ 2 , τ 2 ) exp [ i ( ω 1 t 1 ω 2 t 2 ) ] d t 1 d t 2 .
After tedious calculation and considering the specific Equations (2) and (3), one obtains
W 0 ( ρ 1 , τ 1 ; ρ 2 , τ 2 ) = S ( ρ 1 , ρ 2 ) W ( ω 1 , ω 2 ) ,
where
W ( ω 1 , ω 2 ) = T 0 2 π 2 n n ! Ω 0 exp [ ( ( ω 1 ω 0 ) 2 + ( ω 2 ω 0 ) 2 Ω 0 2 ) ] exp [ ( ω 1 ω 2 ) 2 2 Ω c 2 ] × k = 0 n 2 k k ! ( 1 + Ω c 2 / Ω 0 2 ) n k ( n k ) 2 H 2 ( n k ) [ ω 1 + ω 2 2 ω 0 2 Ω 0 ] ,
and the pulse duration T0, temporal coherence length Tc, spectral bandwidth Ω and the spectral coherence width Ωc are related by
Ω 0 2 4 = 1 T 0 2 + 1 T c 2 ,   Ω c = T c T 0 Ω 0 .
On the one hand, Equations (A2) and (A3) guarantee that the mutual coherence function of the SCTLGSM pulsed beam is a non-negative kernel and is therefore physically realizable [56]. On the other hand, they show how such a pulsed beam can be realized by transformations of an incoherent field E i ( ρ , τ ) with the mutual coherence function
Γ i ( ρ , τ ; ρ , τ ) = E i * ( ρ , τ ) E i ( ρ , τ ) = p S ( ρ ) p T ( τ ) δ ( ρ ρ ) δ ( τ τ ) .
The field E i ( ρ , τ ) can be obtained from a perfectly spatially and temporally incoherent source, e.g., a broadband light-emitting diode, by spatial and temporal amplitude modulation with the modulation profiles p S ( ρ ) and p T ( τ ) respectively.
In order to obtain the SCTLGSM pulsed beam, we need to pass the field E i ( ρ , τ ) through a spatial and temporal imaging system with the transfer function H ( ρ , ε ) h ( τ , v ) so that the output field becomes
E 0 ( ρ , τ ) = + H ( ρ , ε ) h ( τ , v ) E i ( ε , v ) d 2 ε d v ,
which, according to Equations (A2), (A3) and (A12), will give us a field with the desired mutual coherence function Γ 0 ( ρ 1 , τ 1 ; ρ 2 , τ 2 ) =   E 0 * ( ρ 1 , τ 1 ) E 0 ( ρ 2 , τ 2 ) , given by Equation (A1). The spatial transfer function, Equation (A4), can be regarded as a product of three terms (up to an irrelevant constant factor) [12]: (i) the transfer function corresponding to diffraction at a distance L, H diff ( ρ , ε ) = exp [ i k 2 L ( ρ ε ) 2 ] , (ii) the phase factor of a thin lens with the focal distance f, H lens ( ρ ) = exp [ i k 2 f ρ 2 ] , and (iii) the transmission coefficient of a Gaussian amplitude filter with the intensity standard deviation w 0 , H filter ( ρ ) = exp [ ρ 2 4 w 0 2 ] , where k is the wave number corresponding to the carrier frequency and L = f = δ σ k / 2 . The above considerations are valid for the fields of bandwidth much less than the carrier frequency.
In a similar way, the temporal transfer function, Equation (A5), can be regarded as a product of three terms corresponding to standard components of a temporal imaging system [57] (again, up to a constant factor): (i) the transfer function corresponding to propagation through a dispersive medium with the total group delay dispersion (GDD) D, h disp ( τ , v ) = exp [ i 2 D ( τ v ) 2 ] , (ii) the phase factor of a time lens with the focal GDD D f , h lens ( τ ) = exp [ i 2 D f τ 2 ] , and (iii) the transmission coefficient of a Gaussian amplitude modulator with the intensity standard deviation T 0 / 2 , h mod ( τ ) = exp [ τ 2 T 0 2 ] , where D = D f = T c / 2 .
Thus, a combination of elementary transformations described above gives a possible way of generation of a SCTLGSM pulsed beam.
In the Figure A1, we present an explanatory figure for the source of SCTLGSM pulsed beam with spatial cosine-Gaussian and temporal Laguerre–Gaussian correlations and show the spectral density profiles of the source.
Figure A1. (a) Illustrating figure of the SCTLGSM pulsed beam propagating in water. (b) Spatial degree of coherence of the pulsed beam source as a function of ( x 2 x 1 ) / δ . Solid black: m = 0; Dashed red: m = 1; Dash-dotted green: m = 2; Dotted blue: m = 3. (c) Temporal degree of coherence of the pulsed beam source as a function of ( τ 2 τ 1 ) / T c . Solid black: n = 0; Dashed red: n = 1; Dash-dotted green: n = 2; Dotted blue: n = 3. (d) Spectral density of the pulsed beam source as a function of (ωω0)/Ω0. Solid black: Ωc0 =1; Dashed red: Ωc0 = 1; Dash-dotted green: Ωc0 = 0.5; Dotted blue: Ωc0 = 0.1.
Figure A1. (a) Illustrating figure of the SCTLGSM pulsed beam propagating in water. (b) Spatial degree of coherence of the pulsed beam source as a function of ( x 2 x 1 ) / δ . Solid black: m = 0; Dashed red: m = 1; Dash-dotted green: m = 2; Dotted blue: m = 3. (c) Temporal degree of coherence of the pulsed beam source as a function of ( τ 2 τ 1 ) / T c . Solid black: n = 0; Dashed red: n = 1; Dash-dotted green: n = 2; Dotted blue: n = 3. (d) Spectral density of the pulsed beam source as a function of (ωω0)/Ω0. Solid black: Ωc0 =1; Dashed red: Ωc0 = 1; Dash-dotted green: Ωc0 = 0.5; Dotted blue: Ωc0 = 0.1.
Photonics 08 00102 g0a1

References

  1. Korotkova, O. Random Light Beams: Theory and Applications; CRC Press: Boca Raton, FL, USA, 2013. [Google Scholar]
  2. Cai, Y.J.; Chen, Y.H.; Wang, F. Generation and propagation of partially coherent beams with nonconventional correlation functions: A review [Invited]. J. Opt. Soc. Am. A 2014, 31, 2083–2096. [Google Scholar] [CrossRef] [PubMed]
  3. Cai, Y.J.; Wang, F.; Chen, Y.H.; Liu, X.L. Recent advances in partially coherent beams with prescribed degrees of coherence. Sci. Found. China 2017, 25, 38–60. [Google Scholar]
  4. Cai, Y.J.; Chen, Y.H.; Yu, J.Y.; Liu, X.L.; Liu, L. Generation of Partially Coherent Beams. Prog. Opt. 2017, 62, 157–223. [Google Scholar]
  5. Wan, L.P.; Zhao, D.M. Controllable rotating Gaussian Schell-model beams. Opt. Lett. 2019, 44, 735–738. [Google Scholar] [CrossRef] [PubMed]
  6. Zhang, J.P.; Wang, J.; Huang, H.K.; Wang, H.Y.; Zhu, S.J.; Li, Z.H.; Lu, J. Propagation Characteristics of a Twisted Cosine-Gaussian Correlated Radially Polarized Beam. Appl. Sci. 2018, 8, 1485. [Google Scholar] [CrossRef] [Green Version]
  7. Zheng, S.M.; Yuan, C.F.; Ji, X.L.; Wang, T. Ring-shaped twisted Gaussian Schell-model array beams. J. Opt. Soc. Am. A 2020, 37, 444–449. [Google Scholar] [CrossRef]
  8. Hyde IV, M.W.; Bose-Pillai, S.; Voelz, D.G.; Xiao, X. A fast and efficient method for producing partially coherent sources. J. Opt. 2017, 19, 025601. [Google Scholar] [CrossRef]
  9. Wang, J.; Yang, S.; Guo, M.J.; Feng, Z.F.; Li, J.H. Change in phase singularities of a partially coherent Gaussian vortex beam propagating in a GRIN fiber. Opt. Express 2020, 28, 4661–4673. [Google Scholar] [CrossRef]
  10. Mei, Z.R.; Korotkova, O. Cosine-Gaussian Schell-model sources. Opt. Lett. 2013, 38, 2578–2580. [Google Scholar] [CrossRef]
  11. Chen, Y.H.; Gu, J.X.; Wang, F.; Cai, Y.J. Self-splitting properties of a Hermite-Gaussian correlated Schell-model beam. Phys. Rev. A 2015, 91, 013823. [Google Scholar] [CrossRef]
  12. Liang, C.H.; Wang, F.; Liu, X.L.; Cai, Y.J.; Korotkova, O. Experimental generation of cosine-Gaussian-correlated Schell-model beams with rectangular symmetry. Opt. Lett. 2014, 39, 769–772. [Google Scholar] [CrossRef] [PubMed]
  13. Lajunen, H.; Saastamoinen, T. Propagation characteristics of partially coherent beams with spatially varying correlations. Opt. Lett. 2011, 36, 4104–4106. [Google Scholar] [CrossRef] [PubMed]
  14. Chen, Y.H.; Ponomarenko, S.A.; Cai, Y.J. Self-steering partially coherent beams. Sci. Rep. 2017, 7, 39957. [Google Scholar] [CrossRef] [Green Version]
  15. Sahin, S.; Korotkova, O. Light sources generating far fields with tunable flat profiles. Opt. Lett. 2012, 37, 2970–2972. [Google Scholar] [CrossRef]
  16. Liang, C.H.; Khosravi, R.; Liang, X.; Kacerovská, B.; Monfared, Y.E. Standard and elegant higher-order Laguerre-Gaussian correlated Schell-model beams. J.Opt. 2019, 21, 085607. [Google Scholar] [CrossRef]
  17. Korotkova, O. Random sources for rectangular far fields. Opt. Lett. 2014, 39, 64–67. [Google Scholar] [CrossRef]
  18. Wan, L.P.; Zhao, D.M. Optical coherence grids and their propagation characteristics. Opt. Express 2018, 26, 2168–2180. [Google Scholar] [CrossRef]
  19. Chen, Y.H.; Wang, F.; Zhao, C.L.; Cai, Y.J. Experimental demonstration of a Laguerre-Gaussian correlated Schell-model vortex beam. Opt. Express 2014, 22, 5826–5838. [Google Scholar] [CrossRef]
  20. Gbur, G. Partially coherent beam propagation in atmospheric turbulence [Invited]. J. Opt. Soc. Am. A 2014, 31, 2038–2045. [Google Scholar] [CrossRef] [Green Version]
  21. Singh, M.; Lajunen, H.; Tervo, J.; Turunen, J. Imaging with partially coherent light: Elementary-field approach. Opt. Express 2015, 23, 28132–28140. [Google Scholar] [CrossRef]
  22. Wu, G.F.; Cai, Y.J. Detection of a semirough target in turbulent atmosphere by a partially coherent beam. Opt. Lett. 2011, 36, 1939–1941. [Google Scholar] [CrossRef]
  23. Korotkova, O.; Chen, X. Phase structuring of the complex degree of coherence. Opt. Lett. 2018, 43, 4727–4730. [Google Scholar] [CrossRef] [PubMed]
  24. Korotkova, O.; Chen, X.; Setälä, T. Electromagnetic Schell-model beams with arbitrary complex correlation states. Opt. Lett. 2019, 44, 4945–4948. [Google Scholar] [CrossRef] [PubMed]
  25. Pääkkönen, P.; Turunen, J.; Vahimaa, P.; Friberg, A.T.; Wyrowski, F. Partially coherent Gaussian pulses. Opt. Commun. 2002, 204, 53–58. [Google Scholar] [CrossRef]
  26. Lin, Q.; Wang, L.G.; Zhu, S.Y. Partially coherent light pulse and its propagation. Opt. Commun. 2003, 219, 65–70. [Google Scholar] [CrossRef]
  27. Turunen, J. Low-coherence laser beams. In Laser Beam Propagation: Generation and Propagation of Customized Light; Forbes, A., Ed.; CRC Press: Boca Raton, FL, USA, 2014. [Google Scholar]
  28. Bourassin-Bouchet, C.; Couprie, M.E. Partially coherent ultrafast spectrography. Nat. Commun. 2015, 6, 6465. [Google Scholar] [CrossRef] [Green Version]
  29. Lajunen, H.; Saastamoinen, T. Non-uniformly correlated partially coherent pulses. Opt. Express 2013, 21, 190–195. [Google Scholar] [CrossRef]
  30. Ding, C.L.; Korotkova, O.; Pan, L.Z. The control of pulse profiles with tunable temporal coherence. Phys. Lett. A 2014, 378, 1687–1690. [Google Scholar] [CrossRef]
  31. Ding, C.L.; Koivurova, M.; Turunen, J.; Pan, L.Z. Temporal self-splitting of optical pulses. Phys. Rev. A 2018, 97, 053838. [Google Scholar] [CrossRef] [Green Version]
  32. Liu, X.Y.; Zhao, D.M.; Tian, K.H.; Pan, W.Q.; Zhang, K.W. Gaussian temporal modulation for the behavior of multi-sinc Schell-model pulses in dispersive media. Opt. Commun. 2018, 416, 160–165. [Google Scholar] [CrossRef]
  33. Tang, M.M.; Zhao, D.M.; Zhu, Y.B.; Ang, L.-K. Electromagnetic sinc Schell-model pulses in dispersive media. Phys. Lett. A 2016, 380, 794–797. [Google Scholar] [CrossRef]
  34. Ding, C.L.; Korotkova, O.; Zhao, D.M.; Li, D.L.; Zhao, Z.G.; Pan, L.Z. Propagation of temporal coherence gratings in dispersive medium with a chirper. Opt. Express 2020, 28, 7463–7474. [Google Scholar] [CrossRef]
  35. Zhang, Y.T.; Ding, C.L.; Hyde, M.W., IV; Korotkova, O. Non-stationary pulses with complex-valued temporal degree of coherence. J.Opt. 2020, 22, 105607. [Google Scholar] [CrossRef]
  36. Ding, C.L.; Korotkova, O.; Zhang, Y.T.; Pan, L.Z. Cosine-Gaussian correlated Schell-model pulsed beams. Opt. Express 2014, 22, 931–942. [Google Scholar] [CrossRef] [PubMed]
  37. Torres-Company, V.; Mínguez-Vega, G.; Lancis, J.; Friberg, A.T. Controllable generation of partially coherent light pulses with direct space-to-time pulse shaper. Opt. Lett. 2007, 32, 1608–1610. [Google Scholar] [CrossRef] [PubMed]
  38. Christov, I.P. Propagation of Partially Coherent Light Pulses. Opt. Acta Int. J. Opt. 1986, 33, 63–72. [Google Scholar] [CrossRef]
  39. Lajunen, H.; Vahimaa, P.; Tervo, J. Theory of spatially and spectrally partially coherent pulses. J. Opt. Soc. Am. A 2005, 22, 1536–1545. [Google Scholar] [CrossRef]
  40. Torres-Company, V.; Lajunen, H.; Lancis, J.; Friberg, A.T. Ghost interference with classical partially coherent light pulses. Phys. Rev. A 2008, 77, 043811. [Google Scholar] [CrossRef] [Green Version]
  41. Shirai, T.; Setälä, T.; Friberg, A.T. Temporal ghost imaging with classical non-stationary pulsed light. J. Opt. Soc. Am. B 2010, 27, 2549–2555. [Google Scholar] [CrossRef]
  42. Huang, D.; Swanson, E.A.; Lin, C.P.; Schuman, J.S.; Stinson, W.G.; Chang, W.; Hee, M.R.; Flotte, T.; Gregory, K.; Puliafito, C.A. Optical coherence tomography. Science 1991, 254, 1178–1181. [Google Scholar] [CrossRef] [Green Version]
  43. Bertolotti, M.; Ferrari, A.; Sereda, L. Coherence properties of nonstationary polychromatic light sources. J. Opt. Soc. Am. B 1995, 12, 341–347. [Google Scholar] [CrossRef]
  44. Sereda, L.; Bertolotti, M.; Ferrari, A. Coherence properties of nonstationary light wave fields. J. Opt. Soc. Am. A 1998, 15, 695–705. [Google Scholar] [CrossRef]
  45. Dijaili, S.P.; Dienes, A.; Smith, J.S. ABCD matrices for dispersive pulse propagation. IEEE J. Quantum Elect. 1990, 26, 1158–1164. [Google Scholar] [CrossRef]
  46. Mandel, L.; Wolf, E. Optical Coherence and Quantum Optics; Cambridge University: Cambridge, UK, 1995. [Google Scholar]
  47. Luo, M.L.; Zhao, D.M. Characterizing the polarization and cross-polarization of electromagnetic vortex pulses in the space-time and space-frequency domain. Opt. Express 2015, 23, 4153–4162. [Google Scholar] [CrossRef]
  48. Turunen, J.; Friberg, A.T. Matrix representation of Gaussian Schell-model beams in optical systems. Opt. Laser Technol. 1986, 18, 259–267. [Google Scholar] [CrossRef]
  49. Lajunen, H.; Turunen, J.; Vahimaa, P.; Tervo, J.; Wyrowski, F. Spectrally partially coherent pulse trains in dispersive media. Opt. Commun. 2005, 255, 12–22. [Google Scholar] [CrossRef]
  50. Daimon, M.; Masumura, A. Measurement of the refractive index of distilled water from the near-infrared region to the ultraviolet region. Appl. Opt. 2007, 46, 3811–3820. [Google Scholar] [CrossRef] [PubMed]
  51. Lu, W.; Liu, L.R.; Sun, J.F. Influence of temperature and salinity fluctuations on propagation behaviour of partially coherent beams in oceanic turbulence. J. Opt. A Pure Appl. Opt. 2006, 8, 1052. [Google Scholar] [CrossRef]
  52. Pope, R.M.; Fry, E.S. Absorption spectrum (380–700 nm) of pure water. II. Integrating cavity measurements. Appl. Opt. 1997, 36, 8710–8723. [Google Scholar] [CrossRef]
  53. Schnebelin, C.; Cassagne, C.; de Araújo, C.B.; Boudebs, G. Measurements of the third- and fifth-order optical nonlinearities of water at 532 and 1064 nm using the D4σ method. Opt. Lett. 2014, 39, 5046–5049. [Google Scholar] [CrossRef] [Green Version]
  54. Sedov, E.V.; Redyuk, A.A.; Fedoruk, M.P.; Gelash, A.A.; Frumin, L.L.; Turitsyn, S.K. Soliton content in the standard optical OFDM signal. Opt. Lett. 2018, 43, 5985–5988. [Google Scholar] [CrossRef]
  55. Wang, M.; Zhang, J.; Wang, J. Propagation characteristics of non-uniformly Sinc-correlated blue-green laser beam through oceanic turbulence. Infrared Laser Eng. 2020, 49, 20190370. [Google Scholar] [CrossRef]
  56. Gori, F.; Santarsiero, M. Devising genuine spatial correlation functions. Opt. Lett. 2007, 32, 3531–3533. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Patera, G.; Horoshko, D.B.; Kolobov, M.I. Space-time duality and quantum temporal imaging. Phys. Rev. A 2018, 98, 053815. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Evolution of spatiotemporal intensity I(x, y, t, z) (Equation (20)) of SCTLGSM pulsed beam as a function of propagation distance z and horizontal coordinate x with different beam orders: (a) m = 0, (b) m = 1, (c) m = 2 and (d) m = 3.
Figure 1. Evolution of spatiotemporal intensity I(x, y, t, z) (Equation (20)) of SCTLGSM pulsed beam as a function of propagation distance z and horizontal coordinate x with different beam orders: (a) m = 0, (b) m = 1, (c) m = 2 and (d) m = 3.
Photonics 08 00102 g001
Figure 2. Evolution of spatiotemporal intensity I(x, y, t, z) (Equation (20)) of spatial cosine-Gaussian and temporal Laguerre–Gaussian correlated Schell-model (SCTLGSM) pulsed beam as a function of propagation distance z and horizontal coordinate x with different time (a) t = 2 ps and (b) t = 5 ps, and for m = 1, s = 0.
Figure 2. Evolution of spatiotemporal intensity I(x, y, t, z) (Equation (20)) of spatial cosine-Gaussian and temporal Laguerre–Gaussian correlated Schell-model (SCTLGSM) pulsed beam as a function of propagation distance z and horizontal coordinate x with different time (a) t = 2 ps and (b) t = 5 ps, and for m = 1, s = 0.
Photonics 08 00102 g002
Figure 3. Density plots for the normalized spatiotemporal intensity I(x, y, t, z) (Equation (20)) of SCTLGSM pulsed beam in the xy plane at some propagation distances: (a) z = 0, (b) z = 30 m, (c) z = 40 m, (d) z = 60 m and for t = 0, s = 0, m = 1.
Figure 3. Density plots for the normalized spatiotemporal intensity I(x, y, t, z) (Equation (20)) of SCTLGSM pulsed beam in the xy plane at some propagation distances: (a) z = 0, (b) z = 30 m, (c) z = 40 m, (d) z = 60 m and for t = 0, s = 0, m = 1.
Photonics 08 00102 g003
Figure 4. Density plots for the normalized spatiotemporal intensity I(x, y, t, z) (Equation (20)) of SCTLGSM pulsed beam in the xy plane at some propagation distances (a) z = 0, (b) z = 30 m, (c) z = 40 m, (d) z = 60 m and for t = 0, s = 0, m = 2.
Figure 4. Density plots for the normalized spatiotemporal intensity I(x, y, t, z) (Equation (20)) of SCTLGSM pulsed beam in the xy plane at some propagation distances (a) z = 0, (b) z = 30 m, (c) z = 40 m, (d) z = 60 m and for t = 0, s = 0, m = 2.
Photonics 08 00102 g004
Figure 5. Evolution of spatiotemporal intensity I(x, y, t, z) (Equation (20)) of SCTLGSM pulsed beam as a function of propagation distance z and time t with different beam orders (a) n = 0, (b) n =1, (c) n =2 and (d) n =3.
Figure 5. Evolution of spatiotemporal intensity I(x, y, t, z) (Equation (20)) of SCTLGSM pulsed beam as a function of propagation distance z and time t with different beam orders (a) n = 0, (b) n =1, (c) n =2 and (d) n =3.
Photonics 08 00102 g005
Figure 6. Evolution of spatiotemporal intensity I(x, y, t, z) (Equation (20)) of SCTLGSM pulsed beam as a function of propagation distance z and time t with different chirp coefficients (a) s = 0, (b) s =−0.4 ps−2, (c) s =−0.8 ps−2, (d) s =−1.2 ps−2 and for n = 0, m = 0.
Figure 6. Evolution of spatiotemporal intensity I(x, y, t, z) (Equation (20)) of SCTLGSM pulsed beam as a function of propagation distance z and time t with different chirp coefficients (a) s = 0, (b) s =−0.4 ps−2, (c) s =−0.8 ps−2, (d) s =−1.2 ps−2 and for n = 0, m = 0.
Photonics 08 00102 g006
Figure 7. Evolution of spatiotemporal intensity I(x, y, t, z) (Equation (20)) of SCTLGSM pulsed beam as a function of propagation distance z and time t with different chirp coefficients (a) s = 0, (b) s =-0.4 ps−2, (c) s =-0.8 ps−2, (d) s =−1.2 ps−2 and for n = 1, m = 0.
Figure 7. Evolution of spatiotemporal intensity I(x, y, t, z) (Equation (20)) of SCTLGSM pulsed beam as a function of propagation distance z and time t with different chirp coefficients (a) s = 0, (b) s =-0.4 ps−2, (c) s =-0.8 ps−2, (d) s =−1.2 ps−2 and for n = 1, m = 0.
Photonics 08 00102 g007
Figure 8. Density plots for the normalized spatiotemporal intensity I(x, y, t, z) (Equation (20)) of SCTLGSM pulsed beam in the xt plane at some propagation distances (a) z = 0, (b) z =20 m, (c) z =30 m, (d) z =40 m, (e) z =50 m and for s = −0.4ps−2, n = 0, m = 0. (f) Normalized intensity I(x, y, t, z) (Equation (22)) of SCTLGSM pulsed beam as a function of propagation distance z for different chirp coefficients s = 0, −0.4ps−2, −0.8ps−2, −1.2ps−2.
Figure 8. Density plots for the normalized spatiotemporal intensity I(x, y, t, z) (Equation (20)) of SCTLGSM pulsed beam in the xt plane at some propagation distances (a) z = 0, (b) z =20 m, (c) z =30 m, (d) z =40 m, (e) z =50 m and for s = −0.4ps−2, n = 0, m = 0. (f) Normalized intensity I(x, y, t, z) (Equation (22)) of SCTLGSM pulsed beam as a function of propagation distance z for different chirp coefficients s = 0, −0.4ps−2, −0.8ps−2, −1.2ps−2.
Photonics 08 00102 g008
Figure 9. Density plots for the normalized spatiotemporal intensity I(x, y, t, z) (Equation (20)) of SCTLGSM pulsed beam in the xt plane at some propagation distances (a) z = 0, (b) z =20 m, (c) z = 30 m, (d) z = 40 m, (e) z = 50 m and for s = −0.4, n = 1, m = 1. (f) Normalized intensity I(x, y, t, z) (Equation (22)) of SCTLGSM pulsed beam as a function of propagation distance z for different chirp coefficients s = 0, −0.4 ps−2, −0.8 ps−2 and −1.2 ps−2.
Figure 9. Density plots for the normalized spatiotemporal intensity I(x, y, t, z) (Equation (20)) of SCTLGSM pulsed beam in the xt plane at some propagation distances (a) z = 0, (b) z =20 m, (c) z = 30 m, (d) z = 40 m, (e) z = 50 m and for s = −0.4, n = 1, m = 1. (f) Normalized intensity I(x, y, t, z) (Equation (22)) of SCTLGSM pulsed beam as a function of propagation distance z for different chirp coefficients s = 0, −0.4 ps−2, −0.8 ps−2 and −1.2 ps−2.
Photonics 08 00102 g009
Figure 10. Density plots for the normalized spatiotemporal intensity I(x, y, t, z) (Equation (20)) of SCTLGSM pulsed beam in the xt plane for different beam orders (a) m = n = 0, (b) m = n = 1, (c) m = n = 2, (d) m = n = 3 at z = 50 m with s = −0.4ps−2.
Figure 10. Density plots for the normalized spatiotemporal intensity I(x, y, t, z) (Equation (20)) of SCTLGSM pulsed beam in the xt plane for different beam orders (a) m = n = 0, (b) m = n = 1, (c) m = n = 2, (d) m = n = 3 at z = 50 m with s = −0.4ps−2.
Photonics 08 00102 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ding, C.; Korotkova, O.; Horoshko, D.; Zhao, Z.; Pan, L. Evolution of Spatiotemporal Intensity of Partially Coherent Pulsed Beams with Spatial Cosine-Gaussian and Temporal Laguerre–Gaussian Correlations in Still, Pure Water. Photonics 2021, 8, 102. https://doi.org/10.3390/photonics8040102

AMA Style

Ding C, Korotkova O, Horoshko D, Zhao Z, Pan L. Evolution of Spatiotemporal Intensity of Partially Coherent Pulsed Beams with Spatial Cosine-Gaussian and Temporal Laguerre–Gaussian Correlations in Still, Pure Water. Photonics. 2021; 8(4):102. https://doi.org/10.3390/photonics8040102

Chicago/Turabian Style

Ding, Chaoliang, Olga Korotkova, Dmitri Horoshko, Zhiguo Zhao, and Liuzhan Pan. 2021. "Evolution of Spatiotemporal Intensity of Partially Coherent Pulsed Beams with Spatial Cosine-Gaussian and Temporal Laguerre–Gaussian Correlations in Still, Pure Water" Photonics 8, no. 4: 102. https://doi.org/10.3390/photonics8040102

APA Style

Ding, C., Korotkova, O., Horoshko, D., Zhao, Z., & Pan, L. (2021). Evolution of Spatiotemporal Intensity of Partially Coherent Pulsed Beams with Spatial Cosine-Gaussian and Temporal Laguerre–Gaussian Correlations in Still, Pure Water. Photonics, 8(4), 102. https://doi.org/10.3390/photonics8040102

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