2,022
Views
15
CrossRef citations to date
0
Altmetric
Review Articles

Optical spatial shock waves in nonlocal nonlinear media

, , , , &
Article: 1662733 | Received 28 May 2019, Accepted 27 Aug 2019, Published online: 16 Oct 2019

Figures & data

Figure 1. Phase chirp u(x) (a,b), and amplitude ρ(x,z) (c,d) for transverse dimensionality D=1 and different values of z, as indicated. (a,c) are obtained by Eqs. (18) with ϵ=103. (b,d) are simulations of the result of the system (17) with D=1,α=0,χ=1,σ2=5.

Figure reprinted with permission from [Citation29]. Copyright 2007 by the American Physical Society.

Figure 1. Phase chirp u(x) (a,b), and amplitude ρ(x,z) (c,d) for transverse dimensionality D=1 and different values of z, as indicated. (a,c) are obtained by Eqs. (18) with ϵ=10−3. (b,d) are simulations of the result of the system (17) with D=1,α=0,χ=−1,σ2=5.Figure reprinted with permission from [Citation29]. Copyright 2007 by the American Physical Society.

Figure 2. Pictorial representation of an energy landscape. When the system is in the proximity of a local maximum it obeys the RHO Hamiltonian, in figure HˆRO. In the proximity of the minimum the system obeys the Hamiltonian of a harmonic oscillator, in figure HˆHO. The two Hamiltonians are explicitly written in the two the corresponding text boxes, with the related dynamical systems and the discrete eingenvalues. Insets show the transverse profiles of the respective eigenfunctions, bounded on right-hand side for the harmonic oscillator, unbounded on the left-hand side for the RHO.

Reprinted by permission from Macmillan Publishers Ltd. from [Citation38]. Copyright 2015.

Figure 2. Pictorial representation of an energy landscape. When the system is in the proximity of a local maximum it obeys the RHO Hamiltonian, in figure HˆRO. In the proximity of the minimum the system obeys the Hamiltonian of a harmonic oscillator, in figure HˆHO. The two Hamiltonians are explicitly written in the two the corresponding text boxes, with the related dynamical systems and the discrete eingenvalues. Insets show the transverse profiles of the respective eigenfunctions, bounded on right-hand side for the harmonic oscillator, unbounded on the left-hand side for the RHO.Reprinted by permission from Macmillan Publishers Ltd. from [Citation38]. Copyright 2015.

Figure 3. (a) ||fn(x)||2 in Eq. (23) for increasing even order n; (b) corresponding phase chirps xArgfn(x); (c) weights pn(0) [Eq. (27)] of the GV expansion of a Gaussian wave packet.

reprinted with permission from [Citation37]. Copyright 2015 by the American Physical Society.

Figure 3. (a) ||fn−(x)||2 in Eq. (23) for increasing even order n; (b) corresponding phase chirps ∂xArgfn−(x); (c) weights pn(0) [Eq. (27)] of the GV expansion of a Gaussian wave packet.Figure reprinted with permission from [Citation37]. Copyright 2015 by the American Physical Society.

Figure 4. (a) Numerical solution of Eq. (19) with p=104 and σ2=10; (b) projection on GVs for increasing order n for α=0.3 and γ=8; continuous lines are from Eq. (19), dots are from Eq. (27); (c) as in panel (b) for γ=24.

Figure reprinted with permission from [Citation37]. Copyright 2015 by the American Physical Society.

Figure 4. (a) Numerical solution of Eq. (19) with p=104 and σ2=10; (b) projection on GVs for increasing order n for α=0.3 and γ=8; continuous lines are from Eq. (19), dots are from Eq. (27); (c) as in panel (b) for γ=24.Figure reprinted with permission from [Citation37]. Copyright 2015 by the American Physical Society.

Figure 5. Experimental transverse intensity profiles of an initial Gaussian beam propagating in a thermal medium. Measurements are performed for varying input power P=πW02I0. Insets show the 2D output patterns.

Figure reprinted with permission from [Citation29]. Copyright 2007 by the American Physical Society.

Figure 5. Experimental transverse intensity profiles of an initial Gaussian beam propagating in a thermal medium. Measurements are performed for varying input power P=πW02I0. Insets show the 2D output patterns.Figure reprinted with permission from [Citation29]. Copyright 2007 by the American Physical Society.

Figure 6. (a) Experimental setup. Authors of [Citation38,Citation40] collected the transmitted and fluorescence images of the laser beam propagating in RhB samples. Two types of launching lenses L1 were used: a cylindrical and a spherical, for the 1D and 2D experiments, respectively. The top fluorescence image of the propagating beam was collected by a microscope placed above the RhB samples. The second lens is spherical and was used to collect the transverse output profile. (b, c) Top-view intensity distribution as obtained from 2D experiment (b) and numerical simulations (c). Respectively experimental (d) and numerical (e) sections of the images (b) and (c) taken at z=0.2 (red), 0.6 (green) and 0.9 mm (blue).

Reprinted by permission from Macmillan Publishers Ltd. from [Citation38]. Copyright 2015.

Figure 6. (a) Experimental setup. Authors of [Citation38,Citation40] collected the transmitted and fluorescence images of the laser beam propagating in RhB samples. Two types of launching lenses L1 were used: a cylindrical and a spherical, for the 1D and 2D experiments, respectively. The top fluorescence image of the propagating beam was collected by a microscope placed above the RhB samples. The second lens is spherical and was used to collect the transverse output profile. (b, c) Top-view intensity distribution as obtained from 2D experiment (b) and numerical simulations (c). Respectively experimental (d) and numerical (e) sections of the images (b) and (c) taken at z=0.2 (red), 0.6 (green) and 0.9 mm (blue).Reprinted by permission from Macmillan Publishers Ltd. from [Citation38]. Copyright 2015.

Figure 7. (a) Observed intensity decay at different laser powers, obtained by slicing along X0.1mm the top-view intensity distribution the propagation direction (see the yellow line in ). (b) Numerically calculated decays in the conditions of panel A. (c) Peak region of the experimental curve at P=450mW. The superposition of the first two exponential decays unveils the presence of two GVs, the fundamental state, n=0 (slowly decaying) and the first excited state, n=2 (fastly decaying). (d) Decay rates vs P for the fundamental state, Γ0 (filled circles) and the excited state, Γ2, (triangles).

Reprinted by permission from Macmillan Publishers Ltd. from [Citation38]. Copyright 2015.

Figure 7. (a) Observed intensity decay at different laser powers, obtained by slicing along X≃0.1mm the top-view intensity distribution the propagation direction (see the yellow line in Figure 6b). (b) Numerically calculated decays in the conditions of panel A. (c) Peak region of the experimental curve at P=450mW. The superposition of the first two exponential decays unveils the presence of two GVs, the fundamental state, n=0 (slowly decaying) and the first excited state, n=2 (fastly decaying). (d) Decay rates vs P for the fundamental state, Γ0 (filled circles) and the excited state, Γ2, (triangles).Reprinted by permission from Macmillan Publishers Ltd. from [Citation38]. Copyright 2015.

Figure 8. (a, b) CCD image of the light beam at laser powers P=2W and 4W, respectively; the bottom panels show the normalized intensity profile at the maximum waist along Y=0. (c) Analytical solution obtained by Eq. (37) changing Gaussianly the power P in the y direction. (d) As in (c) but for higher powers; the bottom panels show the slice of panel (c) and (d) at y=0, i.e., Eq. (37) square modulus for W=1.5 and γ12 and γ40, respectively. (e) Log-scale normalized intensity as a function of power, as obtained by slicing along Y a region in panel (b). The slopes of the straight lines give the GV decay rates (γ1=8±0.4 and γ2=1.6±0.1). Their quantized ratio is 5.0±0.4 as expected from theory [Citation37]. (f) Intensity oscillations for different power values. (g) Measured oscillations period T as a function of power; continuous line is the fit function TP4, as expected by the theory; the inset shows the same curve of (g) with P1/4 as abscissa axis.

Reprinted with permission from [Citation40]. Copyright 2016 Optical Society of America.

Figure 8. (a, b) CCD image of the light beam at laser powers P=2W and 4W, respectively; the bottom panels show the normalized intensity profile at the maximum waist along Y=0. (c) Analytical solution obtained by Eq. (37) changing Gaussianly the power P in the y direction. (d) As in (c) but for higher powers; the bottom panels show the slice of panel (c) and (d) at y=0, i.e., Eq. (37) square modulus for W=1.5 and γ≃12 and γ≃40, respectively. (e) Log-scale normalized intensity as a function of power, as obtained by slicing along Y a region in panel (b). The slopes of the straight lines give the GV decay rates (γ1=−8±0.4 and γ2=−1.6±0.1). Their quantized ratio is 5.0±0.4 as expected from theory [Citation37]. (f) Intensity oscillations for different power values. (g) Measured oscillations period T as a function of power; continuous line is the fit function T∝P4, as expected by the theory; the inset shows the same curve of (g) with P1/4 as abscissa axis.Reprinted with permission from [Citation40]. Copyright 2016 Optical Society of America.

Figure 9. (a, b) Trajectories x(z) and (b, d) phase space (x,v), respectively, with disorder strength η=0 and η=0.1. z varies from z=0 to z=3.

Figure reprinted with permission from [Citation32]. Copyright 2012 by the American Physical Society.

Figure 9. (a, b) Trajectories x(z) and (b, d) phase space (x,v), respectively, with disorder strength η=0 and η=0.1. z varies from z=0 to z=3.Figure reprinted with permission from [Citation32]. Copyright 2012 by the American Physical Society.

Figure 10. Transverse intensity patterns for different input power P and silica spheres concentration c: (a) P=5mW, c=0w/w, (b)P=400mW, c=0w/w, (c) P=5mW, c=0.017w/w, (d) P=400mW, c=0.017w/w, (e) P=5mW, c=0.030w/w, (f) P=400mW, c=0.030w/w. White 1D curves show the measured section of the intensity profiles vs X.

Figure reprinted with permission from [Citation32]. Copyright 2012 by the American Physical Society.

Figure 10. Transverse intensity patterns for different input power P and silica spheres concentration c: (a) P=5mW, c=0w/w, (b)P=400mW, c=0w/w, (c) P=5mW, c=0.017w/w, (d) P=400mW, c=0.017w/w, (e) P=5mW, c=0.030w/w, (f) P=400mW, c=0.030w/w. White 1D curves show the measured section of the intensity profiles vs X.Figure reprinted with permission from [Citation32]. Copyright 2012 by the American Physical Society.

Figure 11. Far-field intensity profiles at the output of the silica aerogel for Pin ranging from 1mW to 1W, and input beam waist w0 ranging from 43μm to 1.4 mm. Images in the second and third rows correspond to the same incident laser power and beam size, but different positions of the incident laser beam.

Reprinted with permission from [Citation35]. Copyright 2014 Optical Society of America.

Figure 11. Far-field intensity profiles at the output of the silica aerogel for Pin ranging from 1mW to 1W, and input beam waist w0 ranging from 43μm to 1.4 mm. Images in the second and third rows correspond to the same incident laser power and beam size, but different positions of the incident laser beam.Reprinted with permission from [Citation35]. Copyright 2014 Optical Society of America.

Figure 12. Output beam intensity transverse profiles, coming out from a 2mm long M-Cresol/Nylon solution. Input power varies: (a) Pin=2μW, (b) Pin=5mW, (c) Pin=10mW, (d) Pin=20mW.

Reprinted with permission from [Citation36]. Copyright 2014 Optical Society of America.

Figure 12. Output beam intensity transverse profiles, coming out from a 2mm long M-Cresol/Nylon solution. Input power varies: (a) Pin=2μW, (b) Pin=5mW, (c) Pin=10mW, (d) Pin=20mW.Reprinted with permission from [Citation36]. Copyright 2014 Optical Society of America.

Figure 13. Output beam waist for varying hemoglobin concentration and input power. (a) Detected beam diameter as function of input power through the hemoglobin solutions for four different concentrations (Hb1-Hb4): 2.4, 5.1, 8.6, and 15.0 million cells per mL. Nonlinear self-focusing of the beam occurs around 100mW for high concentrations of hemoglobin, but it subsequently expands into thermal defocusing rings at high powers. (b-e) Output beam transverse intensity profiles for (b) self-trapped beam at high concentration and low power, (c) DSW at high concentration and high power, (d) self-trapped beam at low concentration and low power, (e) DSW at low concentration and high power.

Reprinted by permission from Macmillan Publishers Ltd. from [Citation43]. Copyright 2019.

Figure 13. Output beam waist for varying hemoglobin concentration and input power. (a) Detected beam diameter as function of input power through the hemoglobin solutions for four different concentrations (Hb1-Hb4): 2.4, 5.1, 8.6, and 15.0 million cells per mL. Nonlinear self-focusing of the beam occurs around 100mW for high concentrations of hemoglobin, but it subsequently expands into thermal defocusing rings at high powers. (b-e) Output beam transverse intensity profiles for (b) self-trapped beam at high concentration and low power, (c) DSW at high concentration and high power, (d) self-trapped beam at low concentration and low power, (e) DSW at low concentration and high power.Reprinted by permission from Macmillan Publishers Ltd. from [Citation43]. Copyright 2019.

Figure 14. Experimental setup. A CW laser beam, λ=532nm, is sent through a 4-f telescope. A ground-glass plate G, placed in the midst of the telescope (on the focus of the first lens), generates a speckle pattern. The incoherent beam impinges the samples, a cylindrical tube filled with a solution of methanol and graphene nanoscale flakes, with waist W0=2.3mm, while the initial coherence length λc0 is controlled by changing the beam size on G. A CDD camera detects the output.

Reprinted by permission from Macmillan Publishers Ltd. from [Citation39]. Copyright 2015.

Figure 14. Experimental setup. A CW laser beam, λ=532nm, is sent through a 4-f telescope. A ground-glass plate G, placed in the midst of the telescope (on the focus of the first lens), generates a speckle pattern. The incoherent beam impinges the samples, a cylindrical tube filled with a solution of methanol and graphene nanoscale flakes, with waist W0=2.3mm, while the initial coherence length λc0 is controlled by changing the beam size on G. A CDD camera detects the output.Reprinted by permission from Macmillan Publishers Ltd. from [Citation39]. Copyright 2015.

Figure 15. Experimental observation at short-range regime. (a, b) Experimental beam profiles of the output intensity recorded at power a P=0.05W, (b) P=2.50W. (c, d) Zooms on details of (b) that evidence the development of shocklets. (e, f) Numerical simulations of NLSE equation, and (g, h) corresponding zooms.

Reprinted by permission from Macmillan Publishers Ltd. from [Citation39]. Copyright 2015.

Figure 15. Experimental observation at short-range regime. (a, b) Experimental beam profiles of the output intensity recorded at power a P=0.05W, (b) P=2.50W. (c, d) Zooms on details of (b) that evidence the development of shocklets. (e, f) Numerical simulations of NLSE equation, and (g, h) corresponding zooms.Reprinted by permission from Macmillan Publishers Ltd. from [Citation39]. Copyright 2015.

Figure 16. Experimental observation at long-range regime. (a-c) One-dimensional intensity transverse profiles along x at y=0. (d-f) Two-dimensional intensity profiles. The asymmetry in the lower part of the beam is due to convection within the sample. (g-i) Numerical simulations of NLSE. (j-l) Experimental and (m-o) numerical spectrograms: the Z-shaped distortion reveals a dramatic coherence degradation on the annular boundaries of the beam (the coherence length decreases at the shock front), while a significant coherence enhancement occurs in the internal region of the beam. Input beam power: (a, d, g, j, m) P=0.05W, (b, e, h, k, n) P=1.25W, (c, f, i, l, o) P=2.50W.

Reprinted by permission from Macmillan Publishers Ltd. from [Citation39]. Copyright 2015.

Figure 16. Experimental observation at long-range regime. (a-c) One-dimensional intensity transverse profiles along x at y=0. (d-f) Two-dimensional intensity profiles. The asymmetry in the lower part of the beam is due to convection within the sample. (g-i) Numerical simulations of NLSE. (j-l) Experimental and (m-o) numerical spectrograms: the Z-shaped distortion reveals a dramatic coherence degradation on the annular boundaries of the beam (the coherence length decreases at the shock front), while a significant coherence enhancement occurs in the internal region of the beam. Input beam power: (a, d, g, j, m) P=0.05W, (b, e, h, k, n) P=1.25W, (c, f, i, l, o) P=2.50W.Reprinted by permission from Macmillan Publishers Ltd. from [Citation39]. Copyright 2015.