共查询到20条相似文献，搜索用时 31 毫秒
1.
R. H. Miller 《Celestial Mechanics and Dynamical Astronomy》1988,45(1-3):19-26
Computation and a wealth of new observational techniques have reinvigorated dynamical studies of galaxies and star clusters. These objects are examples of the gravitationaln-body problem withn in the range from a few hundred to 1011. Relaxation effects dominate at the low end and are completely negligible at the high end. The gravitationaln-body problem is chaotic, and the principal challenge in doing physics where that problem is involved (whether computationally or with analytic theory) is to ensure that chaos has not vitiated the results. Enforcing a Liouville theorem accomplishes this with collision-free large-n problems, but equivalent recipes are not in common use for smallern. We describe some important insights and discoveries that have come from computation in stellar dynamics, discuss chaos briefly, and indicate the way the physics that comes up in different astronomical contexts is addressed in numerical methods currently in use. Graphics is a vital part of any computational approach. The long range prospects are very promising for continued high scientific productivity in stellar dynamics. 相似文献
2.
Luis A. Aguilar 《Celestial Mechanics and Dynamical Astronomy》1987,41(1-4):3-37
Summary A brief review of previous work and the present situation in the problem of formation of elliptical galaxies via dissipationless collapse are presented, as well as the results of a new set of numerical experiments. It is shown that collapses started from cold initial conditions are different from warmer collapses, due to the presence of a dynamical instability associated with radial orbits. This instability leads to triaxial final configurations, regardless of the initial amount of random kinetic energy, rotational kinetic energy, or shape of the initial conditions, as long as2T/W0.1, whereT is the total (rotational plus thermal) kinetic energy andW is the potential energy of the initial conditions. Warmer initial conditions preserve their initial shape, or become oblate if initially rotating. Cold initial conditions produce equilibrium systems with realistic density profiles, as opposed to collapses from warmer conditions that result in core-halo profiles, unlike the observed surface brightness profiles of elliptical galaxies. Although the same cold collapses that result in triaxial shapes produce realistic density profiles, it is shown that these two effects are not directly connected: cold collapses simulated with anN-body code that enforces spherical symmetry result in realistic density profiles too.Invited Talk at the 18th meeting of the Division on Dynamical Astronomy of the American Astronomical Society held in Braintree, Massachussetts, May 7–9, 1987. 相似文献
3.
Euaggelos E. Zotos 《New Astronomy》2011,16(6):391-401
In the present article, we present a new gravitational galactic model, describing motion in elliptical as well as in disk galaxies, by suitably choosing the dynamical parameters. Moreover, a new dynamical parameter, the S(g) spectrum, is introduced and used, in order to detect islandic motion of resonant orbits and the evolution of the sticky regions. We investigate the regular or chaotic character of motion, with emphasis in the different dynamical models and make an extensive study of the sticky regions of the system. We use the classical method of the Poincaré r ? pr phase plane and the new dynamical parameter of the S(g) spectrum. The L.C.E is used, in order to make an estimation of the degree of chaos in our galactic model. In both cases, the numerical calculations, suggest that our new model, displays a wide variety of families of regular orbits, compared to other galactic models. In addition to the regular motion, this new model displays also chaotic regions. Furthermore, the extent of the chaotic regions increases, as the value of the flatness parameter b of the model increases. Moreover, our simulations indicate, that the degree of chaos in elliptical galaxies, is much smaller than that in dense disk galaxies. In both cases numerical calculations show, that the degree of chaos increases linearly, as the flatness parameter b increases. In addition, a linear relationship between the critical value of angular momentum and the b parameter if found, in both cases (elliptical and disk galaxies). Some theoretical arguments to support the numerical outcomes are presented. Comparison with earlier work is also made. 相似文献
4.
SeveralN-body experiments were performed in order to simulate the dynamical behaviour of systems of galaxies gravitationally dominated by a massive dark background.We discuss mass estimates from the dynamics of the luminous component (M
VT) under the influence of such a background, assuming a constant dark/luminous mass ratio (M
D/M
L) and plausible physical conditions. We extend in this way previous studies (Smith, 1980, 1984) about the dependence ofM
VT on the relative distributions of dark and luminous matter (Limber, 1959). We found that the observed ratio of the virial theorem mass to luminosity (M
VT/L) in systems of galaxies of different sizes could be the result of different stages of their post-virialisation evolution as was previously suggested by White and Rees (1978) and Barnes (1983). This evolution is mainly the result of the dynamical friction that dark matter exerts on the luminous component. Thus our results give support to the idea that compact groups of galaxies are dynamically more evolved than large clusters, which is expected from the hierarchical clustering picture for the formation of such structures. 相似文献
5.
The method for measuring the dynamical age of n-body systems proposed by S. von HOERNER (1976) is applied to clusters of galaxies. The modulus of evolution W characterizing quantitatively the dynamical state of the cluster was determined for several objects. It turns out that the Virgo cluster is one of the dynamically youngest clusters. 相似文献
6.
The dynamical evolution of theN-body systems (N=6) has been studied by numerical simulations. The double and triple subsystems isolated from other bodies have been revealed. The processes of their formation, evolution, and disruption have been followed. At the initial stage of evolution, a common collapse of the system takes place. After that some fast triple approaches of single bodies are possible. At further stages, some physically connected triple subsystems are formed which are disrupted by the intruders or due to their own dynamical unstability. The dissipation mechanisms formulated by van Albada (1968) and Aarseth (1973) have been confirmed. At the final stage, double or hierarchical triple systems are formed. 相似文献
7.
The oldest open clusters in our Galaxy set the lower limit to the age of the Galactic Disk (9–10 Gyr). Although they appear
to be very rich now, it is clear that their primordial populations were much larger. Often considered as transitional objects,
these populous open clusters show structural differences with respect to globular clusters so their dynamics and characteristic
evolutionary time scales can also be different. On the other hand, their large membership lead to different dynamical evolution
as compared with average open clusters. In this paper, the differential features of the evolution of rich open clusters are
studied using N-body simulations, including several of the largest (104 stars) published direct collisional N-body calculations so far, which were performed on a CRAY YMP. The disruption rate of rich open clusters is analysed in detail
and the effect of the initial spatial distribution of the stars in the cluster on its dynamics is studied. The results show
that cluster life-time depends on this initial distribution, decreasing when it is more concentrated. The effect of stellar
evolution on the dynamical evolution of rich clusters is an important subject that also has been considered here. We demonstrate
that the cluster's life-expectancy against evaporation increases because of mass loss by evolving high-mass stars.
This revised version was published online in July 2006 with corrections to the Cover Date. 相似文献
8.
We analyze the relationship between the mass of a spherical component and the minimum possible thickness of stable stellar disks. This relationship for real galaxies allows the lower limit on the dark halo mass to be estimated (the thinner the stable stellar disk is, the more massive the dark halo must be). In our analysis, we use both theoretical relations and numerical N-body simulations of the dynamical evolution of thin disks in the presence of spherical components with different density profiles and different masses. We conclude that the theoretical relationship between the thickness of disk galaxies and the mass of their spherical components is a lower envelope for the model data points. We recommend using this theoretical relationship to estimate the lower limit for the dark halo mass in galaxies. The estimate obtained turns out to be weak. Even for the thinnest galaxies, the dark halo mass within four exponential disk scale lengths must be more than one stellar disk mass. 相似文献
9.
N. Voglis M. Harsoula CH. Efthymiopoulos 《Celestial Mechanics and Dynamical Astronomy》2000,78(1-4):265-278
As it is known, a good number of galaxies are observed to have counterrotating cores. A popular scenario to explain the formation of such galaxies is based on a secondary process of merging of galaxies with their satellites, or gas infall, or merger events between galaxies. An alternative mechanism, proposed by Voglis et al., 1991, and by Harsoula and Voglis 1998, could also be responsible for the formation of these galaxies directly from cosmological initial conditions (direct scenario). The novel mechanism was demonstrated by using quiet cosmological initial conditions in N-body simulations. In the present paper we extend our N-body simulations using clumpy initial conditions and show that this mechanism still works to create counterrotating galaxies. Counterrotation is a result of the considerable amount of memory of initial conditions surviving for times comparable to the Hubble time, despite the large degree of instability of individual orbits and the dramatic redistribution and mixing of the particles in phase space. We show, for example, that the particles remember, in a statistical sense, not only their distance from the center of mass (memory of energy), but also the initial orientation of their position relative to the direction of an external tidal field, which determines the sign and the amount of angular momentum that is transferred to the particles of the system. 相似文献
10.
Bing Zhang Rosemary F. G. Wyse Massimo Stiavelli Joseph Silk 《Monthly notices of the Royal Astronomical Society》2002,332(3):647-675
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. 相似文献
11.
G. Contopoulos N. Voglis C. Kalapotharakos 《Celestial Mechanics and Dynamical Astronomy》2002,83(1-4):191-204
We construct and compare two different self-consistent N-body equilibrium configurations of galactic models. The two systems have their origin in cosmological initial conditions selected so that the radial orbit instability appears in one model and gives an E5 type elliptical galaxy, but not in the other that gives an E1 type. We examine their phase spaces using uniformly distributed orbits of test particles in the resulting potential and compare with the distribution of the orbits of the real particles in the two systems. The main types of orbits in both cases are box, tube and chaotic orbits. One main conclusion is that the orbits of the test particles in the 3-dimensional potential are foliated in a way quite close to the foliation of invariant tori in a 2-dimensional potential. The real particles describe orbits having similar foliation. However, their distribution is far from being uniform. The difference between the two models of equilibrium is realized mainly by different balances of the populations of real particles in box and tube orbits. 相似文献
12.
《New Astronomy》2021
We present an N-body code called Taichi for galactic dynamics and controlled numerical experiments. The code includes two high-order hierarchical multipole expansion methods: the Barnes-Hut (BH) tree and the fast multipole method (FMM). For the time integration, the code can use either a conventional adaptive KDK or a Hamiltonian splitting integrator. The combination of FMM and the Hamiltonian splitting integrator leads to a momentum-conserving N-body scheme with individual time steps. We find Taichi performs well in the typical applications in galactic dynamics. In the isolated and interacting galaxies tests, the momentum conserving scheme produces the same result as a conventional BH tree code. But for similar force accuracies, FMM significantly speeds up the simulations compared to the monopole BH tree. In the cold collapse test, we find the inner structure after relaxation can be sensitive to the force accuracies. Taichi is ready to incorporate special treatment of close encounters thanks to the Hamiltonian splitting integrator, suitable for studying dynamics around central massive bodies. 相似文献
13.
Mercury's capture into the 3:2 spin–orbit resonance can be explained as a result of its chaotic orbital dynamics. One major objective of MESSENGER and BepiColombo spatial missions is to accurately measure Mercury's rotation and its obliquity in order to obtain constraints on internal structure of the planet. Analytical approaches at the first-order level using the Cassini state assumptions give the obliquity constant or quasi-constant. Which is the obliquity's dynamical behavior deriving from a complete spin–orbit motion of Mercury simultaneously integrated with planetary interactions? We have used our SONYR model (acronym of Spin–Orbit N-bodY Relativistic model) integrating the spin–orbit N-body problem applied to the Solar System (Sun and planets). For lack of current accurate observations or ephemerides of Mercury's rotation, and therefore for lack of valid initial conditions for a numerical integration, we have built an original method for finding the libration center of the spin–orbit system and, as a consequence, for avoiding arbitrary amplitudes in librations of the spin–orbit motion as well as in Mercury's obliquity. The method has been carried out in two cases: (1) the spin–orbit motion of Mercury in the 2-body problem case (Sun–Mercury) where an uniform precession of the Keplerian orbital plane is kinematically added at a fixed inclination (S2K case), (2) the spin–orbit motion of Mercury in the N-body problem case (Sun and planets) (Sn case). We find that the remaining amplitude of the oscillations in the Sn case is one order of magnitude larger than in the S2K case, namely 4 versus 0.4 arcseconds (peak-to-peak). The mean obliquity is also larger, namely 1.98 versus 1.80 arcminutes, for a difference of 10.8 arcseconds. These theoretical results are in a good agreement with recent radar observations but it is not excluded that it should be possible to push farther the convergence process by drawing nearer still more precisely to the libration center. We note that the dynamically driven spin precession, which occurs when the planetary interactions are included, is more complex than the purely kinematic case. Nevertheless, in such a N-body problem, we find that the 3:2 spin–orbit resonance is really combined to a synchronism where the spin and orbit poles on average precess at the same rate while the orbit inclination and the spin axis orientation on average decrease at the same rate. As a consequence and whether it would turn out that there exists an irreducible minimum of the oscillation amplitude, quasi-periodic oscillations found in Mercury's obliquity should be to geometrically understood as librations related to these synchronisms that both follow a Cassini state. Whatever the open question on the minimal amplitude in the obliquity's oscillations and in spite of the planetary interactions indirectly acting by the solar torque on Mercury's rotation, Mercury remains therefore in a stable equilibrium state that proceeds from a 2-body Cassini state. 相似文献
14.
We have performed a series ofN-body experiments including the effects of a massive dominant background which follows Schuster's density law in order to simulate clusters of galaxies in which a smoothly distributed dark component is present. The existence of this background is inferred from the weak luminosity segregation observed in clusters which, however, show several characteristics of well-relaxed systems. The comparison of the velocity dispersion profiles of three clusters of galaxies (Coma, Perseus, and Virgo) with those obtained in the numerical experiments allows us to place some constraints on both the distribution and amount of distributed dark material in these clusters. The profiles are rather insensitive to variations in the ratio of the background mass to the mass attached to galaxies (M
b/Mg), but exhibit a strong dependence on their relative concentration. We conclude that highly concentrated background models are not consistent with observations. We find a maximum value for the ratio of the gravitational radius of the galaxies and the background (R
g/Rb) (approximately 0.6) and using previous results (Navarroet al., 1986) we conclude that virial theorem masses underestimate the total mass (M
b+M
g) of the clusters. As a final result, we derive a minimum value for theM
b/Mg ratio. All these conclusions could apply in general if Coma, Perseus, and Virgo constitute a fair sample of the rich clusters of galaxies in the Universe. 相似文献
15.
Simon Portegies Zwart Steve McMillan Stefan Harfst Derek Groen Michiko Fujii Breanndán Ó Nualláin Evert Glebbeek Douglas Heggie James Lombardi Piet Hut Vangelis Angelou Sambaran Banerjee Houria Belkus Tassos Fragos John Fregeau Evghenii Gaburov Rob Izzard Mario Jurić Stephen Justham Andrea Sottoriva Marcel Zemp 《New Astronomy》2009,14(4):369-378
We present MUSE, a software framework for combining existing computational tools for different astrophysical domains into a single multiphysics, multiscale application. MUSE facilitates the coupling of existing codes written in different languages by providing inter-language tools and by specifying an interface between each module and the framework that represents a balance between generality and computational efficiency. This approach allows scientists to use combinations of codes to solve highly coupled problems without the need to write new codes for other domains or significantly alter their existing codes. MUSE currently incorporates the domains of stellar dynamics, stellar evolution and stellar hydrodynamics for studying generalized stellar systems. We have now reached a “Noah’s Ark” milestone, with (at least) two available numerical solvers for each domain. MUSE can treat multiscale and multiphysics systems in which the time- and size-scales are well separated, like simulating the evolution of planetary systems, small stellar associations, dense stellar clusters, galaxies and galactic nuclei. In this paper we describe three examples calculated using MUSE: the merger of two galaxies, the merger of two evolving stars, and a hybrid N-body simulation. In addition, we demonstrate an implementation of MUSE on a distributed computer which may also include special-purpose hardware, such as GRAPEs or GPUs, to accelerate computations. The current MUSE code base is publicly available as open source at http://muse.li. 相似文献
16.
We have compared the frequency distribution of the dynamical observedquantity log (V
z
2
r
p), for a sample of 46 pairs of elliptical galaxies, to the distribution of this quantity obtained from numerical simulations
of pairs of galaxies. From such an analysis, where we have considered the structure of the galaxies and its influence in the
orbital evolution of the system, we have obtained the characteristic mass and the mass-luminosity ratio for the sample. Our
results show that the hypothesis of point-mass in elliptical orbits is, for this sample, an approximation as good as the model
that takes into account the structure of the galaxies. The statistical method used here gives an estimate of a more reliable
mass, it minimizes the contamination of spurious pairs and it considers adequately the contribution of the physical pairs.
We have obtained a characteristic mass to the 46 elliptical pairs of 1.68 × 1012 ± 7.01 ×1011 M⊙ with M/L = 17.6 ± 7.3 (H
0 = 60 km s-1Mpc-1).
This revised version was published online in July 2006 with corrections to the Cover Date. 相似文献
17.
We present a new method for fast numerical integration of close binaries inN-body systems. The basic idea is to slow down the motion of the binary artificially, which makes a faster numerical integration possible but still maintains correct treatment of secular and long-period effects on the motion. We discuss the general principle, with application to close binaries inN-body codes and in the chain regularization. 相似文献
18.
We present a time-transformed leapfrog scheme combined with the extrapolation method to construct an integrator for orbits in N-body systems with large mass ratios. The basic idea can be used to transform any second-order differential equation into a form which may allow more efficient numerical integration. When applied to gravitating few-body systems this formulation permits extremely close two-body encounters to be considered without significant loss of accuracy. The new scheme has been implemented in a direct N-body code for simulations of super-massive binaries in galactic nuclei. In this context relativistic effects may also be included. 相似文献
19.
R. H. Miller 《Celestial Mechanics and Dynamical Astronomy》1985,37(3):307-324
Numerical experiments undertaken to investigate the longevity and behavior of dark-lane elliptical galaxies are described. This is dynamically the same problem as a disk galaxy in a massive halo. Spiral galaxies are disks from a dynamical point of view. A disk of particles embedded in a self-consistent galaxy provides the basic model used for the experiments. This model is applicable to ordinary disk galaxies if the disk is interpreted as the visible galaxy and the galaxy is interpreted as the massive halo thought to be present around disk galaxies. Fully three-dimensional fully self-consistentn-body computer programs that can handle 100,000 particles are used for the experiments. The background galaxy is oblate, and the disk is inclined to the axis of the oblate galaxy, so the disk precesses differentially to produce a warp. A surprising result is that the galaxy center shifted, leaving the disk center orbiting around the galaxy center. This produces interesting phenomena reminiscent of observations in the region of the Galactic center. 相似文献
20.
A. G. Bedregal A. Aragón-Salamanca M. R. Merrifield N. Cardiel 《Monthly notices of the Royal Astronomical Society》2008,387(2):660-676
Using high signal-to-noise ratio VLT/FORS2 long-slit spectroscopy, we have studied the properties of the central stellar populations and dynamics of a sample of S0 galaxies in the Fornax cluster. The central absorption-line indices in these galaxies correlate well with the central velocity dispersions (σ0 ) in accordance with what previous studies found for elliptical galaxies. However, contrary to what it is usually assumed for cluster ellipticals, the observed correlations seem to be driven by systematic age and α-element abundance variations, and not changes in overall metallicity. We also found that the observed scatter in the index–σ0 relations can be partially explained by the rotationally supported nature of these systems. Indeed, even tighter correlations exist between the line indices and the maximum circular velocity of the galaxies. This study suggests that the dynamical mass is the physical property driving these correlations, and for S0 galaxies such masses have to be estimated assuming a large degree of rotational support. The observed trends imply that the most massive S0s have the shortest star formation time-scales and the oldest stellar populations. 相似文献