Search NASA⌕ Search

Engineering topics

Mcmillan, Steve

Publications and source records attributed to Mcmillan, Steve.

Building a better leapfrog

In stellar dynamical computer simulations, as well as other types of simulations using particles, time step size is often held constant in order to guarantee a high degree of energy conservation. In many applications, allowing the time step size to change in time can offer a great saving in computational cost, but variable-size time steps usually imply a substantial degradation in energy conservation. We present a meta-algorithm' for choosing time steps in such a way as to guarantee time symmetry in any integration scheme, thus allowing vastly improved energy conservation for orbital calculations with variable time steps. We apply the algorithm to the familiar leapfrog scheme, and generalize to higher order integration schemes, showing how the stability properties of the fixed-step leapfrog scheme can be extended to higher order, variable-step integrators such as the Hermite method. We illustrate the remarkable properties of these time-symmetric integrators for the case of a highly eccentric elliptical Kepler orbit and discuss applications to more complex problems.

Hut, Piet↗

Binary-binary interactions and the formation of the PSR B1620-26 triple system in M4

The hierarchical triple system containing the millisecond pulsar PSR B1620-26 in M4 is the first 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 pre-existing binary millisecond pulsar. This scenario would have the advantage of reconciling the approximately 10(exp 9) yr timing age of the pulsar with the much shorter lifetime of the triple system in the core of M4.

Rasio, Frederic A.↗

Star cluster evolution with primordial binaries. 3: Effect of the Galactic tidal field

We present the results of N-body simulations of tidally limited star clusters with an initial population of 0%-20% binaries. We find that (1) if enough binaries are initially present, the binary fraction may fall to a minimum value, then increase at late times; (2) the cluster evaporation timescale is quite insensitive to the details of the initial binary distribution; (3) the cluster core radius stabilizes at a few percent of the half-mass radius when binaries are present, just as in the case of isolated clusters; and (4) there may be a marked difference between the spatial distribution of low-energy and high-energy binaries as the cluster evolves. Specifically, the spatial distribution of the lower energy systems is often substantially more extended than that of the more tightly bound pairs. At no time are our simulated clusters well described by simple dynamical models that neglect the close coupling between the binding energies and the center-of-mass energies of the binaries they contain.

Mcmillan, Steve↗

Stellar black holes in globular clusters

The recent discovery of large populations of millisec pulsars associated with neutron stars in globular clusters indicates that several hundred stellar black holes of about 10 solar masses each can form within a typical cluster. While, in clusters of high central density, the rapid dynamical evolution of the black-hole population leads to an ejection of nearly all holes on a short timescale, systems of intermediate density may involve a normal star's capture by one of the surviving holes to form a low-mass X-ray binary. One or more such binaries may be found in the globular clusters surrounding our galaxy.

Kulkarni, S. R.↗

Binaries in globular clusters

Recent observations have shown that globular clusters contain a substantial number of binaries most of which are believed to be primordial. We discuss different successful optical search techniques, based on radial-velocity variables, photometric variables, and the positions of stars in the color-magnitude diagram. In addition, we review searches in other wavelengths, which have turned up low-mass X-ray binaries and more recently a variety of radio pulsars. On the theoretical side, we give an overview of the different physical mechanisms through which individual binaries evolve. We discuss the various simulation techniques which recently have been employed to study the effects of a primordial binary population, and the fascinating interplay between stellar evolution and stellar dynamics which drives globular-cluster evolution.

Hut, Piet↗

The evolution of a primordial binary population in a globular cluster

A simple model is presented for the evolution of a primordial binary population in a globular cluster. Monte Carlo simulations are given for an initial population of 50,000 binaries against a fixed background population of 500,000 single stars in a tidally truncated cluster model. Individual histories of all binaries are followed through mass segregation, scattering recoil, escape from the cluster, or coalescence. It is found that most binaries are destroyed by binary-binary interactions, with the rest escaping in the point-mass approximation. In a more realistic model, the majority of the rest merge. At any instant, most of the remaining binaries are drifting in toward the center before their first strong encounter. A typical binary spends most of its active life in or near the cluster core. The few binaries which receive a recoil sufficient to place them in the halo past the half-mass radius remain there long enough to make a significant contribution to the radial binary distribution.

Hut, Piet↗