-
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.
-
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.
-
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.
-
The lives and deaths of star clusters near the Galactic center
Authors:
Simon Portegies Zwart,
Junichiro Makino,
Stephen McMillan,
Piet Hut
Abstract:
We study the evolution and observability of young, compact star clusters near the Galactic center, such as the Arches and Quintuplet systems. The clusters are modeled by integrating the equations of motion of all stars while accounting for the internal evolution of stars and binaries, as well as the effect of the Galactic tidal field. We find that clusters within 150 pc of the Galactic center di…
▽ More
We study the evolution and observability of young, compact star clusters near the Galactic center, such as the Arches and Quintuplet systems. The clusters are modeled by integrating the equations of motion of all stars while accounting for the internal evolution of stars and binaries, as well as the effect of the Galactic tidal field. We find that clusters within 150 pc of the Galactic center dissolve within ~55 Myr, but their projected densities drop below the background density in the direction of the Galactic center within only a few Myr, effectively making these clusters undetectable after that time. Detailed observations of the Arches cluster, taken at face value, suggest that its mass function is unusually flat and that the cluster contains an overabundance of stars more massive than 20 Msun. Our dynamical analysis, however, shows that the observed characteristics of the Arches cluster are consistent with a perfectly normal initial mass function. The observed anomalies are then caused by a combination of observational selection effects and the dynamical evolution of the cluster. We calibrate the current parameters of the Arches cluster using a normal initial mass function and conclude that the cluster is more massive than 40000 Msun, has a half mass radius of about 0.35 pc and is located between 50 and 90 pc from the Galactic center.
△ Less
Submitted 30 August, 2001; v1 submitted 14 February, 2001;
originally announced February 2001.
-
Immersive 4D Interactive Visualization of Large-Scale Simulations
Authors:
Peter Teuben,
Piet Hut,
Stuart Levy,
Jun Makino,
Steve McMillan,
Simon Portegies Zwart,
Mike Shara,
Carter Emmart
Abstract:
In dense clusters a bewildering variety of interactions between stars can be observed, ranging from simple encounters to collisions and other mass-transfer encounters. With faster and special-purpose computers like GRAPE, the amount of data per simulation is now exceeding 1TB. Visualization of such data has now become a complex 4D data-mining problem, combining space and time, and finding intere…
▽ More
In dense clusters a bewildering variety of interactions between stars can be observed, ranging from simple encounters to collisions and other mass-transfer encounters. With faster and special-purpose computers like GRAPE, the amount of data per simulation is now exceeding 1TB. Visualization of such data has now become a complex 4D data-mining problem, combining space and time, and finding interesting events in these large datasets. We have recently starting using the virtual reality simulator, installed in the Hayden Planetarium in the American Museum for Natural History, to tackle some of these problem. This work (http://www.astro.umd.edu/nemo/amnh/) reports on our first ``observations'', modifications needed for our specific experiments, and perhaps field ideas for other fields in science which can benefit from such immersion. We also discuss how our normal analysis programs can be interfaced with this kind of visualization.
△ Less
Submitted 18 January, 2001;
originally announced January 2001.
-
Runaway Collisions in Star Clusters
Authors:
Simon F. Portegies Zwart,
Junichiro Makino,
Stephen L. W. McMillan,
Piet Hut
Abstract:
We study the occurrence of physical collisions between stars in young and compact star cluster. The calculations are performed on the GRAPE-4 with the starlab software environment which include the dynamical evolution and the nuclear evolution of all stars and binaries. The selection of the initial conditions is based on existing and well observed star clusters, such as R136 in the 30 Doradus re…
▽ More
We study the occurrence of physical collisions between stars in young and compact star cluster. The calculations are performed on the GRAPE-4 with the starlab software environment which include the dynamical evolution and the nuclear evolution of all stars and binaries. The selection of the initial conditions is based on existing and well observed star clusters, such as R136 in the 30 Doradus region in the Large Magellanic Cloud and the Arches and Quintuplet star clusters in the vicinity of the Galactic center.
Collisions between stars occurred rather frequently in our models. At any time a single star dominates the collision history of the system. The collision rate of this runaway merger scales with the initial relaxation time of the cluster and is independent on other cluster parameters, such as the initial mass function or the initial density profile of the cluster. Subsequent encounters result in a steady grow in mass of the coagulating star, until it escapes or explodes in a supernova. The collision rate in these models is about 0.00022 collisions per star per Myr for a cluster with an initial relaxation time of 1 Myr.
△ Less
Submitted 11 December, 2000;
originally announced December 2000.
-
How many young star clusters exist in the Galactic center?
Authors:
Simon Portegies Zwart,
Junichiro Makino,
Stephen McMillan,
Piet Hut
Abstract:
We study the evolution and observability of young compact star clusters within about 200pc of the Galactic center. Calculations are performed using direct N-body integration on the GRAPE-4, including the effects of both stellar and binary evolution and the external influence of the Galaxy. The results of these detailed calculations are used to calibrate a simplified model applicable over a wider…
▽ More
We study the evolution and observability of young compact star clusters within about 200pc of the Galactic center. Calculations are performed using direct N-body integration on the GRAPE-4, including the effects of both stellar and binary evolution and the external influence of the Galaxy. The results of these detailed calculations are used to calibrate a simplified model applicable over a wider range of cluster initial conditions. We find that clusters within 200 pc from the Galactic center dissolve within about 70 Myr. However, their projected densities drop below the background density in the direction of the Galactic center within 20 Myr, effectively making these clusters undetectable after that time. Clusters farther from the Galactic center but at the same projected distance are more strongly affected by this selection effect, and may go undetected for their entire lifetimes. Based on these findings, we conclude that the region within 200 pc of the Galactic center could easily harbor some 50 clusters with properties similar to those of the Arches or the Quintuplet systems.
△ Less
Submitted 24 October, 2000; v1 submitted 30 August, 2000;
originally announced August 2000.
-
Star cluster ecology IVa: Dissection of an open star cluster---photometry
Authors:
Simon F. Portegies Zwart,
Stephen L. W. McMillan,
Piet Hut,
Junichiro Makino
Abstract:
The evolution of star clusters is studied using N-body simulations in which the evolution of single stars and binaries are taken self-consistently into account. Initial conditions are chosen to represent relatively young Galactic open clusters, such as the Pleiades, Praesepe and the Hyades. The calculations include a realistic mass function, primordial binaries and the external potential of the…
▽ More
The evolution of star clusters is studied using N-body simulations in which the evolution of single stars and binaries are taken self-consistently into account. Initial conditions are chosen to represent relatively young Galactic open clusters, such as the Pleiades, Praesepe and the Hyades. The calculations include a realistic mass function, primordial binaries and the external potential of the parent Galaxy. Our model clusters are generally significantly flattened in the Galactic tidal field, and dissolve before deep core collapse occurs. The binary fraction decreases initially due to the destruction of soft binaries, but increases later because lower mass single stars escape more easily than the more massive binaries. At late times, the cluster core is quite rich in giants and white dwarfs. There is no evidence for preferential evaporation of old white dwarfs, on the contrary the formed white dwarfs are likely to remain in the cluster. Stars tend to escape from the cluster through the first and second Lagrange points, in the direction of and away from the Galactic center. Mass segregation manifests itself in our models well within an initial relaxation time. As expected, giants and white dwarfs are much more strongly affected by mass segregation than main-sequence stars. Open clusters are dynamically rather inactive. However, the combined effect of stellar mass loss and evaporation of stars from the cluster potential drives its dissolution on a much shorter timescale than if these effects are neglected. The often-used argument that a star cluster is barely older than its relaxation time and therefore cannot be dynamically evolved is clearly in error for the majority of star clusters.
△ Less
Submitted 11 May, 2000;
originally announced May 2000.
-
Star cluster ecology III: Runaway collisions in young compact star clusters
Authors:
Simon Portegies Zwart,
Junichiro Makino,
Stephen L. W. McMillan,
Piet Hut
Abstract:
The evolution of young compact star clusters is studied using N-body simulations in which both stellar evolution and physical collisions between stars are taken into account. The initial conditions are chosen to represent R136, a compact star cluster in the 30 Doradus region of the Large Magellanic Cloud. The present runs do not include the effects of primordial binaries. We find that physical c…
▽ More
The evolution of young compact star clusters is studied using N-body simulations in which both stellar evolution and physical collisions between stars are taken into account. The initial conditions are chosen to represent R136, a compact star cluster in the 30 Doradus region of the Large Magellanic Cloud. The present runs do not include the effects of primordial binaries. We find that physical collisions between stars in these models are frequent, and that the evolution of the most massive stars and the dynamical evolution of the cluster are closely coupled. In all cases, a single star grows steadily in mass through mergers with other stars, forming a very massive (>100 Msun) star in less than 3-4 Myr. The growth rate of this runaway merger is much larger than estimates based on simple cross-section arguments, mainly because the star is typically found in the core and tends to form binaries with other massive stars there. The runaway is ``rejuvenated'' by each new collision, and its lifetime is extended considerably as a consequence. Observationally, such a star will appear in the Hertzsprung-Russell diagram as a blue straggler. When the runaway forms a black hole, the binary in which it is found is usually dissociated.
△ Less
Submitted 26 February, 1999; v1 submitted 1 December, 1998;
originally announced December 1998.
-
High-Performance Special-Purpose Computers in Science
Authors:
Toshiyuki Fukushige,
Piet Hut,
Jun Makino
Abstract:
The next decade will be an exciting time for computational physicists. After 50 years of being forced to use standardized commercial equipment, it will finally become relatively straightforward to adapt one's computing tools to one's own needs. The breakthrough that opens this new era is the now wide-spread availability of programmable chips that allow virtually every computational scientist to…
▽ More
The next decade will be an exciting time for computational physicists. After 50 years of being forced to use standardized commercial equipment, it will finally become relatively straightforward to adapt one's computing tools to one's own needs. The breakthrough that opens this new era is the now wide-spread availability of programmable chips that allow virtually every computational scientist to design his or her own special-purpose computer.
△ Less
Submitted 26 November, 1998;
originally announced November 1998.
-
Astrophysics on the GRAPE Family of Special Purpose Computers
Authors:
Piet Hut,
Jun Makino
Abstract:
The GRAPE-4, the world's fastest computer in 1995-1997, has produced some major scientific results, through a wide diversity of large-scale simulations in astrophysics. Applications have ranged from planetary formation, through the evolution of star clusters and galactic nuclei, to the formation of galaxies and clusters of galaxies.
The GRAPE-4, the world's fastest computer in 1995-1997, has produced some major scientific results, through a wide diversity of large-scale simulations in astrophysics. Applications have ranged from planetary formation, through the evolution of star clusters and galactic nuclei, to the formation of galaxies and clusters of galaxies.
△ Less
Submitted 26 November, 1998;
originally announced November 1998.
-
Stellar Dynamics of Dense Stellar Systems
Authors:
Piet Hut
Abstract:
Stellar dynamics is almost unreasonably well suited for an implementation in terms of special-purpose hardware. Unlike the case of molecular dynamics, stellar dynamics deals exclusively with a long-range force, gravity, which leads to a computational cost scaling as the square of the number of stars involved. While special tricks can lead to a reduction of this cost from $\sim N^2$ to…
▽ More
Stellar dynamics is almost unreasonably well suited for an implementation in terms of special-purpose hardware. Unlike the case of molecular dynamics, stellar dynamics deals exclusively with a long-range force, gravity, which leads to a computational cost scaling as the square of the number of stars involved. While special tricks can lead to a reduction of this cost from $\sim N^2$ to $\sim N\log N$ in the case of very large particle numbers, such tricks are not suitable for all areas within stellar dynamics. When a stellar system is close to equilibrium, and has a very high density, it still pays to compute all interactions on a star by star basis, even for $N=10^5$. Any $cN\log N$ approach would either gloss over the subtle net effects of near-canceling interactions, driving the evolution of such a system, or would carry a prohibitively large coefficient $c$. This paper presents a brief introduction to the stellar dynamics of dense stellar systems, aimed at researchers using special purpose computers in other branches of physics.
△ Less
Submitted 8 April, 1998;
originally announced April 1998.
-
On the dissolution of evolving star clusters
Authors:
Simon Portegies Zwart,
Piet Hut,
Junichiro Makino,
Stephen McMillan
Abstract:
Using direct N-body simulations which include both the evolution of single stars and the tidal field of the parent galaxy, we study the dynamical evolution of globular clusters and rich open clusters. We compare our results with other N-body simulations and Fokker-Planck calculations. Our simulations, performed on the GRAPE-4, employ up to 32,768 stars. The results are not in agreement with Fokk…
▽ More
Using direct N-body simulations which include both the evolution of single stars and the tidal field of the parent galaxy, we study the dynamical evolution of globular clusters and rich open clusters. We compare our results with other N-body simulations and Fokker-Planck calculations. Our simulations, performed on the GRAPE-4, employ up to 32,768 stars. The results are not in agreement with Fokker-Planck models, in the sense that the lifetimes of stellar systems derived using the latter are an order of magnitude smaller than those obtained in our simulations. For our standard run, Fokker-Plank calculations obtained a lifetime of 0.28 Gyr, while our equivalent N-body calculations find about 4 Gyr. The principal reason for the discrepancy is that a basic assumption of the Fokker-Plank approach is not valid for typical cluster parameters. The stellar evolution timescale is comparable to the dynamical timescale, and therefore the assumption of dynamical equilibrium leads to an overestimate of the dynamical effects of mass loss. Our results suggest that the region in parameter space for which Fokker-Planck studies of globular cluster evolution, including the effects of both stellar evolution and the galactic tidal field, are valid is limited. The discrepancy is largest for clusters with short lifetimes.
△ Less
Submitted 9 May, 1998; v1 submitted 9 March, 1998;
originally announced March 1998.
-
The Equilibrium-Tide Model for Tidal Friction
Authors:
Peter P. Eggleton,
Ludmila G. Kiseleva,
Piet Hut
Abstract:
We derive from first principles equations governing (a) the quadrupole tensor of a star distorted by both rotation and the presence of a companion in a possibly eccentric orbit, (b) a functional form for the dissipative force of tidal friction, based on the concept that the rate of energy loss from a time-dependent tide should be a positive-definite function of the rate of change of the quadrupo…
▽ More
We derive from first principles equations governing (a) the quadrupole tensor of a star distorted by both rotation and the presence of a companion in a possibly eccentric orbit, (b) a functional form for the dissipative force of tidal friction, based on the concept that the rate of energy loss from a time-dependent tide should be a positive-definite function of the rate of change of the quadrupole tensor as seen in the frame which rotates with the star, and (c) the equations governing the rates of change of the magnitude and direction of the stellar rotation, and the orbital period and eccentricity, based on the concept of the Laplace-Runge-Lenz vector. Our analysis leads relatively simply to a closed set of equations, valid for arbitrary inclination of the stellar spin to the orbit. The results are equivalent to classical results based on the rather less clear principle that the tidal bulge lags behind the line of centres by some time determined by the rate of dissipation; our analysis gives the effective lag time as a function of the dissipation rate and the quadrupole moment.
We discuss briefly some possible applications of the formulation.
△ Less
Submitted 26 January, 1998;
originally announced January 1998.
-
HOP: A New Group-Finding Algorithm for N-body Simulations
Authors:
Daniel J. Eisenstein,
Piet Hut
Abstract:
We describe a new method (HOP) for identifying groups of particles in N-body simulations. Having assigned to every particle an estimate of its local density, we associate each particle with the densest of the N_hop particles nearest to it. Repeating this process allows us to trace a path, within the particle set itself, from each particle in the direction of increasing density. The path ends whe…
▽ More
We describe a new method (HOP) for identifying groups of particles in N-body simulations. Having assigned to every particle an estimate of its local density, we associate each particle with the densest of the N_hop particles nearest to it. Repeating this process allows us to trace a path, within the particle set itself, from each particle in the direction of increasing density. The path ends when it reaches a particle that is its own densest neighbor; all particles reaching the same such particle are identified as a group. Combined with an adaptive smoothing kernel for finding the densities, this method is spatially adaptive, coordinate-free, and numerically straight-forward. One can proceed to process the output by truncating groups at a particular density contour and combining groups that share a (possibly different) density contour. While the resulting algorithm has several user-chosen parameters, we show that the results are insensitive to most of these, the exception being the outer density cutoff of the groups.
△ Less
Submitted 15 December, 1997;
originally announced December 1997.
-
Smooth Particle Hydrodynamics: Models, Applications, and Enabling Technologies
Authors:
Piet Hut,
Lars Hernquist,
George Lake,
Jun Makino,
Steve McMillan,
Thomas Sterling
Abstract:
We present the results from a two-day study in which we discussed various implementations of Smooth Particle Hydrodynamics (SPH), one of the leading methods used across a variety of areas of large-scale astrophysical simulations. In particular, we evaluated the suitability of designing special hardware extensions, to further boost the performance of the high-end general purpose computers current…
▽ More
We present the results from a two-day study in which we discussed various implementations of Smooth Particle Hydrodynamics (SPH), one of the leading methods used across a variety of areas of large-scale astrophysical simulations. In particular, we evaluated the suitability of designing special hardware extensions, to further boost the performance of the high-end general purpose computers currently used for those simulations. We considered a range of hybrid architectures, consisting of a mix of custom LSI and reconfigurable logic, combining the extremely high throughput of Special-Purpose Devices (SPDs) with the flexibility of reconfigurable structures, based on Field Programmable Gate Arrays (FPGAs). The main findings of our workshop consist of a clarification of the decomposition of the computational requirements, together with specific estimates for cost/performance improvements that can be obtained at each stage in this decomposition, by using enabling hardware technology to accelerate the performance of general purpose computers.
△ Less
Submitted 20 October, 1997;
originally announced October 1997.
-
Aspects of SU(N_c) Gauge Theories in the Limit of Small Number of Colors
Authors:
Stanley J. Brodsky,
Patrick Huet
Abstract:
We investigate properties of the color space of $SU(N_c)$ gauge theories in the limit of small number of colors $(N_c \to 0)$ and large number of flavors. More generally, we introduce a rescaling of $α_s$ and $n_f$ which assigns a finite limit to colored quantities as $N_c \to 0$, which reproduces their known large-$N_c$ limit, and which expresses them as an analytic function of $N_c^2$ for arbi…
▽ More
We investigate properties of the color space of $SU(N_c)$ gauge theories in the limit of small number of colors $(N_c \to 0)$ and large number of flavors. More generally, we introduce a rescaling of $α_s$ and $n_f$ which assigns a finite limit to colored quantities as $N_c \to 0$, which reproduces their known large-$N_c$ limit, and which expresses them as an analytic function of $N_c^2$ for arbitrary value of $N_c$. The vanishing-$N_c$ limit has an Abelian character and is also the small-$N_c$ limit of $[U(1)]^{N_c-1}$. This limit does not have an obvious quantum field theory interpretation; however, it provides practical consistency checks on QCD perturbative quantities by comparing them to their QED counterparts. Our analysis also describes the two-dimensional topological structure involved in the interpretation of the small $N_c$-limit in color space.
△ Less
Submitted 11 August, 1997; v1 submitted 31 July, 1997;
originally announced July 1997.
-
Star Cluster Ecology II: Binary evolution with single-star encounters
Authors:
Simon F. Portegies Zwart,
Piet Hut,
Stephen L. W. McMillan,
Frank Verbunt
Abstract:
Three-body effects greatly complicate stellar evolution. We model the effects of encounters of binaries with single stars, based on parameters chosen from conditions prevalent in the cores of globular clusters. For our three-body encounters, we start with a population of primordial binaries, and choose incoming stars from the evolving single star populations of Portgies Zwart, Hut and Verbunt (1…
▽ More
Three-body effects greatly complicate stellar evolution. We model the effects of encounters of binaries with single stars, based on parameters chosen from conditions prevalent in the cores of globular clusters. For our three-body encounters, we start with a population of primordial binaries, and choose incoming stars from the evolving single star populations of Portgies Zwart, Hut and Verbunt (1997). In addition, we study the formation of new binaries through tidal capture among the single stars. In subsequent papers in this series, we will combine stellar evolution with the dynamics of a full N-body system.
△ Less
Submitted 10 June, 1997;
originally announced June 1997.
-
Gravitational Thermodynamics
Authors:
Piet Hut
Abstract:
The gravitational N-body problem, for $N>2$, is the oldest unsolved problem in mathematical physics. Some of the most ideal examples that can be found in nature are globular star clusters, with $N \sim 10^6$. In this overview, I discuss six types of fundamental sources of unpredictability, each of which poses a different challenge to attempts to determine the long-term behavior of these systems,…
▽ More
The gravitational N-body problem, for $N>2$, is the oldest unsolved problem in mathematical physics. Some of the most ideal examples that can be found in nature are globular star clusters, with $N \sim 10^6$. In this overview, I discuss six types of fundamental sources of unpredictability, each of which poses a different challenge to attempts to determine the long-term behavior of these systems, governed by a peculiar type of thermodynamics.
△ Less
Submitted 29 April, 1997;
originally announced April 1997.
-
The GRAPE-4, a Teraflops Stellar Dynamics Computer
Authors:
Piet Hut
Abstract:
Recently, special-purpose computers have surpassed general-purpose computers in the speed with which large-scale stellar dynamics simulations can be performed. Speeds up to a Teraflops are now available, for simulations in a variety of fields, such as planetary formation, star cluster dynamics, galactic nuclei, galaxy interactions, galaxy formation, large scale structure, and gravitational lensi…
▽ More
Recently, special-purpose computers have surpassed general-purpose computers in the speed with which large-scale stellar dynamics simulations can be performed. Speeds up to a Teraflops are now available, for simulations in a variety of fields, such as planetary formation, star cluster dynamics, galactic nuclei, galaxy interactions, galaxy formation, large scale structure, and gravitational lensing. Future speed increases for special-purpose computers will be even more dramatic: a Petaflops version, tentatively named the GRAPE-6, could be built within a few years, whereas general-purpose computers are expected to reach this speed somewhere in the 2010-2015 time frame. Boards with a handful of chips from such a machine could be made available to individual astronomers. Such a board, attached to a fast workstation, will then deliver Teraflops speeds on a desktop, around the year 2000.
△ Less
Submitted 28 April, 1997;
originally announced April 1997.
-
Time Symmetrization Meta-Algorithms
Authors:
Piet Hut,
Yoko Funato,
Eiichiro Kokubo,
Junichiro Makino,
Steve McMillan
Abstract:
We present two types of meta-algorithm that can greatly improve the accuracy of existing algorithms for integrating the equations of motion of dynamical systems. The first meta-algorithm takes an integrator that is time-symmetric only for constant time steps, and ensures time-symmetry even in the case of varying time steps. The second meta-algorithm can be applied to any self-starting integratio…
▽ More
We present two types of meta-algorithm that can greatly improve the accuracy of existing algorithms for integrating the equations of motion of dynamical systems. The first meta-algorithm takes an integrator that is time-symmetric only for constant time steps, and ensures time-symmetry even in the case of varying time steps. The second meta-algorithm can be applied to any self-starting integration scheme to create time-symmetry, both for constant and for variable time steps, even if the original scheme was not time-symmetric. Our meta-algorithms are most effective for Hamilton systems or systems with periodic solutions. If the system is not Hamiltonian (for example, if some dissipative force exists), our methods are still useful so long as the dissipation is small.
△ Less
Submitted 28 April, 1997;
originally announced April 1997.
-
GRAPE-6: A Petaflops Prototype
Authors:
Piet Hut,
Jeffrey M. Arnold,
Junichiro Makino,
Stephen L. W. McMillan,
Thomas L. Sterling
Abstract:
We present the outline of a research project aimed at designing and constructing a hybrid computing system that can be easily scaled up to petaflops speeds. As a first step, we envision building a prototype which will consist of three main components: a general-purpose, programmable front end, a special-purpose, fully hardwired computing engine, and a multi-purpose, reconfigurable system. The dr…
▽ More
We present the outline of a research project aimed at designing and constructing a hybrid computing system that can be easily scaled up to petaflops speeds. As a first step, we envision building a prototype which will consist of three main components: a general-purpose, programmable front end, a special-purpose, fully hardwired computing engine, and a multi-purpose, reconfigurable system. The driving application will be a suite of particle-based large-scale simulations in various areas of physics. The prototype system will achieve performance in the $\sim 50 - 100$ teraflops range for a broad class of applications in this area. The combination of a hardwired petaflops-class computational engine and a front end with sustained speed on the order of 10 gigaflops can produce extremely high performance, but only for the limited class of problems in which there exists a single bottleneck with computing cost dominating the total. While the calculation for which the Grape-4 (our system's immediate predecessor) was designed is a prime example of such a problem, in many other applications the primary computational bottleneck, while still related to an inverse-square (gravitational, Coulomb, etc.) force, requires less than 99% of the computing power. Although the remainder of the CPU time is typically dominated by just one secondary bottleneck, its nature varies greatly from problem to problem. It is not cost-effective to attempt to design custom chips for each new problem that arises. FPGA-based systems can restore the balance, guaranteeing scalability from the teraflops to the petaflops domain, while still retaining significant flexibility. (abbreviated abstract)
△ Less
Submitted 17 April, 1997;
originally announced April 1997.
-
Star Cluster Ecology I: A Cluster Core with Encounters between Single Stars
Authors:
Simon F. Portegies Zwart,
Piet Hut,
Frank Verbunt
Abstract:
Star clusters with a high central density contain an ecological network of evolving binaries, affected by interactions with passing stars, while in turn affecting the energy budget of the cluster as a whole by giving off binding energy. This is the first paper in a series aimed at providing the tools for increasingly realistic simulations of these ecological networks.
Star clusters with a high central density contain an ecological network of evolving binaries, affected by interactions with passing stars, while in turn affecting the energy budget of the cluster as a whole by giving off binding energy. This is the first paper in a series aimed at providing the tools for increasingly realistic simulations of these ecological networks.
△ Less
Submitted 9 January, 1997;
originally announced January 1997.
-
Long Range Physics in a Hot Non-Abelian Plasma
Authors:
Patrick Huet,
D. T. Son
Abstract:
We derive a set of equations describing the real time dynamics of modes with spatial momentum of order g^2 T in a high temperature gauge theory, where g is the coupling constant and T is the temperature. This dynamics is stochastic in nature. Important implications for baryon number violation at high temperature and for the physics at the electroweak phase transition, are discussed.
We derive a set of equations describing the real time dynamics of modes with spatial momentum of order g^2 T in a high temperature gauge theory, where g is the coupling constant and T is the temperature. This dynamics is stochastic in nature. Important implications for baryon number violation at high temperature and for the physics at the electroweak phase transition, are discussed.
△ Less
Submitted 2 January, 1997; v1 submitted 5 October, 1996;
originally announced October 1996.
-
Merger Rate of Equal-Mass Spherical Galaxies
Authors:
Junichiro Makino,
Piet Hut
Abstract:
We present cross sections and reaction rates for merging to occur during encounters of equal-mass spherical galaxies. As an application, we determine the rate of galaxy merging in clusters of galaxies. We present results for two types of Plummer models (a full and a truncated one), two King models and the Hernquist model. Cross sections are determined on the basis of a large number ($\sim 500$)…
▽ More
We present cross sections and reaction rates for merging to occur during encounters of equal-mass spherical galaxies. As an application, we determine the rate of galaxy merging in clusters of galaxies. We present results for two types of Plummer models (a full and a truncated one), two King models and the Hernquist model. Cross sections are determined on the basis of a large number ($\sim 500$) of simulations of galaxy encounters, using the 10-Gigaflops GRAPE 3A special-purpose computer. We characterize the overall merger rate of galaxies in a galaxy cluster by a single number, derived from our cross sections by an integration over galaxy encounter velocities in the limit of a constant density in velocity space. For small clusters, where the cluster velocity dispersion may not significantly exceed the internal velocity dispersion of the individual galaxies, this constant-density approximation may not be valid. For those cases, we present separate results, based on integrations of our cross sections over Maxwellian velocity distributions. Finally, tidal effects from the cluster potential as well as from neighboring galaxies may prevent a barely bound galaxy pair from spiraling in after their first encounter. We give a quantitative estimate of the resulting reduction in the actual merger rate, due to these tidal interactions.
△ Less
Submitted 27 March, 1997; v1 submitted 25 August, 1996;
originally announced August 1996.
-
CPT Violation from Planck Scale Physics
Authors:
Patrick Huet
Abstract:
I address the phenomenology of CPT violation in the neutral kaon system under the assumption that it originates from Planck scale physics. This assumption opens the door to a new set of CPT violating parameters whose phenomenology is distinct from the $Δ$ parameter usually considered in the Hamiltonian. The origin of these parameters reflects a possible departure from a $S$-matrix evolution. Exi…
▽ More
I address the phenomenology of CPT violation in the neutral kaon system under the assumption that it originates from Planck scale physics. This assumption opens the door to a new set of CPT violating parameters whose phenomenology is distinct from the $Δ$ parameter usually considered in the Hamiltonian. The origin of these parameters reflects a possible departure from a $S$-matrix evolution. Existing bounds on CPT violation are near the expected range based on naive dimensional analysis. This provides a strong incentive to pursue the quest of CPT violation in near-future kaon experiments. (Talk given at the Workshop on K Physics, Orsay, France, May 30 - June 4, 1996.)
△ Less
Submitted 29 July, 1996; v1 submitted 25 July, 1996;
originally announced July 1996.
-
The Role of Binaries in the Dynamical Evolution of the Core of a Globular Cluster
Authors:
Piet Hut
Abstract:
The size of the core is one of the main diagnostics of the evolutionary state of a globular cluster. Much has been learned over the last few years about the behavior of the core radius during and after core collapse, under a variety of different conditions related to the presence or absence of large numbers of binaries. An overview is presented of the basic physical principles that can be used t…
▽ More
The size of the core is one of the main diagnostics of the evolutionary state of a globular cluster. Much has been learned over the last few years about the behavior of the core radius during and after core collapse, under a variety of different conditions related to the presence or absence of large numbers of binaries. An overview is presented of the basic physical principles that can be used to estimate the core radius. Four different situations are discussed, and expressions are presented for the ratio $r_c/r_h$ of core radius to half mass radius. The regimes are: deep collapse in the absence of primordial binaries; steady post-collapse evolution after primordial binaries have been burned up; chaotic post-collapse evolution under the same conditions; and post-collapse evolution in the presence of primordial binaries. In addition, modifications to all of these cases are indicated for the more realistic situation where effects of the galactic tidal field are taken into account.
△ Less
Submitted 3 May, 1996;
originally announced May 1996.
-
Time-Symmetrized Kustaanheimo-Stiefel Regularization
Authors:
Y. Funato,
P. Hut,
S. McMillan,
J. Makino
Abstract:
In this paper we describe a new algorithm for the long-term numerical integration of the two-body problem, in which two particles interact under a Newtonian gravitational potential. Although analytical solutions exist in the unperturbed and weakly perturbed cases, numerical integration is necessary in situations where the perturbation is relatively strong. Kustaanheimo--Stiefel (KS) regularizati…
▽ More
In this paper we describe a new algorithm for the long-term numerical integration of the two-body problem, in which two particles interact under a Newtonian gravitational potential. Although analytical solutions exist in the unperturbed and weakly perturbed cases, numerical integration is necessary in situations where the perturbation is relatively strong. Kustaanheimo--Stiefel (KS) regularization is widely used to remove the singularity in the equations of motion, making it possible to integrate orbits having very high eccentricity. However, even with KS regularization, long-term integration is difficult, simply because the required accuracy is usually very high. We present a new time-integration algorithm which has no secular error in either the binding energy or the eccentricity, while allowing variable stepsize. The basic approach is to take a time-symmetric algorithm, then apply an implicit criterion for the stepsize to ensure strict time reversibility. We describe the algorithm in detail and present the results of numerical tests involving long-term integration of binaries and hierarchical triples. In all cases studied, we found no systematic error in either the energy or the angular momentum. We also found that its calculation cost does not become higher than those of existing algorithms. By contrast, the stabilization technique, which has been widely used in the field of collisional stellar dynamics, conserves energy very well but does not conserve angular momentum.
△ Less
Submitted 4 April, 1996;
originally announced April 1996.
-
Binary--single-star scattering VII. Hard Binary Exchange Cross Sections for Arbitrary Mass Ratios: Numerical Results and Semi-Analytic Fits
Authors:
D. C. Heggie,
P. Hut,
S. L. W. McMillan
Abstract:
We present the first comprehensive fitting formula for exchange reactions of arbitrary mass ratios. In a comparison with numerical results, this expression is shown to be accurate in the hard binary limit to within 25\% for most mass ratios. The result will be useful in forming quantitative estimates for the branching ratios of various exchange reactions in astrophysical applications. For exampl…
▽ More
We present the first comprehensive fitting formula for exchange reactions of arbitrary mass ratios. In a comparison with numerical results, this expression is shown to be accurate in the hard binary limit to within 25\% for most mass ratios. The result will be useful in forming quantitative estimates for the branching ratios of various exchange reactions in astrophysical applications. For example, it can be used to construct quantitative formation scenarios for unusual objects in globular clusters, such as binaries containing a pulsar.
△ Less
Submitted 2 April, 1996;
originally announced April 1996.
-
Binary--single-star scattering VI. Automatic Determination of Interaction Cross Sections
Authors:
Stephen L. W. McMillan,
Piet Hut
Abstract:
Scattering encounters between binaries and single stars play a central role in determining the dynamical evolution of a star cluster. In addition, three-body scattering can give rise to many interesting exceptional objects: merging can produce blue stragglers; exchange can produce binaries containing millisecond pulsars in environments quite different from those in which the pulsars were spun up…
▽ More
Scattering encounters between binaries and single stars play a central role in determining the dynamical evolution of a star cluster. In addition, three-body scattering can give rise to many interesting exceptional objects: merging can produce blue stragglers; exchange can produce binaries containing millisecond pulsars in environments quite different from those in which the pulsars were spun up; various types of X-ray binaries can be formed, and their activity can be either shut off or triggered as a result of triple interactions.
To date, all published results on three-body scattering have relied on human guidance for determining the correct parameter range for the envelope within which to perform Monte--Carlo scattering experiments. In this paper, we describe the first fully automatic determination of cross sections and reaction rates for binary--single-star scattering. Rather than relying on human inspection of pilot calculations, we have constructed a feedback system that ensures near-optimal coverage of parameter space while guaranteeing completeness. We illustrate our approach with a particular example, in which we describe the results of a three-body encounter between three main-sequence stars of different masses. We provide total cross sections, as well as branching ratios for the various different types of two-body mergers, three-body mergers, and exchanges, both non-resonant and resonant. The companion paper in this series, Paper VII, provides a full survey of unequal-mass three-body scattering for hard binaries in the point-mass limit.
△ Less
Submitted 2 April, 1996;
originally announced April 1996.
-
Exchange Cross Sections for Hard Binaries
Authors:
Douglas C. Heggie,
Piet Hut,
Steven L. W. McMillan
Abstract:
We present results on the exchange cross section for the interaction between a hard binary and a field of single stars, for arbitrary masses. The results are based partly on extensive numerical scattering experiments, and partly on analytic estimates of the mass-dependence of the cross section. They can be used to estimate the rate of exchange in an arbitrary mixture of masses, provided that the…
▽ More
We present results on the exchange cross section for the interaction between a hard binary and a field of single stars, for arbitrary masses. The results are based partly on extensive numerical scattering experiments, and partly on analytic estimates of the mass-dependence of the cross section. They can be used to estimate the rate of exchange in an arbitrary mixture of masses, provided that the binary is hard.
△ Less
Submitted 15 March, 1996;
originally announced March 1996.
-
Overlap formulation of Majorana--Weyl fermions
Authors:
Patrick Huet,
Rajamani Narayanan,
Herbert Neuberger
Abstract:
An overlap method for regularizing Majorana--Weyl fermions interacting with gauge fields is presented. A mod(2) index is introduced in relation to the anomalous violation of a discrete global chiral symmetry. Most of the paper is restricted to 2 dimensions but generalizations to 2+8k dimensions should be straightforward.
An overlap method for regularizing Majorana--Weyl fermions interacting with gauge fields is presented. A mod(2) index is introduced in relation to the anomalous violation of a discrete global chiral symmetry. Most of the paper is restricted to 2 dimensions but generalizations to 2+8k dimensions should be straightforward.
△ Less
Submitted 29 February, 1996;
originally announced February 1996.
-
The Role of Binaries in the Dynamical Evolution of Globular Clusters
Authors:
Piet Hut
Abstract:
Three important developments are vastly increasing our understanding of the role of binaries in the dynamical evolution of globular clusters. From the observational side, the Hubble Space Telescope has shown us detailed pictures of the densest areas in post-collapse cluster cores. From the computational side, the Grape-4 special-purpose hardware is now allowing us to model small globular cluster…
▽ More
Three important developments are vastly increasing our understanding of the role of binaries in the dynamical evolution of globular clusters. From the observational side, the Hubble Space Telescope has shown us detailed pictures of the densest areas in post-collapse cluster cores. From the computational side, the Grape-4 special-purpose hardware is now allowing us to model small globular clusters on a star-by-star basis, and has already given us the first direct evidence of the occurrence of gravothermal oscillations in such systems. From the theoretical astrophysics side, integrated simulations are now becoming feasible that combine stellar dynamics with stellar evolution and hydrodynamics. Given these three developments, we can expect the current rapid progress in our understanding of globular cluster evolution to continue at an even higher rate during the foreseeable future. In this review an outline is given of the current status of globular cluster simulations, and the expected progress over the next five years.
△ Less
Submitted 29 February, 1996;
originally announced February 1996.
-
Dark Matter in Globular Clusters
Authors:
Douglas C. Heggie,
Piet Hut
Abstract:
We first review reasons why dark matter is an interesting issue in connection with star clusters. Next we consider to what extent the presence of dark matter is consistent with their dynamics and structure. We review various model-dependent and model-independent methods which have been applied to two well studied clusters, NGC 6397 and 47 Tuc. We suggest that about half of the mass in each objec…
▽ More
We first review reasons why dark matter is an interesting issue in connection with star clusters. Next we consider to what extent the presence of dark matter is consistent with their dynamics and structure. We review various model-dependent and model-independent methods which have been applied to two well studied clusters, NGC 6397 and 47 Tuc. We suggest that about half of the mass in each object is still unobserved, possibly in the form of a mixture of low-mass stars and white dwarfs.
△ Less
Submitted 24 November, 1995;
originally announced November 1995.
-
Electroweak Baryogenesis in Supersymmetric Models
Authors:
Patrick Huet,
Ann E. Nelson
Abstract:
The baryon density which may be produced during the electroweak phase transition in supersymmetric models is computed, taking into account the previously neglected effects of transport, strong and weak anomalous fermion number violation, thermal scattering, and a new method for computing \cp\ violating processes during the transition. We can account for the observed baryon asymmetry, provided ne…
▽ More
The baryon density which may be produced during the electroweak phase transition in supersymmetric models is computed, taking into account the previously neglected effects of transport, strong and weak anomalous fermion number violation, thermal scattering, and a new method for computing \cp\ violating processes during the transition. We can account for the observed baryon asymmetry, provided new \cp-violating phases are greater than $\sim 10^{-(2-4)}$, and some superpartners are light enough to be relevant during the transition, which takes place at a temperature of (50-100) GeV. In one case, light superpartners are the top squarks and the charginos and/or the neutralinos; in another case the top squarks and both Higgs doublets are light. Our calculation is easily extended to the case of a general two Higgs model, where we find sufficient baryogenesis provided that a certain combination of parameters in the Higgs potential leads to a \cp\ violating space dependent phase in the top quark mass of order $ 10^{-3}$.
△ Less
Submitted 16 January, 1996; v1 submitted 1 July, 1995;
originally announced June 1995.
-
CP Violation and Electroweak Baryogenesis in Extensions of the Standard Model
Authors:
Patrick Huet,
Ann E. Nelson
Abstract:
We develop a new and general method to calculate the effects of \cp\ violation from extensions of the standard model on the mechanism of electroweak baryogenesis. We illustrate its applicability in the framework of two-higgs doublet models.
We develop a new and general method to calculate the effects of \cp\ violation from extensions of the standard model on the mechanism of electroweak baryogenesis. We illustrate its applicability in the framework of two-higgs doublet models.
△ Less
Submitted 20 June, 1995; v1 submitted 28 April, 1995;
originally announced April 1995.
-
Binary-Binary Interactions and the Formation of the PSR B1620-26 Triple System in M4
Authors:
F. A. Rasio,
S. McMillan,
P. Hut
Abstract:
The hierarchical triple system containing the millisecond pulsar PSR B1620-26 in M4 is the first triple star system ever detected in a globular cluster. Such systems should form in globular clusters as a result of dynamical interactions between binaries. We propose that the triple system containing PSR B1620-26 formed through an exchange interaction between a wide primordial binary and a {\it pr…
▽ More
The hierarchical triple system containing the millisecond pulsar PSR B1620-26 in M4 is the first triple star system ever detected in a globular cluster. Such systems should form in globular clusters as a result of dynamical interactions between binaries. We propose that the triple system containing PSR B1620-26 formed through an exchange interaction between a wide primordial binary and a {\it pre-existing\/} binary millisecond pulsar. This scenario would have the advantage of reconciling the $\sim10^9\,$yr timing age of the pulsar with the much shorter lifetime of the triple system in the core of M4.
△ Less
Submitted 22 September, 1994;
originally announced September 1994.
-
Electroweak Baryogenesis and the Standard Model
Authors:
Patrick Huet
Abstract:
Electroweak baryogenesis is addressed within the context of the standard model of particle physics. Although the minimal standard model has the means of fulfilling the three Sakharov's conditions, it falls short to explaining the making of the baryon asymmetry of the universe. In particular, it is demonstrated that the phase of the CKM mixing matrix is an insufficient source of {\it CP} violatio…
▽ More
Electroweak baryogenesis is addressed within the context of the standard model of particle physics. Although the minimal standard model has the means of fulfilling the three Sakharov's conditions, it falls short to explaining the making of the baryon asymmetry of the universe. In particular, it is demonstrated that the phase of the CKM mixing matrix is an insufficient source of {\it CP} violation. The shortcomings of the standard model could be bypassed by enlarging the symmetry breaking sector and adding a new source of {\it CP} violation.
△ Less
Submitted 14 June, 1994; v1 submitted 13 June, 1994;
originally announced June 1994.
-
Testing CPT and Quantum Mechanics in the K_0-Kbar_0 System
Authors:
Patrick Huet
Abstract:
I report a recent study done in collaboration with M.E. Peskin, on the time dependence of a kaon beam propagating according to a generalization of quantum mechanics due to Ellis, Hagelin, Nanopoulos and Srednicki, in which CP- and CPT-violating signatures arise from the evolution of pure states to mixed states. Constraints on the magnitude of its parameters are established on the basis of existi…
▽ More
I report a recent study done in collaboration with M.E. Peskin, on the time dependence of a kaon beam propagating according to a generalization of quantum mechanics due to Ellis, Hagelin, Nanopoulos and Srednicki, in which CP- and CPT-violating signatures arise from the evolution of pure states to mixed states. Constraints on the magnitude of its parameters are established on the basis of existing experimental data. New facilities such as phi-factories are shown to be particularly adequate to study this generalization from \qm\ and to disentangle its parameters from other CPT violating perturbations of the kaon system. ( To appear in Proceedings of the First International Conference on Phenomenology of Unification from Present to Future, Roma, ITALY -- March 23-26, 1994.)
△ Less
Submitted 12 May, 1994;
originally announced May 1994.
-
Electroweak Baryogenesis and Standard Model CP Violation
Authors:
Patrick Huet,
Eric Sather
Abstract:
We analyze the mechanism of electroweak baryogenesis proposed by Farrar and Shaposhnikov in which the phase of the CKM mixing matrix is the only source of $CP$ violation. This mechanism is based on a phase separation of baryons via the scattering of quasiparticles by the wall of an expanding bubble produced at the electroweak phase transition. In agreement with the recent work of Gavela, Hernánd…
▽ More
We analyze the mechanism of electroweak baryogenesis proposed by Farrar and Shaposhnikov in which the phase of the CKM mixing matrix is the only source of $CP$ violation. This mechanism is based on a phase separation of baryons via the scattering of quasiparticles by the wall of an expanding bubble produced at the electroweak phase transition. In agreement with the recent work of Gavela, Hernández, Orloff and Pène, we conclude that QCD damping effects reduce the asymmetry produced to a negligible amount. We interpret the damping as quantum decoherence. We compute the asymmetry analytically. Our analysis reflects the observation that only a thin, outer layer of the bubble contributes to the coherent scattering of the quasiparticles. The generality of our arguments rules out any mechanism of electroweak baryogenesis that does not make use of a new source of $CP$ violation.
△ Less
Submitted 19 April, 1994;
originally announced April 1994.
-
Violation of CPT and Quantum Mechanics in the K0--K0bar System
Authors:
Patrick Huet,
Michael E. Peskin
Abstract:
We reconsider the model of quantum mechanics violation in the $\ko$--$\kob$ system, due to Ellis, Hagelin, Nanopoulos, and Srednicki, in which \cp- and \cpt-violating signatures arise from the evolution of pure states into mixed states. We present a formalism for computing time-dependent asymmetries in this model and show that present data constrains its parameters significantly. In the future,…
▽ More
We reconsider the model of quantum mechanics violation in the $\ko$--$\kob$ system, due to Ellis, Hagelin, Nanopoulos, and Srednicki, in which \cp- and \cpt-violating signatures arise from the evolution of pure states into mixed states. We present a formalism for computing time-dependent asymmetries in this model and show that present data constrains its parameters significantly. In the future, this model will be put to very stringent tests at a $φ$ factory. We present the theory of these tests and show the relation between particular $φ$ decay correlations and the parameters of quantum mechanics violation.
△ Less
Submitted 19 August, 1994; v1 submitted 9 March, 1994;
originally announced March 1994.
-
Blue Stragglers as Tracers of Globular Cluster Evolution
Authors:
Piet Hut
Abstract:
Blue stragglers are natural phenomena in star clusters. They originate through mass transfer in isolated binaries, as well as through encounters between two or more stars, in a complex interplay between stellar dynamics and stellar evolution. While this interplay cannot be modeled quantitatively at present, we will be able to do so in one or two years time. With this prospect, the present paper…
▽ More
Blue stragglers are natural phenomena in star clusters. They originate through mass transfer in isolated binaries, as well as through encounters between two or more stars, in a complex interplay between stellar dynamics and stellar evolution. While this interplay cannot be modeled quantitatively at present, we will be able to do so in one or two years time. With this prospect, the present paper is written largely as a preview.
△ Less
Submitted 8 March, 1993;
originally announced March 1993.
-
Hydrodynamic Stability Analysis of Burning Bubbles in Electroweak Theory and in QCD
Authors:
P. Huet,
K. Kajantie,
R. G. Leigh,
B. -H. Liu,
L. McLerran
Abstract:
Assuming that the electroweak and QCD phase transitions are first order, upon supercooling, bubbles of the new phase appear. These bubbles grow to macroscopic sizes compared to the natural scales associated with the Compton wavelengths of particle excitations. They propagate by burning the old phase into the new phase at the surface of the bubble. We study the hydrodynamic stability of the burni…
▽ More
Assuming that the electroweak and QCD phase transitions are first order, upon supercooling, bubbles of the new phase appear. These bubbles grow to macroscopic sizes compared to the natural scales associated with the Compton wavelengths of particle excitations. They propagate by burning the old phase into the new phase at the surface of the bubble. We study the hydrodynamic stability of the burning and find that for the velocities of interest for cosmology in the electroweak phase transition, the shape of the bubble wall is stable under hydrodynamic perturbations. Bubbles formed in the cosmological QCD phase transition are found to be a borderline case between stability and instability.
△ Less
Submitted 10 December, 1992; v1 submitted 4 December, 1992;
originally announced December 1992.
-
Aspects of the Electroweak Phase Transition
Authors:
Patrick Huet
Abstract:
Presented at the 1992 Meeting of the DPF, Fermilab. The electroweak phase transition is reviewed in light of some recent developments. Emphasis is on the issue whether the transition is first or second order and its possible role in the generation of the baryon asymmetry of the universe.
Presented at the 1992 Meeting of the DPF, Fermilab. The electroweak phase transition is reviewed in light of some recent developments. Emphasis is on the issue whether the transition is first or second order and its possible role in the generation of the baryon asymmetry of the universe.
△ Less
Submitted 7 December, 1992; v1 submitted 3 December, 1992;
originally announced December 1992.
-
Towards the Theory of Cosmological Phase Transitions
Authors:
M. Dine,
R. Leigh,
P. Huet,
A. Linde,
D. Linde
Abstract:
We discuss recent progress (and controversies) in the theory of finite temperature phase transitions. This includes the structure of the effective potential at a finite temperature, the infrared problem in quantum statistics of gauge fields, the theory of formation of critical and subcritical bubbles and the theory of bubble wall propagation.
We discuss recent progress (and controversies) in the theory of finite temperature phase transitions. This includes the structure of the effective potential at a finite temperature, the infrared problem in quantum statistics of gauge fields, the theory of formation of critical and subcritical bubbles and the theory of bubble wall propagation.
△ Less
Submitted 3 March, 1992;
originally announced March 1992.
-
Comments on the Electroweak Phase Transition
Authors:
M. Dine,
R. G. Leigh,
P. Huet,
A. D. Linde,
D. A. Linde
Abstract:
We report on an investigation of various problems related to the theory of the electroweak phase transition. This includes a determination of the nature of the phase transition, a discussion of the possible role of higher order radiative corrections and the theory of the formation and evolution of the bubbles of the new phase. We find in particular that no dangerous linear terms appear in the ef…
▽ More
We report on an investigation of various problems related to the theory of the electroweak phase transition. This includes a determination of the nature of the phase transition, a discussion of the possible role of higher order radiative corrections and the theory of the formation and evolution of the bubbles of the new phase. We find in particular that no dangerous linear terms appear in the effective potential. However, the strength of the first order phase transition is 2/3 times less than what follows from the one-loop approximation. This rules out baryogenesis in the minimal version of the electroweak theory.
△ Less
Submitted 24 February, 1992;
originally announced March 1992.