No Arabic 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 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 vequiv(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 address the turbulent fragmentation scenario for the origin of the stellar initial mass function (IMF), using a large set of numerical simulations of randomly driven supersonic MHD turbulence. The turbulent fragmentation model successfully predicts the main features of the observed stellar IMF assuming an isothermal equation of state without any stellar feedback. As a test of the model, we focus on the case of a magnetized isothermal gas, neglecting stellar feedback, while pursuing a large dynamic range in both space and timescales covering the full spectrum of stellar masses from brown dwarfs to massive stars. Our simulations represent a generic 4 pc region within a typical Galactic molecular cloud, with a mass of 3000 Msun and an rms velocity 10 times the isothermal sound speed and 5 times the average Alfven velocity, in agreement with observations. We achieve a maximum resolution of 50 au and a maximum duration of star formation of 4.0 Myr, forming up to a thousand sink particles whose mass distribution closely matches the observed stellar IMF. A large set of medium-size simulations is used to test the sink particle algorithm, while larger simulations are used to test the numerical convergence of the IMF and the dependence of the IMF turnover on physical parameters predicted by the turbulent fragmentation model. We find a clear trend toward numerical convergence and strong support for the model predictions, including the initial time evolution of the IMF. We conclude that the physics of isothermal MHD turbulence is sufficient to explain the origin of the IMF.
Spatially resolved spectroscopy from SDSS-IV MaNGA survey has revealed a class of quiescent, relatively common early-type galaxies, termed red geysers, that possibly host large scale active galactic nuclei driven winds. Given their potential importance in maintaining low level of star formation at late times, additional evidence confirming that winds are responsible for the red geyser phenomenon is critical. In this work, we present follow-up observations with the Echellette Spectrograph and Imager (ESI) at the Keck telescope of two red geysers (z$<$0.1) using multiple long slit positions to sample different regions of each galaxy. Our ESI data with a spectral resolution (R) $sim$ 8000 improves upon MaNGAs resolution by a factor of four, allowing us to resolve the ionized gas velocity profiles along the putative wind cone with an instrumental resolution of $rm sigma = 16~km~s^{-1}$. The line profiles of H$alpha$ and [NII]$rm lambda 6584$ show asymmetric shapes that depend systematically on location $-$ extended blue wings on the red-shifted side of the galaxy and red wings on the opposite side. We construct a simple wind model and show that our results are consistent with geometric projections through an outflowing conical wind oriented at an angle towards the line of sight. An alternative hypothesis that assigns the asymmetric pattern to beam-smearing of a rotating, ionized gas disk does a poor job matching the line asymmetry profiles. While our study features just two sources, it lends further support to the notion that red geysers are the result of galaxy-scale winds.
We carried out synthetic observations of interstellar atomic hydrogen at 21cm wavelength by utilizing the magneto-hydrodynamical numerical simulations of the inhomogeneous turbulent interstellar medium (ISM) Inoue and Inutsuka (2012). The cold neutral medium (CNM) shows significantly clumpy distribution with a small volume filling factor of 3.5%, whereas the warm neutral medium (WNM) distinctly different smooth distribution with a large filling factor of 96.5%. In projection on the sky, the CNM exhibits highly filamentary distribution with a sub-pc width, whereas the WNM shows smooth extended distribution. In the HI optical depth the CNM is dominant and the contribution of the WNM is negligibly small. The CNM has an area covering factor of 30% in projection, while the WNM has a covering factor of 70%. This causes that the emission-absorption measurements toward radio continuum compact sources tend to sample the WNM with a probability of 70%, yielding smaller HI optical depth and smaller HI column density than those of the bulk HI gas. The emission-absorption measurements, which are significantly affected by the small-scale large fluctuations of the CNM properties, are not suitable to characterize the bulk HI gas. Larger-beam emission measurements which are able to fully sample the HI gas will provide a better tool for that purpose, if a reliable proxy for hydrogen column density, possibly dust optical depth and gamma rays, is available.
Submillimetre-luminous galaxies at high-redshift are the most luminous, heavily star-forming galaxies in the Universe, and are characterised by prodigious emission in the far-infrared at 850 microns (S850 > 5 mJy). They reside in halos ~ 10^13Msun, have low gas fractions compared to main sequence disks at a comparable redshift, trace complex environments, and are not easily observable at optical wavelengths. Their physical origin remains unclear. Simulations have been able to form galaxies with the requisite luminosities, but have otherwise been unable to simultaneously match the stellar masses, star formation rates, gas fractions and environments. Here we report a cosmological hydrodynamic galaxy formation simulation that is able to form a submillimetre galaxy which simultaneously satisfies the broad range of observed physical constraints. We find that groups of galaxies residing in massive dark matter halos have rising star formation histories that peak at collective rates ~ 500-1000 Msun/yr at z=2-3, by which time the interstellar medium is sufficiently enriched with metals that the region may be observed as a submillimetre-selected system. The intense star formation rates are fueled in part by a reservoir gas supply enabled by stellar feedback at earlier times, not through major mergers. With a duty cycle of nearly a gigayear, our simulations show that the submillimetre-luminous phase of high-z galaxies is a drawn out one that is associated with significant mass buildup in early Universe proto-clusters, and that many submillimetre-luminous galaxies are actually composed of numerous unresolved components (for which there is some observational evidence).
We consider the dynamics in and near galaxy clusters. Gas, dark matter and galaxies are presently falling into the clusters between approximately 1 and 5 virial radii. At very large distances, beyond 10 virial radii, all matter is following the Hubble flow, and inside the virial radius the matter particles have on average zero radial velocity. The cosmological parameters are imprinted on the infall profile of the gas, however, no method exists, which allows a measurement of it. We consider the results of two cosmological simulations (using the numerical codes RAMSES and Gadget) and find that the gas and dark matter radial velocities are very similar. We derive the relevant dynamical equations, in particular the generalized hydrostatic equilibrium equation, including both the expansion of the Universe and the cosmological background. This generalized gas equation is the main new contribution of this paper. We combine these generalized equations with the results of the numerical simulations to estimate the contribution to the measured cluster masses from the radial velocity: inside the virial radius it is negligible, and inside two virial radii the effect is below 40%, in agreement the earlier analyses for DM. We point out how the infall velocity in principle may be observable, by measuring the gas properties to distance of about two virial radii, however, this is practically not possible today.