On the origin of high-velocity runaway stars

We explore the hypothesis that some high-velocity runaway stars attain their peculiar velocities in the course of exchange encounters between hard massive binaries and a very massive star (either an ordinary  50–100 M_⊙  star or a more massive one, formed through runaway mergers of ordinary stars in the core of a young massive star cluster). In this process, one of the binary components becomes gravitationally bound to the very massive star, while the second one is ejected, sometimes with a high speed. We performed three-body scattering experiments and found that early B-type stars (the progenitors of the majority of neutron stars) can be ejected with velocities of  ≳200–400 km s⁻¹  (typical of pulsars), while  3–4 M_⊙  stars can attain velocities of  ≳300–400 km s⁻¹  (typical of the bound population of halo late B-type stars). We also found that the ejected stars can occasionally attain velocities exceeding the Milky Ways's escape velocity.     

Multimass spherical structure models for N-body simulations

We present a simple and efficient method to set up spherical structure models for N -body simulations with a multimass technique. This technique reduces by a substantial factor the computer run time needed in order to resolve a given scale as compared to single-mass models. It therefore allows to resolve smaller scales in N -body simulations for a given computer run time. Here, we present several models with an effective resolution of up to  1.68 × 10⁹  particles within their virial radius which are stable over cosmologically relevant time-scales. As an application, we confirm the theoretical prediction by Dehnen that in mergers of collisionless structures like dark matter haloes always the cusp of the steepest progenitor is preserved. We model each merger progenitor with an effective number of particles of approximately 10⁸ particles. We also find that in a core–core merger the central density approximately doubles whereas in the cusp–cusp case the central density only increases by approximately 50 per cent. This may suggest that the central regions of flat structures are better protected and get less energy input through the merger process.     

The formation of a bound star cluster: from the Orion nebula cluster to the Pleiades

Direct N -body calculations are presented of the formation of Galactic clusters using GasEx , which is a variant of the code Nbody6 . The calculations focus on the possible evolution of the Orion nebula cluster (ONC) by assuming that the embedded OB stars explosively drove out 2/3 of its mass in the form of gas about 0.4 Myr ago. A bound cluster forms readily and survives for 150 Myr despite additional mass loss from the large number of massive stars, and the Galactic tidal field. This is the very first time that cluster formation is obtained under such realistic conditions. The cluster contains about 1/3 of the initial 10⁴ stars, and resembles the Pleiades cluster to a remarkable degree, implying that an ONC-like cluster may have been a precursor of the Pleiades. This scenario predicts the present expansion velocity of the ONC, which will be measurable by upcoming astrometric space missions. These missions should also detect the original Pleiades members as an associated expanding young Galactic-field subpopulation. The results arrived at here suggest that Galactic clusters form as the nuclei of expanding OB associations.
The results have wide implications, also for the formation of globular clusters and the Galactic-field and halo stellar populations. In view of this, the distribution of binary orbital periods and the mass function within and outside the model ONC and Pleiades is quantified, finding consistency with observational constraints. Advanced mass segregation is evident in one of the ONC models. The calculations show that the primordial binary population of both clusters could have been much the same as is observed in the Taurus–Auriga star-forming region. The computations also demonstrate that the binary proportion of brown dwarfs is depleted significantly for all periods, whereas massive stars attain a high binary fraction.     

N-body simulations of the Magellanic stream

A suite of high-resolution N -body simulations of the Magellanic Clouds–Milky Way system are presented and compared directly with newly available data from the H  i Parkes All-Sky Survey (HIPASS). We show that the interaction between Small Magellanic Clouds (SMC) and Large Magellanic Clouds results in both a spatial and kinematical bifurcation of both the stream and the leading arm. The spatial bifurcation of the stream is readily apparent in the HIPASS data, and the kinematical bifurcation is also tentatively identified. This bifurcation provides strong support for the tidal disruption origin for the Magellanic stream. A fiducial model for the Magellanic Clouds (MCs) is presented upon completion of an extensive parameter survey of the potential orbital configurations of the MCs and the viable initial boundary conditions for the disc of the SMC. The impact of the choice of these critical parameters upon the final configurations of the stream and leading arm is detailed.     

Post-Newtonian N-body simulations

We report on the first fully consistent conventional cluster simulation which includes terms up to the third-order post-Newtonian approximation. Numerical problems for treating extremely energetic binaries orbiting a single massive object are circumvented by employing the special 'wheel-spoke' regularization method of Zare which has not been used in large- N simulations before. Idealized models containing   N = 1 × 10⁵  particles of mass  1 M_⊙  with a central black hole (BH) of  300 M_⊙  have been studied on GRAPE-type computers. An initial half-mass radius of   r _h≃ 0.1  pc is sufficiently small to yield examples of relativistic coalescence. This is achieved by significant binary shrinkage within a density cusp environment, followed by the generation of extremely high eccentricities which are induced by Kozai cycles and/or resonant relaxation. More realistic models with white dwarfs and 10 times larger half-mass radii also show evidence of general relativity effects before disruption. An experimentation with the post-Newtonian terms suggests that reducing the time-scales for activating the different orders progressively may be justified for obtaining qualitatively correct solutions without aiming for precise predictions of the final gravitational radiation wave form. The results obtained suggest that the standard loss-cone arguments underestimate the swallowing rate in globular clusters containing a central BH.     

Multimass schemes for collisionless N-body simulations

We present a general scheme for constructing Monte Carlo realizations of equilibrium, collisionless galaxy models with known distribution function (DF) f ₀. Our method uses importance sampling to find the sampling DF f _s that minimizes the mean-square formal errors in a given set of projections of the DF f ₀. The result is a multimass N -body realization of the galaxy model in which 'interesting' regions of phase space are densely populated by lots of low-mass particles, increasing the effective N there, and less interesting regions by fewer, higher mass particles.
As a simple application, we consider the case of minimizing the shot noise in estimates of the acceleration field for an N -body model of a spherical Hernquist model. Models constructed using our scheme easily yield a factor of ∼100 reduction in the variance at the central acceleration field when compared to a traditional equal-mass model with the same number of particles. When evolving both models with a real N -body code, the diffusion coefficients in our model are reduced by a similar factor. Therefore, for certain types of problems, our scheme is a practical method for reducing the two-body relaxation effects, thereby bringing the N -body simulations closer to the collisionless ideal.     

N-body simulations for testing the stability of triaxial galaxies in MOND

We perform a stability test of triaxial models in Modified Newtonian Dynamics (MOND) using N -body simulations. The triaxial models considered here have densities that vary with   r ⁻¹  in the centre and   r ⁻⁴  at large radii. The total mass of the model varies from 10⁸ to  10¹⁰ M_⊙  , representing the mass scale of dwarfs to medium-mass elliptical galaxies, respectively, from deep MOND to quasi-Newtonian gravity. We build triaxial galaxy models using the Schwarzschild technique, and evolve the systems for 200 Keplerian dynamical times (at the typical length-scale of 1.0 kpc). We find that the systems are virial overheating, and in quasi-equilibrium with the relaxation taking approximately 5 Keplerian dynamical times (1.0 kpc). For all systems, the change of the inertial (kinetic) energy is less than 10 per cent (20 per cent) after relaxation. However, the central profile of the model is flattened during the relaxation and the (overall) axis ratios change by roughly 10 per cent within 200 Keplerian dynamical times (at 1.0 kpc) in our simulations. We further find that the systems are stable once they reach the equilibrium state.     

The nearest star of spectral type O3: a component of the multiple system HD 150136

From radial velocities determined in high signal-to-noise digital spectra, we report the discovery that the brightest component of the binary system HD 150136 is of spectral type O3. We also present the first double-lined orbital solution for this binary. Our radial velocities confirm the previously published spectroscopic orbital period of 2.6 d. He  ii absorptions appear double at quadratures, but single lines of N  v and N  iv visible in our spectra define a radial velocity orbit of higher semi-amplitude for the primary component than do the He  ii lines. From our orbital analysis, we obtain minimum masses for the binary components of 27 and  18 M_⊙  . The neutral He absorptions apparently do not follow the orbital motion of any of the binary components, thus they most probably arise in a third star in the system.     

Application of a probabilistic neural network in radial‐velocity curve analysis of the spectroscopic binary stars HD 152218, HD 143511, HD 27149, and ER Vul

We use an Artificial Neural Network (ANN) to derive the orbital parameters of spectroscopic binary stars. Using measured radial velocity data of four double‐lined spectroscopic binary systems HD 152218, HD 143511, HD 27149, and ER Vul, we find corresponding orbital and spectroscopic elements. Our numerical results are in good agreement with those obtained by others using more traditional methods (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The separation of the stars in the binary nucleus of the planetary nebula Abell 35

Using the Planetary Camera on board the Hubble Space Telescope , we have measured the projected separation of the binary components in the nucleus of the planetary nebula Abell 35 to be larger than 0.08 arcsec but less than 0.14 arcsec. The system has been imaged in three filters centred at 2950, 3350 and 5785 Å. The white dwarf primary star responsible for ionizing the nebula is half as bright as its companion in the 2950-Å filter, causing the source to be visibly elongated. The 3350-Å setting, on the other hand, shows no elongation as a result of the more extreme flux ratio. The F300W data allow the determinination of the projected separation of the binary. At the minimum distance of 160 pc to the system, our result corresponds to 18 ± 5 au. This outcome is consistent with the wind accretion induced rapid rotation hypothesis, but cannot be reconciled with the binary having emerged from a common-envelope phase.     

Accurate estimators of power spectra in N-body simulations

A method to rapidly estimate the Fourier power spectrum of a point distribution is presented. This method relies on a Taylor expansion of the trigonometric functions. It yields the Fourier modes from a number of fast Fourier transforms (FFTs), which is controlled by the order N of the expansion and by the dimension D of the system. In three dimensions, for the practical value   N = 3  , the number of FFTs required is 20.
We apply the method to the measurement of the power spectrum of a periodic point distribution that is a local Poisson realization of an underlying stationary field. We derive an explicit analytic expression for the spectrum, which allows us to quantify – and correct for – the biases induced by discreteness and by the truncation of the Taylor expansion, and to bound the unknown effects of aliasing of the power spectrum. We show that these aliasing effects decrease rapidly with the order N . For   N = 3  , they are expected to be, respectively, smaller than  ∼10⁻⁴  and 0.02 at half the Nyquist frequency and at the Nyquist frequency of the grid used to perform the FFTs. The only remaining significant source of errors is reduced to the unavoidable cosmic/sample variance due to the finite size of the sample.
The analytical calculations are successfully checked against a cosmological N -body experiment. We also consider the initial conditions of this simulation, which correspond to a perturbed grid. This allows us to test a case where the local Poisson assumption is incorrect. Even in that extreme situation, the third-order Fourier–Taylor estimator behaves well, with aliasing effects restrained to at most the per cent level at half the Nyquist frequency.
We also show how to reach arbitrarily large dynamic range in Fourier space (i.e. high wavenumber), while keeping statistical errors in control, by appropriately 'folding' the particle distribution.     

A search for wide visual companions of exoplanet host stars: The Calar Alto Survey

We have carried out a search for co‐moving stellar and substellar companions around 18 exoplanet host stars with the infrared camera MAGIC at the 2.2 m Calar Alto telescope, by comparing our images with images from the all sky surveys 2MASS, POSS I and II. Four stars of the sample namely HD80606, 55 Cnc, HD46375 and BD–10°3166, are listed as binaries in the Washington Visual Double Star Catalogue (WDS). The binary nature of HD80606, 55 Cnc, and HD46375 is confirmed with both astrometry as well as photometry, thereby the proper motion of the companion of HD46375 was determined here for the first time.We derived the companion masses as well as the longterm stability regions for additional companions in these three binary systems. We can rule out further stellar companions around all stars in the sample with projected separations between 270AU and 2500AU, being sensitive to substellar companions with masses down to ∼60 M_Jup (S /N = 3). Furthermore we present evidence that the two components of the WDS binary BD–10°3166 are unrelated stars, i.e this system is a visual pair. The spectrophotometric distance of the primary (a K0 dwarf) is ∼67 pc, whereas the presumable secondary BD–10°3166B (a M4 to M5 dwarf) is located at a distance of 13 pc in the foreground. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)