Barnett relaxation, first described by Purcell in 1979, appears to play a major role in the alignment of grains with the interstellar magnetic field. In 1999, Lazarian and Draine proposed that Barnett relaxation and its relative, nuclear relaxation, can induce grains to flip. If this thermal flipping is rapid, then the dynamical effect of torques that are fixed relative to the grain body can be greatly reduced. To date, detailed studies of Barnett relaxation have been confined to grains exhibiting dynamic symmetry. In 2003, Weingartner argued that internal relaxation cannot induce flips in any grains, whether they exhibit dynamic symmetry or not. In this work, we develop approximate expressions for the dissipation rate and diffusion coefficient for Barnett relaxation. We revisit the issue of internally induced thermal flipping, finding that it cannot occur for grains with dynamic symmetry but does occur for grains lacking dynamic symmetry.

Collisions of gas particles with a drifting grain give rise to a mechanical torque on the grain. Recent work by Lazarian & Hoang showed that mechanical torques might play a significant role in aligning helical grains along the interstellar magnetic field direction, even in the case of subsonic drift. We compute the mechanical torques on 13 different irregular grains and examine their resulting rotational dynamics, assuming steady rotation about the principal axis of greatest moment of inertia. We find that the alignment efficiency in the subsonic drift regime depends sensitively on the grain shape, with more efficient alignment for shapes with a substantial mechanical torque even in the case of no drift. The alignment is typically more efficient for supersonic drift. A more rigorous analysis of the dynamics is required to definitively appraise the role of mechanical torques in grain alignment.

The degree to which interstellar grains align with respect to the interstellar magnetic field depends on disaligning as well as aligning mechanisms. For decades, it was assumed that disalignment was due primarily to the random angular impulses a grain receives when colliding with gas-phase atoms. Recently, a new disalignment mechanism has been considered, which may be very potent for a grain that has a time-varying electric dipole moment and drifts across the magnetic field. We provide quantitative estimates of the disalignment times for silicate grains with size > approximately 0.1 micron. These appear to be shorter than the time-scale for alignment by radiative torques, unless the grains contain superparamagnetic inclusions.

In interstellar dust grains, internal processes dissipate rotational kinetic energy. The dissipation is accompanied by thermal fluctuations, which transfer energy from the vibrational modes to rotation. Together, these processes are known as internal relaxation. For the past several years, internal relaxation has been thought to give rise to thermal flipping, with profound consequences for grain alignment theory. I show that thermal flipping is not possible in the limit that the inertia tensor does not vary with time.

Radiative torques, due to the absorption and scattering of starlight, are thought to play a major role in the alignment of grains with the interstellar magnetic field. The absorption of radiation also gives rise to recoil torques, associated with the photoelectric effect and photodesorption. The recoil torques are much more difficult to model and compute than the direct radiative torque. Here, we consider the relatively simple case of a spheroidal grain. Given our best estimates for the photoelectric yield and other relevant grain physical properties, we find that the recoil torques contribute at the 10% level or less compared with the direct radiative torque. We recommend that the recoil torques not be included in models of radiation-driven grain alignment at this time. However, additional experimental characterization of the surface properties and photoelectric yield for sub-micron grains is needed to better quantify the magnitude of these torques.

We present the first systematic investigation of the [NeV] (14um/24um) and [SIII] (18um/33um) infrared line flux ratios, traditionally used to estimate the density of the ionized gas, in a sample of 41 Type 1 and Type 2 active galactic nuclei (AGNs) observed with the Infrared Spectrograph on board Spitzer. The majority of galaxies with both [NeV] lines detected have observed [NeV] line flux ratios consistent with or below the theoretical low density limit, based on calculations using currently available collision strengths and ignoring absorption and stimulated emission. We find that Type 2 AGNs have lower line flux ratios than Type 1 AGNs and that all of the galaxies with line flux ratios below the low density limit are Type 2 AGNs. We argue that differential infrared extinction to the [NeV] emitting region due to dust in the obscuring torus is responsible for the ratios below the low density limit and we suggest that the ratio may be a tracer of the inclination angle of the torus to our line of sight. Because the temperature of the gas, the amount of extinction, and the effect of absorption and stimulated emission on the line ratios are all unknown, we are not able to determine the electron densities associated with the [NeV] line flux ratios for the objects in our sample. We also find that the [SIII] emission from the galaxies in our sample is extended and originates primarily in star forming regions. Since the emission from low-ionization species is extended, any analysis using line flux ratios from such species obtained from slits of different sizes is invalid for most nearby galaxies.

Several mechanisms have been proposed to explain the alignment of grains with the interstellar magnetic field, including paramagnetic dissipation, radiative torques, and supersonic gas-grain streaming. These must compete with disaligning processes, including randomly directed torques arising from collisions with gas atoms. I describe a novel disalignment mechanism for grains that have a time-varying electric dipole moment and that drift across the magnetic field. Depending on the drift speed, this mechanism may yield a much shorter disalignment timescale than that associated with random gas atom impacts. For suprathermally rotating grains, the new disaligning process may be more potent for carbonaceous dust than for silicate dust. This could result in efficient alignment for silicate grains but poor alignment for carbonaceous grains.

Photoelectric emission from dust plays an important role in grain charging and gas heating. To date, detailed models of these processes have focused primarily on grains exposed to soft radiation fields. We provide new estimates of the photoelectric yield for neutral and charged carbonaceous and silicate grains, for photon energies exceeding 20 eV. We include the ejection of electrons from both the band structure of the material and the inner shells of the constituent atoms, as well as Auger and secondary electron emission. We apply the model to estimate gas heating rates in planetary nebulae and grain charges in the outflows of broad absorption line quasars. For these applications, secondary emission can be neglected; the combined effect of inner shell and Auger emission is small, though not always negligible. Finally, we investigate the survivability of dust entrained in quasar outflows. The lack of nuclear reddening in broad absorption line quasars may be explained by sputtering of grains in the outflows.

Ever since the pioneering study of Spitzer, it has been widely recognized that grains play an important role in the heating and cooling of photo-ionized environments. A detailed code is necessary to model grains in a photo-ionized medium. In this paper we will use the spectral synthesis code Cloudy for this purpose. The grain model of this code has recently undergone a comprehensive upgrade. Part of this is the newly developed hybrid grain charge model, which will be described in detail. This model allows discrete charge states of very small grains to be modelled accurately while simultaneously avoiding the overhead of fully resolving the charge distribution of large grains, thus making the model both accurate and computationally efficient. A comprehensive comparison with the fully resolved charge state models of Weingartner & Draine (2001) shows that the agreement is very satisfactory for realistic size distributions. The effect of the grain size distribution on the line emission from photo-ionized regions is studied by taking standard models for an H II region and a planetary nebula and adding a dust component to the models with varying grain size distributions. A comparison of the models shows that varying the size distribution has a dramatic effect on the emitted spectrum, and affects the ionization balance as well. These results clearly demonstrate that the grain size distribution is an important parameter in photo-ionization models.

Observations of infrared emission from AGN show that grains exist over a broad range of distances from the central object, extending to the point where they are destroyed by sublimation. These ~1000 K grains produce much of the 1 micron continuum. In this region closest to the central object there must be a gaseous component associated with the hot grains. This paper employs a state of the art grain model and shows that the gas must be very hot, with temperatures in the neighborhood of 10^6 K. The dusty component has a covering factor of roughly 50% and so this region also reprocesses much of the total x-ray emission. Our explicit models of the IR through x-ray spectral energy distributions allow the x-ray component to be predicted from IR observations. We are creating a grid of such predictions and will make them available as an XSPEC add-in, allowing this spectral component to be included in quantitative modeling of AGN spectra.

In this paper we present results obtained with the new grain code in Cloudy which underline the strong effect of photo-electric heating by grains in photo-ionized regions. We study the effect that the distribution of grain sizes has on the magnitude of the effect, and show that this effect is nothing short of dramatic. This makes the grain size distribution an important parameter in modeling of photo-ionized regions such as H II regions and planetary nebulae.

Micron size extrasolar dust particles have been convincingly detected by satellites. Larger extrasolar meteoroids (5-35 microns) have most likely been detected by ground based radar at Arecibo and New Zealand. We present estimates of the minimum detectable particle sizes and collecting areas for both radar systems. We show that particles larger than about 10 microns can propagate for tens of parsecs through the interstellar medium, opening up the possibility that ground based radar systems can detect AGB stars, young stellar objects such as T Tauri stars, and debris disks around Vega-like stars. We provide analytical and numerical estimates of the ejection velocity in the case of a debris disk interacting with a Jupiter mass planet. We provide rough estimates of the flux of large micrometeoroids from all three classes of sources. Current radar systems are unlikely to detect significant numbers of meteors from debris disks such as Beta Pictoris. However, we suggest improvements to radar systems that should allow for the detection of multiple examples of all three classes.

Hubble Space Telescope images of MCG-6-30-15 show a dust lane crossing the galaxy just below the nucleus. In this paper, we argue that this dust lane is responsible for the observed reddening of the nuclear emission and the Fe I edge hinted at in the Chandra spectrum of MCG-6-30-15. We further suggest that the gas within the dust lane can comprise much of the low ionization component (i.e., the one contributing the O VII edge) of the observed warm absorber. Moreover, placing the warm absorbing material at such distances (hundreds of pc) can account for the small outflow velocities of the low ionization absorption lines as well as the constancy of the O VIII edge. Photoionization models of a dusty interstellar gas cloud (with a column appropriate for the reddening toward MCG-6-30-15) using a toy Seyfert 1 spectral energy distribution show that it is possible to obtain a significant O VII edge (\tau~0.2) if the material is ~150 pc from the ionizing source. For MCG-6-30-15, such a distance is consistent with the observed dust lane. The current data on MCG-6-30-15 is unable to constrain the dust composition within the warm absorber. Astronomical silicate is a viable candidate, but there are indications of a very low O abundance in the dust, which is inconsistent with a silicate origin. If true, this may indicate that there were repeated cycles of grain destruction and growth from shocks in the interstellar medium of MCG-6-30-15. Pure iron grains are an unlikely dust constituent due to the limit on their abundance in the Galaxy, yet they cannot be ruled out. The high column densities inferred from the highly ionized zone of the warm absorber implies that this gas is dust-free.

In the previous papers in this series, we found that radiative torques can play a major role in the alignment of grains with the interstellar magnetic field. Since the radiative torques can drive the grains to suprathermal rotational speeds, in previous work we made the simplifying assumption that the grain principal axis of greatest moment of inertia is always parallel to the grain angular momentum. This enabled us to describe many of the features of the grain dynamics. However, this assumption fails when the grains enter periods of thermal rotation, which occur naturally in the radiative torque alignment scenario. In the present paper, we relax this assumption and explore the consequences for the grain dynamics. We develop a treatment to follow the grain dynamics including thermal fluctuations and ``thermal flipping'', and show results for one illustrative example. By comparing with a treatment without thermal fluctuations, we see that inclusion of thermal fluctuations can lead to qualitative changes in the grain dynamics. In a future installment in this series, we will use the more complete dynamical treatment developed here to perform a systematic study of grain alignment by radiative torques.

Recently, Maiolino et al. (2001a, A&A, 365, 28) constructed a sample of active galactic nuclei for which both the reddening E(B-V) and the column density N_H to the nucleus could be determined. For most of the galaxies in their sample, they found that E(B-V)/N_H is substantially smaller than for the diffuse ISM of our Galaxy. They asserted that either the dust-to-gas ratio is lower than in the Galaxy or that the grains are so large that they do not extinct or redden efficiently in the optical. We show that there is no systematic increase in E(B-V) with N_H for the Maiolino et al. (2001a) galaxies, which suggests that the X-ray absorption and optical extinction occur in distinct media. Maiolino et al. (2001b, A&A, 365, 37) suggested that the observed lines of sight for the Maiolino et al. (2001a) galaxies pass through the ``torus'' that obscures the broad line region and nuclear continuum in Seyfert 2 galaxies and argued that the torus grains are larger than Galactic grains. There is no reason to believe that the lines of sight for these galaxies pass through the torus, since the observed column densities are lower than those typically observed in Seyfert 2 galaxies. We suggest instead that the X-ray absorption occurs in material located off the torus and/or accretion disk while the optical extinction occurs in material located beyond the torus. The X-ray absorbing material could either be dust-free or could contain large grains that do not extinct efficently in the optical. There is no conclusive evidence that the grains in active galactic nuclei are systematically larger than those in the diffuse ISM of our Galaxy.

Ever since the pioneering study of Spitzer, it has been widely recognized that grains play an important role in the heating and cooling of photo-ionized environments. This includes the diffuse ISM, as well as H II regions, planetary nebulae, and photo-dissociation regions. A detailed code is necessary to model grains in a photo-ionized medium since the interactions of grains with their environment include a host of microphysical processes, and their importance can only be judged by performing a complete simulation. In this paper we will use the spectral synthesis code Cloudy for this purpose. A comprehensive upgrade of the grain model has been recently incorporated in Cloudy, and certain aspects of this upgrade will be discussed. Special emphasis will be on the new grain charge model. We will consider in detail the physics of grains in both ionized and neutral environments, and will present a calculation of photo-electric heating and collisional cooling rates for a range of physical conditions and grain materials and for a range of grain sizes (including a realistic size distribution). We conclude with a brief discussion of the problems currently hampering progress in this field. The new grain model will be used to model the silicate emission in the Ney-Allen nebula, and will help us better understand the nature of the grains in that part of the Orion complex.

With the high-resolution spectroscopy now available in the optical and satellite UV, it is possible to determine the neutral/ionized column density ratios for several different elements in a single cloud. Assuming ionization equilibrium for each element, one can make several independent determinations of the electron density. For the clouds for which such an analysis has been carried out, these different estimates disagree by large factors, suggesting that some process (or processes) besides photoionization and radiative recombination might play an important role in the ionization balance. One candidate process is collisions of ions with dust grains. Making use of recent work quantifying the abundances of polycyclic aromatic hydrocarbon molecules and other grains in the interstellar medium, as well as recent models for grain charging, we estimate the grain-assisted ion recombination rates for several astrophysically important elements. We find that these rates are comparable to the rates for radiative recombination for conditions typical of the cold neutral medium. Including grain-assisted ion recombination in the ionization equilibrium analysis leads to increased consistency in the various electron density estimates for the gas along the line of sight to 23 Orionis. However, not all of the discrepancies can be eliminated in this way; we speculate on some other processes that might play a role. We also note that grain-assisted recombination of H+ and He+ leads to significantly lower electron fractions than usually assumed for the cold neutral medium.

We construct size distributions for carbonaceous and silicate grain populations in different regions of the Milky Way, LMC, and SMC. The size distributions include sufficient very small carbonaceous grains (including polycyclic aromatic hydrocarbon molecules) to account for the observed infrared and microwave emission from the diffuse interstellar medium. Our distributions reproduce the observed extinction of starlight, which varies depending upon the interstellar environment through which the light travels. As shown by Cardelli, Clayton & Mathis in 1989, these variations can be roughly parameterized by the ratio of visual extinction to reddening, R_V. We adopt a fairly simple functional form for the size distribution, characterized by several parameters. We tabulate these parameters for various combinations of values for R_V and b_C, the C abundance in very small grains. We also find size distributions for the line of sight to HD 210121, and for sightlines in the LMC and SMC. For several size distributions, we evaluate the albedo and scattering asymmetry parameter, and present model extinction curves extending beyond the Lyman limit.

Grains exposed to anisotropic radiation fields are subjected to forces due to the asymmetric photon-stimulated ejection of particles. These forces act in addition to the ``radiation pressure'' due to absorption and scattering. Here we model the forces due to photoelectron emission and the photodesorption of adatoms. The ``photoelectric'' force depends on the ambient conditions relevant to grain charging. We find that it is comparable to the radiation pressure when the grain potential is relatively low and the radiation spectrum is relatively hard. The calculation of the ``photodesorption'' force is highly uncertain, since the surface physics and chemsitry of grain materials are poorly understood at present. For our simple yet plausible model, the photodesorption force dominates the radiation pressure for grains with size >~0.1 micron exposed to starlight from OB stars. We find that the anisotropy of the interstellar radiation field is ~10% in the visible and ultraviolet. We estimate size-dependent drift speeds for grains in the cold and warm neutral media and find that micron-sized grains could potentially be moved across a diffuse cloud during its lifetime.

We model the photoelectric emission from and charging of interstellar dust and obtain photoelectric gas heating efficiencies as a function of grain size and the relevant ambient conditions. Using realistic grain size distributions, we evaluate the net gas heating rate for various interstellar environments, and find less heating for dense regions characterized by R_V=5.5 than for diffuse regions with R_V=3.1. We provide fitting functions which reproduce our numerical results for photoelectric heating and recombination cooling for a wide range of interstellar conditions. In a separate paper we will examine the implications of these results for the thermal structure of the interstellar medium. Finally, we investigate the potential importance of photoelectric heating in H II regions, including the warm ionized medium. We find that photoelectric heating could be comparable to or exceed heating due to photoionization of H for high ratios of the radiation intensity to the gas density. We also find that photoelectric heating by dust can account for the observed variation of temperature with distance from the galactic midplane in the warm ionized medium.

We consider the depletion of elements from the interstellar gas onto a population of very small dust grains. Adopting a grain model in which of order 4% of the cosmic C abundance is in grains with radii <= 10 Angstroms, we find that the rate of accretion onto these grains is adequately fast to account for the observed large depletions of elements like Ti, without invoking unreasonably high rates of mass transfer between interstellar phases or low grain destruction rates. If these grains are composed of arene rings, then only a limited number of metal atoms can be locked up in them. The depletion would be quenched when this limit is reached, unless there is a mechanism for transferring the metals to larger grains and refreshing the very small grain population, for example by grain coagulation and shattering in the diffuse ISM. If Fe depletes onto the very small grains, then for reasonable coagulation rates there is at least one metal atom per five C atoms in the very small grain population. Furthermore, approximately 60% of the cosmic Fe is associated with the carbonaceous grain population. It is unclear whether this scenario is compatible with observations. However, if there is another population of very small grains, with a large capacity for holding Fe atoms, it might be the sink for the most heavily depleted elements.

We discuss the forces on grains exposed to anisotropic radiation fields, including the usual ``radiation pressure'' force and also recoil forces due to photoemitted electrons, photodesorbed hydrogen atoms, and hydrogen molecules which form on the grain surface. We show that these forces can lead to grain dynamics in photodissociation regions which result in enhanced dust-to-gas ratios. Since the gas heating is probably dominated by photoelectric emission from dust, this might explain the unusually high gas temperatures inferred from ISO observations of molecular hydrogen in photodissociation regions.

Radiative torques on irregular dust grains, in addition to producing superthermal rotation, play a direct dynamical role in the alignment of interstellar dust with the local magnetic field. The equations governing the orientation of spinning, precessing grains are derived; H_2 formation torques and paramagnetic dissipation are included in the dynamics. Stationary solutions (constant alignment angle and spin rate) are found; these solutions may be stable ("attractors") or unstable ("repellors"). The equations of motion are numerically integrated for three exemplary irregular grain geometries, exposed to anisotropic radiation with the spectrum of interstellar starlight. The resulting "trajectory maps" are classified as "noncyclic", "semicyclic", or "cyclic", with examples of each given. We find that radiative torques result in rapid grain alignment, even in the absence of paramagnetic dissipation. It appears that radiative torques due to starlight can account for the observed alignment of interstellar grains with the Galactic magnetic field.

Irregular dust grains are subject to radiative torques when irradiated by interstellar starlight. It is shown how these radiative torques may be calculated using the discrete dipole approximation. Calculations are carried out for one irregular grain geometry, and three different grain sizes. It is shown that radiative torques can play an important dynamical role in spinup of interstellar dust grains, resulting in rotation rates which may exceed even those expected from H_2 formation on the grain surface. Because the radiative torque on an interstellar grain is determined by the overall grain geometry rather than merely the state of the grain surface, the resulting superthermal rotation is expected to be long-lived. By itself, long-lived superthermal rotation would permit grain alignment by normal paramagnetic dissipation on the "Davis-Greenstein" timescale. However, radiative torques arising from anisotropy of the starlight background can act directly to alter the grain alignment on much shorter timescales, and are therefore central to the process of interstellar grain alignment. Radiative torques depend strongly on the grain size, measured by a_eff, the radius of a sphere of equal volume. In diffuse clouds, radiative torques dominate the torques due to H2 formation for a_eff=0.2micron grains, but are relatively unimportant for a_eff<0.05micron grains. This may provide a natural explanation for the observation that $a_eff>0.1 micron grains in diffuse clouds are aligned, while there is little alignment of a_eff < 0.05 micron grains. We show that radiative torques are ineffective at producing superthermal rotation within quiescent dark clouds, but can be very effective in star-forming regions such as the M17 molecular cloud.

Recently (Brown \& Bethe 1994) it was suggested that most stars with main sequence mass in the range of about $18 - 30 M_{\odot}$ explode, returning matter to the Galaxy, and then go into low-mass ($\geq 1.5 M_{\odot}$) black holes. Even more massive main-sequence stars would, presumably, chiefly g o into high-mass ($\sim 10 M_{\odot}$) black holes. The Brown-Bethe estimates gave approximately $5 \times {10}^{8}$ low-mass black holes in the Galaxy. A pressing question, which we attempt to answer here, is why, with the possible exception of the compact objects in SN1987A and 4U\,1700--37, none of these have been seen. We address this question in three parts. Firstly, black holes are generally ``seen'' only in binaries, by the accretion of matter from a companion star. High mass black holes are capable of accreting more matter than low-mass black holes, so there is a selection effect favoring them. This, in itself, would not be sufficient to show why low-mass black holes have not been seen, since neutron stars (of nearly the same mass) are seen in abundance. Secondly, and this is our main point, the primary star in a binary ---the first star to evolve--- loses its hydrogen envelope by transfer of matter to the secondary and loss into space, and the resulting ``naked'' helium star evolves differently than a helium core, which is at least initially covered by the hydrogen envelope in a massive main-sequence star. We show that primary stars in binaries can end up as neutron stars even if their initial mass substantially exceeds the mass limit for neutron star formation from single stars ($\sim 18 M_{\odot}$). An example is 4U\,1223--62, in which we suggest that the initial primary mass exceeded $35 M_{\odot}$, yet X-ray pulsations