Statistical mechanics of collisionless orbits. II. Structure of halos

In this paper, we present the density, \rho, velocity dispersion, \sigma, and \rho/\sigma^3 profiles of isotropic systems which have the energy distribution, N(E)\propto[\exp(\phi_0-E)-1], derived in Paper I. This distribution, dubbed DARKexp, is the most probable final state of a collisionless self-gravitating system, which is relaxed in terms of particle energies, but not necessarily in terms of angular momentum. We compare the DARKexp predictions with the results obtained using the extended secondary infall model (ESIM). The ESIM numerical scheme is optimally suited for the purpose because (1) it relaxes only through energy redistribution, leaving shell/particle angular momenta unaltered, and (2) being a shell code with radially increasing shell thickness it has very good mass resolution in the inner halo, where the various theoretical treatments give different predictions. The ESIM halo properties, and especially their energy distributions, are very well fit by DARKexp, implying that the techniques of statistical mechanics can be used to explain the structure of relaxed self-gravitating systems.Comment: 17 pages, 8 figure

Statistical mechanics of collisionless orbits. I. Origin of central cusps in dark-matter halos

We present an equilibrium statistical mechanical theory of collisionless self-gravitational systems with isotropic velocity distributions. Compared to existing standard theories, we introduce two changes: (1) the number of possible microstates is computed in energy (orbit) space rather than phase space and (2) low occupation numbers are treated more appropriately than using Stirling's approximation. Combined, the two modifications predict that the relaxed parts of collisionless self-gravitating systems, such as dark-matter halos, have a differential energy distribution N(E) ~ [exp(phi_0 - E) - 1], dubbed "DARKexp". Such systems have central power-law density cusps rho(r) ~ r^-1, which suggests a statistical mechanical origin of cusps in simulated dark-matter halos.Comment: 6 pages, 2 figure

Statistical mechanics of collisionless orbits. III. Comparison with N-body simulations

We compare the DARKexp differential energy distribution, N(E) \propto \exp(\phi_0-E)-1, obtained from statistical mechanical considerations, to the results of N-body simulations of dark matter halos. We first demonstrate that if DARKexp halos had anisotropic velocity distributions similar to those of N-body simulated halos, their density and energy distributions could not be distinguished from those of isotropic DARKexp halos. We next carry out the comparison in two ways, using (1) the actual energy distribution extracted from simulations, and (2) N-body fitting formula for the density distribution as well as N(E) computed from the density using the isotropic Eddington formula. Both the methods independently agree that DARKexp N(E) with \phi_0\approx 4-5 is an excellent match to N-body N(E). Our results suggest (but do not prove) that statistical mechanical principles of maximum entropy can be used to explain the equilibrated final product of N-body simulations.Comment: 17 pages, 7 figures; ApJ, in pres

Dynamics of merging: Post-merger mixing and relaxation of an Illustris galaxy

During the merger of two galaxies, the resulting system undergoes violent relaxation and seeks stable equilibrium. However, the details of this evolution are not fully understood. Using Illustris simulation, we probe two physically related processes, mixing and relaxation. Though the two are driven by the same dynamics---global time-varying potential for the energy, and torques caused by asymmetries for angular momentum---we measure them differently. We define mixing as the redistribution of energy and angular momentum between particles of the two merging galaxies. We assess the degree of mixing as the difference between the shapes of their N(E)s, and their N(L^2)s. We find that the difference is decreasing with time, indicating mixing. To measure relaxation, we compare N(E) of the newly merged system to N(E) of a theoretical prediction for relaxed collisionless systems, DARKexp, and witness the system becoming more relaxed, in the sense that N(E) approaches DARKexp N(E). Because the dynamics driving mixing and relaxation are the same, the timescale is similar for both. We measure two sequential timescales: a rapid, 1 Gyr phase after the initial merger, during which the difference in N(E) of the two merging halos decreases by ~80%, followed by a slow phase, when the difference decreases by ~50% over ~8.5 Gyrs. This is a direct measurement of the relaxation timescale. Our work also draws attention to the fact that when a galaxy has reached Jeans equilibrium it may not yet have reached a fully relaxed state given by DARKexp, in that it retains information about its past history. This manifests itself most strongly in stars being centrally concentrated. We argue that it is particularly difficult for stars, and other tightly bound particles, to mix because they have less time to be influenced by the fluctuating potential, even across multiple merger events.Comment: accepted for publication in JCA

Ubiquity of density slope oscillations in the central regions of galaxy and cluster-sized systems

One usually thinks of a radial density profile as having a monotonically changing logarithmic slope, such as in NFW or Einasto profiles. However, in two different classes of commonly used systems, this is often not the case. These classes exhibit non-monotonic changes in their density profile slopes which we call oscillations for short. We analyze these two unrelated classes separately. Class 1 consists of systems that have density oscillations and that are defined through their distribution function $f(E)$, or differential energy distribution $N(E)$, such as isothermal spheres, King profiles, or DARKexp, a theoretically derived model for relaxed collisionless systems. Systems defined through $f(E)$ or $N(E)$ generally have density slope oscillations. Class 1 system oscillations can be found at small, intermediate, or large radii but we focus on a limited set of Class 1 systems that have oscillations in the central regions, usually at $\log(r/r_{-2})\lesssim -2$, where $r_{-2}$ is the largest radius where $d\log(\rho)/d\log(r)=-2$. We show that the shape of their $N(E)$ can roughly predict the amplitude of oscillations. Class 2 systems which are a product of dynamical evolution, consist of observed and simulated galaxies and clusters, and pure dark matter halos. Oscillations in the density profile slope seem pervasive in the central regions of Class 2 systems. We argue that in these systems, slope oscillations are an indication that a system is not fully relaxed. We show that these oscillations can be reproduced by small modifications to $N(E)$ of DARKexp. These affect a small fraction of systems' mass and are confined to $\log(r/r_{-2})\lesssim 0$. The size of these modifications serves as a potential diagnostic for quantifying how far a system is from being relaxed.Comment: accepted by the Journal of Cosmology and Astroparticle Physics (JCAP

Non-universality of dark-matter halos: cusps, cores, and the central potential

Dark-matter halos grown in cosmological simulations appear to have central NFW-like density cusps with mean values of $d\log\rho/d\log r \approx -1$, and some dispersion, which is generally parametrized by the varying index $\alpha$ in the Einasto density profile fitting function. Non-universality in profile shapes is also seen in observed galaxy clusters and possibly dwarf galaxies. Here we show that non-universality, at any given mass scale, is an intrinsic property of DARKexp, a theoretically derived model for collisionless self-gravitating systems. We demonstrate that DARKexp - which has only one shape parameter, $\phi_0$ - fits the dispersion in profile shapes of massive simulated halos as well as observed clusters very well. DARKexp also allows for cored dark-matter profiles, such as those found for dwarf spheroidal galaxies. We provide approximate analytical relations between DARKexp $\phi_0$, Einasto $\alpha$, or the central logarithmic slope in the Dehnen-Tremaine analytical $\gamma$-models. The range in halo parameters reflects a substantial variation in the binding energies per unit mass of dark-matter halos.Comment: ApJ, in press, 10 pages, 7 figure

Calibration of the Fundamental Plane Zero-point in the Leo-I Group and an Estimate of the Hubble Constant

We derive new effective radii and total magnitudes for 5 E and S0 galaxies in the Leo-I group from wide-field CCD images. These are used in conjunction with recent literature velocity data to construct the fundamental plane (FP) of the Leo-I group. The rms scatter that we find is only 6 % in distance. The zero point of this relation provides a calibration of the FP as a distance indicator and directly determines the angular diameter distance ratio between the Leo-I group and more distant clusters. In the language of Jerjen and Tammann (1993) we determine a cosmic velocity of the Leo-I group of 757+-68 km/s relative to the Coma cluster, or 796+-57 km/s relative to a frame of 9 clusters. Combining this velocity with the Cepheid distance to M96, a member of Leo-I, we find the Hubble constant to be H_0=67+-8 km/s/Mpc or H_0=70+-7 km/s/Mpc for each case. The distance we obtain for the Coma cluster itself (108+-12 Mpc) is in good agreement with a number of other recent estimates.Comment: 19 pages, LaTeX, includes aaspp4.sty and 3 eps figures. To appear in ApJ. Also available at http://www.nordita.dk/~jens/preprints.htm

Effective interactions and shell model studies of heavy tin isotopes

We present results from large-scale shell-model calculations of even and odd tin isotopes from 134Sn to 142}Sn with a shell-model space defined by the 1f7/2,2p3/2,0h9/2,2p1/2,1f5/2,0i13/2 single-particle orbits. An effective two-body interaction based on modern nucleon-nucleon interactions is employed. The shell-model results are in turn analyzed for their pairing content using a generalized seniority approach. Our results indicate that a pairing-model picture captures a great deal of the structure and the correlations of the lowest lying states for even and odd isotopes.Comment: 7 pages, revtex latex style, submitted to PR

Effective Interaction Techniques for the Gamow Shell Model

We apply a contour deformation technique in momentum space to the newly developed Gamow shell model, and study the drip-line nuclei 5He, 6He and 7He. A major problem in Gamow shell-model studies of nuclear many-body systems is the increasing dimensionality of many-body configurations due to the large number of resonant and complex continuum states necessary to reproduce bound and resonant state energies. We address this problem using two different effective operator approaches generalized to the complex momentum plane. These are the Lee-Suzuki similarity transformation method for complex interactions and the multi-reference perturbation theory method. The combination of these two approaches results in a large truncation of the relevant configurations compared with direct diagonalization. This offers interesting perspectives for studies of weakly bound systems.Comment: 18 pages, 17 figs, Revtex