Non-planetary bodies provide valuable insight into our current under- standing of planetary formation and evolution. Although these objects are challeng- ing to detect and characterize, the potential information to be drawn from them has motivated various searches through a number of techniques. Here, we briefly review the current status in the search of moons, rings, comets, and trojans in exoplanet systems and suggest what future discoveries may occur in the near future.

The inability of primordial baryonic density fluctuations, as observed in the cosmic microwave background (CMB), to grow into the present day astronomical structures is well established, under Newtonian and Einsteinian gravity. It is hence customary to assume the existence of an underlying dark matter component with density fluctuations, $\Delta(M)$, having amplitudes much larger than what CMB observations imply for the baryons. This is in fact one of the recurrent arguments used in support of the dark matter hypothesis. In this letter we prove that the same extended theory of gravity which has been recently shown to accurately reproduce gravitational lensing observations, in absence of any dark matter, and which in the low velocity regime converges to a MONDian force law, implies a sufficiently amplified self-gravity to allow purely baryonic fluctuations with amplitudes in accordance with CMB constraints to naturally grow into the $z=0$ astrophysical structures detected. The linear structure formation scenario which emerges closely resembles the standard concordance cosmology one, as abundantly calibrated over the last decade to match multiple observational constraints at various redshifts. However, in contrast with what occurs in the concordance cosmology, this follows not from a critical dependence on initial conditions and the fine tuning of model parameters, but from the rapid convergence of highly arbitrary initial conditions onto a well defined $\Delta(M,z)$ attractor solution.

That the stellar halo of the Milky Way has a density profile which to first approximation satisfies $\rho \propto r^{-3}$ has been known for a long time. More recently, it has become clear that M31 also has such an extended stellar halo, which approximately follows the same radial scaling. Studies of distant galaxies have revealed the same phenomenology. Also, we now know that the density profiles of the globular cluster systems of our Galaxy and Andromeda to first approximation follow $\rho \propto r^{-3}$, $\Sigma \propto R^{-2}$ in projection. Recently, diffuse populations of stars have been detected spherically surrounding a number of Galactic globular clusters, extending much beyond the Newtonian tidal radii, often without showing any evidence of tidal features. Within the standard Newtonian and GR scenario, numerous and diverse particular explanations have been suggested, individually tailored to each of the different classes of systems described above. Here we show that in a MONDian gravity scenario, any isothermal tenuous halo of tracer particles forming a small perturbation surrounding a spherically symmetric mass distribution, will have an equilibrium configuration which to first approximation satisfies a $\rho \propto r^{-3}$ scaling.

Elliptical galaxies are systems where dark matter is usually less necessary to explain observed dynamics than in the case of spiral galaxies, however there are some instances where Newtonian gravity and the observable mass are insufficient to explain their observed structure and kinematics. Such is the case of NGC 4649, a massive elliptical galaxy in the Virgo cluster for which recent studies report a high fraction of dark matter, 0.78 at $4R_e$. However this galaxy has been studied within the MOND hypothesis, where a good agreement with the observed values of velocity dispersion is found. Using a MONDian gravity force law, here we model this galaxy as a self-consistent gravitational equilibrium dynamical system. This force law reproduces the MOND phenomenology in the $a<a_{0}$ regime, and reduces to the Newtonian case when $a>a_{0}$. Within the MONDian $a<a_{0}$ scales, centrifugal equilibrium or dispersion velocities become independent of radius, and show a direct proportionality to the fourth root of the total baryonic mass, $V^{4}\propto(M G a_{0})$. We find that the recent detailed observations of the surface brightness profile and the velocity dispersion profile for this galaxy are consistent with the phenomenology expected in MONDian theories of modified gravity, without the need of invoking the presence of any hypothetical dark matter.

Over the past couple of years, a number of observational studies have confirmed the flattening of the radial velocity dispersion profiles for stars in various nearby globular clusters. As the projected radial coordinate is increased, a radius appears beyond which, the measured velocity dispersion ceases to drop and settles at a fixed value, $\sigma_{\infty}$. Under Newtonian gravity, this is explained by invoking tidal heating from the overall Milky Way potential on the outer, more loosely bound stars, of the globular clusters in question. From the point of view of modified gravity theories, such an outer flattening is expected on crossing the critical acceleration threshold $a_{0}$, beyond which, a transition to MONDian dynamics is expected, were equilibrium velocities cease to be a function of distance. In this paper we attempt to sort out between the above competing explanations, by looking at their plausibility in terms of an strictly empirical approach. We determine Newtonian tidal radii using masses accurately calculated through stellar population modelling, and hence independent of any dynamical assumptions, distances, size and orbital determinations for a sample of 16 globular clusters. We show that their Newtonian tidal radii at perigalacticon are generally larger that the radii at which the flattening in the velocity dispersion profiles occurs, by large factors of 4, on average. While this point makes the Newtonian tidal explanation suspect, it is found that the radii at which the flattening is observed on average correlate with the radii where the $a_{0}$ threshold is crossed, and that $\sigma_{\infty}$ values scale with the fourth root of the total masses, all features predicted under modified gravity theories.

Assuming Newton's gravity and GR to be valid at all scales leads to the dark matter hypothesis as a requirement demanded by the observed dynamics and measured baryonic content at galactic and extragalactic scales. Alternatively, modified gravity scenarios where a change of regime appears at acceleration scales $a<a_{0}$ have been proposed. This modified regime at $a<a_{0}$ will generically be characterised by equilibrium velocities which become independent of distance. Here we identify a critical test in this debate and we propose its application to samples of wide binary stars. Since for $1 M_{\odot}$ systems the acceleration drops below $a_{0}$ at scales of around 7000 AU, a statistical survey of wide binaries with relative velocities and separations reaching $10^{4}$ AU and beyond should prove useful to the above debate. We apply the proposed test to the best currently available data. Results show a constant upper limit to the relative velocities in wide binaries which is independent of separation for over three orders of magnitude, in analogy with galactic flat rotation curves in the same $a<a_{0}$ acceleration regime. Our results are suggestive of a breakdown of Kepler's third law beyond $a \approx a_{0}$ scales, in accordance with generic predictions of modified gravity theories designed not to require any dark matter at galactic scales and beyond.

In the context of theories of gravity modified to account for the observed dynamics of galactic systems without the need to invoke the existence of dark matter, a prediction often appears regarding low acceleration systems: wherever $a$ falls below $a_{0}$ one should expect a transition from the classical to the modified gravity regime.This modified gravity regime will be characterised by equilibrium velocities which become independent of distance, and which scale with the fourth root of the total baryonic mass, $V^{4} \propto M$. The two above conditions are the well known flat rotation curves and Tully-Fisher relations of the galactic regime. Recently however, a similar phenomenology has been hinted at, at the outskirts of Galactic globular clusters, precisely in the region where $a<a_{0}$. Radial profiles of the projected velocity dispersion have been observed to stop decreasing along Keplerian expectations, and to level off at constant values beyond the radii where $a<a_{0}$. We have constructed gravitational equilibrium dynamical models for a number of globular clusters for which the above gravitational anomaly has been reported, using a modified Newtonian force law which yields equilibrium velocities equivalent to MOND. We find models can be easily constructed having an inner Newtonian region and an outer modified gravity regime, which reproduce all observational constraints, surface brightness profiles, total masses and line of sight velocity dispersion profiles. Through the use of detailed single stellar population models tuned individually to each of the globular clusters in question, we derive estimates of the total masses for these systems. Interestingly, we find that the asymptotic values of the velocity dispersion profiles are consistent with scaling with the fourth root of the total masses, as expected under modified gravity scenarios.

Using dimensional analysis techniques we present an extension of Newton's gravitational theory built under the assumption that Milgrom's acceleration constant is a fundamental quantity of nature. The gravitational force converges to Newton's gravity and to a MOND-like description in two different mass and length regimes. It is shown that a modification on the force sector (and not in the dynamical one as MOND does) is more convenient and can reproduce and predict different phenomena usually ascribed to dark matter at the non-relativistic level.

Modified gravity scenarios where a change of regime appears at acceleration scales $a<a_{0}$ have been proposed. Since for $1 M_{\odot}$ systems the acceleration drops below $a_{0}$ at scales of around 7000 AU, a statistical survey of wide binaries with relative velocities and separations reaching $10^{4}$ AU and beyond should prove useful to the above debate. We apply the proposed test to the best currently available data. Results show a constant upper limit to the relative velocities in wide binaries which is independent of separation for over three orders of magnitude, in analogy with galactic flat rotation curves in the same $a<a_{0}$ acceleration regime. Our results are suggestive of a breakdown of Kepler's third law beyond $a \approx a_{0}$ scales, in accordance with generic predictions of modified gravity theories designed not to require any dark matter at galactic scales and beyond.