Source author record

Eric I. Barnes

Eric I. Barnes appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA astro-ph.CO astro-ph astro-ph.IM cond-mat.stat-mech physics.atom-ph physics.gen-ph

Catalog footprint

What is connected

9works

7topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2016arXiv

Comparing Extended System Interactions with Motions in Softened Potentials

Using an $N$-body evolution code that does not rely on softened potentials, I have created a suite of unbound interacting cluster pair simulations. The motions of the centers of mass of the clusters have been tracked and compared to the trajectories of point masses interacting via one of four different softened potential prescriptions. I find that the relationship between the impact parameter of the cluster interaction and the point-mass softening length that best approximates each cluster's center-of-mass motion depends on the adopted prescription. In general, the range of allowed softening lengths grows roughly linearly with the impact parameter, but zero softening is acceptable in the majority of situations. In an $N$-body simulation that adopts a fixed softening length, such relationships lead to the possibility of two-body effects, like dynamical friction, being either larger or smaller than the corresponding cluster situation. Further consideration of more specific $N$-body situations leads estimating that a very small fraction of point-mass encounters experience two-body effects significantly different than those of equivalent clusters.

preprint2015arXiv

Entropy Production in Collisionless Systems. III. Results from Simulations

The equilibria formed by the self-gravitating, collisionless collapse of simple initial conditions have been investigated for decades. We present the results of our attempts to describe the equilibria formed in $N$-body simulations using thermodynamically-motivated models. Previous work has suggested that it is possible to define distribution functions for such systems that describe maximum entropy states. These distribution functions are used to create radial density and velocity distributions for comparison to those from simulations. A wide variety of $N$-body code conditions are used to reduce the chance that results are biased by numerical issues. We find that a subset of initial conditions studied lead to equilibria that can be accurately described by these models, and that direct calculation of the entropy shows maximum values being achieved.

preprint2013arXiv

Dynamics of One-dimensional Self-gravitating Systems Using Hermite-Legendre Polynomials

The current paradigm for understanding galaxy formation in the universe depends on the existence of self-gravitating collisionless dark matter. Modeling such dark matter systems has been a major focus of astrophysicists, with much of that effort directed at computational techniques. Not surprisingly, a comprehensive understanding of the evolution of these self-gravitating systems still eludes us, since it involves the collective nonlinear dynamics of many-particle systems interacting via long-range forces described by the Vlasov equation. As a step towards developing a clearer picture of collisionless self-gravitating relaxation, we analyze the linearized dynamics of isolated one-dimensional systems near thermal equilibrium by expanding their phase space distribution functions f(x,v) in terms of Hermite functions in the velocity variable, and Legendre functions involving the position variable. This approach produces a picture of phase-space evolution in terms of expansion coefficients, rather than spatial and velocity variables. We obtain equations of motion for the expansion coefficients for both test-particle distributions and self-gravitating linear perturbations of thermal equilibrium. N-body simulations of perturbed equilibria are performed and found to be in excellent agreement with the expansion coefficient approach over a time duration that depends on the size of the expansion series used.

preprint2012arXiv

Collisionless dynamics in Globular Clusters

Since globular clusters (GCs) are old, low-N systems their dynamics is widely believed to be fully dominated by collisional two-body processes, and their surface brightness profiles are fit by King models. However, for many GCs, especially those with HST-resolved central regions, and `extra-tidal' features, King models provide poor fits. We suggest that this is partly because collisionless dynamics is also important and contribute to shaping the cluster properties. We show using time-scale and length-scale arguments that except for the very centers of clusters, collisionless dynamics should be more important than collisional. We then fit 38 GCs analyzed by Noyola and Gebhardt (2006) with (collisional) King and (collisionless) DARKexp models over the full available radial range, and find that the latter provide a better fit to 29 GCs; for six of these the fit is at least ~5x better in term of rms. DARKexp models are theoretically derived maximum entropy equilibrium states of self-gravitating collisionless systems and have already been shown to fit the results of dark matter N-body simulations. (We do not attempt fits with ad hoc fitting functions.)

preprint2012arXiv

Entropy Production in Collisionless Systems. II. Arbitrary Phase-Space Occupation Numbers

We present an analysis of two thermodynamic techniques for determining equilibria of self-gravitating systems. One is the Lynden-Bell entropy maximization analysis that introduced violent relaxation. Since we do not use the Stirling approximation which is invalid at small occupation numbers, our systems have finite mass, unlike Lynden-Bell's isothermal spheres. (Instead of Stirling, we utilize a very accurate smooth approximation for $\ln{x!}$.) The second analysis extends entropy production extremization to self-gravitating systems, also without the use of the Stirling approximation. In addition to the Lynden-Bell (LB) statistical family characterized by the exclusion principle in phase-space, and designed to treat collisionless systems, we also apply the two approaches to the Maxwell-Boltzmann (MB) families, which have no exclusion principle and hence represent collisional systems. We implicitly assume that all of the phase-space is equally accessible. We derive entropy production expressions for both families, and give the extremum conditions for entropy production. Surprisingly, our analysis indicates that extremizing entropy production rate results in systems that have maximum entropy, in both LB and MB statistics. In other words, both thermodynamic approaches lead to the same equilibrium structures.

preprint2011arXiv

Entropy Production in Collisionless Systems. I. Large Phase-Space Occupation Numbers

Certain thermal non-equilibrium situations, outside of the astrophysical realm, suggest that entropy production extrema, instead of entropy extrema, are related to stationary states. In an effort to better understand the evolution of collisionless self-gravitating systems, we investigate the role of entropy production and develop expressions for the entropy production rate in two particular statistical families that describe self-gravitating systems. From these entropy production descriptions, we derive the requirements for extremizing the entropy production rate in terms of specific forms for the relaxation function in the Boltzmann equation. We discuss some implications of these relaxation functions and point to future work that will further test this novel thermodynamic viewpoint of collisionless relaxation.

preprint2009arXiv

The Radial Orbit Instability in Collisionless N-Body Simulations

Using a suite of self-gravitating, collisionless N-body models, we systematically explore a parameter space relevant to the onset and behavior of the radial orbit instability (ROI), whose strength is measured by the systemic axis ratios of the models. We show that a combination of two initial conditions, namely the velocity anisotropy and the virial ratio, determines whether a system will undergo ROI and exactly how triaxial the system will become. A third initial condition, the radial shape of the density profile, plays a smaller, but noticeable role. Regarding the dynamical development of the ROI, the instability a) begins after systems collapse to their most compact configuration and b) evolves fastest when a majority of the particles have radially anisotropic orbits while there is a lack of centrally-concentrated isotropic orbits. We argue that this is further evidence that self-reinforcing torques are the key to the onset of the ROI. Our findings support the idea that a separate orbit instability plays a role in halting the ROI.

preprint2006arXiv

Density Profiles of Collisionless Equilibria. II. Anisotropic Spherical Systems

It has long been realized that dark matter halos formed in cosmological N-body simulations are characterized by density profiles rho(r) that, when suitably scaled, have similar shapes. Additionally, combining the density and velocity dispersion profiles sigma(r), each of which have decidedly nonpower-law shapes, leads to quantity rho/(sigma^3) that is a power-law in radius over 3 orders of magnitude in radius. Halos' velocity anisotropy profiles beta(r) vary from isotropic near the centers of halos to quite radially anisotropic near the virial radius. Finally, there appears to be a nearly linear correlation between beta and the logarithmic density slope gamma for a wide variety of halos. This work is part of a continuing investigation of the above interrelationships and their origins using analytical and semi-analytical techniques. Our finding suggest that the nearly linear beta--gamma relationship is not just another expression of scale-free rho/(sigma^3) behavior. We also note that simultaneously reproducing density and anisotropy profiles like those found in simulations requires beta(r) and gamma(r) to have similar shapes, leading to nearly linear beta--gamma correlations. This work suggests that the beta--gamma and power-law rho/(sigma^3) relations have distinct physical origins.

preprint2001arXiv

Alternative Fourier Expansions for Inverse Square Law Forces

Few-body problems involving Coulomb or gravitational interactions between pairs of particles, whether in classical or quantum physics, are generally handled through a standard multipole expansion of the two-body potentials. We discuss an alternative based on a compact, cylindrical Green's function expansion that should have wide applicability throughout physics. Two-electron "direct" and "exchange" integrals in many-electron quantum systems are evaluated to illustrate the procedure which is more compact than the standard one using Wigner coefficients and Slater integrals.