BRIDGE: A Direct-tree Hybrid N-body Algorithm for Fully Self-consistent Simulations of Star Clusters and their Parent Galaxies

We developed a new direct-tree hybrid N-body algorithm for fully self-consistent N-body simulations of star clusters in their parent galaxies. In such simulations, star clusters need high accuracy, while galaxies need a fast scheme because of the large number of the particles required to model it. In our new algorithm, the internal motion of the star cluster is calculated accurately using the direct Hermite scheme with individual timesteps and all other motions are calculated using the tree code with second-order leapfrog integrator. The direct and tree schemes are combined using an extension of the mixed variable symplectic (MVS) scheme. Thus, the Hamiltonian corresponding to everything other than the internal motion of the star cluster is integrated with the leapfrog, which is symplectic. Using this algorithm, we performed fully self-consistent N-body simulations of star clusters in their parent galaxy. The internal and orbital evolutions of the star cluster agreed well with those obtained using the direct scheme. We also performed fully self-consistent N-body simulation for large-N models ($N=2\times 10^6$). In this case, the calculation speed was seven times faster than what would be if the direct scheme was used.Comment: 12 pages, 13 figures, Accepted for PAS

Evolution of Star Clusters near the Galactic Center: Fully Self-consistent N-body Simulations

We have performed fully self-consistent $N$-body simulations of star clusters near the Galactic center (GC). Such simulations have not been performed because it is difficult to perform fast and accurate simulations of such systems using conventional methods. We used the Bridge code, which integrates the parent galaxy using the tree algorithm and the star cluster using the fourth-order Hermite scheme with individual timestep. The interaction between the parent galaxy and the star cluster is calculate with the tree algorithm. Therefore, the Bridge code can handle both the orbital and internal evolutions of star clusters correctly at the same time. We investigated the evolution of star clusters using the Bridge code and compared the results with previous studies. We found that 1) the inspiral timescale of the star clusters is shorter than that obtained with "traditional" simulations, in which the orbital evolution of star clusters is calculated analytically using the dynamical friction formula and 2) the core collapse of the star cluster increases the core density and help the cluster survive. The initial conditions of star clusters is not so severe as previously suggested.Comment: 19 pages, 19 figures, accepted for publication in Ap

Evolution of Massive Black Hole Binaries

We present the result of large-scale N-body simulations of the stellar-dynamical evolution of a massive black-hole binary at the center of a spherical galaxy. We focus on the dependence of the hardening rate on the relaxation timescale of the parent galaxy. A simple theoretical argument predicts that a binary black hole creates the ``loss cone'' around it. Once the loss cone is formed, the hardening rate is determined by the rate at which field stars diffuse into the loss cone. Therefore the hardening timescale becomes proportional to the relaxation timescale. Recent N-body simulations, however, have failed to confirm this theory and various explanations have been proposed. By performing simulations with sufficiently large N (up to $10^6$) for sufficiently long time, we found that the hardening rate does depend on N. Our result is consistent with the simple theoretical prediction that the hardening timescale is proportional to the relaxation timescale. This dependence implies that most massive black hole binaries are unlikely to merge within the Hubble time through interaction with field stars and gravitational wave radiation alone.Comment: Reviced version accepted for publication in ApJ. Scheduled to appear in the February 10, 2004 issu

On relaxation processes in collisionless mergers

We analyze N-body simulations of halo mergers to investigate the mechanisms responsible for driving mixing in phase-space and the evolution to dynamical equilibrium. We focus on mixing in energy and angular momentum and show that mixing occurs in step-like fashion following pericenter passages of the halos. This makes mixing during a merger unlike other well known mixing processes such as phase mixing and chaotic mixing whose rates scale with local dynamical time. We conclude that the mixing process that drives the system to equilibrium is primarily a response to energy and angular momentum redistribution that occurs due to impulsive tidal shocking and dynamical friction rather than a result of chaotic mixing in a continuously changing potential. We also analyze the merger remnants to determine the degree of mixing at various radii by monitoring changes in radius, energy and angular momentum of particles. We confirm previous findings that show that the majority of particles retain strong memory of their original kinetic energies and angular momenta but do experience changes in their potential energies owing to the tidal shocks they experience during pericenter passages. Finally, we show that a significant fraction of mass (~ 40%) in the merger remnant lies outside its formal virial radius and that this matter is ejected roughly uniformly from all radii outside the inner regions. This highlights the fact that mass, in its standard virial definition, is not additive in mergers. We discuss the implications of these results for our understanding of relaxation in collisionless dynamical systems.Comment: Version accepted for Publication in Astrophysical Journal, March 20, 2007, v685. Minor changes, latex, 14 figure

On the effectiveness of mixing in violent relaxation

Relaxation processes in collisionless dynamics lead to peculiar behavior in systems with long-range interactions such as self-gravitating systems, non-neutral plasmas and wave-particle systems. These systems, adequately described by the Vlasov equation, present quasi-stationary states (QSS), i.e. long lasting intermediate stages of the dynamics that occur after a short significant evolution called "violent relaxation". The nature of the relaxation, in the absence of collisions, is not yet fully understood. We demonstrate in this article the occurrence of stretching and folding behavior in numerical simulations of the Vlasov equation, providing a plausible relaxation mechanism that brings the system from its initial condition into the QSS regime. Area-preserving discrete-time maps with a mean-field coupling term are found to display a similar behaviour in phase space as the Vlasov system.Comment: 10 pages, 11 figure

Sequential actions of Rab5 and Rab7 regulate endocytosis in the Xenopus oocyte

To explore the role of GTPases in endocytosis, we developed an assay using Xenopus oocytes injected with recombinant proteins to follow the uptake of the fluid phase marker HRP. HRP uptake was inhibited in cells injected with GTPγ S or incubated with aluminum fluoride, suggesting a general role for GTPases in endocytosis. Injection of Rab5 into oocytes, as well as Rab5:Q79L, a mutant with decreased GTPase activity, increased HRP uptake. Injection of Rab5:S34N, the dominant-negative mutant, inhibited HRP uptake. Injection of N-ethylmaleimide-sensitive factor (NSF) stimulated HRP uptake, and ATPase-defective NSF mutants inhibited HRP uptake when coinjected with Rab5:Q79L, confirming a requirement for NSF in endocytosis. Surprisingly, injection of Rab7:WT stimulated both uptake and degradation/activation of HRP. The latter appears to be due to enhanced transport to a late endosomal/prelysosomal degradative compartment that is monensin sensitive. Enhancement of uptake by Rab7 appears to function via an Rab5-sensitive pathway in oocytes since the stimulatory effect of Rab7 was blocked by coinjection of Rab5:S34N. Stimulation of uptake by Rab5 was blocked by Rab5:S34N but not by Rab7:T22N. Our results suggest that Rab7, while functioning downstream of Rab5, may be rate limiting for endocytosis in oocytes

The Orbit, Mass, and Albedo of Transneptunian Binary 1999 RZ253

We have observed 1999 RZ253 with the Hubble Space Telescope at seven separate epochs and have fit an orbit to the observed relative positions of this binary. Two orbital solutions have been identified that differ primarily in the inclination of the orbit plane. The best fit corresponds to an orbital period, P=46.263 +0.006/-0.074 days, semimajor axis a=4,660 +/-170 km and orbital eccentricity e=0.460 +/-0.013 corresponding to a system mass m=3.7 +/-0.4 x10^18 kg. For a density of rho = 1000 kg m^-3 the albedo at 477 nm is p = 0.12 +/-0.01, significantly higher than has been commonly assumed for objects in the Kuiper Belt. Multicolor, multiepoch photometry shows this pair to have colors typical for the Kuiper belt with a spectral gradient of 0.35 per 100 nm in the range between 475 and 775 nm. Photometric variations at the four epochs we observed were as large as 12 +/-3% but the sampling is insufficient to confirm the existence of a lightcurve

Towards a framework for validating the matching between notifications and scents in olfactory in-car interaction

Olfactory notifications have been proven to have a positive impact on drivers. This has motivated the use of scents to convey driving-relevant information. Research has proposed the use of such scents as lemon, peppermint, lavender and rose for in-car notifications. However, there is no framework to identify which scent is the most suitable for every application scenario. In this paper, we propose an approach for validating a matching between scents and driving-relevant notifications. We suggest a study in which the olfactory modality is compared with a puff of clean air, visual, auditory, and tactile stimuli while performing the same driving task. For the data analysis, we suggest recording the lane deviation, speed, time required to recover from the error, as well as the perceived liking and comfort ratings. Our approach aims to help automotive UI designers make better decisions about choosing the most suitable scent, as well as possible alternative modalities

The Origin of the Brightest Cluster Galaxies

Most clusters and groups of galaxies contain a giant elliptical galaxy in their centres which far outshines and outweighs normal ellipticals. The origin of these brightest cluster galaxies is intimately related to the collapse and formation of the cluster. Using an N-body simulation of a cluster of galaxies in a hierarchical cosmological model, we show that galaxy merging naturally produces a massive, central galaxy with surface brightness and velocity dispersion profiles similar to observed BCG's. To enhance the resolution of the simulation, 100 dark halos at $z=2$ are replaced with self-consistent disk+bulge+halo galaxy models following a Tully-Fisher relation using 100000 particles for the 20 largest galaxies and 10000 particles for the remaining ones. This technique allows us to analyze the stellar and dark matter components independently. The central galaxy forms through the merger of several massive galaxies along a filament early in the cluster's history. Galactic cannibalism of smaller galaxies through dynamical friction over a Hubble time only accounts for a small fraction of the accreted mass. The galaxy is a flattened, triaxial object whose long axis aligns with the primordial filament and the long axis of the cluster galaxy distribution agreeing with observed trends for galaxy-cluster alignment.Comment: Revised and accepted in ApJ, 25 pages, 10 figures, online version available at http://www.cita.utoronto.ca/~dubinski/bcg

Relaxation of a Collisionless System and the Transition to a New Equilibrium Velocity Distribution

In this paper, we present our conclusions from the numerical study of the collapse of a destabilized collisionless stellar system. We use both direct integration of the Vlasov-Poisson equations and an N-body tree code to obtain our results, which are mutually confirmed. We find that spherical and moderately nonspherical collapse configurations evolve to new equilibrium configurations in which the velocity distribution approaches a Gaussian form, at least in the central regions. The evolution to this state has long been an open question, and in this work we are able to clarify the process responsible and to support predictions made from statistical considerations (Lynden-Bell 1967; Nakamura 2000). The simulations of merging N-body systems show a transition to a Gaussian velocity distribution that is increasingly suppressed as the initial separation of centres is increased. Possible reasons for this are discussed.Comment: 25 pages, LaTeX. Accepted for publication in Ap