We test if Refracted Gravity (RG) can describe the dynamics of disk galaxies without resorting to dark matter. RG is a classical theory of gravity where the standard Poisson equation is modified by the gravitational permittivity, $\epsilon$, a universal monotonic function of the local mass density. We use the rotation curves and the vertical velocity dispersions of 30 galaxies in the DiskMass Survey (DMS) to determine $\epsilon$. RG describes the kinematic profiles with mass-to-light ratios in agreement with SPS models, and disk thicknesses in agreement with observations, once observational biases are considered. Our results rely on setting the three free parameters of RG for each galaxy. However, we show that the differences of these parameters from galaxy to galaxy could be ascribed to statistical fluctuations. We adopt an approximate method to find a single set of parameters that might properly describe the kinematics of the entire sample and suggest the universality of $\epsilon$. We finally show that the RG models of the individual rotation curves can only partly describe the radial acceleration relation (RAR). Evidently, the RG models underestimate the observed accelerations by 0.1-0.3 dex at low Newtonian accelerations. Another problem is the strong correlations, at largely more than 5$\sigma$, between the residuals of the RAR models and three radially-dependent properties of galaxies, whereas the DMS data show considerably less significant correlations, at more than 4$\sigma$, for only two of them. These correlations might originate the non-null intrinsic scatter of the RG models, at odds with the observed intrinsic scatter of galaxy samples, different from DMS, which is consistent with 0. Further studies are required to assess if these discrepancies in the RAR originate from the DMS sample, which might not be ideal for deriving the RAR, or if they are genuine failures of RG.

The Andromeda galaxy is observed to have a system of two large dwarf ellipticals and ~13 smaller satellite galaxies that are currently co-rotating in a thin plane, in addition to 2 counter-rotating satellite galaxies. We explored the consistency of those observations with a scenario where the majority of the co-rotating satellite galaxies originated from a subhalo group, where NGC 205 was the host and the satellite galaxies occupied dark matter sub-subhalos. We ran N-body simulations of a close encounter between NGC 205 and M31. In the simulations, NGC 205 was surrounded by massless particles to statistically sample the distribution of the sub-subhalos expected in a subhalo that has a mass similar to NGC 205. We made Monte Carlo samplings and found that, using a set of reference parameters, the probability of producing a thinner distribution of sub-subhalos than the observed NGC 205 + 15 smaller satellites (thus including the 2 counter-rotators, but excluding M32) increased from <1e-8 for the initial distribution to ~0.01 at pericentre. The probability of the simulated sub-subhalos occupying the locations of the observed co-rotating satellites in the line of sight velocity versus projected on-sky distance plane is at most 0.002 for 11 out of 13 satellites. Increasing the mass of M31 and the extent of the initial distribution of sub-subhalos gives a maximum probability of 0.004 for all 13 co-rotating satellites, but the probability of producing the thinness would drop to ~ 0.001.

The DiskMass survey recently provided measurements of the vertical velocity dispersions of disk stars in a sample of nearly face-on galaxies. By setting the disk scale-heights to be equal to those of edge-on galaxies with similar scale-lengths, it was found that these disks must be sub-maximal, with surprisingly low K-band mass-to-light ratios of the order of $M_\star/L_K \simeq 0.3 M_\odot/L_\odot$. This study made use of a simple relation between the disk surface density and the measured velocity dispersion and scale height of the disk, neglecting the shape of the rotation curve and the dark matter contribution to the vertical force, which can be especially important in the case of sub-maximal disks. Here, we point out that these simplifying assumptions led to an overestimation of the stellar mass-to-light ratios. Relaxing these assumptions, we compute even lower values than previously reported for the mass-to-light ratios, with a median $M_\star/L_K \simeq 0.18 M_\odot/L_\odot$, where 14 galaxies have $M_\star/L_K < 0.11$. Invoking prolate dark matter halos made only a small difference to the derived $M_\star/L_K$, although extreme prolate halos ($q>1.5$ for the axis ratios of the potential) might help. The cross-terms in the Jeans equation are also generally negligible. These deduced K-band stellar mass-to-light ratios are even more difficult to reconcile with stellar population synthesis models than the previously reported ones.

The Modified Newtonian Dynamics (MOND) paradigm generically predicts that the external gravitational field in which a system is embedded can produce effects on its internal dynamics. In this communication, we first show that this External Field Effect can significantly improve some galactic rotation curves fits by decreasing the predicted velocities of the external part of the rotation curves. In modified gravity versions of MOND, this External Field Effect also appears in the Solar System and leads to a very good way to constrain the transition function of the theory. A combined analysis of the galactic rotation curves and Solar System constraints (provided by the Cassini spacecraft) rules out several classes of popular MOND transition functions, but leaves others viable. Moreover, we show that LISA Pathfinder will not be able to improve the current constraints on these still viable transition functions.

This article explores the agreement between the predictions of Modified Newtonian Dynamics (MOND) and the rotation curves and stellar velocity dispersion profiles measured by the DiskMass Survey. A bulge-disk decomposition was made for each of the thirty published galaxies, and a MOND Poisson solver was used to simultaneously compute, from the baryonic mass distributions, model rotation curves and vertical velocity dispersion profiles, which were compared to the measured values. The two main free parameters, the stellar disk's mass-to-light ratio ($M/L$) and its exponential scale-height ($h_z$), were estimated by Markov Chain Monte Carlo modelling. The average best-fit K-band stellar mass-to-light ratio was $M/L \simeq 0.55 \pm 0.15$. However, to match the DiskMass Survey data, the vertical scale-heights would have to be in the range $h_z=200$ to $400$ pc which is a factor of two lower than those derived from observations of edge-on galaxies with a similar scale-length. The reason is that modified gravity versions of MOND characteristically require a larger $M/L$ to fit the rotation curve in the absence of dark matter and therefore predict a stronger vertical gravitational field than Newtonian models. It was found that changing the MOND acceleration parameter, the shape of the velocity dispersion ellipsoid, the adopted vertical distribution of stars, as well as the galaxy inclination, within any realistic range, all had little impact on these results.

This is a brief rebuttal to arXiv:1502.03821, which claims to provide the first observational proof of dark matter interior to the solar circle. We point out that this result is not new, and can be traced back at least a quarter century.

We present the stellar and gaseous kinematics of an Sb galaxy, NGC 3223, with the aim of determining the vertical and radial stellar velocity dispersion as a function of radius, which can help to constrain disk heating theories. Together with the observed NIR photometry, the vertical velocity dispersion is also used to determine the stellar mass-to-light (M/L) ratio, typically one of the largest uncertainties when deriving the dark matter distribution from the observed rotation curve. We find a vertical-to-radial velocity dispersion ratio of sigma_z/sigma_R=1.21+-0.14, significantly higher than expectations from known correlations, and a weakly-constrained Ks-band stellar M/L ratio in the range 0.5-1.7, at the high end of (but consistent with) the predictions of stellar population synthesis models. Such a weak constraint on the stellar M/L ratio, however, does not allow us to securely determine the dark matter density distribution. To achieve this, either a statistical approach or additional data (e.g. integral-field unit) are needed.

In recent years multi-wavelength observations have shown the presence of substructures related to merging events in a high fraction of galaxy clusters. Clusters can be roughly grouped into two categories -- relaxed and non-relaxed -- and a proper characterisation of the dynamical state of these systems is of crucial importance both for astrophysical and cosmological studies. In this paper we investigate the use of a number of morphological parameters (Gini, $M_{20}$, Concentration, Asymmetry, Smoothness, Ellipticity and Gini of the second order moment, $G_{M}$) introduced to automatically classify clusters as relaxed or dynamically disturbed systems. We apply our method to a sample of clusters at different redshifts extracted from the {\it Chandra} archive and we investigate possible correlations between morphological parameters and other X-ray gas properties. We conclude that a combination of the adopted parameters is a very useful tool to properly characterise the X-ray cluster morphology. According to our results three parameters -- Gini, $M_{20}$ and Concentration -- are very promising for identifying cluster mergers. The Gini coefficient is a particularly powerful tool, especially at high redshift, being independent from the choice of the position of the cluster centre. We find that high Gini ($>$ 0.65), high Concentration ($>$ 1.55) and low $M_{20}$ ($<$ -2.0) values are associated with relaxed clusters, while low Gini ($<$ 0.4), low Concentration ($<$ 1.0) and high $M_{20}$ ($>$ -1.4) characterise dynamically perturbed systems. We also estimate the X-ray cluster morphological parameters in the case of {\it radio loud} clusters. In excellent agreement with previous analyses we confirm that diffuse intracluster radio sources are associated with major mergers.

In this article we continue to test cosmological models centred on Modified Newtonian Dynamics (MOND) with light sterile neutrinos, which could in principle be a way to solve the fine-tuning problems of the standard model on galaxy scales while preserving successful predictions on larger scales. Due to previous failures of the simple MOND cosmological model, here we test a speculative model where the modified gravitational field is produced only by the baryons and the sterile neutrinos produce a purely Newtonian field (hence Modified Baryonic Dynamics). We use two component cosmological simulations to separate the baryonic N-body particles from the sterile neutrino ones. The premise is to attenuate the over-production of massive galaxy cluster halos which were prevalent in the original MOND plus light sterile neutrinos scenario. Theoretical issues with such a formulation notwithstanding, the Modified Baryonic Dynamics model fails to produce the correct amplitude for the galaxy cluster mass function for any reasonable value of the primordial power spectrum normalisation.

The classical dwarf spheroidals (dSphs) provide a critical test for Modified Newtonian Dynamics (MOND) because they are observable satellite galactic systems with low internal accelerations and low, but periodically varying, external acceleration. This varying external gravitational field is not commonly found acting on systems with low internal acceleration. Using Jeans modelling, Carina in particular has been demonstrated to require a V-band mass-to-light ratio greater than 5, which is the nominal upper limit for an ancient stellar population. We run MOND N-body simulations of a Carina-like dSph orbiting the Milky Way to test if dSphs in MOND are stable to tidal forces over the Hubble time and if those same tidal forces artificially inflate their velocity dispersions and therefore their apparent mass-to-light ratio. We run many simulations with various initial total masses for Carina, and Galactocentric orbits (consistent with proper motions), and compare the simulation line of sight velocity dispersions (losVDs) with the observed losVDs of Walker et al. (2007). We find that the dSphs are stable, but that the tidal forces are not conducive to artificially inflating the losVDs. Furthermore, the range of mass-to-light ratios that best reproduces the observed line of sight velocity dispersions of Carina is 5.3 to 5.7 and circular orbits are preferred to plunging orbits. Therefore, some tension still exists between the required mass-to-light ratio for the Carina dSph in MOND and those expected from stellar population synthesis models. It remains to be seen whether a careful treatment of the binary population or triaxiality might reduce this tension.

We use our Modified Newtonian Dynamics (MOND) cosmological particle-mesh N-body code to investigate the feasibility of structure formation in a framework involving MOND and light sterile neutrinos in the mass range 11 - 300 eV: always assuming that \Omega_{\nu_s}=0.225 for H_o=72 \kms Mpc^{-1}. We run a suite of simulations with variants on the expansion history, cosmological variation of the MOND acceleration constant, different normalisations of the power spectrum of the initial perturbations and interpolating functions. Using various box sizes, but typically with ones of length 256 Mpc/h, we compare our simulated halo mass functions with observed cluster mass functions and show that (i) the sterile neutrino mass must be larger than 30 eV to account for the low mass (M_{200}<10^{14.6} solar masses) clusters of galaxies in MOND and (ii) regardless of sterile neutrino mass or any of the variations we mentioned above, it is not possible to form the correct number of high mass (M_{200}>10^{15.1} solar masses) clusters of galaxies: there is always a considerable over production. This means that the ansatz of considering the weak-field limit of MOND together with a component of light sterile neutrinos to form structure from z ~ 200 fails. If MOND is the correct description of weak-field gravitational dynamics, it could mean that subtle effects of the additional fields in covariant theories of MOND render the ansatz inaccurate, or that the gravity generated by light sterile neutrinos (or by similar hot dark matter particles) is different from that generated by the baryons.

We examine the formation of clusters of galaxies in numerical simulations of a QUMOND cosmogony with massive sterile neutrinos. Clusters formed in these exploratory simulations develop higher velocities than those found in {\Lambda}CDM simulations. The bulk motions of clusters attain about 1000 km/s by low redshift, comparable to observations whereas {\Lambda}CDM simulated clusters tend to fall short. Similarly, high pairwise velocities are common in cluster-cluster collisions like the Bullet cluster. There is also a propensity for the most massive clusters to be larger in QUMOND and to appear earlier than in {\Lambda}CDM, potentially providing an explanation for 'pink elephants' like El Gordo. However, it is not obvious that the cluster mass function can be recovered.

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.

The bulge dominated galaxy NGC 7814 provides one of the strongest dynamical tests possible for Modified Newtonian Dynamics (MOND). Spitzer 3.6 micron photometry fixes the bulge parameterisation and strongly constrains the properties of the sub-dominant stellar disk. Furthermore, the distance is known to better than 5 percent, virtually eliminating it as a free parameter. The rotation curve is easily measured, since the H I (and stellar) disks are edge on, and both the receding and approaching sides agree very well. In this paper we explore the agreement between the model and observed rotation curves in MOND given that the only two free parameters available are the mass-to-light ratios of the bulge and disk. We use a grid based MOND Poisson solver that accurately solves for the MOND gravity and produces our model rotation curves from a given mass distribution. The input to the Poisson solver is a 3D distribution of N particles which is generated from modelling the observed distribution of stars and gas in the galaxy. By ensuring a superior fit to the radial surface brightness profile than previous works, by virtue of a double Sersic fit to the bulge, we were able to produce excellent fits to the rotation curve with typical values for both mass-to-light ratios. We conclude that the model rotation curve of a mass distribution in MOND is extremely sensitive to the bulge-disk decomposition and even slight deviation from the observed mass distribution can produce large differences in the model rotation curve.

General Relativity is able to describe the dynamics of galaxies and larger cosmic structures only if most of the matter in the Universe is dark, namely it does not emit any electromagnetic radiation. Intriguingly, on the scale of galaxies, there is strong observational evidence that the presence of dark matter appears to be necessary only when the gravitational field inferred from the distribution of the luminous matter falls below an acceleration of the order of 10^(-10) m/s^2. In the standard model, which combines Newtonian gravity with dark matter, the origin of this acceleration scale is challenging and remains unsolved. On the contrary, the full set of observations can be neatly described, and were partly predicted, by a modification of Newtonian dynamics, dubbed MOND, that does not resort to the existence of dark matter. On the scale of galaxy clusters and beyond, however, MOND is not as successful as on the scale of galaxies, and the existence of some dark matter appears unavoidable. A model combining MOND with hot dark matter made of sterile neutrinos seems to be able to describe most of the astrophysical phenomenology, from the power spectrum of the cosmic microwave background anisotropies to the dynamics of dwarf galaxies. Whether there exists a yet unknown covariant theory that contains General Relativity and Newtonian gravity in the weak field limit, and MOND as the ultra-weak field limit is still an open question.

Within the framework of modified Newtonian dynamics (MOND) we investigate the kinematics of two dwarf spiral galaxies belonging to very different environments, namely KK 246 in the Local Void and Holmberg II in the M81 group. A mass model of the rotation curve of KK 246 is presented for the first time, and we show that its observed kinematics are consistent with MOND. We re-derive the outer rotation curve of Holmberg II, by modelling its HI data cube, and find that its inclination should be closer to face-on than previously derived. This implies that Holmberg II has a higher rotation velocity in its outer parts, which, although not very precisely constrained, is consistent with the MOND prediction.

Several scenarios have been suggested to explain the phase-space distribution of the Milky Way (MW) satellite galaxies in a disc of satellites (DoS). To quantitatively compare these different possibilities, a new method analysing angular momentum directions in modelled data is presented. It determines how likely it is to find sets of angular momenta as concentrated and as close to a polar orientation as is observed for the MW satellite orbital poles. The method can be easily applied to orbital pole data from different models. The observed distribution of satellite orbital poles is compared to published angular momentum directions of subhalos derived from six cosmological state-of-the-art simulations in the Aquarius project. This tests the possibility that filamentary accretion might be able to naturally explain the satellite orbits within the DoS. For the most likely alignment of main halo and MW disc spin, the probability to reproduce the MW satellite orbital pole properties turns out to be less than 0.5 per cent in Aquarius models. Even an isotropic distribution of angular momenta has a higher likelihood to produce the observed distribution. The two Via Lactea cosmological simulations give results similar to the Aquarius simulations. Comparing instead with numerical models of galaxy-interactions gives a probability of up to 90 per cent for some models to draw the observed distribution of orbital poles from the angular momenta of tidal debris. This indicates that the formation as tidal dwarf galaxies in a single encounter is a viable, if not the only, process to explain the phase-space distribution of the MW satellite galaxies.

Here we present a new particle-mesh galactic N-body code that uses the full multigrid algorithm for solving the modified Poisson equation of the Quasi Linear formulation of Modified Newtonian Dynamics (QUMOND). A novel approach for handling the boundary conditions using a refinement strategy is implemented and the accuracy of the code is compared with analytical solutions of Kuzmin disks. We then employ the code to compute the predicted rotation curves for a sample of five spiral galaxies from the THINGS sample. We generated static N-body realisations of the galaxies according to their stellar and gaseous surface densities and allowed their distances, mass-to-light ratios (M/L) and both the stellar and gas scale-heights to vary in order to estimate the best fit parameters. We found that NGC 3621, NGC 3521 and DDO 154 are well fit by MOND using expected values of the distance and M/L. NGC 2403 required a moderately larger $M/L$ than expected and NGC 2903 required a substantially larger value. The surprising result was that the scale-height of the dominant baryonic component was well constrained by the rotation curves: the gas scale-height for DDO 154 and the stellar scale-height for the others. In fact, if the suggested stellar scale-height (one-fifth the stellar scale-length) was used in the case of NGC 3621 and NGC 3521 it would not be possible to produce a good fit to the inner rotation curve. For each of the four stellar dominated galaxies, we calculated the vertical velocity dispersions which we found to be, on the whole, quite typical compared with observed stellar vertical velocity dispersions of face on spirals. We conclude that modelling the gas scale-heights of the gas rich dwarf spiral galaxies will be vital in order to make precise conclusions about MOND.

The highly organised distribution of satellite galaxies surrounding the Milky Way is a serious challenge to the concordance cosmological model. Perhaps the only remaining solution, in this framework, is that the dwarf satellite galaxies fall into the Milky Way's potential along one or two filaments, which may or may not plausibly reproduce the observed distribution. Here we test this scenario by making use of the proper motions of the Fornax, Sculptor, Ursa Minor and Carina dwarf spheroidals, and trace their orbits back through several variations of the Milky Way's potential and account for dynamical friction. The key parameters are the proper motions and total masses of the dwarf galaxies. Using a simple model we find no tenable set of parameters that can allow Fornax to be consistent with filamentary infall, mainly because the 1 sigma error on its proper motion is relatively small. The other three must walk a tightrope between requiring a small pericentre (less than 20 kpc) to lose enough orbital energy to dynamical friction and avoiding being tidally disrupted. We then employed a more realistic model with host halo mass accretion, and found that the four dwarf galaxies must have fallen in at least 5 Gyrs ago. This time interval is longer than organised distribution is expected to last before being erased by the randomisation of the satellite orbits.

We present a new Particle-Mesh cosmological N-body code for accurately solving the modified Poisson equation of the Quasi Linear formulation of MOND. We generate initial conditions for the Angus (2009) cosmological model, which is identical to LCDM except that the cold dark matter is switched for a single species of thermal sterile neutrinos. We set the initial conditions at z=250 for a (512 Mpc/h)^3 box with 256^3 particles and we evolve them down to z=0. We clearly demonstrate the necessity of MOND for developing the large scale structure in a hot dark matter cosmology and contradict the naive expectation that MOND cannot form galaxy clusters. We find that the correct order of magnitude of X-ray clusters (with T_X > 4.5 keV) can be formed, but that we overpredict the number of very rich clusters and seriously underpredict the number of lower mass clusters. The latter is a shortcoming of the resolution of our simulations, whereas we suggest that the over production of very rich clusters might be prevented by incorporating a MOND acceleration constant that varies with redshift and an expansion history that cannot be described by the usual Friedmann models. We present evidence that suggests the density profiles of our simulated clusters are compatible with those of observed X-ray clusters in MOND. It remains to be seen if the low mass end of the cluster mass function can be reproduced and if the high densities of dark matter in the central 20 kpc of groups and clusters of galaxies, measured in the MOND framework, can be achieved. As a last test, we computed the relative velocity between pairs of halos within 10 Mpc and find that pairs with velocities larger than 3000 km/s like the bullet cluster, can form without difficulty.

Modified Newtonian Dynamics (MOND) has been shown to be able to fit spiral galaxy rotation curves as well as giving a theoretical foundation for empirically determined scaling relations, such as the Tully - Fisher law, without the need for a dark matter halo. As a complementary analysis, one should investigate whether MOND can also reproduce the dynamics of early - type galaxies (ETGs) without dark matter. As a first step, we here show that MOND can indeed fit the observed central velocity dispersion $\sigma_0$ of a large sample of ETGs assuming a simple MOND interpolating functions and constant anisotropy. We also show that, under some assumptions on the luminosity dependence of the Sersic n parameter and the stellar M/L ratio, MOND predicts a fundamental plane for ETGs : a log - linear relation among the effective radius $R_{eff}$, $\sigma_0$ and the mean effective intensity $\langle I_e \rangle$. However, we predict a tilt between the observed and the MOND fundamental planes.

Using the caustic method, we identify the member stars of five dwarf spheroidal (dSph) galaxies of the Milky Way, the smallest dark matter (DM) dominated systems in the Universe. After our interloper rejection, we compute line-of-sight velocity dispersion profiles that are substantially smoother than previous results. Moreover, two dSphs have line-of-sight velocity dispersions 20% smaller than previous calculations suggested. Our Jeans modelling confirms that the DM content interior to 300pc is roughly constant with satellite luminosity. Finally, if we assume that MOND provides the true law of gravity, our identification of interlopers implies that four dSphs have mass-to-light ratios in agreement with stellar population synthesis models, whereas Carina still has a mass-to-light ratio a factor of two too large and remains a problem for MOND.

We use dwarf spheroidal galaxies as a tool to attempt to put precise lower limits on the mass of the dark matter particle, assuming it is a sterile neutrino. We begin by making cored dark halo fits to the line of sight velocity dispersions as a function of projected radius (taken from Walker et al. 2007) for six of the Milky Way's dwarf spheroidal galaxies. We test Osipkov-Merritt velocity anisotropy profiles, but find that no benefit is gained over constant velocity anisotropy. In contrast to previous attempts, we do not assume any relation between the stellar velocity dispersions and the dark matter ones, but instead we solve directly for the sterile neutrino velocity dispersion at all radii by using the equation of state for a partially degenerate neutrino gas (which ensures hydrostatic equilibrium of the sterile neutrino halo). This yields a 1:1 relation between the sterile neutrino density and velocity dispersion, and therefore gives us an accurate estimate of the Tremaine-Gunn limit at all radii. By varying the sterile neutrino particle mass, we locate the minimum mass for all six dwarf spheroidals such that the Tremaine-Gunn limit is not exceeded at any radius (in particular at the centre). We find sizeable differences between the ranges of feasible sterile neutrino particle mass for each dwarf, but interestingly there exists a small range 270-280eV which is consistent with all dSphs at the 1-$\sigma$ level.

Certain types of globular clusters have the very important property that the predictions for their kinematics in the Newtonian and modified Newtonian dynamics (MOND) contexts are divergent. Here, we caution the recent claim that the stellar kinematics data (using 17 stars) of the globular cluster Palomar 14 are inconsistent with MOND. We compare the observations to the theoretical predictions using a Kolmogorov-Smirnov test, which is appropriate for small samples. We find that, with the currently available data, the MOND prediction for the velocity distribution can only be excluded with a very low confidence level, clearly insufficient to claim that MOND is falsified.

(Abridged) Modified Newtonian Dynamics (MOND) can fit a broad range of galaxy kinematic data, but struggles with clusters of galaxies. MONDian clusters need dark matter, and here we test the 11 eV sterile neutrino - used to fit the first three acoustic peaks of the cosmic microwave background - by investigating their equilibrium distributions in 30 groups and clusters over a wide range of temperatures. Knowing both the sterile neutrino density and velocity dispersion tells us the Tremaine-Gunn phase-space limit at all radii. We find that all 30 systems serendipitously reach the Tremaine-Gunn limit by the centre, which means a portion of the dynamical mass must always be covered by the brightest cluster galaxy. Interestingly, the typical fitted K-band mass-to-light ratio is unity and at most 1.2, which is very consistent - although leaving no margin for error - with stellar population synthesis models. Amidst the sample there are several special cases including the Coma cluster (for which dark matter was first proposed), NGC 720 (where geometrical evidence for dark matter was found) and the bullet cluster (where dark matter - of some kind - in clusters was directly proven to exist). We demonstrate that 11 eV sterile neutrinos are unlikely to influence spiral galaxy rotation curves, as they don't influence even some very massive early-types (NGC 4125 and NGC 6482). Finally, we conclude that it is intriguing that the minimum mass of sterile neutrino particle that can match the cosmic microwave background is identical to the minimum mass found here to be consistent with equilibrium configurations of MONDian clusters of galaxies.

We re-investigate the old problem of the survival of the five globular clusters orbiting the Fornax dwarf galaxy in both standard and modified Newtonian dynamics. For the first time in the history of the topic, we use accurate mass models for the Fornax dwarf, obtained through Jeans modelling of the recently published line of sight velocity dispersion data, and we are also not resigned to circular orbits for the globular clusters. Previously conceived problems stem from fixing the starting distances of the globulars to be less than half the tidal radius. We relax this constraint since there is absolutely no evidence for it and show that the dark matter paradigm, with either cusped or cored dark matter profiles, has no trouble sustaining the orbits of the two least massive globular clusters for a Hubble time almost regardless of their initial distance from Fornax. The three most massive globulars can remain in orbit as long as their starting distances are marginally outside the tidal radius. The outlook for modified Newtonian dynamics is also not nearly as bleak as previously reported. Although dynamical friction inside the tidal radius is far stronger in MOND, outside dynamical friction is negligible due to the absence of stars. This allows highly radial orbits to survive, but more importantly circular orbits at distances more than 85% of Fornax's tidal radius to survive indefinitely. The probability of the globular clusters being on circular orbits at this distance compared with their current projected distances is discussed and shown to be plausible. Finally, if we ignore the presence of the most massive globular (giving it a large line of sight distance) we demonstrate that the remaining four globulars can survive within the tidal radius for the Hubble time with perfectly sensible orbits.

If a single sterile neutrino exists such that $m_{\nu_s}\sim11eV$, it can serendipitously solve all outstanding issues of the Modified Newtonian Dynamics. With it one can explain the dark matter of galaxy clusters without influencing individual galaxies, match the angular power spectrum of the cosmic microwave background and potentially fit the matter power spectrum. This model is flat with $\Omega_{\nu_s}\sim0.23$ and the usual baryonic and dark energy components, thus the Universe has the same expansion history as the $\lcdm$ model and only differs at the galactic scale where the Modified Dynamics outperforms $\lcdm$ significantly.

We take the line of sight velocity dispersions as functions of radius for 8 Milky Way dwarf spheroidal galaxies and use Jeans analysis to calculate the mass-to-light ratios (M/L) in Modified Newtonian Dynamics (MOND). Using the latest structural parameters, distances and variable velocity anisotropy, we find 6/8 dwarfs have sensible M/L using only the stellar populations. Sextans and Draco, however, have M/L=9.2_{-3.0}^{+5.3} and 43.9_{-19.3}^{+29.0} respectively, which poses a problem. Apart from the need for Sextans' integrated magnitude to be reviewed, we propose tidal effects intrinsic to MOND, testable with numerical simulations, but fully orbit dependant, which are disrupting Draco. The creation of the Magellanic Stream is also re-addressed in MOND, the scenario being the stream is ram pressure stripped from the SMC as it crosses the LMC.

Although very successful in explaining the observed conspiracy between the baryonic distribution and the gravitational field in spiral galaxies without resorting to dark matter (DM), the modified Newtonian dynamics (MOND) paradigm still requires DM in X-ray bright systems. Here, to get a handle on the distribution and importance of this DM, and thus on its possible form, we deconstruct the mass profiles of 26 X-ray emitting systems in MOND, with temperatures ranging from 0.5 to 9 keV. Initially we compute the MOND dynamical mass as a function of radius, then subtract the known gas mass along with a component of galaxies which includes the cD galaxy with $M/L_K=1$. Next we test the compatibility of the required DM with ordinary massive neutrinos at the experimental limit of detection ($m_{\nu}=2$ eV), with density given by the Tremaine-Gunn limit. Even by considering that the neutrino density stays constant and maximal within the central 100 or 150 kpc (which is the absolute upper limit of a possible neutrino contribution there), we show that these neutrinos can never account for the required DM within this region. The natural corollary of this finding is that, whereas clusters (T $\ga$ 3 keV) might have most of their mass accounted for if ordinary neutrinos have a 2 eV mass, groups (T $\lsim$ 2 keV) cannot be explained by a 2 eV neutrino contribution. This means that, for instance, cluster baryonic dark matter (CBDM, Milgrom 2007) or even sterile neutrinos would present a more satisfactory solution to the problem of missing mass in MOND X-ray emitting systems.

We summarize the status of Modified Newtonian Dynamics (MOND) in galaxy clusters. The observed acceleration is typically larger than the acceleration threshold of MOND in the central regions, implying that some dark matter is necessary to explain the mass discrepancy there. A plausible resolution of this issue is that the unseen mass in MOND is in the form of ordinary neutrinos with masses just below the experimentally detectable limit. In particular, we show that the lensing mass reconstructions of the clusters 1E0657-56 (the bullet cluster) and Cl0024+17 (the ring) do not pose a new challenge to this scenario. However, the mass discrepancy for cool X-ray emitting groups, in which neutrinos cannot cluster, pose a more serious problem, meaning that dark baryons could present a more satisfactory solution to the problem of unseen mass in MOND clusters.

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

We consider the orbit of the bullet cluster 1E 0657-56 in both CDM and MOND using accurate mass models appropriate to each case in order to ascertain the maximum plausible collision velocity. Impact velocities consistent with the shock velocity (~ 4700km/s) occur naturally in MOND. CDM can generate collision velocities of at most ~ 3800km/s, and is only consistent with the data provided that the shock velocity has been substantially enhanced by hydrodynamical effects.

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 develop a new method to predict the density associated with weak lensing maps of (un)relaxed clusters in a range of theories interpolating between GR and MOND (General Relativity and Modified Newtonian Dynamics). We apply it to fit the lensing map of the bullet merging cluster 1E0657-56, in order to constrain more robustly the nature and amount of collisionless matter in clusters {\it beyond} the usual assumption of spherical equilibrium (Pointecouteau & Silk 2005) and the validity of GR on cluster scales (Clowe et al. 2006). Strengthening the proposal of previous authors we show that the bullet cluster is dominated by a collisionless -- most probably non-baryonic -- component in GR as well as in MOND, a result consistent with the dynamics of many X-ray clusters. Our findings add to the number of known pathologies for a purely baryonic MOND, including its inability to fit the latest data from the Wilkinson Microwave Anisotropy Probe. A plausible resolution of all these issues and standard issues of Cold Dark Matter with galaxy rotation curves is the "marriage" of MOND with ordinary hot neutrinos of 2eV. This prediction is just within the GR-independent maximum of neutrino mass from current $\beta$-decay experiments, and is falsifiable by the Karlsruhe Tritium Neutrino (KATRIN) experiment by 2009. Issues of consistency with strong lensing arcs and the large relative velocity of the two clusters comprising the bullet cluster are also addressed.

We point out that the usual self-similarity in cold dark matter models is broken by encounters with individual normal galactic stars on sub-pc scale. Tidal heating and stripping must have redefined the density and velocity structures of the population of the Earth-mass dark matter halos, which are likely to have been the first bound structures to form in the Universe. The disruption rate depends strongly on {\it galaxy types} and the orbital distribution of the microhalos; in the Milky Way, stochastic radial orbits are destroyed first by stars in the triaxial bulge, microhalos on non-planar retrograde orbits with large pericenters and/or apocenters survive the longest. The final microhalo distribution in the {\it solar neighborhood} is better described as a superposition of filamentry microstreams rather than as a set of discrete spherical clumps in an otherwise homogeneous medium. We discuss its important consequences to our detections of microhalos by direct recoil signal and indirect annihilation signal.

A special purpose N-body simulation has been built to understand the tidal heating of the smallest dark matter substructures ($10^{-6}\msun$ and 0.01pc) from the grainy potential of the Milky Way due to individual stars in the disk and the bulge. To test the method we first run simulations of single encounters of microhalos with an isolated star, and compare with analytical predictions of the dark particle bound fraction as a function of impact parameter. We then follow the orbits of a set of microhalos in a realistic flattened Milky Way potential. We concentrate on (detectable) microhalos passing near the Sun with a range of pericenter and apocenter. Stellar perturbers near the orbital path of a microhalo would exert stochachstic impulses, which we apply in a Monte Carlo fashion according to the Besancon model for the distribution of stars of different masses and ages in our Galaxy. Also incorporated are the usual pericenter tidal heating and disk-shocking heating. We give a detailed diagnosis of typical microhalos and find microhalos with internal tangential anisotropy are slightly more robust than the ones with radial anisotropy. In addition, the dark particles generally go through of a random walk in velocity space and diffuse out of the microhalos. We show that the typical destruction time scales are strongly correlated with the stellar density averaged along a microhalo's orbit over the age of the stellar disk. We also present the morphology of a microhalo at several epochs which may hold the key to dark matter detections.

A proper test of Modified Newtonian Dynamics (MOND) in systems of non-trivial geometries depends on modelling subtle differences in several versions of its postulated theories. This is especially important for lensing and dynamics of barely virialised galaxy clusters with typical gravity of scale $\sim a_0 \sim 1\AA{\rm s}^{-2}$. The original MOND formula, the classical single field modification of the Poisson equation, and the multi-field general relativistic theory of Bekenstein (TeVeS) all lead to different predictions as we stray from spherical symmetry. In this paper, we study a class of analytical MONDian models for a system with a semi-Hernquist baryonic profile. After presenting the analytical distribution function of the baryons in spherical limits, we develop orbits and gravitational lensing of the models in non-spherical geometries. In particular, we can generate a multi-centred baryonic system with a weak lensing signal resembling that of the merging galaxy cluster 1E 0657-56 with a bullet-like light distribution. We finally present analytical scale-free highly non-spherical models to show the subtle differences between the single field classical MOND theory and the multi-field TeVeS theory.