We address the question of identifying the long-term (secular) stability regions in the semi-major axis-eccentricity projected phase space of the Sun-Jupiter planar circular restricted three-body problem in the domains i) below the curve of apsis equal to the planet’s orbital radius (ensuring protection from collisions) and ii) above that curve. This last domain contains several Jupiter’s crossing trajectories. We discuss the structure of the numerical stability map in the (a,e) plane in relation to manifold dynamics. We also present a closed-form perturbation theory for particles with non-crossing highly eccentric trajectories exterior to the planet’s trajectory. Starting with a multipole expansion of the barycentric Hamiltonian, our method carries out a sequence of normalizations by Lie series in closed-form and without relegation. We discuss the applicability of the method as a criterion for estimating the boundary of the domain of regular motion.

We have done N-body simulations with N up to 106, with the aim to determine whether fluctuations in the force field of a globular cluster are caused by nearby or distant encounters. We find that distant encounters are insignificant, in agreement with Chandrasekhar’s expectations, contrary to general opinion.

The origin of spiral structure in disks of galaxies remains an open question. One of the theories predicts that two-armed, grand design spiral arms originate from tidal interactions with another body. Using N-body simulations we find that a Milky Way-like galaxy can develop spiral arms due to tidal force from a cluster-size dark matter halo.

A major difficulty hampering the accuracy of UV/optical star formation rate tracers is the effect of interstellar dust, absorbing and scattering light produced by both young and old stellar populations (SPs). Although empirically calibrated corrections or energy balance SED fitting are often used for fast de-reddening of galaxy stellar emission, eventually only radiative transfer calculations can provide self-consistent predictions of galaxy model spectra, taking into account important factors such as galaxy inclination, different morphological components, non-local heating of the dust and scattered radiation. In addition, dust radiative transfer can be used to determine the fraction of monochromatic dust emission powered by either young or old SPs. This calculation needs to take into account the different response of the dust grains to the UV and optical radiation field, depending on the grain size and composition. We determined the dust heating fractions, on both global and local scales, for a high-resolution galaxy model by using our 3D ray-tracing dust radiative transfer code “DART-Ray”. We show the results obtained using this method and discuss the consequences for star formation rate indicators.

The star formation history of galaxies is a complex process usually considered to be stochastic in nature, for which we can only give average descriptions such as the color-density relation. In this work we follow star-forming gas particles in a hydrodynamical N-body simulation back in time in order to study their initial spatial configuration. By keeping record of the time when a gas particle started forming stars we can produce Lagrangian gas-star isochrone surfaces delineating the surfaces of accreting gas that begin producing stars at different times. These surfaces form a complex a network of filaments in Eulerian space from which galaxies accrete cold gas. Lagrangian accretion surfaces are closely packed inside dense regions, intersecting each other, and as a result galaxies inside proto-clusters stop accreting gas early, naturally explaining the color dependence on density. The process described here has a purely gravitational / geometrical origin, arguably operating at a more fundamental level than complex processes such as AGN and supernovae, and providing a conceptual origin for the color-density relation.

The Cosmic Web is a complicated highly-entangled geometrical object. Remarkably it has formed from practically Gaussian initial conditions, which may be regarded as the simplest departure from exactly uniform universe in purely deterministic mapping. The full complexity of the web is revealed neither in configuration no velocity spaces considered separately. It can be fully appreciated only in six-dimensional (6D) phase space. However, studies of the phase space is complicated by the fact that every projection of it on a three-dimensional (3D) space is multivalued and contained caustics. In addition phase space is not a metric space that complicates studies of geometry. We suggest to use Lagrangian submanifold i.e., x = x(q), where both x and q are 3D vectors instead of the phase space for studies the complexity of cosmic web in cosmological N-body dark matter simulations. Being fully equivalent in dynamical sense to the phase space it has an advantage of being a single valued and also metric space.

We present the current standing of an investigation into the structure of the Milky Way. We use smoothed particle hydrodynamics (SPH) to simulate the ISM gas in the Milky Way under the effect of a number of different gravitational potentials representing the spiral arms and nuclear bars, both fixed and time-dependent. The gas is subject to ISM cooling and chemistry, enabling us to track the CO and HI density. We use a 3D grid-based radiative transfer code to simulate the emission from the SPH output, allowing for the construction of synthetic longitude-velocity maps as viewed from the Earth. By comparing these maps with the observed emission in CO and HI from the Milky Way ([Dame et al. 2001, Kalberla et al. 2005]), we can infer the arm/bar geometry that provides a best fit to our Galaxy. By doing so we aim to answer key questions concerning the morphology of the Milky Way such as the number of the spiral arms, the pattern speeds of the bar(s) and arms, the pitch angle of the arms and shape of the bar(s).

The growing number of detected planets in binary star systems requires methods for a quick assessment of possible habitability of planets in such environments. We offer an analytic method to determine habitability of a terrestrial planet in binary star systems. In this context we give an answer to the most important question: Do the radiative and gravitational perturbations of a secondary influence the extent of the habitable zone (HZ)?

After we have defined the borders of the HZ, we will show the dynamical behaviour of a terrestrial planet in the HZ when adding a Jupiter to the system. In such a system the HZ shows signs of mean motion resonances and secular resonances, depending on the architecture of the planetary system.

Observations of about 60 binary star systems hosting exoplanets indicate the necessity of stability studies of planetary motion in such multi-stellar systems. For wide binary systems with separations between hundreds and thousands of AU, the results from single-star systems may be applicable but, in tight double stars systems, we have to take the stellar interactions into account which influences the planetary motion significantly.

This review discusses the different types of planetary motion in double stars and the stability of the planets for different binary configurations. An application to the most famous tight binary system (γ Cephei) is also shown. Finally, we analyze the habitability from the dynamical point of view in such systems, where we discuss the motion of terrestrial-like planets in the so-called habitable zone.

The evolution and amplification of large-scale magnetic fields in galaxies is investigated by means of high resolution simulations of interacting galaxies. The goal of our project is to consider in detail the role of gravitational interaction of galaxies for the fate of magnetic fields. Since the tidal interaction up to galaxy merging is a basic ingredient of cold-dark matter (CDM) structure formation models we think that our simulations will give important clues for the interplay of galactic dynamics and magnetic fields.

We discuss plans for a new joint effort between observers and theorists to understand the formation of the Milky Way halo back to the first epochs of chemical evolution. New models based on high-resolution N-body simulations coupled to simple models of Galactic chemical evolution show that surviving stars from the epoch of the first galaxies remain in the Milky Way today and should bear the nucleosynthetic imprint of the first stars. We investigate the key physical influences on the formation of stars in the first galaxies and how they appear today, including the relationship between cosmic reionization and surviving Milky Way stars. These models also provide a physically motivated picture of the formation of the Milky Ways “outer halo,” which has been identified from recent large samples of stars from SDSS. The next steps are to use these models to guide rigorous gas simulations of Milky Way formation, including its disk, and to gradually build up the fully detailed theoretical “Virtual Galaxy” that is demanded by the coming generation of massive Galactic stellar surveys.

We present a direct N-body simulation modeling the evolution of the old (7 Gyr) open cluster NGC 188. This is the first N-body open cluster simulation whose initial binary population is directly defined by observations of a specific open cluster: M35 (150 Myr). We compare the simulated color–magnitude diagram at 7 Gyr to that of NGC 188, and discuss the blue stragglers produced in the simulation. We compare the solar-type main-sequence binary period and eccentricity distributions of the simulation to detailed observations of similar binaries in NGC 188. These results demonstrate the importance of detailed observations in guiding N-body open cluster simulations. Finally, we discuss the implications of a few discrepancies between the NGC 188 model and observations and suggest a few methods for bringing N-body open cluster simulations into better agreement with observations.

We discuss the dynamical evolution of minor planetary bodies in the outer regions of planetary systems around the progenitors of DZ white dwarfs. We show that during the planetary-nebula phase of these stars, mass loss can lead to the expansion of all planetary bodies. The orbital eccentricity of the minor bodies, as relics of planetesimals, may be largely excited by the perturbation due to both gas drag effects and nearby gas giant planets. Some of these bodies migrate toward the host star, while others are scattered out of the planetary system. The former have modest probability of being captured by the sweeping secular resonances of giant planets, and induced to migrate toward the host star. When they venture close to their host stars, their orbits are tidally circularized so that they form compact disks where they may undergo further collisionally driven evolution. During the subsequent post main sequence evolution of their host stars, this process may provide an avenue which continually channels heavy elements onto the surface of the white dwarfs. We suggest that this scenario provides an explanation for the recently discovered Calcium line variation in G29-38.

Massive star clusters in the Magellanic Clouds are observed to follow a striking trend in size with age – older clusters exhibit a much greater spread in core radius than do younger clusters, which are generally compact. We present results from realistic N-body modelling of massive star clusters, aimed at investigating a dynamical origin for the radius-age trend. We find that stellar-mass black holes, formed as remnants of the most massive stars in a cluster, can constitute a dynamically important population. If retained, these objects rapidly form a dense core where interactions are common, resulting in the scattering of black holes into the cluster halo, and the ejection of black holes from the cluster. These processes heat the stellar component, resulting in prolonged core expansion of a magnitude matching the observations. Core expansion at early times does not result from the action of black holes, but can be reproduced by the effects of rapid mass-loss due to stellar evolution in a primordially mass segregated cluster.

We study the formation, growth, and co-evolution of single and multiple supermassive black holes (SMBHs) and compact objects like neutron stars, white dwarfs, and stellar mass black holes in galactic nuclei and star clusters, focusing on the role of stellar dynamics. In this paper we focus on one exemplary topic out of a wider range of work done, the study of orbital parameters of binary black holes in galactic nuclei (binding energy, eccentricity, relativistic coalescence) as a function of initial parameters. In some cases the classical evolution of black hole binaries in dense stellar systems drives them to surprisingly high eccentricities, which is very exciting for the emission of gravitational waves and relativistic orbit shrinkage. Such results are interesting to the emerging field of gravitational wave astronomy, in relation to a number of ground and space based instruments designed to measure gravitational waves from astrophysical sources (VIRGO, Geo600, LIGO, LISA). Our models self-consistently cover the entire range from Newtonian dynamics to the relativistic coalescence of SMBH binaries.

External destruction of protoplanetary discs acts mainly due to two mechanisms: gravitational drag by stellar encounters and evaporation by stellar winds and radiation. It is an important question whether any of these mechanisms is important in the stellar evolution process. We focus on the effect of stellar encounters and investigate if there are any observables that could trace this mechanism in young stellar clusters. An analysis of observational data of the Orion Nebula Cluster (ONC) and accompanying n-body simulations both provide evidence for encounters of star-disc systems in the ONC, eventually leading to substantial disk disruption.

About a fifth of the exoplanetary systems that have been discovered contain a so-called hot-Jupiter – a giant planet orbiting within 0.1 AU of the central star. Since these stars are typically of the F/G spectral type, the orbits of any terrestrial planets in their habitable zones at ~1 AU should be dynamically stable. However, because hot-Jupiters are thought to have formed in the outer regions of a protoplanetary disc, and to have then migrated through the terrestrial planet zone to their final location, it is uncertain whether terrestrial planets can actually grow and be retained in these systems. In this paper we review attempts to answer this question. Initial speculations, based on the assumption that migrating giant planets will clear planet-forming material from their swept zone, all concluded that hot-Jupiter systems should lack terrestrial planets. We show that this assumption may be incorrect, for when terrestrial planet formation and giant planet migration are simulated simultaneously, abundant solid material is predicted to remain from which terrestrial planet growth can resume.

We investigate the secular resonances for massless small bodies and Earth-like planets in several planetary systems. We further compare the results with those of Solar System. For example, in the GJ 876 planetary system, we show that the secular resonances ν1 and ν2 (respectively, resulting from the inner and outer giant planets) can excite the eccentricities of the Earth-like planets with orbits 0.21≤ a <0.50 AU and eject them out of the system in a short timescale. However, in a dynamical sense, the potential zones for the existence of Earth-like planets are in the area 0.50≤ a ≤1.00 AU, and there exist all stable orbits last up to 105 yr with low eccentricities. For other systems, e.g., 47 UMa, we also show that the Habitable Zones for Earth-like planets are related to both secular resonances and mean motion resonances in the systems.

More than 300 000 artificial debris particles with diameter larger than 1 cm are orbiting the Earth. The space debris population is similar to the asteroid belt, since it is subject to a process of high-velocity mutual collisions that affects the long-term evolution of its size distribution. The near–Earth space can be divided in three major regions where orbital debris is of concern: Low Earth Orbits (LEOs), below about 2000 km, Geosynchronous Orbits (GEOs), at an altitude of about 36000 km and the Medium Earth Orbits (MEOs) in between. The issues are in principle the same in the three regions, nevertheless they require different approaches and solutions. The space debris are composed by several different populations according to their source and their orbital region. A description of the nature and dynamics of the different populations in the low, medium and high orbital regimes is given. The impact risk posed by these debris is then briefly outlined.

The long term evolution of the whole debris population can be studied with computer models allowing the simulation of all the known source and sinks mechanisms. One of these codes is described and the evolution of the debris environment, over the next 100 years, under different traffic scenarios is shown, pointing out the possible measures to mitigate the growth of the orbital debris population.To search for other articles by the author(s) go to: http://adsabs.harvard.edu/abstract_service.html