Instantaneous axial force of a high-order Bessel vortex beam of acoustic waves incident upon a rigid movable sphere

The present investigation examines the instantaneous force resulting from the interaction of an acoustical high-order Bessel vortex beam (HOBVB) with a rigid sphere. The rigid sphere case is important in fluid dynamics applications because it perfectly simulates the interaction of instantaneous sound waves in a reduced gravity environment with a levitated spherical liquid soft drop in air. Here, a closed-form solution for the instantaneous force involving the total pressure field as well as the Bessel beam parameters is obtained for the case of progressive, stationary and quasi-stationary waves. Instantaneous force examples for progressive waves are computed for both a fixed and a movable rigid sphere. The results show how the instantaneous force per unit cross-sectional surface and unit pressure varies versus the dimensionless frequency ka (k is the wave number in the fluid medium and a is the sphere’s radius), the half-cone angle β and the order m of the HOBVB. It is demonstrated here that the instantaneous force is determined only for (m, n) = (0, 1) (where n is the partial-wave number), and vanishes for m > 0 because of symmetry. In addition, the instantaneous force and normalized amplitude velocity results are computed and compared with those of a rigid immovable (fixed) sphere. It is shown that they differ significantly for ka values below 5. The proposed analysis may be of interest in the analysis of instantaneous forces on spherical particles for particle manipulation, filtering, trapping and drug delivery. The presented solutions may also serve as a method for comparison to other solutions obtained by strictly numerical or asymptotic approaches.     

Axial and transverse acoustic radiation forces on a fluid sphere placed arbitrarily in Bessel beam standing wave tweezers

The axial and transverse radiation forces on a fluid sphere placed arbitrarily in the acoustical field of Bessel beams of standing waves are evaluated. The three-dimensional components of the time-averaged force are expressed in terms of the beam-shape coefficients of the incident field and the scattering coefficients of the fluid sphere using a partial-wave expansion (PWE) method. Examples are chosen for which the standing wave field is composed of either a zero-order (non-vortex) Bessel beam, or a first-order Bessel vortex beam. It is shown here, that both transverse and axial forces can push or pull the fluid sphere to an equilibrium position depending on the chosen size parameter ka   (where k$k$ is the wave-number and a$a$ the sphere’s radius). The corresponding results are of particular importance in biophysical applications for the design of lab-on-chip devices operating with Bessel beams standing wave tweezers. Moreover, potential investigations in acoustic levitation and related applications in particle rotation in a vortex beam may benefit from the results of this study.     

Acoustic scattering of a high-order Bessel beam by an elastic sphere

The exact analytical solution for the scattering of a generalized (or “hollow”) acoustic Bessel beam in water by an elastic sphere centered on the beam is presented. The far-field acoustic scattering field is expressed as a partial wave series involving the scattering angle relative to the beam axis and the half-conical angle of the wave vector components of the generalized Bessel beam. The sphere is assumed to have isotropic elastic material properties so that the nth partial wave amplitude for plane wave scattering is proportional to a known partial-wave coefficient. The transverse acoustic scattering field is investigated versus the dimensionless parameter ka(k is the wave vector, a radius of the sphere) as well as the polar angle θ for a specific dimensionless frequency and half-cone angle β. For higher-order generalized beams, the acoustic scattering vanishes in the backward (θ = π) and forward (θ = 0) directions along the beam axis. Moreover it is possible to suppress the excitation of certain resonances of an elastic sphere by appropriate selection of the generalized Bessel beam parameters.     

Acoustic radiation force of high-order Bessel beam standing wave tweezers on a rigid sphere

Background and objective

Particle manipulation using the acoustic radiation force of Bessel beams is an active field of research. In a previous investigation, [F.G. Mitri, Acoustic radiation force on a sphere in standing and quasi-standing zero-order Bessel beam tweezers, Annals of Physics 323 (2008) 1604–1620] an expression for the radiation force of a zero-order Bessel beam standing wave experienced by a sphere was derived. The present work extends the analysis of the radiation force to the case of a high-order Bessel beam (HOBB) of positive order m having an angular dependence on the phase ?.

Method

The derivation for the general expression of the force is based on the formulation for the total acoustic scattering field of a HOBB by a sphere [F.G. Mitri, Acoustic scattering of a high-order Bessel beam by an elastic sphere, Annals of Physics 323 (2008) 2840–2850; F.G. Mitri, Equivalence of expressions for the acoustic scattering of a progressive high order Bessel beam by an elastic sphere, IEEE Transactions on Ultrasonics Ferroelectrics and Frequency Control 56 (2009) 1100–1103] to derive the general expression for the radiation force function Y_Jm,st(ka,β,m), which is the radiation force per unit characteristic energy density and unit cross-sectional surface. The radiation force function is expressed as a generalized partial wave series involving the half-cone angle β of the wave-number components and the order m of the HOBB.

Results

Numerical results for the radiation force function of a first and a second-order Bessel beam standing wave incident upon a rigid sphere immersed in non-viscous water are computed. The rigid sphere calculations for Y_Jm,st(ka,β,m) show that the force is generally directed to a pressure node when m is a positive even integer number (i.e. Y_Jm,st(ka,β,m)>0), whereas the force is generally directed toward a pressure antinode when m is a positive odd integer number (i.e. Y_Jm,st(ka,β,m)<0).

Conclusion

An expression is derived for the radiation force on a rigid sphere placed along the axis of an ideal non-diffracting HOBB of acoustic standing (or stationary) waves propagating in an ideal fluid. The formulation includes results of a previous work done for a zero-order Bessel beam standing wave (m = 0). The proposed theory is of particular interest essentially due to its inherent value as a canonical problem in particle manipulation using the acoustic radiation force of a HOBB standing wave on a sphere. It may also serve as the benchmark for comparison to other solutions obtained by strictly numerical or asymptotic approaches.     

Conversion of high-order Laguerre-Gaussian beams into Bessel beams of increased, reduced or zeroth order by use of a helical axicon

We propose a new method for transformation of a Laguerre-Gaussian beam of azimuthal index l and radial index n = 0 (LG_l,0) into a vortex, diverging or nondiverging Bessel beam, which can have increased or decreased phase singularity order, or into a zeroth order Bessel beam, by means of a helical axicon. The Bessel beam divergence or nondivergence depends upon the waist position of the input Laguerre-Gaussian beam, regarding the plane where the helical axicon is situated.The expressions for the amplitude and the intensity distribution of the diffracted wave field, in the process of Fresnel diffraction, are deduced using the stationary phase method. The theoretical analysis for the vortex radius and the maximum propagation distance of the Bessel beams obtained is presented.     

Acoustic beam interaction with a rigid sphere: The case of a first-order non-diffracting Bessel trigonometric beam

Mathematical expressions for the acoustic scattering, instantaneous (linear), and time-averaged (nonlinear) forces resulting from the interaction of a new type of Bessel beam, termed here a first-order non-diffracting Bessel trigonometric beam (FOBTB) with a sphere, are derived. The beam is termed “trigonometric” because of the dependence of its phase on the cosine function. The FOBTB is regarded as a superposition of two equi-amplitude first-order Bessel vortex (helicoidal) beams having a unit positive and negative order (known also as topological charge), respectively. The FOBTB is non-diffracting, possesses an axial null, a geometric phase, and has an azimuthal phase that depends on cos(?±?₀), where ?₀ is an initial arbitrary phase angle. Beam rotation around its wave propagation axis can be achieved by varying ?₀. The 3D directivity patterns are computed, and the resulting modifications of the scattering are illustrated for a rigid sphere centered on the beam's axis and immersed in water. Moreover, the backward and forward acoustic scattering by a sphere vanish for all frequencies. The present paper will shed light on the novel scattering properties of an acoustical FOBTB by a sphere that may be useful in particle manipulation and entrapment, non-destructive/medical imaging, and may be extended to other potentially useful applications in optics and electromagnetism.     

Acoustic radiation force on a sphere in standing and quasi-standing zero-order Bessel beam tweezers

Starting from the exact acoustic scattering from a sphere immersed in an ideal fluid and centered along the propagation axis of a standing or quasi-standing zero-order Bessel beam, explicit partial-wave representations for the radiation force are derived. A standing or a quasi-standing acoustic field is the result of propagating two equal or unequal amplitude zero-order Bessel beams, respectively, along the same axis but in opposite sense. The Bessel beam is characterized by the half-cone angle β of its plane wave components, such that β = 0 represents a plane wave. It is assumed here that the half-cone angle β for each of the counter-propagating acoustic Bessel beams is equal. Fluid, elastic and viscoelastic spheres immersed in water are treated as examples. Results indicate the capability of manipulating spherical targets based on their mechanical and acoustical properties. This condition provides an impetus for further designing acoustic tweezers operating with standing or quasi-standing Bessel acoustic waves. Potential applications include particle manipulation in micro-fluidic lab-on-chips as well as in reduced gravity environments.     

Interaction of a high-order Bessel beam with a submerged spherical ultrasound contrast agent shell - Scattering theory

Background and objective

Acoustic scattering properties of ultrasound contrast agents are useful in extending existing or developing new techniques for biomedical imaging applications. A useful first step in this direction is to investigate the acoustic scattering of a new class of acoustic beams, known as helicoidal high-order Bessel beams, to improve the understanding of their scattering characteristics by an ultrasound contrast agent, which at present is very limited.

Method

The transverse acoustic scattering of a commercially available albuminoidal ultrasound contrast agent shell filled with air or a denser gas such as perfluoropropane and placed in a helicoidal Bessel beam of any order is examined numerically. The shell is assumed to possess an outer radius a = 3.5 microns and a thickness of ∼105 nm. Moduli of the total and resonance transverse acoustic scattering form functions are numerically evaluated in the bandwidth 0 < ka? 3, which corresponds to a frequency bandwidth of 0-205 MHz that covers a wide range of applications for imaging with contrast agents. Particular attention is paid to the shell’s material, the content of its interior hollow region and the fluid surrounding its exterior. The contrast agent shell is assumed to be immersed in an ideal compressible fluid so the viscous corrections are not considered. Analytical equations are derived and numerical calculations of the total and resonance form functions are performed with particular emphasis on the effect of varying the half-cone angle, the order of the helicoidal Bessel beam as well as the fluid that fills the interior hollow space.

Results and conclusion

It is shown that shell wave resonance modes can be excited on an encapsulated micro-bubble. The forward and backscattering vanish for a helicoidal high-order Bessel beam. Additionally, the fluid filling the inner core affects the shell’s response significantly. Moreover, there is no monopole contribution to the axial scattering of a helicoidal Bessel beam of order m ? 1 so that the dynamics of contrast agents would be significantly altered. The main finding of the present theory is the suppression or enhancement for a particular resonance that may be used to advantage in imaging with ultrasound contrast agents for clinical applications.     

Propagation dynamics and experimental observation of quadratic vortex beam

The concept of a quadratic vortex beam is proposed, in which phase term of the beam is given by exp(i mθ²). The phase of the quadratic vortex beam increases with azimuthal angle nonlinearly. This change in phase produces several unexpected effects. Unlike the circularly symmetric beam spot of normal vortex beams, the intensity distribution of the quadratic vortex beam is shown to be asymmetric. The phase singularities will shift in the transverse beam plane on propagation.     

Mechanism of the quasi-zero axial acoustic radiation force experienced by elastic and viscoelastic spheres in the field of a quasi-Gaussian beam and particle tweezing

The present analysis investigates the (axial) acoustic radiation force induced by a quasi-Gaussian beam centered on an elastic and a viscoelastic (polymer-type) sphere in a nonviscous fluid. The quasi-Gaussian beam is an exact solution of the source free Helmholtz wave equation and is characterized by an arbitrary waist w₀ and a diffraction convergence length known as the Rayleigh range z_R. Examples are found where the radiation force unexpectedly approaches closely to zero at some of the elastic sphere’s resonance frequencies for kw₀ ? 1 (where this range is of particular interest in describing strongly focused or divergent beams), which may produce particle immobilization along the axial direction. Moreover, the (quasi)vanishing behavior of the radiation force is found to be correlated with conditions giving extinction of the backscattering by the quasi-Gaussian beam. Furthermore, the mechanism for the quasi-zero force is studied theoretically by analyzing the contributions of the kinetic, potential and momentum flux energy densities and their density functions. It is found that all the components vanish simultaneously at the selected ka values for the nulls. However, for a viscoelastic sphere, acoustic absorption degrades the quasi-zero radiation force.     

Spatially-variable retardation plate for efficient generation of radially- and azimuthally-polarized beams

We investigate conversion of a linearly-polarized Gaussian beam to a radially- or an azimuthally-polarized doughnut (0, 1)^∗ Laguerre-Gaussian (LG) beams, performed with a spatially-variable retardation (SVR) plate. The SVR plate is composed of eight sectors of a λ/2 retardation plate, each one with different orientation of the to crystal’s slow axis. The analysis reveal that nearly-pure radially- or azimuthally-polarized LG(01)^∗ beam with M² = 2.2 can be obtained, while the transformation efficiency is 89.6%. In the experiments, performed with Nd:YAG laser, we transformed a Gaussian beam with M² = 1.3 to a radially- and azimuthally-polarized (0, 1)^∗ Laguerre-Gaussian beams with M² = 2.5. We carefully characterized the polarization state of the obtained radially- and azimuthally-polarized beams, measuring Stokes parameters. The polarization purity of the obtained beams, calculated from the measured data, was as high as 96%.     

Generation of rotating intensity blades by superposing optical vortex beams

A kind of optical beam with controllable rotating intensity blades is generated by coaxially superposing two optical vortex beams with frequency difference (Δω), different topological charges (m₁ and m₂) and equal amplitude. It is shown theoretically that the number of the blades is determined by the subtraction of topological charges (m₁ − m₂) and the angular velocity of the rotating pattern is equal to Δω / (m₁ − m₂). In our experiment, the rotating beams were generated by two optical vortex beams with opposite topological charges, where the frequency difference is acquired by a rotating plate glass. The results are quite in accordance with the simulations.     

Gegenbauer expansion to model the incident wave-field of a high-order Bessel vortex beam in spherical coordinates

The aim of this short communication is to report that Gegenbauer’s (partial-wave) expansion, that may be used (under some specific conditions) to represent the incident field of an acoustical (or optical) high-order Bessel beam (HOBB) in spherical coordinates, anticipates earlier expressions for undistorted waves. The incident wave-field is written in terms of the spherical Bessel function of the first kind, the gamma function as well as the Gegenbauer or ultraspherical functions given in terms of the associated Legendre functions when the order m of the HOBB is an integer number. Expressions for high-order and zero-order Bessel beams as well as for plane progressive waves reported in prior works can be deduced from Gegenbauer’s partial-wave expansion by appropriate choice of the beams’ parameters. Hence the value of this note becomes historical. In addition, Gegenbauer’s expansion in spherical coordinates may be used to advantage to model the wave-field of a fractional HOBB at the origin (i.e. z = 0).     

Trapping Rayleigh particles using highly focused higher-order radially polarized beams

The optical trapping characteristics of highly focused higher-order radially polarized beams (R-TEM_p1*) acting on a Rayleigh particle are studied theoretically. Numerical results show that as the order p of beam increases and the numerical aperture NA_o of the objective decreases, the axial trap distance increases but the trap depth and maximum restoring force decreases. In a limit of NA_o = 1, three higher-order R-TEM_p1* beams of p = 1, 2, 3, like the fundamental lowest-order radially polarized beam of p = 0, can three-dimensionally trap a particle to the focus but the axial trap stiffness decreases with the increase of p. When NA_o = 0.95, the focus is still a stable trap point for the two beams of p = 0 and 1 but it becomes an unstable trap point for the two beams of p = 2 and 3. The trap stability is also discussed for higher-order radially polarized beam illumination.     

Axial acoustic radiation force of progressive cylindrical diverging waves on a rigid and a soft cylinder immersed in an ideal compressible fluid

Background and motivation

Previous works investigating the radiation force of diverging spherical progressive waves incident upon spherical particles have demonstrated the direction of reversal of the force when the particle is subjected to a curved wave-front. In this communication, the analysis is extended to the case of diverging cylindrical progressive waves incident upon a rigid or a soft cylinder in a non-viscous fluid with explicit calculations for the radiation force function (which is the radiation force per unit energy density and unit cross-sectional surface) not shown in [F.G. Mitri, Ultrasonics 50 (2010) 620-627].

Method

A closed-form solution presented previously in [F.G. Mitri, Ultrasonics 50 (2010) 620-627] is used to plot the radiation force function with particular emphasis on the difference from the results of incident plane progressive waves versus the size parameter ka (k is the wave number and a is the cylinder’s radius) and the distance of the cylinder from the acoustic source r₀.

Results

Radiation force function calculations for the rigid cylinder unexpectedly reveal that under specific conditions determined by the frequency of the acoustic field, the radius of the cylinder, as well as the distance to the acoustic source, the force becomes attractive (negative force). In addition, the numerical results show that the radiation force on a rigid cylinder does not generally obey the inverse-distance law with respect to the distance from the source.

Conclusion and potential applications

These results suggest that it may be possible, under specific conditions, to pull a cylindrical structure back toward the acoustic source using progressive cylindrical diverging waves. They may also provide a means to predict the radiation force required to manipulate non-destructively a single cylindrical structure. Potential applications include the design of a new generation of acoustic tweezers operating using a single beam of progressive waves (in contrast to the traditional version of acoustical tweezers in which an acoustic standing wave field is produced using two counter-propagating acoustic fields) for investigations in the field of flow cytometry, particle manipulation and entrapment.     

Effect of the spiral phase element on the radial-polarization (0, 1) LG beam

Radially-polarized beams can be strongly amplified without significant birefringent-induced aberrations. However, radially-polarized beam is a high-order beam, and therefore has to be transformed into a fundamental Gaussian beam for reduction the beam-propagation factor M². In effort to transform the radially-polarized beam to a nearly-Gaussian beam, we consider effect of a spiral phase element (SPE) on the Laguerre-Gaussian (LG) (0, 1)^∗ beam with radial polarization, and compare this with the case when the input beam is a LG (0, 1)^∗ beam with spiral phase and uniform or random polarization. The LG (0, 1)^∗ beam with radial polarization, despite its identity in intensity profile to the beam with spiral phase, has distinctly different properties when interacting with the SPE. With the SPE and spatial filter, we transformed the radially-polarized (0, 1)^∗ mode with M² = 2.8 to a nearly-Gaussian beam with M² = 1.7. Measured transformation efficiency was 50%, and the beam brightness P/(M²)² was practically unchanged. The SPE affects polarization state of the radially-polarized beam, leading to appearance of spin angular momentum in the beam center at the far-field.     

Observation of nonstationary effects in saturation spectroscopy

We have studied the nonstationary effects in saturated absorption spectroscopy of the ⁸⁷Rb D₂ line. Varying the size of the σ⁺ polarized pump laser beam, we observed saturated absorption spectra for the σ^± polarized probe beam. For equal polarizations of the pump and probe beams, we found that the resonance signal for the F_g = 1 → F_e = 2 line, and the crossover lines between F_g = 1 → F_e = 2 and F_g = 1 → F_e = 1 (and 0) lines increased to a greater extent than the others. This observation can be understood from the calculated time evolution of the populations of the ground-state sublevels by means of a rate equation model. We also compared experimental data for other conditions with the calculated results. We found good agreement between the calculated results and the data.     

Scattering of Mathieu beams by a rigid sphere

Based on the recent results on the scattering of Bessel beams by a sphere and using the Whittaker integral, the scattering by a rigid sphere centred on a Mathieu beam is derived. The scattering field is expressed as a partial wave series involving the scattering angles relative to the beam axis and Mathieu beam parameters. Some numerical calculations are performed and it is shown that the illumination of a rigid sphere by a Mathieu beam produces asymmetrical scattering as a function of scattering angles θ and ?. The geometrical properties of the scattering Mathieu beam are noted.     

Optimal phase steps of multi-level spiral phase plates

We analyze vortex properties of the optical beams generated by a multi-level spiral phase plate (MLSPP) and find that this kind of beams consists of vortex components with topological charges of L − kM, where k is zero or any integer, M is the level number and L is the intrinsic topological charge of the MLSPP. We proved that the orbital angular momentum of the beams generated by a MLSPP reaches its maximum only if the phase steps (or the ratio of M to L) satisfies some special conditions.