Polar ring galaxies are ideal objects with which to study the three-dimensional shapes of galactic gravitational potentials since two rotation curves can be measured in two perpendicular planes. Observational studies have uncovered systematically larger rotation velocities in the extended polar rings than in the associated host galaxies. In the dark matter context, this can only be explained through dark halos that are systematically flattened along the polar rings. Here, we point out that these objects can also be used as very effective tests of gravity theories, such as those based on Milgromian dynamics (MOND). We run a set of polar ring models using both Milgromian and Newtonian dynamics to predict the expected shapes of the rotation curves in both planes, varying the total mass of the system, the mass of the ring with respect to the host, as well as the size of the hole at the center of the ring. We find that Milgromian dynamics not only naturally leads to rotation velocities being typically higher in the extended polar rings than in the hosts, as would be the case in Newtonian dynamics without dark matter, but that it also gets the shape and amplitude of velocities correct. Milgromian dynamics thus adequately explains this particular property of polar ring galaxies.

We analyze a sample of 23 supermassive elliptical galaxies (central velocity dispersion larger than 330 km s-1), drawn from the SDSS. For each object, we estimate the dynamical mass from the light profile and central velocity dispersion, and compare it with the stellar mass derived from stellar population models. We show that these galaxies are dominated by luminous matter within the radius for which the velocity dispersion is measured. We find that the sizes and stellar masses are tightly correlated, with Re ~ M*^{1.1}$, making the mean density within the de Vaucouleurs radius a steeply declining function of M*: rho_e ~ M*^{-2.2}. These scalings are easily derived from the virial theorem if one recalls that this sample has essentially fixed (but large) sigma_0. In contrast, the mean density within 1 kpc is almost independent of M*, at a value that is in good agreement with recent studies of z ~ 2 galaxies. The fact that the mass within 1 kpc has remained approximately unchanged suggests assembly histories that were dominated by minor mergers -- but we discuss why this is not the unique way to achieve this. Moreover, the total stellar mass of the objects in our sample is typically a factor of ~ 5 larger than that in the high redshift (z ~ 2) sample, an amount which seems difficult to achieve. If our galaxies are the evolved objects of the recent high redshift studies, then we suggest that major mergers were required at z > 1.5, and that minor mergers become the dominant growth mechanism for massive galaxies at z < 1.5.

Recently, the distribution of velocity dispersion as far as 400kpc around red isolated galaxies was derived from statistical studies of satellites in the SDSS (Klypin & Prada 2009). This could help to constrain dark matter models at intermediate scales. We compare the predictions of different DM distributions, LCDM with NFW or cored profiles, and also modified gravity models, with observations. It is shown how the freedom in the various parameters (radial distribution of satellites, velocity anisotropy, external field effect), prevents to disentangle the models, which all can give pretty good fits to the data. In all cases, realistic radial variations of velocity anisotropy are used for the satellites, and a constant stellar-mass to light ratio for the host galaxies.

We review galaxy formation and dynamics under the MOND hypothesis of modified gravity, and compare to similar galaxies in Newtonian dynamics with dark matter. The aim is to find peculiar predictions both to discriminate between various hypotheses, and to make the theory progress through different constraints, touching the interpolation function, or the fundamental acceleration scale. Galaxy instabilities, forming bars and bulges at longer term, evolve differently in the various theories, and help to bring constraints, together with the observations of bar frequency. Dynamical friction and the predicted merger rate could be a sensitive test of theories. The different scenarios of galaxy formation are compared within the various theories and observations.

We consider for the first time the implications on the modified gravity MOND model of galaxies, of the presence of dark baryons, under the form of cold molecular gas in galaxy discs. We show that MOND models of rotation curves are still valid and universal, but the critical acceleration a0 separating the Newtonian and MONDian regimes has a lower value. We quantify this modification, as a function of the scale factor c between the total gas of the galaxy and the measured atomic gas. The main analysis concerns 43 resolved rotation curves and allows us to find the best pair (a0 = 0.96 10e-10 m.s-2, c = 3), which is also compatible to the one obtained from a second method by minimizing the scatter in the baryonic Tully-Fisher relation.

The stability of spiral galaxies is compared in modified Newtonian Dynamics (MOND) and Newtonian dynamics with dark matter (DM). We extend our previous simulations that involved pure stellar discs without gas, to deal with the effects of gas dissipation and star formation. We also vary the interpolating function between the MOND and Newtonian regime. Bar formation is compared in both dynamics, from initial conditions identical in visible component. One first result is that the MOND galaxy evolution is not affected by the choice of the mu-function, it develops bars with the same frequency and strength. The choice of the mu-function significantly changes the equivalent DM models, in changing the dark matter to visible mass ratio and, therefore, changing the stability. The introduction of gas shortens the timescale for bar formation in the DM model, but is not significantly shortened in the MOND model. As a consequence, with gas, the MOND and DM bar frequency histograms are now more similar than without gas. The thickening of the plane occurs through vertical resonance with the bar and peanut formation, and even more quickly with gas. Since the mass gets more concentrated with gas, the radius of the peanut is smaller, and the appearance of the pseudo-bulge is more boxy. The bar strength difference is moderated by saturation, and feedback effects, like the bar weakening or destruction by gas inflow due to gravity torques. Averaged over a series of models representing the Hubble sequence, the MOND models have still more bars, and stronger bars, than the equivalent DM models, better fitting the observations. Gas inflows driven by bars produce accumulations at Lindblad resonances, and MOND models can reproduce observed morphologies quite well, as was found before in the Newtonian dynamics.

We compare N-body simulations performed in MOND with analogs in Newtonian gravity with dark matter (DM). We have developed a code which solves the Poisson equation in both gravity models. It is a grid solver using adaptive mesh refinement techniques, allowing us to study isolated galaxies as well as interacting galaxies. Galaxies in MOND are found to form bars faster and stronger than in the DM model. In Newton dynamics, it is difficult to reproduce the observed high frequency of strong bars, while MOND appears to fit better the observations. Galaxy interactions and mergers, such as the Antennae, are also simulated with Newton and MOND dynamics. In the latter, dynamical friction is much weaker, and merging time-scales are longer. The formation of tidal dwarf galaxies in tidal tails are also compared in MOND and Newton+DM models.

We investigate how different models that have been proposed for solving the dark matter problem can fit the velocity dispersion observed around elliptical galaxies, on either a small scale (~ 20kpc) with stellar tracers, such as planetary nebulae, or large scale (~ 200kpc) with satellite galaxies as tracers. Predictions of Newtonian gravity, either containing pure baryonic matter, or embedded in massive cold dark matter (CDM) haloes, are compared with predictions of the modified gravity of MOND. The standard CDM model has problems on a small scale, and the Newtonian pure baryonic model has difficulties on a large scale, while a fit with MOND is possible on both scales.

The LCDM model is the most commonly admitted to describe our Universe. In spite of a great success with regard to the large scale structure formation, some problems are still unresolved at galactic scales. Alternative scenarios have to be explored such as modified gravity. We have developed an N-body code able to solve in a self consistent way the galactic dynamics in MOND. The first version of the code consists in solving the modified Poisson equation on a uniform Cartesian grid to derive the gravitational force on each particle. With it, we study the evolution of isolated galaxies, like the bar instability, the angular momentum transfer, etc. Galaxies in MOND are found to form stronger bars, faster than in Newtonian dynamics with dark matter. In a second step, we implement an adaptive mesh refinement technique in the code, allowing to run more contrasted simulations on larger scales, like interacting galaxies. During an interaction, the dynamical friction forces are less important in MOND, and merging times are longer than in DM models. The different morphologies of interacting galaxies in the two models are discussed. All simulations are performed in both frameworks of modified gravity and Newtonian gravity with dark matter with equivalent initial conditions.

The recent SDSS measured velocity distribution of satellite galaxies has been modelled in the context of MOND. We show that even when the extra constraint of adhering to the projected satellite number density profile is added, the two line of sight (los) velocity dispersion profiles presented in Klypin & Prada (2007) can be matched simply with a radially varying anisotropy. Interestingly, the anisotropies required to fit the los velocity dispersions are remarkably similar to the anisotropies generated by dissipationless collapse simulations in MOND. The mass-to-light ratios of the two host galaxies used are sensible and positivity of the distribution function is satisfied.

Within the cold dark matter (CDM) framework tidal dwarf galaxies (TDGs) cannot contain dark matter, so the recent results by Bournaud et al. (2007) that 3 rotating TDGs do show significant evidence for being dark matter dominated is inconsistent with the current concordance cosmological theory unless yet another dark matter component is postulated. We confirm that the TDG rotation curves are consistent with Newtonian dynamics only if either an additional dark matter component is postulated, or if all 3 TDGs happen to be viewed nearly edge-on, which is unlikely given the geometry of the tidal debris. We also find that the observed rotation curves are very naturally explained without any free parameters within the modified Newtonian dynamics (MOND) framework if inclinations are adopted as derived by Bournaud et al. We explore different inclination angles and two different assumptions about the external field effect. The results do not change significantly, and we conclude therefore that Newtonian dynamics has severe problems while MOND does exceedingly well in explaining the observed rotation curves of the 3 TDGs studied by Bournaud et al.

The self-binding energy and stability of a galaxy in MOND-based gravity are curiously decreasing functions of its center of mass acceleration towards neighbouring mass concentrations. A tentative indication of this breaking of the Strong Equivalence Principle in field galaxies is the RAVE-observed escape speed in the Milky Way. Another consequence is that satellites of field galaxies will move on nearly Keplerian orbits at large radii (100 - 500 kpc), with a declining speed below the asymptotically constant naive MOND prediction. But consequences of an environment-sensitive gravity are even more severe in clusters, where member galaxies accelerate fast: no more Dark-Halo-like potential is present to support galaxies, meaning that extended axisymmetric disks of gas and stars are likely unstable. These predicted reappearance of asymptotic Keplerian velocity curves and disappearance of "stereotypic galaxies" in clusters are falsifiable with targeted surveys.

We compare N-body simulations of isolated galaxies performed in both frameworks of modified Newtonian dynamics (MOND) and Newtonian gravity with dark matter (DM). We have developed a multigrid code able to efficiently solve the modified Poisson equation derived from the Lagrangian formalism AQUAL. We take particular care of the boundary conditions that are a crucial point in MOND. The 3-dimensional dynamics of initially identical stellar discs is studied in both models. In Newtonian gravity the live DM halo is chosen to fit the rotation curve of the MOND galaxy. For the same value of the Toomre parameter (Q_T), galactic discs in MOND develop a bar instability sooner than in the DM model. In a second phase the MOND bars weaken while the DM bars continue to grow by exchanging angular momentum with the halo. The bar pattern speed evolves quite differently in the two models: there is no dynamical friction on the MOND bars so they keep a constant pattern speed while the DM bars slow down significantly. This affects the position of resonance like the corotation and the peanut. The peanut lobes in the DM model move radially outward while they keep the same position in MOND. Simulations of (only stellar) galaxies of different types on the Hubble sequence lead to a statistical bar frequency that is closer to observations for the MOND than the DM model.