No Arabic abstract
In cold molecular clouds submillimetre emission lines are excited by the ambient radiation field. The pumping is dominated by the cosmic microwave background (CMB). It is usual in molecular line radiative transfer modelling to simply assume that this is the only incident radiation field. In this paper, a molecular line transport code and a dust radiative transfer code are used to explore the effects of the inclusion of a full interstellar radiation field (ISRF) on a simple test molecular cloud. It is found that in many galactic situations, the shape and strength of the line profiles that result are robust to variations in the ISRF and thus that in most cases, it is safe to adopt the CMB radiation field for the molecular line transport calculations. However, we show that in two examples, the inclusion of a plausible radiation field can have a significant effect on the the line profiles. Firstly, in the vicinity of an embedded massive star, there will be an enhanced far infared component to the radiation field. Secondly, for molecular clouds at large redshift, the CMB temperature increases and this of course also alters the radiation field. In both of these cases, the line profiles are weakened significantly compared to a cloud exposed to a standard radiation field. Therefore this effect should be accounted for when investigating prestellar cores in massive star forming regions and when searching for molecular clouds at high redshift.
The evolution of star-forming core analogues undergoing inside-out collapse is studied with a multi-point chemodynamical model which self-consistently computes the abundance distribution of chemical species in the core. For several collapse periods the output chemistry of infall tracer species such as HCO+, CS, and N2H+, is then coupled to an accelerated Lambda-iteration radiative transfer code, which predicts the emerging molecular line profiles using two different input gas/dust temperature distributions. We investigate the sensitivity of the predicted spectral line profiles and line asymmetry ratios to the core temperature distribution, the time-dependent model chemistry, as well as to ad hoc abundance distributions. The line asymmetry is found to be strongly dependent on the adopted chemical abundance distribution. In general, models with a warm central region show higher values of blue asymmetry in optically thick HCO+ and CS lines than models with a starless core temperature profile. We find that in the formal context of Shu-type inside-out infall, and in the absence of rotation or outflows, the relative blue asymmetry of certain HCO+ and CS transitions is a function of time and, subject to the foregoing caveats, can act as a collapse chronometer. The sensitivity of simulated HCO+ line profiles to linear radial variations, subsonic or supersonic, of the internal turbulence field is investigated in the separate case of static cores.
Optical stellar polarimetry in the Perseus molecular cloud direction is known to show a fully mixed bi-modal distribution of position angles across the cloud (Goodman et al. 1990). We study the Gaia trigonometric distances to each of these stars and reveal that the two components in position angles trace two different dust clouds along the line of sight. One component, which shows a polarization angle of -37.6 deg +/- 35.2 deg and a higher polarization fraction of 2.0 +/- 1.7%, primarily traces the Perseus molecular cloud at a distance of 300 pc. The other component, which shows a polarization angle of +66.8 deg +/- 19.1 deg and a lower polarization fraction of 0.8 +/- 0.6%, traces a foreground cloud at a distance of 150 pc. The foreground cloud is faint, with a maximum visual extinction of < 1 mag. We identify that foreground cloud as the outer edge of the Taurus molecular cloud. Between the Perseus and Taurus molecular clouds, we identify a lower-density ellipsoidal dust cavity with a size of 100 -- 160 pc. This dust cavity locates at l = 170 deg, b = -20 deg, and d = 240 pc, which corresponds to an HI shell generally associated with the Per OB2 association. The two-component polarization signature observed toward the Perseus molecular cloud can therefore be explained by a combination of the plane-of-sky orientations of the magnetic field both at the front and at the back of this dust cavity.
Two aspects of filamentary molecular cloud evolution are addressed: (1) Exploring analytically the role of the environment for the evolution of filaments demonstrates that considering them in isolation (i.e. just addressing the fragmentation stability) will result in unphysical conclusions about the filaments properties. Accretion can also explain the observed decorrelation between FWHM and peak column density. (2) Free-fall accretion onto finite filaments can lead to the characteristic fans of infrared-dark clouds around star-forming regions. The fans may form due to tidal forces mostly arising at the ends of the filaments, consistent with numerical models and earlier analytical studies.
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 use N-body simulations to investigate the radial dependence of the density and velocity dispersion in cold dark matter (CDM) halos. In particular, we explore how closely Q rho/sigma^3, a surrogate measure of the phase-space density, follows a power-law in radius. Our study extends earlier work by considering, in addition to spherically-averaged profiles, local Q-estimates for individual particles, Q_i; profiles based on the ellipsoidal radius dictated by the triaxial structure of the halo, Q_i(r); and by carefully removing substructures in order to focus on the profile of the smooth halo, Q^s. The resulting Q_i^s(r) profiles follow closely a power law near the center, but show a clear upturn from this trend near the virial radius, r_{200}. The location and magnitude of the deviations are in excellent agreement with the predictions from Bertschingers spherical secondary-infall similarity solution. In this model, Q propto r^{-1.875} in the inner, virialized regions, but departures from a power-law occur near r_{200} because of the proximity of this radius to the location of the first shell crossing - the shock radius in the case of a collisional fluid. Particles there have not yet fully virialized, and so Q departs from the inner power-law profile. Our results imply that the power-law nature of $Q$ profiles only applies to the inner regions and cannot be used to predict accurately the structure of CDM halos beyond their characteristic scale radius.