adda
adda copied to clipboard
Optimize calculation of Csca and g
Use weighting function for integration of scattered fields to obtain Csca
and g. Take into account at least the diffraction peak. This is similar to
what is currently done in DDSCAT.
Original issue reported on code.google.com by yurkin
on 28 Nov 2008 at 6:53
- Blocked on: #138
It should be noted that this issue is significant only for large particles. And
these particles usually lead to large number of iterations. So even a very
large number of angles used for integration (e.g. Jmax=10) is still relatively
fast, i.e. its computation time is comparable or smaller than the time to get
internal fields.
Original comment by yurkin
on 13 Jun 2011 at 5:00
[deleted comment]
This can be quickly solved through spherical-harmonics expansion of the
scattered fields.
Original comment by yurkin
on 4 Jul 2013 at 10:00
- Now blocked on: #138
Also, certain optimizations are pending for the surface mode. There, the upper
and lower hemispheres should be treated separately. Just separating the range
of theta into two should have some effect (due to bad behavior of Romberg
integration for step functions).
Original comment by yurkin
on 4 Feb 2014 at 5:04
Concerning the spherical-harmonics expansion of the internal field. Csca is then a sum of squares of coefficients of multipoles. The latter is well-known, but interestingly it also follows from Csca ~ P*.G(I).P [1], since the G(I) is diagonal in the spherical-harmonics basis [2].
- Moskalensky A.E. and Yurkin M.A. Energy budget and optical theorem for scattering of source-induced fields, Phys. Rev. A 99, 053824 (2019).
- Moskalensky A.E. and Yurkin M.A. A point electric dipole: from basic optical properties to the fluctuation-dissipation theorem, Rev. Phys. 6, 100047 (2021).