首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 437 毫秒
1.
We consider the use of N -body simulations for studying the evolution of rich star clusters (i.e. globular clusters).The dynamical processes included in this study are restricted to gravitational (point-mass) interactions, the steady tidal field of a galaxy, and instantaneous mass loss resulting from stellar evolution. With evolution driven by these mechanisms, it is known that clusters fall roughly into two broad classes: those that dissipate promptly in the tidal field, as a result of mass loss; and those that survive long enough for their evolution to become dominated by two-body relaxation.
The time-scales of the processes we consider scale in different ways with the number of stars in the simulation, and the main aim of the paper is to suggest how the scaling of a simulation should be done so that the results are representative of the evolution of a 'real' cluster. We investigate three different ways of scaling time. One of these is appropriate to the first type of cluster, i.e. those that dissipate rapidly; similarly, a second scaling is appropriate only to the second (relaxation-dominated) type. We also develop a hybrid scaling, which is a satisfactory compromise for both types of cluster. Finally we present evidence that the widely used Fokker–Planck method produces models that are in good agreement with N -body models of those clusters that are relaxation-dominated, at least for N -body models with several thousand particles, but that the Fokker–Planck models evolve too fast for clusters that dissipate promptly.  相似文献   

2.
We present the first post-core-collapse models of initially rotating star clusters, using the numerical solution of an orbit-averaged 2D Fokker–Planck equation. Based on the code developed by Einsel & Spurzem, we have improved the speed and the stability and included the steady three-body binary heating source. We have confirmed that rotating clusters, whether they are in a tidal field or not, evolve significantly faster than non-rotating ones. Consequences for the observed shapes, density distribution and kinematic properties of young and old star clusters are discussed. The results are compared with gaseous and 1D Fokker–Planck models in the non-rotating case.  相似文献   

3.
In this paper we investigate the gravothermal instability of spherical stellar systems endowed with a radially anisotropic velocity distribution. We focus our attention on the effects of anisotropy on the conditions for the onset of instability and in particular we study the dependence of the spatial structure of critical models on the amount of anisotropy present in a system. The investigation has been carried out by the method of linear series which has already been used in the past to study the gravothermal instability of isotropic systems._   We consider models described by King, Wilson and Woolley–Dickens distribution functions. In the case of King and Woolley–Dickens models, our results show that, for quite a wide range of the amount of anisotropy in the system, the critical value of the concentration of the system (defined as the ratio of the tidal to the King core radius of the system) is approximately constant and equal to the corresponding value for isotropic systems. Only for very anisotropic systems does the critical value of the concentration start to change and it decreases significantly as the anisotropy increases and penetrates the inner parts of the system. For Wilson models the decrease of the concentration of critical models is preceded by an intermediate regime in which critical concentration increases, reaches a maximum and then starts to decrease. The critical value of the central potential always decreases as the anisotropy increases.  相似文献   

4.
Using the standard dynamical theory of spherical systems, we calculate the properties of spherical galaxies and clusters whose density profiles obey the universal form first obtained in high-resolution cosmological N -body simulations by Navarro, Frenk & White (NFW). We adopt three models for the internal kinematics: isotropic velocities, constant anisotropy and increasingly radial OsipkovMerritt anisotropy. Analytical solutions are found for the radial dependence of the mass, gravitational potential, velocity dispersion, energy and virial ratio and we test their variability with the concentration parameter describing the density profile and amount of velocity anisotropy. We also compute structural parameters, such as half-mass radius, effective radius and various measures of concentration. Finally, we derive projected quantities, the surface mass density and line-of-sight as well as aperture-velocity dispersion, all of which can be directly applied in observational tests of current scenarios of structure formation. On the mass scales of galaxies, if constant mass-to-light is assumed, the NFW surface density profile is found to fit HubbleReynolds laws well. It is also well fitted by Sérsic R 1/ m laws, for     but in a much narrower range of m and with much larger effective radii than are observed. Assuming in turn reasonable values of the effective radius, the mass density profiles imply a mass-to-light ratio that increases outwards at all radii.  相似文献   

5.
We present a new method of constraining the mass and velocity anisotropy profiles of galaxy clusters from kinematic data. The method is based on a model of the phase-space density, which allows the anisotropy to vary with radius between two asymptotic values. The characteristic scale of transition between these asymptotes is fixed and tuned to a typical anisotropy profile resulting from cosmological simulations. The model is parametrized by two values of anisotropy, at the centre of the cluster and at infinity, and two parameters of the NFW density profile, the scale radius and the scale mass. In order to test the performance of the method in reconstructing the true cluster parameters, we analyse mock kinematic data for 20 relaxed galaxy clusters generated from a cosmological simulation of the standard Λ cold dark matter model. We use Bayesian methods of inference and the analysis is carried out following the Markov Chain Monte Carlo approach. The parameters of the mass profile are reproduced quite well, but we note that the mass is typically underestimated by 15 per cent, probably due to the presence of small velocity substructures. The constraints on the anisotropy profile for a single cluster are in general barely conclusive. Although the central asymptotic value is determined accurately, the outer one is subject to significant systematic errors caused by substructures at large clustercentric distance. The anisotropy profile is much better constrained if one performs joint analysis of at least a few clusters. In this case, it is possible to reproduce the radial variation of the anisotropy over two decades in radius inside the virial sphere.  相似文献   

6.
We study the mass distribution in six nearby  ( z < 0.06)  relaxed Abell clusters of galaxies A0262, A0496, A1060, A2199, A3158 and A3558. Given the dominance of dark matter in galaxy clusters, we approximate their total density distribution by the Navarro, Frenk & White (NFW) formula characterized by virial mass and concentration. We also assume that the anisotropy of galactic orbits is reasonably well described by a constant and that galaxy distribution traces that of the total density. Using the velocity and position data for 120–420 galaxies per cluster we calculate, after removal of interlopers, the profiles of the lowest order even velocity moments, dispersion and kurtosis. We then reproduce the velocity moments by jointly fitting the moments to the solutions of the Jeans equations. Including the kurtosis in the analysis allows us to break the degeneracy between the mass distribution and anisotropy and constrain the anisotropy as well as the virial mass and concentration. The method is tested in detail on mock data extracted from the N -body simulations of dark matter haloes. We find that the best-fitting Galactic orbits are remarkably close to isotropic in most clusters. Using the fitted pairs of mass and concentration parameters for the six clusters, we conclude that the trend of decreasing concentration for higher masses found in the cosmological N -body simulations is consistent with the data. By scaling the individual cluster data by mass, we combine them to create a composite cluster with 1465 galaxies and perform a similar analysis on such sample. The estimated concentration parameter then lies in the range  1.5 < c < 14  and the anisotropy parameter in the range  −1.1 < β < 0.5  at the 95 per cent confidence level.  相似文献   

7.
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.  相似文献   

8.
We examine the properties of dark matter haloes within a rich galaxy cluster using a high-resolution simulation that captures the cosmological context of a cold dark matter universe. The mass and force resolution permit the resolution of 150 haloes with circular velocities larger than 80 km s−1 within the cluster virial radius of 2 Mpc (with Hubble constant H 0 = 50 km s−1 Mpc−1). This enables an unprecedented study of the statistical properties of a large sample of dark matter haloes evolving in a dense environment. The cumulative fraction of mass attached to these haloes varies from close to zero per cent at 200 kpc to 13 per cent at the virial radius. Even at this resolution the overmerging problem persists; haloes that pass within 100–200 kpc of the cluster centre are tidally disrupted. Additional substructure is lost at earlier epochs within the massive progenitor haloes. The median ratio of apocentric to pericentric radii is 6:1, so that the orbital distribution is close to isotropic, circular orbits are rare and radial orbits are common. The orbits of haloes are unbiased with respect to both position within the cluster and the orbits of the smooth dark matter background, and no velocity bias is detected. The tidal radii of surviving haloes are generally well-fitted using the simple analytic prediction applied to their orbital pericentres. Haloes within clusters have higher concentrations than those in the field. Within the cluster, halo density profiles can be modified by tidal forces and individual encounters with other haloes that cause significant mass loss —'galaxy harassment'. Mergers between haloes do not occur inside the cluster virial radius.  相似文献   

9.
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.  相似文献   

10.
In the construction of multimass King–Michie models of globular clusters, an approximated central energy equipartition between stars of different mass is usually imposed by scaling the velocity parameter of each mass class inversely with the stellar mass, as if the distribution function were isothermal. In this paper, this 'isothermal approximation' has been checked and its consequences on the model parameters studied by a comparison with models including central energy equipartition correctly. It is found that, under the isothermal approximation, the 'temperatures' of a pair of components can differ to a non-negligible amount for low concentration distributions. It is also found that, in general, this approximation leads to a significantly reduced mass segregation in comparison with that given under the exact energy equipartition at the centre. As a representative example, an isotropic three-component model fitting a given projected surface brightness and line-of-sight velocity dispersion profiles is discussed. In this example, the isothermal approximation gives a cluster envelope much more concentrated (central dimensionless potential   W = 3.3  ) than under the true equipartition  ( W = 5.9 × 10−2)  , as well as a higher mass function logarithmic slope. As a consequence, the inferred total mass (and then the global mass-to-light ratio) is a factor of 1.4 times lower than the correct value and the amount of mass in heavy dark remnants is 3.3 times smaller. Under energy equipartition, the fate of stars having a mass below a certain limit is to escape from the system. This limit is derived as a function of the mass and W of the component of giant and turn-off stars.  相似文献   

11.
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 0=3). Models with larger W 0 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.  相似文献   

12.
It is shown that the cuspy density distributions observed in the cores of elliptical galaxies can be realized by dissipationless gravitational collapse. The initial models consist of power-law density spheres such as ρ ∝ r −1 with anisotropic velocity dispersions. Collapse simulations are carried out by integrating the collisionless Boltzmann equation directly, on the assumption of spherical symmetry. From the results obtained, the extent of constant density cores, formed through violent relaxation, decreases as the velocity anisotropy increases radially, and practically disappears for extremely radially anisotropic models. As a result, the relaxed density distributions become more cuspy with increasing radial velocity anisotropy. It is thus concluded that the velocity anisotropy could be a key ingredient for the formation of density cusps in a dissipationless collapse picture. The velocity dispersions increase with radius in the cores according to the nearly power-law density distributions. The power-law index, n , of the density profiles, defined as ρ ∝ r − n , changes from n ≈2.1 at intermediate radii to a shallower power than n ≈2.1 toward the centre. This density bend can be explained from our postulated local phase-space constraint that the phase-space density accessible to the relaxed state is determined at each radius by the maximum phase-space density of the initial state.  相似文献   

13.
We extend the formalism for the calculation of the relativistic corrections to the Sunyaev–Zel'dovich effect for clusters of galaxies and include the multiple scattering effects in the isotropic approximation. We present the results of the calculations by the Fokker–Planck expansion method as well as by the direct numerical integration of the collision term of the Boltzmann equation. The multiple scattering contribution is found to be very small compared with the single scattering contribution. For high-temperature galaxy clusters of     the ratio of both the contributions is −0.2 per cent in the Wien region. In the Rayleigh–Jeans region the ratio is −0.03 per cent. Therefore the multiple scattering contribution is safely neglected for the observed galaxy clusters.  相似文献   

14.
We study the distribution of dark matter in dwarf spheroidal galaxies by modelling the moments of their line-of-sight velocity distributions. We discuss different dark matter density profiles, both cuspy and possessing flat density cores. The predictions are made in the framework of standard dynamical theory of two-component (stars and dark matter) spherical systems with different velocity distributions. We compare the predicted velocity dispersion profiles to observations in the case of Fornax and Draco dwarfs. For isotropic models the dark haloes with cores are found to fit the data better than those with cusps. Anisotropic models are studied by fitting two parameters, dark mass and velocity anisotropy, to the data. In this case all profiles yield good fits, but the steeper the cusp of the profile, the more tangential is the velocity distribution required to fit the data. To resolve this well-known degeneracy of density profile versus velocity anisotropy, we obtain predictions for the kurtosis of the line-of-sight velocity distribution for models found to provide best fits to the velocity dispersion profiles. It turns out that profiles with cores typically yield higher values of kurtosis which decrease more steeply with distance than the cuspy profiles, which will allow us to discriminate between the profiles once the kurtosis measurements become available. We also show that with present quality of the data the alternative explanation of velocity dispersions in terms of Modified Newtonian Dynamics cannot yet be ruled out.  相似文献   

15.
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.  相似文献   

16.
We study the effect of contamination by interlopers in kinematic samples of galaxy clusters. We demonstrate that without the proper removal of interlopers the inferred parameters of the mass distribution in the cluster are strongly biased towards higher mass and lower concentration. The interlopers are removed using two procedures previously shown to work most efficiently on simulated data. One is based on using the virial mass estimator and calculating the maximum velocity available to cluster members and the other relies on the ratio of the virial and projected mass estimators. We illustrate the performance of the methods in detail using the example of A576, a cluster with a strong uniform background contamination, and compare the case of A576 to 15 other clusters with different degree of contamination. We model the velocity dispersion and kurtosis profiles obtained for the cleaned data samples of these clusters solving the Jeans equations to estimate the mass, concentration and anisotropy parameter. We present the mass–concentration relation for the total sample of 22 clusters.  相似文献   

17.
A membership catalogue for Praesepe was compiled and split into four mass bins. A contour plot indicates the presence of a subcluster some 3 pc from the centre of the cluster, of approximately 30 M⊙. A tidally truncated King profile was fitted to the remainder of the cluster and the tidal radius is found to be 12.1 pc; the mass of the cluster (excluding the subcluster) is 630 M⊙. From the calculated velocity dispersions we find that the cluster appears to have too much kinetic energy and should be rapidly disintegrating. X-ray data suggest that there may be an age spread between the main core stars and the subcluster stars. This leads us to the conclusion that Praesepe is two merging clusters.  相似文献   

18.
We use oblate axisymmetric dynamical models including dark haloes to determine the orbital structure of intermediate mass to massive early-type galaxies in the Coma galaxy cluster. We find a large variety of orbital compositions. Averaged over all sample galaxies the unordered stellar kinetic energy in the azimuthal and the radial direction are of the same order, but they can differ by up to 40 per cent in individual systems. In contrast, both for rotating and non-rotating galaxies the vertical kinetic energy is on average smaller than in the other two directions. This implies that even most of the rotating ellipticals are flattened by an anisotropy in the stellar velocity dispersions. Using three-integral axisymmetric toy models, we show that flattening by stellar anisotropy maximizes the entropy for a given density distribution. Collisionless disc merger remnants are radially anisotropic. The apparent lack of strong radial anisotropy in observed early-type galaxies implies that they may not have formed from mergers of discs unless the influence of dissipational processes was significant.  相似文献   

19.
A revision of Stodoíkiewicz's Monte Carlo code is used to simulate evolution of star clusters. The new method treats each superstar as a single star and follows the evolution and motion of all individual stellar objects. The first calculations for isolated, equal-mass N -body systems with three-body energy generation according to Spitzer's formulae show good agreement with direct N -body calculations for N  = 2000, 4096 and 10 000 particles. The density, velocity, mass distributions, energy generation, number of binaries, etc., follow the N -body results. Only the number of escapers is slightly too high compared with N -body results, and there is no level-off anisotropy for advanced post-collapse evolution of Monte Carlo models as is seen in N -body simulations for N  ≤ 2000. For simulations with N  > 10 000 gravothermal oscillations are clearly visible. The calculations of N   2000, 4096, 10 000, 32 000 and 100 000 models take about 2, 6, 20, 130 and 2500 h, respectively. The Monte Carlo code is at least 105 times faster than the N -body one for N  = 32 768 with special-purpose hardware. Thus it becomes possible to run several different models to improve statistical quality of the data and run individual models with N as large as 100 000. The Monte Carlo scheme can be regarded as a method which lies in the middle between direct N -body and Fokker–Planck models and combines most advantages of both methods.  相似文献   

20.
A membership catalogue for the Pleiades is divided into four mass bins, and a tidally truncated King profile is fitted to each bin with good agreement with the data. The tidal radius of the cluster is found to be 13.1 pc, and the total mass of the cluster down to the stellar limit is calculated to be 735 M. The spread of stars in each bin, as well as the relaxation and crossing times, shows the Pleiades to be an approximately relaxed cluster with equilibrium between the density and velocity distributions. The cluster kinetic energy and binding energy are consistent with the virial theorem, indicating no large unseen population of brown dwarfs. However, the 1 σ errors in the cluster parameters provide an upper limit to the mass of any brown dwarf population of 131 M, which would show up in deep CCD surveys as ≤ 5.5 brown dwarfs per 10 × 10 arcminute field in the cluster centre.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号