Prediction of radiation pressure force exerted on moving particles by the two-level skeletonization

A fast full-wave method for computing radiation pressure force (RPF) exerted by shaped light beams on moving particles is presented. The problem of evaluating RPF exerted on a moving particle by a single excitation beam is converted into that of computing RPF's exerted on a static particle by multiple beams. The discretization of different beams leads to distinct right hand sides (RHS's) for the matrix system. To avoid solving each RHS by the brute-force manner, the algorithm conducts low-rank decomposition on the excitation matrix consisting of all RHS's to figure out the so-called skeleton light beams by interpolative decomposition (ID). The peak memory requirement of the skeletonization is a bottle-neck if the particle is large. A two-level skeletonization scheme is proposed to solve this problem. Some numerical experiments on arbitrarily shaped homogeneous particles are performed to illustrate the performance and capability of the developed method.