Self Similar Spherical Collapse Revisited: a Comparison between Gas and Dark Matter Dynamics

We reconsider the collapse of cosmic structures in an Einstein-de Sitter Universe, using the self similar initial conditions of Fillmore & Goldreich (1984). We first derive a new approximation to describe the dark matter dynamics in spherical geometry, that we refer to the "fluid approach". This method enables us to recover the self-similarity solutions of Fillmore & Goldreich for dark matter. We derive also new self-similarity solutions for the gas. We thus compare directly gas and dark matter dynamics, focusing on the differences due to their different dimensionalities in velocity space. This work may have interesting consequences for gas and dark matter distributions in large galaxy clusters, allowing to explain why the total mass profile is always steeper than the X-ray gas profile. We discuss also the shape of the dark matter density profile found in N-body simulations in terms of a change of dimensionality in the dark matter velocity space. The stable clustering hypothesis has been finally considered in the light of this analytical approach.Comment: 14 pages, 2 figures, accepted for publication in The Astrophysical Journa

Long-Term Evolution of Massive Black Hole Binaries. III. Binary Evolution in Collisional Nuclei

[Abridged] In galactic nuclei with sufficiently short relaxation times, binary supermassive black holes can evolve beyond their stalling radii via continued interaction with stars. We study this "collisional" evolutionary regime using both fully self-consistent N-body integrations and approximate Fokker-Planck models. The N-body integrations employ particle numbers up to 0.26M and a direct-summation potential solver; close interactions involving the binary are treated using a new implementation of the Mikkola-Aarseth chain regularization algorithm. Even at these large values of N, two-body scattering occurs at high enough rates in the simulations that they can not be simply scaled to the large-N regime of real galaxies. The Fokker-Planck model is used to bridge this gap; it includes, for the first time, binary-induced changes in the stellar density and potential. The Fokker-Planck model is shown to accurately reproduce the results of the N-body integrations, and is then extended to the much larger N regime of real galaxies. Analytic expressions are derived that accurately reproduce the time dependence of the binary semi-major axis as predicted by the Fokker-Planck model. Gravitational wave coalescence is shown to occur in <10 Gyr in nuclei with velocity dispersions below about 80 km/s. Formation of a core results from a competition between ejection of stars by the binary and re-supply of depleted orbits via two-body scattering. Mass deficits as large as ~4 times the binary mass are produced before coalescence. After the two black holes coalesce, a Bahcall-Wolf cusp appears around the single hole in one relaxation time, resulting in a nuclear density profile consisting of a flat core with an inner, compact cluster, similar to what is observed at the centers of low-luminosity spheroids.Comment: 21 page

Exact Quantum Solutions of Extraordinary N-body Problems

The wave functions of Boson and Fermion gases are known even when the particles have harmonic interactions. Here we generalise these results by solving exactly the N-body Schrodinger equation for potentials V that can be any function of the sum of the squares of the distances of the particles from one another in 3 dimensions. For the harmonic case that function is linear in r^2. Explicit N-body solutions are given when U(r) = -2M \hbar^{-2} V(r) = \zeta r^{-1} - \zeta_2 r^{-2}. Here M is the sum of the masses and r^2 = 1/2 M^{-2} Sigma Sigma m_I m_J ({\bf x}_I - {\bf x}_J)^2. For general U(r) the solution is given in terms of the one or two body problem with potential U(r) in 3 dimensions. The degeneracies of the levels are derived for distinguishable particles, for Bosons of spin zero and for spin 1/2 Fermions. The latter involve significant combinatorial analysis which may have application to the shell model of atomic nuclei. For large N the Fermionic ground state gives the binding energy of a degenerate white dwarf star treated as a giant atom with an N-body wave function. The N-body forces involved in these extraordinary N-body problems are not the usual sums of two body interactions, but nor are forces between quarks or molecules. Bose-Einstein condensation of particles in 3 dimensions interacting via these strange potentials can be treated by this method.Comment: 24 pages, Latex. Accepted for publication in Proceedings of the Royal Societ

Variation of the Dependence of the Transient Process Duration on the Initial Conditions in Systems with Discrete Time

Dependence of the transient process duration on the initial conditions is considered in one- and two-dimensional systems with discrete time, representing a logistic map and the Eno map, respectively.Comment: 4 pages, 2 figure

Partial suppression of the radial orbit instability in stellar systems

It is well known that the simple criterion proposed originally by Polyachenko and Shukhman (1981) for the onset of the radial orbit instability, although being generally a useful tool, faces significant exceptions both on the side of mildly anisotropic systems (with some that can be proved to be unstable) and on the side of strongly anisotropic models (with some that can be shown to be stable). In this paper we address two issues: Are there processes of collisionless collapse that can lead to equilibria of the exceptional type? What is the intrinsic structural property that is responsible for the sometimes noted exceptional stability behavior? To clarify these issues, we have performed a series of simulations of collisionless collapse that start from homogeneous, highly symmetrized, cold initial conditions and, because of such special conditions, are characterized by very little mixing. For these runs, the end-states can be associated with large values of the global pressure anisotropy parameter up to 2K_r/K_T \approx 2.75. The highly anisotropic equilibrium states thus constructed show no significant traces of radial anisotropy in their central region, with a very sharp transition to a radially anisotropic envelope occurring well inside the half-mass radius (around 0.2 r_M). To check whether the existence of such almost perfectly isotropic "nucleus" might be responsible for the apparent suppression of the radial orbit instability, we could not resort to equilibrium models with the above characteristics and with analytically available distribution function; instead, we studied and confirmed the stability of configurations with those characteristics by initializing N-body approximate equilibria (with given density and pressure anisotropy profiles) with the help of the Jeans equations.Comment: 26 pages, 9 figures, accepted for publication in The Astrophysical Journa

A Phase-Space Approach to Collisionless Stellar Systems Using a Particle Method

A particle method for reproducing the phase space of collisionless stellar systems is described. The key idea originates in Liouville's theorem which states that the distribution function (DF) at time t can be derived from tracing necessary orbits back to t=0. To make this procedure feasible, a self-consistent field (SCF) method for solving Poisson's equation is adopted to compute the orbits of arbitrary stars. As an example, for the violent relaxation of a uniform-density sphere, the phase-space evolution which the current method generates is compared to that obtained with a phase-space method for integrating the collisionless Boltzmann equation, on the assumption of spherical symmetry. Then, excellent agreement is found between the two methods if an optimal basis set for the SCF technique is chosen. Since this reproduction method requires only the functional form of initial DFs but needs no assumptions about symmetry of the system, the success in reproducing the phase-space evolution implies that there would be no need of directly solving the collisionless Boltzmann equation in order to access phase space even for systems without any special symmetries. The effects of basis sets used in SCF simulations on the reproduced phase space are also discussed.Comment: 16 pages w/4 embedded PS figures. Uses aaspp4.sty (AASLaTeX v4.0). To be published in ApJ, Oct. 1, 1997. This preprint is also available at http://www.sue.shiga-u.ac.jp/WWW/prof/hozumi/papers.htm

Two-dimensional dissipative maps at chaos threshold: Sensitivity to initial conditions and relaxation dynamics

The sensitivity to initial conditions and relaxation dynamics of two-dimensional maps are analyzed at the edge of chaos, along the lines of nonextensive statistical mechanics. We verify the dual nature of the entropic index for the Henon map, one ($q_{sen}<1$) related to its sensitivity to initial conditions properties, and the other, graining-dependent ($q_{rel}(W)>1$), related to its relaxation dynamics towards its stationary state attractor. We also corroborate a scaling law between these two indexes, previously found for $z$-logistic maps. Finally we perform a preliminary analysis of a linearized version of the Henon map (the smoothed Lozi map). We find that the sensitivity properties of all these $z$-logistic, Henon and Lozi maps are the same, $q_{sen}=0.2445...$Comment: Communication at NEXT2003, Second Sardinian International Conference on News and Expectations in Thermostatistics, Villasimius (Cagliari) Italy, 21-28 September 2003. Submitted to Physica A. Elsevier Latex, 8 pages, 6 eps figure

Two-dimensional maps at the edge of chaos: Numerical results for the Henon map

The mixing properties (or sensitivity to initial conditions) of two-dimensional Henon map have been explored numerically at the edge of chaos. Three independent methods, which have been developed and used so far for the one-dimensional maps, have been used to accomplish this task. These methods are (i)measure of the divergence of initially nearby orbits, (ii)analysis of the multifractal spectrum and (iii)computation of nonextensive entropy increase rates. The obtained results strongly agree with those of the one-dimensional cases and constitute the first verification of this scenario in two-dimensional maps. This obviously makes the idea of weak chaos even more robust.Comment: 4 pages, 3 figure

Analytical solutions of the lattice Boltzmann BGK model

Analytical solutions of the two dimensional triangular and square lattice Boltzmann BGK models have been obtained for the plain Poiseuille flow and the plain Couette flow. The analytical solutions are written in terms of the characteristic velocity of the flow, the single relaxation time $\tau$ and the lattice spacing. The analytic solutions are the exact representation of these two flows without any approximation.Comment: 10 pages, no postscript figure provide

Infinite ergodic theory and Non-extensive entropies

We bring into account a series of result in the infinite ergodic theory that we believe that they are relevant to the theory of non-extensive entropie