Protostars: Forges of cosmic rays?

Galactic cosmic rays (CR) are particles presumably accelerated in supernova remnant shocks that propagate in the interstellar medium up to the densest parts of molecular clouds, losing energy and their ionisation efficiency because of the presence of magnetic fields and collisions with molecular hydrogen. Recent observations hint at high levels of ionisation and at the presence of synchrotron emission in protostellar systems, which leads to an apparent contradiction. We want to explain the origin of these CRs accelerated within young protostars as suggested by observations. Our modelling consists of a set of conditions that has to be satisfied in order to have an efficient CR acceleration through diffusive shock acceleration. We analyse three main acceleration sites, then we follow the propagation of these particles through the protostellar system up to the hot spot region. We find that jet shocks can be strong accelerators of CR protons, which can be boosted up to relativistic energies. Other promising acceleration sites are protostellar surfaces, where shocks caused by impacting material during the collapse phase are strong enough to accelerate CR protons. In contrast, accretion flow shocks are too weak to efficiently accelerate CRs. Though CR electrons are weakly accelerated, they can gain a strong boost to relativistic energies through re-acceleration in successive shocks. We suggest a mechanism able to accelerate both CR protons and electrons through the diffusive shock acceleration mechanism, which can be used to explain the high ionisation rate and the synchrotron emission observed towards protostellar sources. The existence of an internal source of energetic particles can have a strong and unforeseen impact on the ionisation of the protostellar disc, on the star and planet formation processes, and on the formation of pre-biotic molecules.Comment: 22 pages, 15 figures, accepted by Astronomy and Astrophysic

Cosmic ray moderation of the thermal instability

We apply the Hermite-Bieler theorem in the analysis of the effect of cosmic rays on the thermal stability of an initially uniform, static background. The cosmic rays were treated in a fluid approximation and the diffusion coefficient was assumed to be constant in time and space. The inclusion of cosmic rays does not alter the criterion for the thermal stability of a medium subjected to isobaric perturbations. It does alter the criteria for the stability of a medium perturbed by small amplitude sound waves. In the limit of a high background cosmic ray pressure to thermal pressure ratio, the instability in response to high frequency sound waves is suppressed

Galactic interstellar filaments as probed by LOFAR and Planck

Recent Low Frequency Array (LOFAR) observations at 115-175 MHz of a field at medium Galactic latitudes (centered at the bright quasar 3C196) have shown striking filamentary structures in polarization that extend over more than 4 degrees across the sky. In addition, the Planck satellite has released full sky maps of the dust emission in polarization at 353GHz. The LOFAR data resolve Faraday structures along the line of sight, whereas the Planck dust polarization maps probe the orientation of the sky projected magnetic field component. Hence, no apparent correlation between the two is expected. Here we report a surprising, yet clear, correlation between the filamentary structures, detected with LOFAR, and the magnetic field orientation, probed by the Planck satellite. This finding points to a common, yet unclear, physical origin of the two measurements in this specific area in the sky. A number of follow-up multi- frequency studies are proposed to shed light on this unexpected finding.Comment: 6 pages, 4 figures, accepted for publication in MNRAS Letter

A Survey of Extragalactic Faraday Rotation at High Galactic Latitude: The Vertical Magnetic Field of the Milky Way towards the Galactic Poles

We present a study of the vertical magnetic field of the Milky Way towards the Galactic poles, determined from observations of Faraday rotation toward more than 1000 polarized extragalactic radio sources at Galactic latitudes |b| > 77 degs, using the Westerbork Radio Synthesis Telescope and the Australia Telescope Compact Array. We find median rotation measures (RMs) of 0.0 +/- 0.5 rad/m^2 and +6.3 +/- 0.7 rad/m^2 toward the north and south Galactic poles, respectively, demonstrating that there is no coherent vertical magnetic field in the Milky Way at the Sun's position. If this is a global property of the Milky Way's magnetism, then the lack of symmetry across the disk rules out pure dipole or quadrupole geometries for the Galactic magnetic field. The angular fluctuations in RM seen in our data show no preferred scale within the range ~ 0.1 to 25 degs. The observed standard deviation in RM of ~ 9 rad/m^2 then implies an upper limit of ~1microGauss on the strength of the random magnetic field in the warm ionized medium at high Galactic latitudes.Comment: 38 pages, 7 figures, 3 tables Accepted for publication in ApJ, Electronic versions of Tables 1 and 2 are available via email from the first autho

Two-dimensional AMR simulations of colliding flows

Colliding flows are a commonly used scenario for the formation of molecular clouds in numerical simulations. Due to the thermal instability of the warm neutral medium, turbulence is produced by cooling. We carry out a two-dimensional numerical study of such colliding flows in order to test whether statistical properties inferred from adaptive mesh refinement (AMR) simulations are robust with respect to the applied refinement criteria. We compare probability density functions of various quantities as well as the clump statistics and fractal dimension of the density fields in AMR simulations to a static-grid simulation. The static grid with 2048^2 cells matches the resolution of the most refined subgrids in the AMR simulations. The density statistics is reproduced fairly well by AMR. Refinement criteria based on the cooling time or the turbulence intensity appear to be superior to the standard technique of refinement by overdensity. Nevertheless, substantial differences in the flow structure become apparent. In general, it is difficult to separate numerical effects from genuine physical processes in AMR simulations.Comment: 6 pages, 6 figures, submitted to A&

Galaxies in box: A simulated view of the interstellar medium

We review progress in the development of physically realistic three dimensional simulated models of the galaxy.We consider the scales from star forming molecular clouds to the full spiral disc. Models are computed using hydrodynamic (HD) or magnetohydrodynamic (MHD) equations and may include cosmic ray or tracer particles. The range of dynamical scales between the full galaxy structure and the turbulent scales of supernova (SN) explosions and even cloud collapse to form stars, make it impossible with current computing tools and resources to resolve all of these in one model. We therefore consider a hierarchy of models and how they can be related to enhance our understanding of the complete galaxy.Comment: Chapter in Large Scale Magnetic Fields in the Univers

Accurate positions for the ULXs NGC 7319-X4 and NGC 5474-X1 and limiting magnitudes for their optical counterparts

In this paper we report accurate Chandra positions for two ultraluminous X-ray sources: NGC 7319-X4 at Right Ascension (RA) = 339.02917(2) deg, Declination (Dec) = 33.97476(2) deg and NGC 5474-X1 at RA = 211.24859(3) deg, Dec = 53.63584(3) deg. We perform bore-sight corrections on the Chandra X-ray Satellite observations of these sources to get to these accurate positions of the X-ray sources and match these positions with archival optical data from the Wide Field and Planetary Camera 2 on board the Hubble Space Telescope. We do not find the optical counterparts: the limiting absolute magnitudes of the observations in the WFPC2 standard magnitude system are B = -7.9, V = -8.7 and I = -9.3 for NGC 7319-X4 and U = -6.4 for NGC 5474-X1. We report on the X-ray spectral properties and we find evidence for X-ray variability in NGC 5474-X1. Finally, we briefly discuss several options for the nature of these ULXs.Comment: 8 pages, 4 figures, accepted for publication in MNRA

A search for linear polarization in the active galactic nucleus 3C 84 at 239 and 348 GHz

We report a search for linear polarization in the active galactic nucleus (AGN) 3C 84 (NGC 1275) at observed frequencies of 239 GHz and 348 GHz, corresponding to rest-frame frequencies of 243 GHz and 354 GHz. We collected polarization data with the IRAM Plateau de Bure Interferometer via Earth rotation polarimetry. We do not detect linear polarization. Our analysis finds 3-sigma upper limits on the degree of polarization of 0.5% and 1.9% at 239 GHz and 348 GHz, respectively. We regard the influence of Faraday conversion as marginal, leading to expected circular polarizations <0.3%. Assuming depolarization by a local Faraday screen, we constrain the rotation measure, as well as the fluctuations therein, to be 10^6 rad/m^2. From this we estimate line-of-sight magnetic field strengths of >100 microG. Given the physical dimensions of 3C 84 and its observed structure, the Faraday screen appears to show prominent small-scale structure, with \DeltaRM > 10^6 rad/m^2 on projected spatial scales <1 pc.Comment: 7 pages, 4 figures. Accepted by MNRA

On the Origin of Cosmic Magnetic Fields

We review the literature concerning how the cosmic magnetic fields pervading nearly all galaxies actually got started. some observational evidence involves the chemical abundance of the light elements Be and B, while another one is based on strong magnetic fields seen in high red shift galaxies. Seed fields, whose strength is of order 10^{-20} gauss, easily sprung up in the era preceding galaxy formation. Several mechanisms are proposed to amplify these seed fields to microgauss strengths. The standard mechanism is the Alpha-Omega dynamo theory. It has a major difficulty that makes unlikely to provide the sole origin. The difficulty is rooted in the fact that the total flux is constant. This implies that flux must be removed from the galactic discs. This requires that the field and flux be separated, for otherwise interstellar mass must be removed from the deep galactic gravitational and then their strength increased by the alpha omega theory.Comment: 90 pages and 6 figures; accepted for publication in Reports of Progress in Physics as an invited revie