Our Galaxy and the nearby Andromeda Galaxy (M31) form a bound system, even though the relative velocity vector of M31 is currently not well constrained. Their orbital motion is highly dependent on the initial conditions, but all the reliable scenarios imply a first close approach in the next 3-5 Gyrs. In our study, we simulate this interaction via direct N-body integration, using the HiGPUs code. Our aim is to investigate the dependence of the time of the merger on the physical and dynamical properties of the system. Finally, we study the dynamical evolution of the two Supermassive Black Holes placed in the two galactic centers, with the future aim to achieve a proper resolution to follow their motion until they form a tight binary system.

According to our current knowledge about physical and dynamical properties of the Milky Way-M31 system, it seems likely that these two galaxies will collide and eventually merge in a time very sensitive to initial conditions. Using the HiGPUs code, we have performed several numerical simulations to study the dynamics of the system, trying to define the role of indeterminacy in the present day observed relative velocities of the two galaxies and the time of the merger. At the same time, we have followed the dynamics of the two massive black holes sitting in the galactic centers, to check (within the space and time resolution limits of our simulation) their relative motion upon the completion of the galaxies merger process.

Only a small number of galactic open clusters survives for longer than few hundred million years. Longer lifetimes are routinely explained in term of larger initial masses, particularly quiet orbits, and off-plane birth-places. We derive in this work the actual mass of NGC 4337, one of the few open clusters in the Milky Way inner disk that managed to survive for about 1.5 Gyr. We derive its mass in two different ways. First, we exploit an unpublished photometric data set in the UBVI passbands to estimate -using star counts- the cluster luminosity profile, and luminosity and mass function, and hence its actual mass both from the luminosity profile and from the mass function.This data-set is also used to infer crucial cluster parameters, as the cluster half-mass radius and distance. Second, we make use of a large survey of cluster star radial velocities to derive dynamical estimates for the cluster mass. Under the assumption of virial equilibrium and neglecting the external gravitational field leads to values for the mass significantly larger than those obtained by mean of observed density distribution or with the mass function but still marginally compatible with the inferred values of the invisible mass in form of both low mass stars or remnants of high mass stars in the cluster. Finally, we derive the cluster initial mass by computing the mass loss experienced by the cluster during its lifetime, and adopting the various estimates of the actual mass.

Horizontal Branch stars belong to an advanced stage in the evolution of the oldest stellar galactic population, occurring either as field halo stars or grouped in globular clusters. The discovery of multiple populations in these clusters, that were previously believed to have single populations gave rise to the currently accepted theory that the hottest horizontal branch members (the blue hook stars, which had late helium-core flash ignition, followed by deep mixing) are the progeny of a helium-rich "second generation" of stars. It is not known why such a supposedly rare event (a late flash followed by mixing) is so common that the blue hook of {\omega} Cen contains \sim 30% of horizontal branch stars 10 , or why the blue hook luminosity range in this massive cluster cannot be reproduced by models. Here we report that the presence of helium core masses up to \sim 0.04 solar masses larger than the core mass resulting from evolution is required to solve the luminosity range problem. We model this by taking into account the dispersion in rotation rates achieved by the progenitors, whose premain sequence accretion disc suffered an early disruption in the dense environment of the cluster's central regions where second-generation stars form. Rotation may also account for frequent late-flash-mixing events in massive globular clusters.

The application of high end computing to astrophysical problems, mainly in the galactic environment, is under development since many years at the Dep. of Physics of Sapienza Univ. of Roma. The main scientific topic is the physics of self gravitating systems, whose specific subtopics are: i) celestial mechanics and interplanetary probe transfers in the solar system; ii) dynamics of globular clusters and of globular cluster systems in their parent galaxies; iii) nuclear clusters formation and evolution; iv) massive black hole formation and evolution; v) young star cluster early evolution. In this poster we describe the software and hardware computational resources available in our group and how we are developing both software and hardware to reach the scientific aims above itemized.

In the last decade, observational studies have shown the existence of tidal streams in the outer part of many galactic globular clusters. The most striking examples of clusters with well defined tidal tails are represented by Palomar 5 and NGC 5466 (both observed in the framework of the Sloan Digital Sky Survey), which show structures elongated for 4 kpc and 1 kpc in length, respectively. Unfortunately, most of the observational studies about globular clusters (GCs) do not cover such a large field of the sky as the SDSS does. In this framework, by mean of a parallel, adaptive tree-code, we performed detailed N-body simulations of GCs moving in a realistic three-components (bulge, disk and halo) Milky Way potential, in order to clarify whether and to what extent tails in the clusters outer regions (few tidal radii) are tracers of the local orbits and, also, if some kind of correlation exists among the cluster orbital phase and the orientation of such streams.

We present the results of detailed N-body simulations regarding the interaction of four massive globular clusters in the central region of a triaxial galaxy. The systems undergo a full merging event, producing a sort of 'Super Star Cluster' (SSC) whose features are close to those of a superposition of the individual initial mergers. In contrast with other similar simulations, the resulting SSC structural parameters are located along the observed scaling relations of globular clusters. These findings seem to support the idea that a massive SSC may have formed in early phases of the mother galaxy evolution and contributed to the growth of a massive nucleus.

The main topic of this paper is the investigation of the modes of interaction of globular clusters (GCs) moving in the inner part of a galaxy. This is tackled by means of high-resolution N-body simulations, whose first results are presented in this article. Our simulations dealt with primordial very massive (order of 10^7 solar masses) GCs that were able to decay, because of dynamical friction, into the inner regions of triaxial galaxies on a time much shorter than their internal relaxation time. To check the disruptive role of both tidal forces and GC-GC collisions, we maximised the tidal interaction considering GCs on quasi-radial orbits. The available CPU resources allowed us to follow 8 oscillations of the GCs along their orbits and the main findings are: i) clusters with an initial high enough King concentration parameter (c>=1.2), preserve up to 50% of their initial mass; ii) the inner density distribution of the survived clusters keep a King model profile; iii) GC-GC collisions have a negligible effect with respect to that caused by the passage through the galactic center; iv) the orbital energy dissipation due to the tidal interaction is of the same order of that caused by dynamical friction; v) complex sub-structures like "ripples" and "clumps" formed, as observed around real clusters. These findings support the validity of the hypothesis of merging of GCs in the galactic central region, with modes that deserve further careful investigations.

Astrocomp is a joint project, developed by the INAF-Astrophysical Observatory of Catania, University of Roma La Sapienza and Enea. The project has the goal of providing the scientific community of a web-based user-friendly interface which allows running parallel codes on a set of high-performance computing (HPC) resources, without any need for specific knowledge about parallel programming and Operating Systems commands. Astrocomp provides, also, computing time on a set of parallel computing systems, available to the authorized user. At present, the portal makes a few codes available, among which: FLY, a cosmological code for studying three-dimensional collisionless self-gravitating systems with periodic boundary conditions; ATD, a parallel tree-code for the simulation of the dynamics of boundary-free collisional and collisionless self-gravitating systems and MARA, a code for stellar light curves analysis. Other codes are going to be added to the portal.

We present some results of numerical simulations of a globular cluster orbiting in the central region of a triaxial galaxy on a set of 'loop' orbits. Tails start forming after about a quarter of the globular cluster orbital period and develop, in most cases, along the cluster orbit, showing clumpy substructures as observed, for example, in Palomar 5. If completely detectable, clumps can contain about 7,000 solar masses each, i.e. about 10% of the cluster mass at that epoch. The morphology of tails and clumps and the kinematical properties of stars in the tails are studied and compared with available observational data. Our finding is that the stellar velocity dispersion tends to level off at large radii, in agreement to that found for M15 and Omega Centauri.

We study gravitational clustering of mass points in three dimensions with random initial positions and periodic boundary conditions (no expansion) by numerical simulations. Correlation properties are well defined in the system and a sort of thermodynamic limit can be defined for the transient regime of cluste ring. Structure formation proceeds along two paths: (i) fluid-like evolution of density perturbations at large scales and (ii) shift of the granular (non fluid) properties from small to large scales. The latter mechanism finally dominates at all scales and it is responsible for the self-similar characteristics of the clustering.