No Arabic abstract
We present a novel and flexible tensor approach to computing the effect of a time-dependent tidal field acting on a stellar system. The tidal forces are recovered from the tensor by polynomial interpolation in time. The method has been implemented in a direct-summation stellar dynamics integrator (NBODY6) and test-proved through a set of reference calculations: heating, dissolution time and structural evolution of model star clusters are all recovered accurately. The tensor method is applicable to arbitrary configurations, including the important situation where the background potential is a strong function of time. This opens up new perspectives in stellar population studies reaching to the formation epoch of the host galaxy or galaxy cluster, as well as for star-burst events taking place during the merger of large galaxies. A pilot application to a star cluster in the merging galaxies NGC 4038/39 (the Antennae) is presented.
We present a method to couple N-body star cluster simulations to a cosmological tidal field, using the Astrophysical Multipurpose Software Environment. We apply this method to star clusters embedded in the CosmoGrid dark matter-only LambdaCDM simulation. Our star clusters are born at z = 10 (corresponding to an age of the Universe of about 500 Myr) by selecting a dark matter particle and initializing a star cluster with 32,000 stars on its location. We then follow the dynamical evolution of the star cluster within the cosmological environment. We compare the evolution of star clusters in two Milky-Way size haloes with a different accretion history. The mass loss of the star clusters is continuous irrespective of the tidal history of the host halo, but major merger events tend to increase the rate of mass loss. From the selected two dark matter haloes, the halo that experienced the larger number of mergers tends to drive a smaller mass loss rate from the embedded star clusters, even though the final masses of both haloes are similar. We identify two families of star clusters: native clusters, which become part of the main halo before its final major merger event, and the immigrant clusters, which are accreted upon or after this event; native clusters tend to evaporate more quickly than immigrant clusters. Accounting for the evolution of the dark matter halo causes immigrant star clusters to retain more mass than when the z=0 tidal field is taken as a static potential. The reason for this is the weaker tidal field experienced by immigrant star clusters before merging with the larger dark matter halo.
Based on recent findings of a formation mechanism of substructure in tidal tails by Kuepper, Macleod & Heggie (2008) we investigate a more comprehensive set of N-body models of star clusters on orbits about a Milky-Way-like potential. We find that the predicted epicyclic overdensities arise in any tidal tail no matter which orbit the cluster follows as long as the cluster lives long enough for the overdensities to build up. The distance of the overdensities along the tidal tail from the cluster centre depends for circular orbits only on the mass of the cluster and the strength of the tidal field, and therefore decreases monotonically with time, while for eccentric orbits the orbital motion influences the distance, causing a periodic compression and stretching of the tails and making the distance oscillate with time. We provide an approximation for estimating the distance of the overdensities in this case. We describe an additional type of overdensity which arises in extended tidal tails of clusters on eccentric orbits, when the acceleration of the tidal field on the stellar stream is no longer homogeneous. Moreover, we conclude that a pericentre passage or a disk shock is not the direct origin of an overdensity within a tidal tail. Escape due to such tidal perturbations does not take place immediately after the perturbation but is rather delayed and spread over the orbit of the cluster. All observable overdensities are therefore of the mentioned two types. In particular, we note that substructured tidal tails do not imply the existence of dark-matter sub-structures in the haloes of galaxies.
In interacting galaxies, strong tidal forces disturb the global morphology of the progenitors and give birth to the long stellar, gaseous and dusty tails often observed. In addition to this destructive effect, tidal forces can morph into a transient, protective setting called compressive mode. Such modes then shelter the matter in their midst by increasing its gravitational binding energy. This thesis focuses on the study of this poorly known regime by quantifying its properties thanks to numerical and analytical tools applied to a spectacular merging system of two galaxies, commonly known as the Antennae galaxies. N-body simulations of this pair yield compressive modes in the regions where observations reveal a burst of star formation. Furthermore, characteristic time- and energy scales of these modes match well those of self-gravitating substructures such as star clusters and tidal dwarf galaxies. These results suggest that the compressive modes of tidal fields plays an important role in the formation and evolution of young clusters, at least in a statistical sense, over a lapse of ~10 million years. Preliminary results from simulations of stellar associations highlight the importance of embedding the clusters in the evolving background galaxies to account precisely for their morphology and internal evolution.
Until now it has been impossible to observationally measure how star cluster scale height evolves beyond 1Gyr as only small samples have been available. Here we establish a novel method to determine the scale height of a cluster sample using modelled distributions and Kolmogorov-Smirnov tests. This allows us to determine the scale height with a 25% accuracy for samples of 38 clusters or more. We apply our method to investigate the temporal evolution of cluster scale height, using homogeneously selected sub-samples of Kharchenko et al. (MWSC), Dias et al. (DAML02), WEBDA, and Froebrich et al. (FSR). We identify a linear relationship between scale height and log(age/yr) of clusters, considerably different from field stars. The scale height increases from about 40pc at 1Myr to 75pc at 1Gyr, most likely due to internal evolution and external scattering events. After 1Gyr, there is a marked change of the behaviour, with the scale height linearly increasing with log(age/yr) to about 550pc at 3.5Gyr. The most likely interpretation is that the surviving clusters are only observable because they have been scattered away from the mid-plane in their past. A detailed understanding of this observational evidence can only be achieved with numerical simulations of the evolution of cluster samples in the Galactic Disk. Furthermore, we find a weak trend of an age-independent increase in scale height with galactocentric distance. There are no significant temporal or spatial variations of the cluster distribution zero point. We determine the Suns vertical displacement from the Galactic Plane as $Z_odot=18.5pm1.2$pc.
Clusters of galaxies, filled with hot magnetized plasma, are the largest bound objects in existence and an important touchstone in understanding the formation of structures in our Universe. In such clusters, thermal conduction follows field lines, so magnetic fields strongly shape the clusters thermal history; that some have not since cooled and collapsed is a mystery. In a seemingly unrelated puzzle, recent observations of Virgo cluster spiral galaxies imply ridges of strong, coherent magnetic fields offset from their centre. Here we demonstrate, using three-dimensional magnetohydrodynamical simulations, that such ridges are easily explained by galaxies sweeping up field lines as they orbit inside the cluster. This magnetic drape is then lit up with cosmic rays from the galaxies stars, generating coherent polarized emission at the galaxies leading edges. This immediately presents a technique for probing local orientations and characteristic length scales of cluster magnetic fields. The first application of this technique, mapping the field of the Virgo cluster, gives a startling result: outside a central region, the magnetic field is preferentially oriented radially as predicted by the magnetothermal instability. Our results strongly suggest a mechanism for maintaining some clusters in a non-cooling-core state.