The mass-velocity and intensity-velocity relations in jet-driven molecular outflows

We use numerical simulations to examine the mass-velocity and intensity-velocity relations in the CO J=2-1 and H$_2$ S(1)1-0 lines for jet-driven molecular outflows. Contrary to previous expectations, we find that the mass-velocity relation for the swept-up gas is a single power-law, with a shallow slope $\simeq -1.5$ and no break to a steeper slope at high velocities. An analytic bowshock model with no post-shock mixing is shown to reproduce this behaviour very well. We show that molecular dissociation and the temperature dependence of the line emissivity are both critical in defining the shape of the line profiles at velocities above $\sim$ 20 km s$^{-1}$. In particular, the simulated CO J=2-1 intensity-velocity relation does show a break in slope, even though the underlying mass distribution does not. These predicted CO profiles are found to compare remarkably well with observations of molecular outflows, both in terms of the slopes at low and high velocities and in terms of the range of break velocities at which the change in slope occurs. Shallower slopes are predicted at high velocity in higher excitation lines, such as H$_2$ S(1)1-0. This work indicates that, in jet-driven outflows, the CO J=2-1 intensity profile reflects the slope of the underlying mass-velocity distribution only at velocities $\le$ 20 km/s, and that higher temperature tracers are required to probe the mass distribution at higher speed.Comment: 6 pages, 8 figures. Accepted for publication in Astronomy and Astrophysic

An explicit scheme for multifluid magnetohydrodynamics

When modeling astrophysical fluid flows, it is often appropriate to discard the canonical magnetohydrodynamic approximation thereby freeing the magnetic field to diffuse with respect to the bulk velocity field. As a consequence, however, the induction equation can become problematic to solve via standard explicit techniques. In particular, the Hall diffusion term admits fast-moving whistler waves which can impose a vanishing timestep limit. Within an explicit differencing framework, a multifluid scheme for weakly ionised plasmas is presented which relies upon a new approach to integrating the induction equation efficiently. The first component of this approach is a relatively unknown method of accelerating the integration of parabolic systems by enforcing stability over large compound timesteps rather than over each of the constituent substeps. This method, Super Time Stepping, proves to be very effective in applying a part of the Hall term up to a known critical value. The excess of the Hall term above this critical value is then included via a new scheme for pure Hall diffusion.Comment: 8 pages; 4 figures; accepted by MNRAS; minor corrections to equations; addition of appendi

Multifluid magnetohydrodynamic turbulent decay

It is generally believed that turbulence has a significant impact on the dynamics and evolution of molecular clouds and the star formation which occurs within them. Non-ideal magnetohydrodynamic effects are known to influence the nature of this turbulence. We present the results of a suite of 512-cubed resolution simulations of the decay of initially super-Alfvenic and supersonic fully multifluid MHD turbulence. We find that ambipolar diffusion increases the rate of decay of the turbulence while the Hall effect has virtually no impact. The decay of the kinetic energy can be fitted as a power-law in time and the exponent is found to be -1.34 for fully multifluid MHD turbulence. The power spectra of density, velocity and magnetic field are all steepened significantly by the inclusion of non-ideal terms. The dominant reason for this steepening is ambipolar diffusion with the Hall effect again playing a minimal role except at short length scales where it creates extra structure in the magnetic field. Interestingly we find that, at least at these resolutions, the majority of the physics of multifluid turbulence can be captured by simply introducing fixed (in time and space) resistive terms into the induction equation without the need for a full multifluid MHD treatment. The velocity dispersion is also examined and, in common with previously published results, it is found not to be power-law in nature.Comment: 16 pages, 15 figures, Accepted for publication in Ap

A New Time-Dependent Finite Difference Method for Relativistic Shock Acceleration

We present a new approach to calculate the particle distribution function about relativistic shocks including synchrotron losses using the method of lines with an explicit finite difference scheme. A steady, continuous, one dimensional plasma flow is considered to model thick (modified) shocks, leading to a calculation in three dimensions plus time, the former three being momentum, pitch angle and position. The method accurately reproduces the expected power law behaviour in momentum at the shock for upstream flow speeds ranging from 0.1c to 0.995c (1 < \Gamma < 10). It also reproduces approximate analytical results for the synchrotron cutoff shape for a non-relativistic shock, demonstrating that the loss process is accurately represented. The algorithm has been implemented as a hybrid OpenMP--MPI parallel algorithm to make efficient use of SMP cluster architectures and scales well up to many hundreds of CPUs.Comment: Accepted for publication in MNRA

Asteroidal differentiation processes deduced from ultramafic achondrite ureilite meteorites

Ureilite meteorites are partial melt residues of an asteroid-sized object. They record the differentiation process that transformed many asteroids during the earliest stages of solar system formation