next up previous contents
Next: Memory Requirements Up: User Guide for Previous: Choice of Iterative Algorithm

 

Calculation of Radiative Force and Torque

In addition to solving the scattering problem for a dipole array, DDSCAT.5a\ can compute the three-dimensional force tex2html_wrap_inline3477 and torque tex2html_wrap_inline3479 exerted on this array by the incident and scattered radiation fields. This calculation is carried out, after solving the scattering problem, provided DOTORQ has been specified in ddscat.par. For each incident polarization mode, the results are given in terms of dimensionless efficiency vectors tex2html_wrap_inline3481 and tex2html_wrap_inline3217, defined by



where tex2html_wrap_inline3477 and tex2html_wrap_inline3479 are the time-averaged force and torque on the dipole array, tex2html_wrap_inline3377 is the wavenumber in vacuo, and tex2html_wrap_inline3491 is the time-averaged energy density for an incident plane wave with amplitude tex2html_wrap_inline3493. The radiation pressure efficiency vector can be written

where tex2html_wrap_inline3495 is the direction of propagation of the incident radiation, and the vector g is the mean direction of propagation of the scattered radiation:

where tex2html_wrap_inline3497 is the element of solid angle in scattering direction tex2html_wrap_inline3499, and tex2html_wrap_inline3117 is the differential scattering cross section.

Equations for the evaluation of the radiative force and torque are derived by Draine & Weingartner (1996). It is important to note that evaluation of tex2html_wrap_inline3481 and tex2html_wrap_inline3217 involves averaging over scattering directions to evaluate the linear and angular momentum transport by the scattered wave. This evaluation requires appropriate choices of the parameters ICTHM and IPHM - see §21.



Bruce Draine
Thu Aug 10 09:34:16 EDT 2000