Showing 1–50 of 79 results for author: Huet, P

Light Curves and Colors of the Ejecta from Dimorphos after the DART Impact

Authors: Ariel Graykowski, Ryan A. Lambert, Franck Marchis, Dorian Cazeneuve, Paul A. Dalba, Thomas M. Esposito, Daniel O'Conner Peluso, Lauren A. Sgro, Guillaume Blaclard, Antonin Borot, Arnaud Malvache, Laurent Marfisi, Tyler M. Powell, Patrice Huet, Matthieu Limagne, Bruno Payet, Colin Clarke, Susan Murabana, Daniel Chu Owen, Ronald Wasilwa, Keiichi Fukui, Tateki Goto, Bruno Guillet, Patrick Huth, Satoshi Ishiyama , et al. (19 additional authors not shown)

Abstract: On 26 September 2022 the Double Asteroid Redirection Test (DART) spacecraft impacted Dimorphos, a satellite of the asteroid 65803 Didymos. Because it is a binary system, it is possible to determine how much the orbit of the satellite changed, as part of a test of what is necessary to deflect an asteroid that might threaten Earth with an impact. In nominal cases, pre-impact predictions of the orbit… ▽ More On 26 September 2022 the Double Asteroid Redirection Test (DART) spacecraft impacted Dimorphos, a satellite of the asteroid 65803 Didymos. Because it is a binary system, it is possible to determine how much the orbit of the satellite changed, as part of a test of what is necessary to deflect an asteroid that might threaten Earth with an impact. In nominal cases, pre-impact predictions of the orbital period reduction ranged from ~8.8 - 17.2 minutes. Here we report optical observations of Dimorphos before, during and after the impact, from a network of citizen science telescopes across the world. We find a maximum brightening of 2.29 $\pm$ 0.14 mag upon impact. Didymos fades back to its pre-impact brightness over the course of 23.7 $\pm$ 0.7 days. We estimate lower limits on the mass contained in the ejecta, which was 0.3 - 0.5% Dimorphos' mass depending on the dust size. We also observe a reddening of the ejecta upon impact. △ Less

Submitted 9 March, 2023; originally announced March 2023.

Comments: Accepted by Nature

Citizen Science Astronomy with a Network of Small Telescope: The Launch and Deployment of JWST

Authors: R. A. Lambert, F. Marchis, F., J. Asencio, G. Blaclard, L. A. Sgro, J. D. Giorgini, P. Plavchan, T. White, A. Verveen, T. Goto, P. Kuossari, N. Sethu, M. A. Loose, S. Will, K. Sibbernsen, J. W. Pickering, J. Randolph, K. Fukui, P. Huet, B. Guillet, O. Clerget, S. Stahl, N. Yoblonsky, M. Lauvernier , et al. (32 additional authors not shown)

Abstract: We present a coordinated campaign of observations to monitor the brightness of the James Webb Space Telescope (JWST) as it travels toward the second Earth-Sun Lagrange point and unfolds using the network ofUnistellar digital telescopes. Those observations collected by citizen astronomers across the world allowed us to detect specific phases such as the separation from the booster, glare due to a c… ▽ More We present a coordinated campaign of observations to monitor the brightness of the James Webb Space Telescope (JWST) as it travels toward the second Earth-Sun Lagrange point and unfolds using the network ofUnistellar digital telescopes. Those observations collected by citizen astronomers across the world allowed us to detect specific phases such as the separation from the booster, glare due to a change of orientation after a maneuver, the unfurling of the sunshield, and deployment of the primary mirror. After deployment of the sunshield on January 6 2022, the 6-h lightcurve has a significant amplitude and shows small variations due to the artificial rotation of the space telescope during commissionning. These variations could be due to the deployment of the primary mirror or some changes in orientation of the space telescope. This work illustrates the power of a worldwide array of small telescopes, operated by citizen astronomers, to conduct large scientific campaigns over a long timeframe. In the future, our network and others will continue to monitor JWST to detect potential degradations to the space environment by comparing the evolution of the lightcurve. △ Less

Submitted 9 July, 2022; originally announced July 2022.

Comments: 15 pages, 13 figures and 2 tables, SPIE Ground-based and Airborne Telescopes IX, AS22 SPIE Astronomical Telescopes + Instrumentation, 12182-144

The MICA Experiment: Astrophysics in Virtual Worlds

Authors: S. G. Djorgovski, Piet Hut, Rob Knop, Giuseppe Longo, Steve McMillan, Enrico Vesperini, Ciro Donalek, Matthew Graham, Asish Mahabal, Franz Sauer, Charles White, Crista Lopes

Abstract: We describe the work of the Meta-Institute for Computational Astrophysics (MICA), the first professional scientific organization based in virtual worlds. MICA was an experiment in the use of this technology for science and scholarship, lasting from the early 2008 to June 2012, mainly using the Second Life and OpenSimulator as platforms. We describe its goals and activities, and our future plans. W… ▽ More We describe the work of the Meta-Institute for Computational Astrophysics (MICA), the first professional scientific organization based in virtual worlds. MICA was an experiment in the use of this technology for science and scholarship, lasting from the early 2008 to June 2012, mainly using the Second Life and OpenSimulator as platforms. We describe its goals and activities, and our future plans. We conducted scientific collaboration meetings, professional seminars, a workshop, classroom instruction, public lectures, informal discussions and gatherings, and experiments in immersive, interactive visualization of high-dimensional scientific data. Perhaps the most successful of these was our program of popular science lectures, illustrating yet again the great potential of immersive VR as an educational and outreach platform. While the members of our research groups and some collaborators found the use of immersive VR as a professional telepresence tool to be very effective, we did not convince a broader astrophysics community to adopt it at this time, despite some efforts; we discuss some possible reasons for this non-uptake. On the whole, we conclude that immersive VR has a great potential as a scientific and educational platform, as the technology matures and becomes more broadly available and accepted. △ Less

Submitted 28 January, 2013; originally announced January 2013.

Comments: 10 pages, 8 figures; invited paper for the refereed proc. of the SLACTIONS 2012 conference

Journal ref: Ref. proc. of the SLACTIONS 2012 International Research Conference on Virtual Worlds,ceds. L. Morgado, Y. Sivan, A.M. Maia, et al., pp. 49-58 (2012). Vila Real, PT: UTAD, ISBN 978-989-704-102-0, avail. at http://hdl.handle.net/10348/2305

Few-Body Modes of Binary Formation in Core Collapse

Authors: Ataru Tanikawa, Douglas C. Heggie, Piet Hut, Junichiro Makino

Abstract: At the moment of deepest core collapse, a star cluster core contains less than ten stars. This small number makes the traditional treatment of hard binary formation, assuming a homogeneous background density, suspect. In a previous paper, we have found that indeed the conventional wisdom of binary formation, based on three-body encounters, is incorrect. Here we refine that insight, by further diss… ▽ More At the moment of deepest core collapse, a star cluster core contains less than ten stars. This small number makes the traditional treatment of hard binary formation, assuming a homogeneous background density, suspect. In a previous paper, we have found that indeed the conventional wisdom of binary formation, based on three-body encounters, is incorrect. Here we refine that insight, by further dissecting the subsequent steps leading to hard binary formation. For this purpose, we add some analysis tools in order to make the study less subjective. We find that the conventional treatment does remain valid for direct three-body scattering, but fails for resonant three-body scattering. Especially democratic resonance scattering, which forms an important part of the analytical theory of three-body binary formation, takes too much space and time to be approximated as being isolated, in the context of a cluster core around core collapse. We conclude that, while three-body encounters can be analytically approximated as isolated, subsequent strong perturbations typically occur whenever those encounters give rise to democratic resonances. We present analytical estimates postdicting our numerical results. If we only had been a bit more clever, we could have predicted this qualitative behaviour. △ Less

Submitted 20 November, 2013; v1 submitted 20 August, 2012; originally announced August 2012.

Comments: 44 pages, 18 figures, 1 table

PSDF: Particle Stream Data Format for N-Body Simulations

Authors: Will M. Farr, Jeff Ames, Piet Hut, Junichiro Makino, Steve McMillan, Takayuki Muranushi, Koichi Nakamura, Keigo Nitadori, Simon Portegies Zwart

Abstract: We present a data format for the output of general N-body simulations, allowing the presence of individual time steps. By specifying a standard, different N-body integrators and different visualization and analysis programs can all share the simulation data, independent of the type of programs used to produce the data. Our Particle Stream Data Format, PSDF, is specified in YAML, based on the same… ▽ More We present a data format for the output of general N-body simulations, allowing the presence of individual time steps. By specifying a standard, different N-body integrators and different visualization and analysis programs can all share the simulation data, independent of the type of programs used to produce the data. Our Particle Stream Data Format, PSDF, is specified in YAML, based on the same approach as XML but with a simpler syntax. Together with a specification of PSDF, we provide background and motivation, as well as specific examples in a variety of computer languages. We also offer a web site from which these examples can be retrieved, in order to make it easy to augment existing codes in order to give them the option to produce PSDF output. △ Less

Submitted 9 January, 2012; originally announced January 2012.

Comments: 5 pages; submitted to New Astronomy

Journal ref: New Astronomy 17 (2012), pp. 520--523

Unexpected Formation Modes of the First Hard Binary in Core Collapse

Authors: Ataru Tanikawa, Piet Hut, Junichiro Makino

Abstract: The conventional wisdom for the formation of the first hard binary in core collapse is that three-body interactions of single stars form many soft binaries, most of which are quickly destroyed, but eventually one of them survives. We report on direct N-body simulations to test these ideas, for the first time. We find that both assumptions are often incorrect: 1) quite a few three-body interactions… ▽ More The conventional wisdom for the formation of the first hard binary in core collapse is that three-body interactions of single stars form many soft binaries, most of which are quickly destroyed, but eventually one of them survives. We report on direct N-body simulations to test these ideas, for the first time. We find that both assumptions are often incorrect: 1) quite a few three-body interactions produce a hard binary from scratch; 2) and in many cases there are more than three bodies directly and simultaneously involved in the production of the first binary. The main reason for the discrepancies is that the core of a star cluster, at the first deep collapse, contains typically only five or so stars. Therefore, the homogeneous background assumption, which still would be reasonable for, say, 25 stars, utterly breaks down. There have been some speculations in this direction, but we demonstrate this result here explicitly, for the first time. △ Less

Submitted 19 July, 2011; originally announced July 2011.

Comments: 15 pages, 11 figures, 1 table, submitted to New Astronomy

AstroSim: Collaborative Visualization of an Astrophysics Simulation in Second Life

Authors: Arturo Nakasone, Helmut Prendinger, Simon Holland, Piet Hut, Jun Makino, Kenichi Miura

Abstract: We introduce AstroSim, a Second Life based prototype application for synchronous collaborative visualization targeted at astronomers. We introduce AstroSim, a Second Life based prototype application for synchronous collaborative visualization targeted at astronomers. △ Less

Submitted 30 November, 2009; originally announced November 2009.

Comments: For more information about AstroSim, visit: http://www.prendingerlab.net/globallab/?page_id=17. You can download the source code of the project at: http://sourceforge.net/projects/astrosim

Journal ref: IEEE Comp.Graph.App.29:69-81,2009

Exploring the Use of Virtual Worlds as a Scientific Research Platform: The Meta-Institute for Computational Astrophysics (MICA)

Authors: S. G. Djorgovski, P. Hut, S. McMillan, E. Vesperini, R. Knop, W. Farr, M. J. Graham

Abstract: We describe the Meta-Institute for Computational Astrophysics (MICA), the first professional scientific organization based exclusively in virtual worlds (VWs). The goals of MICA are to explore the utility of the emerging VR and VWs technologies for scientific and scholarly work in general, and to facilitate and accelerate their adoption by the scientific research community. MICA itself is an exp… ▽ More We describe the Meta-Institute for Computational Astrophysics (MICA), the first professional scientific organization based exclusively in virtual worlds (VWs). The goals of MICA are to explore the utility of the emerging VR and VWs technologies for scientific and scholarly work in general, and to facilitate and accelerate their adoption by the scientific research community. MICA itself is an experiment in academic and scientific practices enabled by the immersive VR technologies. We describe the current and planned activities and research directions of MICA, and offer some thoughts as to what the future developments in this arena may be. △ Less

Submitted 20 July, 2009; originally announced July 2009.

Comments: 15 pages, to appear in the refereed proceedings of "Facets of Virtual Environments" (FaVE 2009), eds. F. Lehmann-Grube, J. Sablating, et al., ICST Lecture Notes Ser., Berlin: Springer Verlag (2009); version with full resolution color figures is available at http://www.mica-vw.org/wiki/index.php/Publications

An Experiment in Using Virtual Worlds for Scientific Visualization of Self-Gravitating Systems

Authors: Will M. Farr, Piet Hut, Jeff Ames, Adam Johnson

Abstract: In virtual worlds, objects fall straight down. By replacing a few lines of code to include Newton's gravity, virtual world software can become an N-body simulation code with visualization included where objects move under their mutual gravitational attraction as stars in a cluster. We report on our recent experience of adding a gravitational N-body simulator to the OpenSim virtual world physics… ▽ More In virtual worlds, objects fall straight down. By replacing a few lines of code to include Newton's gravity, virtual world software can become an N-body simulation code with visualization included where objects move under their mutual gravitational attraction as stars in a cluster. We report on our recent experience of adding a gravitational N-body simulator to the OpenSim virtual world physics engine. OpenSim is an open-source, virtual world server that provides a 3D immersive experience to users who connect using the popular "Second Life" client software from Linden Labs. With the addition of the N-body simulation engine, which we are calling NEO, short for N-Body Experiments in OpenSim, multiple users can collaboratively create point-mass gravitating objects in the virtual world and then observe the subsequent gravitational evolution of their "stellar" system. We view this work as an experiment examining the suitability of virtual worlds for scientific visualization, and we report on future work to enhance and expand the prototype we have built. We also discuss some standardization and technology issues raised by our unusual use of virtual worlds. △ Less

Submitted 28 September, 2009; v1 submitted 7 May, 2009; originally announced May 2009.

Comments: 10 pages, 4 figures. As accepted by Journal of Virtual Worlds Research

Gravitational N-body Simulations

Authors: M. Trenti, P. Hut

Abstract: Gravitational N-body simulations, that is numerical solutions of the equations of motions for N particles interacting gravitationally, are widely used tools in astrophysics, with applications from few body or solar system like systems all the way up to galactic and cosmological scales. In this article we present a summary review of the field highlighting the main methods for N-body simulations a… ▽ More Gravitational N-body simulations, that is numerical solutions of the equations of motions for N particles interacting gravitationally, are widely used tools in astrophysics, with applications from few body or solar system like systems all the way up to galactic and cosmological scales. In this article we present a summary review of the field highlighting the main methods for N-body simulations and the astrophysical context in which they are usually applied. △ Less

Submitted 24 June, 2008; originally announced June 2008.

Comments: 13 pages, invited refereed review for the Scholarpedia Encyclopedia of Astrophysics, advanced undergraduate level. Available online at http://www.scholarpedia.org/article/N-body_simulations

Virtual Laboratories and Virtual Worlds

Authors: Piet Hut

Abstract: Since we cannot put stars in a laboratory, astrophysicists had to wait till the invention of computers before becoming laboratory scientists. For half a century now, we have been conducting experiments in our virtual laboratories. However, we ourselves have remained behind the keyboard, with the screen of the monitor separating us from the world we are simulating. Recently, 3D on-line technology… ▽ More Since we cannot put stars in a laboratory, astrophysicists had to wait till the invention of computers before becoming laboratory scientists. For half a century now, we have been conducting experiments in our virtual laboratories. However, we ourselves have remained behind the keyboard, with the screen of the monitor separating us from the world we are simulating. Recently, 3D on-line technology, developed first for games but now deployed in virtual worlds like Second Life, is beginning to make it possible for astrophysicists to enter their virtual labs themselves, in virtual form as avatars. This has several advantages, from new possibilities to explore the results of the simulations to a shared presence in a virtual lab with remote collaborators on different continents. I will report my experiences with the use of Qwaq Forums, a virtual world developed by a new company (see http://www.qwaq.com) △ Less

Submitted 11 December, 2007; originally announced December 2007.

Comments: 10 pages, 2 figures, Conference proceedings for IAUS246 'Dynamical Evolution of Dense Stellar Systems', ed. E. Vesperini (Chief Editor), M. Giersz, A. Sills, Capri, Sept. 2007

Modeling Dense Stellar Systems

Authors: Piet Hut, Shin Mineshige, Douglas C. Heggie, Junichiro Makino

Abstract: Black holes and neutron stars present extreme forms of matter that cannot be created as such in a laboratory on Earth. Instead, we have to observe and analyze the experiments that are ongoing in the Universe. The most telling observations of black holes and neutron stars come from dense stellar systems, where stars are crowded close enough to each other to undergo frequent interactions. It is th… ▽ More Black holes and neutron stars present extreme forms of matter that cannot be created as such in a laboratory on Earth. Instead, we have to observe and analyze the experiments that are ongoing in the Universe. The most telling observations of black holes and neutron stars come from dense stellar systems, where stars are crowded close enough to each other to undergo frequent interactions. It is the interplay between black holes, neutron stars and other objects in a dense environment that allows us to use observations to draw firm conclusions about the properties of these extreme forms of matter, through comparisons with simulations. The art of modeling dense stellar systems through computer simulations forms the main topic of this review. △ Less

Submitted 29 July, 2007; originally announced July 2007.

Comments: 23 pages, to appear in Prog. Theor. Phys

Journal ref: Prog.Theor.Phys.118:187-209,2007

Predictions for Triple Stars with and without a Pulsar in Star Clusters

Authors: M. Trenti, S. Ransom, P. Hut, D. C. Heggie

Abstract: Though about 80 pulsar binaries have been detected in globular clusters so far, no pulsar has been found in a triple system in which all three objects are of comparable mass. Here we present predictions for the abundance of such triple systems, and for the most likely characteristics of these systems. Our predictions are based on an extensive set of more than 500 direct simulations of star clust… ▽ More Though about 80 pulsar binaries have been detected in globular clusters so far, no pulsar has been found in a triple system in which all three objects are of comparable mass. Here we present predictions for the abundance of such triple systems, and for the most likely characteristics of these systems. Our predictions are based on an extensive set of more than 500 direct simulations of star clusters with primordial binaries, and a number of additional runs containing primordial triples. Our simulations employ a number N_{tot} of equal mass stars from N_{tot}=512 to N_{tot}=19661 and a primordial binary fraction from 0-50%. In addition, we validate our results against simulations with N=19661 that include a mass spectrum with a turn-off mass at 0.8 M_{sun}, appropriate to describe the old stellar populations of galactic globular clusters. Based on our simulations, we expect that typical triple abundances in the core of a dense cluster are two orders of magnitude lower than the binary abundances, which in itself already suggests that we don't have to wait too long for the first comparable-mass triple with a pulsar to be detected. △ Less

Submitted 3 April, 2008; v1 submitted 29 May, 2007; originally announced May 2007.

Comments: 11 pages, minor changes to match MNRAS accepted version

The core radius of a star cluster containing a massive black hole

Authors: Douglas C. Heggie, Piet Hut, Shin Mineshige, Jun Makino, Holger Baumgardt

Abstract: (Simplified) We present a theoretical framework which establishes how the core radius of a star cluster varies with the mass of an assumed central black hole. Our result is that the ratio of core to half-mass radius varies as the 3/4 power of the ratio of the black hole to cluster mass, when the system is well relaxed. The theory compares favourably with a number of simulations of this problem,… ▽ More (Simplified) We present a theoretical framework which establishes how the core radius of a star cluster varies with the mass of an assumed central black hole. Our result is that the ratio of core to half-mass radius varies as the 3/4 power of the ratio of the black hole to cluster mass, when the system is well relaxed. The theory compares favourably with a number of simulations of this problem, which extend to black hole masses of order 10% of the cluster mass. Though strictly limited as yet to clusters with stars of equal mass, our conclusion strengthens the view that clusters with large core radii are the most promising candidates in which to find a massive black hole. △ Less

Submitted 12 April, 2007; v1 submitted 30 November, 2006; originally announced November 2006.

Comments: 4 pages, 2 figures, to appear in PASJ Letters, PASJ Vol. 59, No. 3. (to be issued on June 25, 2007). Minor revisions

Neutron Stars and Black Holes in Star Clusters

Authors: F. A. Rasio, H. Baumgardt, A. Corongiu, F. D'Antona, G. Fabbiano, J. M. Fregeau, K. Gebhardt, C. O. Heinke, P. Hut, N. Ivanova, T. J. Maccarone, S. M. Ransom, N. A. Webb

Abstract: This article was co-authored by all invited speakers at the Joint Discussion on `Neutron Stars and Black Holes in Star Clusters,' which took place during the IAU General Assembly in Prague, Czech Republic, on August 17 and 18, 2006. Each section presents a short summary of recent developments in a key area of research, incorporating the main ideas expressed during the corresponding panel discuss… ▽ More This article was co-authored by all invited speakers at the Joint Discussion on `Neutron Stars and Black Holes in Star Clusters,' which took place during the IAU General Assembly in Prague, Czech Republic, on August 17 and 18, 2006. Each section presents a short summary of recent developments in a key area of research, incorporating the main ideas expressed during the corresponding panel discussion at the meeting. Our meeting, which had close to 300 registered participants, was broadly aimed at the large community of astronomers around the world working on the formation and evolution of compact objects and interacting binary systems in dense star clusters, such as globular clusters and galactic nuclei. Great advances have occurred in this field during the past few years, including the introduction of fundamentally new theoretical paradigms for the formation and evolution of compact objects in binaries as well as countless new discoveries by astronomers that have challenged many accepted models. Some of the highlights include: a nearly complete census of all the millisecond pulsars in 47 Tuc; first detections of many new radio pulsars in other clusters, particularly Terzan 5; detailed studies of X-ray binary populations and their luminosity functions in many galaxies and extragalactic globular clusters; increasing evidence for intermediate-mass black holes in clusters and greatly improved theoretical understanding of their possible formation processes. △ Less

Submitted 20 November, 2006; originally announced November 2006.

Comments: 29 pages, proceedings of IAU JD06, to appear in Highlights of Astronomy, Volume 14, ed. K. A. van der Hucht

Star Clusters with Primordial Binaries: III. Dynamical Interaction between Binaries and an Intermediate Mass Black Hole

Authors: M. Trenti, E. Ardi, S. Mineshige, P. Hut

Abstract: We present the first study of the dynamical evolution of an isolated star cluster that combines a significant population of primordial binaries with the presence of a central black hole. We use equal-mass direct N-body simulations, with N ranging from 4096 to 16384 and a primordial binary ratio of 0-10%; the black hole mass is about one percent of the total mass of the cluster. The evolution of… ▽ More We present the first study of the dynamical evolution of an isolated star cluster that combines a significant population of primordial binaries with the presence of a central black hole. We use equal-mass direct N-body simulations, with N ranging from 4096 to 16384 and a primordial binary ratio of 0-10%; the black hole mass is about one percent of the total mass of the cluster. The evolution of the binary population is strongly influenced by the presence of the black hole, which gives the cluster a large core with a central density cusp. Starting from a variety of initial conditions (Plummer and King models), we first encounter a phase, that last approximately 10 half-mass relaxation times, in which binaries are disrupted faster compared to analogous simulations without a black hole. Subsequently, however, binary disruption slows down significantly, due to the large core size. The dynamical interplay between the primordial binaries and the black hole thus introduces new features with respect to the scenarios investigated so far, where the influence of the black hole and of the binaries have been considered separately. A large core to half mass radius ratio appears to be a promising indirect evidence for the presence of a intermediate-mass black hole in old globular clusters. △ Less

Submitted 11 October, 2006; originally announced October 2006.

Comments: 11 pages, 11 figures, accepted for publication in MNRAS

Modeling Dense Stellar Systems: Background

Authors: Piet Hut

Abstract: I provide some background about recent efforts made in modeling dense stellar systems, within the context of the MODEST initiative. During the last four years, we have seen more than fifteen MODEST workshops, with an attendance between twenty and a hundred participants, and topics ranging from very specialized discussions to rather general overviews. I provide some background about recent efforts made in modeling dense stellar systems, within the context of the MODEST initiative. During the last four years, we have seen more than fifteen MODEST workshops, with an attendance between twenty and a hundred participants, and topics ranging from very specialized discussions to rather general overviews. △ Less

Submitted 7 October, 2006; originally announced October 2006.

Comments: to appear in the proceedings for JD06, IAU Highlights of Astronomy, Volume 14, ed. K. van der Hucht

Virtual Laboratories

Authors: Piet Hut

Abstract: At the frontier of most areas in science, computer simulations play a central role. The traditional division of natural science into experimental and theoretical investigations is now completely outdated. Instead, theory, simulation, and experimentation form three equally essential aspects, each with its own unique flavor and challenges. Yet, education in computational science is still lagging f… ▽ More At the frontier of most areas in science, computer simulations play a central role. The traditional division of natural science into experimental and theoretical investigations is now completely outdated. Instead, theory, simulation, and experimentation form three equally essential aspects, each with its own unique flavor and challenges. Yet, education in computational science is still lagging far behind, and the number of text books in this area is minuscule compared to the many text books on theoretical and experimental science. As a result, many researchers still carry out simulations in a haphazard way, without properly setting up the computational equivalent of a well equipped laboratory. The art of creating such a virtual laboratory, while providing proper extensibility and documentation, is still in its infancy. A new approach is described here, Open Knowledge, as an extension of the notion of Open Source software. Besides open source code, manuals, and primers, an open knowledge project provides simulated dialogues between code developers, thus sharing not only the code, but also the motivations behind the code. △ Less

Submitted 7 October, 2006; originally announced October 2006.

Comments: to appear in Prog. Theor. Phys

Journal ref: Prog.Theor.Phys.Suppl.164:38-53,2007

Threat Mitigation: The Gravity Tractor

Authors: Russell Schweickart, Clark Chapman, Dan Durda, Piet Hut

Abstract: The Gravity Tractor (GT) is a fully controlled asteroid deflection concept using the mutual gravity between a robotic spacecraft and an asteroid to slowly accelerate the asteroid in the direction of the "hovering" spacecraft. Based on early warning, provided by ground tracking and orbit prediction, it would be deployed a decade or more prior to a potential impact. Ion engines would be utilized f… ▽ More The Gravity Tractor (GT) is a fully controlled asteroid deflection concept using the mutual gravity between a robotic spacecraft and an asteroid to slowly accelerate the asteroid in the direction of the "hovering" spacecraft. Based on early warning, provided by ground tracking and orbit prediction, it would be deployed a decade or more prior to a potential impact. Ion engines would be utilized for both the rendezvous with the asteroid and the towing phase. Since the GT does not dock with or otherwise physically contact the asteroid during the deflection process there is no requirement for knowledge of the asteroid's shape, composition, rotation state or other "conventional" characteristics. The GT would first reduce the uncertainty in the orbit of the asteroid via Earth tracking of its radio transponder while station keeping with the asteroid. If, after analysis of the more precise asteroid orbit a deflection is indeed indicated, the GT would "hover" above the surface of the asteroid in the direction of the required acceleration vector for a duration adequate to achieve the desired velocity change. The orbit of the asteroid is continuously monitored throughout the deflection process and the end state is known in real time. The performance envelope for the GT includes most NEOs which experience close gravitational encounters prior to impact and those below 150-200 meters in diameter on a direct Earth impact trajectory. △ Less

Submitted 15 August, 2006; originally announced August 2006.

Comments: Submitted to NASA Workshop on Near-Earth Objects, Vail, Colorado, June 2006

Threat Mitigation: The Asteroid Tugboat

Authors: Russell Schweickart, Clark Chapman, Dan Durda, Piet Hut

Abstract: The Asteroid Tugboat (AT) is a fully controlled asteroid deflection concept using a robotic spacecraft powered by a high efficiency, electric propulsion system (ion or plasma) which docks with and attaches to the asteroid, conducts preliminary operations, and then thrusts continuously parallel to the asteroid velocity vector until the desired velocity change is achieved. Based on early warning,… ▽ More The Asteroid Tugboat (AT) is a fully controlled asteroid deflection concept using a robotic spacecraft powered by a high efficiency, electric propulsion system (ion or plasma) which docks with and attaches to the asteroid, conducts preliminary operations, and then thrusts continuously parallel to the asteroid velocity vector until the desired velocity change is achieved. Based on early warning, provided by ground tracking and orbit prediction, it would be deployed a decade or more prior to a potential impact. On completion of the initial rendezvous with the near-Earth object (NEO) the AT would first reduce the uncertainty in the orbit of the asteroid via Earth tracking of its radio transponder while it is station keeping with the asteroid. If on analysis of tracking data a deflection is required the AT would execute a reconnaissance phase collecting and processing information about the physical characteristics of the asteroid to support subsequent operations. The AT would then dock at the appropriate pole (i.e. on the spin axis), attach to the asteroid surface, and initiate a NEO reorientation maneuver. Following completion of the NEO reorientation the AT would initiate the deflection phase by thrusting continuously parallel to the asteroid velocity vector until the resultant target orbit is achieved. The orbit of the asteroid is continuously monitored throughout the deflection process and the end state is known in real time. If one assumes a nuclear-electric propulsion (NEP) system similar to that formerly under development in the recently canceled Prometheus Program, the AT would be capable of deflecting threatening NEOs up to 800 meters in diameter or more. △ Less

Submitted 15 August, 2006; originally announced August 2006.

Comments: Submitted to NASA Workshop on Near-Earth Objects, Vail, Colorado, June 2006

Threat Characterization: Trajectory Dynamics

Authors: Russell Schweickart, Clark Chapman, Dan Durda, Bill Bottke, David Nesvorny, Piet Hut

Abstract: Given a primary interest in "mitigation of the potential hazard" of near-Earth objects impacting the Earth, the subject of characterization takes on an aspect not normally present when considering asteroids as abstract bodies. Many deflection concepts are interested in the classic geophysical characteristics of asteroids when considering the physical challenge of modifying their orbits in order… ▽ More Given a primary interest in "mitigation of the potential hazard" of near-Earth objects impacting the Earth, the subject of characterization takes on an aspect not normally present when considering asteroids as abstract bodies. Many deflection concepts are interested in the classic geophysical characteristics of asteroids when considering the physical challenge of modifying their orbits in order to cause them to subsequently miss an impact with Earth. Yet for all deflection concepts there are characteristics of the threat which overwhelm these traditional factors. For example, a close gravitational encounter with Earth some years or decades prior to impact can reduce the velocity change necessary for deflection by several orders of magnitude if the deflection precedes the close encounter (or encounters). Conversely this "benefit" comes at a "price"; a corresponding increase in the accuracy of tracking required to determine the probability of impact. Societal issues, both national and international, also characterize the NEO deflection process and these may strongly contend with the purely technical issues normally considered. Therefore critical factors not normally considered must be brought into play as one characterizes the threat of NEO impacts. △ Less

Submitted 15 August, 2006; originally announced August 2006.

Comments: Submitted to NASA Workshop on Near-Earth Objects, Vail, Colorado, June 2006

Dynamical Formation of Close Binaries in Globular Clusters II: Cataclysmic Variables

Authors: David Pooley, Piet Hut

Abstract: We answer the long-standing question of which production mechanism is responsible for the cataclysmic variables (CVs) in globular clusters. Arguments have been given that range from mostly primordial presence to a significant contribution of later dynamical formation in close stellar encounters. We conclude, based on a thorough analysis of a homogeneous {\it Chandra} data set, that the majority… ▽ More We answer the long-standing question of which production mechanism is responsible for the cataclysmic variables (CVs) in globular clusters. Arguments have been given that range from mostly primordial presence to a significant contribution of later dynamical formation in close stellar encounters. We conclude, based on a thorough analysis of a homogeneous {\it Chandra} data set, that the majority of CVs in globulars has a dynamical origin. △ Less

Submitted 1 May, 2006; originally announced May 2006.

Comments: submitted to ApJL; comments welcome

Journal ref: Astrophys.J.646:L143-L146,2006

A Time-Symmetric Block Time-Step Algorithm for N-Body Simulations

Authors: Junichiro Makino, Piet Hut, Murat Kaplan, Hasan Saygin

Abstract: The method of choice for integrating the equations of motion of the general N-body problem has been to use an individual time step scheme. For the sake of efficiency, block time steps have been the most popular, where all time step sizes are smaller than a maximum time step size by an integer power of two. We present the first successful attempt to construct a time-symmetric integration scheme,… ▽ More The method of choice for integrating the equations of motion of the general N-body problem has been to use an individual time step scheme. For the sake of efficiency, block time steps have been the most popular, where all time step sizes are smaller than a maximum time step size by an integer power of two. We present the first successful attempt to construct a time-symmetric integration scheme, based on block time steps. We demonstrate how our scheme shows a vastly better long-time behavior of energy errors, in the form of a random walk rather than a linear drift. Increasing the number of particles makes the improvement even more pronounced. △ Less

Submitted 26 May, 2006; v1 submitted 18 April, 2006; originally announced April 2006.

Journal ref: New Astron. 12 (2006) 124-133

Star Clusters with Primordial Binaries: II. Dynamical Evolution of Models in a Tidal Field

Authors: M. Trenti, D. C. Heggie, P. Hut

Abstract: [abridged] We extend our analysis of the dynamical evolution of simple star cluster models, in order to provide comparison standards that will aid in interpreting the results of more complex realistic simulations. We augment our previous primordial-binary simulations by introducing a tidal field, and starting with King models of different central concentrations. We present the results of N-body… ▽ More [abridged] We extend our analysis of the dynamical evolution of simple star cluster models, in order to provide comparison standards that will aid in interpreting the results of more complex realistic simulations. We augment our previous primordial-binary simulations by introducing a tidal field, and starting with King models of different central concentrations. We present the results of N-body calculations of the evolution of equal-mass models, starting with primordial binary fractions of 0 - 100 %, and N values from 512 to 16384. We also attempt to extrapolate some of our results to the larger number of particles that are necessary to model globular clusters. We characterize the steady-state `deuterium main sequence' phase in which primordial binaries are depleted in the core in the process of `gravitationally burning'. In this phase we find that the ratio of the core to half-mass radius, r_c/r_h, is similar to that measured for isolated systems. In addition to the generation of energy due to hardening and depletion of the primordial binary population, the overall evolution of the star clusters is driven by a competing process: the tidal disruption of the system. We find that the depletion of primordial binaries before tidal dissolution of the system is possible only if the initial number is below 0.05 N, in the case of a King model with W_0=7 and N=4096 (which is one of our longest living models). We compare our findings, obtained by means of direct N-body simulations but scaled, where possible, to larger N, with similar studies carried out by means of Monte Carlo methods. △ Less

Submitted 6 October, 2006; v1 submitted 17 February, 2006; originally announced February 2006.

Comments: 15 pages, 18 figures, matches MNRAS accepted version, some sections reorganized but no major changes

Journal ref: Mon.Not.Roy.Astron.Soc.374:344-356,2007

Star Clusters with Primordial Binaries: I. Dynamical Evolution of Isolated Models

Authors: D. C. Heggie, M. Trenti, P. Hut

Abstract: In order to interpret the results of complex realistic star cluster simulations, which rely on many simplifying approximations and assumptions, it is essential to study the behavior of even more idealized models, which can highlight the essential physical effects and are amenable to more exact methods. With this aim, we present the results of N-body calculations of the evolution of equal-mass mo… ▽ More In order to interpret the results of complex realistic star cluster simulations, which rely on many simplifying approximations and assumptions, it is essential to study the behavior of even more idealized models, which can highlight the essential physical effects and are amenable to more exact methods. With this aim, we present the results of N-body calculations of the evolution of equal-mass models, starting with primordial binary fractions of 0 - 100 %, with values of N ranging from 256 to 16384. This allows us to extrapolate the main features of the evolution to systems comparable in particle number with globular clusters. In this range, we find that the steady-state `deuterium main sequence' is characterized by a ratio of the core radius to half-mass radius that follows qualitatively the analytical estimate by Vesperini & Chernoff (1994), although the N dependence is steeper than expected. Interestingly, for an initial binary fraction f greater than 10%, the binary heating in the core during the post collapse phase almost saturates (becoming nearly independent of f), and so little variation in the structural properties is observed. Thus, although we observe a significantly lower binary abundance in the core with respect to the Fokker-Planck simulations by Gao et al. (1991), this is of little dynamical consequence. At variance with the study of Gao et al. (1991), we see no sign of gravothermal oscillations before 150 halfmass relaxation times. At later times, however, oscillations become prominent. We demonstrate the gravothermal nature of these oscillations. △ Less

Submitted 17 February, 2006; originally announced February 2006.

Comments: 14 pages, 22 figures, MNRAS accepted

Journal ref: Mon.Not.Roy.Astron.Soc.368:677-689,2006

Dense Stellar Systems as Laboratories for Fundamental Physics

Authors: Piet Hut

Abstract: Galactic nuclei and globular clusters act as laboratories in which nature experiments with normal stars, neutron stars and black holes, through collisions and through the formation of bound states, in the form of binaries. The main difference with the usual Earth-based laboratories is that we cannot control the experiments. Instead, we have no choice but to create virtual laboratories on Earth,… ▽ More Galactic nuclei and globular clusters act as laboratories in which nature experiments with normal stars, neutron stars and black holes, through collisions and through the formation of bound states, in the form of binaries. The main difference with the usual Earth-based laboratories is that we cannot control the experiments. Instead, we have no choice but to create virtual laboratories on Earth, in order to simulate all the relevant physics in large-scale computational experiments. This implies a realistic treatment of stellar dynamics, stellar evolution, and stellar hydrodynamics. Each of these three fields has its own legacy codes, workhorses that are routinely used to simulate star clusters, stars, and stellar collisions, respectively. I outline the main steps that need to be taken in order to embed and where needed transform these legacy codes in order to produce a far more modular and robust environment for modeling dense stellar systems. The time is right to do so: within a few years computers will reach the required speed, in the Petaflops range, to follow a star cluster with a million stars for ten billion years, while resolving the internal binary and multiple star motions. By that time simulation software will be the main bottleneck in our ability to analyze dense stellar systems. Only through full-scale simulations will we be able to critically test our understanding of the `microphysics' of stellar collisions and their aftermath, in a direct comparison with observations. △ Less

Submitted 11 January, 2006; originally announced January 2006.

Comments: review paper, presented at A Life With Stars (Conference in Honor of Ed van den Heuvel), Amsterdam, August, 2005

Three-Body Affairs in the Outer Solar System

Authors: Yoko Funato, Junichiro Makino, Piet Hut, Eiichiro Kokubo, Daisuke Kinoshita

Abstract: Recent observations (Burnes2002,Veillet2002,Margot2002a) have revealed an unexpectedly high binary fraction among the Trans-Neptunian Objects (TNOs) that populate the Kuiper Belt. The TNO binaries are strikingly different from asteroid binaries in four respects (Veillet2002): their frequency is an order of magnitude larger, the mass ratio of their components is closer to unity, and their orbits… ▽ More Recent observations (Burnes2002,Veillet2002,Margot2002a) have revealed an unexpectedly high binary fraction among the Trans-Neptunian Objects (TNOs) that populate the Kuiper Belt. The TNO binaries are strikingly different from asteroid binaries in four respects (Veillet2002): their frequency is an order of magnitude larger, the mass ratio of their components is closer to unity, and their orbits are wider and highly eccentric. Two explanations have been proposed for their formation, one assuming large numbers of massive bodies (Weidenschilling2002), and one assuming large numbers of light bodies (Goldreich2002). We argue that both assumptions are unwarranted, and we show how TNO binaries can be produced from a modest number of intermediate-mass bodies of the type predicted by the gravitational instability theory for the formation of planetesimals (Goldreich and Ward1973). We start with a TNO binary population similar to the asteroid binary population, but subsequently modified by three-body exchange reactions, a process that is far more efficient in the Kuiper belt, because of the much smaller tidal perturbations by the Sun. Our mechanism can naturally account for all four characteristics that distinguish TNO binaries from main-belt asteroid binaries. △ Less

Submitted 3 December, 2005; originally announced December 2005.

Comments: Published in 2003 in the proceedings of the 35th Symposium on Celestial Mechanics. 8 pages

Journal ref: In Proceedings of the 35th Symposium on Celestial Mechanics, eds. E. Kokubo, H. Arakida, and T. Yamamoto. Tokyo, Japan, 2003

The ecology of star clusters and intermediate mass black holes in the Galactic bulge

Authors: Simon Portegies Zwart, Holger Baumgardt, Stephen L. W. McMillan, Junichiro Makino, Piet Hut, Toshi Ebisuzaki

Abstract: We simulate the inner 100pc of the Milky-Way Galaxy to study the formation and evolution of the population of star clusters and intermediate mass black holes. For this study we perform extensive direct N-body simulations of the star clusters which reside in the bulge, and of the inner few tenth of parsecs of the super massive black hole in the Galactic center. In our N-body simulations the dynam… ▽ More We simulate the inner 100pc of the Milky-Way Galaxy to study the formation and evolution of the population of star clusters and intermediate mass black holes. For this study we perform extensive direct N-body simulations of the star clusters which reside in the bulge, and of the inner few tenth of parsecs of the super massive black hole in the Galactic center. In our N-body simulations the dynamical friction of the star cluster in the tidal field of the bulge are taken into account via (semi)analytic soluations. The N-body calculations are used to calibrate a (semi)analytic model of the formation and evolution of the bulge. We find that about 10% of the clusters born within 100pc of the Galactic center undergo core collapse during their inward migration and form intermediate-mass black holes (IMBHs) via runaway stellar merging. After the clusters dissolve, these IMBHs continue their inward drift, carrying a few of the most massive stars with them. We predict that region within about 10 parsec of the SMBH is populated by about 50IMBHs of some 1000Msun. Several of these are expected to be accompanied still by some of the most massive stars from the star cluster. We also find that within a few milliparsec of the SMBH there is a steady population of several IMBHs. This population drives the merger rate between IMBHs and the SMBH at a rate of about one per 10Myr, sufficient to build the accumulate majority of mass of the SMBH. Mergers of IMBHs with SMBHs throughout the universe are detectable by LISA, at a rate of about two per week. △ Less

Submitted 5 December, 2005; v1 submitted 14 November, 2005; originally announced November 2005.

Comments: 26 pages, ApJ in press. Old title: The Galactic center welcomes black hole immigrants

A New Time-Symmetric Block Time-Step Algorithm for N-Body Simulations

Authors: Murat Kaplan, Hasan Saygin, Piet Hut, Jun Makino

Abstract: Time-symmetric integration schemes share with symplectic schemes the property that their energy errors show a much better behavior than is the case for generic integration schemes. Allowing adaptive time steps typically leads to a loss of symplecticity. In contrast, time symmetry can be easily maintained, at least for a continuous choice of time step size. In large-scale N-body simulations, howe… ▽ More Time-symmetric integration schemes share with symplectic schemes the property that their energy errors show a much better behavior than is the case for generic integration schemes. Allowing adaptive time steps typically leads to a loss of symplecticity. In contrast, time symmetry can be easily maintained, at least for a continuous choice of time step size. In large-scale N-body simulations, however, one often uses block time steps, where all time steps are forced to take on values as powers of two. This greatly facilitates parallelization, and hence code efficiency. Straightforward implementation of time-symmetry, translated to block time steps, faces significant hurdles. For example, iteration can lead to oscillatory behavior, and even when such behavior is suppressed, energy errors show a linear drift in time. We present an approach that circumvents these problems. △ Less

Submitted 10 November, 2005; originally announced November 2005.

Comments: 14 pages, 3 figures, submitted for publication in New Astronomy

Performance Tuning of N-Body Codes on Modern Microprocessors: I. Direct Integration with a Hermite Scheme on x86_64 Architecture

Authors: Keigo Nitadori, Junichiro Makino, Piet Hut

Abstract: The main performance bottleneck of gravitational N-body codes is the force calculation between two particles. We have succeeded in speeding up this pair-wise force calculation by factors between two and ten, depending on the code and the processor on which the code is run. These speedups were obtained by writing highly fine-tuned code for x86_64 microprocessors. Any existing N-body code, running… ▽ More The main performance bottleneck of gravitational N-body codes is the force calculation between two particles. We have succeeded in speeding up this pair-wise force calculation by factors between two and ten, depending on the code and the processor on which the code is run. These speedups were obtained by writing highly fine-tuned code for x86_64 microprocessors. Any existing N-body code, running on these chips, can easily incorporate our assembly code programs. In the current paper, we present an outline of our overall approach, which we illustrate with one specific example: the use of a Hermite scheme for a direct N^2 type integration on a single 2.0 GHz Athlon 64 processor, for which we obtain an effective performance of 4.05 Gflops, for double precision accuracy. In subsequent papers, we will discuss other variations, including the combinations of N log N codes, single precision implementations, and performance on other microprocessors. △ Less

Submitted 2 November, 2005; originally announced November 2005.

Comments: 32 pages, 2 figures

Journal ref: New Astron.12:169-181,2006

Primordial binaries and intermediate mass black holes in globular clusters

Authors: M. Trenti, E. Ardi, S. Mineshige, P. Hut

Abstract: We present the first study of the dynamical evolution of a star cluster that combines a significant population of primordial binaries with the presence of a central black hole. We use direct N-body simulations, with a black hole mass of about one percent of the total mass of the cluster. The evolution of the binary population is strongly influenced by the presence of the black hole, which give… ▽ More We present the first study of the dynamical evolution of a star cluster that combines a significant population of primordial binaries with the presence of a central black hole. We use direct N-body simulations, with a black hole mass of about one percent of the total mass of the cluster. The evolution of the binary population is strongly influenced by the presence of the black hole, which gives the cluster a large core with a central density cusp. Starting from a Plummer profile, we first encounter a phase, that last approximately 10 half-mass relaxation times, in which binaries are disrupted faster compared to analogous simulations without a black hole. Subsequently, however, binary disruption slows down significantly, due to the large core size. The dynamical interplay between the primordial binaries and the black hole thus introduces new features with respect to the scenarios investigated so far, where the influence of the black hole and of the binaries have been considered separately. Specifically, the pattern of binary destruction by an intermediate-mass black hole may leave a fingerprint that could be detected observationally. △ Less

Submitted 24 August, 2005; originally announced August 2005.

Comments: 4 pages, 2 figures, submitted to the Astrophysical Journal Letters, uses emulateapj.cls

Black Holes in Massive Star Clusters

Authors: Steve McMillan, Holger Baumgardt, Simon Portegies Zwart, Piet Hut, Junichiro Makino

Abstract: Close encounters and physical collisions between stars in young dense clusters can result in new channels for stellar evolution, and may lead to the formation of very massive stars and black holes via runaway merging. We present some details of this process, using the results of N-body simulations and simple analytical estimates to place limits on the cluster parameters for which it expected to… ▽ More Close encounters and physical collisions between stars in young dense clusters can result in new channels for stellar evolution, and may lead to the formation of very massive stars and black holes via runaway merging. We present some details of this process, using the results of N-body simulations and simple analytical estimates to place limits on the cluster parameters for which it expected to occur. For small clusters, the mass of the runaway is effectively limited by the total number of high-mass stars in the system. For larger clusters, the runaway mass is determined by the fraction of stars that can mass-segregate to the cluster core while still on the main sequence. In typical cases, the result is in the range commonly cited for intermediate-mass black holes. This mechanism may therefore have important implications for the formation of massive black holes and black-hole binaries in dense cluster cores. △ Less

Submitted 6 November, 2004; originally announced November 2004.

Comments: 8 pages, 2 figures; to appear in "Formation and Evolution of Massive Young Star Clusters," eds. H.J.G.L.M. Lamers, A. Nota & L.J. Smith

Which Globular Clusters contain Intermediate-mass Black Holes?

Authors: Holger Baumgardt, Junichiro Makino, Piet Hut

Abstract: It has been assumed that intermediate-mass black holes (IMBHs) in globular clusters can only reside in the most centrally concentrated clusters, with a so-called `core-collapsed' density profile. While this would be a natural guess, it is in fact wrong. We have followed the evolution of star clusters containing IMBHs with masses between 125 \le M_{BH} \le 1000 M_{\odot} through detailed N-body s… ▽ More It has been assumed that intermediate-mass black holes (IMBHs) in globular clusters can only reside in the most centrally concentrated clusters, with a so-called `core-collapsed' density profile. While this would be a natural guess, it is in fact wrong. We have followed the evolution of star clusters containing IMBHs with masses between 125 \le M_{BH} \le 1000 M_{\odot} through detailed N-body simulations, and we find that a cluster with an IMBH, in projection, appears to have a relatively large `core' with surface brightness only slightly rising toward the center. This makes it highly unlikely that any of `core-collapsed' clusters will harbor an IMBH. On the contrary, the places to look for an IMBH are those clusters that can be fitted well by medium-concentration King models. The velocity dispersion of the visible stars in a globular cluster with an IMBH is nearly constant well inside the apparent core radius. For a cluster of mass M_C containing an IMBH of mass M_{BH}, the influence of the IMBH becomes significant only at a fraction 2.5 M_{BH}/M_C of the half-mass radius, deep within the core, where it will affect only a small number of stars. In conclusion, observational detection of an IMBH may be possible, but will be challenging. △ Less

Submitted 25 October, 2004; originally announced October 2004.

Comments: 13 pages, 3 figures, Accepted ofr publication in ApJ (scheduled for February 2005)

Journal ref: Astrophys.J. 620 (2005) 238-243

The formation of massive black holes through collision runaway in dense young star clusters

Authors: Simon Portegies Zwart, Holger Baumgardt, Piet Hut, Junichiro Makino, Stephen McMillan

Abstract: A luminous X-ray source is associated with a cluster (MGG-11) of young stars \~200pc from the center of the starburst galaxy M82. The properties of the X-ray source are best explained by a black hole with a mass of at least 350Msun, which is intermediate between stellar-mass and supermassive black holes. A nearby but somewhat more massive star cluster (MGG-9) shows no evidence of such an interme… ▽ More A luminous X-ray source is associated with a cluster (MGG-11) of young stars \~200pc from the center of the starburst galaxy M82. The properties of the X-ray source are best explained by a black hole with a mass of at least 350Msun, which is intermediate between stellar-mass and supermassive black holes. A nearby but somewhat more massive star cluster (MGG-9) shows no evidence of such an intermediate mass black hole, raising the issue of just what physical characteristics of the clusters can account for this difference. Here we report numerical simulations of the evolution and the motions of stars within the clusters, where stars are allowed to mergers with each other. We find that for MGG-11 dynamical friction leads to the massive stars sinking rapidly to the center of the cluster to participate in a runaway collision, thereby producing a star of 800-3000Msun, which ultimately collapses to an black hole of intermediate mass. No such runaway occurs in the cluster MGG-9 because the larger cluster radius leads to a mass-segregation timescale a factor of five longer than for MGG-11. △ Less

Submitted 26 February, 2004; originally announced February 2004.

Comments: Accepted for publication in Nature (Including supplementary information)

Journal ref: Nature428:724,2004

The formation of Kuiper-belt Binaries through Exchange Reactions

Authors: Yoko Funato, Junichiro Makino, Piet Hut, Eiichiro Kokubo, Daisuke Kinoshita

Abstract: Recent observations have revealed an unexpectedly high binary fraction among the Trans-Neptunian Objects (TNOs) that populate the Kuiper-belt. The discovered binaries have four characteristics they comprise a few percent of the TNOs, the mass ratio of their components is close to unity, their internal orbits are highly eccentric, and the orbits are more than 100 times wider than the primary's ra… ▽ More Recent observations have revealed an unexpectedly high binary fraction among the Trans-Neptunian Objects (TNOs) that populate the Kuiper-belt. The discovered binaries have four characteristics they comprise a few percent of the TNOs, the mass ratio of their components is close to unity, their internal orbits are highly eccentric, and the orbits are more than 100 times wider than the primary's radius. In contrast, theories of binary asteroid formation tend to produce close, circular binaries. Therefore, a new approach is required to explain the unique characteristics of the TNO binaries. Two models have been proposed. Both, however, require extreme assumptions on the size distribution of TNOs. Here we show a mechanism which is guaranteed to produces binaries of the required type during the early TNO growth phase, based on only one plausible assumption, namely that initially TNOs were formed through gravitational instabilities of the protoplanetary dust layer. △ Less

Submitted 13 February, 2004; originally announced February 2004.

Comments: 12pages, 4 figures

MODEST: modeling stellar evolution and (hydro)dynamics

Authors: Piet Hut

Abstract: Simulations of dense stellar systems currently face two major hurdles, one astrophysical and one computational. The astrophysical problem lies in the fact that several major stages in binary evolution, such as common envelope evolution, are still poorly understood. The best we can do in these cases is to parametrize our ignorance, in a way that is reminiscent of the introduction of a mixing leng… ▽ More Simulations of dense stellar systems currently face two major hurdles, one astrophysical and one computational. The astrophysical problem lies in the fact that several major stages in binary evolution, such as common envelope evolution, are still poorly understood. The best we can do in these cases is to parametrize our ignorance, in a way that is reminiscent of the introduction of a mixing length to describe convection in a single star, or an alpha parameter in modeling an accretion disk. The hope is that by modeling a whole star cluster in great detail, and comparing the results to the wealth of observational data currently available, we will be able to constrain the parameters that capture the unknown physics. The computational problem is one of composition: while we have accurate computer codes for modeling stellar dynamics, stellar hydrodynamics, and stellar evolution, we currently have no good way to put all this knowledge together in a single software environment. A year ago, a loosely-knit organization was founded to address these problems, MODEST for MOdeling DEnse STellar systems, with nine working groups and a series of meetings that are held every half year. This report reviews the first year of this initiative. Much more detail can be found on the MODEST web site http://www.manybody.org/modest.html . △ Less

Submitted 15 September, 2003; originally announced September 2003.

Comments: 4 pages, review talk presented at JD-11, IAU General Assembly XXV, Sydney, July 2003, to be included in Highlights of Astronomy, Vol. 13

Dynamical Formation of Close Binary Systems in Globular Clusters

Authors: D. Pooley, W. H. G. Lewin, S. F. Anderson, H. Baumgardt, A. V. Filippenko, B. M. Gaensler, L. Homer, P. Hut, V. M. Kaspi, B. Margon, S. McMillan, S. Portegies Zwart, M. van der Klis, F. Verbunt

Abstract: We know from observations that globular clusters are very efficient catalysts in forming unusual short-period binary systems or their offspring, such as low-mass X-ray binaries (LMXBs; neutron stars accreting matter from low-mass stellar companions), cataclysmic variables (CVs; white dwarfs accreting matter from stellar companions), and millisecond pulsars (MSPs; rotating neutron stars with spin… ▽ More We know from observations that globular clusters are very efficient catalysts in forming unusual short-period binary systems or their offspring, such as low-mass X-ray binaries (LMXBs; neutron stars accreting matter from low-mass stellar companions), cataclysmic variables (CVs; white dwarfs accreting matter from stellar companions), and millisecond pulsars (MSPs; rotating neutron stars with spin periods of a few ms). Although there has been little direct evidence, the overabundance of these objects in globular clusters has been attributed by numerous authors to the high densities in the cores, which leads to an increase in the formation rate of exotic binary systems through close stellar encounters. Many such close binary systems emit X-radiation at low luminosities (L_x < 10^{34} erg/s) and are being found in large numbers through observations with the Chandra X-ray Observatory. Here we present conclusive observational evidence for a link between the number of close binaries observed in X-rays in a globular cluster and the stellar encounter rate of the cluster. We also make an estimate of the total number of LMXBs in globular clusters in our Galaxy. △ Less

Submitted 30 April, 2003; originally announced May 2003.

Comments: 11 pages, 1 b&w figure, 1 color figure

Journal ref: Astrophys.J.591:L131-L134,2003

Kuiper-belt Binary Formation through Exchange Reactions

Authors: Yoko Funato, Junichiro Makino, Piet Hut, Eiichiro Kokubo, Daisuke Kinoshita

Abstract: Recent observations (Burnes et al 2002,Veillet et al 2002, Margot et al 2002a) have revealed an unexpectedly high binary fraction among the Trans-Neptunian Objects (TNOs) that populate the Kuiper Belt. The TNO binaries are strikingly different from asteroid binaries in four respects (Veillet et al 2002): their frequency is an order of magnitude larger, the mass ratio of their components is close… ▽ More Recent observations (Burnes et al 2002,Veillet et al 2002, Margot et al 2002a) have revealed an unexpectedly high binary fraction among the Trans-Neptunian Objects (TNOs) that populate the Kuiper Belt. The TNO binaries are strikingly different from asteroid binaries in four respects (Veillet et al 2002): their frequency is an order of magnitude larger, the mass ratio of their components is closer to unity, and their orbits are wider and highly eccentric. Two explanations have been proposed for their formation, one assuming large numbers of massive bodies (Weidenschilling 2002), and one assuming large numbers of light bodies (Goldreich et al 2002). We argue that both assumptions are unwarranted, and we show how TNO binaries can be produced from a modest number of intermediate-mass bodies of the type predicted by the gravitational instability theory for the formation of planetesimals (Goldreich and Ward 1973). We start with a TNO binary population similar to the asteroid binary population, but subsequently modified by three-body exchange reactions, a process that is far more efficient in the Kuiper belt, because of the much smaller tidal perturbations by the Sun. Our mechanism can naturally account for all four characteristics that distinguish TNO binaries from main-belt asteroid binaries. △ Less

Submitted 28 March, 2003; v1 submitted 5 March, 2003; originally announced March 2003.

Comments: 11 pages, 3 Figures

Modest-2: A Summary

Authors: A. Sills, S. Deiters, P. Eggleton, M. Freitag, M. Giersz, D. Heggie, J. Hurley, P. Hut, N. Ivanova, R. Klessen, P. Kroupa, J. Lombardi, S. McMillan, S. Portegies Zwart, H. Zinnecker

Abstract: This is a summary paper of MODEST-2, a workshop held at the Astronomical Institute ``Anton Pannekoek'' in Amsterdam, 16-17 December 2002. MODEST is a loose collaboration of people interested in MOdelling DEnse STellar systems, particularly those interested in modelling these systems using all the available physics (stellar dynamics, stellar evolution, hydrodynamics and the interplay between the… ▽ More This is a summary paper of MODEST-2, a workshop held at the Astronomical Institute ``Anton Pannekoek'' in Amsterdam, 16-17 December 2002. MODEST is a loose collaboration of people interested in MOdelling DEnse STellar systems, particularly those interested in modelling these systems using all the available physics (stellar dynamics, stellar evolution, hydrodynamics and the interplay between the three) by defining interfaces between different codes. In this paper, we summarize 1) the main advances in this endeavour since MODEST-1; 2) the main science goals which can be and should be addressed by these types of simulations; and 3) the most pressing theoretical and modelling advances that we identified. △ Less

Submitted 11 March, 2003; v1 submitted 23 January, 2003; originally announced January 2003.

Comments: Accepted by New Astronomy

Journal ref: New Astron. 8 (2003) 605-628

A Dynamical Model for the Globular Cluster G1

Authors: Holger Baumgardt, Junichiro Makino, Piet Hut, Steve McMillan, Simon Portegies Zwart

Abstract: We present a comparison between the observational data on the kinematical structure of G1 in M31, obtained with the Hubble WFPC2 and STIS instruments, and the results of dynamical simulations carried out using the special-purpose computer GRAPE-6. We have obtained good fits for models starting with single cluster King-model initial conditions and even better fits when starting our simulations wi… ▽ More We present a comparison between the observational data on the kinematical structure of G1 in M31, obtained with the Hubble WFPC2 and STIS instruments, and the results of dynamical simulations carried out using the special-purpose computer GRAPE-6. We have obtained good fits for models starting with single cluster King-model initial conditions and even better fits when starting our simulations with a dynamically constructed merger remnant of two star clusters. In the latter case, the results from our simulations are in excellent agreement with the observed profiles of luminosity, velocity dispersion, rotation, and ellipticity. We obtain a mass-to-light ratio of $M/L = 4.0 \pm 0.4$ and a total cluster mass of $M=(8\pm 1)\times 10^6 M_\odot$. Given that our dynamical model can fit all available observational data very well, there seems to be no need to invoke the presence of an intermediate-mass black hole in the center of G1. △ Less

Submitted 15 April, 2003; v1 submitted 23 January, 2003; originally announced January 2003.

Comments: 12 pages, 3 figures, ApJL accepted, minor changes

Journal ref: Astrophys.J. 589 (2003) L25-L28

Parameters of core-collapse

Authors: Holger Baumgardt, Douglas C. Heggie, Piet Hut, Junichiro Makino

Abstract: This paper considers the phenomenon of deep core collapse in collisional stellar systems, with stars of equal mass. The collapse takes place on some multiple, $ξ^{-1}$, of the central relaxation time, and produces a density profile in which $ρ\propto r^{-α}$, where $α$ is a constant. The parameters $α$ and $ξ$ have usually been determined from simplified models, such as gas and Fokker-Planck mod… ▽ More This paper considers the phenomenon of deep core collapse in collisional stellar systems, with stars of equal mass. The collapse takes place on some multiple, $ξ^{-1}$, of the central relaxation time, and produces a density profile in which $ρ\propto r^{-α}$, where $α$ is a constant. The parameters $α$ and $ξ$ have usually been determined from simplified models, such as gas and Fokker-Planck models, often with the simplification of isotropy. Here we determine the parameters directly from N-body simulations carried out using the newly completed GRAPE-6. △ Less

Submitted 22 January, 2003; v1 submitted 10 January, 2003; originally announced January 2003.

Comments: 5 pages, 3 figures, accepted by MNRAS. Reference added to Table 1

Journal ref: Mon.Not.Roy.Astron.Soc.341:247,2003

Star cluster ecology V: Dissection of an open star cluster--spectroscopy

Authors: Simon Portegies Zwart, Piet Hut, Steve McMillan, Jun Makino

Abstract: We have modeled in detail the evolution of rich open star clusters such as the Pleiades, Praesepe and Hyades, using simulations that include stellar dynamics as well as the effects of stellar evolution. The dynamics is modeled via direct N-body integration, while the evolution of single stars and binaries is followed through the use of fitting formulae and recipes. The feedback of stellar and bi… ▽ More We have modeled in detail the evolution of rich open star clusters such as the Pleiades, Praesepe and Hyades, using simulations that include stellar dynamics as well as the effects of stellar evolution. The dynamics is modeled via direct N-body integration, while the evolution of single stars and binaries is followed through the use of fitting formulae and recipes. The feedback of stellar and binary evolution on the dynamical evolution of the stellar system is taken into account self-consistently. Our model clusters dissolve in the tidal field of the Galaxy in a time span on the order of a billion years. The rate of mass loss is rather constant, about 1Msun per million years. The binary fraction at first is nearly constant in time, then increases slowly near the end of a cluster's lifetime. For clusters which are more than about 100 million years old the fractions of stars in the form of binaries, giants and collision products in the inner few core radii are considerably higher than in the outer regions, beyond the cluster's half mass radius. When stars with masses >2Msun escape from the cluster, they tend to do so with velocities higher than average.... <Abstract abbreviated, read manuscript for rest> △ Less

Submitted 3 January, 2003; originally announced January 2003.

Comments: 17 pages, submitted to MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 351 (2004) 473

On the central structure of M15

Authors: H. Baumgardt, P. Hut, J. Makino, S. McMillan, S. Portegies Zwart

Abstract: We present a detailed comparison between the latest observational data on the kinematical structure of the core of M15, obtained with the Hubble STIS and WFPC2 instruments, and the results of dynamical simulations carried out using the special-purpose GRAPE-6 computer. The observations imply the presence of a significant amount of dark matter in the cluster core. In our dynamical simulations, ne… ▽ More We present a detailed comparison between the latest observational data on the kinematical structure of the core of M15, obtained with the Hubble STIS and WFPC2 instruments, and the results of dynamical simulations carried out using the special-purpose GRAPE-6 computer. The observations imply the presence of a significant amount of dark matter in the cluster core. In our dynamical simulations, neutron stars and/or massive white dwarfs concentrate to the center through mass segregation, resulting in a sharp increase in $M/L$ toward the center. While consistent with the presence of a central black hole, the Hubble data can also be explained by this central concentration of stellar-mass compact objects. The latter interpretation is more conservative, since such remnants result naturally from stellar evolution, although runaway merging leading to the formation of a black hole may also occur for some range of initial conditions. We conclude that no central massive object is required to explain the observational data, although we cannot conclusively exclude such an object at the level of $\sim500-1000$ solar masses. Our findings are unchanged when we reduce the assumed neutron-star retention fraction in our simulations from 100% to 0%. △ Less

Submitted 22 November, 2002; v1 submitted 5 October, 2002; originally announced October 2002.

Comments: 12 pages, 4 figures. Accepted for publication in ApJL

Journal ref: Astrophys.J.582:L21-L24,2003

MODEST-1: Integrating Stellar Evolution and Stellar Dynamics

Authors: Piet Hut, Michael M. Shara, Sverre J. Aarseth, Ralf S. Klessen, James C. Lombardi Jr., Junichiro Makino, Steve McMillan, Onno R. Pols, Peter J. Teuben, Ronald F. Webbink

Abstract: We summarize the main results from MODEST-1, the first workshop on MOdeling DEnse STellar systems. Our goal is to go beyond traditional population synthesis models, by introducing dynamical interactions between single stars, binaries, and multiple systems. The challenge is to define and develop a software framework to enable us to combine in one simulation existing computer codes in stellar evol… ▽ More We summarize the main results from MODEST-1, the first workshop on MOdeling DEnse STellar systems. Our goal is to go beyond traditional population synthesis models, by introducing dynamical interactions between single stars, binaries, and multiple systems. The challenge is to define and develop a software framework to enable us to combine in one simulation existing computer codes in stellar evolution, stellar dynamics, and stellar hydrodynamics. With this objective, the workshop brought together experts in these three fields, as well as other interested astrophysicists and computer scientists. We report here our main conclusions, questions and suggestions for further steps toward integrating stellar evolution and stellar (hydro)dynamics. △ Less

Submitted 31 October, 2002; v1 submitted 15 July, 2002; originally announced July 2002.

Comments: workshop summary, to appear in New Astronomy, 36 pages, 3 figures; this slightly revised version mainly contains more references

Journal ref: New Astron. 8 (2003) 337-370

Long-term evolution of isolated N-body sytems

Authors: Holger Baumgardt, Piet Hut, Douglas C. Heggie

Abstract: We report results of N-body simulations of isolated star clusters, performed up to the point where the clusters are nearly completely dissolved. Our main focus is on the post-collapse evolution of these clusters. We find that after core collapse, isolated clusters evolve along nearly a single sequence of models whose properties are independent of the initial density profile and particle number.… ▽ More We report results of N-body simulations of isolated star clusters, performed up to the point where the clusters are nearly completely dissolved. Our main focus is on the post-collapse evolution of these clusters. We find that after core collapse, isolated clusters evolve along nearly a single sequence of models whose properties are independent of the initial density profile and particle number. Due to the slower expansion of high-N clusters, relaxation times become almost independent of the particle number after several core collapse times, at least for the particle range of our study. As a result, the dissolution times of isolated clusters exhibit a surprisingly weak dependence on N. We find that most stars escape due to encounters between single stars inside the half-mass radius of the cluster. Encounters with binaries take place mostly in the cluster core and account for roughly 15% of all escapers. Encounters between single stars at intermediate radii are also responsible for the build up of a radial anisotropic velocity distribution in the halo. For clusters undergoing core oscillations, escape due to binary stars is efficient only when the cluster center is in a contracted phase. Our simulations show that it takes about 10^5 N-body time units until the global anisotropy reaches its maximum value. The anisotropy increases with particle number and it seems conceivable that isolated star clusters become vulnerable to radial orbit instabilities for large enough N. However, no indication for the onset of such instabilities was seen in our runs. △ Less

Submitted 1 July, 2002; v1 submitted 15 June, 2002; originally announced June 2002.

Comments: 14 pages, 20 figures, MNRAS in press, V2: Order of authors changed in author-field

Journal ref: Mon.Not.Roy.Astron.Soc. 336 (2002) 1069

The Starlab Environment for Dense Stellar Systems

Authors: Piet Hut

Abstract: Traditionally, a simulation of a dense stellar system required choosing an initial model, running an integrator, and analyzing the output. Almost all of the effort went into writing a clever integrator that could handle binaries, triples and encounters between various multiple systems efficiently. Recently, the scope and complexity of these simulations has increased dramatically, for three reaso… ▽ More Traditionally, a simulation of a dense stellar system required choosing an initial model, running an integrator, and analyzing the output. Almost all of the effort went into writing a clever integrator that could handle binaries, triples and encounters between various multiple systems efficiently. Recently, the scope and complexity of these simulations has increased dramatically, for three reasons: 1) the sheer size of the data sets, measured in Terabytes, make traditional `awking and grepping' of a single output file impractical; 2) the addition of stellar evolution data brings qualitatively new challenges to the data reduction; 3) increased realism of the simulations invites realistic forms of `SOS': Simulations of Observations of Simulations, to be compared directly with observations. We are now witnessing a shift toward the construction of archives as well as tailored forms of visualization including the use of virtual reality simulators and planetarium domes, and a coupling of both with budding efforts in constructing virtual observatories. This review describes these new trends, presenting Starlab as the first example of a full software environment for realistic large-scale simulations of dense stellar systems. △ Less

Submitted 25 April, 2002; originally announced April 2002.

Comments: 12 pages, 1 figure, to appear in "Astrophysical Supercomputing Using Particle Simulations", eds. P. Hut and J. Makino, IAU Symposium 208, Tokyo, July 2001

Theory in a Virtual Observatory

Authors: Peter Teuben, Dave DeYoung, Piet Hut, Stuart Levy, Jun Makino, Steve McMillan, Simon Portegies Zwart, Shawn Slavin

Abstract: During the last couple of years, observers have started to make plans for a Virtual Observatory, as a federation of existing data bases, connected through levels of software that enable rapid searches, correlations, and various forms of data mining. We propose to extend the notion of a Virtual Observatory by adding archives of simulations, together with interactive query and visualization capabi… ▽ More During the last couple of years, observers have started to make plans for a Virtual Observatory, as a federation of existing data bases, connected through levels of software that enable rapid searches, correlations, and various forms of data mining. We propose to extend the notion of a Virtual Observatory by adding archives of simulations, together with interactive query and visualization capabilities, as well as ways to simulate observations of simulations in order to compare them with observations. For this purpose, we have already organized two small workshops, earlier in 2001, in Tucson and Aspen. We have also provided concrete examples of theory data, designed to be federated with a Virtual Observatory. These data stem from a project to construct an archive for our large-scale simulations using the GRAPE-6 (a 32-Teraflops special purpose computer for stellar dynamics). We are constructing interfaces by which remote observers can observe these simulations. In addition, these data will enable detailed comparisons between different simulations. △ Less

Submitted 26 November, 2001; originally announced November 2001.

Comments: 4 pages, ADASS-XI conference, http://bima.astro.umd.edu/nemo/nvo/

Orbital Divergence and Relaxation in the Gravitational N-Body Problem

Authors: P. Hut, D. C. Heggie

Abstract: One of the fundamental aspects of statistical behaviour in many-body systems is exponential divergence of neighbouring orbits, which is often discussed in terms of Liapounov exponents. Here we study this topic for the classical gravitational N-body problem. The application we have in mind is to old stellar systems such as globular star clusters, where N~10^6, and so we concentrate on spherical,… ▽ More One of the fundamental aspects of statistical behaviour in many-body systems is exponential divergence of neighbouring orbits, which is often discussed in terms of Liapounov exponents. Here we study this topic for the classical gravitational N-body problem. The application we have in mind is to old stellar systems such as globular star clusters, where N~10^6, and so we concentrate on spherical, centrally concentrated systems with total energy E<0. Hitherto no connection has been made between the time scale for divergence (denoted here by t_e) and the two-body relaxation time scale (t_r), even though both may be calculated by consideration of two-body encounters. In this paper we give a simplified model showing that divergence in phase space is initially roughly exponential, on a timescale proportional to the crossing time (defined as a mean time for a star to cross from one side of the system to another). In this phase t_e is much less than t_r, if N is not too small (i.e. N is much more than about 30). After several e-folding times, the model shows that the divergence slows down. Thereafter the divergence (measured by the energies of the bodies) varies with time as t^{1/2}, on a timescale nearly proportional to the familiar two-body relaxation timescale, i.e. t_e ~ t_r in this phase. These conclusions are illustrated by numerical results. △ Less

Submitted 31 July, 2002; v1 submitted 1 November, 2001; originally announced November 2001.

Comments: 19 pages, 3 figures, accepted by J. Stat. Phys; new figures, numerical results and references; improved text; fig.1 at low resolution in this version

Missing Link Found? --- The ``runaway'' path to supermassive black holes

Authors: Toshikazu Ebisuzaki, Junichiro Makino, Takeshi Go Tsuru, Yoko Funato, Simon Portegies Zwart, Piet Hut, Steve McMillan, Satoki Matsushita, Hironori Matsumoto, Ryohei Kawabe

Abstract: Observations of stellar kinematics, gas dynamics and masers around galactic nuclei have now firmly established that many galaxies host central supermassive black holes (SMBHs) with masses in the range $10^6 \sim 10^9$M$_{\odot}$. However, how these SMBHs formed is not well understood. One reason for this situation is the lack of observations of intermediate-mass BHs (IMBHs), which could bridge t… ▽ More Observations of stellar kinematics, gas dynamics and masers around galactic nuclei have now firmly established that many galaxies host central supermassive black holes (SMBHs) with masses in the range $10^6 \sim 10^9$M$_{\odot}$. However, how these SMBHs formed is not well understood. One reason for this situation is the lack of observations of intermediate-mass BHs (IMBHs), which could bridge the gap between stellar-mass BHs and SMBHs. Recently, this missing link (i.e., an IMBH) has been found in observations made by the ASCA and the Chandra of the central region of the starburst galaxy M82 \citep{MT99, PG99, MT01, Ka01}. Subsequent observations by SUBARU have revealed that this IMBH apparently coincides with a young compact star cluster. Based on these findings, we suggest a new formation scenario for SMBHs. In this scenario, IMBHs first form in young compact star clusters through runaway merging of massive stars. While these IMBHs are forming, the host star clusters sink toward the galactic nucleus through dynamical friction, and upon evaporation deposit their IMBHs near the galactic center. The IMBHs then form binaries and eventually merge via gravitational radiation, forming an SMBH. △ Less

Submitted 14 June, 2001; originally announced June 2001.

Comments: 12 pages 2 figuresm submitted to Astrophysical Journal, Letters

Panel Discussion on Observing Simulations and Simulating Observations

Authors: Piet Hut, Adrienne Cool, Charles Bailyn, Steve McMillan, Mario Livio, Mike Shara

Abstract: N-body simulations of star cluster evolution have reached a high degree of realism, by incorporating more and more elements of stellar dynamics, stellar evolution, and hydrodynamics. At the end of this conference, six participants discussed how to present the increasingly realistic data from star cluster simulations in a way that is most useful for a direct comparison with observations. N-body simulations of star cluster evolution have reached a high degree of realism, by incorporating more and more elements of stellar dynamics, stellar evolution, and hydrodynamics. At the end of this conference, six participants discussed how to present the increasingly realistic data from star cluster simulations in a way that is most useful for a direct comparison with observations. △ Less

Submitted 11 May, 2001; originally announced May 2001.

Comments: 7 pages, LaTeX, to appear in the proceedings for the conference on Stellar Collisions, Mergers, and their Consequences, held in New York, June 2000, ed.: M. Shara [ASP Conference Series]