No Arabic abstract
We present a new systematic way of setting up galactic gas disks based on the assumption of detailed hydrodynamic equilibrium. To do this, we need to specify the density distribution and the velocity field which supports the disk. We first show that the required circular velocity has no dependence on the height above or below the midplane so long as the gas pressure is a function of density only. The assumption of disks being very thin enables us to decouple the vertical structure from the radial direction. Based on that, the equation of hydrostatic equilibrium together with the reduced Poisson equation leads to two sets of second-order non-linear differential equation, which are easily integrated to set-up a stable disk. We call one approach `density method and the other one `potential method. Gas disks in detailed balance are especially suitable for investigating the onset of the gravitational instability. We revisit the question of global, axisymmetric instability using fully three-dimensional disk simulations. The impact of disk thickness on the disk instability and the formation of spontaneously induced spirals is studied systematically with or without the presence of the stellar potential. In our models, the numerical results show that the threshold value for disk instability is shifted from unity to 0.69 for self-gravitating thick disks and to 0.75 for combined stellar and gas thick disks. The simulations also show that self-induced spirals occur in the correct regions and with the right numbers as predicted by the analytic theory.
Star formation has long been known to be an inefficient process, in the sense that only a small fraction $epsilon_{rm ff}$ of the mass of any given gas cloud is converted to stars per cloud free-fall time. However, developing a successful theory of star formation will require measurements of both the mean value of $epsilon_{rm ff}$ and its scatter from one molecular cloud to another. Because $epsilon_{rm ff}$ is measured relative to the free-fall time, such measurements require accurate determinations of cloud volume densities. Efforts to measure the volume density from two-dimensional projected data, however, have thus far relied on treating molecular clouds as simple uniform spheres, while their real shapes are likely filamentary and their density distributions far from uniform. The resulting uncertainty in the true volume density is likely one of the major sources of error in observational estimates of $epsilon_{rm ff}$. In this paper, we use a suite of simulations of turbulent, magnetized, radiative, self-gravitating star-forming clouds to examine whether it is possible to obtain more accurate volume density estimates and thereby reduce this error. We create mock observations from simulations, and show that current analysis methods relying on the spherical assumption likely yield ~ 0.26 dex underestimations and ~ 0.51 dex errors in volume density estimates, corresponding to a ~ 0.13 dex overestimation and a ~ 0.25 dex scatter in $epsilon_{rm ff}$, comparable to the scatter in observed cloud samples. We build a predictive model that uses information accessible in two-dimensional measurements -- most significantly the Gini coefficient of the surface density distribution -- to estimate volume density with ~ 0.3 dex less scatter. We test our method on a recent observation of the Ophiuchus cloud, and show that it successfully reduces the $epsilon_{rm ff}$ scatter.
Carbon monoxide (CO) is the best tracer of Galactic molecular hydrogen (H2). Its lowest rotational emission lines are in the radio regime and thanks to Galactic rotation emission at different distances is Doppler shifted. For a given gas flow model the observed spectra can thus be deprojected along the line of sight to infer the gas distribution. We use the CO line survey of Dame et al. (2001) to reconstruct the three-dimensional density of H2. We consider the deprojection as a Bayesian variational inference problem. The posterior distribution of the gas densities allows us to estimate both the mean and uncertainty of the reconstructed density. Unlike most of the previous attempts, we take into account the correlations of gas on a variety of scales which allows curing some of the well-known pathologies, like fingers-of-god effects. Both gas flow models that we adopt incorporate a Galactic bar which induces radial motions in the inner few kiloparsecs and thus offers spectral resolution towards the Galactic centre. We compare our gas maps with those of earlier studies and characterise their statistical properties, e.g. the radial profile of the average surface mass density. We have made our three-dimensional gas maps and their uncertainties available to the community at https://dx.doi.org/10.5281/zenodo.4405437 .
Aims. We investigate the spatial distribution of a collection of absorbing gas clouds, some associated with the dense, massive star-forming core NGC6334 I, and others with diffuse foreground clouds. For the former category, we aim to study the dynamical properties of the clouds in order to assess their potential to feed the accreting protostellar cores. Methods. We use spectral imaging from the Herschel SPIRE iFTS to construct a map of HF absorption at 243 micron in a 6x3.5 arcmin region surrounding NGC6334 I and I(N). Results. The combination of new, spatially fully sampled, but spectrally unresolved mapping with a previous, single-pointing, spectrally resolved HF signature yields a 3D picture of absorbing gas clouds in the direction of NGC6334. Toward core I, the HF equivalent width matches that of the spectrally resolved observation. The distribution of HF absorption is consistent with three of the seven components being associated with this dense star-forming envelope. For two of the remaining four components, our data suggest that these clouds are spatially associated with the larger scale filamentary star-forming complex. Our data also implies a lack of gas phase HF in the envelope of core I(N). Using a simple description of adsorption onto and desorption from dust grain surfaces, we show that the overall lower temperature of the envelope of source I(N) is consistent with freeze-out of HF, while it remains in the gas phase in source I. Conclusions. We use the HF molecule as a tracer of column density in diffuse gas (n(H) ~ 10^2 - 10^3 cm^-3), and find that it may uniquely trace a relatively low density portion of the gas reservoir available for star formation that otherwise escapes detection. At higher densities prevailing in protostellar envelopes (>10^4 cm^-3), we find evidence of HF depletion from the gas phase under sufficiently cold conditions.
This paper considers gravitational perturbations in geometrically thin disks with rotation curves dominated by a central object, but with substantial contributions from magnetic pressure and tension. The treatment is general, but the application is to the circumstellar disks that arise during the gravitational collapse phase of star formation. We find the dispersion relation for spiral density waves in these generalized disks and derive the stability criterion for axisymmetric $(m=0)$ disturbances (the analog of the Toomre parameter $Q_T$) for any radial distribution of the mass-to-flux ratio $lambda$. The magnetic effects work in two opposing directions: on one hand, magnetic tension and pressure stabilize the disk against gravitational collapse and fragmentation; on the other hand, they also lower the rotation rate making the disk more unstable. For disks around young stars the first effect generally dominates, so that magnetic fields allow disks to be stable for higher surface densities and larger total masses. These results indicate that magnetic fields act to suppress the formation of giant planets through gravitational instability. Finally, even if gravitational instability can form a secondary body, it must lose an enormous amount of magnetic flux in order to become a planet; this latter requirement represents an additional constraint for planet formation via gravitational instability and places a lower limit on the electrical resistivity.
We examine the gas circulation near a gap opened by a giant planet in a protoplanetary disk. We show with high resolution 3D simulations that the gas flows into the gap at high altitude over the mid-plane, at a rate dependent on viscosity. We explain this observation with a simple conceptual model. From this model we derive an estimate of the amount of gas flowing into a gap opened by a planet with Hill radius comparable to the scale-height of a layered disk (i. e. a disk with viscous upper layer and inviscid midplane). Our estimate agrees with modern MRI simulations(Gressel et al., 2013). We conclude that gap opening in a layered disk can not slow down significantly the runaway gas accretion of Saturn to Jupiter-mass planets.