Hyperfast pulsars as the remnants of massive stars ejected from young star clusters

Recent proper motion and parallax measurements for the pulsar PSR B1508+55 indicate a transverse velocity of  ∼1100 km s⁻¹  , which exceeds earlier measurements for any neutron star. The spin-down characteristics of PSR B1508+55 are typical for a non-recycled pulsar, which implies that the velocity of the pulsar cannot have originated from the second supernova disruption of a massive binary system. The high velocity of PSR B1508+55 can be accounted for by assuming that it received a kick at birth or that the neutron star was accelerated after its formation in the supernova explosion. We propose an explanation for the origin of hyperfast neutron stars based on the hypothesis that they could be the remnants of a symmetric supernova explosion of a high-velocity massive star which attained its peculiar velocity (similar to that of the pulsar) in the course of a strong dynamical three- or four-body encounter in the core of dense young star cluster. To check this hypothesis, we investigated three dynamical processes involving close encounters between: (i) two hard massive binaries, (ii) a hard binary and an intermediate-mass black hole (IMBH) and (iii) a single stars and a hard binary IMBH. We find that main-sequence O-type stars cannot be ejected from young massive star clusters with peculiar velocities high enough to explain the origin of hyperfast neutron stars, but lower mass main-sequence stars or the stripped helium cores of massive stars could be accelerated to hypervelocities. Our explanation for the origin of hyperfast pulsars requires a very dense stellar environment of the order of  10⁶– 10⁷ stars pc⁻³  . Although such high densities may exist during the core collapse of young massive star clusters, we caution that they have never been observed.     

Stellar dynamics in young clusters: the formation of massive runaways and very massive runaway mergers

In the present paper we combine an N-body code that simulates the dynamics of young dense stellar systems with a massive star evolution handler that accounts in a realistic way for the effects of stellar wind mass loss. We discuss two topics.

1. The formation and the evolution of very massive stars (with masses >120 M_⊙) is followed in detail. These very massive stars are formed in the cluster core as a consequence of the successive (physical) collisions of the 10–20 most massive stars in the cluster (this process is known as ‘runaway merging’). The further evolution is governed by stellar wind mass loss during core hydrogen and core helium burning (the WR phase of very massive stars). Our simulations reveal that, as a consequence of runaway merging in clusters with solar and supersolar values, massive black holes can be formed, but with a maximum mass ≈70 M_⊙. In low-metallicity clusters, however, it cannot be excluded that the runaway-merging process is responsible for pair-instability supernovae or for the formation of intermediate-mass black holes with a mass of several 100 M_⊙.
2. Massive runaways can be formed via the supernova explosion of one of the components in a binary system (the Blaauw scenario), or via dynamical interaction of a single star and a binary or between two binaries in a star cluster. We explore the possibility that the most massive runaways (e.g. ζ Pup, λ Cep, BD+43°3654) are the product of the collision and merger of two or three massive stars.

     

Close encounters in young stellar clusters: implications for planetary systems in the solar neighbourhood

The stars that populate the solar neighbourhood were formed in stellar clusters. Through N -body simulations of these clusters, we measure the rate of close encounters between stars. By monitoring the interaction histories of each star, we investigate the singleton fraction in the solar neighbourhood. A singleton is a star which formed as a single star, has never experienced any close encounters with other stars or binaries, or undergone an exchange encounter with a binary. We find that, of the stars which formed as single stars, a significant fraction is not singletons once the clusters have dispersed. If some of these stars had planetary systems, with properties similar to those of the Solar System, the planets' orbits may have been perturbed by the effects of close encounters with other stars or the effects of a companion star within a binary. Such perturbations can lead to strong planet–planet interactions which eject several planets, leaving the remaining planets on eccentric orbits. Some of the single stars exchange into binaries. Most of these binaries are broken up via subsequent interactions within the cluster, but some remain intact beyond the lifetime of the cluster. The properties of these binaries are similar to those of the observed binary systems containing extrasolar planets. Thus, dynamical processes in young stellar clusters will alter significantly any population of Solar System-like planetary systems. In addition, beginning with a population of planetary systems exactly resembling the Solar System around single stars, dynamical encounters in young stellar clusters may produce at least some of the extrasolar planetary systems observed in the solar neighbourhood.     

On the formation of massive stars

We present a model for the formation of massive ( M ≳10 M⊙) stars through accretion-induced collisions in the cores of embedded dense stellar clusters. This model circumvents the problem of accreting on to a star whose luminosity is sufficient to reverse the infall of gas. Instead, the central core of the cluster accretes from the surrounding gas, thereby decreasing its radius until collisions between individual components become sufficient. These components are, in general, intermediate-mass stars that have formed through accretion on to low-mass protostars. Once a sufficiently massive star has formed to expel the remaining gas, the cluster expands in accordance with this loss of mass, halting further collisions. This process implies a critical stellar density for the formation of massive stars, and a high rate of binaries formed by tidal capture.     

Ratios of star cluster core and half-mass radii: a cautionary note on intermediate-mass black holes in star clusters

There is currently much interest in the possible presence of intermediate-mass black holes (IMBHs) in the cores of globular clusters (GCs). Based on theoretical arguments and simulation results it has previously been suggested that a large core radius – or particularly a large ratio of the core radius to half-mass radius – is a promising indicator for finding such a black hole (BH) in a star cluster. In this study N -body models of 100 000 stars with and without primordial binaries are used to investigate the long-term structural evolution of star clusters. Importantly, the simulation data are analysed using the same processes by which structural parameters are extracted from observed star clusters. This gives a ratio of the core and half-mass (or half-light) radii that are directly comparable to the Galactic GC sample. As a result, it is shown that the ratios observed for the bulk of this sample can be explained without the need for an IMBH. Furthermore, it is possible that clusters with large core to half-light radius ratios harbour a BH binary (comprising stellar mass BHs) rather than a single massive BH. This work does not rule out the existence of IMBHs in the cores of at least some star clusters.     

The role of tidal interactions in star formation

Nearly all of the initial angular momentum of the matter that goes into each forming star must somehow be removed or redistributed during the formation process. The possible transport mechanisms and the possible fates of the excess angular momentum are discussed, and it is argued that transport processes in discs are probably not sufficient by themselves to solve the angular momentum problem, while tidal interactions with other stars in forming binary or multiple systems are likely to be of very general importance in redistributing angular momentum during the star formation process. Most, if not all, stars probably form in binary or multiple systems, and tidal torques in these systems can transfer much of the angular momentum from the gas around each forming star to the orbital motions of the companion stars. Tidally generated waves in circumstellar discs may contribute to the overall redistribution of angular momentum. Stars may gain much of their mass by tidally triggered bursts of rapid accretion, and these bursts could account for some of the most energetic phenomena of the earliest stages of stellar evolution, such as jet-like outflows. If tidal interactions are indeed of general importance, planet-forming discs may often have a more chaotic and violent early evolution than in standard models, and shock heating events may be common. Interactions in a hierarchy of subgroups may play a role in building up massive stars in clusters and in determining the form of the upper initial mass function (IMF) . Many of the processes discussed here have analogues on galactic scales, and there may be similarities between the formation of massive stars by interaction-driven accretion processes in clusters and the buildup of massive black holes in galactic nuclei.     

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.     

Hypervelocity collisions of binary stars at the Galactic Centre

Recent surveys have identified seven hypervelocity stars (HVSs) in the halo of the Milky Way. Most of these stars may have originated from the breakup of binary star systems by the nuclear black hole SgrA*. In some instances, the breakup of the binary may lead to a collision between its member stars. We examine the dynamical properties of these collisions by simulating thousands of different binary orbits around SgrA* with a direct N -body integration code. For some orbital parameters, the two stars collide with an impact velocity lower than their escape velocity and may therefore coalesce. It is possible for a coalescing binary to have sufficient velocity to escape the galaxy. Furthermore, some of the massive S-stars near Sgr A* might be the merger remnants of binary systems, however this production method can not account for most of the S-stars.     

Monte Carlo simulations of star clusters – II. Tidally limited, multimass systems with stellar evolution

A revision of Stodółkiewicz's Monte Carlo code is used to simulate evolution of large star clusters. The new method treats each superstar as a single star and follows the evolution and motion of all individual stellar objects. A survey of the evolution of N -body systems influenced by the tidal field of a parent galaxy and by stellar evolution is presented. The process of energy generation is realized by means of appropriately modified versions of Spitzer's and Mikkola's formulae for the interaction cross-section between binaries and field stars and binaries themselves. The results presented are in good agreement with theoretical expectations and the results of other methods (Fokker–Planck, Monte Carlo and N -body). The initial rapid mass loss, resulting from stellar evolution of the most massive stars, causes expansion of the whole cluster and eventually leads to the disruption of less bound systems ( W ₀=3). Models with larger W ₀ survive this phase of evolution and then undergo core collapse and subsequent post-collapse expansion, like isolated models. The expansion phase is eventually reversed when tidal limitation becomes important. The results presented are the first major step in the direction of simulating evolution of real globular clusters by means of the Monte Carlo method.     

Unexpected formation modes of the first hard binary in core collapse

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 the assumptions are incorrect in the majority of the cases: (1) quite a few three-body interactions produce a hard binary from scratch; (2) 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.     

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.     

On the mass of dense star clusters in starburst galaxies from spectrophotometry

The mass of unresolved young star clusters derived from spectrophotometric data may well be off by a factor of 2 or more once the migration of massive stars driven by mass segregation is accounted for. We quantify this effect for a large set of cluster parameters, including variations in the stellar initial mass function (IMF), the intrinsic cluster mass, and mean mass density. Gas-dynamical models coupled with the Cambridge stellar evolution tracks allow us to derive a scheme to recover the real cluster mass given measured half-light radius, one-dimensional velocity dispersion and age. We monitor the evolution with time of the ratio of real to apparent mass through the parameter η. When we compute η for rich star clusters, we find non-monotonic evolution in time when the IMF stretches beyond a critical cut-off mass of  25.5 M_⊙  . We also monitor the rise of colour gradients between the inner and outer volume of clusters: we find trends in time of the stellar IMF power indices overlapping well with those derived for the Large Magellanic Cloud cluster NGC 1818 at an age of 30 Myr. We argue that the core region of massive Antennae clusters should have suffered from much segregation despite their low ages. We apply these results to a cluster mass function, and find that the peak of the mass distribution would appear to observers shifted to lower masses by as much as 0.2 dex. The star formation rate derived for the cluster population is then underestimated by from 20 to 50 per cent.     

On the infant weight loss of low- to intermediate-mass star clusters

Star clusters are born in a highly compact configuration, typically with radii of less than about 1 pc roughly independently of mass. Since the star formation efficiency is less than 50 per cent by observation and because the residual gas is removed from the embedded cluster, the cluster must expand. In the process of doing so it only retains a fraction f _st of its stars. To date there are no observational constraints for f _st, although N -body calculations by Kroupa, Aarseth & Hurley suggest it to be about 20–30 per cent for Orion-type clusters. Here we use the data compiled by Testi et al., Testi, Palla & Natta and Testi, Palla & Natta for clusters around young Ae/Be stars and by de Wit et al. and de Wit et al. around young O stars and the study of de Zeeuw et al. of OB associations and combine these measurements with the expected number of stars in clusters with primary Ae/Be and O stars, respectively, using the empirical correlation between maximal stellar mass and star cluster mass of Weidner & Kroupa. We find that   f _st < 50  per cent with a decrease to higher cluster masses/more massive primaries. The interpretation would be that cluster formation is very disruptive. It appears that clusters with a birth stellar mass in the range  10–10³ M_⊙  keep at most 50 per cent of their stars.     

Stellar contents of two young open clusters: NGC 663 and 654

UBVRI CCD photometry in a wide field around two young open clusters, NGC 663 and 654, has been carried out. Hα and polarimetric observations for the cluster NGC 654 have also been obtained. We use the photometric data to construct colour–colour and colour–magnitude diagrams, from which we can investigate the reddening, age, mass and evolutionary states of the stellar contents of the these clusters. The reddening across the cluster regions is found to be variable. There is evidence for anomalous reddening law in both clusters; however, more infrared and polarimetric data are needed to conclude about the reddening law. Both clusters are situated at about a distance of 2.4 kpc. Star formation in both clusters is found to be a continuous process. In the case of NGC 663, star formation seems to have taken place sequentially, in the sense that formation of low-mass stars precedes the formation of most massive stars. Whereas, in the case of NGC 654, formation of low-mass stars did not cease after the formation of most massive stars in the cluster.     

大质量双星系统的非守恒演化

由于大质量双星系统有强大的星风物质损失,因而在研究其结构和演化时必须考虑星风物质损失,动量损失,物质交换以及由以上原因引起的轨道参量的变化,此外,天文观测又证实,一些大质量双星系统中存在星风冲击波,有Ｘ射线辐射以及有致密天体（白矮星,中子星）的存在,因此在研究大质量双星的演化时,又会遇到在星风冲击波理论及其对演化的影响,双星系统何时会演化成为公共外壳的系统,以及双星系统中如果发生超新星爆发,是否会     

Eclipsing Binaries in Open Clusters

The study of detached eclipsing binaries in open clusters can provide stringent tests of theoretical stellar evolutionary models, which must simultaneously fit the masses, radii, and luminosities of the eclipsing stars and the radiative properties of every other star in the cluster. We review recent progress in such studies and discuss two unusually interesting objects currently under analysis. GV Carinae is an A0 m + A8 m binary in the Southern open cluster NGC 3532; its eclipse depths have changed by 0.1 mag between 1990 and 2001, suggesting that its orbit is being perturbed by a relatively close third body. DW Carinae is a high-mass unevolved B1 V + B1 V binary in the very young open cluster Collinder 228, and displays double-peaked emission in the centre of the Hα line which is characteristic of Be stars. We conclude by pointing out that the great promise of eclipsing binaries in open clusters can only be satisfied when both the binaries and their parent clusters are well-observed, a situation which is less common than we would like.     

Neutron star retention and millisecond pulsar production in globular clusters

We investigate the conditions by which neutron star retention in globular clusters is favoured. We find that neutron stars formed in massive binaries are far more likely to be retained. Such binaries are likely to then evolve into contact before encountering other stars, possibly producing a single neutron star after a common envelope phase. A large fraction of the single neutron stars in globular clusters are then likely to exchange into binaries containing moderate-mass main-sequence stars, replacing the lower-mass components of the original systems. These binaries will become intermediate-mass X-ray binaries (IMXBs), once the moderate-mass star evolves off the main sequence, as mass is transferred on to the neutron star, possibly spinning it up in the process. Such systems may be responsible for the population of millisecond pulsars (MSPs) that has been observed in globular clusters. Additionally, the period of mass-transfer (and thus X-ray visibility) in the vast majority of such systems will have occurred 5–10 Gyr ago, thus explaining the observed relative paucity of X-ray binaries today, given the MSP population.     

Numerical simulations of protostellar encounters — III. Non-coplanar disc–disc encounters

It is expected that an average protostar will undergo at least one impulsive interaction with a neighbouring protostar whilst a large fraction of its mass is still in a massive, extended disc. If protostars are formed individually within a cluster before falling together and interacting, there should be no preferred orientation for such interactions. As star formation within clusters is believed to be coeval, it is probable that, during interactions, both protostars possess massive, extended discs.   We have used an SPH code to carry out a series of simulations of non-coplanar disc–disc interactions. We find that non-coplanar interactions trigger gravitational instabilities in the discs, which may then fragment to form new companions to the existing stars. (This is different from coplanar interactions, in which most of the new companion stars form after material in the discs has been swept up into a shock layer, and this then fragments.) The original stars may also capture each other, leading to the formation of a small- N cluster. If every star undergoes a randomly oriented disc–disc interaction, then the outcome will be the birth of many new stars and substellar objects. Approximately two-thirds of the stars will end up in multiple systems.     

A search for luminous Be stars

As Be stars are restricted to luminosity classes III‐V, but early B‐type stars are believed to evolve into supergiants, it is to be expected that the Be phenomenon disappears at some point in the evolution of a moderately massive star, before it reaches the supergiant phase. As a first stage in an attempt to determine the physical reasons of this cessation, a search of the literature has provided a number of candidates to be Be stars with luminosity classes Ib or II. Spectroscopy has been obtained for candidates in a number of open clusters and associations, as well as several other bright stars in those clusters. Among the objects observed, HD 207329 is the best candidate to be a high‐luminosity Be star, as it appears like a fast‐rotating supergiant with double‐peaked emission lines. The lines of HD 229059, in Berkeley 87, also appear morphologically similar to those of Be stars, but there are reasons to suspect that this object is an interacting binary. At slightly lower luminosities, LS I +56°92 (B4 II) and HD 333452 (O9 II), also appear as intrinsically luminous Be stars. Two Be stars in NGC 6913, HD 229221 and HD 229239, appear to have rather higher intrinsic magnitudes than their spectral type (B0.2 III in both cases) would indicate, being as luminous as luminosity class II objects in the same cluster. HD 344863, in NGC 6823, is also a rather early Be star of moderately high luminosity. The search shows that, though high‐luminosity Be stars do exist, they are scarce and, perhaps surprisingly, tend to have early spectral types. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Black holes and core expansion in massive star clusters

In this study we present the results from realistic N -body modelling of massive star clusters in the Magellanic Clouds. We have computed eight simulations with   N ∼ 10⁵  particles; six of these were evolved for at least a Hubble time. The aim of this modelling is to examine in detail the possibility of large-scale core expansion in massive star clusters, and search for a viable dynamical origin for the radius–age trend observed for such objects in the Magellanic Clouds. We identify two physical processes which can lead to significant and prolonged cluster core expansion – mass-loss due to rapid stellar evolution in a primordially mass-segregated cluster, and heating due to a retained population of stellar mass black holes, formed in the supernova explosions of the most massive cluster stars. These two processes operate over different time-scales and during different periods of a cluster's life. The former occurs only at early times and cannot drive core expansion for longer than a few hundred Myr, while the latter typically does not begin until several hundred Myr have passed, but can result in core expansion lasting for many Gyr. We investigate the behaviour of each of these expansion mechanisms under different circumstances – in clusters with varying degrees of primordial mass segregation, and in clusters with varying black hole retention fractions. In combination, the two processes can lead to a wide variety of evolutionary paths on the radius–age plane, which fully cover the observed cluster distribution and hence define a dynamical origin for the radius–age trend in the Magellanic Clouds. We discuss in some detail the implications of core expansion for various aspects of globular cluster research, as well as the possibility of observationally inferring the presence of a significant population of stellar mass black holes in a cluster.