首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 28 毫秒
1.
We perform a stability test of triaxial models in Modified Newtonian Dynamics (MOND) using N -body simulations. The triaxial models considered here have densities that vary with   r −1  in the centre and   r −4  at large radii. The total mass of the model varies from 108 to  1010 M⊙  , representing the mass scale of dwarfs to medium-mass elliptical galaxies, respectively, from deep MOND to quasi-Newtonian gravity. We build triaxial galaxy models using the Schwarzschild technique, and evolve the systems for 200 Keplerian dynamical times (at the typical length-scale of 1.0 kpc). We find that the systems are virial overheating, and in quasi-equilibrium with the relaxation taking approximately 5 Keplerian dynamical times (1.0 kpc). For all systems, the change of the inertial (kinetic) energy is less than 10 per cent (20 per cent) after relaxation. However, the central profile of the model is flattened during the relaxation and the (overall) axis ratios change by roughly 10 per cent within 200 Keplerian dynamical times (at 1.0 kpc) in our simulations. We further find that the systems are stable once they reach the equilibrium state.  相似文献   

2.
We present predictions for the line-of-sight velocity dispersion profiles of dwarf spheroidal galaxies and compare them to observations in the case of the Fornax dwarf. The predictions are made in the framework of standard dynamical theory of spherical systems with different velocity distributions. The stars are assumed to be distributed according to Sérsic laws with parameters fitted to observations. We compare predictions obtained assuming the presence of dark matter haloes (with density profiles adopted from N -body simulations) with those resulting from Modified Newtonian Dynamics (MOND). If the anisotropy of velocity distribution is treated as a free parameter, observational data for Fornax are reproduced equally well by models with dark matter and with MOND. If stellar mass-to-light ratio of 1 M⊙/L⊙ is assumed, the required mass of the dark halo is     , two orders of magnitude larger than the mass in stars. The derived MOND acceleration scale is     . In both cases a certain amount of tangential anisotropy in the velocity distribution is needed to reproduce the shape of the velocity dispersion profile in Fornax.  相似文献   

3.
We calculate the structure of a wake generated by, and the dynamical friction force on, a gravitational perturber travelling through a gaseous medium of uniform density and constant background acceleration   g ext  , in the context of Modified Newtonian Dynamics (MOND). The wake is described as a linear superposition of two terms. The dominant part displays the same structure as the wake generated in the Newtonian gravity scaled up by a factor  μ−1( g ext/ a 0)  , where a 0 is the constant MOND acceleration and μ the interpolating function. The structure of the second term depends greatly on the angle between   g ext  and the velocity of the perturber. We evaluate the dynamical drag force numerically and compare our MOND results with the Newtonian case. We mention the relevance of our calculations to orbit evolution of globular clusters and satellites in a gaseous protogalaxy. Potential differences in the X-ray emission of gravitational galactic wakes in MOND and in Newtonian gravity with a dark halo are highlighted.  相似文献   

4.
We test the modified Newtonian dynamics (MOND) theory with the velocity dispersion profiles of Galactic globular clusters populating the outermost region of the Milky Way halo, where the Galactic acceleration is lower than the characteristic MOND acceleration a 0. For this purpose, we constructed self-consistent, spherical models of stellar systems in MOND, which are the analogues of the Newtonian King models. The models are spatially limited, reproduce well the surface brightness profiles of globular clusters and have velocity dispersion profiles that differ remarkably in shape from the corresponding Newtonian models. We present dynamical models of six globular clusters, which can be used to efficiently test MOND with the available observing facilities. A comparison with recent spectroscopic data obtained for NGC 2419 suggests that the kinematics of this cluster might be hard to explain in MOND.  相似文献   

5.
The stability of the dynamical trajectories of softened spherical gravitational systems is examined, both in the case of the full N -body problem and that of trajectories moving in the gravitational field of non-interacting background particles. In the latter case, for   N 10 000  , some trajectories, even if unstable, had exceedingly long diffusion times, which correlated with the characteristic e-folding time-scale of the instability. For trajectories of   N ≈100 000  systems this time-scale could be arbitrarily large – and thus appear to correspond to regular orbits. For centrally concentrated systems, low angular momentum trajectories were found to be systematically more unstable. This phenomenon is analogous to the well-known case of trajectories in generic centrally concentrated non-spherical smooth systems, where eccentric trajectories are found to be chaotic. The exponentiation times also correlate with the conservation of the angular momenta along the trajectories. For N up to a few hundred, the instability time-scales of N -body systems and their variation with particle number are similar to those of the most chaotic trajectories in inhomogeneous non-interacting systems. For larger N (up to a few thousand) the values of the these time-scales were found to saturate, increasing significantly more slowly with N . We attribute this to collective effects in the fully self-gravitating problem, which are apparent in the time variations of the time-dependent Liapunov exponents. The results presented here go some way towards resolving the long-standing apparent paradoxes concerning the local instability of trajectories. This now appears to be a manifestation of mechanisms driving evolution in gravitational systems and their interactions – and may thus be a useful diagnostic of such processes.  相似文献   

6.
We investigate the mean velocity dispersion and the velocity dispersion profile of stellar systems in modified Newtonian dynamics (MOND), using the N -body code n-mody , which is a particle-mesh-based code with a numerical MOND potential solver developed by Ciotti, Londrillo & Nipoti. We have calculated mean velocity dispersions for stellar systems following Plummer density distributions with masses in the range of 104 to  109 M⊙  and which are either isolated or immersed in an external field. Our integrations reproduce previous analytic estimates for stellar velocities in systems in the deep MOND regime  ( a i, a e≪ a 0)  , where the motion of stars is either dominated by internal accelerations  ( a i≫ a e)  or constant external accelerations  ( a e≫ a i)  . In addition, we derive for the first time analytic formulae for the line-of-sight velocity dispersion in the intermediate regime  ( a i∼ a e∼ a 0)  . This allows for a much-improved comparison of MOND with observed velocity dispersions of stellar systems. We finally derive the velocity dispersion of the globular cluster Pal 14 as one of the outer Milky Way halo globular clusters that have recently been proposed as a differentiator between Newtonian and MONDian dynamics.  相似文献   

7.
The dynamical friction problem is a long-standing dilemma about globular clusters (hereafter GCs) belonging to dwarf galaxies. GCs are strongly affected by dynamical friction in dwarf galaxies, and are presumed to fall into the galactic centre. But, GCs do exist in dwarf galaxies generally. A solution of the problem has been proposed. If dwarf galaxies have a core dark matter halo which has constant density distribution in its centre, the effect of dynamical friction will be weakened considerably, and GCs should be able to survive beyond the age of the Universe. Then, the solution argued that, in a cored dark halo, interaction between the halo and the GC constructs a new equilibrium state, in which a part of the halo rotates along with the GC (corotating state). The equilibrium state can suppress the dynamical friction in the core region. In this study, I tested whether the solution is reasonable and reconsidered why a constant density, core halo suppresses dynamical friction, by means of N -body simulations. As a result, I conclude that the true mechanism of suppressed dynamical friction is not the corotating state, although a core halo can actually suppress dynamical friction on GCs significantly.  相似文献   

8.
We performed N -body simulations of star cluster encounters with Hernquist's TREECODE in a CRAY YMP-2E computer under different initial conditions (relative positions and velocities, cluster sizes, masses and concentration degrees). The total number of particles per simulation ranged from 1024 to 20480. These models are compared with a series of isodensity maps of cluster pairs in the Magellanic Clouds. Evidence is found that during the interactions, transient morphological effects such as an expanded halo, isophotal deformation and isophotal twisting can occur as a result of tidal effects and dynamical friction. The simulations also show that different outcomes are possible depending on the initial parameters: (i) long-standing changes of concentration degree can occur after the collision; (ii) one member can disaggregate; or (iii) the pair can coalesce into a single cluster with a distinct structure compared with the original ones. These simulations can reproduce a wide range of morphological structures in observed cluster pairs.  相似文献   

9.
The dynamical evolution of small stellar groups composed of N=6 components was numerically simulated within the framework of a gravitational N-body problem. The effects of stellar mass loss in the form of stellar wind, dynamical friction against the interstellar medium, and star mergers on the dynamical evolution of the groups were investigated. A comparison with a purely gravitational N-body problem was made. The state distributions at the time of 300 initial system crossing times were analyzed. The parameters of the forming binary and stable triple systems as well as the escaping single and binary stars were studied. The star-merger and dynamical-friction effects are more pronounced in close systems, while the stellar wind effects are more pronounced in wide systems. Star-mergers and stellar wind slow down the dynamical evolution. These factors cause the mean and median semimajor axes of the final binaries as well as the semimajor axes of the internal and external binaries in stable triple systems to increase. Star mergers and dynamical friction in close systems decrease the fraction of binary systems with highly eccentric orbits and the mean component mass ratios for the final binaries and the internal and external binaries in stable triple systems. Star mergers and dynamical friction in close systems increase the fraction of stable triple systems with prograde motions. Dynamical friction in close systems can both increase and decrease the mean velocities of the escaping single stars, depending on the density of the interstellar medium and the mean velocity of the stars in the system.  相似文献   

10.
We investigate the possibility of discriminating between modified Newtonian dynamics (MOND) and Newtonian gravity with dark matter, by studying the vertical dynamics of disc galaxies. We consider models with the same circular velocity in the equatorial plane (purely baryonic discs in MOND and the same discs in Newtonian gravity embedded in spherical dark matter haloes), and we construct their intrinsic and projected kinematical fields by solving the Jeans equations under the assumption of a two-integral distribution function. We find that the vertical velocity dispersion of deep MOND discs can be much larger than in the equivalent spherical Newtonian models. However, in the more realistic case of high surface density discs, this effect is significantly reduced, casting doubt on the possibility of discriminating between MOND and Newtonian gravity with dark matter by using current observations.  相似文献   

11.
We demonstrate that the chaotic nature of N -body systems can lead to macroscopic variations in the evolution of collisionless simulations containing rotationally supported discs. The unavoidable stochasticity that afflicts all simulations generally causes mild differences between the evolution of similar models but, in order to illustrate that this is not always true, we present a case that shows an extreme bimodal divergence. The divergent behaviour occurs in two different types of code, and is independent of all numerical parameters. We identify and give explicit illustrations of several sources of stochasticity, and also show that macroscopic variations in the evolution can originate from differences at the round-off error level. We obtain somewhat more consistent results from simulations in which the halo is set-up with great care compared with those started from more approximate equilibria, but we have been unable to eliminate diverging behaviour entirely because the main sources of stochasticity are intrinsic to the disc. We show that the divergence is only temporary and that halo friction is merely delayed, for a substantial time in some cases. We argue that the delays are unlikely to arise in real galaxies, and that our results do not affect dynamical friction constraints on halo density. Stochastic variations in the evolution are inevitable in all simulations of disc–halo systems, irrespective of how they were created, although their effect is generally far less extreme than we find here. The possibility of divergent behaviour complicates comparison of results from different workers.  相似文献   

12.
The evolution of substructure embedded in non-dissipative dark haloes is studied through N -body simulations of isolated systems, both in and out of initial equilibrium, complementing cosmological simulations of the growth of structure. We determine by both analytic calculations and direct analysis of the N -body simulations the relative importance of various dynamical processes acting on the clumps, such as the removal of material by global tides, clump–clump heating, clump–clump merging and dynamical friction. The ratio of the internal clump velocity dispersion to that of the dark halo is an important parameter; as this ratio approaches a value of unity, heating by close encounters between clumps becomes less important, while the other dynamical processes continue to increase in importance. Our comparison between merging and disruption processes implies that spiral galaxies cannot be formed in a protosystem that contains a few large clumps, but can be formed through the accretion of many small clumps; elliptical galaxies form in a more clumpy environment than do spiral galaxies. Our results support the idea that the central cusp in the density profiles of dark haloes is the consequence of self-limiting merging of small, dense haloes. This implies that the collapse of a system of clumps/substructure is not sufficient to form a cD galaxy, with an extended envelope; plausibly, subsequent accretion of large galaxies is required. The post-collapse system is in general triaxial, with rounder systems resulting from fewer, but more massive, clumps. Persistent streams of material from disrupted clumps can be found in the outer regions of the final system, and at an overdensity of around 0.75, can cover 10 to 30 per cent of the sky.  相似文献   

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

14.
We study the efficiency at which a black hole or dense star cluster spirals in to the Galactic Centre. This process takes place on a dynamical friction time-scale, which depends on the value of the Coulomb logarithm (ln Λ). We determine the accurate value of this parameter using the direct N -body method, a tree algorithm and a particle-mesh technique with up to two million plus one particles. The three different techniques are in excellent agreement. Our measurement for the Coulomb logarithm appears to be independent of the number of particles. We conclude that  ln Λ= 6.6 ± 0.6  for a massive point particle in the inner few parsec of the Galactic bulge. For an extended object, such as a dense star cluster, ln Λ is smaller, with a value of the logarithm argument Λ inversely proportional to the object size.  相似文献   

15.
N -body simulations are widely used to simulate the dynamical evolution of a variety of systems, among them star clusters. Much of our understanding of their evolution rests on the results of such direct N -body simulations. They provide insight in the structural evolution of star clusters, as well as into the occurrence of stellar exotica. Although the major pure N -body codes starlab/kira and nbody4 are widely used for a range of applications, there is no thorough comparison study yet.
Here, we thoroughly compare basic quantities as derived from simulations performed either with starlab/kira or nbody4 .
We construct a large number of star cluster models for various stellar mass function settings (but without stellar/binary evolution, primordial binaries, external tidal fields, etc.), evolve them in parallel with starlab/kira and nbody4 , analyse them in a consistent way and compare the averaged results quantitatively. For this quantitative comparison, we develop a bootstrap algorithm for functional dependencies.
We find an overall excellent agreement between the codes, both for the clusters' structural and energy parameters as well as for the properties of the dynamically created binaries. However, we identify small differences, like in the energy conservation before core collapse and the energies of escaping stars, which deserve further studies.
Our results reassure the comparability and the possibility to combine results from these two major N -body codes, at least for the purely dynamical models (i.e. without stellar/binary evolution) we performed. Further detailed comparison studies for more complex systems, e.g. including stellar/binary evolution, are required.  相似文献   

16.
The time-scale for galaxies within merging dark matter haloes to merge with each other is an important ingredient in galaxy formation models. Accurate estimates of merging time-scales are required for predictions of astrophysical quantities such as black hole binary merger rates, the build-up of stellar mass in central galaxies and the statistical properties of satellite galaxies within dark matter haloes. In this paper, we study the merging time-scales of extended dark matter haloes using N -body simulations. We compare these results to standard estimates based on the Chandrasekhar theory of dynamical friction. We find that these standard predictions for merging time-scales, which are often used in semi-analytic galaxy formation models, are systematically shorter than those found in simulations. The discrepancy is approximately a factor of 1.7 for M sat/ M host≈ 0.1 and becomes larger for more disparate satellite-to-host mass ratios, reaching a factor of ∼3.3 for M sat/ M host≈ 0.01. Based on our simulations, we propose a new, easily implementable fitting formula that accurately predicts the time-scale for an extended satellite to sink from the virial radius of a host halo down to the halo's centre for a wide range of M sat/ M host and orbits. Including a central bulge in each galaxy changes the merging time-scale by ≲10 per cent. To highlight one concrete application of our results, we show that merging time-scales often used in the literature overestimate the growth of stellar mass by satellite accretion by ≈40 per cent, with the extra mass gained in low mass ratio mergers.  相似文献   

17.
The modified Newtonian dynamics (MOND), suggested by Milgrom as an alternative to dark matter, implies that isothermal spheres with a fixed anisotropy parameter should exhibit a near-perfect relation between the mass and velocity dispersion of the form M ∝ σ  4. This is consistent with the observed Faber–Jackson relation for elliptical galaxies: a luminosity–velocity dispersion relation with large scatter. However, the observable global properties of elliptical galaxies comprise a three-parameter family; they lie on a 'fundamental plane' in a logarithmic space consisting of central velocity dispersion, effective radius ( r e) and luminosity. The scatter perpendicular to this plane is significantly less than that about the Faber–Jackson relation. I show here that, in order to match the observed properties of elliptical galaxies with MOND, models must deviate from being strictly isothermal and isotropic; such objects can be approximated by high-order polytropic spheres with a radial orbit anisotropy in the outer regions. MOND imposes boundary conditions on the inner Newtonian regions which restrict these models to a dynamical fundamental plane of the form where the exponents may differ from the Newtonian expectations ( α =2, γ =1). Scatter about this plane is relatively insensitive to the necessary deviations from homology.  相似文献   

18.
The dynamical friction experienced by a body moving in a gaseous medium is different from the friction in the case of a collisionless stellar system. Here we consider the orbital evolution of a gravitational perturber inside a gaseous sphere using three-dimensional simulations, ignoring however self-gravity. The results are analysed in terms of a 'local' formula with the associated Coulomb logarithm taken as a free parameter. For forced circular orbits, the asymptotic value of the component of the drag force in the direction of the velocity is a slowly varying function of the Mach number in the range 1.0–1.6. The dynamical friction time-scale for free decay orbits is typically only half as long as in the case of a collisionless background, which is in agreement with E. C. Ostriker's recent analytic result. The orbital decay rate is rather insensitive to the past history of the perturber. It is shown that, similarly to the case of stellar systems, orbits are not subject to any significant circularization. However, the dynamical friction time-scales are found to increase with increasing orbital eccentricity for the Plummer model, whilst no strong dependence on the initial eccentricity is found for the isothermal sphere.  相似文献   

19.
In this paper, we extend our numerical method for simulating terrestrial planet formation to include dynamical friction from the unresolved debris component. In the previous work, we implemented a rubble pile planetesimal collision model into direct N -body simulations of terrestrial planet formation. The new collision model treated both accretion and erosion of planetesimals but did not include dynamical friction from debris particles smaller than the resolution limit for the simulation. By extending our numerical model to include dynamical friction from the unresolved debris, we can simulate the dynamical effect of debris produced during collisions and can also investigate the effect of initial debris mass on terrestrial planet formation. We find that significant initial debris mass, 10 per cent or more of the total disc mass, changes the mode of planetesimal growth. Specifically, planetesimals in this situation do not go through a runaway growth phase. Instead, they grow concurrently, similar to oligarchic growth. The dynamical friction from the unresolved debris damps the eccentricities of the planetesimals, reducing the mean impact speeds and causing all collisions to result in merging with no mass loss. As a result, there is no debris production. The mass in debris slowly decreases with time. In addition to including the dynamical friction from the unresolved debris, we have implemented particle tracking as a proxy for monitoring compositional mixing. Although there is much less mixing due to collisions and gravitational scattering when dynamical friction of the background debris is included, there is significant inward migration of the largest protoplanets in the most extreme initial conditions (for which the initial mass in unresolved debris is at least equal to the mass in resolved planetesimals).  相似文献   

20.
We apply the modified acceleration law obtained from Einstein gravity coupled to a massive skew symmetric field,   F μνλ  , to the problem of explaining X-ray galaxy cluster masses without exotic dark matter. Utilizing X-ray observations to fit the gas mass profile and temperature profile of the hot intracluster medium (ICM) with King 'β-models', we show that the dynamical masses of the galaxy clusters resulting from our modified acceleration law fit the cluster gas masses for our sample of 106 clusters without the need of introducing a non-baryonic dark matter component. We are further able to show for our sample of 106 clusters that the distribution of gas in the ICM as a function of radial distance is well fitted by the dynamical mass distribution arising from our modified acceleration law without any additional dark matter component. In a previous work, we applied this theory to galaxy rotation curves and demonstrated good fits to our sample of 101 low surface brightness, high surface brightness and dwarf galaxies including 58 galaxies that were fitted photometrically with the single-parameter mass-to-light ratio ( M / L )stars. The results obtained there were qualitatively similar to those obtained using Milgrom's phenomenological Modified Newtonian Dynamics (MOND) model, although the determined galaxy masses were quantitatively different, and MOND does not show a return to Keplerian behaviour at extragalactic distances. The results obtained here are compared to those obtained using Milgrom's phenomenological MOND model which does not fit the X-ray galaxy cluster masses unless an auxiliary dark matter component is included.  相似文献   

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

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