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.

Herschel surveys of Galactic clouds support a paradigm for low-mass star formation in which dense filaments play a crucial role. The detailed fragmentation properties of star-forming filaments remain poorly understood, however, and the validity of the filament paradigm in the high-mass regime is still unclear. To investigate the density/velocity structure of the filament in the high-mass star-forming region NGC6334, we conducted ALMA observations in the 3mm continuum and the N2H+(1-0) line at ~3arcsec resolution. The filament was detected in both tracers. We identified 26 cores at 3mm and 5 velocity-coherent fiber-like features in N2H+ within the filament. The typical length of, and velocity difference between, the fiber-like features of the NGC6334 filament are reminiscent of the properties for the fibers of the low-mass star-forming filament B211/B213. Only 2 or 3 of the 5 velocity-coherent features are well aligned with the filament and may represent genuine, fiber sub-structures. The core mass distribution has a peak at ~10Msun. They can be divided into 7 groups of cores, closely associated with ArTeMiS clumps. The projected separation between cores and the projected spacing between clumps are roughly consistent with the effective Jeans length in the filament and a physical scale of about 4 times the filament width, respectively, suggesting a bimodal filament fragmentation process. Despite being one order of magnitude denser and more massive than the B211/B213 filament, the NGC6334 filament has a similar density/velocity structure. The difference is that the cores in NGC6334 appear to be an order of magnitude denser and more massive than the cores in Taurus. This suggests that dense filaments may evolve and fragment in a similar manner in low- and high-mass star-forming regions, and that the filament paradigm may hold in the intermediate-mass (if not high-mass) star formation regime.

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.

The connection between the pre-stellar core mass function (CMF) and the stellar initial mass function (IMF) lies at the heart of all star formation theories. In this paper, we study the earliest phases of star formation with a series of high-resolution numerical simulations that include the formation of sinks. In particular, we focus on the transition from cores to sinks within a massive molecular filament. We compare the CMF and IMF between magnetized and unmagnetized simulations, and between different resolutions. We find that selecting cores based on their kinematic virial parameter excludes collapsing objects because they host large velocity dispersions. Selecting only the thermally unstable magnetized cores, we observe that their mass-to-flux ratio spans almost two orders of magnitude for a given mass. We also see that, when magnetic fields are included, the CMF peaks at higher core mass values with respect to pure hydrodynamical simulations. Nonetheless, all models produce sink mass functions with a high-mass slope consistent with Salpeter. Finally, we examine the effects of resolution and find that, in isothermal simulations, even models with very high dynamical range fail to converge in the mass function. Our main conclusion is that, although the resulting CMFs and IMFs have similar slopes in all simulations, the cores have slightly different sizes and kinematical properties when a magnetic field is included. However, a core selection based on the mass-to-flux ratio alone is not enough to alter the shape of the CMF, if we do not take thermal stability into account. Finally, we conclude that extreme care should be given to resolution issues when studying sink formation with an isothermal equation of state.

Herschel observations have emphasized the role of molecular filaments in star formation. However, the origin and evolution of these filaments are not yet well understood, partly because of the lack of kinematic information. To examine whether the B211/B213 filament is accreting background gas due to its gravitational potential, we produced a toy accretion model and compared its predictions to the 12CO(1--0) and 13CO(1--0) velocity patterns. We also examined the spatial distributions of Halpha, 857 GHz continuum, and HI emission to search for evidence of large-scale external effects. We estimated the depth of the cloud around the B211/B213 filament to be 0.3--0.7 pc under the assumption that the density of the gas is the same as the 13CO critical density. Compared to a linear extent of >10 pc in the plane of the sky, this suggests that the 3D morphology of the cloud is sheet-like. 12CO and 13CO PV diagrams perpendicular to the filament axis show that the emission from the gas surrounding B211/B213 is redshifted to the northeast of the filament and blueshifted to the southwest, respectively, and that the velocities of both components approach the filament velocity as the line of sight approaches the filament crest. The PV diagrams predicted by our accretion model are in good agreement with the observed 12CO and 13CO PV diagrams, supporting the scenario of mass accretion into the filament proposed by Palmeirim et al. Moreover, inspection of the distribution of the Halpha and 857 GHz emission in the Taurus-California-Perseus region suggests that the B211/B213 filament may have formed as a result of an expanding supershell generated by the Per OB2 association. Based on these results, we propose a scenario in which the B211/B213 filament was initially formed by large-scale compression of HI gas and then is now growing in mass due to the gravitational accretion of ambient cloud molecular gas.

Magnetic fields play a very important role in the evolution of galaxies through their direct impact on star formation and stellar feedback-induced turbulence. However, their co-evolution with these processes has still not been thoroughly investigated, and the possible effect of the initial conditions is largely unknown. This letter presents the first results from a series of high-resolution numerical models, aimed at deciphering the effect of the initial conditions and of stellar feedback on the evolution of the galactic magnetic field in isolated, Milky-Way-like galaxies. The models start with an ordered, either poloidal or toroidal, magnetic field of varying strength, and are evolved with and without supernova feedback. They include a dark matter halo, a stellar and a gaseous disk, as well as the appropriate cooling and heating processes for the interstellar medium. Independently of the initial conditions, the galaxies develop a turbulent velocity field and a random magnetic field component in under 15 Myrs. Supernova feedback is extremely efficient in building a random magnetic field component up to large galactic heights. However, a random magnetic field emerges even in runs without feedback, which points to an inherent instability of the ordered component. Supernova feedback greatly affects the velocity field of the galaxy up to large galactic heights, and helps restructure the magnetic field up to 10 kpc above the disk, independently of the initial magnetic field morphology. On the other hand, the initial morphology of the magnetic field can accelerate the development of a random component at large heights. These effects have important implications for the study of the magnetic field evolution in galaxy simulations.

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.

Large-scale shocks formed by clustered feedback of young OB stars are considered an important source of mechanical energy for the ISM and a trigger of molecular cloud formation. Their interaction sites are locations where kinetic energy and magnetic fields are redistributed between ISM phases. In this work we study the effect of the magnetic field on the expansion and fragmentation of supershells and look for the signatures of supershell collisions on dense structures and on the kinetic and magnetic energy distribution of the ISM. We performed a series of high-resolution, three-dimensional simulations of colliding supershells. The shocks are created by time-dependent feedback and evolve in a diffuse turbulent environment that is either unmagnetized or has different initial magnetic field configurations. In the hydrodynamical situation, the expansion law of the superbubbles is consistent with the radius-time relation that is theoretically predicted for wind-blown bubbles. The supershells fragment over their entire surface into small dense clumps that carry more than half of the total kinetic energy in the volume. However, this is not the case when a magnetic field is introduced, either in the direction of the collision or perpendicular to the collision. In magnetized situations, the shell surfaces are more stable to dynamical instabilities. When the magnetic field opposes the collision, the expansion law of the supershells also becomes significantly flatter than in the hydrodynamical case. Although a two-phase medium arises in all cases, in the MHD simulations the cold phase is limited to lower densities.

The filamentary structure of the molecular interstellar medium and the potential link of this morphology to star formation have been brought into focus recently by high resolution observational surveys. An especially puzzling matter is that local interstellar filaments appear to have the same thickness, independent of their column density. This requires a theoretical understanding of their formation process and the physics that governs their evolution. In this work we explore a scenario in which filaments are dissipative structures of the large-scale interstellar turbulence cascade and ion-neutral friction (also called ambipolar diffusion) is affecting their sizes by preventing small-scale compressions. We employ high-resolution, 3D MHD simulations, performed with the grid code RAMSES, to investigate non-ideal MHD turbulence as a filament formation mechanism. We focus the analysis on the mass and thickness distributions of the resulting filamentary structures. Simulations of both driven and decaying MHD turbulence show that the morphologies of the density and the magnetic field are different when ambipolar diffusion is included in the models. In particular, the densest structures are broader and more massive as an effect of ion-neutral friction and the power spectra of both the velocity and the density steepen at a smaller wavenumber. The comparison between ideal and non-ideal MHD simulations shows that ambipolar diffusion causes a shift of the filament thickness distribution towards higher values. However, none of the distributions exhibit the pronounced peak found in the observed local filaments. Limitations in dynamical range and the absence of self-gravity in these numerical experiments do not allow us to conclude at this time whether this is due to the different filament selection or due to the physics inherent of the filament formation.

Stellar winds and supernova (SN) explosions of massive stars ("stellar feedback") create bubbles in the interstellar medium (ISM) and insert newly produced heavy elements and kinetic energy into their surroundings, possibly driving turbulence. Most of this energy is thermalized and immediately removed from the ISM by radiative cooling. The rest is available for driving ISM dynamics. In this work we estimate the amount of feedback energy retained as kinetic energy when the bubble walls have decelerated to the sound speed of the ambient medium. We show that the feedback of the most massive star outweighs the feedback from less massive stars. For a giant molecular cloud (GMC) mass of 1e5 solar masses (as e.g. found in the Orion GMCs) and a star formation efficiency of 8% the initial mass function predicts a most massive star of approximately 60 solar masses. For this stellar evolution model we test the dependence of the retained kinetic energy of the cold GMC gas on the inclusion of stellar winds. In our model winds insert 2.34 times the energy of a SN and create stellar wind bubbles serving as pressure reservoirs. We find that during the pressure driven phases of the bubble evolution radiative losses peak near the contact discontinuity (CD), and thus, the retained energy depends critically on the scales of the mixing processes across the CD. Taking into account the winds of massive stars increases the amount of kinetic energy deposited in the cold ISM from 0.1% to a few percent of the feedback energy.

Understanding the formation of binary and multiple stellar systems largely comes down to studying the circumstances for the fragmentation of a condensing core during the first stages of the collapse. However, the probability of fragmentation and the number of fragments seem to be determined to a large degree by the initial conditions. In this work we study the fate of the linear perturbations of a homogeneous gas sphere both analytically and numerically. In particular, we investigate the stability of the well-known homologous solution that describes the collapse of a uniform spherical cloud. The difficulty of the mathematical singularity in the perturbation equations is surpassed here by explicitly introducing a weak shock next to the sonic point. In parallel, we perform adaptive mesh refinement (AMR) numerical simulations of the linear stages of the collapse and compared the growth rates obtained by each method. With this combination of analytical and numerical tools, we explore the behavior of both spherically symmetric and non-axisymmetric perturbations. The numerical experiments provide the linear growth rates as a function of the core's initial virial parameter and as a function of the azimuthal wave number of the perturbation. The overlapping regime of the numerical experiments and the analytical predictions is the situation of a cold and large cloud, and in this regime the analytically calculated growth rates agree very well with the ones obtained from the simulations. The use of a weak shock as part of the perturbation allows us to find a physically acceptable solution to the equations for a continuous range of growth rates. The numerical simulations agree very well with the analytical prediction for the most unstable cores, while they impose a limit of a virial parameter of 0.1 for core fragmentation in the absence of rotation.

Dense, star-forming gas is believed to form at the stagnation points of large-scale ISM flows, but observational examples of this process in action are rare. We here present a giant molecular cloud (GMC) sandwiched between two colliding Milky Way supershells, which we argue shows strong evidence of having formed from material accumulated at the collision zone. Combining 12CO, 13CO and C18O(J=1-0) data with new high-resolution, 3D hydrodynamical simulations of colliding supershells, we discuss the origin and nature of the GMC (G288.5+1.5), favoring a scenario in which the cloud was partially seeded by pre-existing denser material, but assembled into its current form by the action of the shells. This assembly includes the production of some new molecular gas. The GMC is well interpreted as non-self-gravitating, despite its high mass (MH2 ~ 1.7 x 10^5 Msol), and is likely pressure confined by the colliding flows, implying that self-gravity was not a necessary ingredient for its formation. Much of the molecular gas is relatively diffuse, and the cloud as a whole shows little evidence of star formation activity, supporting a scenario in which it is young and recently formed. Drip-like formations along its lower edge may be explained by fluid dynamical instabilities in the cooled gas.

This work aims at improving the current understanding of the interaction between H ii regions and turbulent molecular clouds. We propose a new method to determine the age of a large sample of OB associations by investigating the development of their associated H ii regions in the surrounding turbulent medium. Using analytical solutions, one-dimensional (1D), and three-dimensional (3D) simulations, we constrained the expansion of the ionized bubble depending on the turbulent level of the parent molecular cloud. A grid of 1D simulations was then computed in order to build isochrone curves for H ii regions in a pressure-size diagram. This grid of models allowed to date large sample of OB associations and was used on the H ii Region Discovery Survey (HRDS). Analytical solutions and numerical simulations showed that the expansion of H ii regions is slowed down by the turbulence up to the point where the pressure of the ionized gas is in a quasi-equilibrium with the turbulent ram pressure. Based on this result, we built a grid of 1D models of the expansion of H ii regions in a profile based on Larson laws. The 3D turbulence is taken into account by an effective 1D temperature profile. The ages estimated by the isochrones of this grid agree well with literature values of well-known regions such as Rosette, RCW 36, RCW 79, and M16. We thus propose that this method can be used to give ages of young OB associations through the Galaxy such as the HRDS survey and also in nearby extra-galactic sources.

The fate of metals ejected by young OB associations into the Interstellar Medium (ISM) is investigated numerically. In particular, we study the enrichment of the cold gas phase, which is the material that forms molecular clouds. Following previous work, the expansion and collision of two supershells in a diffuse ISM is simulated, in this case also introducing an advected quantity which represents the metals expelled by the young stars. We adopt the simplest possible approach, not differentiating between metals coming from stellar winds and those coming from supernovae. Even though the hot, diffuse phase of the ISM receives a significant amount of metals from the stars, the cold phase is efficiently shielded, with very little metal enrichment. Significant enrichment of the cold ISM will therefore be delayed by at least the cooling time of this hot phase. No variations in cloud metallicity with distance from the OB association or with direction are found, which means that the shell collision does little to enhance the metallicity of the cold clumps. We conclude that the stellar generation that forms out of molecular structures, triggered by shell collisions cannot be significantly enriched.

The expansion and collision of two wind-blown superbubbles is investigated numerically. Our models go beyond previous simulations of molecular cloud formation from converging gas flows by exploring this process with realistic flow parameters, sizes and timescales. The superbubbles are blown by time-dependent winds and supernova explosions, calculated from population synthesis models. They expand into a uniform or turbulent diffuse medium. We find that dense, cold gas clumps and filaments form naturally in the compressed collision zone of the two superbubbles. Their shapes resemble the elongated, irregular structure of observed cold, molecular gas filaments and clumps. At the end of the simulations, between 65 and 80 percent of the total gas mass in our simulation box is contained in these structures. The clumps are found in a variety of physical states, ranging from pressure equilibrium with the surrounding medium to highly under-pressured clumps with large irregular internal motions and structures which are rotationally supported.

Numerical models of disc galaxy formation predict the existence of extended, hot ( T~10^6 K) gas haloes around present day spirals. The X-ray luminosity of these haloes is predicted to increase strongly with galaxy mass. However, searches for their X-ray emission have not been successful so far. We calculate the all sky O VII column density distributions for the haloes of three Milky Way like disc galaxies, resulting from cosmological high-resolution, N-body/gasdynamical simulations. We perform calculations both including the disc gas and without it, so the disc contribution to the column density is quantified. It is found that the column densities estimated for Milky Way-like galaxies are just below the observational upper limit, making a test of the hot halo paradigm likely within observational reach.