Hydromagnetic turbulence in cometary plasmas

An instability associated with the magnetosonic wave driven unstable due to coupling with electron and ion drift modes has been considered as a potential source for driving the hydromagnetic turbulence observed in Giacobini-Zinner (G-Z) Cometary plasma. The instability has good growth rate for propagation perpendicular to plasma inhomogeneities and exists for all wave numbers. The wave period for waves propagating perpendicular to the gradients is about a few times ion-gyroperiod and higher values of plasma beta ( _e lead to stronger instability.     

Magnetic fields in protoplanetary disks

Magnetic fields likely play a key role in the dynamics and evolution of protoplanetary disks. They have the potential to efficiently transport angular momentum by MHD turbulence or via the magnetocentrifugal acceleration of outflows from the disk surface. Magnetically-driven mixing has implications for disk chemistry and evolution of the grain population, and the effective viscous response of the disk determines whether planets migrate inwards or outwards. However, the weak ionisation of protoplanetary disks means that magnetic fields may not be able to effectively couple to the matter. I examine the magnetic diffusivity in a minimum solar nebula model and present calculations of the ionisation equilibrium and magnetic diffusivity as a function of height from the disk midplane at radii of 1 and 5 AU. Dust grains tend to suppress magnetic coupling by soaking up electrons and ions from the gas phase and reducing the conductivity of the gas by many orders of magnitude. However, once grains have grown to a few microns in size their effect starts to wane and magnetic fields can begin to couple to the gas even at the disk midplane. Because ions are generally decoupled from the magnetic field by neutral collisions while electrons are not, the Hall effect tends to dominate the diffusion of the magnetic field when it is able to partially couple to the gas, except at the disk surfaces where the low density of neutrals permits the ions to remain attached to the field lines. For a standard population of 0.1 μm grains the active surface layers have a combined column Σ_active≈2 g cm⁻² at 1 AU; by the time grains have aggregated to 3 μm, Σ_active≈80 g cm⁻². Ionisation in the active layers is dominated by stellar X-rays. In the absence of grains, X-rays maintain magnetic coupling to 10% of the disk material at 1 AU (i.e. Σ_active≈150 g cm⁻²). At 5 AU the Σ_active≈Σ_total once grains have aggregated to 1 μm in size.     

Cavitation by Nonlinear Interaction Between Inertial Alfvén Waves and Magnetosonic Waves in Low Beta Plasmas

This paper presents the model equations governing the nonlinear interaction between dispersive Alfvén wave (DAW) and magnetosonic wave in the low-β plasmas (β≪m _e/m _i; known as inertial Alfvén waves (IAWs); here \upbeta = 8pn₀T /B₀²\upbeta = 8\pi n_{0}T /B_{0}^{2} is thermal to magnetic pressure, n ₀ is unperturbed plasma number density, T(=T _e≈T _i) represents the plasma temperature, and m _e(m _i) is the mass of electron (ion)). This nonlinear dynamical system may be considered as the modified Zakharov system of equations (MZSE). These model equations are solved numerically by using a pseudo-spectral method to study the nonlinear evolution of density cavities driven by IAW. We observed the nonlinear evolution of IAW magnetic field structures having chaotic behavior accompanied by density cavities associated with the magnetosonic wave. The relevance of these investigations to low-β plasmas in solar corona and auroral ionospheric plasmas has been pointed out. For the auroral ionosphere, we observed the density fluctuations of ∼ 0.07n ₀, consistent with the FAST observation reported by Chaston et al. (Phys. Scr. T84, 64, 2000). The heating of the solar corona observed by Yohkoh and SOHO may be produced by the coupling of IAW and magnetosonic wave via filamentation process as discussed here.     

Viscous-resistive ADAF with a general large-scale magnetic field

We have studied the structure of hot accretion flow bathed in a general large-scale magnetic field. We have considered magnetic parameters , where are the Alfvén sound speeds in three direction of cylindrical coordinate (r,φ,z). The dominant mechanism of energy dissipation is assumed to be the magnetic diffusivity due to turbulence and viscosity in the accretion flow. Also, we adopt a more realistic model for kinematic viscosity (ν=αc _s H), with both c _s and H as a function of magnetic field. As a result in our model, the kinematic viscosity and magnetic diffusivity (η=η ₀ c _s H) are not constant. In order to solve the integrated equations that govern the behavior of the accretion flow, a self-similar method is used. It is found that the existence of magnetic resistivity will increase the radial infall velocity as well as sound speed and vertical thickness of the disk. However the rotational velocity of the disk decreases by the increase of magnetic resistivity. Moreover, we study the effect of three components of global magnetic field on the structure of the disk. We found out that the radial velocity and sound speed are Sub-Keplerian for all values of magnetic field parameters, but the rotational velocity can be Super-Keplerian by the increase of toroidal magnetic field. Also, Our numerical results show that all components of magnetic field can be important and have a considerable effect on velocities and vertical thickness of the disk.     

Magnetically-aided evolution of protoplanetary disks

We investigate the global evolution of a turbulent protoplanetary disk incorporating the effects of Maxwell stress due to a large-scale magnetic field permeating the disk. A magnetic field is produced continuously by an dynamo and the resultant Maxwell stress assists the viscous stress in p roviding the means for disk evolution. The most striking feature of magnetized disk evolution is the presence of the surface density bulge located in the magnetic gap, the region of the disk where the degree of ionization is too low to allow for coupli ng between the magnetic field and the gas. The bulge persists for a time of the order of 10⁵–10⁶ yr. The presence and persistence of the surface density bulge may have important implications for the process of planet formation and the overall characteristics of resultant planetary systems.Operated by USRA under contract No. NASW-4574 with NASA.     

Limits for the Firehose Instability in Space Plasmas

Electromagnetic instabilities in high-β plasmas, where β is the ratio of the kinetic plasma energy to the magnetic energy, have a broad range of astrophysical applications. The presence of temperature anisotropies T _∥ /T _⊥ >1 (where ∥ and ⊥ denote directions relative to the background magnetic field) in solar flares and the solar wind is sustained by the observations and robust acceleration mechanisms that heat plasma particles in the parallel direction. The surplus of parallel kinetic energy can excite either the Weibel-like instability (WI) of the ordinary mode perpendicular to the magnetic field or the firehose instability (FHI) of the circularly polarized waves at parallel propagation. The interplay of these two instabilities is examined. The growth rates and the thresholds provided by the kinetic Vlasov – Maxwell theory are compared. The WI is the fastest growing one with a growth rate that is several orders of magnitude larger than that of the FHI. These instabilities are however inhibited by the ambient magnetic field by introducing a temperature anisotropy threshold. The WI admits a larger anisotropy threshold, so that, under this threshold, the FHI remains the principal mechanism of relaxation. The criteria provided here by describing the interplay of the WI and FHI are relevant for the existence of these two instabilities in any space plasma system characterized by an excess of parallel kinetic energy.     

Multiresolution Analysis of Active Region Magnetic Structure and its Correlation with the Mount Wilson Classification and Flaring Activity

Two different multiresolution analyses are used to decompose the structure of active-region magnetic flux into concentrations of different size scales. Lines separating these opposite polarity regions of flux at each size scale are found. These lines are used as a mask on a map of the magnetic field gradient to sample the local gradient between opposite polarity regions of given scale sizes. It is shown that the maximum, average, and standard deviation of the magnetic flux gradient for α,β,β γ, and β γ δ active-regions increase in the order listed, and that the order is maintained over all length scales. Since magnetic flux gradient is strongly linked to active-region activity, such as flares, this study demonstrates that, on average, the Mt. Wilson classification encodes the notion of activity over all length scales in the active-region, and not just those length scales at which the strongest flux gradients are found. Further, it is also shown that the average gradients in the field, and the average length-scale at which they occur, also increase in the same order. Finally, there are significant differences in the gradient distribution, between flaring and non-flaring active regions, which are maintained over all length scales. It is also shown that the average gradient content of active-regions that have large flares (GOES class “M” and above) is larger than that for active regions containing flares of all flare sizes; this difference is also maintained at all length scales. All of the reported results are independent of the multiresolution transform used. The implications for the Mt. Wilson classification of active-regions in relation to the multiresolution gradient content and flaring activity are discussed.     

Interaction of particles with ion cyclotron waves and magnetosonic waves. Observations from GEOS 1 and GEOS 2

Coordinated observations involving ion composition, thermal plasma, energetic particle, and ULF magnetic field data from GEOS 1 and 2 often reveal the presence of electromagnetic ion cyclotron and magnetosonic waves, which are distinguished by their respective polarization characteristics and frequency spectra. The ion cyclotron waves are identified by a magnetic field perturbation that lies in a plane perpendicular to the Earth's magnetic field B₀ and propagate along B₀. They are associated with the abundance of cold He⁺ in the presence of anisotropic pitch angle distributions of ions having energies E > 20 keV, and were observed at frequencies near the He⁺ gyrofrequency. The magnetosonic waves are characterized by a magnetic field perturbation parallel to B₀ and thus seem to be propagating perpendicular to the Earth's magnetic field. They often occur at harmonics (not always including the fundamental) at the proton gyrofrequency and are associated with phase-space-density distributions that peak at energies E ～ 5–30 keV and at a pitch angle of 90°. Such a ring-like distribution is shown to excite instability in the magnetosonic mode near harmonics of the proton gyrofrequency. Magnetosonic waves are associated in other cases with sharp spatial gradients in energetic ion intensity. Such gradients are encountered in the early afternoon sector (as a consequence of the drift shell distortion caused by the convection electric field) and could likewise constitute a source of free energy for plasma instabilities.     

Simulations of jet formation from a magnetized accretion disk

Axisymmetric magnetohydrodynamic (MHD) simulations have been made of the formation of jets from a Keplerian disk threaded by a magnetic field. The disk is treated as a boundary condition, where matter with high specific entropy is ejected with a Keplerian azimuthal speed and a poloidal speed less than the slow magnetosonic velocity, and where boundary conditions on the magnetic fields correspond to a highly conducting disk. Initially, the space above the disk, the corona, is filled with high specific entropy plasma in the thermal equilibrium in the gravitational field of the central object. The initial magnetic field is poloidal and is represented by the superposition of the fields of monopoles located below the plane of the disk.The rotation of the disk twists the initial poloidal magnetic field lines, and this twist propagates into the corona pushing matter into jet-like outflow in a cylindrical region. After the first switch-on wave, which originates during the first rotation period of the inner radius of the disk, the matter outflowing from the disk starts to flow and accelerate in thez-direction owing to both the magnetic and pressure gradient forces. The flow accelerates through the slow magnetosonic and Alfvén surfaces and at larger distances through the fast magnetosonic surface. The flow velocity of the jet is approximately parallel to thez-axis, with the collimation mainly a result of the pinching force of the toroidal magnetic field. The energy flux of the flow increases with increasing magnetic field strength on the disk. Some of the cases studied have been run for long times, 60 rotation periods of the inner radius of the disk, and show indications of approaching a stationary state.     

On the pulsational-radial oscillations instability of isothermal magnetized accretion disk with self-gravity

This is the first paper to consider the effects of both magnetic field and self-gravity on the pulsational instability. Our main new results are that the self-gravity enhances the instability of the magneto-acoustic mode in the outer disk strongly, and also affects the instability in the inner disk, but stabilized the viscous mode. The effect of self-gravity is much greater than that of magnetic field in the outer disk, while the effect of magnetic field on the instability is weaker than that in the previous work's (Wuet al., 1995; Yanget al., 1995), in which the self-gravity has not been considered. Finally, we discuss our results.     

Evolutionary processes in protoplanetary accretion disks: the propagation of axisymmetric shock waves

In this paper we investigate both the global and the local hydrodynamics of axisymmetric accretion disks around young stellar objects under the simultaneous action of viscosity, self-gravity and pressure forces. For simplicity, we take for the global model a polytropic equation of state, make the infinitely thin disk approximation and characterize the surface density and temperature profiles in the disk as power laws in the radial distance r from the protostar. We solve the problem of the general density profile of a Keplerian disk showing that self-gravity could not be an important factor for the fast formation of the rocky cores of giant gaseous planets in our solar system. Under the hypothesis that the unperturbed rotation curve of the disk is nearly Keplerian throughout the radial extent, we can estimate with our polytropic model a lower limit for the resulting masses M_d(r) of stable disks up to 100 AU. These masses are in the range of the so-called minimum mass solar nebular (_d/M_s ≈ 0.01–0.02).By adopting a simplified viscosity model, where the height-integrated turbulent dynamical viscosity ν is a function of the surface density σ like η ∝ σ^Γ, we derive in the local shearing sheet model linearized evolution equations for small density perturbations describing both a diffusion process and the propagation of acoustic density waves. We solve a special initial value problem and calculate the appropriate Green's function. The analytical solutions so obtained describe in the case Γ < 0 the successive formation of quasi-stationary ring-shaped density structures in a disk with a definite mode of maximum instability, whereas in the case Γ > Γ_c the density wave equation describes the propagation of an “overstable” ring-shaped acoustic density wavelet to the outer ranges of the accretion disk. Whereas the group velocity of the wave packet is subsonic, the phase velocities of individual wave crests in the wave packet are supersonic. The mode of maximum instability, the growth rate and the number of growing waves in the wavelet are controlled by Γ and α. Our present knowledge concerning turbulent viscosity in protoplanetary disks is not sufficient to decide whether or not the case Γ > Γ_c is realized.The suggested structuring processes in the linear theory should initiate in the non-linear regime the formation of narrow ring-shaped density shock waves moving through the protoplanetary disk. These non-linear waves could produce extremely spatially and temporally heterogeneous temperature regions in the disk. We speculate that ring-shaped density waves, excited by inner boundary conditions and which have dominated the disk's evolution at early times, are responsible both for the fast growth of dust to planetesimals and at least for the rapid accretion of the rocky cores of giant gaseous planets in the protoplanetary accretion disk (shock wave trigger hypothesis). We derive provisional scaling rules for planetary systems regarding the spacing of orbits as a function of the mass ratio of the protoplanetary disk to the protostar. However, further analytical work and linear as well as nonlinear numerical simulations of density waves excited by inner boundary conditions are needed to consolidate the results and speculations of our linear wave mechanics in the future.     

Self-similar dynamics of a magnetized polytropic gas

In broad astrophysical contexts of large-scale gravitational collapses and outflows and as a basis for various further astrophysical applications, we formulate and investigate a theoretical problem of self-similar magnetohydrodynamics (MHD) for a non-rotating polytropic gas of quasi-spherical symmetry permeated by a completely random magnetic field. Within this framework, we derive two coupled nonlinear MHD ordinary differential equations (ODEs), examine properties of the magnetosonic critical curve, obtain various asymptotic and global semi-complete similarity MHD solutions, and qualify the applicability of our results. Unique to a magnetized gas cloud, a novel asymptotic MHD solution for a collapsing core is established. Physically, the similarity MHD inflow towards the central dense core proceeds in characteristic manners before the gas material eventually encounters a strong radiating MHD shock upon impact onto the central compact object. Sufficiently far away from the central core region enshrouded by such an MHD shock, we derive regular asymptotic behaviours. We study asymptotic solution behaviours in the vicinity of the magnetosonic critical curve and determine smooth MHD eigensolutions across this curve. Numerically, we construct global semi-complete similarity MHD solutions that cross the magnetosonic critical curve zero, one, and two times. For comparison, counterpart solutions in the case of an isothermal unmagnetized and magnetized gas flows are demonstrated in the present MHD framework at nearly isothermal and weakly magnetized conditions. For a polytropic index γ=1.25 or a strong magnetic field, different solution behaviours emerge. With a strong magnetic field, there exist semi-complete similarity solutions crossing the magnetosonic critical curve only once, and the MHD counterpart of expansion-wave collapse solution disappears. Also in the polytropic case of γ=1.25, we no longer observe the trend in the speed-density phase diagram of finding infinitely many matches to establish global MHD solutions that cross the magnetosonic critical curve twice.      

Fracton Excitations as a Driving Mechanism for the Self-Organized Dynamical Structuring in the Solar Wind

Structural properties of the interplanetary magnetic field (IMF) are discussed. Our main interest is concentrated on the dynamical structuring mechanisms associated with the dominant role of the wave processes in the solar wind. We argue that the IMF possibly reveals the self-organized clustering driven by the low-frequency magnetosonic waves. It is shown that the self-organized geometry of the IMF is a fractal, a specific object having a number of unusual topological features; this fractal geometry is self-consistently generated by the allowed magnetosonic modes. To give an accurate treatment of waves on fractals, we propose an unconventional approach based on the wave equation with the generalized, fractional time derivative. The allowed magnetosonic modes are then defined as the generalized "resonance" solutions to the fractional wave equation and termed "fractons", vibrational excitations of fractal objects. We found that the self-organized fractal geometry of the IMF as maintained by the fractons could be described by the value of the Hausdorff fractal dimension D≈ 4/3. Convection of the IMF fractal structures by a spacecraft observer is shown to result in the power-law behavior of the Fourier energy density spectrum of the in situobserved IMF turbulence, P(f) ∝ f ^−α, with the characteristic slope α ≈ 5/3. This revised version was published online in July 2006 with corrections to the Cover Date.     

Resonant Transition Radiation and Solar Radio Bursts

This paper presents general relations for the intensity of the resonant transition radiation (RTR) and their detailed analysis. This analysis shows that the spectrum amplitude of the x-mode at some frequencies for high-energy electrons can grow with the magnetic field increase in some interval from zero value; it can even dominate over that for the o-mode. With further magnetic field increase, the intensity of the RTR x-mode decreases in comparison with the intensity of the o-mode and this decrease is higher for higher velocities of energetic electrons. The polarization of the RTR depends on the velocity of energetic electrons, too. For velocities lower than some velocity limit v<v _i the RTR emission is unpolarized in a broad interval of magnetic field intensities in the radio source. For reasonable values of indices of the power-law distribution functions of energetic electrons, the RTR is broadband in frequencies (df/f≈0.2−0.4). Furthermore, we show various dependencies of the RTR and its spectral characteristics. Assuming the same radio flux of the transition radiation and the gyro-synchrotron one at the Razin frequency, we estimate the limit magnetic field in the radio source of the transition radiation. Then, we analyze possible sources of small-scale inhomogeneities (thermal density fluctuations, Langmuir and ion-sound waves), which are necessary for the transition radiation. Although the small-scale inhomogeneities connected with the Langmuir waves lead to the plasma radiation, which is essentially stronger than RTR, the inhomogeneities of the ion-sound waves are suitable for the RTR without any other radiation. We present the relations describing the RTR for anisotropic distribution functions of fast electrons. We consider the distribution functions of fast electrons in the form of the Legendre polynomials which depend on the pitch-angle. We analyze the influence of the degree of the anisotropy (an increase of the number of terms in the Legendre polynomial) on spectral characteristics of the RTR. A comparison with previous studies is made. As an example of the use of the derived formulas for the RTR, the 24 December 1991 event is studied. It is shown that the observed decimetric burst can be generated by the RTR in the plasma with the density inhomogeneities at the level 〈ΔN ²〉/N ²=2.5⋅10⁻⁵.     

Polar Coronal Structures and the Global Magnetic Field Evolution Through the Cycle

We compare the shape and position of some plasma formations visible in the polar corona with the cyclic evolution of the global magnetic field. The first type of object is polar crown prominences. A two-fold decrease of the height of polar crown prominences was found during their poleward migration from the middle latitudes to the poles before a polar magnetic field reversal. The effect could be assigned to a decrease of the magnetic field scale. The second type of object is the polar plumes, ray like structures that follow magnetic field lines. Tangents to polar ray structures are usually crossed near some point, “a magnetic focus,” below the surface. The distance q between the focus and the center of the solar disk changes from the maximum value about 0.65 R _⊙ at solar minimum activity to the minimum value about 0.45 R _⊙ at solar maximum. At first glance this behaviour seems to be contrary to the dynamics of spherical harmonics of the global magnetic field throughout a cycle. We believe that the problem could be resolved if one takes into account not only scale changes in the global magnetic field but also the phase difference in the cyclic variations of large-scale and small-scale components of the global field.     

Modeling Magnetic Tower Jets in the Laboratory

The twisting of magnetic fields threading an accretion system can lead to the generation on axis of toroidal field loops. As the magnetic pressure increases, the toroidal field inflates, producing a flow. Collimation is due to a background corona, which radially confines this axially growing “magnetic tower”. We investigate the possibility of studying in the laboratory the dynamics, confinement and stability of magnetic tower jets. We present two-dimensional resistive magnetohydrodynamic simulations of radial arrays, which consist of two concentric electrodes connected radially by thin metallic wires. In the laboratory, a radial wire array is driven by a 1 MA current which produces a hot, low density background plasma. During the current discharge a low plasma beta (β < 1), magnetic cavity develops in the background plasma (β is the ratio of thermal to magnetic pressure). This laboratory magnetic tower is driven by the magnetic pressure of the toroidal field and it is surrounded by a shock envelope. On axis, a high density column is produced by the pinch effect. The background plasma has >rsim1, and in the radial direction the magnetic tower is confined mostly by the thermal pressure. In contrast, in the axial direction the pressure rapidly decays and an elongated, well collimated magnetic-jet develops. This is later disrupted by the development of m = 0 instabilities arising in the axial column.     

Studies on some properties of coronal mass ejections based on angular width

The paper investigates the effects of thermal conductivity and non-uniform magnetic field on the gravitational instability of a non-uniformly rotating infinitely extending axisymmetric cylinder in a homogeneous heat conducting medium. The non-uniform rotation and magnetic field are supposed to act along θ and z directions of the cylinder. It is found that the gravitational instability of this general problem is determined by the same criterion as obtained by Dhiman and Dadwal (Astrophys. Space Sci. 325(2):195–200, 2010) for the self-gravitating isothermal medium in the presence of non-uniform rotation and magnetic field with the only difference that adiabatic sound velocity is now replaced by the isothermal sound velocity. It is found that the thermal conductivity has stabilizing effect on the onset of gravitational instability. Further, the stabilizing/destabilizing effect of the non-uniform magnetic field on the gravitational instability of heat conducting medium has been discussed and is illustrated by considering some special forms of the basic magnetic fields.     

Triggered planet formation in young stellar clusters

Conventional planet formation models via coagulation of planetesimals require timescales in the range of several 10 or even 100 Myr in the outer regions of a protoplanetary disk. But according to observational data, the lifetime of a protoplanetary disk is limited to about 6 Myr. Therefore the existence of Uranus and Neptune poses a problem. Planet formation via gravitational instability may be a solution for this discrepancy. We present a parameter study of the possibility of gravitationally triggered disk instability. Using a restricted N‐body model which allows for a survey of an extended parameter space, we show that a passing dwarf star with a mass between 0.1 and 1 M_⊙ can probably induce gravitational instabilities in the pre‐planetary solar disk for prograde passages with minimum separations below 80‐170 AU. Inclined and retrograde encounters lead to similar results but require slightly closer passages. Such encounter distances are quite likely in young moderately massive star clusters. The induced gravitational instabilities may lead to enhanced planetesimal formation in the outer regions of the protoplanetary disk, and could therefore be relevant for the formation of Uranus and Neptune. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A submillimeter view of protoplanetary dust disks

We present some results from our submillimeter single-dish and aperture synthesis imaging surveys of protoplanetary disks using the JCMT, CSO, and Submillimeter Array (SMA) on Mauna Kea, Hawaii. Employing a simple disk model, we simultaneously fit the spectral energy distributions and spatially resolved submillimeter continuum emission from our SMA survey to constrain disk structure properties, including surface density profiles and sizes. The typical disk structure we infer is consistent with a fiducial accretion disk model with a viscosity parameter α≈0.01. Combined with a large, multiwavelength single-dish survey of similar disks, we show how these observations provide evidence for significant grain growth and rapid evolution in the outer regions of disks, perhaps due to an internal photoevaporation process. In addition, we discuss SMA observations of the disks in the Orion Trapezium (proplyds) in the context of disk evolution in a more extreme environment.     

Dust ejection from planetary bodies by temperature gradients: Laboratory experiments

Laboratory experiments show that dusty bodies in a gaseous environment eject dust particles if they are illuminated. We find that even more intense dust eruptions occur when the light source is turned off. We attribute this to a compression of gas by thermal creep in response to the changing temperature gradients in the top dust layers. The effect is studied at a light flux of 13 kW/m² and 1 mbar ambient pressure. The effect is applicable to protoplanetary disks and Mars. In the inner part of protoplanetary disks, planetesimals can be eroded especially at the terminator of a rotating body. This leads to the production of dust which can then be transported towards the disk edge or the outer disk regions. The generated dust might constitute a significant fraction of the warm dust observed in extrasolar protoplanetary disks. We estimate erosion rates of about 1 kg s^?1 for 100 m parent bodies. The dust might also contribute to subsequent planetary growth in different locations or on existing protoplanets which are large enough not to be susceptible to particle loss by light induced ejection. Due to the ejections, planetesimals and smaller bodies will be accelerated or decelerated and drift outward or inward, respectively. The effect might also explain the entrainment of dust in dust devils on Mars, especially at high altitudes where gas drag alone might not be sufficient.