Simulating `Complex' Problems with Quantum Monte Carlo

We present a new quantum Monte Carlo algorithm suitable for generically complex problems, such as systems coupled to external magnetic fields or anyons in two spatial dimensions. We find that the choice of gauge plays a nontrivial role, and can be used to reduce statistical noise in the simulation. Furthermore, it is found that noise can be greatly reduced by approximate cancellations between the phases of the (gauge dependent) statistical flux and the external magnetic flux.Comment: Revtex, 11 pages. 3 postscript files for figures attache

Outcomes and Duration of Tidal Evolution in a Star-Planet-Moon System

We formulated tidal decay lifetimes for hypothetical moons orbiting extrasolar planets with both lunar and stellar tides. Previous work neglected the effect of lunar tides on planet rotation, and are therefore applicable only to systems in which the moon's mass is much less than that of the planet. This work, in contrast, can be applied to the relatively large moons that might be detected around newly-discovered Neptune-mass and super-Earth planets. We conclude that moons are more stable when the planet/moon systems are further from the parent star, the planets are heavier, or the parent stars are lighter. Inclusion of lunar tides allows for significantly longer lifetimes for a massive moon relative to prior formulations. We expect that the semi-major axis of the planet hosting the first detected exomoon around a G-type star is 0.4-0.6 AU and is 0.2-0.4 AU for an M-type star.Comment: Accepted for publication in ApJ, 19 pages, 19 figure

Dynamically-Driven Star Formation In Models Of NGC 7252

We present new dynamical models of the merger remnant NGC 7252 which include star formation simulated according to various phenomenological rules. By using interactive software to match our model with the observed morphology and gas velocity field, we obtain a consistent dynamical model for NGC 7252. In our models, this proto-elliptical galaxy formed by the merger of two similar gas-rich disk galaxies which fell together with an initial pericentric separation of ~2 disk scale lengths approximately 620 Myr ago. Results from two different star formation rules--- density-dependent and shock-induced--- show significant differences in star formation during and after the first passage. Shock-induced star formation yields a prompt and wide-spread starburst at the time of first passage, while density-dependent star formation predicts a more slowly rising and centrally concentrated starburst. A comparison of the distributions and ages of observed clusters with results of our simulations favors shock-induced mechanism of star formation in NGC 7252. We also present simulated color images of our model of NGC 7252, constructed by incorporating population synthesis with radiative transfer and dust attenuation. Overall the predicted magnitudes and colors of the models are consistent with observations, although the simulated tails are fainter and redder than observed. We suggest that a lack of star formation in the tails, reflected by the redder colors, is due to an incomplete description of star formation in our models rather than insufficient gas in the tails.Comment: 11 pages, 9 figures, to be published in MNRA

The stellar content of the infalling molecular clump G286.21+0.17

The early evolution during massive star cluster formation is still uncertain. Observing embedded clusters at their earliest stages of formation can provide insight into the spatial and temporal distribution of the stars and thus probe different star cluster formation models. We present near-infrared imaging of an 8'*13'(5.4pc*8.7pc) region around the massive infalling clump G286.21+0.17(also known as BYF73). The stellar content across the field is determined and photometry is derived in order to { obtain} stellar parameters for the cluster members. We find evidence for some sub-structure (on scales less than a pc diameter) within the region with apparently at least three different sub-clusters associated with the molecular clump based on differences in extinction and disk fractions. At the center of the clump we identify a deeply embedded sub-cluster. Near-infrared excess is detected for 39-44% in the two sub-clusters associated with molecular material and 27% for the exposed cluster. Using the disk excess as a proxy for age this suggests the clusters are very young. The current total stellar mass is estimated to be at least 200 Msun. The molecular core hosts a rich population of pre-main sequence stars. There is evidence for multiple events of star formation both in terms of the spatial distribution within the star forming region and possibly from the disk frequency.Comment: Submitted to A

A sparse octree gravitational N-body code that runs entirely on the GPU processor

We present parallel algorithms for constructing and traversing sparse octrees on graphics processing units (GPUs). The algorithms are based on parallel-scan and sort methods. To test the performance and feasibility, we implemented them in CUDA in the form of a gravitational tree-code which completely runs on the GPU.(The code is publicly available at: http://castle.strw.leidenuniv.nl/software.html) The tree construction and traverse algorithms are portable to many-core devices which have support for CUDA or OpenCL programming languages. The gravitational tree-code outperforms tuned CPU code during the tree-construction and shows a performance improvement of more than a factor 20 overall, resulting in a processing rate of more than 2.8 million particles per second.Comment: Accepted version. Published in Journal of Computational Physics. 35 pages, 12 figures, single colum

The HI Content of Compact Groups of Galaxies

The HI content of Hickson Compact Groups in the southern hemisphere is measured using data from the HI Parkes All Sky Survey (HIPASS), and dedicated observations using the narrowband filter on the Multibeam instrument on the Parkes telescope. The expected HI mass of these groups was estimated using the luminosity, diameter and morphological types of the member galaxies, calibrated from published data. Taking careful account of non-detection limits, the results show that the compact group population that has been detected by these observations has an HI content similar to that of galaxies in the reference field sample. The upper limits for the undetected groups lie within the normal range; improvement of these limits will require a large increase in sensitivity.Comment: 27 pages, 5 figures. Accepted for publication in PAS

Dissipationless Collapse of Spherical Protogalaxies and the Fundamental Plane

Following on from the numerical work of Capelato, de Carvalho & Carlberg (1995, 1997), where dissipationless merger simulations were shown to reproduce the "Fundamental Plane" (FP) of elliptical galaxies, we investigate whether the end products of pure, spherically symmetric, one-component dissipationless {\it collapses} could also reproduce the FP. Past numerical work on collisionless collapses have addressed important issues on the dynamical/structural characteristics of collapsed equilibrium systems. However, the study of collisionless collapse in the context of the nature of the FP has not been satisfactorily addressed yet. Our aim in this paper is to focus our attention on the resulting collapse of simple one-component spherical models with a range of different initial virial coefficients. We find that the characteristic correlations of the models are compatible with virialized, centrally homologous systems. Our results strengthen the idea that merging may be a fundamental ingredient in forming non-homologous objects.Comment: 9 pages, 4 Postscript figures, Astronomy & Astrophysics in press (2002). Abstract placement correcte

Shock-induced star formation in a model of the Mice

Star formation plays an important role in the fate of interacting galaxies. To date, most galactic simulations including star formation have used a density-dependent star formation rule designed to approximate a Schmidt law. Here, I present a new star formation rule which is governed by the local rate of energy dissipation in shocks. The new and old rules are compared using self-consistent simulations of NGC 4676; shock-induced star formation provides a better match to the observations of this system.Comment: Submitted to MNRAS. 12 pages, 5 B&W figures, 2 color figures. See http://www.ifa.hawaii.edu/~barnes/sisfmm.html for a version with high-resolution figures, and http://www.ifa.hawaii.edu/~barnes/research/interaction_models/mice/index.html for animations of the simulation

Incorporating interactive 3-dimensional graphics in astronomy research papers

Most research data collections created or used by astronomers are intrinsically multi-dimensional. In contrast, all visual representations of data presented within research papers are exclusively 2-dimensional. We present a resolution of this dichotomy that uses a novel technique for embedding 3-dimensional (3-d) visualisations of astronomy data sets in electronic-format research papers. Our technique uses the latest Adobe Portable Document Format extensions together with a new version of the S2PLOT programming library. The 3-d models can be easily rotated and explored by the reader and, in some cases, modified. We demonstrate example applications of this technique including: 3-d figures exhibiting subtle structure in redshift catalogues, colour-magnitude diagrams and halo merger trees; 3-d isosurface and volume renderings of cosmological simulations; and 3-d models of instructional diagrams and instrument designs.Comment: 18 pages, 7 figures, submitted to New Astronomy. For paper with 3-dimensional embedded figures, see http://astronomy.swin.edu.au/s2plot/3dpd