Christian M. Boily

The great disk of Milky-Way satellites and cosmological sub-structures

We show that the shape of the observed distribution of Milky Way (MW) satellites is inconsistent with their being drawn from a cosmological sub-structure population with a confidence of 99.5 per cent. Most of the MW satellites therefore cannot be related to dark-matter dominated satellites.

Local-Group tests of dark-matter concordance cosmology - Towards a new paradigm for structure formation

Predictions of the concordance cosmological model (CCM) of the structures in the environment of large spiral galaxies are compared with observed properties of Local Group galaxies. Five new, most probably irreconcilable problems are uncovered: 1) A wide variety of published CCM models consistently predict some form of relation between dark-matter-mass and luminosity for the Milky Way (MW) satellite galaxies, but none is observed. 2) The mass function of luminous sub-haloes predicted by the CCM contains too few satellites with dark matter (DM) mass ≈10 7 Mwithin their innermost 300 pc than in the case of the MW satellites. 3) The Local Group galaxies and data from extragalactic surveys indicate there is a correlation between bulge-mass and the number of luminous satellites that is not predicted by the CCM. 4) The 13 new ultra-faint MW satellites define a disc-of-satellites (DoS) that is virtually identical to the DoS previously found for the 11 classical MW satellites, implying that most of the 24 MW satellites are correlated in phase-space. 5) The occurrence of two MW-type DM halo masses hosting MW-like galaxies is unlikely in the CCM. However, the properties of the Local Group galaxies provide information leading to a solution of the above problems. The DoS and bulge-satellite correlation suggest that dissipational events forming bulges are related to the processes forming phase-space correlated satellite populations. These events are well known to occur since in galaxy encounters energy and angular momentum are expelled in the form of tidal tails, which can fragment to form populations of tidal-dwarf galaxies (TDGs) and associated star clusters. If Local Group satellite galaxies are to be interpreted as TDGs then the substructure predictions of the CCM are internally in conflict. All findings thus suggest that the CCM does not account for the Local Group observations and that therefore existing as well as new viable alternatives have to be further explored. These are discussed and natural solutions for the above problems emerge.

On the mass function of star clusters

Clusters that form in total 10 3 <N <10 5 stars (type II clusters) lose their gas within a dynamical time as a result of the photo-ionising flux from O stars. Sparser (type I) clusters get rid of their residual gas on a timescale longer or comparable to the nominal crossing time and thus evolve approximately adiabatically. This is also true for massive embedded clusters (type III) for which the velocity dispersion is larger than the sound speed of the ionised gas. On expelling their residual gas, type I and III clusters are therefore expected to lose a smaller fraction of their stellar component than type II clusters. We outline the effect this has on the transformation of the mass function of embedded clusters (ECMF), which is directly related to the mass function of star-cluster-forming molecular cloud cores, to the “initial” MF of bound gas-free star clusters (ICMF). The resulting ICMF has, for a featureless power-law ECMF, a turnover near 10 4.5 M⊙ and a peak near 10 3 M⊙. The peak lies around the initial masses of the Hyades, Praesepe and Pleiades clusters. We also find that the entire Galactic population II stellar spheroid can be generated if star formation proceeded via embedded clusters distributed like a power-law MF with exponent 0.9< � < 2.6.

SUPERBOX – an efficient code for collisionless galactic dynamics

We present SUPERBOX, a particle-mesh code with high resolution sub-grids and an NGP (nearest grid point) force-calculation scheme based on the second derivatives of the potential. SUPERBOX implements a fast low-storage FFT-algorithm, giving the possibility to work with millions of particles on desk-top computers. Test calculations show energy and angular momentum conservation to one part in 10(5) per crossing-time. The effects of grid and numerical relaxation remain negligible, even when these calculations cover a Hubble-time of evolution. As the sub-grids follow the trajectories of individual galaxies, the code allows a highly resolved treatment of interactions in clusters of galaxies, such as high-velocity encounters between elliptical galaxies and the tidal disruption of dwarf galaxies. Excellent agreement is obtained in a comparison with a direct-summation N-body code running on special-purpose GRAPE 3 hardware. The orbital decay of satellite galaxies due to dynamical friction obtained with SUPERBOX agrees with Chandrasekhars treatment when the Coulomb logarithm ln Lambda approximate to 1.5

Low-velocity streams in the solar neighbourhood caused by the Galactic bar

We find that a steady state bar induces transient features at low velocities in the solar neighborhood velocity distribution due to the initial response of the disc, following the formation of the bar. We associate these velocity streams with two quasi-periodic orbital families, librating around the stable x1(1) and x1(2) orbits near the bar’s outer Lindblad resonance (OLR). In a reference frame moving with the bar, these otherwise stationary orbits precess on a timescale dependent on the strength of the bar, consistent with predictions from a simple Hamiltonian model for the resonance. This behavior allows the two orbital families to reach the solar neighborhood and manifest themselves as clumps in the u-v plane moving away from (x1(2)), and toward (x1(1)) the Galactic center. Depending on the bar parameters and time since its formation, this model is consistent with the Pleiades and Coma Berenices, or Pleiades and Sirius moving groups seen in the Hipparcos stellar velocity distribution, if the Milky Way bar angle is 30 � . �0 . 45 � and its pattern speed is b=0 = 1:82±0:07, where 0 is the angular velocity of the local standard of rest (LSR). Since the process is recurrent, we can achieve a good match about every six LSR rotations. However, to be consistent with the fraction of stars in the Pleiades, we estimate that the Milky Way bar formed � 2 Gyr ago. This model argues against a common dynamical origin for the Hyades and Pleiades moving groups. Subject headings: stellar dynamics, Galactic bar, solar neighborhood

Evolution of star clusters in arbitrary tidal fields

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 starburst 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.

On the mass of dense star clusters in starburst galaxies from spectrophotometry

The mass of unresolved young star clusters derived from spectrophotometric data may well be off by a factor of 2 or more once the migration of massive stars driven by mass segregation is accounted for. We quantify this effect for a large set of cluster parameters, including variations in the stellar initial mass function (IMF), the intrinsic cluster mass, and mean mass density. Gas-dynamical models coupled with the Cambridge stellar evolution tracks allow us to derive a scheme to recover the real cluster mass given measured half-light radius, one-dimensional velocity dispersion and age. We monitor the evolution with time of the ratio of real to apparent mass through the parameter η. When we compute η for rich star clusters, we find non-monotonic evolution in time when the IMF stretches beyond a critical cut-off mass of 25.5 M� . We also monitor the rise of colour gradients between the inner and outer volume of clusters: we find trends in time of the stellar IMF power indices overlapping well with those derived for the Large Magellanic Cloud cluster NGC 1818 at an age of 30 Myr. We argue that the core region of massive Antennae clusters should have suffered from much segregation despite their low ages. We apply these results to a cluster mass function, and find that the peak of the mass distribution would appear to observers shifted to lower masses by as much as 0.2 dex. The star formation rate derived for the cluster population is then underestimated by from 20 to 50 per cent.

One moment in time - modeling star formation in the antennae

We present a new high-resolution N-body/smoothed particle hydrodynamics simulation of an encounter of two gas-rich disk galaxies that closely matches the morphology and kinematics of the interacting Antennae galaxies (NGC 4038/39). The simulation includes radiative cooling, star formation, and feedback from Type II supernovae. The large-scale morphology and kinematics are determined by the internal structure and the orbit of the progenitor disks. The properties of the central region, in particular the starburst in the overlap region, only match the observations for a very short time interval (Δt ≈ 20 Myr) after the second encounter. This indicates that the Antennae galaxies are in a special phase only about 40 Myr after the second encounter and 50 Myr before their final collision. This is the only phase in the simulation when a gas-rich overlap region between the nuclei is forming accompanied by enhanced star formation. The star formation rate as well as the recent star formation history in the central region agree well with observational estimates. For the first time, this new model explains the distributed extra-nuclear star formation in the Antennae galaxies as a consequence of the recent second encounter. The proposed model predicts that the Antennae are in a later merger stage than the Mice (NGC 4676) and would therefore lose their first place in the classical Toomre sequence. (Less)

Satellite decay in flattened dark matter haloes

We carry out a set of self-consistent N-body calculations to compare the decay rates of satellite dwarf galaxies orbiting a disc galaxy embedded in a dark matter halo (DMH). We consider both spherical and oblate axisymmetric DMHs of aspect ratio qh = 0.6. The satellites are given different initial orbital inclinations, orbital periods and mass. The live flattened DMHs with embedded discs and bulges are set-up using a new fast algorithm, MaGalie (Boily, Kroupa & Penarrubia 2001). We find that the range of survival times of satellites within a flattened DMH becomes � 100% larger than the same satellites within a spherical DMH. In the oblate DMH, satellites on polar orbits have the longest survival time, whereas satellites on coplanar prograde orbits are destroyed most rapidly. The orbital plane of a satellite tilts as a result of anisotropic dynamical friction, causing the satellites orbit to align with the plane of symmetry of the DMH. Polar orbits are not subjected to alignment. Therefore the decay of a satellites in an axisymmetric DMH may provide a natural explanation for the observed lack of satellites within 0 30 ◦ of their host galaxys disc (Holmberg 1969; Zaritsky & Gonzalez 1999). The computations furthermore indicate that the evolution of the orbital eccen- tricity e is highly dependent of its initial value e(t = 0) and the DMHs shape. We also discuss some implications of flattened DMHs for satellite debris streams.

Efficient N-body realisations of axisymmetric galaxies and halos

We present an efficient method for building equilibrium multi-component galaxies with non-spherical haloes and bulges. The gist of our approach is to tailor the velocity ellipsoid directly to the geometry of the mass distribution. Thus we avoid computing the anisotropic velocity dispersions which leads to large savings in the cpu budget. The computational time of the algorithm for triaxial equilibria scales linearly with the number of particles, N. The approximate solution to the velocity field causes structural relaxation: tests with N = 50,000 show that fluctuations of the inertia tensor (not exceeding the 10 per cent level) disappear within one half of a revolution at the half-mass radius. At later times equilibrium properties settle to values close to those sought from the initial conditions. A disc component is then added as described by Hernquist (1993). Incorporating the above algorithm to his code BuildGal, test runs show that the stability of the disc against vertical heating is not substantially modified by using our method. The code, MaGali, is made generally available.Abstract We present an efficient method for building equilibrium multi-component galaxies with non-spherical haloes and bulges. The gist of our approach is to tailor the velocity ellipsoid directly to the geometry of the mass distribution. Thus we avoid computing the anisotropic velocity dispersions which leads to large savings in the cpu budget. The computational time of the algorithm for triaxial equilibria scales linearly with the number of particles, N. The approximate solution to the velocity field causes structural relaxation: tests with N=50,000 show that fluctuations of the inertia tensor (not exceeding the 10% level) disappear within one half of a revolution at the half-mass radius. At later times equilibrium properties settle to values close to those sought from the initial conditions. A disc component is then added as described by Hernquist [ApJS, 86 (1993) 389]. Incorporating the above algorithm to his code buildgal , test runs show that the stability of the disc against vertical heating is not substantially modified using our method. The code, M a G a L ie , is made generally available (contact the corresponding author).