Magnetically Controlled Spasmodic Accretion During Star Formation. II. Results

The problem of the late accretion phase of the evolution of an axisymmetric, isothermal magnetic disk surrounding a forming star has been formulated in a companion paper. The "central sink approximation" is used to circumvent the problem of describing the evolution inside the opaque central region for densities greater than 10^11 cm^-3 and radii smaller than a few AUs. Only the electrons are assumed to be attached to the magnetic field lines, and the effects of both negatively and positively charged grains are accounted for. After a mass of 0.1 solar mass accumulates in the central cell (forming star), a series of magnetically driven outflows and associated outward propagating shocks form in a quasi-periodic fashion. As a result, mass accretion onto the protostar occurs in magnetically controlled bursts. We refer to this process as spasmodic accretion. The shocks propagate outward with supermagnetosonic speeds. The period of dissipation and revival of the outflow decreases in time, as the mass accumulated in the central sink increases. We evaluate the contribution of ambipolar diffusion to the resolution of the magnetic flux problem of star formation during the accretion phase, and we find it to be very significant although not sufficient to resolve the entire problem yet. Ohmic dissipation is completely negligible in the disk during this phase of the evolution. The protostellar disk is found to be stable against interchange-like instabilities, despite the fact that the mass-to-flux ratio has temporary local maxima.Comment: Astrophysical Journal, in press. 29 pages, 13 figure

Demonstration of magnetic field tomography with starlight polarization towards a diffuse sightline of the ISM

The availability of large datasets with stellar distance and polarization information will enable a tomographic reconstruction of the (plane-of-the-sky-projected) interstellar magnetic field in the near future. We demonstrate the feasibility of such a decomposition within a small region of the diffuse ISM. We combine measurements of starlight (R-band) linear polarization obtained using the RoboPol polarimeter with stellar distances from the second Gaia data release. The stellar sample is brighter than 17 mag in the R band and reaches out to several kpc from the Sun. HI emission spectra reveal the existence of two distinct clouds along the line of sight. We decompose the line-of-sight-integrated stellar polarizations to obtain the mean polarization properties of the two clouds. The two clouds exhibit significant differences in terms of column density and polarization properties. Their mean plane-of-the-sky magnetic field orientation differs by 60 degrees. We show how our tomographic decomposition can be used to constrain our estimates of the polarizing efficiency of the clouds as well as the frequency dependence of the polarization angle of polarized dust emission. We also demonstrate a new method to constrain cloud distances based on this decomposition. Our results represent a preview of the wealth of information that can be obtained from a tomographic map of the ISM magnetic field.Comment: 25 pages, 14 figures, published in ApJ, data appear in journa

Scaling Relations of Dwarf Galaxies without Supernova-Driven Winds

Nearby dwarf galaxies exhibit tight correlations between their global stellar and dynamical properties, such as circular velocity, mass-to-light ratio, stellar mass, surface brightness, and metallicity. Such correlations have often been attributed to gas or metal-rich outflows driven by supernova energy feedback to the interstellar medium. We use high-resolution cosmological simulations of high-redshift galaxies with and without energy feedback, as well as analytic modeling, to investigate whether the observed correlations can arise without supernova-driven outflows. We find that the simulated dwarf galaxies exhibit correlations similar to those observed as early as z~10, regardless of whether supernova feedback is included. We also show that the correlations can be well reproduced by our analytic model that accounts for realistic gas inflow but assumes no outflows, and star formation rate obeying the Kennicutt-Schmidt law with a critical density threshold. We argue that correlations in simulated galaxies arise due to the increasingly inefficient conversion of gas into stars in low-mass dwarf galaxies rather than supernova-driven outflows. We also show that the decrease of the observed effective yield in low-mass objects, often used as an indicator of gas and metal outflows, can be reasonably reproduced in our simulations without outflows. We show that this trend can arise if a significant fraction of metals in small galaxies is spread to the outer regions of the halo outside the stellar extent via mixing. In this case the effective yield can be significantly underestimated if only metals within the stellar radius are taken into account. Measurements of gas metallicity in the outskirts of gaseous disks of dwarfs would thus provide a key test of such explanation.Comment: accepted for publication in ApJ, 19 pages, 12 figures, uses emulateapj

Comparison of prestellar core elongations and large-scale molecular cloud structures in the Lupus 1 region

Turbulence and magnetic fields are expected to be important for regulating molecular cloud formation and evolution. However, their effects on sub-parsec to 100 parsec scales, leading to the formation of starless cores, are not well understood. We investigate the prestellar core structure morphologies obtained from analysis of the Herschel-SPIRE 350 mum maps of the Lupus I cloud. This distribution is first compared on a statistical basis to the large-scale shape of the main filament. We find the distribution of the elongation position angle of the cores to be consistent with a random distribution, which means no specific orientation of the morphology of the cores is observed with respect to the mean orientation of the large-scale filament in Lupus I, nor relative to a large-scale bent filament model. This distribution is also compared to the mean orientation of the large-scale magnetic fields probed at 350 mum with the Balloon-borne Large Aperture Telescope for Polarimetry during its 2010 campaign. Here again we do not find any correlation between the core morphology distribution and the average orientation of the magnetic fields on parsec scales. Our main conclusion is that the local filament dynamics---including secondary filaments that often run orthogonally to the primary filament---and possibly small-scale variations in the local magnetic field direction, could be the dominant factors for explaining the final orientation of each core

Impact of Supernova feedback on the Tully-Fisher relation

Recent observational results found a bend in the Tully-Fisher Relation in such a way that low mass systems lay below the linear relation described by more massive galaxies. We intend to investigate the origin of the observed features in the stellar and baryonic Tully-Fisher relations and analyse the role played by galactic outflows on their determination. Cosmological hydrodynamical simulations which include Supernova feedback were performed in order to follow the dynamical evolution of galaxies. We found that Supernova feedback is a fundamental process in order to reproduce the observed trends in the stellar Tully-Fisher relation. Simulated slow rotating systems tend to have lower stellar masses than those predicted by the linear fit to the massive end of the relation, consistently with observations. This feature is not present if Supernova feedback is turned off. In the case of the baryonic Tully-Fisher relation, we also detect a weaker tendency for smaller systems to lie below the linear relation described by larger ones. This behaviour arises as a result of the more efficient action of Supernovae in the regulation of the star formation process and in the triggering of powerful galactic outflows in shallower potential wells which may heat up and/or expel part of the gas reservoir.Comment: 10 pages, 9 figures, accepted for publication in A&

RoboPol: First season rotations of optical polarization plane in blazars

We present first results on polarization swings in optical emission of blazars obtained by RoboPol, a monitoring program of an unbiased sample of gamma-ray bright blazars specially designed for effective detection of such events. A possible connection of polarization swing events with periods of high activity in gamma rays is investigated using the dataset obtained during the first season of operation. It was found that the brightest gamma-ray flares tend to be located closer in time to rotation events, which may be an indication of two separate mechanisms responsible for the rotations. Blazars with detected rotations have significantly larger amplitude and faster variations of polarization angle in optical than blazars without rotations. Our simulations show that the full set of observed rotations is not a likely outcome (probability $\le 1.5 \times 10^{-2}$) of a random walk of the polarization vector simulated by a multicell model. Furthermore, it is highly unlikely ($\sim 5 \times 10^{-5}$) that none of our rotations is physically connected with an increase in gamma-ray activity.Comment: 16 pages, 9 figure