The SFR and IMF of the galactic disk

There is a long term dynamical heating of stellar populations with age observed in the age – velocity dispersion – relation (AVR). This effect allows a determination of the star formation history SFR(t) from local kinematical data of main sequence stars. Using a self-consistent disk model for the vertical structure of the disk, we find from the kinematics of the stars in the solar neighbourhood that the SFR shows a moderate star burst about 10 Gyr ago followed by a continuous decline to the present day value consistent with the observed number of OB stars. The gravitational potential of the gas component and of the Dark Matter Halo is included and the effect of chemical enrichment, finite lifetime of the stars and mass loss of the stellar component are taken into account. The scale heights for main sequence stars together with the SFR is then used to determine constistently the IMF from the observed local luminosity function. The main new result is that the power law break in the present day mass function (PDMF) around 1 M _⊙ is entirely due to evolutionary effects of the disk and does not appear in the IMF. This revised version was published online in August 2006 with corrections to the Cover Date.     

The tidal stripping of satellites

We present an improved analytic calculation for the tidal radius of satellites and test our results against N -body simulations.
The tidal radius in general depends upon four factors: the potential of the host galaxy, the potential of the satellite, the orbit of the satellite and the orbit of the star within the satellite . We demonstrate that this last point is critical and suggest using three tidal radii to cover the range of orbits of stars within the satellite. In this way we show explicitly that prograde star orbits will be more easily stripped than radial orbits; while radial orbits are more easily stripped than retrograde ones. This result has previously been established by several authors numerically, but can now be understood analytically. For point mass, power-law (which includes the isothermal sphere), and a restricted class of split power-law potentials our solution is fully analytic. For more general potentials, we provide an equation which may be rapidly solved numerically.
Over short times (≲1–2 Gyr ∼1 satellite orbit), we find excellent agreement between our analytic and numerical models. Over longer times, star orbits within the satellite are transformed by the tidal field of the host galaxy. In a Hubble time, this causes a convergence of the three limiting tidal radii towards the prograde stripping radius. Beyond the prograde stripping radius, the velocity dispersion will be tangentially anisotropic.     

The SuperCOSMOS Sky Survey – III. Astrometry

In this, the third in a series of three papers concerning the SuperCOSMOS Sky Survey, we describe the astrometric properties of the data base. We describe the algorithms employed in the derivation of the astrometric parameters of the data, and demonstrate their accuracies by comparison with external data sets using the first release of data, the South Galactic Cap survey. We show that the celestial coordinates, which are tied to the International Celestial Reference Frame via the Tycho–2 reference catalogue, are accurate to better than ±0.2 arcsec at J , R ∼19,18 , rising to ±0.3 arcsec at J , R ∼22,21 , with positional-dependent systematic effects from bright to faint magnitudes at the ∼0.1-arcsec level. The proper motion measurements are shown to be accurate to typically ±10 mas yr⁻¹ at J , R ∼19,18 , rising to ±50 mas yr⁻¹ at J , R ∼22,21 , and are tied to zero using the extragalactic reference frame. We show that the zero-point errors in the proper motions are ≤1 mas yr⁻¹ for R >17 , and are no larger than ∼10 mas yr⁻¹ for R <17 mas yr⁻¹ .     

On water ice formation in interstellar clouds

A model is proposed for the formation of water ice mantles on grains in interstellar clouds. This occurs by direct accretion of monomers from the gas, be they formed by gas or surface reactions. The formation of the first monolayer requires a minimum extinction of interstellar radiation, sufficient to lower the grain temperature to the point where thermal evaporation of monomers is just offset by monomer accretion from the gas. This threshold is mainly determined by the adsorption energy of water molecules on the grain material; for hydrocarbon material, chemical simulation places this energy between 0.5 and 2 kcal mol⁻¹, which sets the (true) visible extinction threshold at a few magnitudes. However, realistic distributions of matter in a cloud will usually add to this an unrelated amount of cloud core extinction, which can explain the large dispersion of observed (apparent) thresholds. Once the threshold is crossed, all available water molecules in the gas are quickly adsorbed, because the grain cools down and the adsorption energy on ice is higher than on bare grain. The relative thickness of the mantle, and, hence, the slope of  τ₃( A _v)  depend only on the available water vapour, which is a small fraction of the oxygen abundance. Chemical simulation was also used to determine the adsorption sites and energies of O and OH on hydrocarbons and study the dynamics of formation of water molecules by surface reactions with gaseous H atoms, as well as their chances to stick in situ.     

On the relation between orbital structure and observed bar morphology

We investigate the morphological relation between the orbits of the central family of periodic orbits ( x ₁ family) and the bar itself using models of test particles moving in a barred potential. We show that different bar morphologies may have as a backbone the same set of x ₁ periodic orbits. We point out that by populating initially axisymmetric stellar discs exponentially with test particles in circular, or almost circular motion, we may end up with a response bar which reveals a shape different in crucial details from that of the individual stable x ₁ orbits. For example, a bar model in which the x ₁ orbits are pure ellipses may have a much more complicated response morphology. This depends on the particular invariant curves around x ₁, which are populated in each model.