Self-action of continuous laser radiation and Pearcey diffraction in a water suspension with light-absorbing particles

.


Introduction
During the past decades, the specific behavior of small particles under the action of light has become a powerful and efficient tool in fundamental and applied optical research. Since the pioneering work by Ashkin et al [1][2][3], light-induced mechanical effects have been widely used for trapping, transportation, orientation and localization of micro-particles (see, e.g., Refs [4][5][6][7][8].), as well as for the study of fine structure of light fields [9][10][11][12][13][14][15][16][17]. Among various further applications, the use of ensembles of suspended optically-driven micro-and nanoparticles has recently been proposed as a means for the control and regulation of light propagation [18,19]. In absorptive media, propagation of inhomogeneous optical radiation induces non-uniform distributions of temperature and, consequently, of many associated optical constants. The corresponding effects are coupled with self-diffraction, self-focusing and emergence of shock waves. Usually, such self-action phenomena accompany the propagation of short powerful laser pulses [20][21][22][23] where, due to high density of the absorbed energy, the medium may be damaged. This is detrimental, for example, in the study of organic and biological objects.
Continuous laser radiation of moderate power can also produce a temperature gradient but normally it is insufficient to generate noticeable non-linear self-action of the optical field. Relatively large thermal non-linearity was observed in the ethanol solutions of the Chinese tea liquid [24][25][26] where the self-induced lens effect was attributed to strong absorption of the He-Ne laser radiation by molecules of chlorophyll and similar herbal dyes [27]. However, in such solutions the non-linear effects were only observable in considerably thick samples (~1 mm).
Recently, we have demonstrated the possibility of extremely high thermal inhomogeneity in a dielectric medium containing suspended absorbing particles [19]. In such systems, remarkable non-linear phenomena may occur within micron-size volumes in the sparing regime preserving integrity of the medium and its constituents. This fact can be used in micro-and nanotechnologies as well as in diagnostics and selective treatment of biological objects. The present work is devoted to the study and realization of the corresponding possibilities.
Usual ways for detection and measurement of weak nonlinearities are based on the wellknown Z scan method [28] in which the sample is placed within a focused Gaussian beam and takes different longitudinal positions near the focal point. Then the quantitative characteristics of the non-linearity are derived from the additional focusing (or defocusing) of the beam which takes place in the sample [23-25]. However, the "strength" of the self-induced inhomogeneity (e.g., the refractive index gradient) inversely depends on the near-focus beam size which cannot be arbitrary small, and this limits the method sensitivity. This limitation can be overcome if the incident field inhomogeneity appears not only due to focusing but also due to additional modulations, e.g. caused by two-wave interference. Emergence of fine periodic structures is especially favorable for the study of weak non-linear effects as it provokes the suitably observable self-diffraction phenomena [19].
In this work, we continue the exploration of the non-linear effects in a disperse medium under the action of a strongly inhomogeneous optical field with moderate power. In view of the above paragraph, the next step is to employ an incident field with a controllable speckle pattern, regular both in the transverse and longitudinal directions. Such situations occur in the coherent optical field formed after a collimated beam passes a thick cylindrical lens. Due to strong aberrations, this focused field is characterized by the special features [29,30]: (i) the intensity distribution is bounded by the caustic surface and (ii) the amplitude and phase distributions contain the well-developed and regular system of bright and dark spots accompanied by the phase singularities, known as the Pearcey diffraction pattern [31].
The methods and background of this work represent a further development of the Z scan method [24,25,28] and the approach employed for investigation of the self-diffraction processes of the two-beam interference field in the same disperse media [19]. The disperse medium (dielectric liquid + suspended particles) is initially characterized by a homogeneous distribution of absorbing particles. Exposure of this medium to an optical field with a spatially inhomogeneous intensity distribution results in inhomogeneous heating, which, in turn, makes the medium refraction index coordinate-dependent. If the disperse medium is situated close to the focal plane of the cylindrical lens, additional modulations of the refraction index in the medium due to the light absorption by the particles becomes strong enough to cause a perturbation of the propagating light beam. As a result, this will induce self-diffraction processes that manifest themselves via the intensity redistribution in the beam leaving the medium. We have performed computer simulation and experimental modeling of this effect for different positions of the sample with respect to the lens focus, and confronting the simulation and experimental data served to measure the non-linear characteristics of the medium.

The Pearcey diffraction pattern
The geometrical conditions of the problem are illustrated by Fig. 1. The action of the cylindrical lens is considered as a passage through a phase transparency with a transmission function is the distance from a point in the lens output plane (P) with the transverse coordinate x to a point in the observation plane (C) with the transverse coordinate ξ, z is the distance between the aperture plane and the observation plane (plane z = 0 coincides with the flat input face of the lens as shown in Fig. 1). The input field is assumed to be a plane wave with unit amplitude, then the aperture function ( )  In contrast to the Kirchhoff integral [32,33], Eq. (1) is applicable to field calculations at small distances between the pupil plane (z = r) and the observation plane. Figure 2 shows the calculated intensity distribution ( , ) I z ξ in the region z > r for a lens aperture 20 m μ and the cylindrical surface radius 10 m r = μ together with the results experimentally observed where the cylindrical lens was formed by a segment of an optical fiber with appropriate radius and refraction index 1.528 m = ; the radiation wavelength is λ = 445 nm. Because of the high spherical aberrations of the lens whose action is illustrated in Fig. 2, its geometric-optics focus does not coincide with the point of maximal energy concentration. Due to the better experimental localizability, the latter is a suitable reference point and from now on we will, conventionally, apply the term 'focal plane' to the transverse plane containing the single axial brightest spot. According to Fig. 2a, this plane corresponds to 28 m z = μ . Figure 2b represents a fragment of the intensity distribution observed in the near-focal region with the help of an oblique photodetector facilitating the visualization of the longitudinal intensity distribution. It clearly reproduces not only caustics but many minor details of the calculated Pearcey pattern.

Influence of the disperse medium near the lens focus
The presence of the cell C with the dispersed medium leads to modification of the standard Pearcey pattern of Fig. 2 in any plane located behind the cell (see Fig. 1). Let the cell be placed in a certain plane z near the focus and the plane of observation be situated at the distance Z from the cylindrical lens, Z > z. In this plane, the field amplitude distribution ( ) , , U z Z ξ ′ distorted due to interaction with the disperse medium can be calculated with the help of the modified diffraction integral (1), , absorption coefficient and d is the thickness of the disperse medium in the cell. In Eq. (2) we suppose the medium layer to be thin enough so that it can be treated as a thin phase transparency with complex amplitude transmission function ( )

Table 1. Values of the disperse medium parameters used in the calculation of the radiation-induced refraction index modification (Eqs. (3) and (4)).
1 Mean size of the pigment particles 0.2 μm 2 Concentration of the pigment particles 0.4 g·cm -3 3 Complex refraction index of the particle coating 1.82 + 0.74i 4 Refraction index of the matrix medium (water) 1.33 5 Absorption coefficient of the medium 7·10 2 cm -1 6 Cell thickness 10 μm (αd = 0.7) 7 Radiation wavelength 0.445 μm (k = 2π/λ = 1.41·10 5 cm -1 ) Specific heat capacity (water) 4.2 kJ⋅kg -1 ⋅K -1 10 Mass density (water) 10 3 kg·m -3 The absorption coefficient α is an effective medium characteristic determined by averaging the light absorption in individual particles, taking into account their distribution and concentration. The pigment particles used in our experiments are close to spherical with mean radius 0.1 m a = μ and contain ~20 nm absorptive coating of carbon-based resin. The complex refraction index of the coating is assumed to be close to that of the atmospheric soot [34]. The parameters of the particles and of the disperse medium used in the experiment are specified in Table 1.
For evaluation of Eq. (2), the key quantity is ( , ) n z ξ Δ . In the first approximation we assume that it is determined completely by the incident field, and that the prevailing mechanism for the refraction index modulation is associated with inhomogeneous heating of the medium due to the light energy absorption [19]. In view of the small thickness of the medium layer compared with the longitudinal scale of the near-focus Pearcey pattern, all points of the medium with the same transverse coordinates experience practically the same field-induced action. In such conditions, the density of the additional energy accumulated in where ( ) dn dT is the refractive index temperature coefficient, C is the specific heat capacity (per unit mass of the medium) and ρ is the medium mass density; for simplicity, here we suppose n Δ to be proportional to q thereby neglecting possible distortion of the absorbed energy distribution due to the heat conductivity. Combining Eqs. (3) and (4) determines the spatial distribution of the induced medium inhomogeneity. According to these equations, the field-induced non-uniform heating only affects the real refraction index thus forming a phase grating. Additionally, due to the gradient optical force [2,9,35] and the photophoresis [7,8,36,37], the suspended particles may concentrate near the intensity minima which contributes to the creation of an amplitude inhomogeneity. However, direct estimates [19] show that at moderate light power, compatible with the medium stability (~40 mW to prevent the water boiling and cavitation effects), variation of the particles' concentration does not exceed 5% and in the actual conditions, influence of the amplitude inhomogeneity is negligible. This permits us to restrict the analysis to the case of a homogeneous particle distribution.
A schematic drawing of the experimental setup is presented in Fig. 3. As a radiation source, the semiconductor laser XJ-A140 was used with wavelength λ = 445 nm and controllable power ≤ 0.3 W. The laser radiation is linearly polarized in the plane orthogonal to the figure plane, which coincides with the focusing plane of the cylindrical lens CL and justifies the scalar approximation of Eqs. (1) and (2). Telescope T consisting of two microobjectives with a pinhole diaphragm in the common focus forms a parallel beam of 10 mm in diameter. Further, the beam is focused by the lens CL onto the silica cell C containing the medium (water) with suspended particles. The lens is made of glass with refraction index m = 1.52 in the shape of a cylindrical segment enclosed between the plane input face and the cylindrical surface of radius r = 1.8 mm; the lens size along the vertical (x) direction equals 3.6 mm and along the y-direction 8 mm. The use of the macroscopic arrangement (in contrast to the micron-size lens discussed in Section 2) is not of principal importance and was motivated by convenience of the cell manipulations. However, in this arrangement we could not measure the 3D diffraction pattern directly but the image of Fig. 2b can still be applied in view of the spatial similarity laws regulating the Pearcey pattern [29]. The lens focus defined as the position of the highest light concentration near the system axis (see Section 2) is located at a distance of f = 5 mm from the lens input plane; the plane of observation is situated at Z = 15 mm (10 mm behind the focus). The cell walls are of thickness 0.5 mm and enclose a 10 m μ layer of the disperse medium. During normal working conditions, the light power reaching the cell is approximately 50 mW. The CCD camera registers the radiation that has passed the cell. Its spatial distribution shows interesting peculiarities that clearly demonstrate the presence of self-diffraction effects (Fig. 4). Here, results of the simulation for different cell positions are compared with the experiment. The cell position is specified by the distance z Δ between the cell central plane and the lens focal plane (see Fig. 3) and varies from 350 m − μ to 50 m μ with 50 m μ steps; each cell position is characterized by a group of three images numbered from (a) to (j) (see Table 2). Table 2. Correspondence between the cell shift with respect to the lens focus and the label of the three-image groups in Fig. 4.  ) and (4), is responsible for the induced refraction index inhomogeneity. Corresponding distribution of n Δ is determined for the parameters of the disperse medium given in Table 1; afterwards, it is substituted into Eq. (2) and the distribution of the field incident onto the CCD camera is simulated (images (ii) in Fig. 4). Patterns marked (iii) show the experimental intensity distributions observed by the CCD camera. It is seen that the calculated profiles (ii) correlate adequately with the experimental distributions (iii) along the axis ' ξ .
During the calculations, values of the coefficient η in Eq. (3) were determined via the fitting procedure aimed at best matching the calculated plots (ii) with the corresponding sections of the measured intensity patterns (iii) of Figs. 4(b)-(h). An example of such a matching is provided by Fig. 5 where the detailed comparison between the simulation data of Fig. 4(e) (ii) and the experimental profile of the image of Fig. 4(e) (iii), taken along y = 0, is illustrated. The results obtained for η due to the best fitting procedure are presented in the second row of Table 2. Fig. 4. Results of the numerical calculations and experimental observations of the focused beam that has passed through the disperse medium situated at different distances from the lens focus ∆z; (і) calculated intensity profile in the near-focus region; (іі) calculated intensity profile behind the focus in the input plane of the CCD camera; (ііі) experimentally observed intensity distributions. Groups of images (a) -(j) correspond to different ∆z (see Table 2). maximal near the axis and gradually decreases to the vertical edges. However, in the patterns observed behind the focus, in many cases, the intensity maxima occur near the edges (Fig. 4  b-g, iii). This is explained by the self-diffraction effects due to which the near-axis energy is redistributed between several maxima along the x-axis. Near the focal-line edges, the intensity is rather small and the radiation propagates in the linear regime with preserving single maximum in the x-direction. Note that the observed intensity distributions in the (xy) plane resemble the Pearcey diffraction patterns in the (xz) plane (see Fig. 2). During the experiment, the laser power approximately of 100 mW (~50 mW in the cell) was maintained; with lowering the power, non-linear self-diffraction effects were rapidly suppressed, while increasing the laser power above this value led to emergence of shock waves, acoustic noise generation and eventually to the medium being destructed.

Conclusion
Following to the way outlined before [19], we have considered additional examples of possible applications of the disperse media containing suspended absorptive micro-and nanoparticles for creation of thermo-induced controllable optical switches and regulators. Carboncoated particles suspended in water efficiently absorb and accumulate the light energy and enable articulate non-linear behavior even being exposed to comparatively weak CW laser radiation with moderate power.
In particular, in this work we present some preliminary results of theoretical and experimental studies of the self-induced laser beam transformations occurring near the focus of a cylindrical lens giving rise to the Pearcey diffraction pattern [31]. This technique can be considered as a generalization of the known Z scan method [24,25,28] enabling increase of sensitivity due to strong spatial inhomogeneity of the incident radiation. It was demonstrated that the self-focusing effect in the absorbing disperse medium takes place because of selfdiffraction of the incident radiation by the phase inhomogeneity caused by the high-gradient Pearcey pattern. Due to the high intensity gradient, the self-diffraction effect occurs in very thin (micron-sized) layers; along with the small consumed power this makes the described scheme promising for integrated optical circuits.
The results obtained as well as the methods developed in this work open up new prospects in detection and measurement of weak optical nonlinearities. They testify for possibilities to observe non-linear optical effects at relatively low energy densities, i.e. in sparing conditions (note that according to the data of Tables 1 and 2, noticeable non-linear effects may occur at quite modest local heating not exceeding 10 K). This may be useful for the study of "delicate" objects, especially for biological investigations dealing with diagnostics of cells, regulation of their activity and manipulation based on the electromagnetic field interaction with matter. The media based on suspensions of strongly absorptive nanoparticles dispersed within a dielectric liquid offer additional advantages associated with the feasibility of preparation and control of the prescribed non-linear properties.