Abstract
We describe an algorithm for constructing N-body realisations of equilibrium stellar systems. The algorithm complements existing orbit-based modelling techniques using linear programming or other optimization algorithms. The equilibria are constructed by integrating an N-body system while slowly adjusting the masses of the particles until the time-averaged density field and other observables converge to a prescribed value. The procedure can be arranged to maximise a linear combination of the entropy of the system and the \\chi^2 statistic for the observables. The equilibria so produced may be useful as initial conditions for N-body simulations or for modelling observations of individual galaxies.
Keywords
Affiliated Institutions
Related Publications
Convergent Algorithm for Sensory Receptive Field Development
An unsupervised developmental algorithm for linear maps is derived which reduces the pixel-entropy (using the measure introduced in previous work) at every update and thus remov...
From NBODY1 to NBODY6: The Growth of an Industry
I review the development of direct N‐body codes at Cambridge over nearly 40 years, highlighting the main stepping stones. The first code (NBODY1) was based on the simple concept...
Linear Model Selection by Cross-Validation
Abstract We consider the problem of selecting a model having the best predictive ability among a class of linear models. The popular leave-one-out cross-validation method, which...
Simulations of dissipative galaxy formation in hierarchically clustering universes - I: Tests of the code
We present tests of a code designed to simulate the evolution of self- gravitating fluids in three dimensions. The code is based on the smoothed-particle hydrodynamics (SPH) tec...
Publication Info
- Year
- 1996
- Type
- article
- Volume
- 282
- Issue
- 1
- Pages
- 223-233
- Citations
- 178
- Access
- Closed
External Links
Social Impact
Social media, news, blog, policy document mentions
Citation Metrics
Cite This
Identifiers
- DOI
- 10.1093/mnras/282.1.223