Search

Hub-filament systems are considered as natural sites for high-mass star formation. Kinematic analysis of the surroundings of hub-filaments is essential to better understand high-mass star formation within such systems. In this work, we present a detailed study of the massive Galactic protocluster IRAS 15394$-$5358, using continuum and molecular line data from the ALMA Three-millimeter Observations of Massive Star-forming Regions (ATOMS) survey. The 3~mm dust continuum map reveals the fragmentation of the massive ($\rm M=843~M_{\odot}$) clump into six cores. The core C-1A is the largest (radius = 0.04~pc), the most massive ($\rm M=157~M_{\odot}$), and lies within the dense central region, along with two smaller cores ($\rm M=7~and~3~M_{\odot}$). The fragmentation process is consistent with the thermal Jeans fragmentation mechanism and virial analysis shows that all the cores have small virial parameter values ($\rm \alpha_{vir}<<2$), suggesting that the cores are gravitationally bound. The mass vs. radius relation indicates that three cores can potentially form at least a single massive star. The integrated intensity map of $\rm H^{13}CO^{+}$ shows that the massive clump is associated with a hub-filament system, where the central hub is linked with four filaments. A sharp velocity gradient is observed towards the hub, suggesting a global collapse where the filaments are actively feeding the hub. We discuss the role of global collapse and the possible driving mechanisms for the massive star formation activity in the protocluster., Comment: 23 pages, 19 figures, accepted for publication in MNRAS

We explore the Schmidt-Kennicutt (SK) relations and the star formation efficiency per free-fall time ($\eff$), mirroring observational studies, in numerical simulations of filamentary molecular clouds undergoing gravitational contraction. We find that {\it a)} collapsing clouds accurately replicate the observed SK relations for galactic clouds and {\it b)} the so-called efficiency per free-fall time ($\eff$) is small and constant in space and in time, with values similar to those found in local clouds. This constancy is a consequence of the similar radial scaling of the free-fall time and the internal mass in density structures with spherically-averaged density profiles near $r^{-2}$. We additionally show that {\it c)} the star formation rate (SFR) increases rapidly in time; {\it d)} the low values of $\eff$ are due to the different time periods over which $\tauff$ and $\tausf$ are evaluated, together with the fast increasing SFR, and {\it e)} the fact that star clusters are significantly denser than the gas clumps from which they form is a natural consequence of the fast increasing SFR, the continuous replenishment of the star-forming gas by the accretion flow, and the near $r^{-2}$ density profile generated by the collapse Finally, we argue that the interpretation of $\eff$ as an efficiency is problematic because its maximum value is not bounded by unity, and because the total gas mass in the clouds is not fixed, but rather depends on the environment where clouds are embedded. In summary, our results show that the SK relation, the typical observed values of $\eff$, and the mass density of clusters arise as a natural consequence of gravitational contraction., Comment: Submitted to MNRAS

We provide a detailed comparison between the ``gravoturbulent'' (GT) and ``global hierarchical collapse'' (GHC) models for molecular clouds and star formation, their respective interpretations of the observational data, the features they share, and suggested tests and observations to discern between them. Also, we clarify common misconceptions in recent literature about the global and hierarchical nature of the GHC scenario, and briefly discuss the evolution of some aspects of both models toward convergence. GT assumes that molecular clouds and their substructures are in approximate virial equilibrium and are in a near-stationary state, interprets the linewidth-derived nonthermal motions exclusively as turbulence, which provides additional pressure against self-gravity. Conversely, GHC assumes that most star-forming molecular clouds and their substructures are part of a continuous gravitationally-driven flow, each accreting from their parent structure. Thus, the clouds and their star formation rate evolve in time. GHC interprets nonthermal motions as a mixture of infall and turbulent components, with the relative importance of the former increasing as the objects become denser and/or more massive. Tests that may provide clues to distinguishing between GT and GHC must take into account that the innermost parts of globally gravitationally bound structures may not locally appear bound, and thus the binding may have to be searched for at the largest scale of the structure., Comment: Submitted to MNRAS. Comments welcome. Animation for Fig. 3 available at https://www.irya.unam.mx/gente/g.gomez/publica/24-GT_vs_GHC-f3.mpg

We investigate the alignment mechanism between the magnetic field and interstellar clouds formed by the collision of warm atomic gas. We find that the magnetic field, initially oriented parallel to the flow, is perturbed by a fast MHD shock, which amplifies magnetic field fluctuations parallel to the shock front. Behind the shock, the compressive downstream velocity field further amplifies the magnetic field component parallel to the shock front. This mechanism causes the magnetic field to become increasingly parallel to the dense layer, and the development of a shear flow around the latter. Furthermore, the bending-mode perturbations on the dense layer are amplified by the non-linear thin-shell instability (NTSI), stretching the density structures formed by the thermal instability, and rendering them parallel to the bent field lines. By extension, we suggest that a tidal stretching velocity gradient such as that produced in gas infalling into a self-gravitating structure must straighten the field lines along the accretion flow, orienting them perpendicular to the density structures. We also find that the upstream superalfv\'enic regime transitions to a transalfv\'enic regime between the shock and the condensation front, and then to a subalfv\'enic regime inside the condensations. Our results provide a feasible physical mechanism for the observed transition from parallel to perpendicular relative orientation of the magnetic field and the density structures as the density structures become increasingly dominated by self-gravity., Comment: 11 pages, 10 figures, 3 animations. Submitted to Astronomy & Astrophysics. Comments are welcome

Only a handful of massive starless core candidates have been discovered so far, but none of them have been fully confirmed. Within the MM1 clump in the filamentary infrared dark cloud G34.43+0.24 that was covered by the ALMA-ATOMS survey at Band 3 ($\sim2\arcsec$, 6000\,au) and the ALMA-QUARKS survey at Band 6 ($\sim 0.3\arcsec$, 900\,au), two prestellar core candidates MM1-C and E1 with masses of 71 and 20 \solarmass~and radii of 2100--4400\,au were discovered. The two cores show no obvious sign of star-formation activities. In particular, MM1-C is a very promising massive prestellar core candidate with a total gas mass of 71\,\solarmass. Within MM1-C, we detected two extremely dense substructures, C1 and C2, as characterized by their high densities of $\rm n_{H_2}\sim 10^{8-9} cm^{-3}$. Moreover, evidence of further fragmentation in C2 was also revealed. We have detected the primordial fragmentation in the earliest stage of massive star formation, and we speculate that MM1-C would be the birthplace of a massive multiple system. However, we cannot fully rule out the possibility that the massive prestellar core MM1-C will just form a cluster of low-mass stars if it undergoes further fragmentation., Comment: 12 pages, 6 figures

The ALMA Survey of Star Formation and Evolution in Massive Protoclusters with Blue Profiles (ASSEMBLE) aims to investigate the process of mass assembly and its connection to high-mass star formation theories in protoclusters in a dynamic view. We observed 11 massive (Mclump>1000 Msun), luminous (Lbol>10,000 Lsun), and blue-profile (infall signature) clumps by ALMA with resolution of 2200-5500 au at 350 GHz (870 um) in continuum and line emission. 248 dense cores were identified, including 106 cores showing protostellar signatures and 142 prestellar core candidates. Compared to early-stage infrared dark clouds (IRDCs) by ASHES, the core mass and surface density within the ASSEMBLE clumps exhibited significant increment, suggesting concurrent core accretion during the evolution of the clumps. The maximum mass of prestellar cores was found to be 2 times larger than that in IRDCs, indicating evolved protoclusters have the potential to harbor massive prestellar cores. The mass relation between clumps and their most massive core (MMCs) is observed in ASSEMBLE but not in IRDCs, which is suggested to be regulated by multiscale mass accretion. The mass correlation between the core clusters and their MMCs has a steeper slope compared to that observed in stellar clusters, which can be due to fragmentation of the MMC and stellar multiplicity. We observe a decrease in core separation and an increase in central concentration as protoclusters evolve. We confirm primordial mass segregation in the ASSEMBLE protoclusters, possibly resulting from gravitational concentration and/or gas accretion., Comment: 37 pages, 13 figures, 5 tables; accepted for publication in ApJS

The formation mechanism of massive stars remains one of the main open problems in astrophysics, in particular the relationship between the mass of the most massive stars, and that of the cores in which they form. Numerical simulations of the formation and evolution of large molecular clouds, within which dense cores and stars form self-consistently, show in general that the cores' masses increase in time, and also that the most massive stars tend to appear later (by a few to several Myr) than lower-mass stars. Here we present an idealized model that incorporates accretion onto the cores as well as onto the stars, in which the core's mass growth is regulated by a ``gravitational choking'' mechanism that does not involve any form of support. This process is of purely gravitational origin, and causes some of the mass accreted onto the core to stagnate there, rather than being transferred to the central stars. Thus, the simultaneous mass growth of the core and of the stellar mass can be computed. In addition, we estimate the mass of the most massive allowed star before its photoionizing radiation is capable of overcoming the accretion flow onto the core. This model constitutes a proof-of-concept for the simultaneous growth of the gas reservoir and the stellar mass, the delay in the formation of massive stars observed in cloud-scale numerical simulations, the need for massive, dense cores in order to form massive stars, and the observed correlation between the mass of the most massive star and the mass of the cluster it resides in. Also, our model implies that by the time massive stars begin to form in a core, a number of low-mass stars are expected to have already formed., Comment: Resubmitted to MNRAS after incorporating reviewer's recommendations. Comments welcome

There is growing evidence that high-mass star formation and hub-filament systems (HFS) are intricately linked. The gas kinematics along the filaments and the forming high-mass star(s) in the central hub are in excellent agreement with the new generation of global hierarchical high-mass star formation models. In this paper, we present an observational investigation of a typical HFS cloud, G310.142+0.758 (G310 hereafter) which reveals unambiguous evidence of mass inflow from the cloud scale via the filaments onto the forming protostar(s) at the hub conforming with the model predictions. Continuum and molecular line data from the ATOMS and MALT90 surveys are used that cover different spatial scales. Three filaments (with total mass $5.7\pm1.1\times 10^3~M_{\odot}$) are identified converging toward the central hub region where several signposts of high-mass star formation have been observed. The hub region contains a massive clump ($1280\pm260~M_{\odot}$) harbouring a central massive core. Additionally, five outflow lobes are associated with the central massive core implying a forming cluster. The observed large-scale, smooth and coherent velocity gradients from the cloud down to the core scale, and the signatures of infall motion seen in the central massive clump and core, clearly unveil a nearly-continuous, multi-scale mass accretion/transfer process at a similar mass infall rate of $\sim 10^{-3}~M_{\odot}~yr^{-1}$ over all scales, feeding the central forming high-mass protostar(s) in the G310 HFS cloud., Comment: Accepted to publish in ApJ. 10 pages with 6 figures and 2 tables

We present ALMA Band-3/7 observations towards "the Heart" of a massive hub-filament system (HFS) SDC335, to investigate its fragmentation and accretion. At a resolution of $\sim0.03$ pc, 3 mm continuum emission resolves two massive dense cores MM1 and MM2, with $383(^{+234}_{-120})$ $M_\odot$ (10-24% mass of "the Heart") and $74(^{+47}_{-24})$ $M_\odot$, respectively. With a resolution down to 0.01 pc, 0.87 mm continuum emission shows MM1 further fragments into six condensations and multi-transition lines of H$_2$CS provide temperature estimation. The relation between separation and mass of condensations at a scale of 0.01 pc favors turbulent Jeans fragmentation where the turbulence seems to be scale-free rather than scale-dependent. We use the H$^{13}$CO$^+$ (1-0) emission line to resolve the complex gas motion inside "the Heart" in position-position-velocity space. We identify four major gas streams connected to large-scale filaments, inheriting the anti-clockwise spiral pattern. Along these streams, gas feeds the central massive core MM1. Assuming an inclination angle of $45(\pm15)^{\circ}$ and a H$^{13}$CO$^+$ abundance of $5(\pm3)\times10^{-11}$, the total mass infall rate is estimated to be $2.40(\pm0.78)\times10^{-3}$ $M_\odot$ yr$^{-1}$, numerically consistent with the accretion rates derived from the clump-scale spherical infall model and the core-scale outflows. The consistency suggests a continuous, near steady-state, and efficient accretion from global collapse, therefore ensuring core feeding. Our comprehensive study of SDC335 showcases the detailed gas kinematics in a prototypical massive infalling clump and calls for further systematic and statistical analyses in a large sample., Comment: 29 pages, 7 figures, Accepted for publication in MNRAS

We investigate the presence of hub-filament systems in a large sample of 146 active proto-clusters, using H$^{13}$CO$^{+}$ J=1-0 molecular line data obtained from the ATOMS survey. We find that filaments are ubiquitous in proto-clusters, and hub-filament systems are very common from dense core scales ($\sim$0.1 pc) to clump/cloud scales ($\sim$1-10 pc). The proportion of proto-clusters containing hub-filament systems decreases with increasing dust temperature ($T_d$) and luminosity-to-mass ratios ($L/M$) of clumps, indicating that stellar feedback from H{\sc ii} regions gradually destroys the hub-filament systems as proto-clusters evolve. Clear velocity gradients are seen along the longest filaments with a mean velocity gradient of 8.71 km s$^{-1}$pc$^{-1}$ and a median velocity gradient of 5.54 km s$^{-1}$pc$^{-1}$. We find that velocity gradients are small for filament lengths larger than $\sim$1~pc, probably hinting at the existence of inertial inflows, although we cannot determine whether the latter are driven by large-scale turbulence or large-scale gravitational contraction. In contrast, velocity gradients below $\sim$1~pc dramatically increase as filament lengths decrease, indicating that the gravity of the hubs or cores starts to dominate gas infall at small scales. We suggest that self-similar hub-filament systems and filamentary accretion at all scales may play a key role in high-mass star formation., Comment: 16 pages

We investigate the effect of numerical magnetic diffusion in magnetohydrodynamic simulations of magnetically supported molecular clouds. To this end, we have performed numerical studies on adaptive mesh isothermal simulations of marginally sub-critical molecular clouds. We find that simulations with low and intermediate resolutions collapse, contrary to what is theoretically expected. However, the simulation with the highest numerical resolution oscillates around an equilibrium state without collapsing. In order to quantify the numerical diffusion of the magnetic field, we ran a second suit of current-sheet simulations in which the numerical magnetic diffusion coefficient can be directly measured, and computed the corresponding diffusion times at various numerical resolutions. On this basis, we propose a criterion for the resolution of magnetic fields in MHD simulations based on requiring that the diffusion time to be larger than the characteristic timescale of the physical process responsible for the dynamic evolution of the structure., Comment: 7 pages, 5 figures, 2 tables. Accepted to MNRAS

We explore how dense filament widths, when measured using different molecular species, may change as a consequence of gas accretion toward the filament. As a gas parcel falls into the filament, it will experience different density, temperature, and extinction values. The rate at which this environment changes will affect differently the abundance of different molecules. So, a molecule that forms quickly will better reflect the local physical conditions a gas parcel experiences than a slower-forming molecule. Since these differences depend on how the respective timescales compare, the different molecular distributions should reflect how rapidly the environment changes, i.e., the accretion rate toward the filament. We find that the filament widths measured from time-dependent abundances for C2H, CO, CN, CS, and C3H2, are sensitive the most to this effect, being those molecules the ones presenting also the wider filament widths. On the contrary, molecules such as N2H+, NH3, H2CO, HNC or CH3OH are not so sensitive to accretion and present the narrowest filament widths. We propose that ratios of filament widths for different tracers could be a useful tool to estimate the accretion rate onto the filament., Comment: Accepted for publication in MNRAS

We study the gravitationally-dominated, accretion-driven evolution of a prestellar core. In our model, as the core's density increases, it remains immersed in a constant-density environment and so it accretes from this environment, increasing its mass and reducing its Jeans length. Assuming a power-law density profile $\rho \propto r^{-p}$, we compute the rate of change of the slope $p$, and show that the value $p=2$ is stationary, and furthermore, an attractor. The radial profile of the Jeans length scales as $r^{p/2}$, implying that, for $p<2$, there is a radius below which the region is smaller than its Jeans length, thus appearing gravitationally stable and in need of pressure confinement, while, in reality, it is part of a larger-scale collapse and is undergoing compression by the infalling material. In this region, the infall speed decreases towards the center, eventually becoming subsonic, thus appearing "coherent", without the need for turbulence dissipation. We present a compilation of observational determinations of density profiles in dense cores and show that the distribution of their slopes peaks at $p \sim 1.7$--1.9, supporting the notion that the profile steepens over time. Finally, we discuss the case of magnetic support in a core in which the field scales as $B \propto \rho^\beta$. For the expected value of $\beta = 2/3$, this implies that the mass to magnetic flux ratio also decreases towards the central parts of the cores, making them appear magnetically supported, while in reality they may be part of larger collapsing supercritical region. We conclude that local signatures of either thermal or magnetic support are not conclusive evidence of stability, that the gravitational instability of a region must be established at the large scales, and that the prestellar stage of collapse is dynamic rather than quasistatic., Comment: Accepted for publication in MNRAS

We present a unified description of the scenario of Global Hierarchical Collapse and fragmentation (GHC) in molecular clouds (MCs), owing to the continuous decrease of the average Jeans mass in the contracting cloud. GHC constitutes a regime of collapses within collapses, in which small-scale collapses begin at later times, but occur on shorter timescales than large-scale ones. The difference in timescales allows for most of the clouds' mass to be dispersed by feedback from the first massive stars, maintaining the global star formation rate low. All scales accrete from their parent structures. The main features of GHC are: star-forming MCs are in an essentially pressureless regime, which produces filaments that accrete onto clumps and cores ("hubs"). The filaments constitute the collapse flow from cloud to hub scales and may approach a quasi-stationary state; the molecular and dense mass fractions of the clouds increase over time; the first (low-mass) stars appear several Myr after global contraction began; more massive stars appear after a few Myr in massive hubs resulting from the collapse of larger scales; the minimum fragment mass may extend well into the brown-dwarf regime; Bondi-Hoyle-Lyttleton accretion occurs at the protostellar and core scales, accounting for a near-Salpeter IMF; the extreme anisotropy of the filamentary network explains the difficulty in detecting large-scale infall signatures; the balance between inertial and gravitationally-driven motions in clumps evolves during the contraction; prestellar cores adopt Bonnor-Ebert-like profiles, but are contracting ever since early times when they may appear to be unbound and to require pressure confinement; stellar clusters develop radial age and mass segregation gradients. Finally, we discuss the incompatibility between supersonic turbulence and the observed scalings in the molecular hierarchy., Comment: Accepted by MNRAS. Sec. 6. 4 has been modified to provide an evolutionary expression for the viral parameter, including turbulent and self-gravitating components. New table comparing various SF scenarios. Abstract abridged

We explore the structure of magnetic field lines in and around filaments in simulations of molecular clouds undergoing global, multi-scale gravitational collapse. In these simulations, filaments are not in a static equilibrium, but are long-lived flow structures that accrete gas from their environment and direct it toward clumps embedded in the filament or at the nodes at the conjunction with other filaments. In this context, the magnetic field is dragged by the collapsing gas, so its structure must reflect the flow that generates the filament. Around the filament, the gas is accreted onto it, and the magnetic lines must then be perpendicular to the filament. As the gas density increases, the gas flow changes direction, becoming almost parallel to the filament, and magnetic lines also tend to align with it. At the spine of the filament, however, magnetic lines become perpendicular again since they must connect to lines on the opposite side of the filament, resulting in "U"-shaped magnetic structures, which tend to be stretched by the longitudinal flow along the filament. Magnetic diffusive processes, however, allow the gas to continue to flow. Assuming a stationary state in which the ram pressure of the flow balances the magnetic tension, the curvature of the field lines is determined by the diffusion rate. We derive an expression relating the curvature of the field lines to the diffusive coefficient, which may be used to observationally determine the nature of the diffusive process., Comment: Submitted to MNRAS

The removal of the interstellar medium (ISM) of disk galaxies through ram pressure stripping (RPS) has been extensively studied in numerous simulations. Nevertheless, the role of magnetic fields (MF) on the gas dynamics in this process has been hardly studied, although the MF influence on the large-scale disk structure is well established. With this in mind, we present a 3D magnetohydrodynamic (MHD) simulation of face-on RPS of a disk galaxy to study the impact of the galactic MF in the gas stripping. The main effect of including a galactic MF is a flared disk. When the intracluster medium (ICM) wind hits this flared disk, oblique shocks are produced at the interaction interface, where the ISM is compressed, generating a gas inflow from large radii towards the central regions of the galaxy. This inflow is observed for $\sim 150$ Myr and may supply the central parts of the galaxy with material for star formation while the outskirts of the disk are being stripped of gas, thus the oblique shocks can induce and enhance the star formation in the remaining disk. We also observed that the MF alters the shape and structure of the swept gas, giving a smooth appearance in the magnetized case and clumpier and filamentary-like morphology in the hydro case. Finally, we estimated the truncation radius expected for our models using the Gunn-Gott criterion and found that is in agreement with the simulations., Comment: 14 pages, 8 figures. Accepted for publication to MNRAS

We analyze the physical properties and energy balance of density enhancements in two SPH simulations of the formation, evolution, and collapse of giant molecular clouds. In the simulations, no feedback is included, so all motions are due either to the initial, decaying turbulence, or to gravitational contraction. We define clumps as connected regions above a series of density thresholds. The resultingfull set of clumps follows the generalized energy-equipartition relation $\sigma_{v}/R^{1/2} \propto \Sigma^{1/2}$, where $\sigma_{v}$ is the velocity dispersion, $R$ is the "radius", and $\Sigma$ is the column density. We interpret this as a natural consequence of gravitational contraction at all scales, rather than virial equilibrium. Nevertheless, clumps with low $\Sigma$ tend to show a large scatter around equipartition. In more than half of the cases, this scatter is dominated by external turbulent compressions that {\it assemble} the clumps, rather than by small-scale random motions that would disperse them. The other half does actually disperse. Moreover, clump sub-samples selected by means of different criteria exhibit different scalings. Sub-samples with narrow $\Sigma$ ranges follow Larson-like relations, although characterized by their respective value of $\Sigma$. Finally, we find that: i) clumps lying in filaments tend to appear sub-virial; ii) high-density cores ($n \ge 10^5$ cm$^3$) that exhibit moderate kinetic energy excesses often contain sink ("stellar") particles, and the excess disappears when the stellar mass is taken into account in the energy balance; iii) cores with kinetic energy excess but no stellar particles are truly in a state of dispersal., Comment: Abstract rephrased. Submitted to ApJ

With the use of a background Milky-Way-like potential model, we performed stellar orbital and magnetohydrodynamic (MHD) simulations. As a first experiment, we studied the gaseous response to a bisymmetric spiral arm potential: the widely employed cosine potential model and a self-gravitating tridimensional density distribution based model called PERLAS. Important differences are noticeable in these simulations, while the simplified cosine potential produces two spiral arms for all cases, the more realistic density based model produces a response of four spiral arms on the gaseous disk, except for weak arms -i.e. close to the linear regime- where a two-armed structure is formed. In order to compare the stellar and gas response to the spiral arms, we have also included a detailed periodic orbit study and explored different structural parameters within observational uncertainties. The four armed response has been explained as the result of ultra harmonic resonances, or as shocks with the massive bisymmetric spiral structure, among other. From the results of this work, and comparing the stellar and gaseous responses, we tracked down an alternative explanation to the formation of branches, based only on the orbital response to a self-gravitating spiral arms model. The presence of features such as branches, might be an indication of transiency of the arms., Comment: 17 pages, 9 figures. Accepted for publication in MNRAS

Motivated by recent numerical simulations of molecular cloud (MC)evolution, in which the clouds engage in global gravitational contraction, and local collapse events culminate significantly earlier than the global collapse, we investigate the growth of density perturbations embedded in a collapsing background, to which we refer as an Inverse Hubble Flow (IHF). We use the standard procedure for the growth of perturbations in a universe that first expands (the usual Hubble Flow) and then recollapses (the IHF). We find that linear density perturbations immersed in an IHF grow faster than perturbations evolving in a static background (the standard Jeans analysis). A fundamental distinction between the two regimes is that, in the Jeans case, the time $\tau_\mathrm{nl}$ for a density fluctuation to become nonlinear increases without limit as its initial value approaches zero, while in the IHF case $\tau_\mathrm{nl} \le \tau_\mathrm{ff}$ always, where $\tau_\mathrm{ff}$ is the free-fall time of the background density. We suggest that this effect, although moderate, implies that small-scale density fluctuations embedded in globally-collapsing clouds must collapse earlier than their parent cloud, regardless of whether the initial amplitude of the fluctuations is moderate or strongly nonlinear, thus allowing the classical mechanism of Hoyle fragmentation to operate in multi-Jeans-mass MCs. More fundamentally, our results show that, contrary to the standard paradigm that fluctuations of all scales grow at the same rate in the linear regime, the hierarchical nesting of the fluctuations of different scales does affect their growth even in the linear stage., Comment: 7 pages, 3 figures; published in MNRAS

We analyze the spiral structure that results in a numerical simulation of a galactic disk with stellar and gaseous components evolving in a potential that includes an axisymmetric halo and bulge. We perform a second simulation without the gas component to observe how it affects the spiral structure in the disk. To quantify this, we use a Fourier analysis and obtain values for the pitch angle and the velocity of the self-excited spiral pattern of the disk. The results show a tighter spiral in the simulation with gaseous component. The spiral structure is consistent with a superposition of waves, each with a constant pattern velocity in given radial ranges.

We report on the filaments that develop self-consistently in a new numerical simulation of cloud formation by colliding flows. As in previous studies, the forming cloud begins to undergo gravitational collapse because it rapidly acquires a mass much larger than the average Jeans mass. Thus, the collapse soon becomes nearly pressureless, proceeding along its shortest dimension first. This naturally produces filaments in the cloud, and clumps within the filaments. The filaments are not in equilibrium at any time, but instead are long-lived flow features, through which the gas flows from the cloud to the clumps. The filaments are long-lived because they accrete from their environment while simultaneously accreting onto the clumps within them; they are essentially the locus where the flow changes from accreting in two dimensions to accreting in one dimension. Moreover, the clumps also exhibit a hierarchical nature: the gas in a filament flows onto a main, central clump, but other, smaller-scale clumps form along the infalling gas. Correspondingly, the velocity along the filament exhibits a hierarchy of jumps at the locations of the clumps. Two prominent filaments in the simulation have lengths ~15 pc, and masses ~600 Msun above density n ~ 10^3 cm-3 (~2x10^3 Msun at n > 50 cm-3). The density profile exhibits a central flattened core of size ~0.3 pc and an envelope that decays as r^-2.5, in reasonable agreement with observations. Accretion onto the filament reaches a maximum linear density rate of ~30 Msun Myr^-1 pc^-1., Comment: Revised to address the referee's comments, submitted to ApJ. See related animations in http://www.crya.unam.mx/~g.gomez/publica.html

We present a numerical study of the evolution of molecular clouds, from their formation by converging flows in the warm ISM, to their destruction by the ionizing feedback of the massive stars they form. We improve with respect to our previous simulations by including a different stellar-particle formation algorithm, which allows them to have masses corresponding to single stars rather than to small clusters, and with a mass distribution following a near-Salpeter stellar IMF. We also employ a simplified radiative-transfer algorithm that allows the stellar particles to feed back on the medium at a rate that depends on their mass and the local density. Our results are as follows: a) Contrary to the results from our previous study, where all stellar particles injected energy at a rate corresponding to a star of ~ 10 Msun, the dense gas is now completely evacuated from 10-pc regions around the stars within 10-20 Myr, suggesting that this feat is accomplished essentially by the most massive stars. b) At the scale of the whole numerical simulations, the dense gas mass is reduced by up to an order of magnitude, although star formation (SF) never shuts off completely, indicating that the feedback terminates SF locally, but new SF events continue to occur elesewhere in the clouds. c) The SF efficiency (SFE) is maintained globally at the ~ 10% level, although locally, the cloud with largest degree of focusing of its accretion flow reaches SFE ~ 30%. d) The virial parameter of the clouds approaches unity before the stellar feedback begins to dominate the dynamics, becoming much larger once feedback dominates, suggesting that clouds become unbound as a consequence of the stellar feedback. e) The erosion of the filaments that feed the star-forming clumps produces chains of isolated dense blobs reminiscent of those observed in the vicinity of the dark globule B68., Comment: Accepted for publication in MNRAS. This version include two more figures and four more pages with a more extended discusion

It is generally assumed that gas in a galactic disk follows closely non self-intersecting periodic stellar orbits. In order to test this common assumption, we have performed MHD simulations of a galactic-like disk under the influence of a spiral galactic potential. We also have calculated the actual orbit of a gas parcel and compared it to stable periodic stellar orbits in the same galactic potential and position. We found that the gaseous orbits approach periodic stellar orbits far from the major orbital resonances only. Gas orbits initialized at a given galactocentric distance but at different azimuths can be different, and scattering is conspicuous at certain galactocentric radii. Also, in contrast to the stellar behaviour, near the 4:1 (or higher order) resonance the gas follows nearly circular orbits, with much shorter radial excursions than the stars. Also, since the gas does not settle into a steady state, the gaseous orbits do not necessarily close on themselves., Comment: 8 pages, accepted for publication in MNRAS

We investigate the collapse of non-spherical substructures, such as sheets and filaments, which are ubiquitous in molecular clouds. Such non-spherical substructures collapse homologously in their interiors but are influenced by an edge effect that causes their edges to be preferentially accelerated. We analytically compute the homologous collapse timescales of the interiors of uniform-density, self-gravitating filaments and find that the homologous collapse timescale scales linearly with the aspect ratio. The characteristic timescale for an edge driven collapse mode in a filament, however, is shown to have a square root dependence on the aspect ratio. For both filaments and circular sheets, we find that selective edge acceleration becomes more important with increasing aspect ratio. In general, we find that lower dimensional objects and objects with larger aspect ratios have longer collapse timescales. We show that estimates for star formation rates, based upon gas densities, can be overestimated by an order of magnitude if the geometry of a cloud is not taken into account., Comment: 10 pages, 2 figures, accepted by ApJ, minor grammatical errors fixed

Molecular clouds often exhibit filamentary or sheet-like shapes. We compute the free-fall time ($\tff$) for finite, uniform, self-gravitating circular sheets and filamentary clouds of small but finite thickness, so that their volume density $\rho$ can still be defined. We find that, for thin sheets, the free-fall time is larger than that of a uniform sphere with the same volume density by a factor proportional to $\sqrt{A}$, where the aspect ratio $A$ is given by $A=R/h$, $R$ being the sheet's radius and $h$ is its thickness. For filamentary clouds, the aspect ratio is defined as $A=L/\calR$, where $L$ is the filament's half length and $\calR$ is its (small) radius, and the modification factor is a more complicated, although in the limit of large $A$ it again reduces to nearly $\sqrt{A}$. We propose that our result for filamentary shapes naturally explains the ubiquitous configuration of clumps fed by filaments observed in the densest structures of molecular clouds. Also, the longer free-fall times for non-spherical geometries in general may contribute towards partially alleviating the "star-formation conundrum", namely, that the star formation rate in the Galaxy appears to be proceeding in a timescale much larger than the total molecular mass in the Galaxy divided by its typical free-fall time. If molecular clouds are in general formed by thin sheets and long filaments, then their relevant free-fall time may have been systematically underestimated, possibly by factors of up to one order of magnitude., Comment: To appear on The Astrophysical Journal

We report a set of numerical experiments aimed at addressing the applicability of competitive accretion to explain the high-mass end of the stellar initial mass function in a sheet geometry with shallow gravitational potential, in contrast to most previous simulations which have assumed formation in a cluster gravitational potential. Our flat cloud geometry is motivated by models of molecular cloud formation due to large-scale flows in the interstellar medium. The experiments consisted of SPH simulations of gas accretion onto sink particles formed rapidly from Jeans-unstable dense clumps placed randomly in the finite sheet. These simplifications allow us to study accretion with a minimum of free parameters, and to develop better statistics on the resulting mass spectra. We considered both clumps of equal mass and gaussian distributions of masses, and either uniform or spatially-varying gas densities. In all cases, the sink mass function develops a power law tail at high masses, with $dN/dlog M \propto M^{-\Gamma}$. The accretion rates of individual sinks follow $\dot{M} \propto M^2$ at high masses; this results in a continual flattening of the slope of the mass function towards an asymptotic form $\Gamma \sim 1$ (where the Salpeter slope is $\Gamma = 1.35$). The asymptotic limit is most rapidly reached when starting from a relatively broad distribution of initial sink masses. In general the resulting upper mass slope is correlated with the maximum sink mass; higher sink masses are found in simulations with flatter upper mass slopes. Although these simulations are of a highly idealized situation, the results suggest that competitive accretion may be relevant in a wider variety of environments than previously considered, and in particular that the upper mass distribution may generally evolve towards a limiting value of $\Gamma \sim 1$., Comment: 20 pages, 12 figures

We numerically investigate the effect of feedback from the ionizing radiation heating from massive stars on the evolution of giant molecular clouds (GMCs) and their star formation efficiency (SFE). We find that the star-forming regions within the GMCs are invariably formed by gravitational contraction. After an initial period of contraction, the collapsing clouds begin forming stars, whose feedback evaporates part of the clouds' mass, opposing the continuing accretion from the infalling gas. The competition of accretion against dense gas consumption by star formation (SF) and evaporation by the feedback, regulates the clouds' mass and energy balance, as well as their SFE. We find that, in the presence of feedback, the clouds attain levels of the SFE that are consistent at all times with observational determinations for regions of comparable SF rates (SFRs). However, we observe that the dense gas mass is larger in general in the presence of feedback, while the total (dense gas + stars) is nearly insensitive to the presence of feedback, suggesting that the total mass is determined by the accretion, while the feedback inhibits mainly the conversion of dense gas to stars. The factor by which the SFE is reduced upon the inclusion of feedback is a decreasing function of the cloud's mass, for clouds of size ~ 10 pc. This naturally explains the larger observed SFEs of massive-star forming regions. We also find that the clouds may attain a pseudo-virialized state, with a value of the virial mass very similar to the actual cloud mass. However, this state differs from true virialization in that the clouds are the center of a large-scale collapse, continuously accreting mass, rather than being equilibrium entities., Comment: Submitted to ApJ (abstract abridged)

We investigate the response of the star formation efficiency (SFE) to the main parameters of simulations of molecular cloud formation by the collision of warm diffuse medium (WNM) cylindrical streams, neglecting stellar feedback and magnetic fields. The parameters we vary are the Mach number of the inflow velocity of the streams, Msinf, the rms Mach number of the initial background turbulence in the WNM, and the total mass contained in the colliding gas streams, Minf. Because the SFE is a function of time, we define two estimators for it, the "absolute" SFE, measured at t = 25 Myr into the simulation's evolution (sfeabs), and the "relative" SFE, measured 5 Myr after the onset of star formation in each simulation (sferel). The latter is close to the "star formation rate per free-fall time" for gas at n = 100 cm^-3. We find that both estimators decrease with increasing Minf, although by no more than a factor of 2 as Msinf increases from 1.25 to 3.5. Increasing levels of background turbulence similarly reduce the SFE, because the turbulence disrupts the coherence of the colliding streams, fragmenting the cloud, and producing small-scale clumps scattered through the numerical box, which have low SFEs. Finally, the SFE is very sensitive to the mass of the inflows, with sferel decreasing from ~0.4 to ~0.04 as the the virial parameter in the colliding streams increases from ~0.15 to ~1.5. This trend is in partial agreement with the prediction by Krumholz & McKee (2005), since the latter lies within the same range as the observed efficiencies, but with a significantly shallower slope. We conclude that the observed variability of the SFE is a highly sensitive function of the parameters of the cloud formation process, and may be the cause of significant scatter in observational determinations., Comment: 19 pages, submitted to MNRAS

Assuming that density waves trigger star formation, and that young stars preserve the velocity components of the molecular gas where they are born, we analyze the effects that non-circular gas orbits have on color gradients across spiral arms. We try two approaches, one involving semi-analytical solutions for spiral shocks, and another with magnetohydrodynamic (MHD) numerical simulation data. We find that, if non-circular motions are ignored, the comparison between observed color gradients and stellar population synthesis models would in principle yield pattern speed values that are systematically too high for regions inside corotation, with the difference between the real and the measured pattern speeds increasing with decreasing radius. On the other hand, image processing and pixel averaging result in systematically lower measured spiral pattern speed values, regardless of the kinematics of stellar orbits. The net effect is that roughly the correct pattern speeds are recovered, although the trend of higher measured $\Omega_p$ at lower radii (as expected when non-circular motions exist but are neglected) should still be observed. We examine the Martinez-Garcia et al. (2009) photometric data and confirm that this is indeed the case. The comparison of the size of the systematic pattern speed offset in the data with the predictions of the semi-analytical and MHD models corroborates that spirals are more likely to end at Outer Lindblad Resonance, as these authors had already found., Comment: 32 pages, 15 figures, accepted to ApJ

We investigate the properties of "star forming regions" in a previously published numerical simulation of molecular cloud formation out of compressive motions in the warm neutral atomic interstellar medium, neglecting magnetic fields and stellar feedback. In this simulation, the velocity dispersions at all scales are caused primarily by infall motions rather than by random turbulence. We study the properties (density, total gas+stars mass, stellar mass, velocity dispersion, and star formation rate) of the cloud hosting the first local, isolated "star formation" event in the simulation and compare them with those of the cloud formed by a later central, global collapse event. We suggest that the small-scale, isolated collapse may be representative of low- to intermediate-mass star-forming regions, while the large-scale, massive one may be representative of massive star forming regions. We also find that the statistical distributions of physical properties of the dense cores in the region of massive collapse compare very well with those from a recent survey of the massive star forming region in the Cygnus X molecular cloud. The star formation efficiency per free-fall time (SFE_ff) of the high-mass SF clump is low, ~0.04. This occurs because the clump is accreting mass at a high rate, not because its specific SFR (SSFR) is low. This implies that a low value of the SFE_ff does not necessarily imply a low SSFR, but may rather indicate a large gas accretion rate. We suggest that a globally low SSFR at the GMC level can be attained even if local star forming sites have much larger values of the SSFR if star formation is a spatially intermittent process, so that most of the mass in a GMC is not participating of the SF process at any given time., Comment: Accepted by ApJ. Revised version, according to exchanges with referee. Original results unchanged. Extensive new discussion on the low global efficiency vs. high local efficiency of star formation. Abstract abridged

Only a few molecular clouds in the Solar Neighborhood exhibit the formation of only low-mass stars. Traditionally, these clouds have been assumed to be supported against more vigorous collapse by magnetic fields. The existence of strong magnetic fields in molecular clouds, however, poses serious problems for the formation of stars and of the clouds themselves. In this {\em Letter}, we review the three-dimensional structure and kinematics of Taurus --the archetype of a region forming only low-mass stars-- as well as its orientation within the Milky way. We conclude that the particularly low star-formation efficiency in Taurus may naturally be explained by tidal forces from the Galaxy, with no need for magnetic regulation or stellar feedback., Comment: Minor changes. 5 pages. Accepted by MNRAS

(Abridged) The gravitational term for clouds and cores entering in the virial theorem is usually assumed to be equal to the gravitational energy, since the contribution to the gravitational force from the mass distribution outside the volume of integration is assumed to be negligible. Such approximation may not be valid in the presence of an important external net potential. In the present work we analyze the effect of an external gravitational field on the gravitational budget of a density structure. Our cases under analysis are (a) a giant molecular cloud (GMC) with different aspect ratios embedded within a galactic net potential, and (b) a molecular cloud core embedded within the gravitational potential of its parent molecular cloud. We find that for roundish GMCs, the tidal tearing due to the shear in the plane of the galaxy is compensated by the tidal compression in the z direction. The influence of the external effective potential on the total gravitational budget of these clouds is relatively small, although not necessarily negligible. However, for more filamentary GMCs, the external effective potential can be dominant and can even overwhelm self-gravity, regardless of whether its main effect on the cloud is to disrupt it or compress it. This may explain the presence of some GMCs with few or no signs of massive star formation, such as the Taurus or the Maddalena's clouds. In the case of dense cores embedded in their parent molecular cloud, we found that the gravitational content due to the external field may be more important than the gravitational energy of the cores themselves. This effect works in the same direction as the gravitational energy, i.e., favoring the collapse of cores. We speculate on the implications of these results for star formation models., Comment: Accepted for publication in MNRAS

(Abridged) We present numerical hydrodynamical simulations of the formation, evolution and gravitational collapse of isothermal molecular cloud cores. A compressive wave is set up in a constant sub-Jeans density distribution of radius r = 1 pc. As the wave travels through the simulation grid, a shock-bounded spherical shell is formed. The inner shock of this shell reaches and bounces off the center, leaving behind a central core with an initially almost uniform density distribution, surrounded by an envelope consisting of the material in the shock-bounded shell, with a power-law density profile that at late times approaches a logarithmic slope of -2 even in non-collapsing cases. The resulting density structure resembles a quiescent core of radius < 0.1 pc, with a Bonnor-Ebert-like (BE-like) profile, although it has significant dynamical differences: it is initially non-self-gravitating and confined by the ram pressure of the infalling material, and consequently, growing continuously in mass and size. With the appropriate parameters, the core mass eventually reaches an effective Jeans mass, at which time the core begins to collapse. Thus, there is necessarily a time delay between the appearance of the core and the onset of its collapse, but this is not due to the dissipation of its internal turbulence as it is often believed. These results suggest that pre-stellar cores may approximate Bonnor-Ebert structures which are however of variable mass and may or may not experience gravitational collapse, in qualitative agreement with the large observed frequency of cores with BE-like profiles., Comment: Accepted for publication in ApJ. Associated mpeg files can be found in http://www.astrosmo.unam.mx/~g.gomez/publica.html

Errors in the kinematic distances, under the assumption of circular gas orbits, were estimated by performing synthetic observations of a model disk galaxy. It was found that the error is < 0.5 kpc for most of the disk when the measured rotation curve was used, but larger if the real rotation curve is applied. In both cases, the error is significantly larger at the positions of the spiral arms. The error structure is such that, when kinematic distances are used to develope a picture of the large scale density distribution, the most significant features of the numerical model are significantly distorted or absent, while spurious structure appears. By considering the full velocity field in the calculation of the kinematic distances, most of the original density structures can be recovered., Comment: Accepted for publication in AJ

We study the stability of proto-planetary disks with vertical velocity gradients in their equilibrium rotation rates; such gradients are expected to develop when dust settles into the midplane. Using a linear stability analysis of a simple three-layer model, we show that the onset of instability occurs at a larger value of the Richardson number, and therefore for a thicker layer, when the effects of Coriolis forces are included. This analysis also shows that even-symmetry (midplane-crossing) modes develop faster than odd-symmetry ones. These conclusions are corroborated by a large number of nonlinear numerical simulations with two different parameterized prescriptions for the initial (continuous) dust distributions. Based on these numerical experiments, the Richardson number required for marginal stability is more than an order of magnitude larger than the traditional 1/4 value. The dominant modes that grow have horizontal wavelengths of several initial dust scale heights, and in nonlinear stages mix solids fairly homogeneously over a comparable vertical range. We conclude that gravitational instability may be more difficult to achieve than previously thought, and that the vertical distribution of matter within the dust layer is likely globally, rather than locally, determined., Comment: Accepted for publication in ApJ

We generated synthetic observations from the four-arm model presented in Gomez & Cox (2004) for the Galactic ISM in the presence of a spiral gravitational perturbation. We found that velocity crowding and diffusion have a strong effect in the l-v diagram. The v-b diagram presents structures at the expected spiral arm velocities, that can be explained by the off-the-plane structure of the arms presented in previous papers of this series. Such structures are observed in the Leiden/Dwingeloo HI survey. The rotation curve, as measured from the inside of the modeled galaxy, shows similarities with the observed one for the Milky Way Galaxy, although it has large deviations from the smooth circular rotation corresponding to the background potential. The magnetic field inferred from a synthetic synchrotron map shows a largely circular structure, but with interesting deviations in the midplane due to distortion of the field from circularity in the interarm regions., Comment: Accepted for publication in ApJ. Better quality figures in http://www.astro.umd.edu/~gomez/publica/3d_galaxy-3.pdf

The extra stiffness that the magnetic field adds to the ISM changes the way it reacts to the presence of a spiral perturbation. At intermediate to high z, the gas shoots up before the arm, flows over, and falls behind it, as it approaches the next arm. This generates a multicell circulation pattern, within each of which the net radial mass flux is positive near the midplane and negative at higher z. The flow distorts the magnetic field lines. In the arm region, the gas flows nearly parallel to the arm, and therefore, the magnetic field adopts a similar pitch angle. Between the arms, the gas flows out in radius, generating a negative pitch angle in the magnetic field. The intensity and direction of the field yield synthetic synchrotron maps that reproduce some features of the synchrotron maps of external galaxies, like the islands of emission and the displacement between the gaseous and synchrotron arms. When comparing the magnitude of the field with the local gas density, two distinctive relations appear, depending on whether the magnetic pressure is dominant. Above the plane, the density structure develops a shape resembling a breaking wave. This structure collapses and rises again with a period of about 60 Myr, similar to that of a vertical oscillation mode. The falling gas plays an important part in the overall hydrostatics, since its deceleration compresses the low z gas, raising the average midplane pressure in the interarm region above that provided by the weight of the material above., Comment: Accepted for publication in ApJ. Better quality figures in http://www.astro.umd.edu/~gomez/publica/3d_galaxy-2.pdf

We show the initial results of our 3D MHD simulations of the flow of the Galactic atmosphere as it responds to a spiral perturbation in the potential. In our standard case, as the gas approaches the arm, there is a downward converging flow that terminates in a complex of shocks just ahead of the midplane density peak. The density maximum slants forward at high z, preceeded by a similarly leaning shock. The latter diverts the flow upward and over the arm, as in a hydraulic jump. Behind the gaseous arm, the flow falls again, generating further secondary shocks as it approaches the lower z material. Structures similar to the high z part of the gaseous arms are found in the interarm region of our two-armed case, while broken arms and low column density bridges are present in the four-armed case. We present three examples of what can be learned from these models., Comment: 33 pages, 17 figures. Accepted for publication in Apj. Better quality images in http://www.journals.uchicago.edu/ApJ/journal/preprints/ApJ55782.preprint.pdf

When modeling the three-dimensional hydrodynamics of interstellar material rotating in a galactic gravitational potential, it is useful to have an analytic expression for gravitational perturbations due to stellar spiral arms. We present such an expression for which changes in the assumed characteristics of the arms can be made easily and the sensitivity of the hydrodynamics to those characteristics examined. This analytic expression also makes it easy to rotate the force field at the pattern angular velocity with little overhead on the calculations., Comment: 23 pages, 15 figures. Accepted for publication in ApJS

We present a list of 109 pulsars with independent distance information compiled from the literature. Since the compilation of Frail & Weisberg, there are 35 pulsars with new distance estimate and 25 pulsars for which the distance or distance uncertainty have been revised. We used this data to fit a smooth, axisymmetric, two disk model of the distribution of galactic electrons. The two exponential model components have mean local midplane densities at the solar circle of 2.03e-2 cm^-3 and 0.71e-2 cm^-3, and scale heights of 1.07 and 0.053 kpc. The thick component shows very little radial variation, while the second has a radial scale length of only a few kiloparsecs. We also examined a model which varies as sech^2(x), rather than exp(-x), in both the radial and vertical direction. We prefer this model with no midplane cusp, but find that the fit parameters essentially describe the same electron distribution. The distances predicted by this distribution have a similar scatter as the more complex model of Taylor & Cordes. We examine the pulsars that deviate strongly from this model. There are two regions of enhanced dispersion measure, one of which correlates well with the Sagittarius-Carina spiral arm. We find that the scatter of the observed dispersion measure from the model is not fit well by either a normal or log-normal distribution of lump sizes, but may be caused instead by the uncertainties in the distances., Comment: 31 pages, 10 embedded figures, submitted to AJ