共查询到20条相似文献，搜索用时 15 毫秒
1.
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. 相似文献
2.
In our previous papers, we showed that at the final phases of the dynamical evolution of an open cluster, an extended population of stars elongated along its Galactic orbit, the stellar tail of the cluster, is formed. The tail stars that escaped from the cluster at different times move in a common orbit with low relative velocities. Experiencing a weak interaction with Galactic field stars, these objects, the relics of open clusters, can exist for a fairly long time. In this paper, we investigate the structures of such stellar tails in the nearest open clusters: Hyades, Pleiades, Praesepe, Alpha Persei, Coma, IC 2391, and IC 2602. To this end, we performed several numerical simulations of the dynamical evolution of these clusters in the tidal field of the Galaxy. Our computations of the dynamical evolution were based on known cluster age estimates and real Galactic orbits. The initial conditions were chosen in such a way that the parameters of the simulated clusters corresponded to their observed parameters. As a result, we obtained models of the stellar tails for the nearest open clusters and estimated such parameters of the tails as their sizes, densities, locations relative to the solar neighborhood, and others. 相似文献
3.
对于星系际弥散恒星的研究是分别从观测、数值模拟和半解析模型这三个方面进行的.现在已经在邻近星系团及中低红移处观测到弥散恒星,甚至在Virgo及Coma星系团中观测到了单个的弥散恒星.观测数据的积累使得人们能够从统计上了解星系际弥散恒星的性质.研究表明星系际弥散恒星围绕着星系团势阱中心呈椭球状对称分布,其在星系团恒星总质... 相似文献
4.
Sverre J. Aarseth 《Celestial Mechanics and Dynamical Astronomy》1999,73(1-4):127-137
This paper concentrates on four key tools for performing star cluster simulations developed during the last decade which are
sufficient to handle all the relevant dynamical aspects. First we discuss briefly the Hermite integration scheme which is
simple to use and highly efficient for advancing the single particles. The main numerical challenge is in dealing with weakly
and strongly perturbed hard binaries. A new treatment of the classical Kustaanheimo-Stiefel two-body regularization has proved
to be more accurate for studying binaries than previous algorithms based on divided differences or Hermite integration. This
formulation employs a Taylor series expansion combined with the Stumpff functions, still with one force evaluation per step,
which gives exact solutions for unperturbed motion and is at least comparable to the polynomial methods for large perturbations.
Strong interactions between hard binaries and single stars or other binaries are studied by chain regularization which ensures
a non-biased outcome for chaotic motions. A new semi-analytical stability criterion for hierarchical systems has been adopted
and the long-term effects on the inner binary are now treated by averaging techniques for cases of interest. These modifications
describe consistent changes of the orbital variables due to large Kozai cycles and tidal dissipation. The range of astrophysical
processes which can now be considered by N-body simulations include tidal capture, circularization, mass transfer by Roche-lobe
overflow as well as physical collisions, where the masses and radii of individual stars are modelled by synthetic stellar
evolution.
This revised version was published online in July 2006 with corrections to the Cover Date. 相似文献
5.
The luminous barred galaxy NGC 6872 is one of the largest spiral galaxies known. Star formation is traced all along the arms,
which extend over more than 100 kpc. The galaxy is likely to be affected by tidal perturbations from the nearby companion
IC 4970 passing by on a low-inclination, prograde orbit. We have mapped the large-scale distribution and kinematics of the
atomic gas in the NGC 6872/IC 4970 system and carried out N-body simulations. Our aim is to investigate the evolution of gas
and stars in such close encounters, examine the influence of a dark matter halo on the length of the tidal tails, and test
models of collisionally induced star formation.
This revised version was published online in August 2006 with corrections to the Cover Date. 相似文献
6.
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. 相似文献
7.
M. Hénon 《Astrophysics and Space Science》1971,13(2):284-299
8.
Kazuhide Akiyama 《Earth, Moon, and Planets》1991,54(3):203-239
Dynamical evolution of globular clusters in the Large Magellanic Cloud (LMC) is investigated by means of N-body simulations; particular attention is paid to time evolution in the ellipticitical figure of globular clusters. The simulations were started with a binary globular cluster. It merged into a single cluster with ellipticity of about 0.3. The simulations were continued until the cluster became rounder due to the effects of two body relaxation and of tidal field of LMC. It is found that the outward angular momentum transport due to the gravothermal contraction makes the inner region rounder; the ellipticity at about the initial half-mass radius (r
h) decreases with the e-folding time of 20 relaxation times. On the other hand, the outer region becomes rounder due to the stripping of stars by the tidal field; the ellipticity at about 3r
h decreases with the e-folding time of 80 crossing times therein, though the time scale depends on the direction of the tidal field relative to the spin of the cluster. These two effects are comparable at about the half-mass radius. Taking account of such theoretical results we reanalyzed observed data for the ellipticity at about the half-mass radius of LMC clusters. We estimated the relaxation time and crossing time for each of the observed clusters, from which we calculated the effective time of getting round of the cluster. We plotted the observed ellipticity of the clusters against their non-dimensional age — i.e., the age normalized by the effective time. We found that observed ellipticity distribution is consistent with our picture. 相似文献
9.
Jarrod R. Hurley † Christopher A. Tout Sverre J. Aarseth Onno R. Pols 《Monthly notices of the Royal Astronomical Society》2001,323(3):630-650
We present a state-of-the-art N -body code which includes a detailed treatment of stellar and binary evolution as well as the cluster dynamics. This code is ideal for investigating all aspects relating to the evolution of star clusters and their stellar populations. It is applicable to open and globular clusters of any age. We use the N -body code to model the blue straggler population of the old open cluster M67. Preliminary calculations with our binary population synthesis code show that binary evolution alone cannot explain the observed numbers or properties of the blue stragglers. On the other hand, our N -body model of M67 generates the required number of blue stragglers and provides formation paths for all the various types found in M67. This demonstrates the effectiveness of the cluster environment in modifying the nature of the stars it contains, and highlights the importance of combining dynamics with stellar evolution. We also perform a series of N =10 000 simulations in order to quantify the rate of escape of stars from a cluster subject to the Galactic tidal field. 相似文献
10.
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. 相似文献
11.
Kenji Bekki 《Astrophysics and Space Science》2001,276(2-4):1033-1036
A significant fraction of clusters of galaxies are observed to have substructure, which implies that merging between clusters
and subclusters is a rather common physical process in cluster formation. It still remains unclear how cluster merging affects
the evolution of cluster member galaxies. We report the results of numerical simulations that show the dynamical evolution
of a gas-rich, late-type spiral in a merger between a small group of galaxies and a cluster. The simulations demonstrate that
the time-dependent tidal gravitational field during merging excites non-axisymmetric structure of the galaxy, subsequently
drives efficient transfer of gas to the central region and finally triggers a secondary starburst. This result provides a
close physical relationship between the emergence of starburst galaxies and the formation of substructure in clusters. We
accordingly interpret post-starburst galaxies located near substructure of the Coma Cluster as one observational example indicating
the global tidal effects of group–cluster merging. Our numerical results further suggest a causal link between the observed
excess of blue galaxies in distant clusters and the cluster virialization process through hierarchical merging of subclusters.
This revised version was published online in July 2006 with corrections to the Cover Date. 相似文献
12.
Nurur Rahman †‡ Janusz Krywult Patrick M. Motl Piotr Flin Sergei F. Shandarin 《Monthly notices of the Royal Astronomical Society》2006,367(2):838-850
We have made a comparative study of morphological evolution in simulated dark matter (DM) haloes and X-ray brightness distribution, and in optical clusters. Samples of simulated clusters include star formation with supernovae feedback, radiative cooling and simulation in the adiabatic limit at three different redshifts, z = 0.0, 0.10 and 0.25. The optical sample contains 208 Abell, Corwin & Olowin (ACO) clusters within redshift, z ≤ 0.25 . Cluster morphology, within 0.5 and 1.0 h −1 Mpc from cluster centre, is quantified by multiplicity and ellipticity.
We find that the distribution of the DM haloes in the adiabatic simulation appears to be more elongated than the galaxy clusters. Radiative cooling brings halo shapes in excellent agreement with observed clusters; however, cooling along with feedback mechanism makes the haloes more flattened.
Our results indicate relatively stronger structural evolution and more clumpy distributions in observed clusters than in the structure of simulated clusters, and slower increase in simulated cluster shapes compared to those in the observed one.
Within z ≤ 0.1 , we note an interesting agreement in the shapes of clusters obtained from the cooling simulations and observation. We also note that the different samples of observed clusters differ significantly in morphological evolution with redshift. We highlight a few possibilities responsible for the discrepancy in morphological evolution of simulated and observed clusters. 相似文献
We find that the distribution of the DM haloes in the adiabatic simulation appears to be more elongated than the galaxy clusters. Radiative cooling brings halo shapes in excellent agreement with observed clusters; however, cooling along with feedback mechanism makes the haloes more flattened.
Our results indicate relatively stronger structural evolution and more clumpy distributions in observed clusters than in the structure of simulated clusters, and slower increase in simulated cluster shapes compared to those in the observed one.
Within z ≤ 0.1 , we note an interesting agreement in the shapes of clusters obtained from the cooling simulations and observation. We also note that the different samples of observed clusters differ significantly in morphological evolution with redshift. We highlight a few possibilities responsible for the discrepancy in morphological evolution of simulated and observed clusters. 相似文献
13.
We performed 2D numerical simulations of oscillatory tidal flexing to study the interrelationship between tidal dissipation (calculated using the Maxwell model) and a heterogeneous temperature structure in Europa’s ice shell. Our 2D simulations show that, if the temperature is spatially uniform, the tidal dissipation rate peaks when the Maxwell time is close to the tidal period, consistent with previous studies. The tidal dissipation rate in a convective plume encased in a different background temperature depends on both the plume and background temperature. At a fixed background temperature, the dissipation increases strongly with plume temperature at low temperatures, peaks, and then decreases with temperature near the melting point when a melting-temperature viscosity of 1013 Pa s is used; however, the peak occurs at significantly higher temperature in this heterogeneous case than in a homogeneous medium for equivalent rheology. For constant plume temperature, the dissipation rate in a plume decreases as the surrounding temperature increases; plumes that are warmer than their surroundings can exhibit enhanced heating not only relative to their surroundings but relative to the Maxwell-model prediction for a homogeneous medium at the plume temperature. These results have important implications for thermal feedbacks in Europa’s ice shell.To self-consistently determine how convection interacts with tidal heating that is correctly calculated from the time-evolving heterogeneous temperature field, we coupled viscoelastic simulations of oscillatory tidal flexing (using Tekton) to long-term simulations of the convective evolution (using ConMan). Our simulations show that the tidal dissipation rate resulting from heterogeneous temperature can have a strong impact on thermal convection in Europa’s ice shell. Temperatures within upwelling plumes are greatly enhanced and can reach the melting temperature under plausible tidal-flexing amplitude for Europa. A pre-existing fracture zone (at least 6 km deep) promotes the concentration of tidal dissipation (up to ∼20 times more than that in the surroundings), leading to lithospheric thinning. This supports the idea that spatially variable tidal dissipation could lead locally to high temperatures, partial melting, and play an important role in the formation of ridges, chaos, or other features. 相似文献
14.
In a previous work (Callegari and Yokoyama, Celest. Mech. Dyn. Astr. 98:5–30, 2007), the main features of the motion of the
pair Enceladus–Dione were analyzed in the frozen regime, i.e., without considering the tidal evolution. Here, the results
of a great deal of numerical simulations of a pair of satellites similar to Enceladus and Dione crossing the 2:1 mean-motion
resonance are shown. The resonance crossing is modeled with a linear tidal theory, considering a two-degrees-of-freedom model
written in the framework of the general three-body planar problem. The main regimes of motion of the system during the passage
through resonance are studied in detail. We discuss our results comparing them with classical scenarios of tidal evolution
of the system. We show new scenarios of evolution of the Enceladus–Dione system through resonance not shown in previous approaches
of the problem. 相似文献
15.
Julia M. Comerford Priyamvada Natarajan 《Monthly notices of the Royal Astronomical Society》2007,379(1):190-200
The properties of clusters of galaxies offer key insights into the assembly process of structure in the universe. Numerical simulations of cosmic structure formation in a hierarchical, dark matter dominated universe suggest that galaxy cluster concentrations, which are a measure of a halo's central density, decrease gradually with virial mass. However, cluster observations have yet to confirm this correlation. The slopes of the run of measured concentrations with virial mass are often either steeper or flatter than that predicted by simulations. In this work, we present the most complete sample of observed cluster concentrations and masses yet assembled, including new measurements for 10 strong-lensing clusters, thereby more than doubling the existing number of strong-lensing concentration estimates. We fit a power law to the observed concentrations as a function of virial mass, and find that the slope is consistent with the slopes found in simulations, though our normalization factor is higher. Observed lensing concentrations appear to be systematically larger than X-ray concentrations, a more pronounced effect than that found in simulations. We also find that at a fixed mass, the bulk of observed cluster concentrations are distributed lognormally, with the exception of a few anomalously high concentration clusters. We examine the physical processes likely responsible for the discrepancy between lensing and X-ray concentrations, and for the anomalously high concentrations in particular. The forthcoming Millennium simulation results will offer the most comprehensive comparison set to our findings of an observed concentration–mass power law relation. 相似文献
16.
J. Donnert K. Dolag H. Lesch E. Müller 《Monthly notices of the Royal Astronomical Society》2009,392(3):1008-1021
We performed cosmological, magnetohydrodynamical simulations to follow the evolution of magnetic fields in galaxy clusters, exploring the possibility that the origin of the magnetic seed fields is galactic outflows during the starburst phase of galactic evolution. To do this, we coupled a semi-analytical model for magnetized galactic winds as suggested by Bertone, Vogt & Enßlin to our cosmological simulation. We find that the strength and structure of magnetic fields observed in galaxy clusters are well reproduced for a wide range of model parameters for the magnetized, galactic winds and do only weakly depend on the exact magnetic structure within the assumed galactic outflows. Although the evolution of a primordial magnetic seed field shows no significant differences to that of galaxy cluster fields from previous studies, we find that the magnetic field pollution in the diffuse medium within filaments is below the level predicted by scenarios with pure primordial magnetic seed field. We therefore conclude that magnetized galactic outflows and their subsequent evolution within the intracluster medium can fully account for the observed magnetic fields in galaxy clusters. Our findings also suggest that measuring cosmological magnetic fields in low-density environments such as filaments is much more useful than observing cluster magnetic fields to infer their possible origin. 相似文献
17.
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. 相似文献
18.
Mirek Giersz 《Monthly notices of the Royal Astronomical Society》2006,371(1):484-494
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. 相似文献
19.
Shogo Inagaki 《Astrophysics and Space Science》1986,118(1-2):367-370
N-body simulations were made to study post-collapse evolution of small-N clusters. The qualitative behaviour of the dynamical evolution is the same as in the Fokker-Planck approximation. The flatter is the mass spectrum, the more violent is the activity of binaries and so are the expansion and contraction of the cluster core.Paper presented at the IAU Third Asian-Pacific Regional Meeting, held in Kyoto, Japan, between 30 September–6 October, 1984. 相似文献