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, $\delta v\equiv(V_{thick}-V_{thin})/\Delta 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.

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.

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.

We report on high-sensitivity and high-angular resolution archival Submillimeter Array (SMA) observations of the large ($\sim$15000 AU) putative circumstellar disk associated with the O-type protostar NGC 7538 S. Observations of the continuum resolve this putative circumstellar disk into five compact sources, with sizes $\sim$ 3000 AU and masses $\sim10 M_\odot$. This confirm the results of recent millimeter observations made with CARMA/BIMA towards this object. However, we find that from most of these compact sources eject collimated bipolar outflows, revealed by our silicon monoxide (SiO {\it J}=5-4) observations and confirm that these sources have a (proto)stellar nature. All outflows are perpendicular to the large and rotating dusty structure. We propose therefore that, rather than being a single massive circumstellar disk, NGC 7538 S could be instead a large and massive contracting or rotating filament that is fragmenting at scales of 0.1 to 0.01 pc to form several B-type stars, via the standard process involving outflows and disks. As in recent high spatial resolution studies of dusty filaments, our observations also suggest that thermal pressure does not seem to be sufficient to support the filament, so that either additional support needs to be invoked, or else the filament must be in the process of collapsing. An SPH numerical simulation of the formation of a molecular cloud by converging warm neutral medium flows produces contracting filaments whose dimensions and spacings between the stars forming within them, as well as their column densities, strongly resemble those observed in the filament reported here.