Revisiting the ABC flow dynamo

The ABC flow is a prototype for fast dynamo action, essential to the origin of magnetic field in large astrophysical objects. Probably the most studied configuration is the classical 1:1:1 flow. We investigate its dynamo properties varying the magnetic Reynolds number Rm. We identify two kinks in the growth rate, which correspond respectively to an eigenvalue crossing and to an eigenvalue coalescence. The dominant eigenvalue becomes purely real for a finite value of the control parameter. Finally we show that even for Rm = 25000, the dominant eigenvalue has not yet reached an asymptotic behaviour. Its still varies very significantly with the controlling parameter. Even at these very large values of Rm the fast dynamo property of this flow cannot yet be established

Irregular Satellites of the Planets: Products of Capture in the Early Solar System

All four giant planets in the Solar system possess irregular satellites, characterized by large, highly eccentric and/or inclined orbits that are distinct from the nearly circular, uninclined orbits of the regular satellites. This difference can be traced directly to different modes of formation. Whereas the regular satellites grew by accretion within circumplanetary disks the irregular satellites were captured from initially heliocentric orbits at an early epoch. Recently, powerful survey observations have greatly increased the number of known irregular satellites, permitting a fresh look at the group properties of these objects and motivating a re-examination of the mechanisms of capture. None of the suggested mechanisms, including gas-drag, pull-down, and three-body capture, convincingly fit the group characteristics of the irregular satellites. The sources of the satellites also remain unidentified.Comment: 51 pages, 17 figures, 5 tables, to appear in ARAA 200

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

Capture of field stars by globular clusters in dense bulge regions

The recent detection of a double Red Giant Branch in the optical color-magnitude diagram (CMD) of the bulge globular cluster HP1 (Ortolani et al. 1997), a more populated metal-poor steep one corresponding to the cluster itself, and another metal-rich curved, led us to explore in the present Letter the possibility of capture of field stars by a globular cluster orbiting in dense bulge regions over several gigayears. Analytical arguments, as well as N-body calculations for a cluster model of 10^5 solar masses in a bulge-like environment, suggest that a significant fraction of cluster stars may consist of captures. Metal-poor globular clusters in the inner bulge, like HP1, contrasting at least in Delta [Fe/H] = 1.0 dex with respect to the surrounding metal-rich stars, are ideal probes to further test the capture scenario. In turn, if this scenario is confirmed, the double RGB of HP1 could provide direct estimates of blanketing amounts, which is fundamental for the photometric calibration of metal-rich stellar populations.Comment: 6 pages, 2 included figures, aas2pp4,sty Latex style. To be published in Astrophysical Journal Letter

Two-Component Fokker-Planck Models for the Evolution of Isolated Globular Clusters

Two-component (normal and degenerate stars) models are the simplest realization of clusters with a mass spectrum because high mass stars evolve quickly into degenerates, while low mass stars remain on the main-sequence for the age of the universe. Here we examine the evolution of isolated globular clusters using two-component Fokker-Planck (FP) models that include heating by binaries formed in tidal capture and in three-body encounters. Three-body binary heating dominates and the postcollapse expansion is self-similar, at least in models with total mass M <= 3 x 10^5 M_\odot, initial half-mass radius r_{h,i} >= 5 pc, component mass ratio m_2/m_1 <= 2, and number ratio N_1/N_2 <= 300 when m_2=1.4 M_\odot. We derive scaling laws for \rho_c, v_c, r_c, and r_h as functions of m_1/m_2, N, M, and time t from simple energy-balance arguments, and these agree well with the FP simulations. We have studied the conditions under which gravothermal oscillations (GTOs) occur. If E_{tot} and E_c are the energies of the cluster and of the core, respectively, and t_{rh} and t_c are their relaxation times, then \epsilon \equiv (E_{tot}/t_{rh})/(E_c/t_{rc}) is a good predictor of GTOs: all models with \epsilon>0.01 are stable, and all but one with \epsilon < 0.01 oscillate. We derive a scaling law for \epsilon against N and m_1/m_2 and compared with our numerical results. Clusters with larger m_2/m_1 or smaller N are stabler.Comment: 15 pages (LaTeX) with 8 figures. To appear in ApJ March 10, 1998 issu

A Two-Temperature Model of the Intracluster Medium

We investigate evolution of the intracluster medium (ICM), considering the relaxation process between the ions and electrons. According to the standard scenario of structure formation, ICM is heated by the shock in the accretion flow to the gravitational potential well of the dark halo. The shock primarily heats the ions because the kinetic energy of an ion entering the shock is larger than that of an electron by the ratio of masses. Then the electrons and ions exchange the energy through coulomb collisions and reach the equilibrium. From simple order estimation we find that the region where the electron temperature is considerably lower than the ion temperature spreads out on a Mpc scale. We then calculate the ion and electron temperature profiles by combining the adiabatic model of two-temperature plasma by Fox & Loeb (1997) with spherically symmetric N-body and hydrodynamic simulations based on three different cosmological models. It is found that the electron temperature is about a half of the mean temperature at radii $\sim$ 1 Mpc. This could lead to an about 50 % underestimation in the total mass contained within $\sim$ 1 Mpc when the electron temperature profiles are used. The polytropic indices of the electron temperature profiles are $\simeq 1.5$ whereas those of mean temperature $\simeq 1.3$ for $r \geq 1$ Mpc. This result is consistent both with the X-ray observations on electron temperature profiles and with some theoretical and numerical predictions about mean temperature profiles.Comment: 20 pages with 6 figures. Accepted for publication in Ap

Dynamical stability analysis of the HD202206 system and constraints to the planetary orbits

Long-term precise Doppler measurements with the CORALIE spectrograph revealed the presence of two massive companions to the solar-type star HD202206. Although the three-body fit of the system is unstable, it was shown that a 5:1 mean motion resonance exists close to the best fit, where the system is stable. We present here an extensive dynamical study of the HD202206 system aiming at constraining the inclinations of the two known companions, from which we derive possible ranges of value for the companion masses. We study the long term stability of the system in a small neighborhood of the best fit using Laskar's frequency map analysis. We also introduce a numerical method based on frequency analysis to determine the center of libration mode inside a mean motion resonance. We find that acceptable coplanar configurations are limited to inclinations to the line of sight between 30 and 90 degrees. This limits the masses of both companions to roughly twice the minimum. Non coplanar configurations are possible for a wide range of mutual inclinations from 0 to 90 degrees, although $\Delta\Omega = 0 [\pi]$ configurations seem to be favored. We also confirm the 5:1 mean motion resonance to be most likely. In the coplanar edge-on case, we provide a very good stable solution in the resonance, whose $\chi^2$ does not differ significantly from the best fit. Using our method to determine the center of libration, we further refine this solution to obtain an orbit with a very low amplitude of libration, as we expect dissipative effects to have dampened the libration.Comment: 14 pages, 18 figure

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

Non-abelian plane waves and stochastic regimes for (2+1)-dimensional gauge field models with Chern-Simons term

An exact time-dependent solution of field equations for the 3-d gauge field model with a Chern-Simons (CS) topological mass is found. Limiting cases of constant solution and solution with vanishing topological mass are considered. After Lorentz boost, the found solution describes a massive nonlinear non-abelian plane wave. For the more complicate case of gauge fields with CS mass interacting with a Higgs field, the stochastic character of motion is demonstrated.Comment: LaTeX 2.09, 13 pages, 11 eps figure

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