Effect of solar radiation on interplanetary dust grains

Let , , be cylindrical polar coordinates in a frame of reference, centered on the Sun, that is aligned with the orbital plane of the dust grain, as described in Section I.1. Let and be the (relative) position and velocity of the grain, respectively. Consider a photon emitted by the Sun. Let and be the photon's energy and momentum, respectively, in the heliocentric frame. Let and be the corresponding quantities in the dust grain's instantaneous rest frame. According to standard relativistic theory (Rindler 1977),

(10.160) | ||||||

and | (10.161) |

where . Assuming that the grain is moving non-relativistically, so that , the previous expressions reduce to

(10.162) | ||||

and | (10.163) |

respectively, to first order in . The corresponding inverse transforms are

(10.164) | ||||

and | (10.165) |

To first order in , the previous four expressions also apply to energy and momentum measured per unit time.

Let and be the electromagnetic energy and momentum, respectively, absorbed by the grain per unit time in its instantaneous rest frame. Let and be the corresponding quantities emitted by the grain. Now, we are assuming that

(10.166) |

In other words, all energy absorbed by the grain is immediately re-emitted. We also assuming that

(10.167) |

In other words, the emitted radiation pattern is isotropic in the grain's instantaneous rest frame, and consequently carries off zero net momentum. It follows that the grain's net rate of energy gain in its instantaneous rest frame is

whereas the corresponding rate of momentum gain is

Transforming to the heliocentric frame, we obtain

(10.170) |

where use has been made of Equations (10.165), (10.168), and (10.169). However, according to Equation (10.163),

(10.171) |

Furthermore, to first order in ,

(10.172) | ||||||

and | (10.173) |

where is the mean energy of a photon emitted by the Sun, the corresponding mean momentum (Rindler 1977), the local photon number density, and the cross-sectional area of the grain (normal to ). It follows that, to first order in ,

(10.174) |

Of course, is the local solar electromagnetic energy flux, where is the

In the heliocentric frame, the net force per unit mass acting on the grain is

(10.175) |

where is the grain's mass. It follows from the previous analysis that

where

(10.177) |

and

(10.178) |

[A more rigorous derivation of the previous two equations is given in Robertson 1937 and Klačka 1993.] Here, is the radiation pressure force

where use has been made of Equations (I.3), (I.7), and (I.10).

Assuming that the total radiation pressure force, (10.176), is small compared to the force of gravitational attraction between the Sun and the dust grain--and can, thus, be treated as a perturbation--the grain's orbit can be modeled as Keplerian ellipse whose six elements evolve slowly in time under the influence of the pressure. The six elements in question are chosen to be the major radius, , the mean anomaly at epoch, , the eccentricity, , the argument of the perigee, , the inclination (to the ecliptic plane), , and the longitude of the ascending node (measured with respect to the vernal equinox), . (See Section 4.12.) The evolution of these elements is governed by the Gauss planetary equations, (I.53)-(I.58). Now, in the perturbative limit, the evolution of the dust grain's orbital elements takes place on a timescale that is much longer than its orbital period. We can concentrate on this evolution, and filter out any relatively short-term oscillations in the elements, by averaging the Gauss planetary equations over an orbital period. A suitable orbit-average operator is

(10.182) |

Here, we have made use of the fact that is a constant of the unperturbed motion. (See Chapter 4.) Note that is the unperturbed mean orbital angular velocity. Moreover, , where is the solar mass. (Of course, we are ignoring the mass of the grain with respect to that of the Sun.) Making use of Equations (10.179)-(10.181), as well as some of the formulae appearing in Appendix I, the orbit-averaged Gauss planetary equations become

It can be seen that the radiation pressure force causes the grain's orbital major radius and eccentricity to both decrease monotonically in time [because the left-hand sides of Equations (10.183) and (10.185) are both negative.] On the other hand, radiation pressure has no effect on the either the location of the grain's perihelion point or the orientation of its orbital plane (because , and are all zero). Finally, given that the right-hand side of Equation (10.184) is negative, the radiation pressure force produces a slight reduction in the grain's mean orbital angular velocity: the modified angular velocity being . Note that, to lowest order in , the radial radiation pressure force, , averages out of the expressions for and . Hence, the reduction in the grain's orbital major radius and eccentricity is due to the Poynting-Robertson drag, , combined with the first-order component of . ( is similar in magnitude to the first-order component of .) On the other hand, the lowest-order component of makes the dominant contribution to the expression for (in fact, we have neglected any first-order contributions to this expression).

In order to be in the perturbative limit, the relative changes in the dust grain's orbital elements induced by the radiation pressure force in an orbital period must all be small. Because [given that , which follows because the grain is moving non-relativistically], this requirement yields , or

(10.189) |

which is equivalent to the requirement that the radiation pressure force be much smaller than the force of gravitational attraction between the grain and the Sun. Let us model the grain as a sphere of radius and mass density . It follows that and . Thus, the previous expression yields

(10.190) |

Given that , , , (Yoder 1995), and assuming that , we obtain

(10.191) |

We conclude that the radiation pressure force is perturbative for ``large'' dust grains whose radii are much greater than a micron. On the other hand, the radiation pressure force--whose dominant component is directed radially outward from the Sun--exceeds the force of gravitational attraction for ``small'' dust grains whose radii are less than, or of order, a micron. We would, therefore, expect radiation pressure to eventually expel small dust grains from the solar system.

For the case of a large dust grain in a circular orbit around the Sun, Equation (10.183) gives

(10.192) |

Hence, the time required for a grain whose orbit has an initial radius to spiral into the Sun is

(10.193) |

Again modeling the grain as a sphere of radius and mass density , we obtain

(10.194) |

Here, is measured in astronomical units, and in meters. Thus, a dust grain of radius , and initial orbital radius 1 AU, takes about years to spiral into the Sun.