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

Non-ideal MHD turbulent decay in molecular clouds

It is well known that non-ideal magnetohydrodynamic effects are important in the dynamics of molecular clouds: both ambipolar diffusion and possibly the Hall effect have been identified as significant. We present the results of a suite of simulations with a resolution of 512-cubed of turbulent decay in molecular clouds incorporating a simplified form of both ambipolar diffusion and the Hall effect simultaneously. The initial velocity field in the turbulence is varied from being super-Alfv\'enic and hypersonic, through to trans-Alfv\'enic but still supersonic. We find that ambipolar diffusion increases the rate of decay of the turbulence increasing the decay from $t^{-1.25}$ to $t^{-1.4}$. The Hall effect has virtually no impact in this regard. The power spectra of density, velocity and the magnetic field are all affected by the non-ideal terms, being steepened significantly when compared with ideal MHD turbulence with exponents. The density power spectra components change from about 1.4 to about 2.1 for the ideal and non-ideal simulations respectively, and power spectra of the other variables all show similar modifications when non-ideal effects are considered. Again, the dominant source of these changes is ambipolar diffusion rather than the Hall effect. There is also a decoupling between the velocity field and the magnetic field at short length scales. The Hall effect leads to enhanced magnetic reconnection, and hence less power, at short length scales. The dependence of the velocity dispersion on the characteristic length scale is studied and found not to be power-law in nature.Comment: 16 pages, 17 figure

Massive molecular outflows at high spatial resolution

We present high-spatial resolution Plateau de Bure Interferometer CO(2-1) and SiO(2-1) observations of one intermediate-mass and one high-mass star-forming region. The intermediate-mass region IRAS20293+3952 exhibits four molecular outflows, one being as collimated as the highly collimated jet-like outflows observed in low-mass star formation sources. Furthermore, comparing the data with additional infrared H2 and cm observations we see indications that the nearby ultracompact HII region triggers a shock wave interacting with the outflow. The high-mass region IRAS19217+1651 exhibits a bipolar outflow as well and the region is dominated by the central driving source. Adding two more sources from the literature, we compare position-velocity diagrams of the intermediate- to high-mass sources with previous studies in the low-mass regime. We find similar kinematic signatures, some sources can be explained by jet-driven outflows whereas other are better constrained by wind-driven models. The data also allow to estimate accretion rates varying from a few times 10^{-5}Msun/yr for the intermediate-mass sources to a few times 10^{-4}Msun/yr for the high-mass source, consistent with models explaining star formation of all masses via accretion processes.Comment: 14 pages text, 4 tables, 8 figures, accepted for Ap

Relativistic Blastwaves and Synchrotron Emission

Relativistic shocks accelerate particles by the first order Fermi mechanism. These particles then emit synchrotron emission in the post shock gas. We have developed a numerical code which integrates the relativistic Euler equations for fluid dynamics with a general equation of state, together with the Liouville equation for the accelerated particles. We present tests of this code and, in addition, we use it to study the gamma ray burst afterglow predicted by the fireball model, along with the hydrodynamics of a relativistic blastwave. We find that, while, broadly speaking, the behaviour of the emission is similar to that already predicted with semi-analytic approaches, the detailed behaviour is somewhat different. The ``breaks'' in the synchrotron spectrum behave differently with time, and the spectrum above the final break is harder than previously expected. These effects are due to the incorporation of the geometry of the (spherical) blastwave, along with relativistic beaming and adiabatic cooling of the energetic particles leading to a mix, in the observed spectrum, between recently injected "uncooled" particles and the older "cooled" population in different parts of the evolving, inhomogeneous flow.Comment: 12 pages, 10 figures, accepted for publication in MNRAS. Expanded discussion in section 5, more tests of the code, and other minor change

Dense Molecular Gas and the Role of Star Formation in the Host Galaxies of Quasi-Stellar Objects

New millimeter-wave CO and HCN observations of the host galaxies of infrared-excess Palomar Green quasi-stellar objects (PG QSOs) previously detected in CO are presented. These observations are designed to assess the validity of using the infrared luminosity to estimate star formation rates of luminous AGN by determining the relative significance of dust-heating by young, massive stars and active galactic nuclei (AGN) in QSO hosts and IRAS galaxies with warm, AGN-like infrared colors. The HCN data show the PG QSO host IZw1 and most of the warm IRAS galaxies to have high L_IR / L'_HCN (>1600) relative to the cool IRAS galaxy population for which the median L_IR / L'_HCN ~ 890(+440,-470). If the assumption is made that the infrared emission from cool IRAS galaxies is reprocessed light from embedded star-forming regions, then high values of L_IR / L'_HCN are likely the result of dust heating by the AGN. Further, if the median ratio of L'_HCN / L'_CO ~ 0.06 observed for Seyfert galaxies and IZw1 is applied to the PG QSOs not detected in HCN, then the derived L_IR / L'_HCN correspond to a stellar contribution to the production of L_IR of ~ 7-39%, and star formation rates ~ 2-37 M_sun/yr are derived for the QSO hosts. Alternatively, if the far-infrared is adopted as the star formation component of the total infrared in cool galaxies, the stellar contributions in QSO hosts to their L_FIR are up to 35% higher than the percentages derived for L_IR. This raises the possibility that the L_FIR in several of the PG QSO hosts, including IZw1, could be due entirely to dust heated by young, massive stars. Finally, there is no evidence that the global HCN emission is enhanced relative to CO in galaxies hosting luminous AGN.Comment: LaTex, 31 pages, including 9 postscript figures, AJ, in press (December 2006

Spatial variation of the physical conditions of molecular gas in galaxies

Multi-line studies of CO-12, CO-13, C-18O, HCN, and HCO(+) at 3 mm, 1.3 mm, and 0.8 mm using the Institute for Radio Astronomy in the Millimeter range (IRAM) 30 m telescope, with the IRAM superconductor insulator superconductor (SIS) receivers and the Max Planck Institute for External Physics (MPE) 350 GHz SIS receiver, show that the densities and temperatures of molecular gas in external galaxies change significantly with position. CO-12 measures the densities and temperature of diffuse interclump molecular gas, but not the bulk of the molecular gas. Simple one-component models, with or without external heating, cannot account for the weakness of the CO-12 J = 3 to 2 line relative to J = 2 to 1 and J = 1 to 0. CO-12 does not trace the bulk of the molecular gas, and optical depth effects obviate a straightforward interpretation of CO-12 data. Instead, researchers turned to the optically thin CO isotopes and other molecular species. Isotopic CO lines measure the bulk of the molecular gas, and HCN and HCO(+) pick out denser regions. Researchers find a warm ridge of gas in IC 342 (Eckart et al. 1989), denser gas in the starburst nucleus of IC 342, and a possible hot-spot in NGC 2903. In IC 342, NGC 2146, and NGC 6764, the CO-13 J = 2 to 1 line is subthermally populated, implying gas densities less than or equal to 10(exp 4) cm(-3)

MHDSTS: a new explicit numerical scheme for simulations of partially ionised solar plasma

The interaction of plasma with magnetic field in the partially ionised solar atmosphere is frequently modelled via a single-fluid approximation, which is valid for the case of a strongly coupled collisional media, such as solar photosphere and low chromosphere. Under the single-fluid formalism the main non-ideal effects are described by a series of extra terms in the generalised induction equation and in the energy conservation equation. These effects are: Ohmic diffusion, ambipolar diffusion, the Hall effect, and the Biermann battery effect. From the point of view of the numerical solution of the single-fluid equations, when ambipolar diffusion or Hall effects dominate can introduce severe restrictions on the integration time step and can compromise the stability of the numerical scheme. In this paper we introduce two numerical schemes to overcome those limitations. The first of them is known as Super Time-Stepping (STS) and it is designed to overcome the limitations imposed when the ambipolar diffusion term is dominant. The second scheme is called the Hall Diffusion Scheme (HDS) and it is used when the Hall term becomes dominant. These two numerical techniques can be used together by applying Strang operator splitting. This paper describes the implementation of the STS and HDS schemes in the single-fluid code Mancha3D. The validation for each of these schemes is provided by comparing the analytical solution with the numerical one for a suite of numerical tests