Matlab code for a multipolar decomposition of optical forces

This dataset supplements Figure 5 from the publication "Multipolar Origin of the Unexpected Transverse
Force Resulting from Two-Wave Interference*" by Karim Achouri, Andrei Kiselev, and Olivier J. F. Martin. Here, we provide the code for the multipolar analysis of the optical force based on the vector spherical harmonic decomposition. Based on the Mie solution, electric fields can be obtained numerically in the far-field by using the software developed by Dr. Karim Achouri https://github.com/kachourim/MieScatteringPEC. This code analyzes the far-field scattered by a perfect electric conductor sphere placed in vacuum for different sphere radii. In the framework of the Maxwell stress tensor, we find the optical force acting on a sphere along the z-direction in the illumination configuration presented in Figure 1*. The use of the vector spherical decomposition allows to observe contributions from different multipolar interactions. In this code,we analyze the force appearing as a result of the interaction between electric dipole aligned along x, px, and electric quadrupole with components xz, Qexz and compare it with the total force along the z-axis appearing as a result of interaction between all multipoles supported by the sphere with given radius.