Probing the physics of interacting galaxies

The morphological and velocity structures in the gaseous (HI and CO) and stellar components of two interacting systems are examined. Both Arp 140 and Arp 104 reveal extended tidal tails in the HI. The Hα and FIR fluxes of Arp 140 yield similar SFR of ∼ 0.8 M_⊙ yr^-1. In contrast the Hα flux of Arp 104 yields a SFR of ∼ 0.05 M_⊙ yr^-1, ∼ 20 times smaller than that obtained from the FIR flux. Spectra were used to examine the changing velocity of atomic and molecular gas in NGC 5218 (Arp 104). The atomic and molecular gas were found to be dynamically similar with comparable velocities and velocity widths across the galaxy; consistent with the two phases responding similarly to the interaction, or enhanced HI to CO conversion in the centre of the galaxy. This revised version was published online in August 2006 with corrections to the Cover Date.     

The efficiency of the spiral-in of a black hole to the Galactic Centre

We study the efficiency at which a black hole or dense star cluster spirals in to the Galactic Centre. This process takes place on a dynamical friction time-scale, which depends on the value of the Coulomb logarithm (ln Λ). We determine the accurate value of this parameter using the direct N -body method, a tree algorithm and a particle-mesh technique with up to two million plus one particles. The three different techniques are in excellent agreement. Our measurement for the Coulomb logarithm appears to be independent of the number of particles. We conclude that  ln Λ= 6.6 ± 0.6  for a massive point particle in the inner few parsec of the Galactic bulge. For an extended object, such as a dense star cluster, ln Λ is smaller, with a value of the logarithm argument Λ inversely proportional to the object size.     

A numerical study of the pre-ejection, magnetically-sheared corona as a free boundary problem

A class of magnetostatic equilibria with axial symmetry outside a unit sphere in the presence of plasma pressure and an r ^–2 gravitational field is constructed. The structure contains a localized current-carrying region confined by a background bipolar potential field, and the shape of the region changes subject to the variation of the electric current. The continuity requirement for the magnetic field and plasma pressures at the outer boundary of the cavity defines a free boundary problem, which is solved numerically using a spectral boundary scheme. The model is then used to study the expansion of the current-carrying region, caused by the buildup of magnetic shear, against the background confining field. The magnetic shear in our model is induced by the loading of an azimuthal field, accompanied by a depletion of plasma density.We show that due to the additional effect of confinement by the dense surrounding plasma, the energy of the magnetic field can exceed the energy of its associated open field, presumably a necessary condition for the onset of coronal mass ejections. (However, the plasma beta of the confining fluid is higher than that in the outer boundary of a realistic helmet-streamer structure.) Furthermore, under the assumption that coronal mass ejections are driven by magnetic buoyancy, the result from our model study lends further support to the notion of a suspended magnetic flux rope in the low-density cavity of a helmet-streamer as a promising pre-ejection configuration.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Non-Newtonian topographic relaxation on Europa

Models of topographic support on Europa by lateral shell thickness variations have previously assumed a Newtonian ice viscosity. Here I show that using a more realistic stress-dependent viscosity gives relaxation times which can be significantly different. Topography of wavelength 100 km cannot be supported by lateral shell thickness variations for ∼50 Myr, unless the shell thickness is <10 km or the ice grain size >10 mm. Shorter wavelength topography would require even thinner shells, but may be supported elastically. Global-scale variations in shell thickness, however, can be supported for geological timescales if the shell thickness is O(10 km).