Pico: Parameters for the Impatient Cosmologist

We present a fast, accurate, robust and flexible method of accelerating parameter estimation. This algorithm, called Pico, can compute the CMB power spectrum and matter transfer function as well as any computationally expensive likelihoods in a few milliseconds. By removing these bottlenecks from parameter estimation codes, Pico decreases their computational time by 1 or 2 orders of magnitude. Pico has several important properties. First, it is extremely fast and accurate over a large volume of parameter space. Furthermore, its accuracy can continue to be improved by using a larger training set. This method is generalizable to an arbitrary number of cosmological parameters and to any range of l-values in multipole space. Pico is approximately 3000 times faster than CAMB for flat models, and approximately 2000 times faster then the WMAP 3 year likelihood code. In this paper, we demonstrate that using Pico to compute power spectra and likelihoods produces parameter posteriors that are very similar to those using CAMB and the official WMAP3 code, but in only a fraction of the time. Pico and an interface to CosmoMC are made publicly available at http://www.astro.uiuc.edu/~bwandelt/pico/.Comment: 9 pages, 10 figures, submitted to ApJ, LaTeX with emulateap

Collimation of astrophysical jets - the role of the accretion disk magnetic field distribution

We have applied axisymmetric MHD simulations to investigate the impact of the accretion disk magnetic flux profile on the jet collimation. Using the ZEUS-3D code modified for magnetic diffusivity, our simulations evolve from an initial hydrostatic equilibrium state in a force-free magnetic field. Considering a power law for the disk poloidal magnetic field profile Bp ~ r^{-mu} and for the disk wind density profile rho ~ r^{-mu_rho} we performed a systematic study over a wide parameter range mu and mu_rho. We find a degree of collimation (ratio of mass flow rates in axial and lateral direction) decreasing for steeper disk magnetic field profiles (increasing mu). Varying the total magnetic flux doesn't change the degree of jet collimation substantially, it only affects the time scale of outflow evolution and the terminal jet speed. As our major result we find a general relation between the collimation degree with the disk wind magnetization power law exponent. Outflows with high collimation degree resulting from a flat disk magnetic field profile tend to be unsteady, producing axially propagating knots as discussed earlier. Depending slightly on the inflow density profile this unsteady behavior sets in for mu < 0.4. We also performed simulations of jet formation with artificially enhanced decay of the toroidal magnetic field in order to investigate the idea of a purely "poloidal collimation" discussed in the literature. These outflows remain weakly collimated and propagate with lower velocity. Thanks to our large numerical grid size (7x14 AU for protostars), we may apply our results to recently observed hints of jet rotation (DG Tau) indicating a relatively flat disk magnetic field profile, mu ~ 0.5. In general, our results are applicable to both stellar and extragalactic sources of MHD jets.Comment: accepted by ApJ, high resolution version under www.mpia-hd.mpg.de/homes/fendt

Relativistic expansion of a magnetized fluid

We study semi-analytical time-dependent solutions of the relativistic magnetohydrodynamic (MHD) equations for the fields and the fluid emerging from a spherical source. We assume uniform expansion of the field and the fluid and a polytropic relation between the density and the pressure of the fluid. The expansion velocity is small near the base but approaches the speed of light at the light sphere where the flux terminates. We find self-consistent solutions for the density and the magnetic flux. The details of the solution depend on the ratio of the toroidal and the poloidal magnetic field, the ratio of the energy carried by the fluid and the electromagnetic field and the maximum velocity it reaches.Comment: 17 pages, 6 figures, accepted by Geophysical and Astrophysical Fluid Dynamic

Ultra-Relativistic Magneto-Hydro-Dynamic Jets in the context of Gamma Ray Bursts

We present a detailed numerical study of the dynamics and evolution of ultrarelativistic magnetohydrodynamic jets in the black hole-disk system under extreme magnetization conditions. We find that Lorentz factors of up to 3000 are achieved and derived a modifiedMichel scaling (Gamma ~ sigma) which allows for a wide variation in the flow Lorentz factor. Pending contamination induced by mass-entrainment, the linear Michel scaling links modulations in the ultrarelativistic wind to variations in mass accretion in the disk for a given magnetization. The jet is asymptotically dominated by the toroidal magnetic field allowing for efficient collimation. We discuss our solutions (jets) in the context of Gamma ray bursts and describe the relevant features such as the high variability in the Lorentz factor and how high collimation angles (~ 0-5 degrees), or cylindrical jets, can be achieved. We isolate a jet instability mechanism we refer to as the "bottle-neck" instability which essentially relies on a high magnetization and a recollimation of the magnetic flux surfaces. The instability occurs at large radii where any dissipation of the magnetic energy into radiation would in principle result in an optically thin emission.Comment: 31 pages, 6 figures. Submitted to ApJ. Higher Quality figures at http://www.capca.ucalgary.ca/paper

Rico: An Accurate Cosmological Recombination Code

We present Rico, a code designed to compute the ionization fraction of the Universe during the epoch of hydrogen and helium recombination with an unprecedented combination of speed and accuracy. This is accomplished by training the machine learning code Pico on the calculations of a multi-level cosmological recombination code which self-consistently includes several physical processes that were neglected previously. After training, Rico is used to fit the free electron fraction as a function of the cosmological parameters. While, for example at low redshifts (z<~900), much of the net change in the ionization fraction can be captured by lowering the hydrogen fudge factor in Recfast by about 3%, Rico provides a means of effectively using the accurate ionization history of the full recombination code in the standard cosmological parameter estimation framework without the need to add new or refined fudge factors or functions to a simple recombination model. Within the new approach presented here it is easy to update Rico whenever a more accurate full recombination code becomes available. Once trained, Rico computes the cosmological ionization history with negligible fitting error in ~10 milliseconds, a speed-up of at least 10^6 over the full recombination code that was used here. Also Rico is able to reproduce the ionization history of the full code to a level well below 0.1%, thereby ensuring that the theoretical power spectra of CMB fluctuations can be computed to sufficient accuracy and speed for analysis from upcoming CMB experiments like Planck. Furthermore it will enable cross-checking different recombination codes across cosmological parameter space, a comparison that will be very important in order to assure the accurate interpretation of future cosmic microwave background data.Comment: 14 pages, 11 figures, submitted to PR

VLBI imaging of a flare in the Crab Nebula: More than just a spot

We report on very long baseline interferometry (VLBI) observations of the radio emission from the inner region of the Crab Nebula, made at 1.6 GHz and 5 GHz after a recent high-energy flare in this object. The 5 GHz data have provided only upper limits of 0.4 milli-Jansky (mJy) on the flux density of the pulsar and 0.4 mJy/beam on the brightness of the putative flaring region. The 1.6 GHz data have enabled imaging the inner regions of the nebula on scales of up to ~40". The emission from the inner "wisps" is detected for the first time with VLBI observations. A likely radio counterpart (designated "C1") of the putative flaring region observed with Chandra and HST is detected in the radio image, with an estimated flux density of $0.5\pm 0.3$\,mJy and a size of 0.2-0.6". Another compact feature ("C2") is also detected in the VLBI image closer to the pulsar, with an estimated flux density of 0.4 +- 0.2 mJy and a size smaller than 0{\farcs}2. Combined with the broad-band SED of the flare, the radio properties of C1 yield a lower limit of ~0.5 mG for the magnetic field and a total minimum energy of 1.2*10^41 ergs vested in the flare (corresponding to using about 0.2% of the pulsar spin-down power). The 1.6 GHz observations provide upper limits for the brightness (0.2 mJy/beam) and total flux density (0.4 mJy) of the optical Knot 1 located at 0.6" from the pulsar. The absolute position of the Crab pulsar is determined, and an estimate of the pulsar proper motion is obtained.Comment: Astronomy & Astrophysics; accepted; 10 pages, 8 figure

Molecular outflows in the young open cluster IC348

We present a wide-field survey of the young open cluster IC348 for molecular H2 outflows. Outflow activity is only found at its south-western limit, where a new subcluster of embedded sources is in an early phase of its formation. If the IC348 cluster had been built up by such subclusters forming at different times, this could explain the large age-spread that Herbig (1998) found for the IC348 member stars. In addition to several compact groups of H2 knots, our survey reveals a large north-south oriented outflow, and we identify the newly discovered far-infrared and mm-object IC348MMS as its source. New deep images in the 1-0 S(1) line of molecular hydrogen trace the HH211 jet and counterjet as highly-collimated chains of knots, resembling the interferometric CO and SiO jets. This jet system appears rotated counter-clockwise by about 3 degrees with respect to the prominent H2 bow shocks. Furthermore, we resolve HH211-mm as a double point-like source in the mm-continuum.Comment: 10 pages, 9 figures, accepted for publication in Ap

Can Protostellar Jets Drive Supersonic Turbulence in Molecular Clouds?

Jets and outflows from young stellar objects are proposed candidates to drive supersonic turbulence in molecular clouds. Here, we present the results from multi-dimensional jet simulations where we investigate in detail the energy and momentum deposition from jets into their surrounding environment and quantify the character of the excited turbulence with velocity probability density functions. Our study include jet--clump interaction, transient jets, and magnetised jets. We find that collimated supersonic jets do not excite supersonic motions far from the vicinity of the jet. Supersonic fluctuations are damped quickly and do not spread into the parent cloud. Instead subsonic, non-compressional modes occupy most of the excited volume. This is a generic feature which can not be fully circumvented by overdense jets or magnetic fields. Nevertheless, jets are able to leave strong imprints in their cloud structure and can disrupt dense clumps. Our results question the ability of collimated jets to sustain supersonic turbulence in molecular clouds.Comment: 33 pages, 18 figures, accepted by ApJ, version with high resolution figures at: http://www.ita.uni-heidelberg.de/~banerjee/publications/jet_paper.pd

The Propagation of Magneto-Centrifugally Launched Jets: I

We present simulations of the propagation of magnetized jets. This work differs from previous studies in that the cross-sectional distributions of the jets's state variables are derived from analytical models for magneto-centrifugal launching. The source is a magnetized rotator whose properties are specfied as boundary conditions. The jets in these simulations are considerably more complex than the ``top-hat''constant density etc. profiles used in previous work. We find that density and magnetic field stratification (with radius) in the jet leads to new behavior including the separation of an inner jet core from a low density collar. We find this {\it jet within a jet} structure, along with the magnetic stresses, leads to propagation behaviors not observed in previous simulation studies. Our methodology allows us to compare MHD jets from different types of sources whose properties could ultimately be derived from the behavior of the propagating jets.Comment: 42 pages, accepted by the Ap