共查询到20条相似文献，搜索用时 51 毫秒
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. 相似文献
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 study the evolution of globular clusters with mass spectra under the influence of the steady Galactic tidal field, including the effects of velocity anisotropy. Similarly to single-mass models, velocity anisotropy develops as the cluster evolves, but the degree of anisotropy is much smaller than in isolated clusters. Except for very early epochs of the cluster evolution, the velocity distributions of nearly all mass components become tangentially anisotropic at the outer parts. We examine how the mass function (MF) changes in time. Specifically, we find that the power-law index of the MF decreases monotonically with the total mass of the cluster, in agreement with previous findings based on isotropic models or N -body studies. This is also consistent with the behaviour of the observed slopes of MFs for a limited number of clusters. We attempt to compare our results with multimass King models, although it is almost impossible to fit the entire density profiles for all mass components. When the MF is fixed, the central densities of individual components show significant differences between Fokker–Planck and King models. We obtain 'best-fitting' multimass King models, for which the central density of individual components as well as the total density distribution agrees with the Fokker–Planck models by adjusting the MF. The MFs obtained in this way closely resemble the MF within the half-mass radius of the Fokker–Planck result. Also, we find that the local MFs predicted by Fokker–Planck calculations vary more rapidly with radius than best-fitting multimass King models. The projected velocity profiles for anisotropic models show significant flattening toward the tidal radius compared with the isotropic model. This is caused by the fact that the tangential velocity dispersion becomes dominant at the outer parts. Such a behaviour of velocity profile appears to be consistent with the observed profiles of the collapsed cluster M15. 相似文献
3.
Mirek Giersz 《Monthly notices of the Royal Astronomical Society》2001,324(1):218-230
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. 相似文献
4.
Mirek Giersz 《Monthly notices of the Royal Astronomical Society》1998,298(4):1239-1248
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. 相似文献
5.
Naoki Itoh Youhei Kawana Satoshi Nozawa Yasuharu Kohyama 《Monthly notices of the Royal Astronomical Society》2001,327(2):567-576
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. 相似文献
6.
M. Di Criscienzo P. Ventura F. D'Antona M. Marconi A. Ruoppo V. Ripepi 《Monthly notices of the Royal Astronomical Society》2008,389(1):325-332
We applied the aton evolutionary code to the computation of detailed grids of standard (non-rotating) and rotating pre-main sequence (PMS) models and computed their adiabatic oscillation spectra, with the aim of exploring the seismic properties of young stars. As, until now, only a few frequencies have been determined for ∼40 PMS stars, the way of approaching the interpretation of the oscillations is not unique. We adopt a method similar to the matching mode method by Guenther and Brown making use, when necessary, also of our rotating evolutionary code to compute the models for PMS stars. The method is described by a preliminary application to the frequency spectrum of two PMS stars (85 and 278) in the young open cluster NGC 6530. For the Star 85, we confirm with self-consistent rotating models, previous interpretation of the data, attributing three close frequencies to the mode n = 4, l = 1 and m = 0 , +1 and −1. For the Star 278, we find a different fit for the frequencies, corresponding to a model within the original error box of the star, and dispute the possibility that this star has a T eff much cooler that the red boundary of the radial instability strip. 相似文献
7.
Douglas C. Heggie Michele Trenti Piet Hut 《Monthly notices of the Royal Astronomical Society》2006,368(2):677-689
In order to interpret the results of complex realistic star cluster simulations, which rely on many simplifying approximations and assumptions, it is essential to study the behaviour of even more idealized models, which can highlight the essential physical effects and are amenable to more exact methods. With this aim, we present the results of N -body calculations of the evolution of equal-mass models, starting with primordial binary fractions of 0–100 per cent, with values of N ranging from 256 to 16 384. This allows us to extrapolate the main features of the evolution to systems comparable in particle number with globular clusters.
In this range, we find that the steady-state 'deuterium main sequence' is characterized by a ratio of the core radius to half-mass radius that follows qualitatively the analytical estimate by Vesperini & Chernoff, although the N dependence is steeper than expected. Interestingly, for an initial binary fraction f greater than 10 per cent, the binary heating in the core during the post-collapse phase almost saturates (becoming nearly independent of f ), and so little variation in the structural properties is observed. Thus, although we observe a significantly lower binary abundance in the core with respect to the Fokker–Planck simulations by Gao et al., this is of little dynamical consequence.
At variance with the study of Gao et al., we see no sign of gravothermal oscillations before 150 half-mass relaxation times. At later times, however, oscillations become prominent. We demonstrate the gravothermal nature of these oscillations. 相似文献
In this range, we find that the steady-state 'deuterium main sequence' is characterized by a ratio of the core radius to half-mass radius that follows qualitatively the analytical estimate by Vesperini & Chernoff, although the N dependence is steeper than expected. Interestingly, for an initial binary fraction f greater than 10 per cent, the binary heating in the core during the post-collapse phase almost saturates (becoming nearly independent of f ), and so little variation in the structural properties is observed. Thus, although we observe a significantly lower binary abundance in the core with respect to the Fokker–Planck simulations by Gao et al., this is of little dynamical consequence.
At variance with the study of Gao et al., we see no sign of gravothermal oscillations before 150 half-mass relaxation times. At later times, however, oscillations become prominent. We demonstrate the gravothermal nature of these oscillations. 相似文献
8.
Rainer Spurzem Holger Baumgardt Peter Berczik Christian Boily Stefan Deiters Pavel Kroupa Panos Patsis 《Astronomische Nachrichten》2005,326(7):589-605
G01 New evidence for a connection between massive black holes and ULX G02 Long‐Term Evolution of Massive Black Hole Binaries G03 NBODY Meets Stellar Population Synthesis G04 N‐body modelling of real globular star clusters G05 Fokker‐Planck rotating models of globular clusters with black hole G06 Observational Manifestation of chaos in spiral galaxies: quantitative analysis and qualitative explanation G07 GRAPE Clusters: Beyond the Million‐Body Problem G08 Orbital decay of star clusters and Massive Black Holes in cuspy galactic nuclei G09 An Edge‐on Disk Galaxy Catalog G10 Complexes of open clusters in the Solar neighborhood G11 Search for and investigation of new stellar clusters using the data from huge stellar catalogues G12 Computing 2D images of 3D galactic disk models G13 Outer Pseudoring in the Galaxy G14 Where are tidal‐dwarf galaxies? G15 Ultra compact dwarf galaxies in nearby clusters G16 Impact of an Accretion Disk on the Structure of a stellar cluster in active galactic nuclei G17 Order and Chaos in the edge‐on profiles of disk galaxies G18 On the stability of OB‐star configurations in the Orion Nebula cluster G19 Older stars captured in young star clusters by cloud collapse G20 General features of the population of open clusters within 1 kpc from the Sun G21 Unstable modes in thin stellar disks G22 From Newton to Einstein – Dynamics of N‐body systems G23 On the relation between the maximum stellar mass and the star cluster mass 相似文献
9.
Gravitational wave emission by coalescing black holes (BHs) kicks the remnant BH with a typical velocity of hundreds of km s−1 . This velocity is sufficiently large to remove the remnant BH from a low-mass galaxy but is below the escape velocity from the Milky Way (MW) galaxy. If central BHs were common in the galactic building blocks that merged to make the MW, then numerous BHs that were kicked out of low-mass galaxies should be freely floating in the MW halo today. We use a large statistical sample of possible merger tree histories for the MW to estimate the expected number of recoiled BH remnants present in the MW halo today. We find that hundreds of BHs should remain bound to the MW halo after leaving their parent low-mass galaxies. Each BH carries a compact cluster of old stars that populated the core of its original host galaxy. Using the time-dependent Fokker–Planck equation, we find that the present-day clusters are ≲1 pc in size, and their central bright regions should be unresolved in most existing sky surveys. These compact systems are distinguishable from globular clusters by their internal (Keplerian) velocity dispersion greater than 100 km s−1 and their high mass-to-light ratio owing to the central BH. An observational discovery of this relic population of star clusters in the MW halo would constrain the formation history of the MW and the dynamics of BH mergers in the early Universe. A similar population should exist around other galaxies and may potentially be detectable in M31 and M33. 相似文献
10.
Dong Lai 《Monthly notices of the Royal Astronomical Society》1999,307(4):1001-1007
Gravitational radiation tends to drive gravity modes in rotating neutron stars to become unstable. For an inviscid star, the instability sets in when the rotation frequency is about 0.7 times the corresponding mode frequency of the non-rotating star. Neutron stars with spin frequencies ≳100 Hz are susceptible to this instability, with a growth time of the order of years. However, it is likely that viscous dissipation suppresses the instability except for a narrow range of temperatures around 109 K. We also show that the viscosity-driven instability of g-modes is absent. 相似文献
11.
Globular clusters rotate significantly, and with the increasing amount of detailed morphological and kinematical data obtained in recent years on galactic globular clusters many interesting features show up. We show how our theoretical evolutionary models of rotating clusters can be used to obtain fits, which at least properly model the overall rotation and its implied kinematics in full 2D detail (dispersions, rotation velocities). Our simplified equal mass axisymmetric rotating model provides detailed two-dimensional kinematical and morphological data for star clusters. The degree of rotation is not dominant in energy, but also non-negligible for the phase-space distribution function, shape and kinematics of clusters. Therefore, the models are well applicable for galactic globular clusters. Since previously published papers on that matter by us made it difficult to do detailed comparisons with observations, we provide a much more comprehensive and easy-to-use set of data here, which uses as entries dynamical age and flattening of observed cluster and then offers a limited range of applicable models in full detail. The method, data structure and some exemplary comparison with observations are presented. Future work will improve modelling and data base to take a central black hole, a mass spectrum and stellar evolution into account. 相似文献
12.
N. Babkovskaia A. Brandenburg J. Poutanen 《Monthly notices of the Royal Astronomical Society》2008,386(2):1038-1044
In an attempt to model the accretion on to a neutron star in low-mass X-ray binaries, we present 2D hydrodynamical models of the gas flow in close vicinity of the stellar surface. First, we consider a gas pressure-dominated case, assuming that the star is non-rotating. For the stellar mass we take M star = 1.4 × 10−2 M⊙ and for the gas temperature T = 5 × 106 K . Our results are qualitatively different in the case of a realistic neutron star mass and a realistic gas temperature of T ≃ 108 K, when the radiation pressure dominates. We show that to get the stationary solution in a latter case, the star most probably has to rotate with the considerable velocity. 相似文献
13.
Holger Baumgardt Douglas C. Heggie Piet Hut Junichiro Makino 《Monthly notices of the Royal Astronomical Society》2003,341(1):247-250
This paper considers the phenomenon of deep core collapse in collisional stellar systems, with stars of equal mass. The collapse takes place on some multiple, ξ−1 , of the central relaxation time, and produces a density profile in which ρ∝ r −α , where α is a constant. The parameters α and ξ have usually been determined from simplified models, such as gas and Fokker–Planck models, often with the simplification of isotropy. Here we determine the parameters directly from N -body simulations carried out using the newly completed GRAPE-6. 相似文献
14.
Yacine Ali-Haïmoud Christopher M. Hirata Clive Dickinson 《Monthly notices of the Royal Astronomical Society》2009,395(2):1055-1078
We present a comprehensive treatment of the spectrum of electric dipole emission from spinning dust grains, updating the commonly used model of Draine & Lazarian. Grain angular velocity distributions are computed using the Fokker–Planck equation; we revisit the drift and diffusion coefficients for the major torques on the grain, including collisions, grain-plasma interactions and infrared emission. We use updated grain optical properties and size distributions. The theoretical formalism is implemented in the companion code, spdust , which is publicly available. The effect of some environmental and grain parameters on the emissivity is shown and analysed. 相似文献
15.
《天文和天体物理学研究(英文版)》2017,(7)
We present a new tool for color-magnitude diagram(CMD) studies, Powerful CMD. This tool is built based on the advanced stellar population synthesis(ASPS) model, in which single stars, binary stars, rotating stars and star formation history have been taken into account. Via Powerful CMD, the distance modulus, color excess, metallicity, age, binary fraction, rotating star fraction and star formation history of star clusters can be determined simultaneously from observed CMDs. The new tool is tested via both simulated and real star clusters. Five parameters of clusters NGC 6362, NGC 6652, NGC 6838 and M67 are determined and compared to other works. It is shown that this tool is useful for CMD studies, in particular for those utilizing data from the Hubble Space Telescope(HST). Moreover, we find that inclusion of binaries in theoretical stellar population models may lead to smaller color excess compared to the case of single-star population models. 相似文献
16.
Non-linear hydrodynamical evolution of rotating relativistic stars: numerical methods and code tests
José A. Font † Nikolaos Stergioulas † Kostas D. Kokkotas † 《Monthly notices of the Royal Astronomical Society》2000,313(4):678-688
We present numerical hydrodynamical evolutions of rapidly rotating relativistic stars, using an axisymmetric, non-linear relativistic hydrodynamics code. We use four different high-resolution shock-capturing (HRSC) finite-difference schemes (based on approximate Riemann solvers) and compare their accuracy in preserving uniformly rotating stationary initial configurations in long-term evolutions. Among these four schemes, we find that the third-order piecewise parabolic method scheme is superior in maintaining the initial rotation law in long-term evolutions, especially near the surface of the star. It is further shown that HRSC schemes are suitable for the evolution of perturbed neutron stars and for the accurate identification (via Fourier transforms) of normal modes of oscillation. This is demonstrated for radial and quadrupolar pulsations in the non-rotating limit, where we find good agreement with frequencies obtained with a linear perturbation code. The code can be used for studying small-amplitude or non-linear pulsations of differentially rotating neutron stars, while our present results serve as testbed computations for three-dimensional general-relativistic evolution codes. 相似文献
17.
We discuss the nature of the various modes of pulsation of superfluid neutron stars using comparatively simple Newtonian models and the Cowling approximation. The matter in these stars is described in terms of a two-fluid model, where one fluid is the neutron superfluid, which is believed to exist in the core and inner crust of mature neutron stars, and the other fluid represents a conglomerate of all other constituents (crust nuclei, protons, electrons, etc.). In our model, we incorporate the non-dissipative interaction known as the entrainment effect, whereby the momentum of one constituent (e.g. the neutrons) carries along part of the mass of the other constituent. We show that there is no independent set of pulsating g-modes in a non-rotating superfluid neutron star core, even though the linearized superfluid equations contain a well-defined (and real-valued) analogue to the so-called Brunt–Väisälä frequency. Instead, what we find are two sets of spheroidal perturbations whose nature is predominately acoustic. In addition, an analysis of the zero-frequency subspace (i.e. the space of time-independent perturbations) reveals two sets of degenerate spheroidal perturbations, which we interpret to be the missing g-modes, and two sets of toroidal perturbations. We anticipate that the degeneracy of all these zero-frequency modes will be broken by the Coriolis force in the case of rotating stars. To illustrate this we consider the toroidal pulsation modes of a slowly rotating superfluid star. This analysis shows that the superfluid equations support a new class of r-modes, in addition to those familiar from, for example, geophysical fluid dynamics. Finally, the role of the entrainment effect on the superfluid mode frequencies is shown explicitly via solutions to dispersion relations that follow from a 'local' analysis of the linearized superfluid equations. 相似文献
18.
The models of non-rotating and rotating 2.31M
\ stars of Population I composition have been calculated, starting at the threshold of stability. A 2.31M
\ star was chosen to compare the results with the observational parameters of the primary component of the well-known detached binary YZ Cassiopeiae. The effects of rotation on the internal structure during the evolution of the star were studied by constructing sequences of axisymmetric rotating models under the assumption that angular momentum was conserved according to a predetermined angular velocity distribution depending on the structure of the star.The first section of this paper deals with effects of rotation on the evolutionary behaviours of the 2.31M
\ star through the pre-Main-Sequence evolution as well as the zero-age Main Sequence.In the second section of this paper, the evolutionary studies have been extended up to near-hydrogen exhaustion phase in order to obtain a theoretical model corresponding to the given mass and radius of the primary component of YZ Cassiopeiae. The theoretical models were found to be in a good agreement with observational parameters. The computed rotating models of the primary of YZ Cassiopeiae indicates that its evolutionary age is 6.01×108 years; and the central hydrogen content 0.183 — which means that about 75% of its original value was depleted. 相似文献
19.
K. S. Wong 《Monthly notices of the Royal Astronomical Society》2009,396(1):439-461
In this paper, we shall study the problem of stability of polar oscillations of a non-rotating Newtonian superfluid star. We find that the stability of the latter system is guaranteed by the positive definiteness of a matrix quantity. We also find conditions which characterize the occurrence of zero-frequency polar modes, which marks the onset of polar instability. We find that the negative definiteness of the mentioned matrix quantity implies instability of the system. We apply our results to show that the polar oscillations of a non-rotating Newtonian superfluid star in the zero-temperature approximation are marginally stable. 相似文献
20.
A. Passamonti B. Haskell N. Andersson 《Monthly notices of the Royal Astronomical Society》2009,396(2):951-963
Using time evolutions of the relevant linearized equations, we study non-axisymmetric oscillations of rapidly rotating and superfluid neutron stars. We consider perturbations of Newtonian axisymmetric background configurations and account for the presence of superfluid components via the standard two-fluid model. Within the Cowling approximation, we are able to carry out evolutions for uniformly rotating stars up to the mass-shedding limit. This leads to the first detailed analysis of superfluid neutron star oscillations in the fast rotation regime, where the star is significantly deformed by the centrifugal force. For simplicity, we focus on background models where the two fluids (superfluid neutrons and protons) corotate, are in β-equilibrium and co-exist throughout the volume of the star. We construct sequences of rotating stars for two analytical model equations of state. These models represent relatively simple generalizations of single fluid, polytropic stars. We study the effects of entrainment, rotation and symmetry energy on non-radial oscillations of these models. Our results show that entrainment and symmetry energy can have a significant effect on the rotational splitting of non-axisymmetric modes. In particular, the symmetry energy modifies the inertial mode frequencies considerably in the regime of fast rotation. 相似文献