We present a combination of theoretical and simulation-based examinations of the role of two-fluid ambipolar drift on molecular line widths. The dissipation provided by ion-neutral interactions can produce a significant difference between the widths of neutral molecules and the widths of ionic species, comparable to the sound speed. We demonstrate that Alfven waves and certain families of magnetosonic waves become strongly damped on scales comparable to the ambipolar diffusion scale. Using the RIEMANN code, we simulate two-fluid turbulence with ionization fractions ranging from 10^{-2} to 10^{-6}. We show that the wave damping causes the power spectrum of the ion velocity to drop below that of the neutral velocity when measured on a relative basis. Following a set of motivational observations by Li & Houde (2008), we produce synthetic line width-size relations that shows a difference between the ion and neutral line widths, illustrating that two-fluid effects can have an observationally detectable role in modifying the MHD turbulence in the clouds.

We present a series of simulations of turbulent stratified protostellar discs with the goal of characterizing the settling of dust throughout a minimum-mass solar nebula. We compare the evolution of both compact spherical grains, as well as highly fractal grains. Our simulations use a shearing-box formulation to study the evolution of dust grains locally within the disc, and collectively our simulations span the entire extent of a typical accretion disc. The dust is stirred by gas that undergoes MRI-driven turbulence. This establishes a steady state scale height for the dust that is different for dust of different sizes. This sedimentation of dust is an important first step in planet formation and we predict that ALMA should be able to observationally verify its existence. When significant sedimentation occurs, the dust will participate in a streaming instability that significantly enhances the dust density. We show that the streaming instability is pervasive in the outer disc. We characterize the scale heights of dust whose size ranges from a few microns to a few centimeters. We find that for spherical grains, a power-law relationship develops for the scale height with grain size, with a slope that is slightly steeper than -1/2. The sedimentation is strongest in the outer disc and increases for large grains. The results presented here show that direct measurements of grain settling can be made by ALMA and we present favorable conditions for observability. The streaming instability should also be directly observable and we provide conditions for directly observing it. We calculate collision rates and growth rates for the dust grains in our simulations of various sizes colliding with other grains, and find that these rates are significantly enhanced through the density enhancement arising from the streaming instability.

(Abridged) In this paper we have used the RIEMANN code for computational astrophysics to study the interaction of a realistic distribution of dust grains with gas in a vertically stratified protostellar accretion disc. The disc was modeled to have the density and temperature of a minimum mass solar nebula, and was driven to a fully-developed turbulence via the magnetorotational instability (MRI). We find that the inclusion of standard dust to gas ratios does not have any significant effect on the MRI even when the dust sediments to the midplane of the accretion disc. The density distribution of the dust reaches a Gaussian profile, and the scale heights for the dust that we derive are shown to be proportional to the reciprocal of the square root of the dust radius. The largest dust shows a strong tendency to settle to the midplane of the accretion disc, and tends to organize itself into elongated clumps of high density. The dynamics of these clumps is shown to be consistent with a streaming instability. The streaming instability is seen to be very vigorous and persistent once it forms. Each stream of high density dust displays a reduced RMS velocity dispersion, and the densest clumpings of large dust are shown to form where the streams intersect. We have also shown that the mean free path and collision time for the dust that participates in the streaming instability is reduced by almost two orders of magnitude relative to the average mean free paths and collision times. We show that some of the large dust in our 10 au simulations should have a propensity for grain coalescence.

We present a semi-implicit method for isothermal two-fluid ion-neutral ambipolar drift that is second-order accurate in space and time. The method has been implemented in the RIEMANN code for astrophysical fluid dynamics. We present four test problems that show the method works and correctly tracks the propagation of MHD waves and the structure of two-fluid C-shocks. The accurate propagation of MHD waves in the two-fluid approximation is shown to be a stringent test of the algorithm. We demonstrate that highly accurate methods are required in order to properly capture the MHD wave behaviour in the presence of ion-neutral friction.

Anisotropic thermal conduction plays an important role in various astrophysical systems. One of the most stringent tests of thermal conduction can be found in supernova remnants. In this paper we study anisotropic thermal conduction and examine the physical nature of the flux of thermal conduction in the classical and saturated limits. We also present a temporally second-order accurate implicit-explicit scheme for the time-update of thermal conduction terms within a numerical MHD scheme. Several simulations of supernova remnants are presented for a range of ISM parameters. The role of thermal conduction in such remnants has been studied. We find that thermal conduction produces cooler temperatures and higher densities in the hot gas bubbles that form in the remnants. The effect of thermal conduction in changing the thermal characteristics of the hot gas bubble increases as the remnant propagates through denser ISMs. Remnants evolving in denser ISMs are shown to make a faster transition to a centre-bright x-ray morphology, with the trend emerging earlier in hard x-rays than in the soft x-rays.

We present a large number of two and a half dimensional simulations of supernova remnants expanding into interstellar media having a range of densities, temperatures and magnetic field strengths. The volume of hot gas produced is strongly dependent on the inclusion of thermal conduction and magnetic fields. The four-volumes and three-areas of hot gas have been catalogued and their dependence on interstellar parameters documented. Simulated line widths of radioactive species ejected by supernovae have also been catalogued.

(Abridged) We present a series of decaying turbulence simulations that represent a cluster-forming clump within a molecular cloud, investigating the role of magnetic fields on the formation of potential star-forming cores. We present an exhaustive analysis of numerical data from these simulations that includes a compilation of all of the distributions of physical properties that characterize bound cores - including their masses, radii, mean densities, angular momenta, spins, magnetizations, and mass-to-flux ratios. We also present line maps of our models that can be compared with observations. Our simulations range between 5-30 Jeans masses of gas, and are representative of molecular cloud clumps with masses between 100-1000 solar masses. The cores have mass-to-flux ratios that are generally less than that of the original cloud, and so a cloud that is initially highly supercritical can produce cores that are slightly supercritical, similar to that seen by Zeeman measurements of molecular cloud cores. Clouds that are initially only slightly supercritical will instead collapse along the field lines into sheets, and the cores that form as these sheets fragment have a different mass spectrum than what is observed. The spin rates of these cores suggests that subsequent fragmentation into multiple systems is likely. The sizes of the bound cores that are produced are typically 0.02-0.2 pc and have densities in the range 10^4-10^5 cm^{-3} in agreement with observational surveys. Finally, our numerical data allow us to test theoretical models of the mass spectrum of cores, such as the turbulent fragmentation picture of Padoan-Nordlund. We find that while this model gets the shape of the core mass spectrum reasonably well, it fails to predict the peak mass in the core mass spectrum.

We present new observations of the fundamental ro-vibrational CO spectrum of V1647 Ori, the young star whose recent outburst illuminated McNeil's Nebula. Previous spectra, acquired during outburst in 2004 February and July, had shown the CO emission lines to be broad and centrally peaked-similar to the CO spectrum of a typical classical T Tauri star. In this paper, we present CO spectra acquired shortly after the luminosity of the source returned to its pre-outburst level (2006 February) and roughly one year later (2006 December and 2007 February). The spectrum taken in 2006 February revealed blue-shifted CO absorption lines superimposed on the previously observed CO emission lines. The projected velocity, column density, and temperature of this outflowing gas was 30 km/s, 3^{+2}_{-1}E18 cm^{-2$, and 700^{+300}_{-100} K, respectively. The absorption lines were not observed in the 2006 December and 2007 February data, and so their strengths must have decreased in the interim by a factor of 9 or more. We discuss three mechanisms that could give rise to this unusual outflow.

We explore the role of anisotropic thermal conduction on the evolution of supernova remnants through interstellar media with a range of densities via numerical simulations. We find that a remnant expanding in a dense environment can produce centre-bright hard x-ray emission within 20 kyr, and centre-bright soft x-ray emission within 60 kyr of the supernova event. In a more tenuous environment, the appearance of a centre-bright structure in hard x-rays is delayed until about 60 kyr. The soft x-ray emission from such a remnant may not become centre bright during its observable lifetime. This can explain the observations that show that mixed-morphology supernova remnants preferentially occur close to denser, molecular environments. Remnants expanding into denser environments tend to be smaller, making it easier for thermal conduction to make larger changes in the temperatures of their hot gas bubbles. We show that the lower temperatures make it very favorable to use high-stage ions as diagnostics of the hot gas bubbles in SNRs. In particular, the distribution of O VIII transitions from shell-bright at early epochs to centre-bright at later epochs in the evolution of an SNR expanding in a dense ISM when the physics of thermal conduction is included.

We explore the importance of anisotropic thermal conduction in the evolution of supernova remnants via numerical simulations. The mean temperature of the bubble of hot gas is decreased by a factor of ~3 compared to simulations without thermal conduction, together with an increase in the mean density of hot gas by a similar factor. Thus, thermal conduction greatly reduces the volume of hot gas produced over the life of the remnant. This underscores the importance of thermal conduction in estimating the hot gas filling fraction and emissivities in high-stage ions in Galactic and proto-galactic ISMs.

We present high-resolution infrared spectra of four YSOs (T Tau N, T Tau S, RNO 91, and HL Tau). The spectra exhibit narrow absorption lines of 12CO, 13CO, and C18O as well as broad emission lines of gas phase12CO. The narrow absorption lines of CO are shown to originate from the colder circumstellar gas. We find that the line of sight gas column densities resulting from the CO absorption lines are much higher than expected for the measured extinction for each source and suggest the gas to dust ratio is measuring the dust settling and/or grain coagulation in these extended disks. We provide a model of turbulence, dust settling and grain growth to explain the results. The techniques presented here allow us to provide some observationally-motivated bounds on accretion disk alpha in protostellar systems.

We present the results of a series of numerical simulations of compressible, self-gravitating hydrodynamic turbulence of cluster-forming clumps in molecular clouds. We examine the role that turbulence has in the formation of gravitationally bound cores, studying the dynamical state, internal structure and bulk properties of these cores. Complex structure in turbulent clumps is formed provided that the damping time of the turbulence, $t_{damp}$ is longer than the gravitational free-fall time $t_{ff}$ in a region. We find a variety of density and infall velocity structures among the cores in the simulation, including cores that resemble the Larson-Penston collapse of an isothermal sphere ($\rho \propto r^{-2}$) as well as cores that resemble the McLaughlin-Pudritz collapse of logatropic spheres ($\rho \propto r^{-1}$). The specific angular momentum profiles range between $j \propto r^{1} - r^{2}$. The masses of the bound cores that form are well-fit by the turbulent mass spectrum of Padoan and Nordlund (2002), while the specific angular momentum distribution can be fit by a broken power law. While our hydrodynamic simulations reproduce many of the observed properties of cores, we find an upper limit for the star formation efficiency (SFE) in clusters of 40-50 per cent.

We develop models for the self-similar collapse of magnetized isothermal cylinders. We find solutions for the case of a fluid with a constant toroidal flux-to-mass ratio (Gamma_phi=constant) and the case of a fluid with a constant gas to magnetic pressure ratio (beta=constant). In both cases, we find that a low magnetization results in density profiles that behave as rho ~ r^{-4} at large radii, and at high magnetization we find density profiles that behave as rho ~ r^{-2}. This density behaviour is the same as for hydrostatic filamentary structures, suggesting that density measurements alone cannot distinguish between hydrostatic and collapsing filaments--velocity measurements are required. Our solutions show that the self-similar radial velocity behaves as v_r ~ r during the collapse phase, and that unlike collapsing self-similar spheres, there is no subsequent accretion (i.e. expansion-wave) phase. We also examine the fragmentation properties of these cylinders, and find that in both cases, the presence of a toroidal field acts to strengthen the cylinder against fragmentation. Finally, the collapse time scales in our models are shorter than the fragmentation time scales. Thus, we anticipate that highly collapsed filaments can form before they are broken into pieces by gravitational fragmentation.

We performed simulations of self-gravitating hydrodynamic turbulence to model the formation of filaments, clumps and cores in molecular clouds. We find that when the mass on the initial computational grid is comparable to the Jeans mass, turbulent pressure is able to prevent gravitational collapse. When the turbulence has damped away sufficiently, gravitational collapse can occur, and the resulting structure closely resembles the pre-singularity collapse of an isothermal sphere of Penston (1969). If several Jeans masses are initially placed on the grid, turbulence may not be sufficient to prevent collapse before turbulence can be significantly damped. In this case, the cores have density structures which are considerably shallower than expected for an isothermal gas, and resemble the solutions for a logatropic equation of state.

We study the evolution of gravitational waves through the preheating era that follows inflation. The oscillating inflaton drives parametric resonant growth of scalar field fluctuations, and although super-Hubble tensor modes are not strongly amplified, they do carry an imprint of preheating. This is clearly seen in the Weyl tensor, which provides a covariant description of gravitational waves.