We study the evolution of the virial parameter, $\alpha$, and the star formation activity of a star forming region in a numerical simulation and two observational samples likely to be at different evolutionary stages, the Pipe and the G14.225 clouds. We consider a numerical simulation of turbulence in the warm atomic gas, in which clouds form by a combination of thermal instability and turbulent compressions. In this simulation we select a region with physical properties similar to those of the Pipe in both its global physical properties and the dense cores it contains, and show that this region evolves to become similar to the G14.225 cloud and its substructure in the course of $\sim$1.5 Myr. Within this region, we follow the evolution of $\alpha$ and the star formation activity of the cores over this time interval, during which they grow in mass by accretion and develop further substructure of higher density. We find that the individual cores evolve by exhibiting first a mild decrease in $\alpha$ followed by a rapid increase when star formation begins. We suggest this evolution is due to an early loss of external compressive kinetic energy followed by an increase of the gravitationally driven motions. Nevertheless, collectively, the ensamble of clumps and cores reproduces the recently observed trend of decreasing $\alpha$ for objects of higher mass. We also find that the star formation rate and star formation efficiency increase monotonically as the region evolves. Thus, we propose that the energy balance, $\alpha$, and the star formation activity determine the evolutionary stage of a star forming region, and that a star forming region like the Pipe may evolve into a G14.225-like stage over the course of a few megayears if it is allowed to accrete from its environment. This evolution is consistent with the recently proposed scenario of hierarchical gravitational collapse.

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.

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.

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.

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.

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^\alpha$, with $\alpha$ 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.

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.

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.

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.

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.

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 $\delta v-R$ plane, and the dependence of the velocity dispersion on the surface density ${\delta v^2/ R} \propto \Sigma$. 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.