The covariant scale-invariant dynamics (SID) theory has recently been proposed as a possible explanation for the observed dynamical discrepancies in galaxies (Maeder & Gueorguiev 2020). SID implies that these discrepancies $-$ commonly attributed to dark matter $-$ arise instead from a non-standard velocity-dependent force that causes two-body near-Keplerian orbits to expand. We show that the predicted expansion of the Earth-Moon orbit is incompatible with lunar laser ranging data at $>200\sigma$. Moreover, SID predicts that the gravitating mass of any object was much smaller in the past. If true, a low-mass red giant star must be significantly older than in standard theory. This would make it much older than the conventional age of the Universe, which however is expected to be similarly old in SID. Moreover, it is not completely clear whether SID truly contains new physics beyond General Relativity, with several previous works arguing that the extra degree of freedom is purely mathematical. We conclude that the SID model is falsified at high significance by observations across a range of scales, even if it is theoretically well formulated.

The gravitational dipole theory of Hadjukovic (2010) is based on the hypothesis that antimatter has a negative gravitational mass and thus falls upwards on Earth. Astrophysically, the model is similar to but more fundamental than Modified Newtonian Dynamics (MOND), with the Newtonian gravity $g_{_N}$ towards an isolated point mass boosted by the factor $\nu = 1 + \left( \alpha/x \right) \tanh \left( \sqrt{x}/\alpha \right)$, where $x \equiv g_{_N}/a_{_0}$ and $a_{_0} = 1.2 \times 10^{-10}$ m/s$^2$ is the MOND acceleration constant. We show that $\alpha$ must lie in the range ${0.4-1}$ to acceptably fit galaxy rotation curves. In the Solar System, this interpolating function implies an extra Sunwards acceleration of ${\alpha a_{_0}}$. This would cause Saturn to deviate from Newtonian expectations by ${7000 \left( \alpha/0.4 \right)}$ km over 15 years, starting from known initial position and velocity on a near-circular orbit. We demonstrate that this prediction should not be significantly altered by the postulated dipole haloes of other planets due to the rather small region in which each planet's gravity dominates over that of the Sun. The orbit of Saturn should similarly be little affected by a possible ninth planet in the outer Solar System and by the Galactic gravity causing a non-spherical distribution of gravitational dipoles several kAU from the Sun. Radio tracking of the Cassini spacecraft orbiting Saturn yields a ${5\sigma}$ upper limit of 160 metres on deviations from its conventionally calculated trajectory. These measurements imply a much more stringent upper limit on $\alpha$ than the minimum required for consistency with rotation curve data. Therefore, no value of $\alpha$ can simultaneously match all available constraints, falsifying the gravitational dipole theory in its current form at extremely high significance.

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 acceleration-dependent factor of $\nu$. Recently, a much more computer-friendly quasi-linear 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 consider the feasibility of testing Newtonian gravity at low accelerations using wide binary (WB) stars separated by $\ge 3$ kAU. These systems probe the accelerations at which galaxy rotation curves unexpectedly flatline, possibly due to Modified Newtonian Dynamics (MOND). We conduct Newtonian and MOND simulations of WBs covering a grid of model parameters in the system mass, semi-major axis, eccentricity and orbital plane. We self-consistently include the external field (EF) from the rest of the Galaxy on the Solar neighbourhood using an axisymmetric algorithm. For a given projected separation, WB relative velocities reach larger values in MOND. The excess is ${\approx 20\%}$ adopting its simple interpolating function, as works best with a range of Galactic and extragalactic observations. This causes noticeable MOND effects in accurate observations of ${\approx 500}$ WBs, even without radial velocity measurements. We show that the proposed Theia mission may be able to directly measure the orbital acceleration of Proxima Centauri towards the 13 kAU-distant $\alpha$ Centauri. This requires an astrometric accuracy of $\approx 1 \, \mu$as over 5 years. We also consider the long-term orbital stability of WBs with different orbital planes. As each system rotates around the Galaxy, it experiences a time-varying EF because this is directed towards the Galactic Centre. We demonstrate approximate conservation of the angular momentum component along this direction, a consequence of the WB orbit adiabatically adjusting to the much slower Galactic orbit. WBs with very little angular momentum in this direction are less stable over Gyr periods. This novel direction-dependent effect might allow for further tests of MOND.

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 non-linear 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 self-consistent 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.

(Abridged) Any viable cosmological model in which galaxies interact predicts the existence of primordial and tidal dwarf galaxies (TDGs). In particular, in the standard model of cosmology ($\Lambda$CDM), according to the dual dwarf galaxy theorem, there must exist both primordial dark matter-dominated and dark matter-free TDGs with different radii. We study the frequency, evolution, and properties of TDGs in a $\Lambda$CDM cosmology. We use the hydrodynamical cosmological Illustris-1 simulation to identify tidal dwarf galaxy candidates (TDGCs) and study their present-day physical properties. We also present movies on the formation of a few galaxies lacking dark matter, confirming their tidal dwarf nature. TDGCs can however also be formed via other mechanisms, such as from ram-pressure-stripped material or, speculatively, from cold-accreted gas. We find 97 TDGCs with $M_{stellar} >5 \times 10^7 M_\odot$ at redshift $z = 0$, corresponding to a co-moving number density of $2.3 \times 10^{-4} h^3 cMpc^{-3}$. The most massive TDGC has $M_{total} = 3.1 \times 10^9 M_\odot$, comparable to that of the Large Magellanic Cloud. TDGCs are phase-space-correlated, reach high metallicities, and are typically younger than dark matter-rich dwarf galaxies. We report for the first time the verification of the dual dwarf theorem in a self-consistent $\Lambda$CDM cosmological simulation. Simulated TDGCs and dark matter-dominated galaxies populate different regions in the radius-mass diagram in disagreement with observations of early-type galaxies. The dark matter-poor galaxies formed in Illustris-1 have comparable radii to observed dwarf galaxies and to TDGs formed in other galaxy-encounter simulations. In Illustris-1, only 0.17% of all selected galaxies with $M_{stellar} = 5 \times 10^7-10^9 M_\odot$ are TDGCs or dark matter-poor dwarf galaxies. The occurrence of NGC 1052-DF2-type objects is discussed.

The observed line-of-sight velocity dispersion $\sigma_{los}$ of the ultra diffuse galaxy Dragonfly 44 (DF44) requires a Newtonian dynamical mass-to-light ratio of $M_{dyn}/L_I=26^{+7}_{-6}$ Solar units. This is well outside the acceptable limits of our stellar population synthesis (SPS) models, which we construct using the integrated galactic initial mass function (IGIMF) theory. Assuming DF44 is in isolation and using Jeans analysis, we calculate $\sigma_{los}$ profiles of DF44 in Milgromian dynamics (MOND) and modified gravity (MOG) theories without invoking dark matter. Comparing with the observed kinematics, the best-fitting MOND model has $M_{dyn}/L_I = 3.6_{-1.2}^{+1.6}$ and a constant orbital anisotropy of $\beta=-0.5_{-1.6}^{+0.4}$. In MOG, we first fix its two theoretical parameters $\alpha$ and $\mu$ based on previous fits to the observed rotation curve data of The HI Nearby Galaxy Survey (THINGS). The DF44 $\sigma_{los}$ profile is best fit with $M_{dyn}/L_I=7.4_{-1.4}^{+1.5}$, larger than plausible SPS values. MOG produces a $\sigma_{los}$ profile for DF44 with acceptable $M_{dyn}/L_I$ and isotropic orbits if $\alpha$ and $\mu$ are allowed to vary. MOND with the canonical $a_0$ can explain DF44 at the $2.40 \sigma$ confidence level (1.66\%) if considering both its observed kinematics and typical star formation histories in an IGIMF context. However, MOG is ruled out at $5.49 \sigma$ ($P$-value of $4.07 \times 10^{-8}$) if its free parameters are fixed at the highest values consistent with THINGS data.

Recently van Dokkum et al. (2018b) reported that the galaxy NGC 1052-DF2 (DF2) lacks dark matter if located at $20$ Mpc from Earth. In contrast, DF2 is a dark-matter-dominated dwarf galaxy with a normal globular cluster population if it has a much shorter distance near $10$ Mpc. However, DF2 then has a high peculiar velocity wrt. the cosmic microwave background of $886$ $\rm{km\,s^{-1}}$, which differs from that of the Local Group (LG) velocity vector by $1298$ $\rm{km\,s^{-1}}$ with an angle of $117 \, ^{\circ}$. Taking into account the dynamical $M/L$ ratio, the stellar mass, half-light radius, peculiar velocity, motion relative to the LG, and the luminosities of the globular clusters, we show that the probability of finding DF2-like galaxies in the lambda cold dark matter ($\Lambda$CDM) TNG100-1 simulation is at most $1.0\times10^{-4}$ at $11.5$ Mpc and is $4.8\times10^{-7}$ at $20.0$ Mpc. At $11.5$ Mpc, the peculiar velocity is in significant tension with the TNG100-1, TNG300-1, and Millennium simulations, but occurs naturally in a Milgromian cosmology. At $20.0$ Mpc, the unusual globular cluster population would challenge any cosmological model. Estimating that precise measurements of the internal velocity dispersion, stellar mass, and distance exist for $100$ galaxies, DF2 is in $2.6\sigma$ ($11.5$ Mpc) and $4.1\sigma$ ($20.0$ Mpc) tension with standard cosmology. Adopting the former distance for DF2 and assuming that NGC 1052-DF4 is at $20.0$ Mpc, the existence of both is in tension at $\geq4.8\sigma$ with the $\Lambda$CDM model. If both galaxies are at $20.0$ Mpc the $\Lambda$CDM cosmology has to be rejected by $\geq5.8\sigma$.

The Galaxy is conventionally thought to be surrounded by a massive dark matter (DM) halo. As the Sun goes through this halo, it excites a DM wake behind it. This local asymmetry in the DM distribution would gravitationally affect the motions of Solar System planets, potentially allowing the DM wake to be detected or ruled out. Hernandez (2019) recently calculated that the DM-induced perturbation to Saturn's position is 252 metres net of the effect on the Sun. No such anomaly is seen in Saturn's motion despite very accurate tracking of the Cassini spacecraft, which orbited Saturn for >13 years. Here, we revisit the calculation of how much Saturn would deviate from Keplerian motion if we fix its position and velocity at some particular time. The DM wake induces a nearly resonant perturbation whose amplitude grows almost linearly with time. We show that the Hernandez (2019) result applies only for an observing duration comparable to the ${\approx 250}$ million year period of the Sun's orbit around the Galaxy. Over a 100 year period, the perturbation to Saturn's orbit amounts to <1 cm, which is quite consistent with existing observations. Even smaller perturbations are expected for the terrestrial planets.

The relative velocity distribution of wide binary (WB) stars is sensitive to the law of gravity at the low accelerations typical of galactic outskirts. I consider the feasibility of this wide binary test using the `line velocity' method. This involves considering only the velocity components along the direction within the sky plane orthogonal to the systemic proper motion of each WB. I apply this technique to the WB sample of Hernandez et. al. (2019), carefully accounting for large-angle effects at one order beyond leading. Based on Monte Carlo trials, the uncertainty in the one-dimensional velocity dispersion is $\approx 100$ m/s when using sky-projected relative velocities. Using line velocities reduces this to $\approx 30$ m/s because these are much less affected by distance uncertainties. My analysis does not support the Hernandez et. al. (2019) claim of a clear departure from Newtonian dynamics beyond a radius of $\approx 10$ kAU, partly because I use $2\sigma$ outlier rejection to clean their sample first. Nonetheless, the uncertainties are small enough that existing WB data are nearly sufficient to distinguish Newtonian dynamics from Modified Newtonian Dynamics. I estimate that $\approx 1000$ WB systems will be required for this purpose if using only line velocities. In addition to a larger sample, it will also be important to control for systematics like undetected companions and moving groups. This could be done statistically. The contamination can be minimized by considering a narrow theoretically motivated range of parameters and focusing on how different theories predict different proportions of WBs in this region.

We consider how the trajectory of an interstellar precursor mission would be affected by the gravity of the Sun in Newtonian and Milgromian dynamics (MOND). The Solar gravity is ${\approx 50\%}$ stronger in MOND beyond a distance of ${\approx 7000}$ astronomical units, the MOND radius of the Sun. A spacecraft travelling at $0.01 \, c$ reaches this distance after 11.1 years. We show that the extra gravity in MOND causes an anomalous deceleration that reduces its radial velocity by ${\approx 3}$ cm/s and the two-way light travel time from the inner Solar System by ${\approx 0.1}$ seconds after 20 years. A distinctive signature of MOND is that the gravity from the Sun is not directly towards it. This is due to the non-linear nature of MOND and the external gravitational field from the rest of the Galaxy, which we self-consistently include in our calculations. As a result, the sky position of the spacecraft would deviate by up to 0.2 mas over 20 years. This deviation is always in the plane containing the spacecraft trajectory and the direction towards the Galactic Centre. By launching spacecraft in different directions, it is possible to test the characteristic pattern of angular deviations expected in MOND. This would minimize the chance that any detected anomalies are caused by other processes like drag from the interstellar medium. Such confounding factors could also be mitigated using an onboard accelerometer to measure non-gravitational forces. We briefly discuss how the gravity theories could be conclusively distinguished using a Cavendish-style active gravitational experiment beyond the Sun's MOND radius.

The wide binary orbit of Proxima Centauri around $\alpha$ Centauri A and B differs significantly between Newtonian and Milgromian dynamics (MOND). By combining previous calculations of this effect with mock observations generated using a Monte Carlo procedure, we show that this prediction can be tested using high precision astrometry of Proxima Centauri. This requires ${\approx 10}$ years of observations at an individual epoch precision of $0.5 \, \mu$as, within the design specifications of the proposed Theia mission. In general, the required duration should scale as the 2/5 power of the astrometric precision. A long-period planet could produce a MOND-like astrometric signal, but only if it has a particular ratio of mass to separation squared and a sky position close to the line segment connecting Proxima Centauri with $\alpha$ Centauri. Uncertainties in perspective effects should be small enough for this test if the absolute radial velocity of Proxima Centauri can be measured to within ${\approx 10}$ m/s, better than the present accuracy of 32 m/s. We expect the required improvement to become feasible using radial velocity zero points estimated from larger samples of close binaries, with the Sun providing an anchor. We demonstrate that possible astrometric microlensing of Proxima Centauri is unlikely to affect the results. We also discuss why it should be possible to find sufficiently astrometrically stable reference stars. Adequately addressing these and other issues would enable a decisive test of gravity in the currently little explored low acceleration regime relevant to the dynamical discrepancies in galactic outskirts.

The ultra-diffuse dwarf galaxy NGC 1052-DF2 (DF2) has ten (eleven) measured globular clusters (GCs) with a line-of-sight velocity dispersion of $\sigma=7.8^{+5.2}_{-2.2}\,$km/s ($\sigma=10.6^{+3.9}_{-2.3}\,$km/s). Our conventional statistical analysis of the original ten GCs gives $\sigma=8.0^{+4.3}_{-3.0}\,$km/s. The overall distribution of velocities agrees well with a Gaussian of this width. Due to the non-linear Poisson equation in MOND, a dwarf galaxy has weaker self-gravity when in close proximity to a massive host. This external field effect is investigated using a new analytic formulation and fully self-consistent live $N$-body models in MOND. Our formulation agrees well with that of Famaey and McGaugh (2012). These new simulations confirm our analytic results and suggest that DF2 may be in a deep-freeze state unique to MOND. The correctly calculated MOND velocity dispersion agrees with our inferred dispersion and that of van Dokkum et al. (2018b) if DF2 is within 150 kpc of NGC 1052 and both are 20 Mpc away. The GCs of DF2 are however significantly brighter and larger than normal GCs, a problem which disappears if DF2 is significantly closer to us. A distance of 10-13 Mpc makes DF2 a normal dwarf galaxy even more consistent with MOND and the 13 Mpc distance reported by Trujillo et. al. (2019). We discuss the similar dwarf DF4, finding good agreement with MOND. We also discuss possible massive galaxies near DF2 and DF4 along with their distances and peculiar velocities, noting that NGC 1052 may lie at a distance near 10 Mpc.

A great challenge in present-day physics is to understand whether the observed internal dynamics of galaxies is due to dark matter matter or due to a modification of the law of gravity. Recently, van Dokkum et al. reported that the ultra-diffuse dwarf galaxy NGC1052-DF2 lacks dark matter, and they claimed that this would -- paradoxically -- be problematic for modified gravity theories like Milgromian dynamics (MOND). However, NGC1052-DF2 is not isolated, so that a valid prediction of its internal dynamics in MOND cannot be made without properly accounting for the external gravitational fields from neighbouring galaxies. Including this external field effect following Haghi et al. shows that NGC1052-DF2 is consistent with MOND.

We consider disk stability in the quasi-linear formulation of MOND (QUMOND), the basis for some $N$-body integrators. We generalize the Toomre criterion for the stability of disks to tightly wound, axisymmetric perturbations. We apply this to a family of thin exponential disks with different central surface densities. By numerically calculating their QUMOND rotation curves, we obtain the minimum radial velocity dispersion required for stability against local self-gravitating collapse. MOND correctly predicts much higher rotation speeds in low surface brightness galaxies (LSBs) than does Newtonian dynamics without dark matter. Newtonian models thus require putative very massive halos, whose inert nature implies they would strongly stabilize the disk. MOND also increases the stability of galactic disks, but in contradistinction to Newtonian gravity, this extra stability is limited to a factor of 2. MOND is thus rather more conducive to the formation of bars and spiral arms. Therefore, observation of such features in LSBs could be problematic for Newtonian galaxy models. This could constitute a crucial discriminating test. We quantitatively account for these facts in QUMOND. We also compare numerical QUMOND rotation curves of thin exponential disks to those predicted by two algebraic expressions commonly used to calculate MOND rotation curves. For the choice that best approximates QUMOND, we find the circular velocities agree to within 1.5% beyond $\approx 0.5$ disk scale lengths, regardless of the central surface density. The other expression can underestimate the rotational speed by up to 12.5% at one scale length, though rather less so at larger radii.

A central problem of contemporary physics is whether the law of gravity is non-Newtonian on galaxy scales. Rodrigues et al. argue that Milgromian gravitation, which solves the flat rotation curve problem without the need for dark matter particles, is ruled out at > 10{\sigma} significance. To a large extent, this conclusion relies on galaxies with very uncertain distances and/or nearly edge-on orientations, where dust obscuration often becomes significant. Applying appropriate quality cuts to the data leaves only a handful of outliers to the predictions of Milgromian gravitation according to the analysis of Rodrigues et al., but even these outliers can be explained with Milgromian gravitation.

The works in this portfolio test the hypothesis that it is not possible to extrapolate the Newtonian inverse square law of gravity from Solar System to galaxy scales. In particular, I look into various tests of Modified Newtonian Dynamics (MOND), which posits a modification below a very low acceleration threshold. Although discrepancies with Newtonian dynamics are indeed observed, they can usually be explained by invoking an appropriate distribution of invisible mass known as dark matter (DM). This leads to the standard cosmological paradigm, $\Lambda$CDM. I consider how it may be distinguished from MOND using collision velocities of galaxy clusters, which should sometimes be much faster in MOND. I focus on measuring these velocities more accurately and conclude that this test ought to be feasible in the near future. For the time being, I look at the much nearer and more accurately observed Local Group (LG) of galaxies. Its main constituents $-$ the Milky Way (MW) and Andromeda (M31) $-$ should have undergone a past close flyby in MOND but not in $\Lambda$CDM. The fast MW-M31 relative motion around the time of their flyby would have allowed them to gravitationally slingshot any passing LG dwarf galaxies out at high speed. I consider whether there is any evidence for such high-velocity galaxies (HVGs). Several candidates are found in two different $\Lambda$CDM models of the LG, one written by a founding figure of the paradigm. The properties of these HVGs are similar to what might be expected in MOND, especially their tendency to lie close to a plane. Being more confident of its validity, I then used MOND to determine the escape velocity curve of the MW over the distance range 8$-$50 kpc, finding reasonable agreement with the latest observations. I finish by discussing possible future directions for MOND research.

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 MW-M31 flyby occurred ${\approx 8}$ Gyr ago. Using the timing argument, we obtain MW-M31 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 best-fitting model yields a significant fraction (0.22) of counter-rotating material, perhaps explaining why Sculptor counter-rotates within the MW satellite plane. In contrast, our model yields no counter-rotating 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 MW-M31 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 MW-M31 barycentre and is oriented so the MW-M31 line is inclined by $16^\circ$ to it. We develop a toy model to constrain the scenario whereby a past MW-M31 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 MW-M31 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 fast-moving MW and M31 gravitationally slingshot test particles outwards at high speeds. The most distant such particles preferentially lie within the MW-M31 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}$).

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 best-fitting 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 fast-moving 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.

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 time-dependent (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 edge-on 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 best-fitting $\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 fast-moving galaxies could plausibly explain why some non-satellite LG galaxies are moving away from us even faster than a pure Hubble flow.

In this short note we investigate the ultra-relativistic collisions of small amplitude oscillons in 1+1 dimensions. Using the amplitude of the oscillons and the inverse relativistic boost factor $\gamma^{-1}$ as the perturbation variables, we analytically calculate the leading order spatial and temporal phase shifts, and the change in the amplitude of the oscillons after the collisions. At leading order, we find that only the temporal phase shift receives a nonzero contribution, and that the collision is elastic. This work is also the first application of the general kinematic framework for understanding ultra-relativistic collisions (arXiv:1308.0606) to intrinsically time-dependent solitons.

I present an analysis of ~13000 stars from RAVE Data Release 4 (DR4) to better understand the formation of the thick disk. The stars I consider are mostly within 1 kpc of the Sun. Based on a Monte-Carlo analysis of the data and comparison with a toy model, I suggest the thick disk formed fairly rapidly when [Mg/Fe] ${\approx}$0.2-0.3. We consider an origin via direct accretion unlikely as the thick disk is fast-rotating and fairly metal-rich for its age. As internal disk heating can not easily reach the high observed velocity dispersions and would likely be more gradual, it also appears implausible. A very early formation of the thick disk (e.g. from star formation in the collapsing primordial gas cloud of the Milky Way) appears unlikely given it would require [Mg/Fe] ${\geq}$ 0.4 at the time the thick disk formed. An origin via tidal interaction with another galaxy is one of the few remaining possibilities. My proposed explanation for the nature of the event which formed the thick disk is a close flyby of Andromeda, whose orbital dynamics in MOND are consistent with my estimate that the event occurred 7-11 Gyr ago. This leads to reasonable thick disk velocity dispersions. It will be important to determine whether the event heated the outer regions of the thin disk uniformly or if substantial portions of it were left largely unaffected. This may distinguish between models (like mine) involving a single close encounter and models with multiple encounters (as common in mergers). The main limiting factor at present seems to be the accuracy of elemental abundances and lack of direct observations at different galactocentric radii to the Sun.