The tidal stripping of satellites

We present an improved analytic calculation for the tidal radius of satellites and test our results against N -body simulations.
The tidal radius in general depends upon four factors: the potential of the host galaxy, the potential of the satellite, the orbit of the satellite and the orbit of the star within the satellite . We demonstrate that this last point is critical and suggest using three tidal radii to cover the range of orbits of stars within the satellite. In this way we show explicitly that prograde star orbits will be more easily stripped than radial orbits; while radial orbits are more easily stripped than retrograde ones. This result has previously been established by several authors numerically, but can now be understood analytically. For point mass, power-law (which includes the isothermal sphere), and a restricted class of split power-law potentials our solution is fully analytic. For more general potentials, we provide an equation which may be rapidly solved numerically.
Over short times (≲1–2 Gyr ∼1 satellite orbit), we find excellent agreement between our analytic and numerical models. Over longer times, star orbits within the satellite are transformed by the tidal field of the host galaxy. In a Hubble time, this causes a convergence of the three limiting tidal radii towards the prograde stripping radius. Beyond the prograde stripping radius, the velocity dispersion will be tangentially anisotropic.     

Single star scattering on an open cluster

We examine the scattering of single stars from an open star cluster. The probability of the capture of a star by a star cluster is dependent on the velocity and mass of the star, and the stars that are not captured experience a velocity change. For low-velocity stars there is an exponential decrease of the capture probability with the initial velocity, and the velocity change decreases almost linearly. For high-velocity stars there is a v ⁻⁶ dependence for the capture probability, and a v ⁻¹ dependence for the velocity change. Analytical estimations, Monte Carlo and full N -body simulations are all in good agreement.     

Scaling of N-body calculations

We report results of collisional N -body simulations aimed at studying the N dependence of the dynamical evolution of star clusters. Our clusters consist of equal-mass stars and are in virial equilibrium. Clusters moving in external tidal fields and clusters limited by a cut-off radius are simulated. Our main focus is to study the dependence of the lifetimes of the clusters on the number of cluster stars and the chosen escape condition.
We find that star clusters in external tidal fields exhibit a scaling problem in the sense that their lifetimes do not scale with the relaxation time. Isolated clusters show a similar problem if stars are removed only after their distance to the cluster centre exceeds a certain cut-off radius. If stars are removed immediately after their energy exceeds the energy necessary for escape, the scaling problem disappears.
We show that some stars that gain the energy necessary for escape are scattered to lower energies before they can leave the cluster. As the efficiency of this process decreases with increasing particle number, it causes the lifetimes not to scale with the relaxation time. Analytic formulae are derived for the scaling of the lifetimes in the different cases.     

On the structure of tidal tails

We examine the longitudinal distribution of the stars escaping from a cluster along tidal tails. Using both theory and simulations, we show that, even in the case of a star cluster in a circular galactic orbit, when the tide is steady, the distribution exhibits maxima at a distance of many tidal radii from the cluster.     

Winds from clusters with non-uniform stellar distributions

We present the analytic and numerical models of the 'cluster wind' resulting from the multiple interactions of the winds ejected by the stars of a dense cluster of massive stars. We consider the case in which the distribution of stars (i.e. the number of stars per unit volume) within the cluster is spherically symmetric, has a power-law radial dependence, and drops discontinuously to zero at the outer radius of the cluster. We carry out comparisons between an analytic model (in which the stars are considered in terms of a spatially continuous injection of mass and energy) and 3D gasdynamic simulations (in which we include 100 stars with identical winds, located in 3D space by statistically sampling the stellar distribution function). From the analytic model, we find that for stellar distributions with steep enough radial dependencies, the cluster wind flow develops a very high central density and a non-zero central velocity, and for steeper dependencies, it becomes fully supersonic throughout the volume of the cluster (these properties are partially reproduced by the 3D numerical simulations). Therefore, the wind solutions obtained for stratified clusters can differ dramatically from the case of a homogeneous stellar distribution (which produces a cluster wind with zero central velocity, and a fully subsonic flow within the cluster radius). Finally, from our numerical simulations, we compute predictions of X-ray emission maps and luminosities, which can be directly compared with observations of cluster wind flows.     

Blue straggler stars in dwarf spheroidal galaxies – II. Sculptor and Fornax

The existence of blue straggler stars (BSSs) in dwarf spheroidal galaxies (dSphs) is still an open question. In fact, many BSS candidates have been observed in the Local Group dSphs, but it is unclear whether they are real BSSs or young stars. Shedding light on the nature of these BSS candidates is crucial in order to understand the star formation history of dSphs. In this paper, we consider BSS candidates in Sculptor and Fornax. In Fornax, there are strong hints that the BSS population is contaminated by young stars, whereas in Sculptor there is no clear evidence of recent star formation. We derive the radial and luminosity distribution of BSS candidates from wide field imaging data extending beyond the nominal tidal radius of these galaxies. The observations are compared with the radial distribution of BSSs expected from dynamical simulations. In Sculptor, the radial distribution of BSS candidates is consistent with that of red horizontal branch (RHB) stars and is in agreement with theoretical expectations for BSSs generated via mass transfer in binaries. On the contrary, in Fornax, the radial distribution of BSS candidates is more concentrated than that of all the considered stellar populations. This result supports the hypothesis that most of BSS candidates in Fornax are young stars, and this is consistent with previous studies.     

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.     

The time-scale of escape from star clusters

In this paper a cluster is modelled as a smooth potential (due to the cluster stars) plus the steady tidal field of the Galaxy. In this model there is a minimum energy below which stars cannot escape. Above this energy, however, the time-scale on which a star escapes varies with the orbital parameters of the star (mainly its energy) in a way which we attempt to quantify, with both theoretical arguments and computer simulations. Within the limitations of the model we show that the time-scale is long enough to complicate the interpretation of full N -body simulations of clusters, and that stars above the escape energy may remain bound to the cluster for about a Hubble time.     

Resonant tidal disruption in galactic nuclei

It has recently been shown by Rauch 38 Tremaine that the rate of angular momentum relaxation in nearly Keplerian star clusters is greatly increased by a process termed 'resonant relaxation'; it was also argued, via a series of scaling arguments, that tidal disruption of stars in galactic nuclei containing massive black holes could be noticeably enhanced by this process. We describe here the results of numerical simulations of resonant tidal disruption which quantitatively test the predictions made by Rauch 38 Tremaine. The simulation method is based on an N -body routine incorporating cloning of stars near the loss cone and a semirelativistic symplectic integration scheme. Normalized disruption rates for resonant and non-resonant nuclei are derived at orbital energies both above and below the critical energy, and the corresponding angular momentum distribution functions are found. The black hole mass above which resonant tidal disruption is quenched by relativistic precession is determined. We also briefly describe the discovery of chaos in the Wisdom–Holman symplectic integrator applied to highly eccentric orbits and propose a modified integration scheme that remains robust under these conditions. We find that resonant disruption rates exceed their non-resonant counterparts by an amount consistent with the predictions; in particular, we estimate the net tidal disruption rate for a fully resonant cluster to be about twice that of its non-resonant counterpart. No significant enhancement in rates is observed outside the critical radius. Relativistic quenching of the effect is found to occur for hole masses M  >  M _Q  = (8 ± 3) × 10⁷  M . The numerical results combined with the observed properties of galactic nuclei indicate that for most galaxies the resonant enhancement to tidal disruption rates will be very small.     

A highly abnormal massive star mass function in the Orion Nebula cluster and the dynamical decay of trapezium systems

The Orion Nebula cluster (ONC) appears to be unusual on two grounds: the observed constellation of the OB stars of the entire ONC and its Trapezium at its centre implies a time-scale problem given the age of the Trapezium, and an initial mass function (IMF) problem for the whole OB star population in the ONC. Given the estimated crossing time of the Trapezium, it ought to have totally dynamically decayed by now. Furthermore, by combining the lower limit of the ONC mass with a standard IMF it emerges that the ONC should have formed at least about 40 stars heavier than  5 M_⊙  while only 10 are observed. Using the N -body experiments we (i) confirm the expected instability of the Trapezium and (ii) show that beginning with a compact OB-star configuration of about 40 stars both the number of observed OB stars after 1 Myr within 1 pc radius and a compact trapezium configuration can be reproduced. These two empirical constraints thus support our estimate of 40 initial OB stars in the cluster. Interestingly, a more-evolved version of the ONC resembles the Upper Scorpius OB association. The N -body experiments are performed with the new C-code catena by integrating the equations of motion using the chain-multiple-regularization method. In addition, we present a new numerical formulation of the IMF.     

The eccentricities of the barium stars

We investigate the eccentricities of barium (Ba  ii ) stars formed via a stellar wind accretion model. We carry out a series of Monte Carlo simulations using a rapid binary evolution algorithm, which incorporates full tidal evolution, mass loss and accretion, and nucleosynthesis and dredge-up on the thermally pulsing asymptotic giant branch. We follow the enhancement of barium in the envelope of the accreting main-sequence companion and dilution into its convective envelope once the star ascends the giant branch.
The observed eccentricities of Ba  ii stars are significantly smaller than those of an equivalent set of normal red giants but are nevertheless non-zero. We show that such a distribution of eccentricities is consistent with a wind accretion model for Ba  ii star production with weak viscous tidal dissipation in the convective envelopes of giant stars. We successfully model the distribution of orbital periods and the number of observed Ba  ii stars. The actual distribution of eccentricities is quite sensitive to the strength of the tides, so that we are able to confirm that this strength is close to, but less than, what is expected theoretically and found with alternative observational tests. Two systems – one very short-period but eccentric, and one long-period and highly eccentric – still lie outside the envelope of our models, and so require a more exotic formation mechanism. All our models, even those which were a good fit to the observed distributions, overproduced the number of high-period barium stars, a problem that could not be solved by some combination of the three parameters: tidal strength, tidal enhancement and wind accretion efficiency.     

Tidal disruption of globular clusters in dwarf galaxies with triaxial dark matter haloes

We use N -body simulations to study the tidal evolution of globular clusters (GCs) in dwarf spheroidal (dSph) galaxies. Our models adopt a cosmologically motivated scenario in which the dSph is approximated by a static Navarro, Frenk & White halo with a triaxial shape. We apply our models to five GCs spanning three orders of magnitude in stellar density and two in mass, chosen to represent the properties exhibited by the five GCs of the Fornax dSph. We show that only the object representing Fornax's least dense GC (F1) can be fully disrupted by Fornax's internal tidal field – the four denser clusters survive even if their orbits decay to the centre of Fornax. For a large set of orbits and projection angles, we examine the spatial and velocity distribution of stellar debris deposited during the complete disruption of an F1-like GC. Our simulations show that such debris appears as shells, isolated clumps and elongated overdensities at low surface brightness (≥26 mag arcsec⁻²), reminiscent of substructure observed in several Milky Way dSphs. Such features arise from the triaxiality of the galaxy potential and do not dissolve in time. The kinematics of the debris depends strongly on the progenitor's orbit. Debris associated with box and resonant orbits does not display stream motions and may appear 'colder'/'hotter' than the dSph's field population if the viewing angle is perpendicular/parallel to the progenitor's orbital plane. In contrast, debris associated with loop orbits shows a rotational velocity that may be detectable out to a few kpc from the galaxy centre. Chemical tagging that can distinguish GC debris from field stars may reveal whether the merger of GCs contributed to the formation of multiple stellar components observed in dSphs.     

Effects of external tidal field on the evolution of the outer regions of multi-mass star clusters

We present N -body simulations (including an initial mass function) of globular clusters in the Galaxy in order to study effects of the tidal field systematically on the properties of the outer parts of globular clusters. Using nbody6 , which correctly takes into account the two-body relaxation, we investigate the development of tidal tails of globular clusters in the Galactic tidal field. For simplicity, we have employed only the spherical components (bulge and halo) of the Galaxy, and ignored the effects of stellar evolution which could have been important in the very early phase of the cluster evolution. The total number of stars in our simulations is about 20 000, which is much smaller than the realistic number of stars. All simulations had been done for several orbital periods in order to understand the development of the tidal tails. In our scaled-down models, the relaxation time is sufficiently short to show the mass segregation effect, but we did not go far enough to see the core collapse, and the fraction of stars lost from the cluster at the end of the simulations is only ∼10 per cent. The radial distribution of extra-tidal stars can be described by a power law with a slope around −3 in surface density. The directions of tidal tails are determined by the orbits and locations of the clusters. We find that the length of tidal tails increases towards the apogalacticon and decreases towards the perigalacticon. This is an anti-correlation with the strength of the tidal field, caused by the fact that the time-scale for the stars to respond to the potential is similar to the orbital time-scale of the cluster. The escape of stars in the tidal tails towards the pericentre could be another reason for the decrease of the length of tidal tails. We find that the rotational angular velocity of tidally induced clusters shows quite different behaviour from that of initially rotating clusters.     

Monte Carlo simulations of star clusters – III. A million-body star cluster

A revision of Stodółkiewicz's Monte Carlo code is used to simulate the evolution of million-body star clusters. The new method treats each superstar as a single star and follows the evolution and motion of all individual stellar objects. The evolution of N -body systems influenced by the tidal field of a parent galaxy and by stellar evolution is presented. All models consist of 1 000 000 stars. 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. During the evolution, the initial mass function (IMF) changes significantly. The local mass function around the half-mass radius closely resembles the actual global mass function. At the late stages of evolution, the mass of the evolved stars inside the core can be as high as 97 per cent of the total mass in this region. For the whole system, the evolved stars can compose up to 75 per cent of the total mass. The evolution of cluster anisotropy strongly depends on initial cluster concentration, IMF and the strength of the tidal field. The results presented are the first step in the direction of simulating the evolution of real globular clusters by means of the Monte Carlo method.     

Numerical simulations of the Hyades dynamics and the nature of the moving Hyades cluster

We present our numerical simulations of the dynamical evolution of the Hyades open cluster. The simulations were performed usinga modified NBODY6 algorithm that included tidal forces and a realistic orbit of the cluster in a gravitational field described by the Miyamoto-Nagai potential. Our goal was to study the nature of movingclu sters. We show that the stars that were earlier cluster members could be later identified within a sphere of 50 pc in diameter around the Sun. The number of such stars for the chosen initial mass and virial radius of the cluster does not exceed ten. The maximum space velocity of these stars relative to the core of the current cluster does not exceed 3 km s^?1. Our numerical simulations confirm the assumption that some of the moving clusters near the Sun could consist of stars that have escaped from open clusters in the course of their dynamical evolution.     

Probing FSR star cluster candidates in bulge/disc directions with 2MASS colour–magnitude diagrams

We analyse 20 star cluster candidates projected mostly in the bulge direction  (|ℓ| < 60°)  . The sample contains all candidates in that sector classified by Froebrich, Scholz & Raftery with quality flags denoting high probability of being star clusters. Bulge contamination in the colour–magnitude diagrams (CMDs) is in general important, while at lower Galactic latitudes disc stars contribute as well. Properties of the candidates are investigated with Two Micron All Sky Survey (2MASS) CMDs and stellar radial density profiles (RDPs) built with field star decontaminated photometry. To uncover the nature of the structures we decontaminate the CMDs from field stars using tools that we previously developed to deal with objects in dense fields. We confirm in all cases excesses in the RDPs with respect to the background level, as expected from the method the candidates were originally selected. CMDs and RDPs taken together revealed six open clusters, five uncertain cases that require deeper observations, while nine objects are possibly field density fluctuations.     

On the dissolution of star clusters in the Galactic Centre – I. Circular orbits

We present N -body simulations of dissolving star clusters close to Galactic Centres. For this purpose, we developed a new N -body program called nbody6gc based on Aarseth's series of N -body codes. We describe the algorithm in detail. We report about the density wave phenomenon in the tidal arms which has been recently explained by Küpper, Macleod & Heggie. Standing waves develop in the tidal arms. The wave knots or clumps develop at the position, where the emerging tidal arm hits the potential wall of the effective potential and is reflected. The escaping stars move through the wave knots further into the tidal arms. We show the consistency of the positions of the wave knots with the theory in Just et al. We also demonstrate a simple method to study the properties of tidal arms. By solving many eigenvalue problems along the tidal arms, we numerically construct a one-dimensional coordinate system whose direction is always along a principal axis of the local tensor of inertia. Along this coordinate system, physical quantities can be evaluated. The half-mass or dissolution times of our models are almost independent of the particle number which indicates that two-body relaxation is not the dominant mechanism leading to the dissolution. This may be a typical situation for many young star clusters. We propose a classification scheme which sheds light on the dissolution mechanism.     

Seismic signatures of strange stars with crust

We study acoustic oscillations (eigenfrequencies, velocity distributions, damping times) of normal crusts of strange stars. These oscillations are very specific because of huge density jump at the interface between the normal crust and the strange matter core. The oscillation problem is shown to be self-similar. For a low (but non-zero) multipolarity l , the fundamental mode (without radial nodes) has a frequency of ∼300 Hz and mostly horizontal oscillation velocity; other pressure modes have frequencies ≳20 kHz and almost radial oscillation velocities. The latter modes are similar to radial oscillations (having approximately the same frequencies and radial velocity profiles). The oscillation spectrum of strange stars with crust differs from the spectrum of neutron stars. If detected, acoustic oscillations would allow one to discriminate between strange stars with crust and neutron stars and constrain the mass and radius of the star.     

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.