We investigate models of the Milky Way disc taking into account simultaneously the bar and a two-armed quasi-static spiral pattern. Away from major resonance overlaps, the mean stellar radial motions in the plane are essentially a linear superposition of the isolated effects of the bar and spirals. Thus, provided the bar is strong enough, even in the presence of spiral arms, these mean radial motions are predominantly affected by the Galactic bar for large scale velocity fluctuations. This is evident when comparing the peculiar line-of-sight velocity power spectrum of our coupled models with bar-only models. However, we show how forthcoming spectroscopic surveys could disentangle bar-only non-axisymmetric models of the Galaxy from models in which spiral arms have a significant amplitude. We also point out that overlaps of low-order resonances are sufficient to enhance stellar churning within the disc, even when the spirals amplitude is kept constant. Nevertheless, for churning to be truly non-local, stronger or (more likely) transient amplitudes would be needed: otherwise the disc is actually mostly unaffected by churning in the present models. Finally, regarding vertical breathing modes, the combined effect of the bar and spirals on vertical motions is a clear non-linear superposition of the isolated effects of both components, significantly superseding the linear superposition of modes produced by each perturber separately, thereby providing an additional effect to consider when analysing the observed breathing mode of the Galactic disc in the extended Solar neighbourhood.

This paper discusses a new method to generate self-coherent initial conditions for young substructured stellar cluster. The expansion of a uniform system allows stellar sub-structures (clumps) to grow from fragmentation modes by adiabatic cooling. We treat the system mass elements as stars, chosen according to a Salpeter mass function, and the time-evolution is performed with a collisional N-body integrator. This procedure allows to create a fully-coherent relation between the clumps' spatial distribution and the underlying velocity field. The cooling is driven by the gravitational field, as in a cosmological Hubble-Lema\^itre flow. The fragmented configuration has a `fractal'-like geometry but with a self-grown velocity field and mass profile. We compare the characteristics of the stellar population in clumps with that obtained from hydrodynamical simulations and find a remarkable correspondence between the two in terms of the stellar content and the degree of spatial mass-segregation. In the fragmented configuration, the IMF power index is ~0.3 lower in clumps in comparison to the field stellar population, in agreement with observations in the Milky Way. We follow in time the dynamical evolution of fully fragmented and sub-virial configurations, and find a soft collapse, leading rapidly to equilibrium (timescale of 1 Myr for a ~ 10^4 Msun system). The low-concentration equilibrium implies that the dynamical evolution including massive stars is less likely to induce direct collisions and the formation of exotic objects. Low-mass stars already ejected from merging clumps are depleted in the end-result stellar clusters, which harbour a top-heavy stellar mass function.

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.

(abridged) 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. However, the Local Group properties provide hints that may lead 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. Such 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 sub-structure predictions of 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.

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 x_1(1) and x_1(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 (x_1(2)), and toward (x_1(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<phi_0<45 [deg] and its pattern speed is Omega_b/Omega_0=1.82\pm 0.07, where Omega_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.

We present a new high-resolution N-body/SPH simulation of an encounter of two gas-rich disk galaxies which closely matches the morphology and kinematics of the interacting Antennae galaxies (NGC 4038/39). The simulation includes radiative cooling, star formation and feedback from SNII. 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 of ~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.

The disruptive effect of galactic tides is a textbook example of gravitational dynamics. However, depending on the shape of the potential, tides can also become fully compressive. When that is the case, they might trigger or strengthen the formation of galactic substructures (star clusters, tidal dwarf galaxies), instead of destroying them. We perform N-body simulations of interacting galaxies to quantify this effect. We demonstrate that tidal compression occurs repeatedly during a galaxy merger, independently of the specific choice of parameterization. With a model tailored to the Antennae galaxies, we show that the distribution of compressive tides matches the locations and timescales of observed substructures. After extending our study to a broad range of parameters, we conclude that neither the importance of the compressive tides (~15% of the stellar mass) nor their duration (~ 10 Myr) are strongly affected by changes in the progenitors' configurations and orbits. Moreover, we show that individual clumps of matter can enter compressive regions several times in the course of a simulation. We speculate that this may spawn multiple star formation episodes in some star clusters, through e.g., enhanced gas retention.

We point out a strong time-evolution of the mass-to-light conversion factor eta commonly used to estimate masses of unresolved star clusters from observed cluster spectro-photometric measures. We present a series of gas-dynamical models coupled with the Cambridge stellar evolution tracks to compute line-of-sight velocity dispersions and half-light radii weighted by the luminosity. We explore a range of initial conditions, varying in turn the cluster mass and/or density, and the stellar population's IMF. We find that eta, and hence the estimated cluster mass, may increase by factors as large as 3 over time-scales of 50 million years. We apply these results to an hypothetic cluster mass distribution function (d.f.) and show that the d.f. shape may be strongly affected at the low-mass end by this effect. Fitting truncated isothermal (Michie-King) models to the projected light profile leads to over-estimates of the concentration parameter c of delta c ~ 0.3 compared to the same functional fit applied to the projected mass density.

Physics lectures always refer to the tides as a disruptive effect. However, tides can also be compressive. When the potential of two galaxies overlap, as happens during a merger, fully compressive tides can develop and have a strong impact on the dynamics of substructures such as star clusters or tidal dwarf galaxies. Using N-body simulations of a large set of mergers, we noticed the importance of these tidal modes at cluster scale. With a model of the Antennae galaxies, we conclude that the positions and timescales of these tidal modes match the actual distribution of young clusters. A detailed study of the statistics of the compressive tides shows a stunning correlation between this purely gravitational effect and the observed properties of the star clusters. In this contribution, we introduce the concept of compressive tide and show its relevance in the special case of the Antennae galaxies. We extend our conclusions to a broad range of parameters and discuss their implications on several critical points such as the infant mortality, multiple star formation epochs in clusters or the age distribution.

The tidal field of galaxies is known generally to be disruptive. However, in the case of galaxy mergers, a compressive mode of tidal wave may develop and last long enough to cocoon the formation of star clusters. Using an N-body simulation of the Antennae galaxies, we derive the positions of these compressive regions and the statistics of their duration. Excellent agreement between the spatial distribution of tides and observed young clusters is found, while the characteristic e-folding times of 10 to 30 Myrs derived for the tidal field compare well with cluster formation time-scales.

In the framework of hierarchical structure formation ellipticals can form from merging of smaller disk galaxies. The nearby interacting 'Antennae' galaxy pair (NGC 4038/39) is one of the best-studied local systems of merging spirals, thus presenting us with an ideal laboratory for the study of galaxy evolution models. The Antennae are believed to be in a state prior to their final encounter with rapid subsequent merging, which puts them in the first position of the Toomre (1977) merger sequence. Here we present first numerical high-resolution, self-consistent, smoothed particle hydrodynamics (SPH) simulations of the Antennae system, including star formation and stellar feedback, and compare our results to VLA HI observations by Hibbard et al. (2001). We are able to obtain a close, but not yet perfect match to the observed morphology and kinematics of the system.

Gravitational tides are widely understood to strip and destroy galactic substructures. In the course of a galaxy merger, however, transient totally compressive tides may develop and prevent star forming regions from dissolving, after they condensed to form clusters of stars. We study the statistics of such compressive modes in an N-body model of the galaxy merger NGC 4038/39 (the Antennae) and show that ~15% of the disc material undergoes compressive tides at pericentre. The spatial distribution of observed young clusters in the overlap and nuclear regions of the Antennae matches surprisingly well the location of compressive tides obtained from simulation data. Furthermore, the statistics of time intervals spent by individual particles embedded in a compressive tide yields a log-normal distribution of characteristic time ~10 Myr, comparable to star cluster formation timescales. We argue that this generic process is operative in galaxy mergers at all redshifts and possibly enhances the formation of star clusters. We show with a model calculation that this process will prevent the dissolution of a star cluster during the formation phase, even for a star formation efficiency as low as ~10%. The transient nature of compressive tides implies that clusters may dissolve rapidly once the tidal field switches to the usual disruptive mode.

We implement an efficient method to quantify time-dependent orbital complexity in gravitational $N$-body simulations. The technique, which we name DWaTIM, is based on a discrete wavelet transform of velocity orbital time series. The wavelet power-spectrum is used to measure trends in complexity continuously in time. We apply the method to the test cases N=3 Pythagorean- and a perturbed N=5 Caledonian configurations. The method recovers the well-known time-dependent complexity of the dynamics in these small-$N$ problems. We then apply the technique to an equal-mass collisional N=256 body simulation ran through core-collapse. We find that a majority of stars evolve on relatively complex orbits up to the time when the first hard binary forms, whereas after core-collapse, less complex orbits are found on the whole as a result of expanding mass shells.

The motion of a black hole about the centre of gravity of its host galaxy induces a strong response from the surrounding stellar population. We treat the case of a harmonic potential analytically and show that half of the stars on circular orbits in that potential shift to an orbit of lower energy, while the other half receive a positive boost and recede to a larger radius. The black hole itself remains on an orbit of fixed amplitude and merely acts as a catalyst for the evolution of the stellar energy distribution function f(E). We show that this effect is operative out to a radius of approx 3 to 4 times the hole's influence radius, R_bh. We use numerical integration to explore more fully the response of a stellar distribution to black hole motion. We consider orbits in a logarithmic potential and compare the response of stars on circular orbits, to the situation of a `warm' and `hot' (isotropic) stellar velocity field. While features seen in density maps are now wiped out, the kinematic signature of black hole motion still imprints the stellar line-of-sight mean velocity to a magnitude ~18% the local root mean-square velocity dispersion sigma.

We study the survival of substructures (clumps) within larger self-gravitating dark matter halos. Building on scaling relations obtained from N-body calculations of violent relaxation, we argue that the tidal field of galaxies and halos can only destroy substructures if spherical symmetry is imposed at formation. We explore other mechanisms that may tailor the number of halo substructures during the course of virialization. Unless the larger halo is built up from a few large clumps, we find that clump-clump encounters are unlikely to homogenize the halo on a dynamical timescale. Phase mixing would proceed faster in the inner parts and allow for the secular evolution of a stellar disk.

We present a purely theoretical study of the morphological evolution of self-gravitating systems formed through the dissipationless collapse of N-point sources. We explore the effects of resolution in mass and length on the growth of triaxial structures formed by an instability triggered by an excess of radial orbits. We point out that as resolution increases, the equilibria shift, from mildly prolate, to oblate. A number of particles N ~= 100000 or larger is required for convergence of axial aspect ratios. An upper bound for the softening, e ~ 1/256, is also identified. We then study the properties of a set of equilibria formed from scale-free cold initial mass distributions, ro ~ r^-g with 0 <= g <= 2. Oblateness is enhanced for initially more peaked structures (larger values of g). We map the run of density in space and find no evidence for a power-law inner structure when g <= 3/2 down to a mass fraction <~0.1 per cent of the total. However, when 3/2 < g <= 2, the mass profile in equilibrium is well matched by a power law of index ~g out to a mass fraction ~ 10 per cent. We interpret this in terms of less-effective violent relaxation for more peaked profiles when more phase mixing takes place at the centre. We map out the velocity field of the equilibria and note that at small radii the velocity coarse-grained distribution function (DF) is Maxwellian to a very good approximation.

The mass of unresolved young star clusters derived from spectro-photometric 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 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 eta. When we compute eta for rich star clusters, we find non-monotonic evolution in time when the IMF stretches beyond a critical cutoff mass of 25.5 solar mass. We also monitor the rise of color 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 LMC 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 (SFR) derived for the cluster population is then underestimated by from 20 to 50 per cent.

We point out a strong time-evolution of the mass-to-light conversion factor \eta commonly used to estimate masses of dense star clusters from observed cluster radii and stellar velocity dispersions. We use a gas-dynamical model coupled with the Cambridge stellar evolution tracks to compute line-of-sight velocity dispersions and half-light radii weighted by the luminosity. Stars at birth are assumed to follow the Salpeter mass function in the range [0.15--17 M_\sun]. We find that $\eta$, and hence the estimated cluster mass, increases by factors as large as 3 over time-scales of 20 million years. Increasing the upper mass limit to $50 M_\sun$ leads to a sharp rise of similar amplitude but in as little as 10 million years. Fitting truncated isothermal (Michie-King) models to the projected light profile leads to over-estimates of the concentration par ameter c of $\delta c\approx 0.3$ compared to the same functional fit applied to the proj ected mass density.

We show that the shape of the observed distribution of Milky Way (MW) satellites is inconsistent with 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.

We subject to an N-body numerical investigation our analysis of Paper I on the survival of stellar clusters undergoing rapid mass loss. We compare analytical tracks of bound mass-fraction {\it vs} star formation efficiency $\epsilon$ to those obtained with N-body integration. We use these to argue that stellar clusters must develop massive cores of high-binding energy if they are to remain bound despite a star formation efficiency as low as 30% or lower suggested by observations. The average local virial ratio $<\sigma^2/|\phi|>$ is introduced to classify bound clusters as function of their critical $\epsilon$ for dissolution. Clusters dissolving at lower $\epsilon$ achieve the lowest ratio. We applied this classification parameter successfully to Michie-King and Hernquist-type distribution functions. The Plummer sphere is exceptional in that it defies this and other classification parameters we tried. The reasons for the discrepancy include less effective energy redistribution during the expansion phase for this case.

We study analytically the disruptive effect of instantaneous gas removal from a cluster containing O stars. We setup an iterative calculation based on the stellar velocity distribution function to compute the fraction of stars that remain bound once the cluster has ejected the gas and is out of equilibrium. We show that the stellar bound fraction is a function of the initial cluster distribution function as well as the star formation efficiency, $\epsilon$, taken constant throughout the cluster. The case of the Plummer sphere is dealt with in greater details. We find for this case that up to ~ 50% of the stars may remain bound when $\epsilon$ assumes values < 1/2, contrary to expectations derived from the virial theorem. The fraction of bound stars is expressed algebraically for polytropic distribution functions.

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} Msun and a peak near 10^3 Msun. 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 < beta < 2.6.

The purpose of this article is to show that when dynamically cold, dissipationless self-gravitating systems collapse, their evolution is a strong function of the symmetry in the initial distribution. We explore with a set of pressure-less homogeneous fluids the time-evolution of ellipsoidal distributions and map the depth of potential achieved during relaxation as function of initial ellipsoid axis ratios. We then perform a series of $N$-body numerical simulations and contrast their evolution with the fluid solutions. We verify an analytic relation between collapse factor ${\cal C}$ and particle number $N$ in spherical symmetry, such that ${\cal C} \propto N^{1/3}$. We sought a similar relation for axisymmetric configurations, and found an empirical scaling relation such that ${\cal C} \propto N^{1/6}$ in these cases. We then show that when mass distributions do not respect spherical- or axial-symmetry, the ensuing gravitational collapse deepens with increasing particle number $N$ but only slowly: 86% of triaxial configurations may collapse by a factor of no more than 40 as $N\to\infty$. For $N\approx 10^5$ and larger, violent relaxation develops fully under the Lin-Mestel-Shu instability such that numerical $N$-body solutions now resolve the different initial morphologies adequately.

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 q_h=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 and Pe\~{n}arrubia 2001). We find that the range of survival times of satellites within a flattened DMH becomes of the order of 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 satellite's 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) degrees of their host galaxy's disc (Holmberg 1969; Zaritsky and Gonz\'alez 1999). The computations furthermore indicate that the evolution of the orbital eccentricity $e$ is highly dependent of its initial value e(t=0) and the DMH's shape. We also discuss some implications of flattened DMHs for satellite debris streams.

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.

We study both analytically and numerically the disruptive effect of instantaneous gas removal from an embedded cluster. We setup a calculation based on the stellar velocity distribution function, to compute the fraction of stars that remain bound once the cluster has ejected the gas and is out of equilibrium. We find tracks of bound mass-fraction vs star formation efficiency similar to those obtained with N-body calculations. We use these to argue that embedded clusters must develop high-binding energy cores if they are to survive as bound clusters despite a star formation rate as low as 20 per cent or lower suggested by observations.

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 Grape3 hardware. The orbital decay of satellite galaxies due to dynamical friction obtained with Superbox agrees with Chandrasekhar's treatment when the Coulomb logarithm is approximately 1.5.

The evolution of rotating, isolated clusters of stars up to core-collapse is investigated with n-body numerical codes. The simulations start off from axisymmetric generalisations of King profiles, with added global angular momentum. In this contribution I report on results obtained for two sets of single-mass cluster simulations. These confirm the more rapid evolution of even mildly-rotating clusters. A model is presented with rotational energy comparable to omega-Centauri's; it reaches core-collapse in less than half the time required for non-rotating model clusters.

1- Simulations of star clusters and galactic nuclei: rotation and other challenges (R. Spurzem & C. Einsel) 2- Relaxation of elliptical galaxies (C. Boily) 3- Simulations with an Hybrid Nbody-Scf code (M. Hemsendorf) These are three two-page contributions to the Rutgers galaxy dynamics held in august, 1998. Each will be included in the proceedings to be published by the Astronomical Society of the Pacific.