Collapse of a Molecular Cloud Core to Stellar Densities: The First Three-Dimensional Calculations

We present results from the first three-dimensional calculations ever to follow the collapse of a molecular cloud core (~ 10^{-18} g cm^{-3}) to stellar densities (> 0.01 g cm^{-3}). The calculations resolve structures over 7 orders of magnitude in spatial extent (~ 5000 AU - 0.1 R_\odot), and over 17 orders of magnitude in density contrast. With these calculations, we consider whether fragmentation to form a close binary stellar system can occur during the second collapse phase. We find that, if the quasistatic core that forms before the second collapse phase is dynamically unstable to the growth of non-axisymmetric perturbations, the angular momentum extracted from the central regions of the core, via gravitational torques, is sufficient to prevent fragmentation and the formation of a close binary during the subsequent second collapse.Comment: ApJ Letters, in press (will appear in Nov 20 issue; available from the ApJ Rapid Release web page). 7 pages, incl. 5 figures. Also available at http://www.mpia-hd.mpg.de/theory/bat

On Tits' Centre Conjecture for Fixed Point Subcomplexes

We give a short and uniform proof of a special case of Tits' Centre Conjecture using a theorem of J-P. Serre and a result from our earlier work. We consider fixed point subcomplexes $X^H$ of the building $X = X(G)$ of a connected reductive algebraic group $G$, where $H$ is a subgroup of $G$.Comment: 4 pages; minor changes, to appear in C. R. Acad. Sci. Paris Ser. I Mat

Simulations of Stellar Collisions Involving Pre-Main Sequence Stars

In this paper, we present the results of smoothed particle hydrodynamic (SPH) simulations of collisions between pre-main sequence stars and a variety of other kinds of stars. Simulations over a range of impact parameters and velocities were performed. We find that pre-main sequence stars tend to ``wrap themselves'' around their impactor. We discuss the probable evolutionary state of products of collisions between pre-main sequence stars and pre-main sequence, main sequence, giant branch, and compact stars. The nature of the collision product does not depend strongly on the impact parameter or the velocity of the collision.Comment: Accepted by Ap

Modelling circumstellar discs with 3D radiation hydrodynamics

We present results from combining a grid-based radiative transfer code with a Smoothed Particle Hydrodynamics code to produce a flexible system for modelling radiation hydrodynamics. We use a benchmark model of a circumstellar disc to determine a robust method for constructing a gridded density distribution from SPH particles. The benchmark disc is then used to determine the accuracy of the radiative transfer results. We find that the SED and the temperature distribution within the disc are sensitive to the representation of the disc inner edge, which depends critically on both the grid and SPH resolution. The code is then used to model a circumstellar disc around a T-Tauri star. As the disc adjusts towards equilibrium vertical motions in the disc are induced resulting in scale height enhancements which intercept radiation from the central star. Vertical transport of radiation enables these perturbations to influence the mid-plane temperature of the disc. The vertical motions decay over time and the disc ultimately reaches a state of simultaneous hydrostatic and radiative equilibrium.Comment: MNRAS accepted; 15 pages; 17 figures, 4 in colou

Stellar Encounters with Massive Star-Disk Systems

The dense, clustered environment in which massive stars form can lead to interactions with neighboring stars. It has been hypothesized that collisions and mergers may contribute to the growth of the most massive stars. In this paper we extend the study of star-disk interactions to explore encounters between a massive protostar and a less massive cluster sibling using the publicly available SPH code GADGET-2. Collisions do not occur in the parameter space studied, but the end state of many encounters is an eccentric binary with a semi-major axis ~ 100 AU. Disk material is sometimes captured by the impactor. Most encounters result in disruption and destruction of the initial disk, and periodic torquing of the remnant disk. We consider the effect of the changing orientation of the disk on an accretion driven jet, and the evolution of the systems in the presence of on-going accretion from the parent core.Comment: 11 pages, 10 figures, accepted to Ap

Astrometric signatures of self-gravitating protoplanetary discs

We use high resolution numerical simulations to study whether gravitational instabilities within circumstellar discs can produce astrometrically detectable motion of the central star. For discs with masses of M_disc = 0.1 M_star, which are permanantly stable against fragmentation, we find that the magnitude of the astrometric signal depends upon the efficiency of disc cooling. Short cooling times produce prominent filamentary spiral structures in the disc, and lead to stellar motions that are potentially observable with future high precision astrometric experiments. For a disc that is marginally unstable within radii of \~10 au, we estimate astrometric displacements of 10-100 microarcsec on decade timescales for a star at a distance of 100 pc. The predicted displacement is suppressed by a factor of several in more stable discs in which the cooling time exceeds the local dynamical time by an order of magnitude. We find that the largest contribution comes from material in the outer regions of the disc and hence, in the most pessimistic scenario, the stellar motions caused by the disc could confuse astrometric searches for low mass planets orbiting at large radii. They are, however, unlikely to present any complications in searches for embedded planets orbiting at small radii, relative to the disc size, or Jupiter mass planets or greater orbiting at large radii.Comment: 6 pages, 9 figures, accepted for publication in MNRA

Formation of Globular Clusters in Galaxy Mergers

We present a high-resolution simulation of globular cluster formation in a galaxy merger. For the first time in such a simulation, individual star clusters are directly identified and followed on their orbits. We quantitatively compare star formation in the merger to that in the unperturbed galaxies. The merging galaxies show a strong starburst, in sharp contrast to their isolated progenitors. Most star clusters form in the tidal features. With a mass range of $5\times10^{5}$--$5\times 10^{6} M_{\odot}$, they are identified as globular clusters. The merger remnant is an elliptical galaxy. Clusters with different mass or age have different radial distributions in the galaxy. Our results show that the high specific frequency and bimodal distribution of metallicity observed in elliptical galaxies are natural products of gas-rich mergers, supporting a merger origin for the ellipticals and their globular cluster systems.Comment: ApJL accepted, version with high quality color images can be found in http://research.amnh.org/~yuexing/astro-ph/0407248.pd

The effect of magnetic fields on star cluster formation

We examine the effect of magnetic fields on star cluster formation by performing simulations following the self-gravitating collapse of a turbulent molecular cloud to form stars in ideal MHD. The collapse of the cloud is computed for global mass-to-flux ratios of infinity, 20, 10, 5 and 3, that is using both weak and strong magnetic fields. Whilst even at very low strengths the magnetic field is able to significantly influence the star formation process, for magnetic fields with plasma beta < 1 the results are substantially different to the hydrodynamic case. In these cases we find large-scale magnetically-supported voids imprinted in the cloud structure; anisotropic turbulent motions and column density structure aligned with the magnetic field lines, both of which have recently been observed in the Taurus molecular cloud. We also find strongly suppressed accretion in the magnetised runs, leading to up to a 75% reduction in the amount of mass converted into stars over the course of the calculations and a more quiescent mode of star formation. There is also some indication that the relative formation efficiency of brown dwarfs is lower in the strongly magnetised runs due to the reduction in the importance of protostellar ejections.Comment: 16 pages, 9 figures, 8 very pretty movies, MNRAS, accepted. Version with high-res figures + movies available from http://www.astro.ex.ac.uk/people/dprice/pubs/mcluster/index.htm

Formation of the First Supermassive Black Holes

We consider the physical conditions under which supermassive black holes could have formed inside the first galaxies. Our SPH simulations indicate that metal-free galaxies with a virial temperature ~10^4 K and with suppressed H2 formation (due to an intergalactic UV background) tend to form a binary black hole system which contains a substantial fraction (>10%) of the total baryonic mass of the host galaxy. Fragmentation into stars is suppressed without substantial H2 cooling. Our simulations follow the condensation of ~5x10^6 M_sun around the two centers of the binary down to a scale of < 0.1pc. Low-spin galaxies form a single black hole instead. These early black holes lead to quasar activity before the epoch of reionization. Primordial black hole binaries lead to the emission of gravitational radiation at redshifts z>10 that would be detectable by LISA.Comment: 11 pages, 9 figures, revised version, ApJ in press (October 10, 2003