Atomic Diagnostics of X-ray Irradiated Protoplanetary Disks

We study atomic line diagnostics of the inner regions of protoplanetary disks with our model of X-ray irradiated disk atmospheres which was previously used to predict observable levels of the NeII and NeIII fine-structure transitions at 12.81 and 15.55mum. We extend the X-ray ionization theory to sulfur and calculate the fraction of sulfur in S, S+, S2+ and sulfur molecules. For the D'Alessio generic T Tauri star disk, we find that the SI fine-structure line at 25.55mum is below the detection level of the Spitzer Infrared Spectrometer (IRS), in large part due to X-ray ionization of atomic S at the top of the atmosphere and to its incorporation into molecules close to the mid-plane. We predict that observable fluxes of the SII 6718/6732AA forbidden transitions are produced in the upper atmosphere at somewhat shallower depths and smaller radii than the neon fine-structure lines. This and other forbidden line transitions, such as the OI 6300/6363AA and the CI 9826/9852AA lines, serve as complementary diagnostics of X-ray irradiated disk atmospheres. We have also analyzed the potential role of the low-excitation fine-structure lines of CI, CII, and OI, which should be observable by SOFIA and Herschel.Comment: Accepted by Ap

A Spherical Model for "Starless" Cores of Magnetic Molecular Clouds and Dynamical Effects of Dust Grains

In the standard picture of isolated star formation, dense ``starless'' cores are formed out of magnetic molecular clouds due to ambipolar diffusion. Under the simplest spherical geometry, I demonstrate that ``starless'' cores formed this way naturally exhibit a large scale inward motion, whose size and speed are comparable to those detected recently by Taffala et al. and Williams et al. in ``starless'' core L1544. My model clouds have a relatively low mass (of order 10 $M_\odot$) and low field strength (of order 10 $\mu$G) to begin with. They evolve into a density profile with a central plateau surrounded by a power-law envelope, as found previously. The density in the envelope decreases with radius more steeply than those found by Mouschovias and collaborators for the more strongly magnetized, disk-like clouds. At high enough densities, dust grains become dynamically important by greatly enhancing the coupling between magnetic field and the neutral cloud matter. The trapping of magnetic flux associated with the enhanced coupling leads, in the spherical geometry, to a rapid assemblage of mass by the central protostar, which exacerbates the so-called ``luminosity problem'' in star formation.Comment: 27 pages, 4 figures, accepted by Ap

TMC-1C: an accreting starless core

We have mapped the starless core TMC-1C in a variety of molecular lines with the IRAM 30m telescope. High density tracers show clear signs of self-absorption and sub-sonic infall asymmetries are present in N2H+ (1-0) and DCO+ (2-1) lines. The inward velocity profile in N2H+ (1-0) is extended over a region of about 7,000 AU in radius around the dust continuum peak, which is the most extended ``infalling'' region observed in a starless core with this tracer. The kinetic temperature (~12 K) measured from C17O and C18O suggests that their emission comes from a shell outside the colder interior traced by the mm continuum dust. The C18O (2-1) excitation temperature drops from 12 K to ~10 K away from the center. This is consistent with a volume density drop of the gas traced by the C18O lines, from ~4x10^4 cm^-3 towards the dust peak to ~6x10^3 cm^-3 at a projected distance from the dust peak of 80" (or 11,000 AU). The column density implied by the gas and dust show similar N2H+ and CO depletion factors (f_D < 6). This can be explained with a simple scenario in which: (i) the TMC-1C core is embedded in a relatively dense environment (H2 ~10^4 cm^-3), where CO is mostly in the gas phase and the N2H+ abundance had time to reach equilibrium values; (ii) the surrounding material (rich in CO and N2H+) is accreting onto the dense core nucleus; (iii) TMC-1C is older than 3x10^5 yr, to account for the observed abundance of N2H+ across the core (~10^-10 w.r.t. H2); and (iv) the core nucleus is either much younger (~10^4 yr) or ``undepleted'' material from the surrounding envelope has fallen towards it in the past 10,000 yr.Comment: 29 pages, including 5 tables and 15 figure

On the Rapid Collapse and Evolution of Molecular Clouds

Stars generally form faster than the ambipolar diffusion time, suggesting that several processes short circuit the delay and promote a rapid collapse. These processes are considered here, including turbulence compression in the outer parts of giant molecular cloud (GMC) cores and GMC envelopes, GMC core formation in an initially supercritical state, and compression-induced triggering in dispersing GMC envelopes. The classical issues related to star formation timescales are addressed: high molecular fractions, low efficiencies, long consumption times for CO and HCN, rapid GMC core disruption and the lack of a stable core, long absolute but short relative timescales with accelerated star formation, and the slow motions of protostars. We consider stimuli to collapse from changes in the density dependence of the ionization fraction, the cosmic ray ionization rate, and various dust properties at densities above ~10^5 cm^{-3}. We favor the standard model of subcritical GMC envelops and suggest they would be long lived if not for disruption by rapid star formation in GMC cores. The lifecycle of GMCs is illustrated by a spiral arm section in the Hubble Heritage image of M51, showing GMC formation, star formation, GMC disruption with lingering triggered star formation, and envelope dispersal. There is no delay between spiral arm dustlanes and star formation; the classical notion results from heavy extinction in the dust lane and triggered star formation during cloud dispersal. Differences in the IMF for the different modes of star formation are considered.Comment: 46 pages, 5 figures, scheduled for ApJ 668, October 20, 200

Molecular ions in L1544. II. The ionization degree

The maps presented in Paper I are here used to infer the variation of the column densities of HCO+, DCO+, N2H+, and N2D+ as a function of distance from the dust peak. These results are interpreted with the aid of a crude chemical model which predicts the abundances of these species as a function of radius in a spherically symmetric model with radial density distribution inferred from the observations of dust emission at millimeter wavelengths and dust absorption in the infrared. Our main observational finding is that the N(N2D+)/N(N2H+) column density ratio is of order 0.2 towards the L1544 dust peak as compared to N(DCO+)/N(HCO+) = 0.04. We conclude that this result as well as the general finding that N2H+ and N2D+ correlate well with the dust is caused by CO being depleted to a much higher degree than molecular nitrogen in the high density core of L1544. Depletion also favors deuterium enhancement and thus N2D+, which traces the dense and highly CO-depleted core nucleus, is much more enhanced than DCO+. Our models do not uniquely define the chemistry in the high density depleted nucleus of L1544 but they do suggest that the ionization degree is a few times 10^{-9} and that the ambipolar diffusion time scale is locally similar to the free fall time. It seems likely that the lower limit which one obtains to ionization degree by summing all observable molecular ions is not a great underestimate of the true ionization degree. We predict that atomic oxygen is abundant in the dense core and, if so, H3O+ may be the main ion in the central highly depleted region of the core.Comment: 31 pages, 8 figures, to be published in Ap

The role of damped Alfven waves on magnetospheric accretion models of young stars

We examine the role of Alfven wave damping in heating the plasma in the magnetic funnels of magnetospheric accretion models of young stars. We study four different damping mechanisms of the Alfven waves: nonlinear, turbulent, viscous-resistive and collisional. Two different possible origins for the Alfven waves are discussed: 1) Alfven waves generated at the surface of the star by the shock produced by the infalling matter; and 2) Alfven waves generated locally in the funnel by the Kelvin-Helmholtz instability. We find that, in general, the damping lengths are smaller than the tube length. Since thermal conduction in the tube is not efficient, Alfven waves generated only at the star's surface cannot heat the tube to the temperatures necessary to fit the observations. Only for very low frequency Alfven waves ~10^{-5} the ion cyclotron frequency, is the viscous-resistive damping length greater than the tube length. In this case, the Alfven waves produced at the surface of the star are able to heat the whole tube. Otherwise, local production of Alfven waves is required to explain the observations. The turbulence level is calculated for different frequencies for optically thin and thick media. We find that turbulent velocities varies greatly for different damping mechanisms, reaching \~100 km s^{-1} for the collisional damping of small frequency waves.Comment: 29 pages, 12 figures, to appear in The Astrophysical Journa

Cosmic-ray ionization of molecular clouds

Low-energy cosmic rays are a fundamental source of ionization for molecular clouds, influencing their chemical, thermal and dynamical evolution. The purpose of this work is to explore the possibility that a low-energy component of cosmic-rays, not directly measurable from the Earth, can account for the discrepancy between the ionization rate measured in diffuse and dense interstellar clouds. We collect the most recent experimental and theoretical data on the cross sections for the production of H2+ and He+ by electron and proton impact, and we discuss the available constraints on the cosmic-ray fluxes in the local interstellar medium. Starting from different extrapolations at low energies of the demodulated cosmic-ray proton and electron spectra, we compute the propagated spectra in molecular clouds in the continuous slowing-down approximation taking into account all the relevant energy loss processes. The theoretical value of the cosmic-ray ionization rate as a function of the column density of traversed matter is in agreement with the observational data only if either the flux of cosmic-ray electrons or of protons increases at low energies. The most successful models are characterized by a significant (or even dominant) contribution of the electron component to the ionization rate, in agreement with previous suggestions. However, the large spread of cosmic-ray ionization rates inferred from chemical models of molecular cloud cores remains to be explained. Available data combined with simple propagation models support the existence of a low-energy component (below about 100 MeV) of cosmic-ray electrons or protons responsible for the ionization of molecular cloud cores and dense protostellar envelopes.Comment: 14 pages, 15 figure

Impact of grain evolution on the chemical structure of protoplanetary disks

We study the impact of dust evolution in a protoplanetary disk around a T Tauri star on the disk chemical composition. For the first time we utilize a comprehensive model of dust evolution which includes growth, fragmentation and sedimentation. Specific attention is paid to the influence of grain evolution on the penetration of the UV field in the disk. A chemical model that includes a comprehensive set of gas phase and grain surface chemical reactions is used to simulate the chemical structure of the disk. The main effect of the grain evolution on the disk chemical composition comes from sedimentation, and, to a lesser degree, from the reduction of the total grain surface area. The net effect of grain growth is suppressed by the fragmentation process which maintains a population of small grains, dominating the total grain surface area. We consider three models of dust properties. In model GS both growth and sedimentation are taken into account. In models A5 and A4 all grains are assumed to have the same size (10(-5) cm and 10(-4) cm, respectively) with constant gas-to-dust mass ratio of 100. Like in previous studies, the "three-layer" pattern (midplane, molecular layer, hot atmosphere) in the disk chemical structure is preserved in all models, but shifted closer to the midplane in models with increased grain size (GS and A4). Unlike other similar studies, we find that in models GS and A4 column densities of most gas-phase species are enhanced by 1-3 orders of magnitude relative to those in a model with pristine dust (A5), while column densities of their surface counterparts are decreased. We show that column densities of certain species, like C2H, HC(2n+1)N (n=0-3), H2O and some other molecules, as well as the C2H2/HCN abundance ratio which are accessible with Herschel and ALMA can be used as observational tracers of early stages of the grain evolution process in protoplanetary disks.Comment: 50 pages, 4 tables, 11 figures, accepted to the Ap

Molecular Hydrogen Emission from Protoplanetary Disks II. Effects of X-ray Irradiation and Dust Evolution

Detailed models for the density and temperature profiles of gas and dust in protoplanetary disks are constructed by taking into account X-ray and ultraviolet (UV) irradiation from a central T Tauri star, as well as dust size growth and settling toward the disk midplane. The spatial and size distributions of dust grains in the disks are numerically computed by solving the coagulation equation for settling dust particles. The level populations and line emission of molecular hydrogen are calculated using the derived physical structure of the disks. X-ray irradiation is the dominant heating source of the gas in the inner disk region and in the surface layer, while the far UV heating dominates otherwise. If the central star has strong X-ray and weak UV radiation, the H2 level populations are controlled by X-ray pumping, and the X-ray induced transition lines could be observable. If the UV irradiation is strong, the level populations are controlled by thermal collisions or UV pumping, depending on the properties of the dust grains in the disks. As the dust particles evolve in the disks, the gas temperature at the disk surface drops because the grain photoelectric heating becomes less efficient, while the UV radiation fields become stronger due to the decrease of grain opacity. This makes the H2 level populations change from local thermodynamic equilibrium (LTE) to non-LTE distributions, which results in changes to the line ratios of H2 emission. Our results suggest that dust evolution in protoplanetary disks could be observable through the H2 line ratios. The emission lines are strong from disks irradiated by strong UV and X-rays and possessing small dust grains; such disks will be good targets in which to observe H2 emission.Comment: 33 pages, accepted for publication in the Astrophysical Journa

Magnetic fields in protoplanetary disks

Magnetic fields likely play a key role in the dynamics and evolution of protoplanetary discs. They have the potential to efficiently transport angular momentum by MHD turbulence or via the magnetocentrifugal acceleration of outflows from the disk surface, and magnetically-driven mixing has implications for disk chemistry and evolution of the grain population. However, the weak ionisation of protoplanetary discs means that magnetic fields may not be able to effectively couple to the matter. I 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. For a standard population of 0.1 micron grains the active surface layers have a combined column of about 2 g/cm^2 at 1 AU; by the time grains have aggregated to 3 microns the active surface density is 80 g/cm^2. In the absence of grains, x-rays maintain magnetic coupling to 10% of the disk material at 1 AU (150 g/cm^2). At 5 AU the entire disk thickness becomes active once grains have aggregated to 1 micron in size.Comment: 11 pages, 11 figs, aastex.cls. Accepted for publication in Astrophysics & Space Science. v3 corrects bibliograph