-
Unveiling two expanding stellar groups formed through violent relaxation in The Lagoon Nebula Cluster
Authors:
A. Bonilla-Barroso,
J. Ballesteros-Paredes,
J. Hernández,
L. Aguilar,
M. Zamora-Avilés
Abstract:
The current kinematic state of young stellar clusters can give clues on their actual dynamical state and origin. In this contribution, we use Gaia DR3 data of the Lagoon Nebula Cluster (LNC) to show that the cluster is composed of two expanding groups, likely formed from different molecular cloud clumps. We find no evidence of massive stars having larger velocity dispersion than low-mass stars or…
▽ More
The current kinematic state of young stellar clusters can give clues on their actual dynamical state and origin. In this contribution, we use Gaia DR3 data of the Lagoon Nebula Cluster (LNC) to show that the cluster is composed of two expanding groups, likely formed from different molecular cloud clumps. We find no evidence of massive stars having larger velocity dispersion than low-mass stars or being spatially segregated across the LNC, as a whole, or within the Primary group. However, the Secondary group, with 1/5th of the stars, exhibits intriguing features. On the one hand, it shows a bipolar nature, with an aspect ratio of $\sim$3:1. In addition, the massive stars in this group exhibit larger velocity dispersion than the low-mass stars, although they are not concentrated towards the center of the group. This suggests that this group may have undergone dynamical relaxation, first, and some explosive event afterward. However, further observations and numerical work have to be performed to confirm this hypothesis. The results of this work suggest that, although stellar clusters may form by the global and hierarchical collapse of their parent clump, still some dynamical relaxation may take place.
△ Less
Submitted 29 February, 2024;
originally announced March 2024.
-
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.
-
The stellar population responsible for a kiloparsec size superbubble seen in the JWST "phantom" images of NGC628
Authors:
Y. D. Mayya,
J. A. Alzate,
L. Lomelí-Núñez,
J. Zaragoza-Cardiel,
V. M. A. Gómez-González,
S. Silich,
D. Fernández-Arenas,
O. Vega,
P. A. Ovando,
L. H. Rodríguez,
D. Rosa-González,
A. Luna,
M. Zamora-Avilés,
F. Rosales-Ortega
Abstract:
We here study the multi-band properties of a kiloparsec-size superbubble in the late-type spiral galaxy NGC628. The superbubble is the largest of many holes seen in the early release images using JWST/MIRI filters that trace the Polycyclic Aromatic Hydrocarbon (PAH) emissions. The superbubble is located in the interarm region ~3 kpc from the galactic center in the south-east direction. The shell s…
▽ More
We here study the multi-band properties of a kiloparsec-size superbubble in the late-type spiral galaxy NGC628. The superbubble is the largest of many holes seen in the early release images using JWST/MIRI filters that trace the Polycyclic Aromatic Hydrocarbon (PAH) emissions. The superbubble is located in the interarm region ~3 kpc from the galactic center in the south-east direction. The shell surrounding the superbubble is detected in HI, CO, and Halpha with an expansion velocity of 12 km/s, and contains as much as 2x10^7 Msun of mass in gas that is mostly in molecular form. We find a clear excess of blue, bright stars inside the bubble as compared to the surrounding disk on the HST/ACS images. These excess blue, bright stars are part of a stellar population of 10^5 Msun mass that is formed over the last 50 Myr in different star formation episodes, as determined from an analysis of color-magnitude diagrams using a Bayesian technique. The mechanical power injected by the massive stars of these populations is sufficient to provide the energy necessary for the expansion of the shell gas. Slow and steady, rather than violent, injection of energy is probably the reason for the maintenance of the shell structure over the kiloparsec scale. The expanding shell is currently the site for triggered star formation as inferred from the JWST 21 micron (F2100W filter) and the Halpha images.
△ Less
Submitted 24 February, 2023;
originally announced February 2023.
-
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.
-
Why most molecular clouds are gravitationally dominated
Authors:
Laura Ramírez-Galeano,
Javier Ballesteros-Paredes,
Rowan Smith,
Vianey Camacho,
Manuel Zamora-Aviles
Abstract:
Observational and theoretical evidence suggests that a substantial population of molecular clouds (MCs) appear to be unbound, dominated by turbulent motions. However, these estimations are made typically via the so-called viral parameter $α_{\rm vir}^{\rm class}$, which is an observational proxy to the virial ratio between the kinetic and the gravitational energy. This parameter intrinsically assu…
▽ More
Observational and theoretical evidence suggests that a substantial population of molecular clouds (MCs) appear to be unbound, dominated by turbulent motions. However, these estimations are made typically via the so-called viral parameter $α_{\rm vir}^{\rm class}$, which is an observational proxy to the virial ratio between the kinetic and the gravitational energy. This parameter intrinsically assumes that MCs are isolated, spherical, and with constant density. However, MCs are embedded in their parent galaxy and thus are subject to compressive and disruptive tidal forces from their galaxy, exhibit irregular shapes, and show substantial substructure. We, therefore, compare the typical estimations of $α_{\rm vir}^{\rm class}$ to a more precise definition of the virial parameter, $α_{\rm vir}^{\rm full}$, which accounts not only for the self-gravity (as $α_{\rm vir}^{\rm class}$), but also for the tidal stresses, and thus, it can take negative (self-gravity) and positive (tides) values. While we recover the classical result that most of the clouds appear to be unbound, having $α_{\rm vir}^{\rm class} > 2$, we show that, with the more detailed definition considering the full gravitational energy, (i) 50\%\ of the total population is gravitationally bound, however, (ii) another 20\%\ is gravitationally dominated, but with tides tearing them apart; (iii) the source of those tides does not come from the galactic structure (bulge, halo, spiral arms), but from the molecular cloud complexes in which clouds reside, and probably (iv) from massive young stellar complexes, if they were present. (v) Finally, our results also suggest that, interstellar turbulence can have, at least partially, a gravitational origin.
△ Less
Submitted 18 June, 2022;
originally announced June 2022.
-
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.
-
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.
-
What is the physics behind the Larson mass-size relation?
Authors:
Javier Ballesteros-Paredes,
Carlos Román-Zúñiga,
Quentin Salomé,
Manuel Zamora-Avilés,
María Jesús Jiménez-Donaire
Abstract:
Different studies have reported a power-law mass-size relation $M \propto R^q$ for ensembles of molecular clouds. In the case of nearby clouds, the index of the power-law $q$ is close to 2. However, for clouds spread all over the Galaxy, indexes larger than 2 are reported. We show that indexes larger than 2 could be the result of line-of-sight superposition of emission that does not belong to the…
▽ More
Different studies have reported a power-law mass-size relation $M \propto R^q$ for ensembles of molecular clouds. In the case of nearby clouds, the index of the power-law $q$ is close to 2. However, for clouds spread all over the Galaxy, indexes larger than 2 are reported. We show that indexes larger than 2 could be the result of line-of-sight superposition of emission that does not belong to the cloud itself. We found that a random factor of gas contamination, between 0.001\%\ and 10\%\ of the line-of-sight, allows to reproduce the mass-size relation with $q \sim 2.2-2.3$ observed in Galactic CO surveys. Furthermore, for dense cores within a single cloud, or molecular clouds within a single galaxy, we argue that, even in these cases, there is observational and theoretical evidence that some degree of superposition may be occurring. However, additional effects may be present in each case, and are briefly discussed. We also argue that defining the fractal dimension of clouds via the mass-size relation is not adequate, since the mass is not {necessarily} a proxy to the area, and the size reported in $M-R$ relations is typically obtained from the square root of the area, rather than from an estimation of the size independent from the area. Finally, we argue that the statistical analysis of finding clouds satisfying the Larson's relations does not mean that each individual cloud is in virial equilibrium.
△ Less
Submitted 6 September, 2019;
originally announced September 2019.
-
Flipping-up the field: gravitational feedback as a mechanism for young clusters dispersal
Authors:
Manuel Zamora-Avilés,
Javier Ballesteros-Paredes,
Jesús Hernández,
Carlos Román-Zúñiga,
Verónica Lora,
Marina Kounkel
Abstract:
Recent analyses of Gaia data have provided direct evidence that most young stellar clusters are in a state of expansion, with velocities of the order of ~0.5 km/s. Traditionally, expanding young clusters have been pictured as entities that became unbound due to the lack of gravitational binding once the gas from the parental cloud that formed the cluster has been expelled by the stellar radiation…
▽ More
Recent analyses of Gaia data have provided direct evidence that most young stellar clusters are in a state of expansion, with velocities of the order of ~0.5 km/s. Traditionally, expanding young clusters have been pictured as entities that became unbound due to the lack of gravitational binding once the gas from the parental cloud that formed the cluster has been expelled by the stellar radiation of the massive stars in the cluster. In the present contribution, we used radiation-magnetohydrodynamic numerical simulations of molecular cloud formation and evolution to understand how stellar clusters form and disperse. We found that the ionising feedback from the newborn massive stars expels the gas from the collapse centre, flipping-up the gravitational potential as a consequence of the mass removal from the inside-out. Since neither the parental clouds nor the formed shells are distributed symmetrically around the HII region, net forces pulling out the stars are present, accelerating them towards the edges of the cavity. We call this mechanism ``gravitational feedback", in which the gravity from the expelled gas appears to be the crucial mechanism producing unbound clusters that expand away from their formation centre in an accelerated way in young stellar clusters. This mechanism naturally explains the "Hubble flow-like" expansion observed in several young clusters.
△ Less
Submitted 9 July, 2019;
originally announced July 2019.
-
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.
-
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.
-
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.
-
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.
-
Are fibres in molecular cloud filaments real objects?
Authors:
Manuel Zamora-Avilés,
Javier Ballesteros-Paredes,
Lee W. Hartmann
Abstract:
We analyse the morphology and kinematics of dense filamentary structures produced in a numerical simulation of a star--forming cloud of $1.4 \times 10^4 \, \Msun$ evolving under their own self--gravity in a magnetized medium. This study is motivated by recent observations of velocity--coherent substructures ("fibres") in star-forming filaments. We find such "fibres" ubiquitously in our simulated f…
▽ More
We analyse the morphology and kinematics of dense filamentary structures produced in a numerical simulation of a star--forming cloud of $1.4 \times 10^4 \, \Msun$ evolving under their own self--gravity in a magnetized medium. This study is motivated by recent observations of velocity--coherent substructures ("fibres") in star-forming filaments. We find such "fibres" ubiquitously in our simulated filament. We found that a fibre in one projection is not necessarily a fibre in another projection, and thus, caution should be taken into account when considering them as real objects. We found that only the densest parts of the filament ($\sim$30\% of the densest volume, which contains $\sim$70\% of the mass) belong to fibres in 2 projections. Moreover, it is quite common that they are formed by separated density enhancements superposed along the line of sight. Observations of fibres can yield insight into the level of turbulent substructure driven by gravity, but care should be taken in interpreting the results given the problem of line of sight superposition. We also studied the morphology and kinematics of the 3D skeleton (spine), finding that subfilaments accrete structured material mainly along the magnetic field lines, which are preferentially perpendicular to the skeleton. The magnetic field is at the same time dragged by the velocity field due to the gravitational collapse.
△ Less
Submitted 12 August, 2017; v1 submitted 4 August, 2017;
originally announced August 2017.
-
Magnetized converging flows towards the hot core in the intermediate/high-mass star-forming region NGC 6334 V
Authors:
Carmen Juárez,
Josep M. Girart,
Manuel Zamora-Avilés,
Ya-Wen Tang,
Patrick M. Koch,
Hauyu Baobab Liu,
Aina Palau,
Javier Ballesteros-Paredes,
Qizhou Zhang,
Keping Qiu
Abstract:
We present Submillimeter Array (SMA) observations at 345 GHz towards the intermediate/high-mass cluster-forming region NGC 6334 V. From the dust emission we spatially resolve three dense condensations, the brightest one presenting the typical chemistry of a hot core. The magnetic field (derived from the dust polarized emission) shows a bimodal converging pattern towards the hot core. The molecular…
▽ More
We present Submillimeter Array (SMA) observations at 345 GHz towards the intermediate/high-mass cluster-forming region NGC 6334 V. From the dust emission we spatially resolve three dense condensations, the brightest one presenting the typical chemistry of a hot core. The magnetic field (derived from the dust polarized emission) shows a bimodal converging pattern towards the hot core. The molecular emission traces two filamentary structures at two different velocities, separated by 2 km/s, converging to the hot core and following the magnetic field distribution. We compare the velocity field and the magnetic field derived from the SMA observations with MHD simulations of star-forming regions dominated by gravity. This comparison allows us to show how the gas falls in from the larger-scale extended dense core (~0.1 pc) of NGC 6334 V towards the higher-density hot core region (~0.02 pc) through two distinctive converging flows dragging the magnetic field, whose strength seems to have been overcome by gravity.
△ Less
Submitted 12 June, 2017;
originally announced June 2017.
-
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.
-
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.
-
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.
-
An Evolutionary Model for Collapsing Molecular Clouds and Their Star Formation Activity
Authors:
Manuel Zamora-Aviles,
Enrique Vázquez-Semadeni,
Pedro Colín
Abstract:
We present an idealized, semi-empirical model for the evolution of gravitationally contracting molecular clouds (MCs) and their star formation rate (SFR) and efficiency (SFE). The model assumes that the instantaneous SFR is given by the mass above a certain density threshold divided by its free-fall time. The instantaneous number of massive stars is computed assuming a Kroupa IMF. These stars feed…
▽ More
We present an idealized, semi-empirical model for the evolution of gravitationally contracting molecular clouds (MCs) and their star formation rate (SFR) and efficiency (SFE). The model assumes that the instantaneous SFR is given by the mass above a certain density threshold divided by its free-fall time. The instantaneous number of massive stars is computed assuming a Kroupa IMF. These stars feed back on the cloud through ionizing radiation, eroding it. The main controlling parameter of the evolution turns out to be the maximum cloud mass, $\Mmax$. This allows us to compare various properties of the model clouds against their observational counterparts. A giant molecular cloud (GMC) model ($\Mmax \sim 10^5 \Msun$) adheres very well to the evolutionary scenario recently inferred by Kawamura et al. (2009) for GMCs in the Large Magellanic Cloud. A model cloud with $\Mmax \approx 2000 \Msun$ evolves in the Kennicutt-Schmidt diagram first passing through the locus of typical low- to-intermediate mass star-forming clouds, and then moving towards the locus of high-mass star-forming ones over the course of $\sim 10$ Myr. Also, the stellar age histograms for this cloud a few Myr before its destruction agree very well with those observed in the $ρ$-Oph stellar association, whose parent cloud has a similar mass, and imply that the SFR of the clouds increases with time. Our model thus agrees well with various observed properties of star-forming MCs, suggesting that the scenario of gravitationally collapsing MCs, with their SFR regulated by stellar feedback, is entirely feasible and in agreement with key observed properties of molecular clouds.
△ Less
Submitted 29 March, 2012; v1 submitted 24 May, 2011;
originally announced May 2011.
-
Gravity or Turbulence? The velocity dispersion-size relation
Authors:
Javier Ballesteros-Paredes,
Lee W. Hartmann,
Enrique Vázquez-Semadeni,
Fabian Heitsch,
Manuel A. Zamora-Avilés
Abstract:
We discuss the nature of the velocity dispersion vs. size relation for molecular clouds. In particular, we add to previous observational results showing that the velocity dispersions in molecular clouds and cores are not purely functions of spatial scale but involve surface gas densities as well. We emphasize that hydrodynamic turbulence is required to produce the first condensations in the progen…
▽ More
We discuss the nature of the velocity dispersion vs. size relation for molecular clouds. In particular, we add to previous observational results showing that the velocity dispersions in molecular clouds and cores are not purely functions of spatial scale but involve surface gas densities as well. We emphasize that hydrodynamic turbulence is required to produce the first condensations in the progenitor medium. However, as the cloud is forming, it also becomes bound, and gravitational accelerations dominate the motions. Energy conservation in this case implies $|E_g| \sim E_k$, in agreement with observational data, and providing an interpretation for two recent observational results: the scatter in the $δv-R$ plane, and the dependence of the velocity dispersion on the surface density ${δv^2/ R} \propto Σ$. We argue that the observational data are consistent with molecular clouds in a state of hierarchical gravitational collapse, i.e., developing local centers of collapse throughout the whole cloud while the cloud itself is collapsing, and making equilibrium unnecessary at all stages prior to the formation of actual stars. Finally, we discuss how this mechanism need not be in conflict with the observed star formation rate.
△ Less
Submitted 8 September, 2010;
originally announced September 2010.