
We determine the escape velocity from the Milky Way (MW) at a range of
Galactocentric radii in the context of Modified Newtonian Dynamics (MOND). Due
to its nonlinear nature, escape is possible if the MW is considered embedded
in a constant external gravitational field (EF) from distant objects. We model
this situation using a fully selfconsistent method based on a direct solution
of the governing equations out to several thousand disk scale lengths. We try
out a range of EF strengths and mass models for the MW in an attempt to match
the escape velocity measurements of Williams et al. (2017).
A reasonable match is found if the EF on the MW is ${\sim 0.03 a_{_0}}$,
towards the higher end of the range considered. Our models include a hot gas
corona surrounding the MW, but our results suggest that this should have a very
low mass of ${\sim 2 \times 10^{10} M_\odot}$ to avoid pushing the escape
velocity too high. Our analysis favours a slightly lower baryonic disk mass
than the ${\sim 7 \times 10^{10} M_\odot}$ required to explain its rotation
curve in MOND. However, given the uncertainties, MOND is consistent with both
the locally measured amplitude of the MW rotation curve and its escape velocity
over Galactocentric distances of 8$$50 kpc.

The standard $\Lambda$CDM paradigm seems to describe cosmology and large
scale structure formation very well. However, a number of puzzling observations
remain on galactic scales. An example is the anisotropic distribution of
satellite galaxies in the Local Group. This has led to suggestions that a
modified gravity theory might provide a better explanation than Newtonian
gravity supplemented by dark matter. One of the leading modified gravity
theories is Modified Newtonian Dynamics (MOND). For an isolated point mass, it
boosts gravity by an accelerationdependent factor of $\nu$.
Recently, a much more computerfriendly quasilinear formulation of MOND
(QUMOND) has become available. We investigate analytically the solution for a
point mass embedded in a constant external field of $\mathbf{g}_{ext}$. We find
that the potential is $\Phi =  ~ \frac{GM \nu_{ext}}{r}\left(1 + \frac{K_0}{2}
\sin^2 \theta \right)$, where $r$ is distance from the mass $M$ which is in an
external field that `saturates' the $\nu$ function at the value $\nu_{ext}$,
leading to a fixed value of $K_0 \equiv \frac{\partial Ln ~ \nu}{\partial Ln ~
g_{ext}}$. In a very weak gravitational field $\left(\left \mathbf{g}_{ext}
\right \ll a_0 \right)$, $K_0 = \frac{1}{2}$. The angle $\theta$ is that
between the external field direction and the direction towards the mass.
Our results are quite close to the more traditional aquadratic Lagrangian
(AQUAL) formulation of MOND. We apply both theories to a simple model of the
Sagittarius tidal stream. We find that they give very similar results, with the
tidal stream seeming to spread slightly further in AQUAL.

We attempt to understand the planes of satellite galaxies orbiting the Milky
Way (MW) and M31 in the context of Modified Newtonian Dynamics (MOND), which
implies a close MWM31 flyby occurred ${\approx 8}$ Gyr ago. Using the timing
argument, we obtain MWM31 trajectories consistent with cosmological initial
conditions and present observations. We adjust the present M31 proper motion
within its uncertainty in order to simulate a range of orbital geometries and
closest approach distances. Treating the MW and M31 as point masses, we follow
the trajectories of surrounding test particle disks, thereby mapping out the
tidal debris distribution.
Around each galaxy, the resulting tidal debris tends to cluster around a
particular orbital pole. We find some models in which these preferred spin
vectors align fairly well with those of the corresponding observed satellite
planes. The radial distributions of material in the simulated satellite planes
are similar to what we observe. Around the MW, our bestfitting model yields a
significant fraction (0.22) of counterrotating material, perhaps explaining
why Sculptor counterrotates within the MW satellite plane. In contrast, our
model yields no counterrotating material around M31. This is testable with
proper motions of M31 satellites.
In our best model, the MW disk is thickened by the flyby 7.65 Gyr ago to a
root mean square height of 0.75 kpc. This is similar to the observed age and
thickness of the Galactic thick disk. Thus, the MW thick disk may have formed
together with the MW and M31 satellite planes during a past MWM31 flyby.

We recently showed that several Local Group (LG) galaxies have much higher
radial velocities (RVs) than predicted by a 3D dynamical model of the standard
cosmological paradigm. Here, we show that 6 of these 7 galaxies define a thin
plane with root mean square thickness of only 101 kpc despite a widest extent
of nearly 3 Mpc, much larger than the conventional virial radius of the Milky
Way (MW) or M31. This plane passes within ${\sim 70}$ kpc of the MWM31
barycentre and is oriented so the MWM31 line is inclined by $16^\circ$ to it.
We develop a toy model to constrain the scenario whereby a past MWM31 flyby
in Modified Newtonian Dynamics (MOND) forms tidal dwarf galaxies that settle
into the recently discovered planes of satellites around the MW and M31. The
scenario is viable only for a particular MWM31 orbital plane. This roughly
coincides with the plane of LG dwarfs with anomalously high RVs.
Using a restricted $N$body simulation of the LG in MOND, we show how the
once fastmoving MW and M31 gravitationally slingshot test particles outwards
at high speeds. The most distant such particles preferentially lie within the
MWM31 orbital plane, probably because the particles ending up with the highest
RVs are those flung out almost parallel to the motion of the perturber. This
suggests a dynamical reason for our finding of a similar trend in the real LG,
something not easily explained as a chance alignment of galaxies with an
isotropic or mildly flattened distribution (probability $= {0.0015}$).

It has recently been proposed, by assuming that dark matter is a superfluid,
that MONDlike effects can be achieved on small scales whilst preserving the
success of $\Lambda$CDM on large scales. Here we aim to provide the first set
of spherical models of galaxy clusters in the context of superfluid dark
matter. We first outline the theoretical structure of the superfluid core and
the surrounding "normal phase" dark halo of quasiparticles in thermal
equlibrium. The latter should encompass the largest part of galaxy clusters.
Here, we set the SfDM transition at the radius where the density and pressure
of the superfluid and normal phase coincides, neglecting the effect of phonons
in the suprefluid core. We then apply the theory to a sample of galaxy
clusters, and directly compare the SfDM predicted mass profiles to data. We
find that the superfluid formulation can reproduce the Xray dynamical mass
profile of clusters. The SfDM fits however display slight underpredictions of
the gravity in the central regions which might be partly related to our
neglecting of the effect of phonons in these regions. We conclude that this
superfluid formulation is successful in describing galaxy clusters, but further
work will be needed to determine whether the parameter choice is consistent
with galaxies. Our model could be made more realistic by exploring
nonsphericity and the SfDM transition condition we impose.

Adopting Schwarzschild's orbitsuperposition technique, we construct a series
of selfconsistent galaxy models, embedded in the external field of galaxy
clusters in the framework of Milgrom's MOdified Newtonian Dynamics. These
models represent relatively massive ellipticals with a Hernquist radial profile
at various distances from the cluster centre. Using $N$body simulations, we
perform a first analysis of these models and their evolution. We find that
selfgravitating axisymmetric density models, even under a weak external field,
lose their symmetry by instability and generally evolve to triaxial
configurations. A kinematic analysis suggests that the instability originates
from both box and nonclassified orbits with low angular momentum. We also
consider a selfconsistent isolated system which is then placed in a strong
external field and allowed to evolve freely. This model, just as the
corresponding equilibrium model in the same external field, eventually settles
to a triaxial equilibrium as well, but has a higher velocity radial anisotropy
and is rounder. The presence of an external field in MOND universe generically
predicts some lopsidedness of galaxy shapes.

We model the cluster A1689 in two modified MOND frameworks (EMOND and what we
call generalised MOND or GMOND) with the aim of determining whether it is
possible to explain the inferred acceleration profile, from gravitational
lensing, without the aid of dark matter. We also compare our result to
predictions from MOG/STVG and Emergent Gravity. By using a baryonic mass model,
we determine the total gravitational acceleration predicted by the modified
gravitational equations and compare the result to NFW profiles of dark matter
studies on A1689 from the literature. Theory parameters are inferred
empirically, with the aid of previous work. We are able to reproduce the
desired acceleration profile of A1689 for GMOND, EMOND and MOG, but not
Emergent Gravity. There is much more work which needs to be conducted with
regards to understanding how the GMOND parameters behave in different
environments. Furthermore, we show that the exact baryonic profile becomes very
important when undertaking modified gravity modelling rather than {\Lambda}CDM.

A sample of Coma cluster ultradiffuse galaxies (UDGs) are modelled in the
context of Extended Modified Newtonian Dynamics (EMOND) with the aim to explain
the large dark matterlike effect observed in these cluster galaxies.
We first build a model of the Coma cluster in the context of EMOND using gas
and galaxy mass profiles from the literature. Then assuming the dynamical mass
of the UDGs satisfies the fundamental manifold of other ellipticals, and that
the UDG stellar masstolight matches their colour, we can verify the EMOND
formulation by comparing two predictions of the baryonic mass of UDGs.
We find that EMOND can explain the UDG mass, within the expected modelling
errors, if they lie on the fundamental manifold of ellipsoids, however, given
that measurements show one UDG lying off the fundamental manifold, observations
of more UDGs are needed to confirm this assumption.

We attempt to fit the observed radial velocities (RVs) of ${\,{\sim}\,}30$
Local Group (LG) galaxies using a 3D dynamical model of it and its immediate
environment within the context of the standard cosmological paradigm,
$\Lambda$CDM. This extends and confirms the basic results of our previous
axisymmetric investigation of the LG (MNRAS, 459, 2237). We find that there
remains a tendency for observed RVs to exceed those predicted by our
bestfitting model. The typical mismatch is slightly higher than in our 2D
model, with a root mean square value of ${\,{\sim}\,}50$ km/s. Our main finding
is that including the 3D distribution of massive perturbing dark matter halos
is unlikely to help greatly with the high velocity galaxy problem. Nonetheless,
the 2D and 3D results differ in several other ways such as which galaxies' RVs
are most problematic and the preferred values of parameters common to both
models.
The anomalously high RVs of several LG dwarfs may be better explained if the
Milky Way (MW) and Andromeda (M31) were once moving much faster than in our
models. This would allow LG dwarfs to gain very high RVs via gravitational
slingshot encounters with a massive fastmoving galaxy. Such a scenario is
possible in some modified gravity theories, especially those which require the
MW and M31 to have previously undergone a close flyby. In a $\Lambda$CDM
context, however, this scenario is not feasible as the resulting dynamical
friction would cause a rapid merger.

Modified Newtonian Dynamics (MOND) is a gravitational framework designed to
explain the astronomical observations in the Universe without the inclusion of
particle dark matter. Modified Newtonian Dynamics, in its current form, cannot
explain the missing mass in galaxy clusters without the inclusion of some extra
mass, be it in the form of neutrinos or nonluminous baryonic matter. We
investigate whether the MOND framework can be generalized to account for the
missing mass in galaxy clusters by boosting gravity in high gravitational
potential regions. We examine and review Extended MOND (EMOND), which was
designed to increase the MOND scale acceleration in high potential regions,
thereby boosting the gravity in clusters. We seek to investigate galaxy cluster
mass profiles in the context of MOND with the primary aim at explaining the
missing mass problem fully without the need for dark matter. Using the
assumption that the clusters are in hydrostatic equilibrium, we can compute the
dynamical mass of each cluster and compare the result to the predicted mass of
the EMOND formalism. We find that EMOND has some success in fitting some
clusters, but overall has issues when trying to explain the mass deficit fully.
We also investigate an empirical relation to solve the cluster problem, which
is found by analysing the cluster data and is based on the MOND paradigm. We
discuss the limitations in the text.

Counter to intuition, the images of an extended galaxy lensed by a moving
galaxy cluster should have slightly different spectra in any metric gravity
theory. This is mainly for two reasons. One relies on the gravitational
potential of a moving lens being timedependent (the $\text{Moving}$
$\text{Cluster}$ $\text{Effect}$, $\text{MCE}$). The other is due to uneven
magnification across the extended, rotating source (the $\text{Differential}$
$\text{Magnification}$ $\text{Effect}$, $\text{DME}$). The time delay between
the images can also cause their redshifts to differ because of cosmological
expansion. This Differential Expansion Effect is likely to be small. Using a
simple model, we derive these effects from first principles.
One application would be to the Bullet Cluster, whose large tangential
velocity may be inconsistent with the $\Lambda CDM$ paradigm. This velocity can
be estimated with complicated hydrodynamic models. Uncertainties with such
models can be avoided using the MCE. We argue that the MCE should be observable
with ALMA.
However, such measurements can be corrupted by the DME if typical spiral
galaxies are used as sources. Fortunately, we find that if detailed spectral
line profiles were available, then the DME and MCE could be distinguished. It
might also be feasible to calculate how much the DME should affect the mean
redshift of each image. Resolved observations of the source would be required
to do this accurately.
The DME is of order the source angular size divided by the Einstein radius
times the redshift variation across the source. Thus, it mostly affects nearly
edgeon spiral galaxies in certain orientations. This suggests that observers
should reduce the DME by careful choice of target, a possibility we discuss in
some detail.

The positions and velocities of galaxies in the Local Group (LG) measure the
gravitational field within it. This is mostly due to the Milky Way (MW) and
Andromeda (M31). We constrain their masses using distance and radial velocity
(RV) measurements of 32 LG galaxies. To do this, we follow the trajectories of
many simulated particles starting on a pure Hubble flow at redshift 9. For each
observed galaxy, we obtain a trajectory which today is at the same position.
Its final velocity is the model prediction for the velocity of that galaxy.
Unlike previous simulations based on spherical symmetry, ours are
axisymmetric and include gravity from Centaurus A. We find the total LG mass is
${4.33^{+0.37}_{0.32}\times{10}^{12}M_\odot}$, with $0.14 \pm 0.07$ of this
being in the MW. We approximately account for IC 342, M81, the Great Attractor
and the Large Magellanic Cloud.
No plausible set of initial conditions yields a good match to the RVs of our
sample of LG galaxies. Observed RVs systematically exceed those predicted by
the bestfitting $\Lambda$CDM model, with a typical disagreement of
${45.1^{+7.0}_{5.7}}$ km/s and a maximum of ${110 \pm 13}$ km/s for DDO 99.
Interactions between LG dwarf galaxies can't easily explain this.
One possibility is a past close flyby of the MW and M31. This arises in some
modified gravity theories but not in $\Lambda$CDM. Gravitational slingshot
encounters of material in the LG with either of these massive fastmoving
galaxies could plausibly explain why some nonsatellite LG galaxies are moving
away from us even faster than a pure Hubble flow.

The TMT Detailed Science Case describes the transformational science that the
Thirty Meter Telescope will enable. Planned to begin science operations in
2024, TMT will open up opportunities for revolutionary discoveries in
essentially every field of astronomy, astrophysics and cosmology, seeing much
fainter objects much more clearly than existing telescopes. Per this
capability, TMT's science agenda fills all of space and time, from nearby
comets and asteroids, to exoplanets, to the most distant galaxies, and all the
way back to the very first sources of light in the Universe.
More than 150 astronomers from within the TMT partnership and beyond offered
input in compiling the new 2015 Detailed Science Case. The contributing
astronomers represent the entire TMT partnership, including the California
Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA),
the National Astronomical Observatories of the Chinese Academy of Sciences
(NAOC), the National Astronomical Observatory of Japan (NAOJ), the University
of California, the Association of Canadian Universities for Research in
Astronomy (ACURA) and US associate partner, the Association of Universities for
Research in Astronomy (AURA).

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.

Following earlier authors, we reexamine constraints on the radial velocity
anisotropy of generic stellar systems using arguments for phase space density
positivity, stability, and separability. It is known that although the majority
of commonly used systems have an maximum anisotropy of less than half of the
logarithmic density slope \emph{i.e.} $\beta < \gamma/2$, there are exceptions
for separable models with large central anisotropy. Here we present a new
exceptional case with abovethreshold anisotropy locally but with an isotropic
center nevertheless. These models are nonseparable and we maintain positivity.
Our analysis suggests that regions of abovethreshold anisotropy are more
related to regions of possible secular instability, which might be observed in
selfconsistent galaxies in a shortlived phase. s an equilibrium solution for
a collisionless system.

We study the evolution of the phasespace of collisionless Nbody systems
under repeated stirrings or perturbations, which has been shown to lead to a
convergence towards a limited group of end states. This socalled attractor was
previously shown to be independent of the initial system and environmental
conditions. However the fundamental reason for its appearance is still unclear.
It has been suggested that the origin of the attractor may be either radial
infall (RI), the radial orbit instability (ROI), or energy exchange which, for
instance, happens during violent relaxation. Here we examine the effects of a
set of controlled perturbations, referred to as `kicks', which act in addition
to the standard collisionless dynamics by allowing prespecified instantaneous
perturbations in phasespace. We first demonstrate that the attractor persists
in the absence of RI and ROI by forcing the system to expand. We then consider
radial velocity kicks in a rigid potential and isotropic velocity kicks, since
there are no energy exchanges in these two recipes of kicks. We find that these
kicks do not lead to the attractor, indicating that the energy exchange in a
dynamic potential is indeed the physical mechanism responsible for the
attractor.

There is growing interest in testing alternative gravity theories using the
subtle gravitational redshifts in clusters of galaxies. However, current models
all neglect a transverse Doppler redshift of similar magnitude, and some models
are not selfconsistent. An equilibrium model would fix the gravitational and
transverse Doppler velocity shifts to be about 6sigma^2/c and 3sigma^2/2c in
order to fit the observed velocity dispersion sigma selfconsistently. This
result comes from the Virial Theorem for a spherical isotropic cluster, and is
insensitive to the theory of gravity. A gravitational redshift signal also does
not directly distinguish between the Einsteinian and f(R) gravity theories,
because each theory requires different dark halo mass function to keep the
clusters in equilibrium. When this constraint is imposed, the gravitational
redshift has no sensitivity to theory. Indeed our Nbody simulations show that
the halo mass function differs in f(R), and that the transverse Doppler effect
is stronger than analytically predicted due to nonequilibrium.

The Local Group timing has been one of the first historical probes of the
missing mass problem. Whilst modern cosmological probes indicate that pure
baryonic dynamics is not sufficient on the largest scales, nearby galaxies and
small galaxy groups persistently obey Milgrom's MOND law, which implies that
dynamics at small scales is possibly entirely predicted by the baryons. Here,
we investigate the Local Group timing in this context of Milgromian dynamics.
Making use of the latest measured proper motions and radial velocities for
Andromeda and the Magellanic clouds, we integrate backwards their orbits by
making use of the Milgromian twobody equation of motion. We find that, with
the currently measured proper motions and radial velocity of M31, MOND would
imply that the Milky Way and M31 first moved apart via Hubble expansion after
birth, but then necessarily got attracted again by the Milgromian gravitational
attraction, and had a past flyby encounter before coming to their present
positions. This encounter would most probably have happened 7 to 11 Gyr ago
(0.8<z<3). The absence of a dark matter halo and its associated dynamical
friction is necessary for such a close encounter not to trigger a merger.
Observational arguments which could exclude or favour such a past encounter
would thus be very important in view of falsifying or vindicating Milgromian
dynamics on the scale of the Local Group. Interestingly, the closest approach
of the encounter is small enough (<55 kpc) to have severe consequences on the
disk dynamics, including perhaps thick disk formation, and on the satellite
systems of both galaxies. The ages of the satellite galaxies and of the young
halo globular clusters, all of which form the vast polar structure around the
Milky Way, are consistent with these objects having been born in this
encounter.

We use Schwarzschild's orbitsuperposition technique to construct
selfconsistent models of the Galactic bar. Using $\chi^2$ minimisation, we
find that the bestfit Galactic bar model has a pattern speed $\Omega_{\rm
p}=60 \rm{km s^{1} kpc^{1}}$, disk mass $\rm{M_{\rm
d}=1.0\times10^{11}M_{\odot}}$ and bar angle $\theta_{\rm bar}=20^{\circ}$ for
an adopted bar mass $\rm{M_{\rm bar}=2\times10^{10}M_{\odot}}$. The model can
reproduce not only the threedimensional and projected density distributions
but also velocity and velocity dispersion data from the BRAVA survey. We also
predict the proper motions in the range $l=[12^{\circ},12^{\circ}]$,
$b=[10^{\circ},10^{\circ}]$, which appear to be higher than observations in
the longitudinal direction. The model is stable within a timescale of 0.5 Gyr,
but appears to deviate from steadystate on longer timescales. Our model can be
further tested by future observations such as those from GAIA.

A remarkably tight relation is observed between the Newtonian gravity sourced
by the baryons and the actual gravity in galaxies of all sizes. This can be
interpreted as the effect of a single effective force law depending on
acceleration. This is however not the case in larger systems with much deeper
potential wells, such as galaxy clusters. Here we explore the possibility of an
effective force law reproducing mass discrepancies in all extragalactic systems
when depending on both acceleration and the depth of the potential well. We
exhibit, at least at a phenomenological level, one such possible construction
in the classical gravitational potential theory. Interestingly, the framework,
dubbed EMOND, is able to reproduce the observed mass discrepancies in both
galaxies and galaxy clusters, and to produce multicenter systems with offsets
between the peaks of gravity and the peaks of the baryonic distribution.

A unique signature of the modified Newtonian dynamics (MOND) paradigm is its
peculiar behavior in the vicinity of the points where the total Newtonian
acceleration exactly cancels. In the Solar System, these are the saddle points
of the gravitational potential near the planets. Typically, such points are
embedded into lowacceleration bubbles where modified gravity theories a` la
MOND predict significant deviations from Newton's laws. As has been pointed out
recently, the EarthSun bubble may be visited by the LISA Pathfinder spacecraft
in the near future, providing a unique occasion to put these theories to a
direct test. In this work, we present a highprecision model of the Solar
System's gravitational potential to determine accurate positions and motions of
these saddle points and study the predicted dynamical anomalies within the
framework of quasilinear MOND. Considering the expected sensitivity of the
LISA Pathfinder probe, we argue that interpolation functions which exhibit a
"faster" transition between the two dynamical regimes have a good chance of
surviving a null result. An example of such a function is the QMOND analog of
the socalled simple interpolating function which agrees well with much of the
extragalactic phenomenology. We have also discovered that several of Saturn's
outermost satellites periodically intersect the SaturnSun bubble, providing
the first example of Solar System objects that regularly undergo the MOND
regime.

We study the evolution of the phasespace of collisionless Nbody systems
under repeated stirrings or perturbations. We find convergence towards a
limited solution group, in accordance with Hansen 2010, that is independent of
the initial system and environmental conditions, paying particular attention to
the assumed gravitational paradigm (Newtonian and MOND). We examine the effects
of changes to the perturbation scheme and in doing so identify a large group of
perturbations featuring radial orbit instability (ROI) which always lead to
convergence. The attractor is thus found to be a robust and reproducible effect
under a variety of circumstances.

We present new radial velocity measurements from the Bulge Radial Velocity
Assay (BRAVA), a large scale spectroscopic survey of Mtype giants in the
Galactic bulge/bar region. The sample of ~4500 new radial velocities, mostly in
the region 10 deg < l < +10 deg and b ~ 6 deg more than doubles the existent
published data set. Our new data extend our rotation curve and velocity
dispersion profile to +20 deg, which is ~2.8 kpc from the Galactic Center. The
new data confirm the cylindrical rotation observed at 6 deg and 8 deg, and
are an excellent fit to the Shen et al. (2010) Nbody bar model. We measure the
strength of the TiO molecular band as a first step towards a metallicity
ranking of the stellar sample, from which we confirm the presence of a vertical
abundance gradient. Our survey finds no strong evidence of previously unknown
kinematic streams. We also publish our complete catalog of radial velocities,
photometry, TiO band strengths, and spectra, which is available at the IRSA
archive: http://irsa.ipac.caltech.edu/ as well as at UCLA:
http://brava.astro.ucla.edu/.

Certain covariant theories of the modified Newtonian dynamics paradigm seem
to require an additional hot dark matter (HDM) component  in the form of
either heavy ordinary neutrinos or more recently light sterile neutrinos (SNs)
with a mass around 11eV  to be relieved of problems ranging from cosmological
scales down to intermediate ones relevant for galaxy clusters. Here we suggest
using gravitational lensing by galaxy clusters to test such a marriage of
neutrino HDM and modified gravity, adopting the framework of
tensorvectorscalar theory (TeVeS). Unlike conventional cold dark matter
(CDM), such HDM is subject to strong phasespace constraints, which allows one
to check cluster lens models inferred within the modified framework for
consistency. Since the considered HDM particles cannot collapse into
arbitrarily dense clumps and only form structures well above the galactic
scale, systems which indicate the need for dark substructure are of particular
interest. As a first example, we study the cluster lens Abell 2390 and its
impressive straight arc with the help of numerical simulations. Based on our
results, we outline a general and systematic approach to model cluster lenses
in TeVeS which significantly reduces the calculation complexity. We further
consider a simple bimodal lens configuration, capable of producing the straight
arc, to demonstrate our approach. We find that such a model is marginally
consistent with the hypothesis of 11eV SNs. Future work including more detailed
and realistic lens models may further constrain the necessary SN distribution
and help to conclusively assess this point. Cluster lenses could therefore
provide an interesting discriminator between CDM and such modified gravity
scenarios supplemented by SNs or other choices of HDM.

Using strong lensing data Milgrom's MOdified Newtonian Dynamics (MOND) or its
covariant TeVeS (TensorVectorScalar Theory) is being examined here as an
alternative to the conventional $\Lambda$CDM paradigm. We examine 10
doubleimage gravitational lensing systems, in which the lens masses have been
estimated by stellar population synthesis models. While mild deviations exist,
we do not find out that strong cases for outliers to the TeVeS theory.