Abstract
We describe an algorithm for constructing N-body realizations 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 maximize a linear combination of the entropy of the system and the χ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.
Sign-in/Register to access full text options 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
