An Advanced N body Model for Interacting Multiple Stellar Systems
Abstract
We construct an advanced model for interacting multiple stellar systems in which we compute all trajectories with a numerical N body integrator, namely the Bulirsch–Stoer from the SWIFT package. We can then derive various observables: astrometric positions, radial velocities, minima timings (TTVs), eclipse durations, interferometric visibilities, closure phases, synthetic spectra, spectral energy distribution, and even complete light curves. We use a modified version of the Wilson–Devinney code for the latter, in which the instantaneous true phase and inclination of the eclipsing binary are governed by the N body integration. If all of these types of observations are at one’s disposal, a joint χ {sup 2} metric and an optimization algorithm (a simplex or simulated annealing) allow one to search for a global minimum and construct very robust models of stellar systems. At the same time, our N body model is free from artifacts that may arise if mutual gravitational interactions among all components are not selfconsistently accounted for. Finally, we present a number of examples showing dynamical effects that can be studied with our code and we discuss how systematic errors may affect the results (and how to prevent this from happening).
 Authors:
 Astronomical Institute of the Charles University, Faculty of Mathematics and Physics, V Holešovičkách 2, CZ18000 Praha 8 (Czech Republic)
 Publication Date:
 OSTI Identifier:
 22661119
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Astrophysical Journal, Supplement Series; Journal Volume: 230; Journal Issue: 2; Other Information: Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 79 ASTROPHYSICS, COSMOLOGY AND ASTRONOMY; ALGORITHMS; DIAGRAMS; ECLIPSE; ENERGY SPECTRA; GRAVITATIONAL INTERACTIONS; INCLINATION; METRICS; OPTIMIZATION; RADIAL VELOCITY; STAR CLUSTERS; STARS; VISIBILITY; VISIBLE RADIATION
Citation Formats
Brož, Miroslav. An Advanced N body Model for Interacting Multiple Stellar Systems. United States: N. p., 2017.
Web. doi:10.3847/15384365/AA7207.
Brož, Miroslav. An Advanced N body Model for Interacting Multiple Stellar Systems. United States. doi:10.3847/15384365/AA7207.
Brož, Miroslav. Thu .
"An Advanced N body Model for Interacting Multiple Stellar Systems". United States.
doi:10.3847/15384365/AA7207.
@article{osti_22661119,
title = {An Advanced N body Model for Interacting Multiple Stellar Systems},
author = {Brož, Miroslav},
abstractNote = {We construct an advanced model for interacting multiple stellar systems in which we compute all trajectories with a numerical N body integrator, namely the Bulirsch–Stoer from the SWIFT package. We can then derive various observables: astrometric positions, radial velocities, minima timings (TTVs), eclipse durations, interferometric visibilities, closure phases, synthetic spectra, spectral energy distribution, and even complete light curves. We use a modified version of the Wilson–Devinney code for the latter, in which the instantaneous true phase and inclination of the eclipsing binary are governed by the N body integration. If all of these types of observations are at one’s disposal, a joint χ {sup 2} metric and an optimization algorithm (a simplex or simulated annealing) allow one to search for a global minimum and construct very robust models of stellar systems. At the same time, our N body model is free from artifacts that may arise if mutual gravitational interactions among all components are not selfconsistently accounted for. Finally, we present a number of examples showing dynamical effects that can be studied with our code and we discuss how systematic errors may affect the results (and how to prevent this from happening).},
doi = {10.3847/15384365/AA7207},
journal = {Astrophysical Journal, Supplement Series},
number = 2,
volume = 230,
place = {United States},
year = {Thu Jun 01 00:00:00 EDT 2017},
month = {Thu Jun 01 00:00:00 EDT 2017}
}

Relativistic particle dynamics for an Nbody interacting system
Claims in a recent Comment by Liou to find new forms for the threebody interactions to order 1/c/sup 2/ in an Nbody relativistic interacting system are shown to arise from a misunderstanding. The proposed forms are encompassed by the solution to this problem obtained previously by the present authors. Several other clarifications are also made. (AIP) 
MOMENT OF INERTIA OF INTERACTING MANYBODY FERMION SYSTEMS
It is shown that the moment of inertia of a noninteracting manybody fermion system moving under periodic boundary conditions has the classical or rigid value when calculated on the "cranking" model of Inglis. By investigating the analogous "pushing" case for the inertial mass we show that the rigid moment can be associated with rigid rotation in spite of apparent surface currents. The effect of particleparticle forces is investigated in the lowest order of perturbation theory. The terms corresponding to a level shift or effective mass are just compensated by other terms and there is no change in the moment ofmore » 
Spectroscopy of Collective Excitations in Interacting LowDimensional ManyBody Systems Using Quench Dynamics
We study the problem of rapid change of the interaction parameter (quench) in a manybody lowdimensional system. It is shown that, measuring the correlation functions after the quench, the information about a spectrum of collective excitations in a system can be obtained. This observation is supported by analysis of several integrable models and we argue that it is valid for nonintegrable models as well. Our conclusions are supplemented by performing exact numerical simulations on finite systems. We propose that measuring the power spectrum in a dynamically split 1D BoseEinsten condensate into two coupled condensates can be used as an experimentalmore » 
Nearestneighbor distributions and tunneling splittings in interacting manybody twolevel boson systems
We study the nearestneighbor distributions of the kbody embedded ensembles of random matrices for n bosons distributed over twodegenerate singleparticle states. This ensemble, as a function of k, displays a transition from harmonicoscillator behavior (k=1) to randommatrixtype behavior (k=n). We show that a large and robust quasidegeneracy is present for a wide interval of values of k when the ensemble is timereversal invariant. These quasidegenerate levels are Shnirelman doublets which appear due to the integrability and timereversal invariance of the underlying classical systems. We present results related to the frequency in the spectrum of these degenerate levels in terms ofmore »