# Electrostrictive Mechanism of Radiation Self-Action in Nanofluids

Authors: Albert Livashvili; Victor Krishtop; Margarita Yakunina
Journal: Advances in Condensed Matter Physics
 (2013)

---

## Abstract

The electrostriction mechanism of beam self-focusing in nanofluids is theoretically investigated. An analytical solution of the diffusion equation, which describes the dynamics of particles in nanofluids, was obtained and studied. Explicit expressions for the nonlinear part of the refractive index and concentration lens focal length are presented. It is shown that there is a limit on the radiation intensity associated with the physical and hydrodynamic characteristics of the phenomena in these processes.


---

## Body

## 1. Introduction

Nanoparticles are being increasingly adopted in new and different areas of science and technology [1–4]. Colloidal suspensions that contain nanoparticles, also known as nanofluids, have found a variety of important applications in modern technologies [5]. For example, magnetic fluids are largely used for polishing optical components [6, 7], suspensions of silica particles in liquid crystals have exhibited extraordinary capabilities for optical storage [8, 9], and artificial media with a high optical nonlinearity were obtained in colloids of submicrometer-size particles in liquids [10–13]. As shown in recent studies [14–21], the liquid-phase environment of dispersed nanoparticles of wide bandgap semiconductors or insulators is very effective for a number of nonlinear optical effects. But the physical mechanisms involved, particularly with nonlinear optical processes in such media, are not entirely clear and require further study.When placed in an electromagnetic field, particles in a microheterogeneous medium, with components that have different refractive indices, are subject to electrostrictive forces that could lead to the appearance of concentration flows [21–23]. Depending on the sign of the polarizability, these microparticles can be pulled in or pushed out of high electric field areas. At the same time that this electrostriction flow phenomenon occurs, thermal diffusion due to the temperature gradient [23] also occurs, but we do not consider thermal effects here. The concentration dynamics of colloidal particles in a periodic light field have been theoretically investigated in [24, 25], and similar problems within the small-perturbation approximation of the concentration were studied in [26, 27]. In this paper, the restrictions of this approximation are removed.Here, using the exact solution of the transformed diffusion equation for particles in the field of a light wave, we explore the dynamics of the particle concentration. The results allow us to find an explicit expression for the nonlinear part of the refractive index and calculate the focal length of the lens concentration formed by the self-interaction effects of the light field.


## 2. Diffusion Equation

The geometry of the problem considered in the paper is as follows. Electromagnetic radiation with some profile of the intensity falls on the cylindrical cell containing the investigated medium (Figure1).Figure 1
The geometry of the problem. 1: cell; 2: liquid medium with nanoparticles.The equation describing the spatial and temporal variations of the nanoparticle mass concentrationC(r) can be written as a diffusion equation in terms of the external field gradient (electrostrictive) forces as follows (we consider the axisymmetric case) [24]:
(1)∂C∂t=D(∂2C∂r2+1r∂C∂r)-div(γC∂I∂r),
where C(r,t)=m0/m, m0 is the total nanoparticle mass, m is the mass of the whole fluid sample, D is the diffusion coefficient, γ=4πβD/c-nkBT, β is the polarizability of the particles, kB is the Boltzmann constant, T is the temperature of the medium, c- is the speed of light in vacuum, n is the effective refractive index of the medium, I=I~0exp(-r2/r02) is the intensity of the light beam, and r0 is the beam radius. Equation (1) is sometimes referred to as the Smoluchowski equation [28].To obtain the exact solution, consider the axial profile of the intensity of the beam in the paraxial Gaussian approximationI=I~0(1-r2/ro2). Then, as it is small, we drop the term proportional to ∇C∇I (comparing ∇I∇C≈δCδI/l02 and C∇2I≈C(δI/l02), where δC and δI are small increments of concentration and intensity, resp., and l0 is a characteristic length. We find that δC<<C, which confirms the validity of the approximation). This is equivalent to neglecting convection, and in some cases, the experiment can be carried out in such a way that the influence convection is minimal. Equation (1), after the relevant changes for the initial and boundary conditions, can then be written as
(2)∂C∂t=D(∂2C∂r2+1r∂C∂r)+φC,(3)C(r,0)=C0,∂C∂r|r=0=0,C(r0,t)=C0,
where φ=4I~0γ/r02 and C0 is the initial particle concentration.The second boundary condition expresses the fact that the functionC(r,t) is finite at the origin and its symmetry axis of the cell. It means that particle flux is constant at the origin. The third boundary condition is due to the fact that, in the paraxial approximation, the concentration perturbations almost do not reach to the boundary.We seek a solution in the form(4)C(r,t)=C~(r,t)exp(φt),
which leads to
(5)∂C~∂t=D(∂2C~∂r2+1r∂C~∂r),C~(r,0)=C0,∂C~∂r|r=0=0,C~(r0,t)=C0exp(-φt).By using the Green’s function of this problem [29], taking into account (4), and carrying out the relevant integration, the exact solution can be written as
(6)C(r,t)=2C0∑n=1∞J0(μnr/r0)J1(μn)μnexp(φ-Dμn2/r02)t+2C0∑n=1∞J0(μnr/r0)J1(μn)μn2(1-φr02/Dμn2)×[(1-exp(φ-Dμn2/r02)t)],
where J0(x) is the Bessel function of order zero and μn are the roots of the equation J0(μ)=0. From the resulting solutions, we can see that, in the case of a positive exponent, C(r,t)→∞ if t→∞, which is physically absurd. A negative sign can be provided in two ways: (1)φ<Dμn2/r02 or (2)φ<0. In the first case, using the explicit form of φ, we obtain a limit on the light intensity I~0<μn2c-nkB/16πβ (we also give an estimate below). In the second case, we must proceed from the explicit form of the polarizability, which for a sphere of radius a0 can be written as [30]
(7)β=εpεp-εlεp+2εla03,
where εp and εl are the dielectric constant of the particles and liquid, respectively. Thus, case (2) is realized when εp<εl. Consequently, the transition to stationary conditions (t→∞) in (6) is possible and can be represented by
(8)Cst(r)=2C0∑n=1∞J0(μnr/r0)J1(μn)(1-φr02/Dμn2)=C0J0(rφ/D)J0(r0φ/D),ifφ>0,(9)Cst(r)=C0I0(rφ/D)I0(r0φ/D),ifφ<0,
where I0(x) is the modified Bessel function of order zero. A series summation has been used to obtain (8) and (9) [31].Note that the vanishing of the denominator in the second term does not lead to divergence as the uncertainty can be easily removed by means of l’Hôpital’s rule. With regard to the convergence of (6) (omitting the details of the analysis), it can be argued that convergence is ensured by the exponents and the condition r/r0<1. Furthermore, for small density perturbations, corresponding expressions for (8) and (9) were found in [26].


## 3. Concentration Optical Nonlinearity and Effects of Self-Action

If we take into account the dependence of the refractive index on the concentration, we can write(10)n(C)≈n(C0)+∂n∂C(C(r,t)-C0).Here,n(C0)=n0 is the unperturbed (linear part) of the refractive index and ∂n/∂C is a constant of the medium. Thus, referring to the above results for the concentration, we can say that the value of n(C) is a result of the nonstationary, nonlinear function of the radiation intensity, which is not observed in the medium with an instant response. This is characteristic of media with a nonlocal (slow) nonlinearity mechanism (in our case, electrostriction).Based on the inequalityCst(r=0)<1 and using (8) or (9), we can obtain a limit on the light intensity, but it is then necessary to solve the irrational inequality C0<J0(r0φ/D). By using the asymptotic behavior of the Bessel functions for small arguments, we get C0<1-r02φ2/D2. The possibility of such representation implies the following condition: r02φ2/D<1. Therefore, we first estimate the argument of the Bessel function by assuming β=0.1a03, a0=2×10-9 m, n=1.3, and T=300 K, and we obtain 4πI~0β/c-nkBT≈5×10-6I~0. That is, up to an intensity of I~0≈1010 W/m2, the argument of the Bessel functions is less than one and a decomposition is possible. Carrying out this expansion for (8) and (9) and substituting the result into (10), we obtain
(11)n(C)≈n(C0)+0.4πa03c-nkBT(1-r2r02)C0I~0∂n∂C,ifΦ>0,(12)n(C)≈n(C0)-0.4πa03c-nkBT(1-r2r02)C0I~0∂n∂C,ifΦ<0.Thus, within the paraxial approximation, the sign of the nonlinear part of the refractive index depends on the sign of the nanoparticle polarizability, that is, the values of the variablesεp and εl. Physical processes are related to the sign of the polarizability clearly presented in [15].Note that the limiting intensity can be obtained by solving a simple hydrodynamic equation describing the motion of a particle subject to two forces: the gradient force exerted by the electric field of the light wavesF∇=4πβ∇I/c-n and the Stokes drag force Fs=6πηa0v, where η is the viscosity of the fluid and ν is the velocity of the particles. Thus, by determining the maximum speed of the particles (we omit the details of the calculation) and restricting the Reynolds number to Re ≪ 1, we get I~0≪15(e4c-n/a02ρ)η2, where ρ is the particle density. Setting η=1×10-3 kg/m·s (water), we get the desired estimate of I-0≪1016 W/m2, which is consistent with our results. Using (11), we estimate the nonlinear part of the refractive index n2 in the case of r=0 as
(13)n2=0.4πa03c-nkBTC0.Takinga0=2.34×10-7 m, C0=10-3, and ∂n/∂C=10-2 gives n2≈2×10-13 m2/W. Note that, in [32], a value of n2≈3×10-13 m2/W which is consistent with our calculations was obtained experimentally for an aqueous suspension of latex particles of the same radius.Self-action, in particular, is the formation of a lens in a layer of the medium. For small particulate concentrations, the effective index of refraction of the medium is proportional to the particle concentration and characterized by the parameter(∂n/∂C).These results allow us to calculate the focal length of the lens concentration formed in a layer of thicknessd as follows:
(14)1F=-d(∂n∂C)(∂2C∂r2),forr=0.Then, as above, except in (8) and (9), where the arguments of the Bessel functions were quite small, after expanding in powers of the ratio Φ/D and the use of (13), we obtain an explicit expression for the focal length of
(15)1F=±φ2DdC0(1-r02φD)(∂n∂C).In this case, the positive sign corresponds to the case withΦ>0, that is, εp>εl, and the negative sign to Φ<0, that is, εp<εl.Let us rewrite (15) in a more convenient form for estimation as follows:
(16)F=r04c-nkBT0.2a03dPC0(∂n∂C).Assuming here thatr0=10-4 m, C0=10-2, n=1.3, d=10-2 m, and P=10-2 W (along with the other quantities above), we then get F≈10-4 m, and so the focus of the concentration lens will be inside the medium.As is well known, there is a discrete spectrum of nonlinear modes in the propagating beam, each of which takes its critical power from the self-focused beam. For example, the lowest axially symmetric mode has a critical power of [33]
(17)Pcr=1.3×λ232π2n2.Forλ=5.3×10-7 m and using the value of n2 derived earlier, we obtain an estimate of Pcr≈3×10-3 W. Knowing the explicit form of n2, it can be argued that the power Pcr does not depend on the width of the light beam.


## 4. Discussion

The approach that we have used lets us find an explicit solution of the spectral problem in (2) and (3), which describes the spatial-time dynamics of nanoparticles in a liquid medium and confirms that there is a response. Some previous theoretical descriptions of the self-action of light beams have relied on equations that assume a slowly changing amplitude and a nonstationary heat equation (a diffusion equation is used here), but obtaining an analytical solution for such systems is extremely difficult; problems involving heat nonlinearity are usually resolved by numerical methods. The explicit expression for the nonlinear part of the index of refraction that we have calculated here is caused by the action of the gradient force of the electric field of the light wave, and we have shown that it depends nonlinearly on the intensity.With the exception of [11, 34, 35], we could not find any publications that have experimentally studied the contribution to the concentration nonlinearity directly associated with electrostriction. In these published studies, the simultaneous contributions of both thermal diffusion and electrostriction were investigated. In two experiments [34, 35], the radiation source had a power of less than 500 mW. We believe that with a higher radiation intensity for particles and a liquid with low absorption, it will be possible to distinguish the electrostriction contribution of the flow because the amount of light energy that is converted to heat is minimized. In our opinion, it is interesting that the intensity limitation appears naturally both in the thermodynamics and hydrodynamics. Note that this kind of limitation is absent in the self-action due to the heat nonlinearity of the medium, although both optical nonlinearity (heat and concentration) mechanisms are inertial.In a previous paper [27], the dynamics of the nanoparticle concentration and the effects of self-action due to both thermal diffusion and electrostriction were described but within a small perturbation approximation. Therefore, for a more detailed study of these processes, it is necessary to solve (2) using (3) with the term that describes the thermal diffusion flux and to assess the contributions of the respective flows. This will be the topic of our future research.


## 5. Conclusion

We have described the spatial-temporal dynamics of nanoparticles in liquid-phase media and the electrostriction response using an approach that provides an exact solution to the spectral problem. Based on these results, we calculated an explicit expression for the nonlinear part of the refractive index due to the action of the gradient force exerted by the electric field of the light wave and an expression for the concentration lens focal length. We showed that this part of the index of refraction depends nonlinearly on the radiation intensity. In our view, the constraints on the value of the radiation intensity are inherent to the framework of thermodynamics and hydrodynamics.



---

*Source: 591087-2013-09-29.xml*