The magnetic fields of spiral galaxies are so strong that they cannot be primordial. Their typical values are over one billion times higher than any value predicted for the early Universe. Explaining this immense growth and incorporating it in galaxy evolution theories is one of the long-standing challenges in astrophysics. So far, the most successful theory for the sustained growth of the galactic magnetic field is the alpha-omega dynamo. This theory predicts a characteristic dipolar or quadrupolar morphology for the galactic magnetic field, which has been observed in external galaxies. However, so far, there has been no direct demonstration of a mean-field dynamo operating in direct, multi-physics simulations of spiral galaxies. We do so in this work. We employ numerical models of isolated, star-forming spiral galaxies that include a magnetized gaseous disk, a dark matter halo, stars, and stellar feedback. Naturally, the resulting magnetic field has a complex morphology that includes a strong random component. Using a smoothing of the magnetic field on small scales, we are able to separate the mean from the turbulent component and analyze them individually. We find that a mean-field dynamo naturally occurs as a result of the dynamical evolution of the galaxy and amplifies the magnetic field by an order of magnitude over half a Gyr. Despite the highly dynamical nature of these models, the morphology of the mean component of the field is identical to analytical predictions. This result underlines the importance of the mean-field dynamo in galactic evolution. Moreover, by demonstrating the natural growth of the magnetic field in a complex galactic environment, it brings us a step closer to understanding the cosmic origin of magnetic fields.

We use N-body simulations to examine whether a characteristic turnaround radius, as predicted from the spherical collapse model in a $\rm {\Lambda CDM}$ Universe, can be meaningfully identified for galaxy clusters, in the presence of full three-dimensional effects. We use The Dark Sky Simulations and Illustris-TNG dark-matter--only cosmological runs to calculate radial velocity profiles around collapsed structures, extending out to many times the virial radius $R_{200}$. There, the turnaround radius can be unambiguously identified as the largest non-expanding scale around a center of gravity. We find that: (a) Indeed, a single turnaround scale can meaningfully describe strongly non-spherical structures. (b) For halos of masses $M_{200}>10^{13}M_\odot$, the turnaround radius $R_{ta}$ scales with the enclosed mass $M_{ta}$ as $M_{ta}^{1/3}$, as predicted by the spherical collapse model. (c) The deviation of $R_{ta}$ in simulated halos from the spherical collapse model prediction is insensitive to halo asphericity. Rather, it is sensitive to the tidal forces due to massive neighbors when such are present. (d) Halos exhibit a characteristic average density within the turnaround scale. This characteristic density is dependent on cosmology and redshift. For the present cosmic epoch and for concordance cosmological parameters ($\Omega_m \sim 0.7$; $\Omega_\Lambda \sim 0.3$) turnaround structures exhibit an average matter density contrast with the background Universe of $\delta \sim 11$. Thus $R_{ta}$ is equivalent to $R_{11}$ -- in a way analogous to defining the "virial" radius as $R_{200}$ -- with the advantage that $R_{11}$ is shown in this work to correspond to a kinematically relevant scale in N-body simulations.

Polarization measurements done using Imaging Polarimeters such as the Robotic Polarimeter are very sensitive to the presence of artefacts in images. Artefacts can range from internal reflections in a telescope to satellite trails that could contaminate an area of interest in the image. With the advent of wide-field polarimetry surveys, it is imperative to develop methods that automatically flag artefacts in images. In this paper, we implement a Convolutional Neural Network to identify the most dominant artefacts in the images. We find that our model can successfully classify sources with 98\% true positive and 97\% true negative rates. Such models, combined with transfer learning, will give us a running start in artefact elimination for near-future surveys like WALOP.

We report the first detection of galactic spiral structure by means of thermal emission from magnetically aligned dust grains. Our 89 $\mu$m polarimetric imaging of NGC 1068 with the High-resolution Airborne Wideband Camera/Polarimeter (HAWC+) on NASA's Stratospheric Observatory for Infrared Astronomy (SOFIA) also sheds light on magnetic field structure in the vicinity of the galaxy's inner-bar and active galactic nucleus (AGN). We find correlations between the 89 $\mu$m magnetic field vectors and other tracers of spiral arms, and a symmetric polarization pattern as a function of the azimuthal angle arising from the projection and inclination of the disk field component in the plane of the sky. The observations can be fit with a logarithmic spiral model with pitch angle of $16.9^{+2.7}_{-2.8}$$^{\circ}$ and a disk inclination of $48\pm2^{\circ}$. We infer that the bulk of the interstellar medium from which the polarized dust emission originates is threaded by a magnetic field that closely follows the spiral arms. Inside the central starburst disk ($<1.6$ kpc), the degree of polarization is found to be lower than for far-infrared sources in the Milky Way, and has minima at the locations of most intense star formation near the outer ends of the inner-bar. Inside the starburst ring, the field direction deviates from the model, becoming more radial along the leading edges of the inner-bar. The polarized flux and dust temperature peak $\sim 3-6$" NE of the AGN at the location of a bow shock between the AGN outflow and the surrounding interstellar medium, but the AGN itself is weakly polarized ($< 1$%) at both 53 and 89 \um.

We propose in this White Paper a concept for a space experiment using cold atoms to search for ultra-light dark matter, and to detect gravitational waves in the frequency range between the most sensitive ranges of LISA and the terrestrial LIGO/Virgo/KAGRA/INDIGO experiments. This interdisciplinary experiment, called Atomic Experiment for Dark Matter and Gravity Exploration (AEDGE), will also complement other planned searches for dark matter, and exploit synergies with other gravitational wave detectors. We give examples of the extended range of sensitivity to ultra-light dark matter offered by AEDGE, and how its gravitational-wave measurements could explore the assembly of super-massive black holes, first-order phase transitions in the early universe and cosmic strings. AEDGE will be based upon technologies now being developed for terrestrial experiments using cold atoms, and will benefit from the space experience obtained with, e.g., LISA and cold atom experiments in microgravity. This paper is based on a submission (v1) in response to the Call for White Papers for the Voyage 2050 long-term plan in the ESA Science Programme. ESA limited the number of White Paper authors to 30. However, in this version (v2) we have welcomed as supporting authors participants in the Workshop on Atomic Experiments for Dark Matter and Gravity Exploration held at CERN: ({\tt https://indico.cern.ch/event/830432/}), as well as other interested scientists, and have incorporated additional material.

We report on polarimetric maps made with HAWC+/SOFIA toward Rho Oph A, the densest portion of the Rho Ophiuchi molecular complex. We employed HAWC+ bands C (89 $\mu$m) and D (154 $\mu$m). The slope of the polarization spectrum was investigated by defining the quantity R_DC = p_D/p_C, where p_C and p_D represent polarization degrees in bands C and D, respectively. We find a clear correlation between R_DC and the molecular hydrogen column density across the cloud. A positive slope (R_DC > 1) dominates the lower density and well illuminated portions of the cloud, that are heated by the high mass star Oph S1, whereas a transition to a negative slope (R_DC < 1) is observed toward the denser and less evenly illuminated cloud core. We interpret the trends as due to a combination of: (1) Warm grains at the cloud outskirts, which are efficiently aligned by the abundant exposure to radiation from Oph S1, as proposed in the radiative torques theory; and (2) Cold grains deep in the cloud core, which are poorly aligned due to shielding from external radiation. To assess this interpretation, we developed a very simple toy model using a spherically symmetric cloud core based on Herschel data, and verified that the predicted variation of R_DC is consistent with the observations. This result introduces a new method that can be used to probe the grain alignment efficiency in molecular clouds, based on the analysis of trends in the far-infrared polarization spectrum.

SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This paper presents an overview of the main science drivers for B-BOP, including high dynamic range polarimetric imaging of the cold interstellar medium (ISM) in both our Milky Way and nearby galaxies. Thanks to a cooled telescope, B-BOP will deliver wide-field 100-350 micron images of linearly polarized dust emission in Stokes Q and U with a resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges comparable to those achieved by Herschel images of the cold ISM in total intensity (Stokes I). The B-BOP 200 micron images will also have a factor ~30 higher resolution than Planck polarization data. This will make B-BOP a unique tool for characterizing the statistical properties of the magnetized interstellar medium and probing the role of magnetic fields in the formation and evolution of the interstellar web of dusty molecular filaments giving birth to most stars in our Galaxy. B-BOP will also be a powerful instrument for studying the magnetism of nearby galaxies and testing galactic dynamo models, constraining the physics of dust grain alignment, informing the problem of the interaction of cosmic rays with molecular clouds, tracing magnetic fields in the inner layers of protoplanetary disks, and monitoring accretion bursts in embedded protostars.

Galactic dust emission is polarized at unexpectedly high levels, as revealed by Planck. The origin of the observed $\simeq 20\%$ polarization fractions can be identified by characterizing the properties of optical starlight polarization in a region with maximally polarized dust emission. We measure the R-band linear polarization of 22 stars in a region with a submillimeter polarization fraction of $\simeq 20%$. A subset of 6 stars is also measured in the B, V and I bands to investigate the wavelength dependence of polarization. We find that starlight is polarized at correspondingly high levels. Through multiband polarimetry we find that the high polarization fractions are unlikely to arise from unusual dust properties, such as enhanced grain alignment. Instead, a favorable magnetic field geometry is the most likely explanation, and is supported by observational probes of the magnetic field morphology. The observed starlight polarization exceeds the classical upper limit of $\left[p_V/E\left(B-V\right)\right]_{\rm max} = 9$%mag$^{-1}$ and is at least as high as 13%mag$^{-1}$ that was inferred from a joint analysis of Planck data, starlight polarization and reddening measurements. Thus, we confirm that the intrinsic polarizing ability of dust grains at optical wavelengths has long been underestimated.

We present the design and performance of RoboPol, a four-channel optical polarimeter operating at the Skinakas Observatory in Crete, Greece. RoboPol is capable of measuring both relative linear Stokes parameters $q$ and $u$ (and the total intensity $I$) in one sky exposure. Though primarily used to measure the polarization of point sources in the R-band, the instrument features additional filters (B, V and I), enabling multi-wavelength imaging polarimetry over a large field of view (13.6' $\times$ 13.6'). We demonstrate the accuracy and stability of the instrument throughout its five years of operation. Best performance is achieved within the central region of the field of view and in the R band. For such measurements the systematic uncertainty is below 0.1% in fractional linear polarization, $p$ (0.05% maximum likelihood). Throughout all observing seasons the instrumental polarization varies within 0.1% in $p$ and within 1$^\circ$ in polarization angle.

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.

We report new polarimetric and photometric maps of the massive star-forming region OMC-1 using the HAWC+ instrument on the Stratospheric Observatory for Infrared Astronomy (SOFIA). We present continuum polarimetric and photometric measurements of this region at 53, 89, 154, and 214 microns at angular resolutions of 5.1, 7.9, 14.0, and 18.7 arcseconds for the four bands, respectively. The photometric maps enable the computation of improved SEDs for the region. We find that at the longer wavelengths, the inferred magnetic field configuration matches the `hourglass' configuration seen in previous studies, indicating magnetically-regulated star formation. The field morphology differs at the shorter wavelengths. The magnetic field inferred at these wavelengths traces the bipolar structure of the explosive Becklin-Neugebauer (BN)/Kleinman-Low (KL) outflow emerging from OMC-1 behind the Orion Nebula. Using statistical methods to estimate the field strength in the region, we find that the explosion dominates the magnetic field near the center of the feature. Farther out, the magnetic field is close to energetic equilibrium with the ejecta and may be providing confinement to the explosion. The correlation between polarization fraction and the local polarization angle dispersion indicates that the depolarization as a function of unpolarized intensity is a result of intrinsic field geometry as opposed to decreases in grain alignment efficiency in denser regions.

We present Far-Infrared polarimetry observations of M82 at 53 and $154~\mu \rm{m}$ and NGC 253 at $89~\mu \rm{m}$, which were taken with HAWC+ in polarimetry mode on the Stratospheric Observatory for Infrared Astronomy (SOFIA). The polarization of M82 at $53~\mu \rm{m}$ clearly shows a magnetic field geometry perpendicular to the disk in the hot dust emission. For M82 the polarization at $154~\mu \rm{m}$ shows a combination of field geometry perpendicular to the disk in the nuclear region, but closer to parallel to the disk away from the nucleus. The fractional polarization at $53~\mu \rm{m}$ $(154~\mu \rm{m})$ ranges from 7% (3%) off nucleus to 0.5% (0.3%) near the nucleus. A simple interpretation of the observations of M82 invokes a massive polar outflow, dragging the field along, from a region $\sim 700$~pc in diameter that has entrained some of the gas and dust, creating a vertical field geometry seen mostly in the hotter $(53~\mu \rm{m})$ dust emission. This outflow sits within a larger disk with a more typical planar geometry that more strongly contributes to the cooler $(154~\mu \rm{m})$ dust emission. For NGC 253, the polarization at $89~\mu \rm{m}$ is dominated by a planar geometry in the tilted disk, with weak indication of a vertical geometry above and below the plane from the nucleus. The polarization observations of NGC 253 at $53~\mu \rm{m}$ were of insufficient S/N for detailed analysis.

PASIPHAE (the Polar-Areas Stellar Imaging in Polarization High-Accuracy Experiment) is an optopolarimetric survey aiming to measure the linear polarization from millions of stars, and use these to create a three-dimensional tomographic map of the magnetic field threading dust clouds within the Milky Way. This map will provide invaluable information for future CMB B-mode experiments searching for inflationary gravitational waves, providing unique information regarding line-of-sight integration effects. Optical polarization observations of a large number of stars at known distances, tracing the same dust that emits polarized microwaves, can map the magnetic field between them. The Gaia mission is measuring distances to a billion stars, providing an opportunity to produce a tomographic map of Galactic magnetic field directions, using optical polarization of starlight. Such a map will not only boost CMB polarization foreground removal, but it will also have a profound impact in a wide range of astrophysical research, including interstellar medium physics, high-energy astrophysics, and evolution of the Galaxy. Taking advantage of the novel technology implemented in our high-accuracy Wide-Area Linear Optical Polarimeters (WALOPs) currently under construction at IUCAA, India, we will engage in a large-scale optopolarimetric program that can meet this challenge: a survey of both northern and southern Galactic polar regions targeted by CMB experiments, covering over 10,000 square degrees, which will measure linear optical polarization of over 360 stars per square degree (over 3.5 million stars, a 1000-fold increase over the state of the art). The survey will be conducted concurrently from the South African Astronomical Observatory in Sutherland, South Africa in the southern hemisphere, and the Skinakas Observatory in Crete, Greece, in the north.

Recent studies of the diffuse parts of molecular clouds have revealed the presence of parallel, ordered low-density filaments termed striations. Flows along magnetic field lines, Kelvin-Helmholtz instabilities and hydromagnetic waves are amongst the various formation mechanisms proposed. Through a synergy of observational, numerical and theoretical analysis, previous studies singled out the hydromagnetic waves model as the only one that can account for the observed properties of striations. Based on the predictions of that model, we develop here a method for measuring the temporal evolution of striations through a combination of molecular and dust continuum observations. Our method allows us to not only probe temporal variations in molecular clouds but also estimate the strength of both the ordered and fluctuating components of the magnetic field projected on the plane-of-the-sky. We benchmark our new method against chemical and radiative transfer effects through two-dimensional magnetohydrodynamic simulations coupled with non-equilibrium chemical modelling and non-local thermodynamic equilibrium line radiative transfer. We find good agreement between theoretical predictions, simulations and observations of striations in the Taurus molecular cloud. We find a value of $\rm{27 \pm 7} ~\rm{\mu G}$ for the plane-of-sky magnetic field, in agreement with previous estimates via the Davis-Chandrasekhar-Fermi method, and a ratio of fluctuating to ordered component of the magnetic field of $\sim$ 10\%.

We conduct a small-scale pathfinding survey designed to identify the average polarization properties of the diffuse ISM locally at the lowest dust content regions. We perform deep optopolarimetric surveys within three $\sim 15' \times 15'$ regions located at $b > 48^\circ$, using the RoboPol instrument. The observed samples of stars are photometrically complete to $\sim$16 mag in the R-band. The selected regions exhibit low dust emission at 353 GHz and low total reddening compared to the majority of high-latitude sightlines. We measure the level of systematic uncertainty for all observing epochs and find it to be 0.1\% in fractional linear polarization, $p$. The majority of individual stellar measurements are non-detections. However, our survey strategy enables us to locate the mean fractional linear polarization $p_{mean}$ in each of the three regions. The region with lowest dust content yields $p_{mean}=(0.054 \pm 0.038) \%$, not significantly different from zero. We find significant detections for the remaining two regions of: $p_{mean}=(0.113 \pm 0.036) \%$ and $p_{mean}=(0.208 \pm 0.044) \%$. Using a Bayesian approach we provide upper limits on the intrinsic spread of the small-scale distributions of $q$ and $u$. At the detected $p_{mean}$ levels, the determination of the systematic uncertainty is critical for the reliability of the measurements. We verify the significance of our detections with statistical tests, accounting for all sources of uncertainty. Using publicly available HI emission data, we identify the velocity components that most likely account for the observed $p_{mean}$ and find their morphologies to be misaligned with the orientation of the mean plane-of-sky magnetic field at a spatial resolution of 10$\arcmin$. We find indications that the standard upper envelope of $p$ with reddening underestimates the maximum $p$ at very low E(B-V) ($\leq 0.01$ mag).

We describe PyRaTE, a new, non-local thermodynamic equilibrium (non-LTE) line radiative transfer code developed specifically for post-processing astrochemical simulations. Population densities are estimated using the escape probability method. When computing the escape probability, the optical depth is calculated towards all directions with density, molecular abundance, temperature and velocity variations all taken into account. A very easy-to-use interface, capable of importing data from simulations outputs performed with all major astrophysical codes, is also developed. The code is written in Python using an `embarrassingly parallel' strategy and can handle all geometries and projection angles. We benchmark the code by comparing our results with those from RADEX (van der Tak et al. 2007) and against analytical solutions and present case studies using hydrochemical simulations. The code is available on GitHub (https://github.com/ArisTr/PyRaTE).

Stars and planets are formed inside dense interstellar molecular clouds, by processes imprinted on the 3-dimensional (3D) morphology of the clouds. Determining the 3D structure of interstellar clouds remains challenging, due to projection effects and difficulties measuring their extent along the line of sight. We report the detection of normal vibrational modes in the isolated interstellar cloud Musca, allowing determination of the 3D physical dimensions of the cloud. Musca is found to be vibrating globally, with the characteristic modes of a sheet viewed edge-on, not a filament as previously supposed. We reconstruct the physical properties of Musca through 3D magnetohydrodynamic simulations, reproducing the observed normal modes and confirming a sheet-like morphology.

We use results of our 3 year polarimetric monitoring program to investigate the previously suggested connection between rotations of the polarization plane in the optical emission of blazars and their gamma-ray flares in the GeV band. The homogeneous set of 40 rotation events in 24 sources detected by {\em RoboPol} is analysed together with the gamma-ray data provided by {\em Fermi}-LAT. We confirm that polarization plane rotations are indeed related to the closest gamma-ray flares in blazars and the time lags between these events are consistent with zero. Amplitudes of the rotations are anticorrelated with amplitudes of the gamma-ray flares. This is presumably caused by higher relativistic boosting (higher Doppler factors) in blazars that exhibit smaller amplitude polarization plane rotations. Moreover, the time scales of rotations and flares are marginally correlated.

Filaments in Herschel molecular cloud images are found to exhibit a "characteristic width". This finding is in tension with spatial power spectra of the data, which show no indication of this characteristic scale. We demonstrate that this discrepancy is a result of the methodology adopted for measuring filament widths. First, we perform the previously used analysis technique on artificial scale-free data, and obtain a peaked width distribution of filament-like structures. Next, we repeat the analysis on three Herschel maps and reproduce the narrow distribution of widths found in previous studies $-$ when considering the average width of each filament. However, the distribution of widths measured at all points along a filament spine is broader than the distribution of mean filament widths, indicating that the narrow spread (interpreted as a "characteristic" width) results from averaging. Furthermore, the width is found to vary significantly from one end of a filament to the other. Therefore, the previously identified peak at 0.1 pc cannot be understood as representing the typical width of filaments. We find an alternative explanation by modelling the observed width distribution as a truncated power-law distribution, sampled with uncertainties. The position of the peak is connected to the lower truncation scale and is likely set by the choice of parameters used in measuring filament widths. We conclude that a "characteristic" width of filaments is not supported by the available data.

The magnetic field (B-field) of the starless dark cloud L1544 has been studied using near-infrared (NIR) background starlight polarimetry (BSP) and archival data in order to characterize the properties of the plane-of-sky B-field. NIR linear polarization measurements of over 1,700 stars were obtained in the H-band and 201 of these were also measured in the K-band. The NIR BSP properties are correlated with reddening, as traced using the RJCE (H-M) method, and with thermal dust emission from the L1544 cloud and envelope seen in Herschel maps. The NIR polarization position angles change at the location of the cloud and exhibit their lowest dispersion of position angles there, offering strong evidence that NIR polarization traces the plane-of-sky B-field of L1544. In this paper, the uniformity of the plane-of-sky B-field in the envelope region of L1544 is quantitatively assessed. This allowed evaluating the approach of assuming uniform field geometry when measuring relative mass-to-flux ratios in the cloud envelope and core based on averaging of the envelope radio Zeeman observations, as in Crutcher et al. (2009). In L1544, the NIR BSP shows the envelope B-field to be significantly non-uniform and likely not suitable for averaging Zeeman properties without treating intrinsic variations. Deeper analyses of the NIR BSP and related data sets, including estimates of the B-field strength and testing how it varies with position and gas density, are the subjects of later papers in this series.

We present average R-band optopolarimetric data, as well as variability parameters, from the first and second RoboPol observing season. We investigate whether gamma- ray--loud and gamma-ray--quiet blazars exhibit systematic differences in their optical polarization properties. We find that gamma-ray--loud blazars have a systematically higher polarization fraction (0.092) than gamma-ray--quiet blazars (0.031), with the hypothesis of the two samples being drawn from the same distribution of polarization fractions being rejected at the 3{\sigma} level. We have not found any evidence that this discrepancy is related to differences in the redshift distribution, rest-frame R-band lu- minosity density, or the source classification. The median polarization fraction versus synchrotron-peak-frequency plot shows an envelope implying that high synchrotron- peaked sources have a smaller range of median polarization fractions concentrated around lower values. Our gamma-ray--quiet sources show similar median polarization fractions although they are all low synchrotron-peaked. We also find that the random- ness of the polarization angle depends on the synchrotron peak frequency. For high synchrotron-peaked sources it tends to concentrate around preferred directions while for low synchrotron-peaked sources it is more variable and less likely to have a pre- ferred direction. We propose a scenario which mediates efficient particle acceleration in shocks and increases the helical B-field component immediately downstream of the shock.

Dust continuum and molecular observations of the low column density parts of molecular clouds have revealed the presence of elongated structures which appear to be well aligned with the magnetic field. These so-called striations are usually assumed to be streams that flow towards or away from denser regions. We perform ideal magnetohydrodynamic (MHD) simulations adopting four models that could account for the formation of such structures. In the first two models striations are created by velocity gradients between ambient, parallel streamlines along magnetic field lines. In the third model striations are formed as a result of a Kelvin-Helmholtz instability perpendicular to field lines. Finally, in the fourth model striations are formed from the nonlinear coupling of MHD waves due to density inhomogeneities. We assess the validity of each scenario by comparing the results from our simulations with previous observational studies and results obtained from the analysis of CO (J = 1 - 0) observations from the Taurus molecular cloud. We find that the first three models cannot reproduce the density contrast and the properties of the spatial power spectrum of a perpendicular cut to the long axes of striations. We conclude that the nonlinear coupling of MHD waves is the most probable formation mechanism of striations.

We present a new set of optical polarization plane rotations in blazars, observed during the third year of operation of RoboPol. The entire set of rotation events discovered during three years of observations is analysed with the aim of determining whether these events are inherent in all blazars. It is found that the frequency of the polarization plane rotations varies widely among blazars. This variation cannot be explained either by a difference in the relativistic boosting or by selection effects caused by a difference in the average fractional polarization. We conclude that the rotations are characteristic of a subset of blazars and that they occur as a consequence of their intrinsic properties.

In diffuse molecular clouds, possible precursors of star-forming clouds, the effect of the magnetic field is unclear. In this work we compare the orientations of filamentary structures in the Polaris Flare, as seen through dust emission by Herschel, to the plane-of-the-sky magnetic field orientation ($\rm B_{pos}$) as revealed by stellar optical polarimetry with RoboPol. Dust structures in this translucent cloud show a strong preference for alignment with $\rm B_{pos}$. 70 % of field orientations are consistent with those of the filaments (within 30$^\circ$). We explore the spatial variation of the relative orientations and find it to be uncorrelated with the dust emission intensity and correlated to the dispersion of polarization angles. Concentrating in the area around the highest column density filament, and in the region with the most uniform field, we infer the $\rm B_{pos}$ strength to be 24 $-$ 120 $\mu$G. Assuming that the magnetic field can be decomposed into a turbulent and an ordered component, we find a turbulent-to-ordered ratio of 0.2 $-$ 0.8, implying that the magnetic field is dynamically important, at least in these two areas. We discuss implications on the 3D field properties, as well as on the distance estimate of the cloud.

We present a 3 mm spectral line and continuum survey of L1451 in the Perseus Molecular Cloud. These observations are from the CARMA Large Area Star Formation Survey (CLASSy), which also imaged Barnard 1, NGC 1333, Serpens Main and Serpens South. L1451 is the survey region with the lowest level of star formation activity---it contains no confirmed protostars. HCO+, HCN, and N2H+ (J=1-0) are all detected throughout the region, with HCO+ the most spatially widespread, and molecular emission seen toward 90% of the area above N(H_2) column densities of 1.9x10^21 cm^-2. HCO+ has the broadest velocity dispersion, near 0.3 km/s on average, compared to ~0.15 km/s for the other molecules, thus representing a range from supersonic to subsonic gas motions. Our non-binary dendrogram analysis reveals that the dense gas traced by each molecule has similar hierarchical structure, and that gas surrounding the candidate first hydrostatic core (FHSC), L1451-mm, and other previously detected single-dish continuum clumps have similar hierarchical structure; this suggests that different sub-regions of L1451 are fragmenting on the pathway to forming young stars. We determined the three-dimensional morphology of the largest detectable dense gas structures to be relatively ellipsoidal compared to other CLASSy regions, which appeared more flattened at largest scales. A virial analysis shows the most centrally condensed dust structures are likely unstable against collapse. Additionally, we identify a new spherical, centrally condensed N2H+ feature that could be a new FHSC candidate. The overall results suggest L1451 is a young region starting to form its generation of stars within turbulent, hierarchical structures.

The stages before the formation of stars in molecular clouds are poorly understood. Insights can be gained by studying the properties of quiescent clouds, such as their magnetic field structure. The plane-of-the-sky orientation of the field can be traced by polarized starlight. We present the first extended, wide-field ($\sim$10 $\rm deg^2$) map of the Polaris Flare cloud in dust-absorption induced optical polarization of background stars, using the RoboPol polarimeter at the Skinakas Observatory. This is the first application of the wide-field imaging capabilities of RoboPol. The data were taken in the R-band and analysed with the automated reduction pipeline of the instrument. We present in detail optimizations in the reduction pipeline specific to wide-field observations. Our analysis resulted in reliable measurements of 641 stars with median fractional linear polarization 1.3%. The projected magnetic field shows a large scale ordered pattern. At high longitudes it appears to align with faint striations seen in the Herschel-SPIRE map of dust emission (250 $\mu m$), while in the central 4-5 deg$^2$ it shows an eddy-like feature. The overall polarization pattern we obtain is in good agreement with large scale measurements by Planck of the dust emission polarization in the same area of the sky.

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.

We present a new method for assessing the intrinsic 3D shape of prestellar cores from molecular column densities. We have employed hydrodynamic simulations of contracting, isothermal cores considering three intrinsic geometries: spherical, cylindrical/filamentary and disk-like. We have coupled our hydrodynamic simulations with non-equilibrium chemistry. We find that a) when cores are observed very elongated (i.e. for aspect ratios $\le$ 0.15) the intrinsic 3D geometry can be probed by their 2D molecular emission maps, since these exhibit significant qualitative morphological differences between cylindrical and disk-like cores. Specifically, if a disk-like core is observed as a filamentary object in dust emission, then it will be observed as two parallel filaments in $\rm{N_2H^{+}}$; b) for cores with higher aspect ratios (i.e. 0.15 $\sim$ 0.9) we define a metric $\Delta$ that quantifies whether a molecular column density profile is centrally peaked, depressed or flat. We have identified one molecule ($\rm{CN}$) for which $\Delta$ as a function of the aspect ratio probes the 3D geometry of the core; and c) for cores with almost circular projections (i.e. for aspect ratios $\sim$ 1), we have identified three molecules ($\rm{OH}$, $\rm{CO}$ and $\rm{H_2CO}$) that can be used to probe the intrinsic 3D shape by close inspection of their molecular column density radial profiles. We alter the temperature and the cosmic-ray ionization rate and demonstrate that our method is robust against the choice of parameters.

We have carried out a search for filamentary structures in the Taurus molecular cloud using $\rm^{13}CO$ line emission data from the FCRAO survey of $\rm \sim100 \, deg^2$. We have used the topological analysis tool, DisPerSe, and post-processed its results to include a more strict definition of filaments that requires an aspect ratio of at least 3:1 and cross section intensity profiles peaked on the spine of the filament. In the velocity-integrated intensity map only 10 of the hundreds of filamentary structures identified by DisPerSe comply with our criteria. Unlike Herschel analyses, which find a characteristic width for filaments of $\rm \sim0.1 \, pc$, we find a much broader distribution of profile widths in our structures, with a peak at 0.4 pc. Furthermore, even if the identified filaments are cylindrical objects, their complicated velocity structure and velocity dispersions imply that they are probably gravitationally unbound. Analysis of velocity channel maps reveals the existence of hundreds of `velocity-coherent' filaments. The distribution of their widths is peaked at lower values (0.2 pc) while the fluctuation of their peak intensities is indicative of stochastic origin. These filaments are suppressed in the integrated intensity map due to the blending of diffuse emission from different velocities. Conversely, integration over velocities can cause filamentary structures to appear. Such apparent filaments can also be traced, using the same methodology, in simple simulated maps consisting of randomly placed cores. They have profile shapes similar to observed filaments and contain most of the simulated cores.

We present measurements of rotations of the optical polarization of blazars during the second year of operation of RoboPol, a monitoring programme of an unbiased sample of gamma-ray bright blazars specially designed for effective detection of such events, and we analyse the large set of rotation events discovered in two years of observation. We investigate patterns of variability in the polarization parameters and total flux density during the rotation events and compare them to the behaviour in a non-rotating state. We have searched for possible correlations between average parameters of the polarization-plane rotations and average parameters of polarization, with the following results: (1) there is no statistical association of the rotations with contemporaneous optical flares; (2) the average fractional polarization during the rotations tends to be lower than that in a non-rotating state; (3) the average fractional polarization during rotations is correlated with the rotation rate of the polarization plane in the jet rest frame; (4) it is likely that distributions of amplitudes and durations of the rotations have physical upper bounds, so arbitrarily long rotations are not realised in nature.

Recent Planck results have shown that radiation from the cosmic microwave background passes through foregrounds in which aligned dust grains produce polarized dust emission, even in regions of the sky with the lowest level of dust emission. One of the most commonly used ways to remove the dust foreground is to extrapolate the polarized dust emission signal from frequencies where it dominates (e.g., ~ 350 GHz) to frequencies commonly targeted by cosmic microwave background experiments (e.g., ~150 GHz). In this paper, we describe an interstellar medium effect that can lead to decorrelation of the dust emission polarization pattern between different frequencies due to multiple contributions along the line of sight. Using a simple 2-cloud model we show that there are two conditions under which this decorrelation can be large: (a) the ratio of polarized intensities between the two clouds changes between the two frequencies; (b) the magnetic fields between the two clouds contributing along a line of sight are significantly misaligned. In such cases, the 350 GHz polarized sky map is not predictive of that at 150 GHz. We propose a possible correction for this effect, using information from optopolarimetric surveys of dichroicly absorbed starlight.

We revisit the relation between magnetic-field strength ($B$) and gas density ($\rho$) for contracting interstellar clouds and fragments (or, cores), which is central in observationally determining the dynamical importance of magnetic fields in cloud evolution and star formation. Recently, it has been claimed that a relation $B \propto \rho^{2/3} $ is statistically preferred over $B \propto \rho^{1/2}$ in molecular clouds, when magnetic field detections and nondetections from Zeeman observations are combined. This finding has unique observational implications on cloud and core geometry: The relation $B \propto \rho^{2/3} $ can only be realized under spherical contraction. However, no indication of spherical geometry can be found for the objects used in the original statistical analysis of the $B-\rho$ relation. We trace the origin of the inconsistency to simplifying assumptions in the statistical model used to arrive at the $B\propto \rho^{2/3}$ conclusion and to an underestimate of observational uncertainties in the determination of cloud and core densities. We show that, when these restrictive assumptions are relaxed, $B \propto \rho^{1/2}$ is the preferred relation for the (self-gravitating) molecular-cloud data, as theoretically predicted four decades ago.

We present observations of N2H+(1-0), HCO+(1-0), and HCN(1-0) toward the Serpens Main molecular cloud from the CARMA Large Area Star Formation Survey (CLASSy). We mapped 150 square arcminutes of Serpens Main with an angular resolution of 7 arcsecs. The gas emission is concentrated in two subclusters (the NW and SE subclusters). The SE subcluster has more prominent filamentary structures and more complicated kinematics compared to the NW subcluster. The majority of gas in the two subclusters has subsonic to sonic velocity dispersions. We applied a dendrogram technique with N2H+(1-0) to study the gas structures; the SE subcluster has a higher degree of hierarchy than the NW subcluster. Combining the dendrogram and line fitting analyses reveals two distinct relations: a flat relation between nonthermal velocity dispersion and size, and a positive correlation between variation in velocity centroids and size. The two relations imply a characteristic depth of 0.15 pc for the cloud. Furthermore, we have identified six filaments in the SE subcluster. These filaments have lengths of 0.2 pc and widths of 0.03 pc, which is smaller than a characteristic width of 0.1 pc suggested by Herschel observations. The filaments can be classified into two types based on their properties. The first type, located in the northeast of the SE subcluster, has larger velocity gradients, smaller masses, and nearly critical mass-per-unit-length ratios. The other type, located in the southwest of the SE subcluster, has the opposite properties. Several YSOs are formed along two filaments which have supercritical mass per unit length ratios, while filaments with nearly critical mass-per-unit-length ratios are not associated with YSOs, suggesting that stars are formed on gravitationally unstable filaments.

Polarization in classical Be stars results from Thomson scattering of the unpolarized light from the Be star in the circumstellar disc. Theory and observations agree that the maximum degree of polarization from isolated Be stars is < 4%. We report on the first optical polarimetric observations of the Be/X-ray binary EXO\,2030+375. We find that the optical (R band) light is strongly linearly polarized with a degee of polarization of 19%, the highest ever measured either in a classical or Be/X-ray binary. We argue that the interstellar medium cannot account for this high polarization degree and that a substantial amount must be intrinsic to the source. We propose that it may result from the alignment of non-spherical ferromagnetic grains in the Be star disc due to the strong neutron star magnetic field.

We report the polarized optical light curve of a gamma-ray burst afterglow obtained using the RoboPol instrument. Observations began 655 seconds after the initial burst of gamma-rays from GRB131030A, and continued uninterrupted for 2 hours. The afterglow displayed a low, constant fractional linear polarization of $p = (2.1 \pm 1.6)\,\%$ throughout, which is similar to the interstellar polarization measured on nearby stars. The optical brightness decay is consistent with a forward-shock propagating in a medium of constant density, and the low polarization fraction indicates a disordered magnetic field in the shock front. This supports the idea that the magnetic field is amplified by plasma instabilities on the shock front. These plasma instabilities produce strong magnetic fields with random directions on scales much smaller than the total observable region of the shock, and the resulting randomly-oriented polarization vectors sum to produce a low net polarization over the total observable region of the shock.

We present details of the CARMA Large Area Star Formation Survey (CLASSy), while focusing on observations of Barnard 1. CLASSy is a CARMA Key Project that spectrally imaged N2H+, HCO+, and HCN (J=1-0 transitions) across over 800 square arcminutes of the Perseus and Serpens Molecular Clouds. The observations have angular resolution near 7" and spectral resolution near 0.16 km/s. We imaged ~150 square arcminutes of Barnard 1, focusing on the main core, and the B1 Ridge and clumps to its southwest. N2H+ shows the strongest emission, with morphology similar to cool dust in the region, while HCO+ and HCN trace several molecular outflows from a collection of protostars in the main core. We identify a range of kinematic complexity, with N2H+ velocity dispersions ranging from ~0.05-0.50 km/s across the field. Simultaneous continuum mapping at 3 mm reveals six compact object detections, three of which are new detections. A new non-binary dendrogram algorithm is used to analyze dense gas structures in the N2H+ position-position-velocity (PPV) cube. The projected sizes of dendrogram-identified structures range from about 0.01-0.34 pc. Size-linewidth relations using those structures show that non-thermal line-of-sight velocity dispersion varies weakly with projected size, while rms variation in the centroid velocity rises steeply with projected size. Comparing these relations, we propose that all dense gas structures in Barnard 1 have comparable depths into the sky, around 0.1-0.2 pc; this suggests that over-dense, parsec-scale regions within molecular clouds are better described as flattened structures rather than spherical collections of gas. Science-ready PPV cubes for Barnard 1 molecular emission are available for download.

The ratio of mass and magnetic flux determines the relative importance of magnetic and gravitational forces in the evolution of molecular clouds and their cores. Its measurement is thus central in discriminating between different theories of core formation and evolution. Here we discuss the effect of chemical depletion on measurements of the mass-to-flux ratio using the same molecule (OH) both for Zeeman measurements of the magnetic field and the determination of the mass of the region. The uncertainties entering through the OH abundance in determining separately the magnetic field and the mass of a region have been recognized in the literature. It has been proposed however that, when comparing two regions of the same cloud, the abundance will in both cases be the same. We show that this assumption is invalid. We demonstrate that when comparing regions with different densities, the effect of OH depletion in measuring changes of the mass-to-flux ratio between different parts of the same cloud can even reverse the direction of the underlying trends (for example, the mass-to-flux ratio may appear to decrease as we move to higher density regions). The systematic errors enter primarily through the inadequate estimation of the mass of the region.

We present the N2H+(J=1-0) map of the Serpens South molecular cloud obtained as part of the CARMA Large Area Star Formation Survey (CLASSy). The observations cover 250 square arcminutes and fully sample structures from 3000 AU to 3 pc with a velocity resolution of 0.16 km/s, and they can be used to constrain the origin and evolution of molecular cloud filaments. The spatial distribution of the N2H+ emission is characterized by long filaments that resemble those observed in the dust continuum emission by Herschel. However, the gas filaments are typically narrower such that, in some cases, two or three quasi-parallel N2H+ filaments comprise a single observed dust continuum filament. The difference between the dust and gas filament widths casts doubt on Herschel ability to resolve the Serpens South filaments. Some molecular filaments show velocity gradients along their major axis, and two are characterized by a steep velocity gradient in the direction perpendicular to the filament axis. The observed velocity gradient along one of these filaments was previously postulated as evidence for mass infall toward the central cluster, but these kind of gradients can be interpreted as projection of large-scale turbulence.

We present first results from RoboPol, a novel-design optical polarimeter operating at the Skinakas Observatory in Crete. The data, taken during the May - June 2013 commissioning of the instrument, constitute a single-epoch linear polarization survey of a sample of gamma-ray - loud blazars, defined according to unbiased and objective selection criteria, easily reproducible in simulations, as well as a comparison sample of, otherwise similar, gamma-ray - quiet blazars. As such, the results of this survey are appropriate for both phenomenological population studies and for tests of theoretical population models. We have measured polarization fractions as low as $0.015$ down to $R$ magnitude of 17 and as low as $0.035$ down to 18 magnitude. The hypothesis that the polarization fractions of gamma-ray - loud and gamma-ray - quiet blazars are drawn from the same distribution is rejected at the $10^{-3}$ level. We therefore conclude that gamma-ray - loud and gamma-ray - quiet sources have different optical polarization properties. This is the first time this statistical difference is demonstrated in optical wavelengths. The polarization fraction distributions of both samples are well-described by exponential distributions with averages of $\langle p \rangle =6.4 ^{+0.9}_{-0.8}\times 10^{-2}$ for gamma-ray--loud blazars, and $\langle p \rangle =3.2 ^{+2.0}_{-1.1}\times 10^{-2}$ for gamma-ray--quiet blazars. The most probable value for the difference of the means is $3.4^{+1.5}_{-2.0}\times 10^{-2}$. The distribution of polarization angles is statistically consistent with being uniform.

We describe the data reduction pipeline and control system for the RoboPol project. The RoboPol project is monitoring the optical $R$-band magnitude and linear polarization of a large sample of active galactic nuclei that is dominated by blazars. The pipeline calibrates and reduces each exposure frame, producing a measurement of the magnitude and linear polarization of every source in the $13'\times 13'$ field of view. The control system combines a dynamic scheduler, real-time data reduction, and telescope automation to allow high-efficiency unassisted observations.

Comparison of linewidths of spectral line profiles of ions and neutral molecules have been recently used to estimate the strength of the magnetic field in turbulent star-forming regions. However, the ion (HCO+) and neutral (HCN) species used in such studies may not be necessarily co-evolving at every scale and density and may thus not trace the same regions. Here, we use coupled chemical/dynamical models of evolving prestellar molecular cloud cores including non-equilibrium chemistry, with and without magnetic fields, to study the spatial distribution of HCO+ and HCN, which have been used in observations of spectral linewidth differences to date. In addition, we seek new ion-neutral pairs that are good candidates for such observations because they have similar evolution and are approximately co-spatial in our models. We identify three such good candidate pairs: HCO+/NO, HCO+/CO, and NO+/NO.

We study the effect that non-equilibrium chemistry in dynamical models of collapsing molecular cloud cores has on measurements of the magnetic field in these cores, the degree of ionization, and the mean molecular weight of ions. We find that OH and CN, usually used in Zeeman observations of the line-of-sight magnetic field, have an abundance that decreases toward the center of the core much faster than the density increases. As a result, Zeeman observations tend to sample the outer layers of the core and consistently underestimate the core magnetic field. The degree of ionization follows a complicated dependence on the number density at central densities up to 10^5 cm^{-3} for magnetic models and 10^6 cm^{-3} in non-magnetic models. At higher central densities the scaling approaches a power-law with a slope of -0.6 and a normalization which depends on the cosmic-ray ionization rate {\zeta} and the temperature T as ({\zeta}T)^1/2. The mean molecular weight of ions is systematically lower than the usually assumed value of 20 - 30, and, at high densities, approaches a value of 3 due to the asymptotic dominance of the H3+ ion. This significantly lower value implies that ambipolar diffusion operates faster.

We combine dynamical and non-equilibrium chemical modeling of evolving prestellar molecular cloud cores, and explore the evolution of molecular abundances in the contracting core. We model both magnetic cores, with varying degrees of initial magnetic support, and non-magnetic cores, with varying collapse delay times. We explore, through a parameter study, the competing effects of various model parameters in the evolving molecular abundances, including the elemental C/O ratio, the temperature, and the cosmic-ray ionization rate. We find that different models show their largest quantitative differences at the center of the core, whereas the outer layers, which evolve slower, have abundances which are severely degenerate among different dynamical models. There is a large range of possible abundance values for different models at a fixed evolutionary stage (central density), which demonstrates the large potential of chemical differentiation in prestellar cores. However, degeneracies among different models, compounded with uncertainties induced by other model parameters, make it difficult to discriminate among dynamical models. To address these difficulties, we identify abundance ratios between particular molecules, the measurement of which would have maximal potential for discrimination among the different models examined here. In particular, we find that the ratios between NH3 and CO; NH2 and CO; NH3 and HCO+ are sensitive to the evolutionary timescale, and that the ratio between HCN and OH is sensitive to the C/O ratio. Finally, we demonstrate that measurements of the central deviation (central depletion or enhancement) of abundances of certain molecules are good indicators of the dynamics of the core.

The tight relation of star formation with molecular gas indicated by observations and assumed in recent models implies that the efficiency with which galaxies convert their gas into stars depends on gas metallicity. This is because the abundance of molecular hydrogen is sensitive to the abundance of dust, which catalyzes the formation of H_2 and helps to shield it from dissociating radiation. In this study we point out that in the absence of significant pre-enrichment by Population III stars forming out of zero metallicity gas, such H_2-based star formation is expected to leave an imprint in the form of bi-modality in the metallicity distribution among dwarf galaxies and in the metallicity distribution of stars within individual galaxies. The bi-modality arises because when gas metallicity (and dust abundance) is low, formation of molecular gas is inefficient, the gas consumption time scale is long, and star formation and metal enrichment proceed slowly. When metallicity reaches a critical threshold value star formation and enrichment accelerate, which leads to rapid increase in both stellar mass and metallicity of galaxies. We demonstrate this process both using a simple analytical model and full cosmological simulations. In contrast, observed metallicity distributions of dwarf galaxies or stars within them are not bi-modal. We argue that this discrepancy points to substantial early stochastic pre-enrichment by population III stars to levels Z ~ 0.01 Z_sun in dense, star forming regions of early galaxies.

We propose a simple analytical method for estimating the central volume density of prestellar molecular cloud cores from their column density profiles. Prestellar cores feature a flat central part of the column density and volume density profiles of the same size indicating the existence of a uniform density inner region. The size of this region is set by the thermal pressure force which depends only on the central volume density and temperature of the core, and can provide a direct measurement of the central volume density. Thus a simple length measurement can immediately yield a central density estimate independent of any dynamical model for the core and without the need for fitting. Using the radius at which the column density is 90% of the central value as an estimate of the size of the flat inner part of the column density profile yields an estimate of the central volume density within a factor of 2 for well resolved cores.

We investigate the effect of the Milky Way's magnetic field in star forming regions using archived 350 micron polarization data on 52 Galactic star formation regions from the Hertz polarimeter module. The polarization angles and percentages for individual telescope beams were combined in order to produce a large-scale average for each source and for complexes of sources. In more than 80% of the sources, we find a meaningful mean magnetic field direction, implying the existence of an ordered magnetic field component at the scale of these sources. The average polarization angles were analyzed with respect to the Galactic coordinates in order to test for correlations between polarization percentage, polarization angle, intensity, and Galactic location. No correlation was found, which suggests that the magnetic field in dense molecular clouds is decoupled from the large-scale Galactic magnetic field. Finally, we show that the magnetic field directions in the complexes are consistent with a random distribution on the sky.

We present results from numerical simulations of the cooling-core cluster A2199 produced by the two-dimensional (2-D) resistive magnetohydrodynamics (MHD) code MACH2. In our simulations we explore the effect of anisotropic thermal conduction on the energy balance of the system. The results from idealized cases in 2-D axisymmetric geometry underscore the importance of the initial plasma density in ICM simulations, especially the near-core values since the radiation cooling rate is proportional to ${n_e}^2$. Heat conduction is found to be non-effective in preventing catastrophic cooling in this cluster. In addition we performed 2-D planar MHD simulations starting from initial conditions deliberately violating both thermal balance and hydrostatic equilibrium in the ICM, to assess contributions of the convective terms in the energy balance of the system against anisotropic thermal conduction. We find that in this case work done by the pressure on the plasma can dominate the early evolution of the internal energy over anisotropic thermal conduction in the presence of subsonic flows, thereby reducing the impact of the magnetic field. Deviations from hydrostatic equilibrium near the cluster core may be associated with transient activity of a central active galactic nucleus and/or remnant dynamical activity in the ICM and warrant further study in three dimensions.

We had recently re-analyzed in a self-consistent way OH-Zeeman observations in four molecular-cloud envelopes and we had shown that, contrary to claims by Crutcher et al., there is no evidence that the mass-to-flux ratio decreases from the envelopes to the cores of these clouds. The key difference between our data analysis and the earlier one by Crutcher et al. is the relaxation of the overly restrictive assumption made by Crutcher et al, that the magnetic field strength is independent of position in each of the four envelopes. In a more recent paper, Crutcher et al. (1) claim that our analysis is not self-consistent, in that it misses a cosine factor, and (2) present new arguments to support their contention that the magnetic-field strength is indeed independent of position in each of the four envelopes. We show that the claim of the missing cosine factor is false, that the new arguments contain even more serious problems than the Crutcher et al. original data analysis, and we present new observational evidence, independent of the OH-Zeeman data, that suggests significant variations in the magnetic-field strength in the four cloud envelopes.

Recent observations of column densities in molecular clouds find lognormal distributions with power-law high-density tails. These results are often interpreted as indications that supersonic turbulence dominates the dynamics of the observed clouds. We calculate and present the column-density distributions of three clouds, modeled with very different techniques, none of which is dominated by supersonic turbulence. The first star-forming cloud is simulated using smoothed particle hydrodynamics (SPH); in this case gravity, opposed only by thermal-pressure forces, drives the evolution. The second cloud is magnetically subcritical with subsonic turbulence, simulated using nonideal MHD; in this case the evolution is due to gravitationally-driven ambipolar diffusion. The third cloud is isothermal, self-gravitating, and has a smooth density distribution analytically approximated with a uniform inner region and an r^-2 profile at larger radii. We show that in all three cases the column-density distributions are lognormal. Power-law tails develop only at late times (or, in the case of the smooth analytic profile, for strongly centrally concentrated configurations), when gravity dominates all opposing forces. It therefore follows that lognormal column-density distributions are generic features of diverse model clouds, and should not be interpreted as being a consequence of supersonic turbulence.

The ambipolar-diffusion theory of star formation predicts the formation of fragments in molecular clouds with mass-to-flux ratios greater than that of the parent-cloud envelope. By contrast, scenarios of turbulence-induced fragmentation do not yield such a robust prediction. Based on this property, Crutcher et al. (2009) proposed an observational test that could potentially discriminate between fragmentation theories. However, the analysis applied to the data severely restricts the discriminative power of the test: the authors conclude that they can only constrain what they refer to as the "idealized" ambipolar-diffusion theory that assumes initially straight-parallel magnetic field lines in the parent cloud. We present an original, self-consistent analysis of the same data taking into account the nonuniformity of the magnetic field in the cloud envelopes, which is suggested by the data themselves, and we discuss important geometrical effects that must be accounted for in using this test. We show quantitatively that the quality of current data does not allow for a strong conclusion about any fragmentation theory. Given the discriminative potential of the test, we urge for more and better-quality data.

We present a novel statistical analysis aimed at deriving the intrinsic shapes and magnetic field orientations of molecular clouds using dust emission and polarization observations by the Hertz polarimeter. Our observables are the aspect ratio of the projected plane-of-the-sky cloud image, and the angle between the mean direction of the plane-of-the-sky component of the magnetic field and the short axis of the cloud image. To overcome projection effects due to the unknown orientation of the line-of-sight, we combine observations from 24 clouds, assuming that line-of-sight orientations are random and all are equally probable. Through a weighted least-squares analysis, we find that the best-fit intrinsic cloud shape describing our sample is an oblate disk with only small degrees of triaxiality. The best-fit intrinsic magnetic field orientation is close to the direction of the shortest cloud axis, with small (~24 deg) deviations toward the long/middle cloud axes. However, due to the small number of observed clouds, the power of our analysis to reject alternative configurations is limited.

We describe a phenomenological model for molecular hydrogen formation suited for applications in galaxy formation simulations, which includes on-equilibrium formation of molecular hydrogen on dust and approximate treatment of both its self-shielding and shielding by dust from the dissociating UV radiation. The model is applicable in simulations in which individual star forming regions - the giant molecular complexes - can be identified (resolution of tens of pc) and their mean internal density estimated reliably, even if internal structure is not resolved. In agreement with previous studies, calculations based on our model show that the transition from atomic to fully molecular phase depends primarily on the metallicity, which we assume is directly related to the dust abundance, and clumpiness of the interstellar medium. The clumpiness simply boosts the formation rate of molecular hydrogen, while dust serves both as a catalyst of molecular hydrogen formation and as an additional shielding from dissociating UV radiation. The upshot is that it is difficult to form fully-shielded giant molecular clouds while gas metallicity is low. However, once the gas is enriched to Z ~ 0.01-0.1 solar, the subsequent star formation and enrichment can proceed at a much faster rate. This may keep star formation efficiency in the low-mass, low-metallicity progenitors of galaxies very low for a certain period of time with the effect similar to a strong "feedback" mechanism. [abridged]

Until recently, many of the dozens of quantitative predictions of the ambipolar-diffusion theory of gravitational fragmentation (or core formation) of molecular clouds have been confirmed by observations and, just as importantly, no prediction has been contradicted by any observation. A recent paper, however, claims that measurements of the variation of the mass-to-flux ratio from envelopes to cores in four clouds {\it de}creases, in direct contrast to a prediction of the theory but in agreement with turbulent fragmentation (in the absence of gravity) and, therefore, the ambipolar-diffusion theory is invalid (Crutcher et al 2008). The paper treats magnetic-field nondetections as if they were detections. We show that the analysis of the data is fundamentally flawed and, moreover, the comparison with the theoretical prediction ignores major geometrical effects, suggested by the data themselves if taken at face value. The magnetic fluxes of the envelopes are also miscalculated. We carry out a proper error analysis and treatment of the nondetections and we show that the claimed measurement of the variation of the mass-to-flux ratio from envelopes to cores is not valid, no contradiction with the ambipolar-diffusion theory can be concluded, and no theory can be tested on the basis of these data.

Using interferometric data from BIMA observations, combined with detailed modeling in Fourier space of the physical structures predicted by models, we constrain the circumstellar envelope parameters for four Class 0 young stellar objects, as well as their embedded circumstellar disks. The envelopes of these objects are still undergoing collapse, and theoretical collapse models can be compared to the observations. Since it has been suggested in a previous study that both the Larson-Penston and Shu similarity solutions underestimate the age of the system, we adopt Tassis & Mouschovias' model of the collapse process, which includes all relevant magnetic fields effects. The results of the model fitting show a good consistency between theory and data; furthermore, no age problem exists since the Tassis & Mouschovias' model is age independent for the first 255 kyr. Although the majority of the continuum dust emission arises from the circumstellar envelopes, these objects have well known outflows, which suggest the presence of circumstellar disks. At the highest resolution, most of the large-scale envelope emission is resolved out by interferometry, but the small-scale residual emission remains, making it difficult to observe only the compact disk component. By modeling the emission of the envelope and subtracting it from the total emission, we constrain the disk masses in our four systems to be comparable to or smaller than the typical disk masses for T Tauri systems.

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.

The surface density of the star formation rate in different galaxies, as well as in different parts of a single galaxy, scales nonlinearly with the surface density of the total gas. This observationally established relation is known as the Kennicutt-Schmidt star formation law. The slope of the star formation law has been shown to change with the density of the gas against which the star formation rate is plotted. This dependence implies a nonlinear scaling between the dense gas and the total gas surface densities within galaxies. Here, we explore a possible interpretation of this scaling as a property of the geometry of the interstellar medium (ISM), and we find that it arises naturally if the topology of the ISM is multifractal. Under the additional assumption that, at very high densities, the star formation timescale is roughly constant, the star formation law itself can also be recovered as a consequence of the multifractal geometry of the ISM. The slope of the scaling depends on the width of the global probability density function (PDF), and is between 1.5 and 1.6 for wide PDFs relevant to high-mass systems, while it is higher for narrower PDFs appropriate for lower-mass dwarf galaxies, in agreement with observations.

We investigate the intrinsic shapes of starless cores in the Orion GMC, using the prestellar core sample of Nutter and Ward-Thompson (2007), which is based on submillimeter SCUBA data. We employ a maximum-likelihood method to reconstruct the intrinsic distribution of ellipsoid axial ratios from the axial ratios of projected plane-of-the-sky core ellipses. We find that, independently of the details of the assumed functional form of the distribution, there is a strong preference for oblate cores of finite thickness. Cores with varying finite degrees of triaxiality are a better fit than purely axisymmetric cores although cores close to axisymmetry are not excluded by the data. The incidence of prolate starless cores in Orion is found to be very infrequent. We also test the consistency of the observed data with a uniform distribution of intrinsic shapes, which is similar to those found in gravoturbulent fragmentation simulations. This distribution is excluded at the 0.1% level. These findings have important implications for theories of core formation within molecular clouds.

We formulate the problem of the formation of magnetically supercritical cores in magnetically subcritical parent molecular clouds, and the subsequent collapse of the cores to high densities, past the detachment of ions from magnetic field lines and into the opaque regime. We employ the six-fluid MHD equations, accounting for the effects of grains (negative, positive and neutral) including their inelastic collisions with other species. We do not assume that the magnetic flux is frozen in any of the charged species. We derive a generalized Ohm's law that explicitly distinguishes between flux advection (and the associated process of ambipolar diffusion) and Ohmic dissipation, in order to assess the contribution of each mechanism to the increase of the mass-to-flux ratio of the central parts of a collapsing core and possibly to the resolution of the magnetic flux problem of star formation. We show how our formulation is related to and can be transformed into the traditional, directional formulation of the generalized Ohm's law, and we derive formulae for the perpendicular, parallel and Hall conductivities entering the latter, which include, for the first time, the effect of inelastic collisions between grains. In addition, we present a general (valid in any geometry) solution for the velocities of charged species as functions of the velocity of the neutrals and of the effective flux velocity (which can in turn be calculated from the dynamics of the system and Faraday's law). The last two sets of formulae can be adapted for use in any general non-ideal MHD code to study phenomena beyond star formation in magnetic clouds. The results, including a detailed parameter study, are presented in two accompanying papers.

In two previous papers we formulated and solved, for a fiducial set of free parameters, the problem of the formation and evolution of a magnetically supercritical core inside a magnetically subcritical parent cloud. In this paper we present a parameter study to assess the sensitivity of the results (1) to the density at which the equation of state becomes adiabatic; (2) to the initial mass-to-flux ratio of the parent cloud; and (3) to ionization by radioactive decay of different nuclei (40K and 26Al) at high densities (number density > 10^12 particles per cubic cm). We find that (1) the results depend only slightly on the density at which the onset of adiabaticity occurs; (2) memory of the initial mass-to-flux ratio is completely lost at late times, which emphasizes the relevance of this work, idependently of the adopted theory of core formation; and (3) the precise source of radioactive ionization alters the degree of attachment of the electrons to the field lines (at high densities), and the relative importance of ambipolar diffusion and Ohmic dissipation in reducing the magnetic flux of the protostar. The value of the magnetic field at the end of the runs is insensitive to the values of the free parameters and in excellent agreement with meteoritic measurements of the protosolar nebula magnetic field. The magnetic flux problem of star formation is resolved for at least strongly magnetic newborn stars. A complete detachment of the magnetic field from the matter is unlikely. The formation of a "magnetic wall" (with an associated magnetic shock) is independent of the assumed equation of state, although the process is enhanced and accelerated by the formation of a central hydrostatic core.

We follow the ambipolar-diffusion--driven formation and evolution of a fragment in a magnetically supported molecular cloud, until a hydrostatic protostellar core forms at its center. This problem was formulated in Paper I. We determine the density, velocity and magnetic field as functions of space and time, and the contribution of ambipolar diffusion and Ohmic dissipation to the resolution of the magnetic flux problem of star formation. The issue of whether the magnetic field ever decouples from the (neutral) matter is also addressed. We also find that the electrons do not decouple from the field lines before thermal ionization becomes important and recouples the magnetic field to the neutral matter. Ohmic dissipation becomes more effective than ambipolar diffusion as a flux reduction mechanism only at the highest densities (a few times 10^12 particles per cubic cm). In the high-density central parts of the core, the magnetic field acquires an almost spatially uniform structure, with a value that, at the end of the calculation (number density ~ 5 times 10^14 particles per cubic cm), is found to be in excellent agreement with meteoritic measurements of magnetic fields in the protosolar nebula. Outside the hydrostatic protostellar core, a concentration of magnetic flux (a "magnetic wall") forms, which gives rise to a magnetic shock. This magnetic shock is the precursor of the repeated shocks found by Tassis & Mouschovias (2005) which cause spasmodic accretion onto the hydrostatic core at later times.

We revisit the problem of the star formation timescale and the ages of molecular clouds. The apparent overabundance of star-forming molecular clouds over clouds without active star formation has been thought to indicate that molecular clouds are "short-lived" and that star formation is "rapid". We show that this statistical argument lacks self-consistency and, even within the rapid star-formation scenario, implies cloud lifetimes of approximately 10 Myr. We discuss additional observational evidence from external galaxies that indicate lifetimes of molecular clouds and a timescale of star formation of approximately 10 Myr . These long cloud lifetimes in conjunction with the rapid (approximately 1 Myr) decay of supersonic turbulence present severe difficulties for the scenario of turbulence-controlled star formation. By contrast, we show that all 31 existing observations of objects for which the linewidth, the size, and the magnetic field strength have been reliably measured are in excellent quantitative agreement with the predictions of the ambipolar-diffusion theory. Within the ambipolar-diffusion-controlled star formation theory the linewidths may be attributed to large-scale non-radial cloud oscillations (essentially standing large-amplitude, long-wavelength Alfven waves), and the predicted relation between the linewidth, the size, and the magnetic field is a natural consequence of magnetic support of self-gravitating clouds.

We formulate the problem of the late accretion phase of the evolution of an isothermal magnetic disk surrounding a forming star. The evolution is described by the six-fluid MHD equations, accounting for the presence of neutrals, atomic and molecular ions, electrons, and neutral, positively, and negatively charged grains. Only the electron fluid is assumed to be attached to the magnetic field, in order to investigate the effect of the detachment of the ions from the magnetic field lines that begins at densities as low as 10^8 cm^-3. 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. In this way, the structure and evolution of the isothermal disk surrounding the forming star can be studied at late times without having to implement the numerically costly radiative transfer required by the physics of the opaque core. The mass and magnetic flux accumulating in the forming star arecalculated, as are their effects on the structure & evolution of the surrounding disk. The numerical method of solution first uses an adaptive grid and later, after a central region a few AU in radius becomes opaque, switches to a stationary but nonuniform grid with a central sink cell. It also involves an implicit time integrator, an advective difference scheme that possesses the transportive property, a second-order difference approximation of forces inside a cell, an integral approximation of the gravitational and magnetic fields, and tensor artificial viscosity that permits an accurate investigation of the formation and evolution of shocks in the neutral fluid.

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.

We re-examine critically the estimates of the duration of different phases of star formation and the lifetimes of molecular clouds, based on the ages of T-Tauri stars, age spreads of stars in clusters, and statistics of pre-stellar cores. We show that all available observational data are consistent with lifetimes of molecular clouds comparable to 10 Myr, as well as with the predictions of the theory of self-initiated, ambipolar-diffusion--controlled star formation. We conclude that there exists no observational support for either "young" molecular clouds or "rapid" star formation.

We present first results from three-dimensional hydrodynamic simulations of the high redshift formation of dwarf galaxies. The simulations use an Eulerian adaptive mesh refinement technique to follow the non-equilibrium chemistry of hydrogen and helium with cosmological initial conditions drawn from a popular Lambda-dominated CDM model. We include the effects of reionization using a uniform radiation field, a phenomenological description of the effect of star formation and, in a separate simulation, the effects of stellar feedback. The results highlight the effects of stellar feedback and photoionization on the baryon content and star formation of galaxies with virial temperatures of approximately 10^4K. Dwarf sized dark matter halos that assemble prior to reionization are able to form stars. Most halos of similar mass that assemble after reionization do not form stars by redshift of three. Dwarf galaxies that form stars show large variations in their gas content because of stellar feedback and photoionization effects. Baryon-to-dark matter mass ratios are found to lie below the cosmic mean as a result of stellar feedback. The supposed substructure problem of CDM is critically assessed on the basis of these results. The star formation histories modulated by radiative and stellar feedbacks are discussed. In addition, metallicities of individual objects are shown to be naturally correlated with their mass-to-light ratios as is also evident in the properties of local dwarf galaxies.