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

The thermodynamics of collapsing molecular cloud cores using smoothed particle hydrodynamics with radiative transfer

We present the results of a series of calculations studying the collapse of molecular cloud cores performed using a three-dimensional smoothed particle hydr odynamics code with radiative transfer in the flux-limited diffusion approximation. The opacities and specific heat capacities are identical for each calculation. However, we find that the temperature evolution during the simulations varies significantly when starting from different initial conditions. Even spherically-symmetric clouds with different initial densities show markedly different development. We conclude that simple barotropic equations of state like those used in some previous calculations provide at best a crude approximation to the thermal behaviour of the gas. Radiative transfer is necessary to obtain accurate temperatures.Comment: 8 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

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

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

Competitive accretion in embedded stellar cluster

We investigate the physics of gas accretion in young stellar clusters. Accretion in clusters is a dynamic phenomenon as both the stars and the gas respond to the same gravitational potential. Accretion rates are highly non-uniform with stars nearer the centre of the cluster, where gas densities are higher, accreting more than others. This competitive accretion naturally results in both initial mass segregation and a spectrum of stellar masses. Accretion in gas-dominated clusters is well modelled using a tidal-lobe radius instead of the commonly used Bondi-Hoyle accretion radius. This works as both the stellar and gas velocities are under the influence of the same gravitational potential and are thus comparable. The low relative velocity that results means that the tidal radius is smaller than the Bondi-Hoyle radius in these systems. In contrast, when the stars dominate the potential and are virialised, the Bondi-Hoyle radius is smaller than the tidal radius and thus Bondi-Hoyle accretion is a better fit to the accretion rates.Comment: 11 pages, 11 figures, MNRAS in pres