Numerical simulations of type III planetary migration – II. Inward migration of massive planets

We present a numerical study of rapid, so-called type III migration for Jupiter-sized planets embedded in a protoplanetary disc. We limit ourselves to the case of inward migration, and study in detail its evolution and physics, concentrating on the structure of the corotation and circumplanetary regions, and processes for stopping migration. We also consider the dependence of the migration behaviour on several key parameters. We perform this study using the results of global, two-dimensional hydrodynamical simulations with adaptive mesh refinement. The initial conditions are chosen to satisfy the condition for rapid inward migration. We find that type III migration can be divided into two regimes, fast and slow. The structure of the co-orbital region, mass accumulation rate and migration behaviour differ between these two regimes. All our simulations show a transition from the fast to the slow regime, ending type III migration well before reaching the star. The stopping radius is found to be larger for more massive planets and less massive discs. A sharp density drop is also found to be an efficient stopping mechanism. In the fast migration regime the migration rate and induced eccentricity are lower for less massive discs, but almost do not depend on planet mass. Eccentricity is damped on the migration time-scale.     

Simulations of spiral galaxies with an active potential: molecular cloud formation and gas dynamics

We describe simulations of the response of a gaseous disc to an active spiral potential. The potential is derived from an N -body calculation and leads to a multi-armed time-evolving pattern. The gas forms long spiral arms typical of grand-design galaxies, although the spiral pattern is asymmetric. The primary difference from a grand-design spiral galaxy, which has a consistent two-/four-armed pattern, is that instead of passing through the spiral arms, gas generally falls into a developing potential minimum and is released only when the local minimum dissolves. In this case, the densest gas is coincident with the spiral potential, rather than offset as in the grand-design spirals. We would therefore expect no offset between the spiral shock and star formation, and no obvious corotation radius. Spurs which occur in grand-design spirals when large clumps are sheared off leaving the spiral arms, are rare in the active, time-evolving spiral reported here. Instead, large branches are formed from spiral arms when the underlying spiral potential is dissolving due to the N -body dynamics. We find that the molecular cloud mass spectrum for the active potential is similar to that for clouds in grand-design calculations, depending primarily on the ambient pressure rather than the nature of the potential. The largest molecular clouds occur when spiral arms collide, rather than by agglomeration within a spiral arm.     

Molecular line mapping of the giant molecular cloud associated with RCW 106 – II. Column density and dynamical state of the clumps

We present a fully sampled C¹⁸O (1–0) map towards the southern giant molecular cloud (GMC) associated with the H  ii region RCW 106, and use it in combination with previous ¹³CO (1–0) mapping to estimate the gas column density as a function of position and velocity. We find localized regions of significant ¹³CO optical depth in the northern part of the cloud, with several of the high-opacity clouds in this region likely associated with a limb-brightened shell around the H  ii region G333.6−0.2. Optical depth corrections broaden the distribution of column densities in the cloud, yielding a lognormal distribution as predicted by simulations of turbulence. Decomposing the ¹³CO and C¹⁸O data cubes into clumps, we find relatively weak correlations between size and linewidth, and a more sensitive dependence of luminosity on size than would be predicted by a constant average column density. The clump mass spectrum has a slope near −1.7, consistent with previous studies. The most massive clumps appear to have gravitational binding energies well in excess of virial equilibrium; we discuss possible explanations, which include magnetic support and neglect of time-varying surface terms in the virial theorem. Unlike molecular clouds as a whole, the clumps within the RCW 106 GMC, while elongated, appear to show random orientations with respect to the Galactic plane.     

Environment and mass dependencies of galactic λ spin parameter: cosmological simulations and observed galaxies compared

We use a sample of galaxies from the Sloan Digital Sky Survey (SDSS) to search for correlations between the λ spin parameter and the environment and mass of galaxies. In order to calculate the total value of λ for each observed galaxy, we employed a simple model of the dynamical structure of the galaxies, which allows a rough estimate of the value of λ using only readily obtainable observables from the luminous galaxies. Use of a large volume-limited sample (upwards of 11 000) allows reliable inferences of mean values and dispersions of λ distributions. We find, in agreement with some N -body cosmological simulations, no significant dependence of λ on the environmental density of the galaxies. For the case of mass, our results show a marked correlation with λ, in the sense that low-mass galaxies present both higher mean values of λ and associated dispersions, than high-mass galaxies. These results provide interesting constrain on the mechanisms of galaxy formation and acquisition of angular momentum, a valuable test for cosmological models.     

Molecular line profiles as diagnostics of protostellar collapse: modelling the 'blue asymmetry' in inside-out infall

The evolution of star-forming core analogues undergoing inside-out collapse is studied with a multipoint chemodynamical model which self-consistently computes the abundance distribution of chemical species in the core. For several collapse periods the output chemistry of infalling tracer species such as HCO⁺, CS and N₂H⁺ is then coupled to an accelerated Λ-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.