Radiation Force Function Research Articles

Optical radiation force and torque of light-sheets on a cylindrical particle near an infinite boundary

This work aims to extend the previous studies on the optical radiation force and torque for a cylindrical dielectric particle to the case near an infinite boundary. Without loss of generality, the two-dimensional cylindrical object is assumed to have an arbitrary shape and be immersed in any light-sheet beam of arbitrary wave front. The partial-wave series expansion method in cylindrical coordinates and the method of images are utilized to derive the exact expressions for the axial and transverse radiation force and torque functions, which are dependent on the beam shape coefficients of the incident light-sheet and the scattering coefficients characterized by the boundary conditions at the surface. Numerical computations are performed for a circular cylindrical particle illuminated by a two-dimensional plane wave and a two-dimensional Gauss beam, respectively, with particular emphases on the size parameter, the refractive index, the particle-boundary distance, the beam waist and the offset shifts. The radiation force will increase in magnitude and reverse its direction under selected conditions. When the absorptive cylindrical particle is shifted off-axially with respect to the beam axis, it will be rotated counterclockwise or clockwise by the radiation torque. The results obtained in this work have potential applications in non-contact particle manipulation and transportation using optical tweezers.

Induced effects of a rigid boundary on a pulsating spherical source

This work develops a comprehensive analytical formalism to study the boundary-induced effects on a pulsating spherical radiating source near a rigid planar boundary in an inviscid fluid. Using the modal multipole expansion method, the image method and the translational addition theorem for spherical wave functions, exact closed-form expressions are derived for the dimensionless axial acoustic radiation force function as well as the radiation, amplification and extinction energy efficiency factors. The source-boundary distance and the dimensionless size parameter are especially addressed in the numerical simulations for spherical sources with different vibration modes. It is shown that the axial radiation force is positive or negative, depending on the choice of parameters. Moreover, cases when the extinction acoustic power completely vanishes are predicted theoretically for dipole and quadrupole vibrations, which can be used to achieve acoustic invisibility in acoustic cloaking devices and acoustic stealth technologies. Combination of more than one vibration modes is also investigated to better simulate the acoustic manipulation in practical applications. This work can help us understand the underlying mechanism in manipulating the active spherical radiating sources in a bounded field.

Sustainable pulling motion of an active scatterer

In this paper, we concern the generation of attractive net motion with respect to the location of external wave source for active spherical carriers. Here we recall that the exerted acoustic radiation force in an acoustic field resulting from an incident wave, a scattered wave, and the radiated wave from the active carrier, can be positive, negative, or zero based on the location of the object in the field. Thus, in a general case, sustainable pulling motion is not guaranteed. In this work, by considering the point that the average of the radiation force for an object active in a single mode, over a complete wavelength, is equal to the radiation force applied to a passive object (which is always positive), we put forward a technique to generate acoustic radiation force with a negative average over a complete wavelength to ensure that the net motion is attractive. The idea here is to simultaneously excite of the object in two modes of breathing and first, with a difference in the phase of excitations. We will also show that by controlling the phase difference, the average of the force exerted on the object can be positive, negative or zero. Moreover, we show that for specific phase differences not only the average of the force is negative but also the force itself never experiences a positive value in a whole wavelength at those phase differences, which this can be a desired state to achieve a perfect negative net motion. The formulation of the paper is developed based on the velocity distribution over the surface of the sphere, but an implementation method with piezoelectric actuators is suggested. Here, by setting the zero-radiation force state (i.e., radiation force function cancellation) as a reasonable criterion for determination of the required amplitude for surface velocity or the operation voltage, we calculate the averaged force vs the frequency and the phase difference. The current work, with the hope of achieving complete maneuverability extends the concept of manipulation of smart carriers and reinforces the literature on single beam acoustic handling techniques.

Axial acoustic radiation force on an elastic spherical shell near an impedance boundary for zero-order quasi-Bessel–Gauss beam* *Project supported by the National Natural Science Foundation of China (Grant Nos. 81527901, 11604361, and 91630309).

Shell structures have increasingly widespread applications in biomedical ultrasound fields such as contrast agents and drug delivery, which requires the precise prediction of the acoustic radiation force under various circumstances to improve the system efficiency. The acoustic radiation force exerted by a zero-order quasi-Bessel–Gauss beam on an elastic spherical shell near an impedance boundary is theoretically and numerically studied in this study. By means of the finite series method and the image theory, a zero-order quasi-Bessel–Gauss beam is expanded in terms of spherical harmonic functions, and the exact solution of the acoustic radiation force is derived based on the acoustic scattering theory. The acoustic radiation force function, which represents the radiation force per unit energy density and per unit cross-sectional surface, is especially investigated. Some simulated results for a polymethyl methacrylate shell and an aluminum shell are provided to illustrate the behavior of acoustic radiation force in this case. The simulated results show the oscillatory property and the negative radiation force caused by the impedance boundary. An appropriate relative thickness of the shell can generate sharp peaks for a polymethyl methacrylate shell. Strong radiation force can be obtained at small half-cone angles and the beam waist only affects the results at high frequencies. Considering that the quasi-Bessel–Gauss beam possesses both the energy focusing property and the non-diffracting advantage, this study is expected to be useful in the development of acoustic tweezers, contrast agent micro-shells, and drug delivery applications.

Optical TM ⇄ TE mode conversion contribution to the radiation force on a cylinder exhibiting rotary polarization in circularly polarized light

Circularly polarized light can induce additional mode conversion when it interacts with a material that rotates the orientation of the plane of polarization, which is known as “optical rotation” or “rotary polarization”. The effect of this additional mode conversion phenomenon is studied in the framework of the optical radiation force theory. The time-averaged force (per-length) acting on an infinitely long perfect electromagnetic conductor (PEMC) circular cylinder is considered. The multipole modal expansion method in cylindrical coordinates is utilized to derive exact series expansions for the components of the longitudinal radiation force per-length (i.e. acting along the direction of wave propagation) applicable to any range of frequencies. Numerical predictions for the radiation force function (which is the radiation force per unit energy density and cross-sectional surface) and its components clearly demonstrate the contribution of the co-polarized and cross-polarized waves. It is shown that the total radiation force is not the mere sum of the TM and TE mode components in a circularly polarized wave-field. There is an additional cross-interference term describing the contributions of TM ⇄ TE mode conversion to the total force that cannot be neglected. The computations show that the interference term can be positive, null or negative as the dimensionless size parameter ka and admittance M of the cylinder vary (where k is the wavenumber in the medium of wave propagation and a is the radius of the cylinder). The results show that the total force is always positive and acts in the direction of wave propagation (i.e., force of repulsion). The analysis is further extended to compute the energy efficiencies, demonstrating the validity of the results from the standpoint of energy conservation. This work generalizes the radiation force formalism for linearly polarized fields by introducing additional terms corresponding to TM ⇄ TE mode conversions in the series expansion for the longitudinal radiation force on cylindrical materials exhibiting rotary polarization in circularly polarized plane progressive waves.

Dynamic oscillatory radiation force in optical heterodyning

In contrast with the optical radiation force induced by conventional lasers of continuous mono-frequency wave-fields, a dynamic or oscillatory mode arises when a beating effect caused by intensity changes over time (i.e., amplitude modulation) occurs. The purpose of this analysis is directed toward examining theoretically the dynamic (oscillatory) radiation force in optical heterodyning, caused by interfering/mixing two electromagnetic/optical plane waves driven at slightly different frequencies. The example of a dielectric cylinder material having a circular geometric cross-section is considered. Based on the integration of Maxwell's stress tensor over the surface of the object in the near-field, the modal expansion method in cylindrical coordinates is used in conjunction with the short-term time average to obtain mathematical series for the longitudinal dynamic optical radiation force per-length (i.e. acting along the direction of wave propagation). The cases of incident TM and TE polarized plane progressive waves are considered. Numerical illustrative results for the dimensionless dynamic radiation force function are performed with particular emphasis on changing the size parameter of the cylinder as well as the normalized difference-frequency and time parameters. The results show that the total force is not the mere sum of the components originated by the primary waves. There is an additional cross-term factor related to the dynamic component which cannot be neglected, suggesting that the radiation force phenomenon using beating optical waves is slowly time-varying (oscillatory). Moreover, a resonance splitting effect in the plots of the dynamic radiation force is observed and discussed. The present analysis has the potential to open a novel method in the development of dynamic/oscillatory optical heterodyne tweezers and tractor beams for particle manipulation and characterization.

HYDRODYNAMIC ASSESMENT ON MOTION PERFORMANCE CHARACTERISTICS OF GENERIC WAVE ENERGY POINT ABSORBERS

This work condenses various modeling techniques for different Point Absorber configurations. An alternating frequency - time domain model is implemented in MatlabFORTRAN in order to compute the displacement, velocities and the power absorbed in the heave mode for both single and multiple body configurations. Coupling of different degrees of freedom are merely contemplated with regard to a single buoy motion. NEMOH and BEMIO solvers are applied in the solution of Newtons second law according to the Boundary Element Methodology. Initially, this Wave-to-Wire model is validated through comparison with previous experimental results for a floating cone cylinder shape (Buldra-FO3). A single, generic, vertical floating cylinder is contemplated then, that responds to the action of the passing regular waves excitation. Later, two equally sized vertical floating cylinders aligned with the incident wave direction are modeled for a variable distance between the bodies. In deep water, we approximate the convolutive radiation force function term through the Prony method. Using for instance triangular or diamond shaped arrays of three and four bodies respectively, the study delves into the interaction effects for regular waves. The results highlight the most efficient layout for maximizing the energy production whilst providing important insights into their performance, revealing displacement amplification-, capture width-ratios and the commonly known park effect.

Optical radiation force circular dichroism spectroscopy

This work introduces the method of circular dichroism spectroscopy in the framework of the electromagnetic/optical radiation force theory. This technique is defined here as the difference in radiation force of left-handed and right-handed circularly polarized electromagnetic waves illuminating an object exhibiting rotary polarization. The example of a lossless material, such as the perfect electromagnetic conductor (PEMC) cylinder having a circular geometric cross-section, is considered. The modal expansion method in cylindrical coordinates is used to obtain exact mathematical series expansions for the longitudinal radiation force per-length (i.e. acting along the direction of wave propagation) considering left-handed and right-handed circularly polarized cylindrically diverging waves emanating from a line source. The case of plane progressive waves is recovered when the source is located far from the cylinder. Numerical illustrative results for the dimensionless radiation force functions as well as the scattering, extinction and absorption energy efficiencies and their co-polarized and cross-polarized components are performed with particular emphasis on the size parameter of the cylinder, the dimensionless distance parameter from the line source, and the admittance of the cylinder. The results reveal that the individual radiation force functions for left-handed and right-handed circularly polarized waves can be negative, zero, or positive depending on the cylinder distance from the source. Moreover, the optical radiation force circular dichroism (ORFCD) and the extinction energy efficiency circular dichroism (EEECD) are positive for a negative admittance of the cylinder, while they reverse sign for a positive admittance. While the EEECD shows some form of symmetry versus admittance sign change, the ORFCD does not. The possibility of achieving invisibility cloaking for a small PEMC cylinder is also predicted. The present ORFCD spectroscopy method is applicable to any cylinder material exhibiting rotary polarization such as chiral, topological insulator, plasma, liquid crystal etc.

Acoustic radiation force in standing and quasi-standing high-order Bessel beams on a multilayered sphere

The acoustic radiation force exerted by a standing and quasi-standing high-order Bessel beam on a multilayered sphere immersed in an ideal fluid is theoretically and numerically studied in this work. By means of the finite series method, the incident beam is expanded in terms of spherical harmonic functions. An analytical expression of the acoustic radiation force is derived based on the sound scattering theory. The dimensionless radiation force function versus ka is also simulated, with a particular emphasis on the relative thickness of each layer as well as the half-cone angle. To better simulate the practical manipulation, the effect of sound absorption is also analyzed. The simulated results of the first three orders reveal that the radius of each layer affects the position of peaks and dips of the curves, and acoustic trapping can be achieved at selected ka. It is also found that the curves exhibit the same general trend but the opposite directions when the order is 1 and 2, respectively. This study is expected to be useful in acoustic manipulation, drug delivery in the field of biomedical ultrasound and material sciences.

Induced radiation force of an optical line source on a cylinder material exhibiting circular dichroism.

The optical radiation force experienced by a cylinder material of circular cross section exhibiting circular dichroism (known also as rotary polarization) in an electric line source illumination is considered. An exact analytical expression for the radiation force (per length) valid for any frequency range is derived assuming an electric line source radiating cylindrically diverging TM-polarized waves without any approximations. The partial-wave series expansion method in cylindrical coordinates utilizing standard Bessel and Hankel functions is used to derive the electric and magnetic field expressions and a dimensionless radiation force function (or efficiency), which depends on the scattering coefficient of the cylinder as well as the distance from the radiating source. To illustrate the analysis, numerical computations for the dimensionless radiation force function for a perfect electromagnetic conductor (PEMC) cylinder are performed with emphasis on its dimensionless size parameter and source distance, which clearly draw attention to the contribution of the cross-polarized scattered waves (resulting from the rotary polarization effect) to the total force. The numerical predictions demonstrate the possibility to pull a circular-shaped cylinder material with rotary polarization toward the illuminating electric line source with TM-polarized waves using a curved wavefront depending on the PEMC material admittance, distance to the source, and size of the cylinder.

Axial acoustic radiation force on a spherical particle in a zero-order Mathieu beam.

In this paper, an expression for the acoustic radiation force on a sphere located on the axis of a zero-order Mathieu beam propagating in an ideal fluid is obtained. The simulation results of different kinds of spherical particles, including rigid spheres, fluid spheres, and elastic spheres, are presented to illustrate the theory. Curves of the acoustic radiation force function versus ka are obtained for different half-cone angles β and ellipticity parameters q of the Mathieu beam. The simulation results reveal that β and q have a direct effect on the acoustic radiation force exerted on a spherical object in an ideal fluid. The acoustic radiation force has negative values for a fluid sphere in a zero-order Mathieu beam when the half-cone angle is β≥60°. The value of the half-cone angle β also affects the direction of the acoustic radiation force for an elastic sphere. When the half-cone angle is set as a constant, the greater the ellipticity parameter is, the slower the decrease in the value of the acoustic radiation force function for a spherical particle. The conclusions of this work may be helpful for the design of acoustic tweezers and may also extend the potential applications of acoustic manipulation technology.

Electromagnetic radiation force on a perfect electromagnetic conductor (PEMC) circular cylinder

Unlike isotropic dielectric objects, the interaction of incident electromagnetic (EM) or optical waves with a material allowing rotary polarization produces coupled polarized internal and scattered fields. The aim of this investigation is to examine theoretically a novel physical effect, which concerns the contributions of the co-polarized and cross-polarized fields to the radiation force (per-length) experienced by an infinitely long perfect electromagnetic conductor (PEMC) cylinder having a circular cross-section and illuminated by TM-polarized plane progressive waves propagating perpendicularly to its axis. The multipole partial-wave series expansion method in cylindrical coordinates is used to derive exact series expansions for the co-polarized and cross-polarized components of the longitudinal radiation force per-length (i.e. acting along the direction of wave propagation). In contrast with the perfect electric or magnetic conductors (PECs or PMCs), or the dielectric cylinder case, numerical illustrative results for the radiation force function (which is the radiation force per unit energy density and cross-sectional surface) clearly demonstrate the contribution of the cross-polarized component of the radiation force for a PEMC cylinder allowing rotary polarization. The results show that the cross-polarized component of the radiation force function can be positive or negative as the dimensionless frequency parameter ka varies (where k is the wavenumber in the medium of wave propagation and a is the radius of the cylinder). Moreover, it vanishes for ka = 0.67946 regardless of the admittance of the PEMC cylinder. Notice that the total force (i.e. the sum of the co-polarized and cross-polarized components) is always repulsive (i.e., positive). It is also verified that the results are in complete agreement with the law of energy conservation applied to scattering. The present analysis generalizes the classical radiation force investigations by introducing extra new terms in the series expansion for the longitudinal radiation force function for cylindrical materials allowing rotary polarization.

Calculation and optimization of acoustic radiation force produced by a two-dimensional transducer array

The beam-shape coefficients of a focused sound field from a two-dimensional transducer array are obtained by using the translational addition theorem for spherical functions, and the acoustic radiation force function of a sphere is given in terms of the beam-shape coefficients and the scattering coefficients. The numerical results of acoustic radiation force function show that when the radius of the sphere is much smaller than the acoustic wavelength, the acoustic radiation force function increases sharply with the increase in the radius, while it begins to decrease as the radius continues to increase. Increasing the vibration frequency of the transducer, the array element spacing and the number of elements will enhance the directivity of the acoustic beam and make the acoustic radiation force distribution much more concentrated at the bottom of the sphere. Through optimization of acoustic radiation force distribution by adjusting the phase distribution of the transducer array, the acoustic radiation force distribution range at the bottom of the sphere becomes wider.

Optical radiation force (per–length) on an electrically conducting elliptical cylinder having a smooth or ribbed surface

The aim of this work is to develop a formal semi-analytical model using the modal expansion method in cylindrical coordinates to calculate the optical/electromagnetic (EM) radiation force-per-length experienced by an infinitely long electrically-conducting elliptical cylinder having a smooth or wavy/corrugated surface in EM plane progressive waves with different polarizations. One of the semi-axes of the elliptical cylinder coincides with the direction of the incident field. The modal matching method is used to determine the scattering coefficients by matrix inversion. Standard cylindrical (Bessel and Hankel) wave functions are used. Simplified expressions leading to exact series expansions for the optical/EM radiation forces assuming either electric (TM) of magnetic (TE) plane wave incidences are provided without any approximations, in addition to integral equations demonstrating the direct relationship of the radiation force with the square of the scattered field magnitude. Numerical computations for the non-dimensional radiation force function are performed for electrically conducting elliptic and circular cylinders having a smooth or ribbed/corrugated surface. Adequate convergence plots confirm the validity and correctness of the method to evaluate the radiation force with no limitation to a particular frequency range (i.e. Rayleigh, Mie or geometrical optics regimes). Emphases are given on the aspect ratio, the non-dimensional size of the cylinder, the corrugation characteristic of its surface, and the polarization of the incident field. The results are particularly relevant in optical tweezers and other related applications in fluid dynamics, where the shape and stability of a cylindrical drop stressed by a uniform external electric/magnetic field are altered. Furthermore, a direct analogy with the acoustical counterpart is discussed.

Electromagnetic binding and radiation force reversal on a pair of electrically conducting cylinders of arbitrary geometrical cross-section with smooth and corrugated surfaces

The electromagnetic (EM) radiation force-per-length exerted on a pair of electrically-conducting cylindrical particles of circular and non-circular cross-sections is examined using a formal semi-analytical method based on boundary matching in cylindrical coordinates. Initially, the scattering coefficients of the particle pair are determined by imposing suitable boundary conditions leading linear systems of equations computed via matrix inversion and numerical integration procedures. Standard cylindrical (Bessel and Hankel) wave functions are used and closed-form expressions for the dimensionless longitudinal and transverse radiation force functions are evaluated assuming either magnetic (TE) or electric (TM) plane wave incidences. Particle pairs with smooth and corrugated surfaces are considered and numerical computations are performed with emphasis on the distance separating their centers of mass, the angle of incidence of the incident illuminating field and the surface roughness. Adequate convergence plots confirm the validity of the method to evaluate the radiation force functions, and the model is adaptable to any frequency range (i.e. Rayleigh, Mie or geometrical optics regimes). The results can find potential applications in optical tweezers and other related applications in fluid dynamics. In addition, the acoustical analogue is discussed.

Acoustic radiation force on a rigid cylinder in an off-axis Gaussian beam near an impedance boundary

The exact equations of the axial and transverse acoustic radiation force functions of a Gaussian beam arbitrarily incident on an infinite rigid cylinder close to an impedance boundary and immersed in an ideal fluid are deduced by expressing the incident wave, the scattering wave and the boundary reflected wave in terms of the cylindrical wave function. The effects of the beam waist, the sound reflection coefficient, the cylinder position and the distance from the impedance boundary on the acoustic radiation force are studied using numerical simulations. The simulation results show that the amplitude of the acoustic radiation force function increases with beam width. Moreover, the values of the acoustic radiation force in both the axial and transverse directions reach those of a plane wave when the beam width is considerably larger than the wavelength of the Gaussian beam. The properties of the impedance boundary and the position of the cylinder in the Gaussian beam have a considerable effect on the magnitude and direction of the force. The simulation results, particularly in the case of a transverse force, indicate the presence of a negative acoustic radiation force that is related to the nondimensional frequency and position of the cylinder in the Gaussian beam.

Acoustic radiation force of a solid elastic sphere immersed in a cylindrical cavity filled with ideal fluid

An expression for the acoustic radiation force function on a solid elastic spherical particle placed in an infinite rigid cylindrical cavity filled with an ideal fluid is deduced when the incident wave is a plane progressive wave propagated along the cylindrical axis. The acoustic radiation force of the spherical particle with different materials was computed to validate the theory. The simulation results demonstrate that the acoustic radiation force changes demonstrably because of the influence of the reflective acoustic wave from the cylindrical cavity. The sharp resonance peaks, which result from the resonance of the fluid-filled cylindrical cavity, appear at the same positions in the acoustic radiation force curve for the spherical particle with different radii and materials. Relative radius, which is the ratio of the sphere radius and the cylindrical cavity radius, has more influence on acoustic radiation force. Moreover, the negative radiation forces, which are opposite to the progressive directions of the plane wave, are observed at certain frequencies.

Acoustic radiation force on a cylindrical particle near a planar rigid boundary

The aim of this investigation is to provide improved mathematical series expansions of the longitudinal and transverse acoustic radiation forces for a rigid cylindrical particle in 2D of arbitrary cross-section located near a planar rigid wall. Incident plane progressive waves with variable angle of incidence are considered in a non-viscous fluid. The multiple scattering effects occurring between the particle and the rigid boundary are described using the partial-wave decomposition in cylindrical coordinates, the method of images and the translational addition theorem. Initially, an effective acoustic field incident on the particle is defined, which includes the primary incident field, the reflected waves from the flat wall and the scattered field from the image object. Subsequently, the incident effective field along with the scattered field from the object are utilized to obtain closed-form mathematical expressions for the longitudinal and transverse radiation force functions, based on a scattering approach in the far-field. The radiation force vector components are formulated in partial-wave series in cylindrical coordinates, which involve the incidence angle, the expansion coefficients of the scatterer and its image, and the distance from the center of mass of the particle to the boundary. Numerical examples for a rigid circular cylinder are considered. Computations for the longitudinal and transverse non-dimensional radiation force functions are performed. Emphasis is given on varying the size of the particle, the incidence angle of the source field and the particle-wall distance. Depending on the particle-wall distance and incidence angle, zero-longitudinal and transverse force components arise, thus, the particle becomes unaffected by the linear momentum transfer. Moreover, pushing or pulling forces between the particle and wall are predicted depending on the particle-wall distance, the incidence angle and size parameter. The results may find possible applications in the development of acousto-fluidic devices, acoustic levitation of particles nearby a boundary, cloaking/invisibility, and underwater acoustics to name a few areas, where most investigations resort initially to numerical simulations to guide the experimental design processes.

Acoustic radiation force induced by two Airy–Gaussian beams on a cylindrical particle**Project supported by the National Key R & D Program, China (Grant No. 2016YFF0203000), the National Natural Science Foundation of China (Grant Nos. 11774167 and 61571222), Fundamental Research Funds for the Central Universities, China (Grant No. 020414380001), State Key Laboratory of Acoustics, Chinese Academy of Sciences (Grant No. SKLA201609),

Based on the angular spectrum decomposition and partial-wave series expansion methods, we investigate the radiation force functions of two Airy–Gaussian (AiG) beams on a cylindrical particle and the motion trajectory of the particle. The simulations show that the particle can be pulled or propelled into either the positive or negative transverse direction by turning the phase difference between the two AiG beams appropriately; and the larger the beam widths of the two AiG beams are, the bigger the radiation force can be obtained to control the particle. In addition, the direction of the accelerated particle can be controlled while the dimensionless frequency bandwidth changes. The results indicate that the phase plays an important role in controlling the direction of the particle, which may provide a theoretical basis for the design of acoustical tweezers and the development of drug delivery.

Pushing, pulling and electromagnetic radiation force cloaking by a pair of conducting cylindrical particles

The present analysis shows that two conducting cylindrical particles illuminated by an axially-polarized electric field of plane progressive waves at arbitrary incidence will attract, repel or become totally cloaked (i.e., invisible to the transfer of linear momentum carried by the incident waves), depending on their sizes, the interparticle distance as well as the angle of incidence of the incident field. Based on the rigorous multipole expansion method and the translational addition theorem of cylindrical wave functions, the electromagnetic (EM) radiation forces arising from multiple scattering effects between a pair of perfectly conducting cylindrical particles of circular cross-sections are derived and computed. An effective incident field on a particular particle is determined first, and used subsequently with its corresponding scattered field to derive the closed-form analytical expressions for the radiation force vector components. The mathematical expressions for the EM radiation force components (i.e. longitudinal and transverse) are exact, and have been formulated in partial-wave series expansions in cylindrical coordinates involving the angle of incidence, the interparticle distance and the expansion coefficients. Numerical examples illustrate the analysis for two perfectly conducting circular cylinders in a homogeneous nonmagnetic medium of wave propagation. The computations for the dimensionless radiation force functions are performed with particular emphasis on varying the angle of incidence, the interparticle distance, and the sizes of the particles. Depending on the interparticle distance and angle of incidence, the cylinders yield total neutrality (or invisibility); they experience no force and become unresponsive to the transfer of the EM linear momentum due to multiple scattering cancellation effects. Moreover, pushing or pulling EM forces between the two cylinders arise depending on the interparticle distance, the angle of incidence and their size parameters. This study provides a complete analytical method and computations for the longitudinal and transverse radiation force components in the multiple scattering of EM plane progressive waves with potential applications in particle manipulation, optically-engineered metamaterials with reconfigurable periodicities and cloaking devices to name a few examples.