Showing 1–50 of 130 results for author: Vázquez-Semadeni, E

The GT and GHC models for molecular clouds compared. Differences, similarities, and myths

Authors: Enrique Vázquez-Semadeni, Aina Palau, Gilberto C. Gómez, Griselda Arroyo-Chávez, Christian Alig, Javier Ballesteros-Paredes, Vianey Camacho, Alejandro González-Samaniego, Andreas Burkert

Abstract: 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… ▽ More 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. △ Less

Submitted 19 August, 2024; originally announced August 2024.

Comments: 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

On the relation between magnetic field strength and gas density in the interstellar medium: A multiscale analysis

Authors: David J. Whitworth, Sundar Srinivasan, Ralph E. Pudritz, Mordecai M. Mac Low, Rowan J. Smith, Aina Palau, Kate Pattle, Gwendoline Eadie, Hector Robinson, Rachel Pillsworth, James Wadsley, Noe Brucy, Ugo Lebreuilly, Patrick Hennebelle, Philipp Girichidis, Fred A. Gent, Jessy Marin, Lylon Sánchez Valido, Vianey Camacho, Ralf S. Klessen, Enrique Vázquez-Semadeni

Abstract: The relation between magnetic field strength B and gas density n in the interstellar medium is of fundamental importance to many areas of astrophysics, from protostellar disks to galaxy evolution. We present and compare Bayesian analyses of the B - n relation for a comprehensive observational data set, as well as a large body of numerical MHD simulations. We extend the original Zeeman relation of… ▽ More The relation between magnetic field strength B and gas density n in the interstellar medium is of fundamental importance to many areas of astrophysics, from protostellar disks to galaxy evolution. We present and compare Bayesian analyses of the B - n relation for a comprehensive observational data set, as well as a large body of numerical MHD simulations. We extend the original Zeeman relation of Crutcher et al. (2010) with a large body of magnetic data that includes 700 observations with the Davis-Chandrasekhar-Fermi method. By using a new multiparameter Bayesian analysis we present a new, more general, time-averaged observational relation: B \propto n^{0.27 \pm 0.017} for n \leq n_0 and B \propto n^{0.54 \pm 0.18} for n \geq n_0, with n_0 = 924^(+145-144) cm^-3. We perform a separate analysis on 19 numerical magnetohydrodynamics simulations that cover a wide range of scales, resolutions, initial conditions, and completed with a variety of codes: arepo, flash, pencil, and ramses. The power law exponents derived from the simulations depend on several physical factors including: dynamo effects, time scales, turbulence, and the initial seed field strength. In particular, early-time simulations where the density, velocity and magnetic fields are unevolved do not match the observational scalings. The simulations that trace the observed density range best, the evolved dwarf galaxy and Milky Way like galaxy simulations, settle into a near consistent exponent of = 0.5 in the dense gas, with variability in the diffuse gas exponent. △ Less

Submitted 25 July, 2024; originally announced July 2024.

Comments: 23 figures, 28 pages, submitted to MNRAS, Comments welcome

Direct observational evidence of multi-epoch massive star formation in G24.47+0.49

Authors: Anindya Saha, Anandmayee Tej, Hong-Li Liu, Tie Liu, Guido Garay, Paul F. Goldsmith, Chang Won Lee, Jinhua He, Mika Juvela, Leonardo Bronfman, Tapas Baug, Enrique Vazquez-Semadeni, Patricio Sanhueza, Shanghuo Li, James O. Chibueze, N. K. Bhadari, Lokesh K. Dewangan, Swagat Ranjan Das, Feng-Wei Xu, Namitha Issac, Jihye Hwang, L. Viktor Toth

Abstract: Using new continuum and molecular line data from the ALMA Three-millimeter Observations of Massive Star-forming Regions (ATOMS) survey and archival VLA, 4.86 GHz data, we present direct observational evidence of hierarchical triggering relating three epochs of massive star formation in a ring-like H II region, G24.47+0.49. We find from radio flux analysis that it is excited by a massive star(s) of… ▽ More Using new continuum and molecular line data from the ALMA Three-millimeter Observations of Massive Star-forming Regions (ATOMS) survey and archival VLA, 4.86 GHz data, we present direct observational evidence of hierarchical triggering relating three epochs of massive star formation in a ring-like H II region, G24.47+0.49. We find from radio flux analysis that it is excited by a massive star(s) of spectral type O8.5V-O8V from the first epoch of star formation. The swept-up ionized ring structure shows evidence of secondary collapse, and within this ring a burst of massive star formation is observed in different evolutionary phases, which constitutes the second epoch. ATOMS spectral line (e.g., HCO$^+$(1-0)) observations reveal an outer concentric molecular gas ring expanding at a velocity of $\sim$ 9 $\rm km\,s^{-1}$, constituting the direct and unambiguous detection of an expanding molecular ring. It harbors twelve dense molecular cores with surface mass density greater than 0.05 $\rm g\,cm^{-2}$, a threshold typical of massive star formation. Half of them are found to be subvirial, and thus in gravitational collapse, making them third epoch of potential massive star-forming sites. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 18 pages, 7 figures, accepted for publication in The Astrophysical Journal Letters

The physical mechanism behind magnetic field alignment in interstellar clouds

Authors: Guido Granda-Muñoz, Enrique Vázquez-Semadeni, Gilberto C. Gómez

Abstract: 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 mag… ▽ More 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énic regime transitions to a transalfvénic regime between the shock and the condensation front, and then to a subalfvénic 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. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 11 pages, 10 figures, 3 animations. Submitted to Astronomy & Astrophysics. Comments are welcome

On the properties and implications of collapse-driven MHD turbulence

Authors: Enrique Vázquez-Semadeni, Yue Hu, Siyao Xu, Rubén Guerrero-Gamboa, Alex Lazarian

Abstract: We numerically investigate the driving of MHD turbulence by gravitational contraction using simulations of an initially spherical, magnetically supercritical cloud core with initially transonic and trans-Alfvénic turbulence. We perform a Helmholtz decomposition of the velocity field, and investigate the evolution of its solenoidal and compressible parts, as well as of the velocity component along… ▽ More We numerically investigate the driving of MHD turbulence by gravitational contraction using simulations of an initially spherical, magnetically supercritical cloud core with initially transonic and trans-Alfvénic turbulence. We perform a Helmholtz decomposition of the velocity field, and investigate the evolution of its solenoidal and compressible parts, as well as of the velocity component along the gravitational acceleration vector, a proxy for the infall component of the velocity field. We find that: 1) In spite of being supercritical, the core first contracts to a sheet perpendicular to the mean field, and the sheet itself collapses. 2) The solenoidal component of the turbulence remains at roughly its initial level throughout the simulation, while the compressible component increases continuously. This implies that turbulence does {\it not} dissipate towards the center of the core. 3) The distribution of simulation cells in the $B$-$ρ$ plane occupies a wide triangular region at low densities, bounded below by the expected trend for fast MHD waves ($B \propto ρ$, applicable for high local Alfvénic Mach number $\Ma$) and above by the trend expected for slow waves ($B \sim$ constant, applicable for low local $\Ma$). At high densities, the distribution follows a single trend $B \propto ρ^{\gamef}$, with $1/2 < \gamef < 2/3$, as expected for gravitational compression. 4) The measured mass-to-magnetic flux ratio $λ$ increases with radius $r$, due to the different scalings of the mass and magnetic flux with $r$. At a fixed radius, $λ$ increases with time due to the accretion of material along field lines. 5) The solenoidal energy fraction is much smaller than the total turbulent component, indicating that the collapse drives the turbulence mainly compressibly, even in directions orthogonal to that of the collapse. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: Resubmitted to MNRAS after first set of reviewer's recommendations. Comments welcome

The ALMA Survey of Star Formation and Evolution in Massive Protoclusters with Blue Profiles (ASSEMBLE): Core Growth, Cluster Contraction, and Primordial Mass Segregation

Authors: Fengwei Xu, Ke Wang, Tie Liu, Mengyao Tang, Neal J. Evans II, Aina Palau, Kaho Morii, Jinhua He, Patricio Sanhueza, Hong-Li Liu, Amelia Stutz, Qizhou Zhang, Xi Chen, Pak Shing Li, Gilberto C. Gómez, Enrique Vázquez-Semadeni, Shanghuo Li, Xiaofeng Mai, Xing Lu, Meizhu Liu, Li Chen, Chuanshou Li, Hongqiong Shi, Zhiyuan Ren, Di Li , et al. (18 additional authors not shown)

Abstract: 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 a… ▽ More 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. △ Less

Submitted 26 September, 2023; originally announced September 2023.

Comments: 37 pages, 13 figures, 5 tables; accepted for publication in ApJS

Multi-scale accretion in dense cloud cores and the delayed formation of massive stars

Authors: Enrique Vázquez-Semadeni, Gilberto C. Gómez, Alejandro González-Samaniego

Abstract: 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… ▽ More 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. △ Less

Submitted 15 February, 2024; v1 submitted 23 June, 2023; originally announced June 2023.

Comments: Resubmitted to MNRAS after incorporating reviewer's recommendations. Comments welcome

Direct observational evidence of the multi-scale, dynamical mass accretion toward a high-mass star forming hub-filament system

Authors: Dongting Yang, Hong-Li Liu, Anandmayee Tej, Tie Liu, Patricio Sanhueza, Sheng-Li Qin, Xing Lu, Ke Wang, Sirong Pan, Feng-Wei Xu, Enrique Vazquez-Semadeni, Shanghuo Li, Gilberto C. Gomez, Aina Palau, Guido Garay, Paul F. Goldsmith, Mika Juvela, Anindya Saha, Leonardo Bronfman, Chang Won Lee, Kenichi Tatematsu, Lokesh Dewangan, Jianwen Zhou, Yong Zhang, Amelia Stutz , et al. (6 additional authors not shown)

Abstract: 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+… ▽ More 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. △ Less

Submitted 17 June, 2023; originally announced June 2023.

Comments: Accepted to publish in ApJ. 10 pages with 6 figures and 2 tables

High-resolution APEX/LAsMA $^{12}$CO and $^{13}$CO (3-2) observation of the G333 giant molecular cloud complex : I. Evidence for gravitational acceleration in hub-filament systems

Authors: J. W. Zhou, F. Wyrowski, S. Neupane, J. S. Urquhart, N. J. Evans II, E. Vázquez-Semadeni, K. M. Menten, Y. Gong, T. Liu

Abstract: Hub-filament systems are suggested to be the birth cradles of high-mass stars and clusters. We apply the FILFINDER algorithm to the integrated intensity maps of the 13CO (3-2) line to identify filaments in the G333 complex, and extract the velocity and intensity along the filament skeleton from moment maps. Clear velocity and density fluctuations are seen along the filaments, allowing us to fit ve… ▽ More Hub-filament systems are suggested to be the birth cradles of high-mass stars and clusters. We apply the FILFINDER algorithm to the integrated intensity maps of the 13CO (3-2) line to identify filaments in the G333 complex, and extract the velocity and intensity along the filament skeleton from moment maps. Clear velocity and density fluctuations are seen along the filaments, allowing us to fit velocity gradients around the intensity peaks. The velocity gradients fitted to the LAsMA data and ALMA data agree with each other over the scales covered by ALMA observations in the ATOMS survey. Changes of velocity gradient with scale indicate a ''funnel'' structure of the velocity field in PPV space, indicative of a smooth, continuously increasing velocity gradient from large to small scales, and thus consistent with gravitational acceleration. The typical velocity gradient corresponding to a 1 pc scale is ~1.6km/s/pc. Assuming free-fall, we estimate a kinematic mass within 1 pc of ~1190 M$_\odot$, which is consistent with typical masses of clumps in the ATLASGAL survey. We find direct evidence for gravitational acceleration from comparison of the observed accelerations to those predicted by free-fall onto dense hubs. On large scales, we find that the inflow may be driven by the larger scale structure, consistent with hierarchical structure in the molecular cloud and gas inflow from large to small scales. The hub-filament structures at different scales may be organized into a hierarchical system extending up to the largest scales probed, through the coupling of gravitational centers at different scales. We argue that the ''funnel'' structure in PPV space can be an effective probe for the gravitational collapse motions in molecular clouds. The large scale gas inflow is driven by gravity, implying that the molecular clouds in G333 complex may be in the state of global gravitational collapse. △ Less

Submitted 21 May, 2023; originally announced May 2023.

Journal ref: A&A 676, A69 (2023)

Simulated observations of star formation regions: infrared evolution of globally collapsing clouds

Authors: Jesús M. Jáquez-Domínguez, Roberto Galván-Madrid, Jacopo Fritz, Manuel Zamora-Avilés, Peter Camps, Gustavo Bruzual, Maarten Baes, Yuxin Lin, Enrique Vázquez-Semadeni

Abstract: The direct comparison between hydrodynamical simulations and observations is needed to improve the physics included in the former and test biases in the latter. Post-processing radiative transfer and synthetic observations are now the standard way to do this. We report on the first application of the \texttt{SKIRT} radiative transfer code to simulations of a star-forming cloud. The synthetic obser… ▽ More The direct comparison between hydrodynamical simulations and observations is needed to improve the physics included in the former and test biases in the latter. Post-processing radiative transfer and synthetic observations are now the standard way to do this. We report on the first application of the \texttt{SKIRT} radiative transfer code to simulations of a star-forming cloud. The synthetic observations are then analyzed following traditional observational workflows. We find that in the early stages of the simulation, stellar radiation is inefficient in heating dust to the temperatures observed in Galactic clouds, thus the addition of an interstellar radiation field is necessary. The spectral energy distribution of the cloud settles rather quickly after $\sim3$ Myr of evolution from the onset of star formation, but its morphology continues to evolve for $\sim8$ Myr due to the expansion of \textsc{Hii} regions and the respective creation of cavities, filaments, and ridges. Modeling synthetic \textit{Herschel} fluxes with 1- or 2-component modified black bodies underestimates total dust masses by a factor of $\sim2$. Spatially-resolved fitting recovers up to about $70\%$ of the intrinsic value. This ``missing mass'' is located in a very cold dust component with temperatures below $10$ K, which does not contribute appreciably to the far-infrared flux. This effect could bias real observations if such dust exists in large amounts. Finally, we tested observational calibrations of the SFR based on infrared fluxes and concluded that they are in agreement when compared to the intrinsic SFR of the simulation averaged over $\sim100$ Myr. △ Less

Submitted 10 April, 2023; originally announced April 2023.

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -XIII. Ongoing triggered star formation within clump-fed scenario found in the massive ($\sim1500$ $\rm M_\odot$) clump

Authors: S. Zhang, K. Wang, T. Liu, A. Zavagno, M. Juvela, H. Liu, A. Tej, A. M. Stutz, S. Li, L. Bronfman, Q. Zhang, P. F. Goldsmith, C. W. Lee, E. Vázquez-Semadeni, K. Tatematsu, W. Jiao, F. Xu, C. Wang, J. -W. Zhou

Abstract: Whether ionization feedback triggers the formation of massive stars is highly debated. Using ALMA 3 mm observations with a spatial resolution of $\sim 0.05$ pc and a mass sensitivity of 1.1 $\rm M_\odot$ beam$^{-1}$ at 20 K, we investigate the star formation and gas flow structures within the ionizing feedback-driven structure, a clump-scale massive ($\gtrsim 1500$ $\rm M_\odot$) bright-rimmed clo… ▽ More Whether ionization feedback triggers the formation of massive stars is highly debated. Using ALMA 3 mm observations with a spatial resolution of $\sim 0.05$ pc and a mass sensitivity of 1.1 $\rm M_\odot$ beam$^{-1}$ at 20 K, we investigate the star formation and gas flow structures within the ionizing feedback-driven structure, a clump-scale massive ($\gtrsim 1500$ $\rm M_\odot$) bright-rimmed cloud (BRC) associated with IRAS 18290-0924. This BRC is bound only if external compression from ionized gas is considered. A small-scale ($\lesssim1$ pc) age sequence along the direction of ionizing radiation is revealed for the embedded cores and protostars, which suggests triggered star formation via radiation-driven implosion (RDI). Furthermore, filamentary gas structures converge towards the cores located in the BRC's center, indicating that these filaments are fueling mass towards cores. The local core-scale mass infall rate derived from H$^{13}$CO$^+$ $J=1-0$ blue profile is of the same order of magnitude as the filamentary mass inflow rate, approximately 1 $\rm M_\odot$ kyr$^{-1}$. A photodissociation region (PDR) covering the irradiated clump surface is detected in several molecules, such as CCH, HCO$^+$, and CS whereas the spatial distribution stratification of these molecules is indistinct. CCH spectra of the PDR possibly indicate a photoevaporation flow leaving the clump surface with a projected velocity of $\sim2$ km s$^{-1}$. Our new observations show that RDI accompanied by a clump-fed process is operating in this massive BRC. Whether this combined process works in other massive BRCs is worth exploring with dedicated surveys. △ Less

Submitted 5 January, 2023; originally announced January 2023.

Comments: 34 pages with 20 figures, Accepted by MNRAS on 2022 December 28

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- XV. Steady Accretion from Global Collapse to Core Feeding in Massive Hub-filament System SDC335

Authors: Feng-Wei Xu, Ke Wang, Tie Liu, Paul F. Goldsmith, Qizhou Zhang, Mika Juvela, Hong-Li Liu, Sheng-Li Qin, Guang-Xing Li, Anandmayee Tej, Guido Garay, Leonardo Bronfman, Shanghuo Li, Yue-Fang Wu, Gilberto C. Gómez, Enrique Vázquez-Semadeni, Ken'ichi Tatematsu, Zhiyuan Ren, Yong Zhang, L. Viktor Toth, Xunchuan Liu, Nannan Yue, Siju Zhang, Tapas Baug, Namitha Issac , et al. (15 additional authors not shown)

Abstract: 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.0… ▽ More 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. △ Less

Submitted 4 January, 2023; originally announced January 2023.

Comments: 29 pages, 7 figures, Accepted for publication in MNRAS

The kinetic and magnetic energy budget of hub-filament systems during the gravitational fragmentation of molecular clouds

Authors: Vianey Camacho, Enrique Vázquez-Semadeni, Manuel Zamora-Avilés, Aina Palau

Abstract: We present a numerical study of the balance between the gravitational (Eg), kinetic (Ek), and magnetic (Em) energies of structures within a hub-filament system in a simulation of the formation and global hierarchical collapse (GHC) of a giant molecular cloud. For structures defined by various density thresholds, and at different evolutionary stages, we investigate the scaling of the virial paramet… ▽ More We present a numerical study of the balance between the gravitational (Eg), kinetic (Ek), and magnetic (Em) energies of structures within a hub-filament system in a simulation of the formation and global hierarchical collapse (GHC) of a giant molecular cloud. For structures defined by various density thresholds, and at different evolutionary stages, we investigate the scaling of the virial parameter, $α$, with mass $M$, and of the Larson ratio, ${\cal{L}}\equivσ_v/R^{1/2}$, with column density $Σ$, where $σ_v$ is the 1D velocity dispersion, and $R$ is an effective radius. We also investigate these scalings for the corresponding magnetic parameters $α_m$ and ${\cal{L}}_m$. Finally, we compare our numerical results with an observational sample of massive clumps. We find that: 1) $α_m$ and ${\cal{L}}_m$ follow similar scalings as their kinetic counterparts, although the ratio Em/Ek decreases as |Eg| increases. 2) The largest objects, defined by the lowest thresholds, tend to appear gravitationally bound (and magnetically supercritical), while their internal substructures tend to appear unbound (and subcritical). This suggests that the latter are being compressed by the infall of their parent structures, and supports earlier suggestions that the measured mass-to-magnetic flux ratio $μ$ decreases inwards in a centrally-peaked cloud under ideal MHD. 3)~The scatter in the $α$-$M$ and ${\cal{L}}$-$Σ$ plots is reduced when Ek and Em are plotted directly against Eg, suggesting that the scatter is due to an ambiguity between mass and size. 4) The clumps in our GHC simulation follow the same trends as the observational sample of massive clumps in the $α$-$M$ and ${\cal{L}}$-$Σ$ diagrams. We conclude that the main controlling parameter of the energy budget in the structures is Eg, with the kinetic and magnetic energies being derived from it. △ Less

Submitted 22 May, 2023; v1 submitted 20 December, 2022; originally announced December 2022.

Comments: Accepted for publication in MNRAS

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- XI. From inflow to infall in hub-filament systems

Authors: Jian-Wen Zhou, Tie Liu, Neal J. Evans II, Guido Garay, Paul F. Goldsmith, Gilberto C. Gomez, Enrique Vazquez-Semadeni, Hong-Li Liu, Amelia M. Stutz, Ke Wang, Mika Juvela, Jinhua He, Di Li, Leonardo Bronfman, Xunchuan Liu, Feng-Wei Xu, Anandmayee Tej, L. K. Dewangan, Shanghuo Li, Siju Zhang, Chao Zhang, Zhiyuan Ren, Kenichi Tatematsu, Pak Shing Li, Chang Won Lee , et al. (15 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 16 June, 2022; originally announced June 2022.

Comments: 16 pages

GASKAP-HI Pilot Survey Science I: ASKAP Zoom Observations of HI Emission in the Small Magellanic Cloud

Authors: N. M. Pingel, J. Dempsey, N. M. McClure-Griffiths, J. M. Dickey, K. E. Jameson, H. Arce, G. Anglada, J. Bland-Hawthorn, S. L. Breen, F. Buckland-Willis, S. E. Clark, J. R. Dawson, H. Dénes, E. M. Di Teodoro, B. -Q. For, Tyler J. Foster, J. F. Gómez, H. Imai, G. Joncas, C. -G. Kim, M. -Y. Lee, C. Lynn, D. Leahy, Y. K. Ma, A. Marchal , et al. (31 additional authors not shown)

Abstract: We present the most sensitive and detailed view of the neutral hydrogen (HI) emission associated with the Small Magellanic Cloud (SMC), through the combination of data from the Australian Square Kilometre Array Pathfinder (ASKAP) and Parkes (Murriyang), as part of the Galactic Australian Square Kilometre Array Pathfinder (GASKAP) pilot survey. These GASKAP-HI pilot observations, for the first time… ▽ More We present the most sensitive and detailed view of the neutral hydrogen (HI) emission associated with the Small Magellanic Cloud (SMC), through the combination of data from the Australian Square Kilometre Array Pathfinder (ASKAP) and Parkes (Murriyang), as part of the Galactic Australian Square Kilometre Array Pathfinder (GASKAP) pilot survey. These GASKAP-HI pilot observations, for the first time, reveal HI in the SMC on similar physical scales as other important tracers of the interstellar medium, such as molecular gas and dust. The resultant image cube possesses an rms noise level of 1.1 K (1.6 mJy/beam) per 0.98 km s$^{-1}$ spectral channel with an angular resolution of 30$''$ ($\sim$10 pc). We discuss the calibration scheme and the custom imaging pipeline that utilizes a joint deconvolution approach, efficiently distributed across a computing cluster, to accurately recover the emission extending across the entire $\sim$25 deg$^2$ field-of-view. We provide an overview of the data products and characterize several aspects including the noise properties as a function of angular resolution and the represented spatial scales by deriving the global transfer function over the full spectral range. A preliminary spatial power spectrum analysis on individual spectral channels reveals that the power-law nature of the density distribution extends down to scales of 10 pc. We highlight the scientific potential of these data by comparing the properties of an outflowing high velocity cloud with previous ASKAP+Parkes HI test observations. △ Less

Submitted 10 December, 2021; v1 submitted 9 November, 2021; originally announced November 2021.

Comments: Accepted for publication in PASA, 34 pages, 18 figures, 5 tables

A resolution criterion based on characteristic timescales for MHD simulations of molecular clouds

Authors: Guido Granda-Muñoz, Enrique Vázquez-Semadeni, Gilberto C. Gómez, Manuel Zamora-Avilés

Abstract: 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… ▽ More 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. △ Less

Submitted 25 January, 2022; v1 submitted 25 October, 2021; originally announced October 2021.

Comments: 7 pages, 5 figures, 2 tables. Accepted to MNRAS

Evolution of the angular momentum during gravitational fragmentation of molecular clouds

Authors: G. Arroyo-Chávez, E. Vázquez-Semadeni

Abstract: We investigate the origin of the observed scaling $j\sim R^{3/2}$ between the specific angular momentum $j$ and the radius $R$ of molecular clouds (MCs) and their their substructures, and of the observed near independence of $β$, the ratio of rotational to gravitational energy, from $R$. To this end, we measure the angular momentum (AM) of sets of particles in an SPH simulation of the formation, c… ▽ More We investigate the origin of the observed scaling $j\sim R^{3/2}$ between the specific angular momentum $j$ and the radius $R$ of molecular clouds (MCs) and their their substructures, and of the observed near independence of $β$, the ratio of rotational to gravitational energy, from $R$. To this end, we measure the angular momentum (AM) of sets of particles in an SPH simulation of the formation, collapse and fragmentation of giant MCs. The sets of SPH particles are defined either as ``clumps'' (connected particle sets), or as lagrangian sets that conform a connected clump only at a certain time $\tdef$. We find that: {\it i)} Clumps evolve along the observed \jR\ relation at all times, {\it ii)} Lagrangian particle sets evolve along the observed relation when the volume containing them also contains a large number of other ``intruder'' particles. Otherwise, they evolve with $j\sim$ cst. {\it iii)} Tracking lagrangian sets to the future, we find that a subset of the SPH particles participates in the collapse, while another disperses away. {\it iv)} Noting that, under AM conservation, $β$ increases during contraction, we suggest that its near independence of radius may arise from the competition between this increase and an increase in the AM exchange rate at higher rotational energy density gradient. {\it v)} We suggest that, if MCs are globally dominated by gravity, the observed \jR\ relation arises because the observational selection of its dense structures amounts to selecting the fragments that have lost AM via Reynolds stresses from their neighbors. △ Less

Submitted 18 June, 2021; originally announced June 2021.

Comments: 20 pages, 15 figures, submitted to ApJ, comments welcome

Filamentary collapse flow in molecular clouds

Authors: Raúl Naranjo-Romero, Enrique Vázquez-Semadeni, Robert M. Loughnane

Abstract: We present idealized numerical simulations of prestellar gravitational collapse of a moderate initial filamentary perturbation with an additional central ellipsoidal enhancement (a core) considering a uniform, and a stratified background, the latter representing flattened clouds. Both simulations maintain the filamentary structure during the collapse, developing a hierarchical accretion flow from… ▽ More We present idealized numerical simulations of prestellar gravitational collapse of a moderate initial filamentary perturbation with an additional central ellipsoidal enhancement (a core) considering a uniform, and a stratified background, the latter representing flattened clouds. Both simulations maintain the filamentary structure during the collapse, developing a hierarchical accretion flow from the cloud to the plane; from there to the filament, and from the filament to the core. The flow changes direction smoothly at every step of the hierarchy, with no density divergence nor a shock developing at the filament's axis during the studied prestellar evolution. The flow drives accretion onto the central core and drains material from the filament, slowing down the growth of the latter. As a consequence, the ratio of the central density of the core to the filament density increases in time, diverging at the time of singularity formation in the core. The stratified simulation produces the best match for observed Plummer-like radial column density profiles of filaments, while the uniform simulation does not produce a flat central density profile. This result supports recent suggestions that MCs may be preferentially flattened structures. We examine the possibility that the filamentary flow might approach a quasi-stationary regime in which the radial accretion onto the filament is balanced by the longitudinal accretion onto the core. A simple argument suggests that such a stationary state may be an attractor for the system. Our simulations, do not attain this stationary stage, but appear to be approaching it during the prestellar stage. △ Less

Submitted 11 January, 2022; v1 submitted 23 December, 2020; originally announced December 2020.

Comments: 15 pages, 16 figures, submitted to MNRAS, category astro-ph.GA Comments are welcome

Density profile evolution during prestellar core collapse: Collapse starts at the large scale

Authors: Gilberto C. Gómez, Enrique Vázquez-Semadeni, Aina Palau

Abstract: 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 $ρ\propto r^{-p}$, we compute the rate of change of the slope $p$, and show that… ▽ More 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 $ρ\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 ρ^β$. For the expected value of $β= 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. △ Less

Submitted 9 February, 2021; v1 submitted 29 September, 2020; originally announced September 2020.

Comments: Accepted for publication in MNRAS

Non-adiabatic turbulence driving during gravitational collapse

Authors: Rubén Guerrero-Gamboa, Enrique Vázquez-Semadeni

Abstract: We investigate the generation of turbulence during the prestellar gravitational contraction of a turbulent spherical core. We define the ratio $g$ of the one-dimensional turbulent velocity dispersion, $σ_\mathrm{1D}$ to the gravitational velocity $v_g$, to then analytically estimate $g$ under the assumptions of a) equipartition or virial equilibrium between the gravitational ($E_g$) and turbulent… ▽ More We investigate the generation of turbulence during the prestellar gravitational contraction of a turbulent spherical core. We define the ratio $g$ of the one-dimensional turbulent velocity dispersion, $σ_\mathrm{1D}$ to the gravitational velocity $v_g$, to then analytically estimate $g$ under the assumptions of a) equipartition or virial equilibrium between the gravitational ($E_g$) and turbulent kinetic ($E_\mathrm{turb}$) energies and b) stationarity of transfer from gravitational to turbulent energy (implying $E_\mathrm{turb}/E_g=$cst). In the equipartition and virial cases, we find $g=\sqrt{1/3}\approx0.58$ and $g=\sqrt{1/6}\approx0.41$, respectively; in the stationary case we find $g=\langle v_\mathrm{rad}\rangle L_d/(4π\sqrt{3}ηRv_g)$, where $η$ is an efficiency factor, $L_d$ is the energy injection scale of the turbulence, and $R$ is the core's radius. Next, we perform AMR simulations of the prestellar collapse of an isothermal, transonic turbulent core at two different resolutions, and a non-turbulent control simulation. We find that the turbulent simulations collapse at the same rate as the non-turbulent one, so that the turbulence generation does not significantly slow down the collapse. We also find that a) the simulations approach near balance between the rates of energy injection from the collapse and of turbulence dissipation; b) $g\approx0.395\pm0.035$, close to the "virial" value (turbulence is $\sim35-40\%$ of non-thermal linewidth); c) the injection scale is $L_d\lesssim R$, and d) the "turbulent pressure" $ρσ_\mathrm{1D}^2$ scales as $\simρ^{1.64}$, an apparently nearly-adiabatic scaling. We propose that this scaling and the nearly virial values of the turbulent velocity dispersion may be reconciled with the non-delayed collapse rate if the turbulence is dissipated as soon as it is generated. △ Less

Submitted 26 September, 2020; v1 submitted 10 September, 2020; originally announced September 2020.

Comments: 14 pages, 6 figures. Accepted for publication in ApJ

From diffuse gas to dense molecular cloud cores

Authors: Javier Ballesteros-Paredes, Philippe André, Patrck Hennebelle, Ralf S. Klessen, Shu-ichiro Inutsuka, J. M. Diederik Kruijssen, Mélanie Chevance, Fumitaka Nakamura, Angela Adamo, Enrique Vázquez-Semadeni

Abstract: Molecular clouds are a fundamental ingredient of galaxies: they are the channels that transform the diffuse gas into stars. The detailed process of how they do it is not completely understood. We review the current knowledge of molecular clouds and their substructure from scales $\sim~$1~kpc down to the filament and core scale. We first review the mechanisms of cloud formation from the warm diffus… ▽ More Molecular clouds are a fundamental ingredient of galaxies: they are the channels that transform the diffuse gas into stars. The detailed process of how they do it is not completely understood. We review the current knowledge of molecular clouds and their substructure from scales $\sim~$1~kpc down to the filament and core scale. We first review the mechanisms of cloud formation from the warm diffuse interstellar medium down to the cold and dense molecular clouds, the process of molecule formation and the role of the thermal and gravitational instabilities. We also discuss the main physical mechanisms through which clouds gather their mass, and note that all of them may have a role at various stages of the process. In order to understand the dynamics of clouds we then give a critical review of the widely used virial theorem, and its relation to the measurable properties of molecular clouds. Since these properties are the tools we have for understanding the dynamical state of clouds, we critically analyse them. We finally discuss the ubiquitous filamentary structure of molecular clouds and its connection to prestellar cores and star formation. △ Less

Submitted 1 June, 2020; originally announced June 2020.

Comments: 73 pages, 11 figures, preprint of a review to appear in Space Science Reviews Topical collection "Star formation". Submitted January 31st 2020; accepted May 17th 2020

Zooming in on Individual Star Formation: Low- and High-mass Stars

Authors: Anna L. Rosen, Stella S. R. Offner, Sarah I. Sadavoy, Asmita Bhandare, Enrique Vázquez-Semadeni, Adam Ginsburg

Abstract: Star formation is a multi-scale, multi-physics problem ranging from the size scale of molecular clouds ($\sim$10s pc) down to the size scales of dense prestellar cores ($\sim$0.1 pc) that are the birth sites of stars. Several physical processes like turbulence, magnetic fields and stellar feedback, such as radiation pressure and outflows, are more or less important for different stellar masses and… ▽ More Star formation is a multi-scale, multi-physics problem ranging from the size scale of molecular clouds ($\sim$10s pc) down to the size scales of dense prestellar cores ($\sim$0.1 pc) that are the birth sites of stars. Several physical processes like turbulence, magnetic fields and stellar feedback, such as radiation pressure and outflows, are more or less important for different stellar masses and size scales. During the last decade a variety of technological and computing advances have transformed our understanding of star formation through the use of multi-wavelength observations, large scale observational surveys, and multi-physics multi-dimensional numerical simulations. Additionally, the use of synthetic observations of simulations have provided a useful tool to interpret observational data and evaluate the importance of various physical processes on different scales in star formation. Here, we review these recent advancements in both high- ($M \gtrsim 8 \, M_{\rm \odot}$) and low-mass star formation. △ Less

Submitted 15 May, 2020; originally announced May 2020.

Comments: 53 pages, 7 figures, To appear in Space Science Reviews (submitted January 31, 2020; accepted May 7, 2020), topical collection on Star formation

The Physics of Star Cluster Formation and Evolution

Authors: Martin G. H. Krause, Stella S. R. Offner, Corinne Charbonnel, Mark Gieles, Ralf S. Klessen, Enrique Vazquez-Semadeni, Javier Ballesteros-Paredes, Philipp Girichidis, J. M. Diederik Kruijssen, Jacob L. Ward, Hans Zinnecker

Abstract: Star clusters form in dense, hierarchically collapsing gas clouds. Bulk kinetic energy is transformed to turbulence with stars forming from cores fed by filaments. In the most compact regions, stellar feedback is least effective in removing the gas and stars may form very efficiently. These are also the regions where, in high-mass clusters, ejecta from some kind of high-mass stars are effectively… ▽ More Star clusters form in dense, hierarchically collapsing gas clouds. Bulk kinetic energy is transformed to turbulence with stars forming from cores fed by filaments. In the most compact regions, stellar feedback is least effective in removing the gas and stars may form very efficiently. These are also the regions where, in high-mass clusters, ejecta from some kind of high-mass stars are effectively captured during the formation phase of some of the low mass stars and effectively channeled into the latter to form multiple populations. Star formation epochs in star clusters are generally set by gas flows that determine the abundance of gas in the cluster. We argue that there is likely only one star formation epoch after which clusters remain essentially clear of gas by cluster winds. Collisional dynamics is important in this phase leading to core collapse, expansion and eventual dispersion of every cluster. We review recent developments in the field with a focus on theoretical work. △ Less

Submitted 15 May, 2020; v1 submitted 2 May, 2020; originally announced May 2020.

Comments: 53 pages, 7 figures, preprint of a review to appear in Space Science Reviews Topical collection "Star formation". Changes in this version: updated comments line to to acceptance of the paper by the journal

Dynamical cloud formation traced by atomic and molecular gas

Authors: H. Beuther, Y. Wang, J. D. Soler, H. Linz, J. Henshaw, E. Vazquez-Semadeni, G. Gomez, S. Ragan, Th. Henning, S. C. O. Glover, M. -Y. Lee, R. Guesten

Abstract: Context: Atomic and molecular cloud formation is a dynamical process. However, kinematic signatures of these processes are still observationally poorly constrained. Methods: Targeting the cloud-scale environment of the prototypical infrared dark cloud G28.3, we employ spectral line imaging observations of the two atomic lines HI and [CI] as well as molecular lines observations in 13CO in the 1--0… ▽ More Context: Atomic and molecular cloud formation is a dynamical process. However, kinematic signatures of these processes are still observationally poorly constrained. Methods: Targeting the cloud-scale environment of the prototypical infrared dark cloud G28.3, we employ spectral line imaging observations of the two atomic lines HI and [CI] as well as molecular lines observations in 13CO in the 1--0 and 3--2 transitions. The analysis comprises investigations of the kinematic properties of the different tracers, estimates of the mass flow rates, velocity structure functions, a Histogram of Oriented Gradients (HOG) study as well as comparisons to simulations. Results: The central IRDC is embedded in a more diffuse envelope of cold neutral medium (CNM) traced by HI self-absorption (HISA) and molecular gas. The spectral line data as well as the HOG and structure function analysis indicate a possible kinematic decoupling of the HI from the other gas compounds. Spectral analysis and position-velocity diagrams reveal two velocity components that converge at the position of the IRDC. Estimated mass flow rates appear rather constant from the cloud edge toward the center. The velocity structure function analysis is consistent with gas flows being dominated by the formation of hierarchical structures. Conclusions: The observations and analysis are consistent with a picture where the IRDC G28 is formed at the center of two converging gas flows. While the approximately constant mass flow rates are consistent with a self-similar, gravitationally driven collapse of the cloud, external compression by, e.g., spiral arm shocks or supernovae explosions cannot be excluded yet. Future investigations should aim at differentiating the origin of such converging gas flows. △ Less

Submitted 16 April, 2020; v1 submitted 14 April, 2020; originally announced April 2020.

Comments: 13 pages, 15 figures, Astronomy & Astrophysics in press, a higher resolution version can be found at http://www.mpia.de/homes/beuther/papers.html

Journal ref: A&A 638, A44 (2020)

The molecular cloud lifecycle

Authors: Mélanie Chevance, J. M. Diederik Kruijssen, Enrique Vazquez-Semadeni, Fumitaka Nakamura, Ralf Klessen, Javier Ballesteros-Paredes, Shu-ichiro Inutsuka, Angela Adamo, Patrick Hennebelle

Abstract: Giant molecular clouds (GMCs) and their stellar offspring are the building blocks of galaxies. The physical characteristics of GMCs and their evolution are tightly connected to galaxy evolution. The macroscopic properties of the interstellar medium propagate into the properties of GMCs condensing out of it, with correlations between e.g. the galactic and GMC scale gas pressures, surface densities… ▽ More Giant molecular clouds (GMCs) and their stellar offspring are the building blocks of galaxies. The physical characteristics of GMCs and their evolution are tightly connected to galaxy evolution. The macroscopic properties of the interstellar medium propagate into the properties of GMCs condensing out of it, with correlations between e.g. the galactic and GMC scale gas pressures, surface densities and volume densities. That way, the galactic environment sets the initial conditions for star formation within GMCs. After the onset of massive star formation, stellar feedback from e.g. photoionisation, stellar winds, and supernovae eventually contributes to dispersing the parent cloud, depositing energy, momentum and metals into the surrounding medium, thereby changing the properties of galaxies. This cycling of matter between gas and stars, governed by star formation and feedback, is therefore a major driver of galaxy evolution. Much of the recent debate has focused on the durations of the various evolutionary phases that constitute this cycle in galaxies, and what these can teach us about the physical mechanisms driving the cycle. We review results from observational, theoretical, and numerical work to build a dynamical picture of the evolutionary lifecycle of GMC evolution, star formation, and feedback in galaxies. △ Less

Submitted 13 April, 2020; originally announced April 2020.

Comments: 52 pages, 7 figures. To appear in Space Science Reviews (submitted February 1, 2020; accepted April 8, 2020), topical collection Star formation

The effect of photoionising feedback on the shaping of hierarchically-forming stellar clusters

Authors: Alejandro González-Samaniego, Enrique Vazquez-Semadeni

Abstract: We use two hydrodynamical simulations (with and without photoionising feedback) of the self-consistent evolution of molecular clouds (MCs) undergoing global hierarchical collapse (GHC), to study the effect of the feedback on the structural and kinematic properties of the gas and the stellar clusters formed in the clouds. During this early stage, the evolution of the two simulations is very similar… ▽ More We use two hydrodynamical simulations (with and without photoionising feedback) of the self-consistent evolution of molecular clouds (MCs) undergoing global hierarchical collapse (GHC), to study the effect of the feedback on the structural and kinematic properties of the gas and the stellar clusters formed in the clouds. During this early stage, the evolution of the two simulations is very similar (implying that the feedback from low mass stars does not affect the cloud-scale evolution significantly) and the star-forming region accretes faster than it can convert gas to stars, causing the instantaneous measured star formation efficiency (SFE) to remain low even in the absence of significant feedback. Afterwards, the ionising feedback first destroys the filamentary supply to star-forming hubs and ultimately removes the gas from it, thus first reducing the star formation (SF) and finally halting it. The ionising feedback also affects the initial kinematics and spatial distribution of the forming stars, because the gas being dispersed continues to form stars, which inherit its motion. In the non-feedback simulation, the groups remain highly compact and do not mix, while in the run with feedback, the gas dispersal causes each group to expand, and the cluster expansion thus consists of the combined expansion of the groups. Most secondary star-forming sites around the main hub are also present in the non-feedback run, implying a primordial rather than triggered nature. We do find one example of a peripheral star-forming site that appears only in the feedback run, thus having a triggered origin. However, this appears to be the exception rather than the rule, although this may be an artifact of our simplified radiative transfer scheme. △ Less

Submitted 23 September, 2020; v1 submitted 28 March, 2020; originally announced March 2020.

Comments: 14 pages, 5 figures, accepted for publication in MNRAS

Simultaneous evolution of the virial parameter and star formation rate in molecular clumps undergoing global hierarchical collapse

Authors: Vianey Camacho, Enrique Vázquez-Semadeni, Aina Palau, Gemma Busquet, Manuel Zamora-Avilés

Abstract: We compare dense clumps and cores in a numerical simulation of molecular clouds (MCs) undergoing global hierarchical collapse (GHC) to observations in two MCs at different evolutionary stages, the Pipe and the G14.225 clouds, to test the ability of the GHC scenario to follow the early evolution of the energy budget and star formation activity of these structures. In the simulation, we select a reg… ▽ More We compare dense clumps and cores in a numerical simulation of molecular clouds (MCs) undergoing global hierarchical collapse (GHC) to observations in two MCs at different evolutionary stages, the Pipe and the G14.225 clouds, to test the ability of the GHC scenario to follow the early evolution of the energy budget and star formation activity of these structures. In the simulation, we select a region that contains cores of sizes and densities similar to the Pipe cores, and find that it evolves through accretion, developing substructure similar to that of G14.225 cloud after $\sim 1.6$ Myr. Within this region, we follow the evolution of the Larson ratio $\mathcal{L} \equiv σ_{\rm v}/R^{1/2}$, where $σ_{\rm v}$ is the velocity dispersion and $R$ is the size, the virial parameter $α$, and the star formation activity of the cores/clumps. In the simulation, we find that as the region evolves: $i)$ its clumps have $\mathcal{L}$ and $α$ values first consistent with those of the Pipe substructures and later with those of G14.225; $ii)$ the individual cores first exhibit a decrease in $α$ followed by an increase when star formation begins; $iii)$ collectively, the ensemble of cores/clumps reproduces the observed trend of lower $α$ for higher-mass objects, and $iv)$ the star formation rate and star formation efficiency increase monotonically. We suggest that this evolution is due to the simultaneous loss of externally-driven compressive kinetic energy and increase of the self-gravity-driven motions. We conclude that the GHC scenario provides a realistic description of the evolution of the energy budget of the clouds' substructure at early times, which occurs simultaneously with an evolution of the star formation activity. △ Less

Submitted 23 September, 2020; v1 submitted 4 February, 2020; originally announced February 2020.

Comments: Accepted for publication in ApJ

Global Hierarchical Collapse In Molecular Clouds. Towards a Comprehensive Scenario

Authors: Enrique Vázquez-Semadeni, Aina Palau, Javier Ballesteros-Paredes, Gilberto C. Gómez, Manuel Zamora-Avilés

Abstract: 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 t… ▽ More 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. △ Less

Submitted 16 October, 2019; v1 submitted 27 March, 2019; originally announced March 2019.

Comments: 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

Structure and Expansion Law of HII Regions in structured Molecular Clouds

Authors: Manuel Zamora-Avilés, Enrique Vázquez-Semadeni, Ricardo F. González, José Franco, Steven N. Shore, Lee W. Hartmann, Javier Ballesteros-Paredes, Robi Banerjee, Bastian Körtgen

Abstract: We present radiation-magnetohydrodynamic simulations aimed at studying evolutionary properties of H\,{\normalsize II} regions in turbulent, magnetised, and collapsing molecular clouds formed by converging flows in the warm neutral medium. We focus on the structure, dynamics and expansion laws of these regions. Once a massive star forms in our highly structured clouds, its ionising radiation eventu… ▽ More We present radiation-magnetohydrodynamic simulations aimed at studying evolutionary properties of H\,{\normalsize II} regions in turbulent, magnetised, and collapsing molecular clouds formed by converging flows in the warm neutral medium. We focus on the structure, dynamics and expansion laws of these regions. Once a massive star forms in our highly structured clouds, its ionising radiation eventually stops the accretion (through filaments) toward the massive star-forming regions. The new over-pressured H\,{\normalsize II} regions push away the dense gas, thus disrupting the more massive collapse centres. Also, because of the complex density structure in the cloud, the H\,{\normalsize II} regions expand in a hybrid manner: they virtually do not expand toward the densest regions (cores), while they expand according to the classical analytical result towards the rest of the cloud, and in an accelerated way, as a blister region, towards the diffuse medium. Thus, the ionised regions grow anisotropically, and the ionising stars generally appear off-centre of the regions. Finally, we find that the hypotheses assumed in standard H\,{\normalsize II}-region expansion models (fully embedded region, blister-type, or expansion in a density gradient) apply simultaneously in different parts of our simulated H\,{\normalsize II} regions, producing a net expansion law ($R \propto t^α$, with $α$ in the range of 0.93-1.47 and a mean value of $1.2 \pm 0.17$) that differs from any of those of the standard models. △ Less

Submitted 17 May, 2019; v1 submitted 4 March, 2019; originally announced March 2019.

Comments: Accepted for publication in MNRAS

Challenges and Techniques for Simulating Line Emission

Authors: Karen P. Olsen, Andrea Pallottini, Aida Wofford, Marios Chatzikos, Mitchell Revalski, Francisco Guzmán, Gergö Popping, Enrique Vázquez-Semadeni, Georgios E. Magdis, Mark L. A. Richardson, Michaela Hirschmann, William J. Gray

Abstract: Modeling emission lines from the millimeter to the UV and producing synthetic spectra is crucial for a good understanding of observations, yet it is an art filled with hazards. This is the proceedings of "Walking the Line", a 3-day conference held in 2018 that brought together scientists working on different aspects of emission line simulations, in order to share knowledge and discuss the methodol… ▽ More Modeling emission lines from the millimeter to the UV and producing synthetic spectra is crucial for a good understanding of observations, yet it is an art filled with hazards. This is the proceedings of "Walking the Line", a 3-day conference held in 2018 that brought together scientists working on different aspects of emission line simulations, in order to share knowledge and discuss the methodology. Emission lines across the spectrum from the millimeter to the UV were discussed, with most of the focus on the interstellar medium, but also some topics on the circumgalactic medium. The most important quality of a useful model is a good synergy with observations and experiments. Challenges in simulating line emission are identified, some of which are already being worked upon, and others that must be addressed in the future for models to agree with observations. Recent advances in several areas aiming at achieving that synergy are summarized here, from micro-physical to galactic and circum-galactic scale. △ Less

Submitted 24 August, 2018; originally announced August 2018.

Comments: 29 pages, 10 figures, 2 tables. Report of the findings at the Walking the Line workshop held in March 2018. Submitted to the MDPI journal Galaxies

Synthetic Spectral Signatures from Isothermal Collapsing Gas and the Interpretation of Infall Profiles

Authors: Robert M. Loughnane, Enrique Vázquez-Semadeni, Raúl Naranjo-Romero

Abstract: We revisit the interpretation of blue-excess molecular lines from dense collapsing cores, considering recent numerical results that suggest prestellar core collapse occurs from the outside-in, and not inside-out. We thus create synthetic molecular-line observations of simulated collapsing, spherically-symmetric, density fluctuations of low initial amplitude, embedded in a uniform, globally gravita… ▽ More We revisit the interpretation of blue-excess molecular lines from dense collapsing cores, considering recent numerical results that suggest prestellar core collapse occurs from the outside-in, and not inside-out. We thus create synthetic molecular-line observations of simulated collapsing, spherically-symmetric, density fluctuations of low initial amplitude, embedded in a uniform, globally gravitationally unstable background, without a turbulent component. The collapsing core develops a flattened, Bonnor-Ebert-like density profile, but with an outside-in radial velocity profile, where the peak infall speeds are at large radii, in lower density gas, with cloud-to-core accretion, and no hydrostatic outer envelope. Using optically thick HCO$^+$ J=1-0 and 3-2 rotational lines, we consider several "typical" beamwidths and use a simple line-fitting model to infer infall speeds from the synthetic profiles similarly to what is done in standard line modeling. We find that the inferred infall speeds are ~ 25-30% of the actual peak infall speed as the largest speeds are downweighted by the low density gas in which they occur, due to the outside-in nature of the actual radial collapse profile. Also, with the N$_2$H$^+$ $J_{F_1F}=1_{01}-0_{12}$ hyperfine line, we investigate the change in the asymmetry parameter, $δv\equiv(V_{thick}-V_{thin})/Δv_{thin}$, during the collapse, finding good agreement with observed values. Finally, the HCO$^+$ J=3-2 lines exhibit extreme $T_b$/$T_r$-ratios like those for evolved cores, for larger beams late in the collapse. Our results suggest that standard dynamical infall reproduces several observed features, but that low-mass core infall speeds are generally undervalued, often interpreted as being subsonic although the actual speeds are supersonic, due to incorrectly assuming an inside-out infall radial velocity profile with a static outer envelope. △ Less

Submitted 26 January, 2022; v1 submitted 14 August, 2018; originally announced August 2018.

Comments: The velocities in the radiative transfer analysis for each snapshot of the simulated collapsing core are a factor of 1.732 (the square root of 3) smaller along each of the three dimensions. This error is currently being rectified

Molecular Cloud Evolution VI. Measuring cloud ages

Authors: Enrique Vázquez-Semadeni, Manuel Zamora-Avilés, Roberto Galván-Madrid, Jan Forbrich

Abstract: In previous contributions, we have presented an analytical model describing the evolution of molecular clouds (MCs) undergoing hierarchical gravitational contraction. The cloud's evolution is characterized by an initial increase in its mass, density, and star formation rate (SFR) and efficiency (SFE) as it contracts, followed by a decrease of these quantities as newly formed massive stars begin to… ▽ More In previous contributions, we have presented an analytical model describing the evolution of molecular clouds (MCs) undergoing hierarchical gravitational contraction. The cloud's evolution is characterized by an initial increase in its mass, density, and star formation rate (SFR) and efficiency (SFE) as it contracts, followed by a decrease of these quantities as newly formed massive stars begin to disrupt the cloud. The main parameter of the model is the maximum mass reached by the cloud during its evolution. Thus, specifying the instantaneous mass and some other variable completely determines the cloud's evolutionary stage. We apply the model to interpret the observed scatter in SFEs of the cloud sample compiled by Lada et al.\ as an evolutionary effect so that, although clouds such as California and Orion A have similar masses, they are in very different evolutionary stages, causing their very different observed SFRs and SFEs. The model predicts that the California cloud will eventually reach a significantly larger total mass than the Orion A cloud. Next, we apply the model to derive estimated ages of the clouds since the time when approximately 25\% of their mass had become molecular. We find ages from $\sim 1.5$ to 27 Myr, with the most inactive clouds being the youngest. Further predictions of the model are that clouds with very low SFEs should have massive atomic envelopes constituting the majority of their gravitational mass, and that low-mass clouds ($M \sim 10^3$-$10^4 \, M_\odot$) end their lives with a mini-burst of star formation, reaching SFRs $\sim 300$-$500\, M_\odot$ Myr$^{-1}$. By this time, they have contracted to become compact ($\sim 1$ pc) massive star-forming clumps, in general embedded within larger GMCs. △ Less

Submitted 16 July, 2019; v1 submitted 17 May, 2018; originally announced May 2018.

Comments: 11 pages, 3 figures, and 2 tables. Published in MNRAS. Journal-ref: MNRAS, 2018, Volume 479, p.3254-3263. ADS link: https://ui.adsabs.harvard.edu/abs/2018MNRAS.479.3254V/abstract

The magnetic field structure in molecular cloud filaments

Authors: Gilberto C. Gomez, Enrique Vazquez-Semadeni, Manuel Zamora-Aviles

Abstract: 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 othe… ▽ More 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. △ Less

Submitted 9 January, 2018; originally announced January 2018.

Comments: Submitted to MNRAS

Gravity or turbulence? IV. Collapsing cores in out-of-virial disguise

Authors: Javier Ballesteros-Paredes, Enrique Vázquez-Semadeni, Aina Palau, Ralf S. Klessen

Abstract: We study the dynamical state of cores by using a simple analytical model, a sample of observational massive cores, and numerical simulations of collapsing massive cores. From the model, we find that, if cores are formed from turbulent compressions, they evolve from small to large column densities, increasing their velocity dispersion as they collapse, while they tend to equipartition between gravi… ▽ More We study the dynamical state of cores by using a simple analytical model, a sample of observational massive cores, and numerical simulations of collapsing massive cores. From the model, we find that, if cores are formed from turbulent compressions, they evolve from small to large column densities, increasing their velocity dispersion as they collapse, while they tend to equipartition between gravity and kinetic energy. From the observational sample, we find that: (a) cores with substantially different column densities in the sample do not follow a Larson-like linewidth-size relation. Instead, cores with higher column densities tend to be located in the upper-left corner of the Larson velocity dispersion-size diagram, a result predicted previously (Ballesteros-Paredes et al. 2011a). (b) The data exhibit cores with overvirial values. Finally, in the simulations we reproduce the behavior depicted by the model and observational sample: cores evolve towards larger velocity dispersions and smaller sizes as they collapse and increase their column density. However, collapsing cores appear to approach overvirial states within a free-fall time. The cause of this apparent excess of kinetic energy is an underestimation of the gravitational energy, due to the assumption that the gravitational energy is given by the energy of an isolated sphere of constant column density. This excess disappears when the gravitational energy is correctly calculated from the actual spatial mass distribution, where inhomogeneities, as well as the potential due to the mass outside of the core, also contribute to the gravitational energy. We conclude that the observed energy budget of cores in surveys is consistent with their non-thermal motions being driven by their self-gravity and in a hierarchical and chaotic collapse. △ Less

Submitted 19 October, 2017; originally announced October 2017.

Comments: 15 pages, 4 figures. Submitted to MNRAS. Comments welcome

Partially Ionized Plasmas in Astrophysics

Authors: Jose Luis Ballester, Igor Alexeev, Manuel Collados, Turlough Downes, Robert F. Pfaff, Holly Gilbert, Maxim Khodachenko, Elena Khomenko, Ildar F. Shaikhislamov, Roberto Soler, Enrique Vazquez-Semadeni, Teimuraz Zaqarashvili

Abstract: Partially ionized plasmas are found across the Universe in many different astrophysical environments. They constitute an essential ingredient of the solar atmosphere, molecular clouds, planetary ionospheres and protoplanetary disks, among other environments, and display a richness of physical effects which are not present in fully ionized plasmas. This review provides an overview of the physics of… ▽ More Partially ionized plasmas are found across the Universe in many different astrophysical environments. They constitute an essential ingredient of the solar atmosphere, molecular clouds, planetary ionospheres and protoplanetary disks, among other environments, and display a richness of physical effects which are not present in fully ionized plasmas. This review provides an overview of the physics of partially ionized plasmas, including recent advances in different astrophysical areas in which partial ionization plays a fundamental role. We outline outstanding observational and theoretical questions and discuss possible directions for future progress. △ Less

Submitted 26 July, 2017; v1 submitted 25 July, 2017; originally announced July 2017.

Comments: 163 pages, submitted to Space Science Reviews

Hierarchical Cluster Assembly in Globally Collapsing Clouds

Authors: Enrique Vazquez-Semadeni, Alejandro Gonzalez-Samaniego, Pedro Colin

Abstract: We discuss the mechanism of cluster formation in a numerical simulation of a molecular cloud (MC) undergoing global hierarchical collapse (GHC). The global nature of the collapse implies that the SFR increases over time. The hierarchical nature of the collapse consists of small-scale collapses within larger-scale ones. The large-scale collapses culminate a few Myr later than the small-scale ones a… ▽ More We discuss the mechanism of cluster formation in a numerical simulation of a molecular cloud (MC) undergoing global hierarchical collapse (GHC). The global nature of the collapse implies that the SFR increases over time. The hierarchical nature of the collapse consists of small-scale collapses within larger-scale ones. The large-scale collapses culminate a few Myr later than the small-scale ones and consist of filamentary flows that accrete onto massive central clumps. The small-scale collapses form clumps that are embedded in the filaments and falling onto the large-scale collapse centers. The stars formed in the early, small-scale collapses share the infall motion of their parent clumps. Thus, the filaments feed both gaseous and stellar material to the massive central clump. This leads to the presence of a few older stars in a region where new protostars are forming, and also to a self-similar structure, in which each unit is composed of smaller-scale sub-units that approach each other and may merge. Because the older stars formed in the filaments share the infall motion of the gas onto the central clump, they tend to have larger velocities and to be distributed over larger areas than the younger stars formed in the central clump. Finally, interpreting the IMF at face-value as a probability distribution implies that massive stars only form once the {\it local} SFR is large enough to sample the IMF up to high masses. In combination with the increase of the SFR, this implies that massive stars tend to appear late in the evolution of the MC, and only in the central massive clumps. We discuss the correspondence of these features with observed properties of young stellar clusters, finding very good qualitative agreement, thus providing support to the scenario of global, hierarchical collapse of MCs, while explaining the origin of the observed cluster structure. △ Less

Submitted 31 October, 2016; originally announced November 2016.

Comments: 17 pages, 12 figures. Submitted to MNRAS. Comments welcome

Energy budget of forming clumps in numerical simulations of collapsing clouds

Authors: Vianey Camacho, Enrique Vázquez-Semadeni, Javier Ballesteros-Paredes, Gilberto C. Gómez, S. Michael Fall, M. Dolores Mata-Chávez

Abstract: 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 resul… ▽ More 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 $σ_{v}/R^{1/2} \propto Σ^{1/2}$, where $σ_{v}$ is the velocity dispersion, $R$ is the "radius", and $Σ$ 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 $Σ$ 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 $Σ$ ranges follow Larson-like relations, although characterized by their respective value of $Σ$. 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. △ Less

Submitted 3 October, 2016; v1 submitted 28 September, 2016; originally announced September 2016.

Comments: Abstract rephrased. Submitted to ApJ

Supernova Feedback in Molecular Clouds: Global Evolution and Dynamics

Authors: Bastian Körtgen, Daniel Seifried, Robi Banerjee, Enrique Vázquez-Semadeni, Manuel Zamora-Avilés

Abstract: We use magnetohydrodynamical simulations of converging warm neutral medium flows to analyse the formation and global evolution of magnetised and turbulent molecular clouds subject to supernova feedback from massive stars. We show that supernova feedback alone fails to disrupt entire, gravitationally bound, molecular clouds, but is able to disperse small--sized (~10 pc) regions on timescales of les… ▽ More We use magnetohydrodynamical simulations of converging warm neutral medium flows to analyse the formation and global evolution of magnetised and turbulent molecular clouds subject to supernova feedback from massive stars. We show that supernova feedback alone fails to disrupt entire, gravitationally bound, molecular clouds, but is able to disperse small--sized (~10 pc) regions on timescales of less than 1 Myr. Efficient radiative cooling of the supernova remnant as well as strong compression of the surrounding gas result in non-persistent energy and momentum input from the supernovae. However, if the time between subsequent supernovae is short and they are clustered, large hot bubbles form that disperse larger regions of the parental cloud. On longer timescales, supernova feedback increases the amount of gas with moderate temperatures (T~300-3000 K). Despite its inability to disrupt molecular clouds, supernova feedback leaves a strong imprint on the star formation process. We find an overall reduction of the star formation efficiency by a factor of 2 and of the star formation rate by roughly factors of 2-4. △ Less

Submitted 31 March, 2016; originally announced March 2016.

Comments: 16 pages, 12 figures (2 in appendix), revised version, submitted to MNRAS

Cluster Assembly in Hierarchically Collapsing Clouds

Authors: Enrique Vazquez-Semadeni, Alejandro Gonzalez-Samaniego, Manuel Zamora-Aviles, Pedro Colin

Abstract: We discuss the mechanism of cluster formation in hierarchically collapsing molecular clouds. Recent evidence, both observational and numerical, suggests that molecular clouds (MCs) may be undergoing global, hierarchical gravitational collapse. The "hierarchical" regime consists of small-scale collapses within larger-scale ones. The latter implies that the star formation rate increases systematical… ▽ More We discuss the mechanism of cluster formation in hierarchically collapsing molecular clouds. Recent evidence, both observational and numerical, suggests that molecular clouds (MCs) may be undergoing global, hierarchical gravitational collapse. The "hierarchical" regime consists of small-scale collapses within larger-scale ones. The latter implies that the star formation rate increases systematically during the early stages of evolution, and occurs via filamentary flows onto "hubs" of higher density, mass, and velocity dispersion, and culminates a few Myr after than the small-scale collapses have started to form stars. In turn, the small-scale collapses occur in clumps embedded in the filaments, and are themselves falling into the larger potential well of the still-ongoing large-scale collapse. The stars formed in the early, small-scale collapses share the infall motion of their parent clumps towards the larger potential trough, so that the filaments feed both gaseous and stellar material to the hubs. This leads to the presence of older stars in a region where new protostars are still forming, to a scale-free or fractal structure of the clusters, in which each unit is composed of smaller-scale ones, and to the eventual merging of the subunits, explaining the observed structural features of open clusters. △ Less

Submitted 23 February, 2016; originally announced February 2016.

Comments: 6 pages, 7 figures, to appear in the proceedings of IAU Symposium 316: Formation, Evolution, and Survival of Massive Star Clusters

Hierarchical gravitational fragmentation. I. Collapsing cores within collapsing clouds

Authors: Raúl Naranjo-Romero, Enrique Vázquez-Semadeni, Robert M. Loughnane

Abstract: We investigate the Hierarchical Gravitational Fragmentation scenario through numerical simulations of the prestellar stages of the collapse of a marginally gravitationally unstable isothermal sphere immersed in a strongly gravitationally unstable, uniform background medium. The core developes a Bonnor-Ebert (BE)-like density profile, while at the time of singularity (the protostar) formation the e… ▽ More We investigate the Hierarchical Gravitational Fragmentation scenario through numerical simulations of the prestellar stages of the collapse of a marginally gravitationally unstable isothermal sphere immersed in a strongly gravitationally unstable, uniform background medium. The core developes a Bonnor-Ebert (BE)-like density profile, while at the time of singularity (the protostar) formation the envelope approaches a singular-isothermal-sphere (SIS)-like $r^-2$ density profile. However, these structures are never hydrostatic. In this case, the central flat region is characterized by an infall speed, while the envelope is characterized by a uniform speed. This implies that the hydrostatic SIS initial condition leading to Shu's classical inside-out solution is not expected to occur, and therefore neither should the inside-out solution. Instead, the solution collapses from the outside-in, naturally explaining the observation of extended infall velocities. The core, defined by the radius at which it merges with the background, has a time-variable mass, and evolves along the locus of the ensemble of observed prestellar cores in a plot of $M/M_{BE}$ vs. $M$, where $M$ is the core's mass and $M_{BE}$ is the critical Bonnor-Ebert mass, spanning the range from the "stable" to the "unstable" regimes, even though it is collapsing at all times. We conclude that the presence of an unstable background allows a core to evolve dynamically from the time when it first appears, even when it resembles a pressure-confined, stable BE-sphere. The core can be thought of as a ram-pressure confined BE-sphere, with an increasing mass due to the accretion from the unstable background. △ Less

Submitted 19 October, 2015; originally announced October 2015.

Comments: 24 pages, 3 figures. Accepted for publication in the Astrophysical Journal

Gravity or turbulence? -III. Evidence of pure thermal Jeans fragmentation at ~0.1 pc scale

Authors: Aina Palau, Javier Ballesteros-Paredes, Enrique Vazquez-Semadeni, Alvaro Sanchez-Monge, Robert Estalella, S. Michael Fall, Luis A. Zapata, Vianey Camacho, Laura Gomez, Raul Naranjo-Romero, Gemma Busquet, Francesco Fontani

Abstract: We combine previously published interferometric and single-dish data of relatively nearby massive dense cores that are actively forming stars to test whether their `fragmentation level' is controlled by turbulent or thermal support. We find no clear correlation between the fragmentation level and velocity dispersion, nor between the observed number of fragments and the number of fragments expected… ▽ More We combine previously published interferometric and single-dish data of relatively nearby massive dense cores that are actively forming stars to test whether their `fragmentation level' is controlled by turbulent or thermal support. We find no clear correlation between the fragmentation level and velocity dispersion, nor between the observed number of fragments and the number of fragments expected when the gravitationally unstable mass is calculated including various prescriptions for `turbulent support'. On the other hand, the best correlation is found for the case of pure thermal Jeans fragmentation, for which we infer a core formation efficiency around 13 per cent, consistent with previous works. We conclude that the dominant factor determining the fragmentation level of star-forming massive dense cores at 0.1 pc scale seems to be thermal Jeans fragmentation. △ Less

Submitted 13 October, 2015; v1 submitted 28 April, 2015; originally announced April 2015.

Comments: accepted in MNRAS

Galactic and Magellanic Evolution with the SKA

Authors: Naomi M. McClure-Griffiths, Snezana Stanimirovic, Claire E. Murray, Di Li, John M. Dickey, Enrique Vazquez-Semadeni, Josh E. G. Peek, Mary Putman, Susan E. Clark, Marc-Antoine Miville-Deschenes, Joss Bland-Hawthorn, Lister Staveley-Smith

Abstract: As we strive to understand how galaxies evolve it is crucial that we resolve physical processes and test emerging theories in nearby systems that we can observe in great detail. Our own Galaxy, the Milky Way, and the nearby Magellanic Clouds provide unique windows into the evolution of galaxies, each with its own metallicity and star formation rate. These laboratories allow us to study with more d… ▽ More As we strive to understand how galaxies evolve it is crucial that we resolve physical processes and test emerging theories in nearby systems that we can observe in great detail. Our own Galaxy, the Milky Way, and the nearby Magellanic Clouds provide unique windows into the evolution of galaxies, each with its own metallicity and star formation rate. These laboratories allow us to study with more detail than anywhere else in the Universe how galaxies acquire fresh gas to fuel their continuing star formation, how they exchange gas with the surrounding intergalactic medium, and turn warm, diffuse gas into molecular clouds and ultimately stars. The $λ$21-cm line of atomic hydrogen (HI) is an excellent tracer of these physical processes. With the SKA we will finally have the combination of surface brightness sensitivity, point source sensitivity and angular resolution to transform our understanding of the evolution of gas in the Milky Way, all the way from the halo down to the formation of individual molecular clouds. △ Less

Submitted 6 January, 2015; originally announced January 2015.

Comments: 25 pages, from "Advancing Astrophysics with the Square Kilometre Array", to appear in Proceedings of Science

Inverse Hubble Flows in Molecular Clouds

Authors: Jesús A. Toalá, Enrique Vázquez-Semadeni, Pedro Colín, Gilberto C. Gómez

Abstract: 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… ▽ More 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 $τ_\mathrm{nl}$ for a density fluctuation to become nonlinear increases without limit as its initial value approaches zero, while in the IHF case $τ_\mathrm{nl} \le τ_\mathrm{ff}$ always, where $τ_\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. △ Less

Submitted 22 January, 2015; v1 submitted 6 November, 2014; originally announced November 2014.

Comments: 7 pages, 3 figures; published in MNRAS

The identification of filaments on far infrared and submillimiter images. Morphology, physical conditions and relation with star formation of filamentary structure

Authors: E. Schisano, K. L. Rygl, S. Molinari, G. Busquet, D. Elia, M. Pestalozzi, D. Polychroni, N. Billot, A. Noriega-Crespo, S. Carey, R. Paladini, T. J. T. Moore, R. Plume, S. C. O. Glover, E. Vazquez-Semadeni

Abstract: Observations of molecular clouds reveal a complex structure, with gas and dust often arranged in filamentary rather than spherical geometries. The associations of pre- and proto- stellar cores with the filaments suggest a direct link with the process of star formation. Any study of the properties of such filaments requires a representative samples from different enviroments and so an unbiased dete… ▽ More Observations of molecular clouds reveal a complex structure, with gas and dust often arranged in filamentary rather than spherical geometries. The associations of pre- and proto- stellar cores with the filaments suggest a direct link with the process of star formation. Any study of the properties of such filaments requires a representative samples from different enviroments and so an unbiased detection method. We developed such an approach using the Hessian matrix of a surface-brightness distribution to identify filaments and determine their physical and morphological properties. After testing the method on simulated, but realistic filaments, we apply the algorithms to column-density maps computed from Herschel observations of the Galactic Plane obtained by the Hi-GAL project. We identified ~500 filaments, in the longitude range of l=216.5 to l=225.5, with lengths from ~1 pc up to ~30 pc and widths between 0.1 pc and 2.5 pc. Average column densities are between 10^20 cm^-2 and 10^22 cm^-2. Filaments include the majority of dense material with N_H_2 > 6x10^21cm^-2. We find that the pre- and proto-stellar compact sources already identified in the same region are mostly associated with filaments. However, surface densities in excess of the expected critical values for high-mass star formation are only found on the filaments, indicating that these structures are necessary to channel material into the clumps. Furthermore, we analyze the gravitational stability of filaments and discuss their relationship with star formation. △ Less

Submitted 17 June, 2014; originally announced June 2014.

Comments: 32 pages, 23 Figures. Accepted for publication in ApJ

Molecular cloud formation as seen in synthetic Hi and molecular gas observations

Authors: Jonathan S. Heiner, Enrique Vázquez-Semadeni, Javier Ballesteros-Paredes

Abstract: We present synthetic Hi and CO observations of a simulation of decaying turbulence in the thermally bistable neutral medium. We first present the simulation, with clouds initially consisting of clustered clumps. Self-gravity causes these clump clusters to form more homogeneous dense clouds. We apply a simple radiative transfer algorithm, and defining every cell with <Av> > 1 as molecular. We then… ▽ More We present synthetic Hi and CO observations of a simulation of decaying turbulence in the thermally bistable neutral medium. We first present the simulation, with clouds initially consisting of clustered clumps. Self-gravity causes these clump clusters to form more homogeneous dense clouds. We apply a simple radiative transfer algorithm, and defining every cell with <Av> > 1 as molecular. We then produce maps of Hi, CO-free molecular gas, and CO, and investigate the following aspects: i) The spatial distribution of the warm, cold, and molecular gas, finding the well-known layered structure, with molecular gas surrounded by cold Hi, surrounded by warm Hi. ii) The velocity of the various components, with atomic gas generally flowing towards the molecular gas, and that this motion is reflected in the frequently observed bimodal shape of the Hi profiles. This conclusion is tentative, because we do not include feedback. iii) The production of Hi self-absorption (HISA) profiles, and the correlation of HISA with molecular gas. We test the suggestion of using the second derivative of the brightness temperature Hi profile to trace HISA and molecular gas, finding limitations. On a scale of ~parsecs, some agreement is obtained between this technique and actual HISA, as well as a correlation between HISA and N(mol). It quickly deteriorates towards sub-parsec scales. iv) The N-PDFs of the actual Hi gas and those recovered from the Hi line profiles, with the latter having a cutoff at column densities where the gas becomes optically thick, thus missing the contribution from the HISA-producing gas. We find that the power-law tail typical of gravitational contraction is only observed in the molecular gas, and that, before the power-law tail develops in the total gas density PDF, no CO is yet present, reinforcing the notion that gravitational contraction is needed to produce this component. (abridged) △ Less

Submitted 25 March, 2014; originally announced March 2014.

Comments: Submitted to MNRAS - comments are welcome. 22 pages, 21 figures

The Milky Way as a Star Formation Engine

Authors: Sergio Molinari, John Bally, Simon Glover, Toby Moore, Alberto Noriega-Crespo, René Plume, Leonardo Testi, Enrique Vázquez-Semadeni, Annie Zavagno, Jean-Philippe Bernard, Peter Martin

Abstract: The cycling of material from the interstellar medium (ISM) into stars and the return of stellar ejecta into the ISM is the engine that drives the "galactic ecology" in normal spirals, a cornerstone in the formation and evolution of galaxies through cosmic time. Major observational and theoretical challenges need to be addressed in determining the processes responsible for converting the low-densit… ▽ More The cycling of material from the interstellar medium (ISM) into stars and the return of stellar ejecta into the ISM is the engine that drives the "galactic ecology" in normal spirals, a cornerstone in the formation and evolution of galaxies through cosmic time. Major observational and theoretical challenges need to be addressed in determining the processes responsible for converting the low-density ISM into dense molecular clouds, forming dense filaments and clumps, fragmenting them into stars, OB associations and bound clusters, and characterizing the feedback that limits the rate and efficiency of star formation. This formidable task can be now effectively attacked thanks to the combination of new global-scale surveys of the Milky Way Galactic Plane from infrared to radio wavelengths, offering the possibility of bridging the gap between local and extragalactic star formation studies. The Herschel, Spitzer and WISE mid to far infrared continuum surveys, complemented by analogue surveys from ground-based facilities in the millimetre and radio wavelengths, enables us to measure the Galactic distribution and physical properties of dust on all scales and in all components of the ISM from diffuse clouds to filamentary complexes and tens of thousands of dense clumps. A complementary suite of spectroscopic surveys in various atomic and molecular tracers is providing the chemical fingerprinting of dense clumps and filaments, as well as essential kinematic information to derive distances and thus transform panoramic data into a 3D representation. The latest results emerging from these Galaxy-scale surveys are reviewed. New insights into cloud formation and evolution, filaments and their relationship to channeling gas onto gravitationally-bound clumps, the properties of these clumps, density thresholds for gravitational collapse, and star and cluster formation rates are discussed. △ Less

Submitted 16 May, 2014; v1 submitted 25 February, 2014; originally announced February 2014.

Comments: Accepted for publication as a chapter in Protostars and Planets VI, University of Arizona Press (2014), eds. H. Beuther, R. Klessen, C. Dullemond, Th. Henning (minor text revisions, plus updated figure 11)

Formation of Molecular Clouds and Global Conditions for Star Formation

Authors: Clare L. Dobbs, Mark R. Krumholz, Javier Ballesteros-Paredes, Alberto D. Bolatto, Yasuo Fukui, Mark Heyer, Mordecai-Mark Mac Low, Eve C. Ostriker, Enrique Vázquez-Semadeni

Abstract: Giant molecular clouds (GMCs) are the primary reservoirs of cold, star-forming molecular gas in the Milky Way and similar galaxies, and thus any understanding of star formation must encompass a model for GMC formation, evolution, and destruction. These models are necessarily constrained by measurements of interstellar molecular and atomic gas, and the emergent, newborn stars. Both observations and… ▽ More Giant molecular clouds (GMCs) are the primary reservoirs of cold, star-forming molecular gas in the Milky Way and similar galaxies, and thus any understanding of star formation must encompass a model for GMC formation, evolution, and destruction. These models are necessarily constrained by measurements of interstellar molecular and atomic gas, and the emergent, newborn stars. Both observations and theory have undergone great advances in recent years, the latter driven largely by improved numerical simulations, and the former by the advent of large-scale surveys with new telescopes and instruments. This chapter offers a thorough review of the current state of the field. △ Less

Submitted 11 December, 2013; originally announced December 2013.

Comments: 25 pages, 7 figures, Accepted for publication as a chapter in Protostars and Planets VI, University of Arizona Press (2014), eds. H. Beuther, R. Klessen, C. Dullemond, Th. Henning

Testing assumptions and predictions of star-formation theories

Authors: Alejandro González-Samaniego, Enrique Vázquez-Semadeni, Ricardo F. González, Jongsoo Kim

Abstract: (Abridged). We present numerical simulations of isothermal, MHD, supersonic turbulence, designed to test various hypotheses frequently assumed in star formation(SF) theories. We consider three simulations, each with a different combination of physical size, rms sonic Mach number, and Jeans parameter, but chosen as to give the same value of the virial parameter and to conform with Larson's scaling… ▽ More (Abridged). We present numerical simulations of isothermal, MHD, supersonic turbulence, designed to test various hypotheses frequently assumed in star formation(SF) theories. We consider three simulations, each with a different combination of physical size, rms sonic Mach number, and Jeans parameter, but chosen as to give the same value of the virial parameter and to conform with Larson's scaling relations. As in the non-magnetic case: we find no simultaneously subsonic and super-Jeans structures in our MHD simulations. We find that the fraction of small-scale super-Jeans structures increases when self gravity is turned on, and that the production of gravitationally unstable dense cores by turbulence alone is very low. This implies that self-gravity is in general necessary not only to induce the collapse of Jeans-unstable cores, but also to form them. We find that denser regions tend to have more negative values of the velocity divergence, implying a net inwards flow towards the regions' centers. We compare the results from our simulations with the predictions from the recent SF theories by Krumholz & McKee, Padoan & Nordlund, and Hennebelle & Chabrier, using the expressions recently provided by Federrath & Klessen. We find that none of these theories reproduces the dependence of the SFEff with Ms observed in our simulations in the MHD case. The SFEff predicted by the theories ranges between half and one order of magnitude larger than what we observe in the simulations in both the HD and the MHD cases. We conclude that the type of flow used in simulations like the ones presented here and assumed in recent SF theories, may not correctly represent the flow within actual clouds, and that theories that assume it does may be missing a fundamental aspect of the flow. We suggest that a more realistic regime may be that of hierarchical gravitational collapse. △ Less

Submitted 28 March, 2014; v1 submitted 23 October, 2013; originally announced October 2013.

Comments: 20 pages, 13 figures, minor changes after final referee report, results unchanged; added references. Accepted for publication in MNRAS

Filaments in Simulations of Molecular Cloud Formation

Authors: Gilberto C. Gomez, Enrique Vazquez-Semadeni

Abstract: 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… ▽ More 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. △ Less

Submitted 11 July, 2014; v1 submitted 28 August, 2013; originally announced August 2013.

Comments: Revised to address the referee's comments, submitted to ApJ. See related animations in http://www.crya.unam.mx/~g.gomez/publica.html

An Evolutionary Model for Collapsing Molecular Clouds and Their Star Formation Activity. II. Mass Dependence of the Star Formation Rate

Authors: Manuel Zamora-Avilés, Enrique Vázquez-Semadeni

Abstract: We discuss the dependence of various properties of the star formation rate (SFR) and efficiency (SFE) in molecular clouds (MCs) on the maximum mass reached by the clouds, based on a previously-published model for MC and SFR evolution in which the clouds were assumed to be undergoing global collapse, and the SFR was controlled by ioniztion feedback. Because the model neglects various other processe… ▽ More We discuss the dependence of various properties of the star formation rate (SFR) and efficiency (SFE) in molecular clouds (MCs) on the maximum mass reached by the clouds, based on a previously-published model for MC and SFR evolution in which the clouds were assumed to be undergoing global collapse, and the SFR was controlled by ioniztion feedback. Because the model neglects various other processes, the results presented are upper limits. We find that clouds with $\Mmax \lesssim 10^4 \Msun$ end their lives with a mini-burst, at which the SFR reaches a peak of $\sim 10^4 ~\Msun \Myr^{-1}$, although its time average is only $\SFRavg \sim \hbox{ a few} \times 10^2 \Msun \Myr^{-1}$. The corresponding efficiencies are $\SFEmax \lesssim $60\%$ and $\SFEavg \lesssim $1\%$. For more massive clouds ($\Mmax \gtrsim 10^5 ~ \Msun$), the SFR first increases and then remains roughly constant for $\sim 10^7$ yr, because the clouds are influenced by the stellar feedback since earlier in their evolution. We find that $\SFRavg$ and $\SFEavg$ are well represented by the fits $\SFRavg \approx 100 (1+\Mmax/2 \times 10^5 ~ \Msun)^{2} ~ \Msun \Myr^{-1}$ and $\SFEavg \approx 0.024 (\Mmax/10^5 ~ \Msun)^{0.28}$, respectively. The massive model clouds follow the SFR-dense gas mass relation obtained by Gao \& Solomon for infrared galaxies, extrapolated down to MCs scales. Low-mass clouds fall above this relation, in agreement with recent observations. An integration of the model-predicted $\SFRavg$ over a Galactic GMC mass spectrum yields a realistic value for the Galactic SFR. Our results reinforce the suggestion that star-forming GMCs may be in global collapse, and still have low net SFRs and SFEs, due to the evaporation of most of the cloud material by the feedback from massive stars. △ Less

Submitted 23 July, 2014; v1 submitted 22 August, 2013; originally announced August 2013.

Comments: Accepted for publication in the ApJ

MSC Class: 85-01