We discuss the in-flight autonomous assembly as the means to build advanced planetary science payloads to explore the outer regions of the solar system. These payloads are robotically constructed from modular parts delivered by a group of smallsats (< 20 kg) which are placed on fast solar system transfer trajectories while being accelerated by solar sail propulsion to velocities of ~10 AU/yr. This concept provides the planetary science community with inexpensive, frequent access to distant regions of the solar system with flexible, reconfigurable instruments and systems that are assembled in flight. It permits faster revisit times, rapid replenishment and technology insertions, longer mission capability with lower costs. It also increases the science capabilities of smallsats via the use of modular, redundant architectures and allows for proliferation of sensing instrumentation throughout the solar system.

We examined the solar gravitational lens (SGL) as the means to produce direct high-resolution, multipixel images of exoplanets. The properties of the SGL are remarkable: it offers maximum light amplification of ~1e11 and angular resolution of ~1e-10 arcsec. A probe with a 1-m telescope in the SGL focal region can image an exoplanet at 30 pc with 10-kilometer resolution on its surface, sufficient to observe seasonal changes, oceans, continents, surface topography. We reached and exceeded all objectives set for our study: We developed a new wave-optical approach to study the imaging of exoplanets while treating them as extended, resolved, faint sources at large but finite distances. We properly accounted for the solar corona brightness. We developed deconvolution algorithms and demonstrated the feasibility of high-quality image reconstruction under realistic conditions. We have proven that multipixel imaging and spectroscopy of exoplanets with the SGL are feasible. We have developed a new mission concept that delivers an array of optical telescopes to the SGL focal region relying on three innovations: i) a new way to enable direct exoplanet imaging, ii) use of smallsats solar sails fast transit through the solar system and beyond, iii) an open architecture to take advantage of swarm technology. This approach enables entirely new missions, providing a great leap in capabilities for NASA and the greater aerospace community. Our results are encouraging as they lead to a realistic design for a mission that will be able to make direct resolved images of exoplanets in our stellar neighborhood. It could allow exploration of exoplanets relying on the SGL capabilities decades, if not centuries, earlier than possible with other extant technologies. The architecture and mission concepts for a mission to the strong interference region of the SGL are promising and should be explored further.

We study the image formation process with the solar gravitational lens (SGL) in the case of an extended, resolved source. An imaging telescope, modeled as a convex lens, is positioned within the image cylinder formed by the light received from the source. In the strong interference region of the SGL, this light is greatly amplified, forming the Einstein ring around the Sun, representing a distorted image of the extended source. We study the intensity distribution within the Einstein ring observed in the focal plane of the convex lens. For any particular telescope position in the image plane, we model light received from the resolved source as a combination of two signals: light received from the directly imaged region of the source and light from the rest of the source. We also consider the case when the telescope points away from the extended source or, equivalently, it observes light from sources in sky positions that are some distance away from the extended source, but still in its proximity. At even larger distances from the optical axis, in the weak interference or geometric optics regions, our approach recovers known models related to microlensing, but now obtained via the wave-optical treatment. We then derive the power of the signal and related photon fluxes within the annulus that contains the Einstein ring of the extended source, as seen by the imaging telescope. We discuss the properties of the deconvolution process, especially its effects on noise in the recovered image. We compare anticipated signals from realistic exoplanetary targets against estimates of noise from the solar corona and estimate integration times needed for the recovery of high-quality images of faint sources. The results demonstrate that the SGL offers a unique, realistic capability to obtain resolved images of exoplanets in our galactic neighborhood.

We study image formation with the solar gravitational lens (SGL). We consider a point source that is positioned at a large but finite distance from the Sun. We assume that an optical telescope is positioned in the image plane, in the focal region of the SGL. We model the telescope as a convex lens and evaluate the intensity distribution produced by the electromagnetic field that forms the image in the focal plane of the convex lens. We first investigate the case when the telescope is located on the optical axis of the SGL or in its immediate vicinity. This is the region of strong interference where the SGL forms an image of a distant source, which is our primary interest. We derive analytic expressions that describe the progression of the image from an Einstein ring corresponding to an on-axis telescope position, to the case of two bright spots when the telescope is positioned some distance away from the optical axis. At greater distances from the optical axis, in the region of weak interference and that of geometric optics, we recover expressions that are familiar from models of gravitational microlensing, but developed here using a wave-optical treatment. We discuss applications of the results for imaging and spectroscopy of exoplanets with the SGL.

We discuss the optical properties of the solar gravitational lens (SGL). We estimate the power of the EM field received by an imaging telescope. Studying the behavior of the EM field at the photometric detector, we develop expressions that describe the received power from a point source as well as from an extended resolved source. We model the source as a disk with uniform surface brightness and study the contribution of blur to a particular image pixel. To describe this process, we develop expressions describing the power received from the directly imaged region of the exoplanet, from the rest of the exoplanet and also the power for off-image pointing. We study the SGL's amplification and its angular resolution in the case of observing an extended source with a modest size telescope. The results can be applied to direct imaging of exoplanets using the SGL.

We present calculations for the anticipated shadow sizes of Sgr A* and the supermassive black hole in the galaxy M87 in the context of the MOG modified theory of gravity (also known as Scalar-Tensor-Vector-Gravity, or STVG). We demonstrate that mass estimates derived from stellar and gas dynamics in the vicinity of these black holes are the Newtonian masses of the black holes even in the MOG theory. Consequently, shadow sizes increase as a function of the key dimensionless MOG parameter alpha that characterizes the variable gravitational coupling coefficient G, and may offer an observational means to distinguish the MOG theory from standard general relativity.

We present a statistical study on the orientation of outflows with respect to large-scale filaments and the magnetic fields. Although filaments are widely observed toward Galactic star-forming regions, the exact role of filaments in star formation is unclear. Studies toward low-mass star-forming regions revealed both preferred and random orientation of outflows respective to the filament long-axes, while outflows in massive star-forming regions mostly oriented perpendicular to the host filaments, and parallel to the magnetic fields at similar physical scales. Here, we explore outflows in a sample of 11 protoclusters in HII regions, a more evolved stage compared to IRDCs, using ALMA CO (3-2) line observations. We identify a total of 105 outflow lobes in these protoclusters. Among the 11 targets, 7 are embedded within parsec-scale filamentary structures detected in $^{13}$CO line and 870 $\mu m$ continuum emissions. The angles between outflow axes and corresponding filaments ($\gamma_\mathrm{Fil}$) do not show any hint of preferred orientations (i.e., orthogonal or parallel as inferred in numerical models) with respect to the position angle of the filaments. Identified outflow lobes are also not correlated with the magnetic fields and Galactic plane position angles. Outflows associated with filaments aligned along the large-scale magnetic fields are also randomly orientated. Our study presents the first statistical results of outflow orientation respective to large-scale filaments and magnetic fields in evolved massive star-forming regions. The random distribution suggests a lack of alignment of outflows with filaments, which may be a result of the evolutionary stage of the clusters.

We investigate the optical properties of the solar gravitational lens (SGL) with respect to an extended source located at a large but finite distance from the Sun. The static, spherically symmetric gravitational field of the Sun is modeled within the first post-Newtonian approximation of the general theory of relativity. We consider the propagation of monochromatic electromagnetic (EM) waves near the Sun. We develop, based on a Mie theory, a vector theory of diffraction that accounts for the refractive properties of the solar gravitational field. The finite distance to a point source can be accounted for using a rotation of the coordinate system to align its polar axis with the axis directed from the point source to the center of the Sun, which we call the optical axis. We determine the EM field and study the key optical properties of the SGL in all four regions formed behind the Sun by an EM wave diffracted by the solar gravity field: the shadow, geometric optics, and weak and strong interference regions. Extended sources can then be represented as collections of point sources. We present the power density of the signal received by a telescope in the image plane. Our discussion concludes with considering the implications for imaging with the SGL.

We study the propagation of electromagnetic (EM) waves in the solar system and develop a Mie theory that accounts for the refractive properties of the free electron plasma in the extended solar corona. We use a generic model for the electron number density distribution and apply the eikonal approximation to find a solution in terms of Debye potentials, which is then used to determine the EM field both within the inner solar system and at large heliocentric distances. As expected, the solution for the EM wave propagating through the solar system is characterized by a plasma-induced phase shift and related change in the light ray's direction of propagation. Our approach quantitatively accounts for these effects, providing a wave-optical treatment for diffraction in the solar plasma. As such, it may be used in practical applications involving big apertures, large interferometric baselines or otherwise widely distributed high-precision astronomical instruments.

We study the optical properties of the solar gravitational lens (SGL) under the combined influence of the static spherically symmetric gravitational field of the Sun---modeled within the first post-Newtonian approximation of the general theory of relativity---and of the solar corona---modeled as a generic, steady-state, spherically symmetric free electron plasma. For this, we consider the propagation of monochromatic electromagnetic (EM) waves near the Sun and develop a Mie theory that accounts for the refractive properties of the gravitational field of the Sun and that of the free electron plasma in the extended solar system. We establish a compact, closed-form solution to the boundary value problem, which extends previously known results into the new regime where gravity and plasma are both present. Relying on the wave-optical approach, we consider three different regions of practical importance for the SGL, including the shadow region directly behind the Sun, the region of geometrical optics and the interference region. We demonstrate that the presence of the solar plasma affects all characteristics of an incident unpolarized light, including the direction of the EM wave propagation, its amplitude and its phase. We show that the presence of the solar plasma leads to a reduction of the light amplification of the SGL and to a broadening of its point spread function. We also show that the wavelength-dependent plasma effect is important at radio frequencies, where it drastically reduces both the amplification factor of the SGL and also its angular resolution. However, for optical and shorter wavelengths, the plasma's contribution to the EM wave is negligibly small, leaving the plasma-free optical properties of the SGL practically unaffected.

We model the velocity dispersion of the ultra-diffuse galaxy NGC 1052-DF2 using Newtonian gravity and modified gravity (MOG). The velocity dispersion predicted by MOG is higher than the Newtonian gravity prediction, but it is fully consistent with the observed velocity dispersion that is obtained from the motion of 10 globular clusters.

We investigate gravitational lensing in the context of the MOG modified theory of gravity. Using a formulation of the theory with no adjustable or fitted parameters, we present the MOG equations of motion for slow, nonrelativistic test particles and for ultrarelativistic test particles, such as rays of light. We demonstrate how the MOG prediction for the bending of light can be applied to astronomical observations. Our investigation first focuses on a small set of strong lensing observations where the properties of the lensing objects are found to be consistent with the predictions of the theory. We also present an analysis of the colliding clusters 1E0657-558 (known also as the Bullet Cluster) and Abell 520; in both cases, the predictions of the MOG theory are in good agreement with observation.

The Solar Gravitational Lens (SGL) allows for major brightness amplification ($\sim 10^{11}$ at wavelength of $1~\mu$m) and extreme angular resolution ($\sim10^{-10}$ arcsec) within a narrow field of view. A meter-class telescope, with a modest coronagraph to block solar light with 1e-6 suppression placed in the focal area of the SGL, can image an exoplanet at a distance of 30 parsec with few kilometer-scale resolution on its surface. Notably, spectroscopic broadband SNR is $\sim 10^{-6}$ in two weeks of integration time, providing this instrument with incredible remote sensing capabilities. A mission capable of exploiting the remarkable optical properties of the SGL allows for direct high-resolution imaging/spectroscopy of a potentially habitable exoplanet. Such missions could allow exploration of exoplanets relying on the SGL capabilities decades, if not centuries, earlier than possible with other extant technologies.

The remarkable optical properties of the solar gravitational lens (SGL) include major brightness amplification (~1e11 at wavelength of 1 um) and extreme angular resolution (~1e-10 arcsec) in a narrow field of view. A mission to the SGL carrying a modest telescope and coronagraph opens up a possibility for direct megapixel imaging and high-resolution spectroscopy of a habitable Earth-like exoplanet at a distance of up to 100 light years. The entire image of such a planet is compressed by the SGL into a region with a diameter of ~1.3 km in the vicinity of the focal line. The telescope, acting as a single pixel detector while traversing this region, can build an image of the exoplanet with kilometer-scale resolution of its surface, enough to see its surface features and signs of habitability. We report here on the results of our initial study of a mission to the deep outer regions of our solar system, with the primary mission objective of conducting direct megapixel high-resolution imaging and spectroscopy of a potentially habitable exoplanet by exploiting the remarkable optical properties of the SGL. Our main goal was to investigate what it takes to operate spacecraft at such enormous distances with the needed precision. Specifically, we studied i) how a space mission to the focal region of the SGL may be used to obtain high-resolution direct imaging and spectroscopy of an exoplanet by detecting, tracking, and studying the Einstein ring around the Sun, and ii) how such information could be used to detect signs of life on another planet. Our results indicate that a mission to the SGL with an objective of direct imaging and spectroscopy of a distant exoplanet is challenging, but possible. We composed a list of recommendations on the mission architectures with risk and return tradeoffs and discuss an enabling technology development program.

Modified gravity (MOG) is a covariant, relativistic, alternative gravitational theory whose field equations are derived from an action that supplements the spacetime metric tensor with vector and scalar fields. Both gravitational (spin 2) and electromagnetic waves travel on null geodesics of the theory's one metric. Despite a recent claim to the contrary, MOG satisfies the weak equivalence principle and is consistent with observations of the neutron star merger and gamma ray burster event GW170817/GRB170817A.

Self-gravitating Bose-Einstein condensates (BEC) have been proposed in various astrophysical contexts, including Bose-stars and BEC dark matter halos. These systems are described by a combination of the Gross-Pitaevskii and Poisson equations (the GPP system). In the analysis of these hypothetical objects, the Thomas-Fermi (TF) approximation is widely used. This approximation is based on the assumption that in the presence of a large number of particles, the kinetic term in the Gross-Pitaevskii energy functional can be neglected, yet it is well known that this assumption is violated near the condensate surface. We also show that the total energy of the self-gravitating condensate in the TF-approximation is positive. The stability of a self-gravitating system is dependent on the total energy being negative. Therefore, the TF-approximation is ill suited to formulate initial conditions in numerical simulations. As an alternative, we offer an approximate solution of the full GPP system.

TianQin is a proposal for a space-borne detector of gravitational waves in the millihertz frequencies. The experiment relies on a constellation of three drag-free spacecraft orbiting the Earth. Inter-spacecraft laser interferometry is used to monitor the distances between the test masses. The experiment is designed to be capable of detecting a signal with high confidence from a single source of gravitational waves within a few months of observing time. We describe the preliminary mission concept for TianQin, including the candidate source and experimental designs. We present estimates for the major constituents of the experiment's error budget and discuss the project's overall feasibility. Given the current level of technology readiness, we expect TianQin to be flown in the second half of the next decade.

Using recently developed nonrelativistic numerical simulation code, we investigate the stability properties of compact astrophysical objects that may be formed due to the Bose-Einstein condensation of dark matter. Once the temperature of a boson gas is less than the critical temperature, a Bose-Einstein condensation process can always take place during the cosmic history of the universe. Due to dark matter accretion, a Bose-Einstein condensed core can also be formed inside massive astrophysical objects such as neutron stars or white dwarfs, for example. Numerically solving the Gross-Pitaevskii-Poisson system of coupled differential equations, we demonstrate, with longer simulation runs, that within the computational limits of the simulation the objects we investigate are stable. Physical properties of a self-gravitating Bose-Einstein condensate are examined both in non-rotating and rotating cases.

Galaxy rotation curves determined observationally out to a radius well beyond the galaxy cores can provide a critical test of modified gravity models without dark matter. The predicted rotational velocity curve obtained from Scalar-Vector-Tensor Gravity (STVG or MOG) is in excellent agreement with data for the Milky Way without a dark matter halo, with a mass of $5\times 10^{10}\,M_{\odot}$. The velocity rotation curve predicted by modified Newtonian Dynamics (MOND) does not agree with the data.

We develop a high-precision model for laser ranging interferometric (LRI) observables of the GRACE Follow-On (GRACE-FO) mission. For this, we study the propagation of an electromagnetic wave in the gravitational field in the vicinity of an extended body, in the post-Newtonian approximation of the general theory of relativity. We present a general relativistic model for the phase of a plane wave that accounts for contributions of all the multipoles of the gravitating body, its angular momentum, as well as the contribution of tidal fields produced by external sources. We develop a new approach to model a coherent signal transmission in the gravitational field of the solar system that relies on a relativistic treatment of the phase. We use this approach to describe high-precision interferometric measurements on GRACE-FO and formulate the key LRI observables, namely the phase and phase rate of a coherent laser link between the two spacecraft. We develop a relativistic model for the LRI-enabled range between the two GRACE-FO spacecraft, accurate to less than 1 nm, and a high-precision model for the corresponding range rate, accurate to better than 0.1 nm/s. We also formulate high-precision relativistic models for the double one-way range (DOWR) and DOWR-enabled range rate observables originally used on GRACE and now studied for interferometric measurements on GRACE-FO. Our formulation justifies the basic assumptions behind the design of the GRACE-FO mission and highlights the importance of achieving nearly circular and nearly identical orbits for the GRACE-FO spacecraft.

We completed the development of simulation code that is designed to study the behavior of a conjectured dark matter galactic halo that is in the form of a Bose-Einstein Condensate (BEC). The BEC is described by the Gross-Pitaevskii equation, which can be solved numerically using the Crank-Nicholson method. The gravitational potential, in turn, is described by Poisson's equation, that can be solved using the relaxation method. Our code combines these two methods to study the time evolution of a self-gravitating BEC. The inefficiency of the relaxation method is balanced by the fact that in subsequent time iterations, previously computed values of the gravitational field serve as very good initial estimates. The code is robust (as evidenced by its stability on coarse grids) and efficient enough to simulate the evolution of a system over the course of 1E9 years using a finer (100x100x100) spatial grid, in less than a day of processor time on a contemporary desktop computer.

The Jordan-Brans-Dicke theory of gravitation, which promotes the gravitational constant to a dynamical scalar field, predicts a value for the Eddington-Robertson post-Newtonian parameter gamma that is significantly different from the general relativistic value of unity. This contradicts precision solar system measurements that tightly constrain gamma around 1. We consider a modification of the theory, in which the scalar field is sourced explicitly by matter. We find that this leads to a modified expression for the gamma-parameter. In particular, a specific choice of the scalar current yields gamma=1, just as in general relativity, while the weak equivalence principle is also satisfied. This result has important implications for theories that mimic Jordan-Brans-Dicke theory in the post-Newtonian limit in the solar system, including our scalar-tensor-vector modified gravity theory (MOG).

We investigate the possibility that the anomalous acceleration of the Pioneer 10 and 11 spacecraft is due to the recoil force associated with an anisotropic emission of thermal radiation off the vehicles. To this end, relying on the project and spacecraft design documentation, we constructed a comprehensive finite-element thermal model of the two spacecraft. Then, we numerically solve thermal conduction and radiation equations using the actual flight telemetry as boundary conditions. We use the results of this model to evaluate the effect of the thermal recoil force on the Pioneer 10 spacecraft at various heliocentric distances. We found that the magnitude, temporal behavior, and direction of the resulting thermal acceleration are all similar to the properties of the observed anomaly. As a novel element of our investigation, we develop a parameterized model for the thermal recoil force and estimate the coefficients of this model independently from navigational Doppler data. We find no statistically significant difference between the two estimates and conclude that once the thermal recoil force is properly accounted for, no anomalous acceleration remains.

Our modified gravity theory (MOG) is a gravitational theory without exotic dark matter, based on an action principle. MOG has been used successfully to model astrophysical phenomena such as galaxy rotation curves, galaxy cluster masses, and lensing. MOG may also be able to account for cosmological observations. We assume that the MOG point source solution can be used to describe extended distributions of matter via an appropriately modified Poisson equation. We use this result to model perturbation growth in MOG and find that it agrees well with the observed matter power spectrum at present. As the resolution of the power spectrum improves with increasing survey size, however, significant differences emerge between the predictions of MOG and the standard LCDM model, as in the absence of exotic dark matter, oscillations of the power spectrum in MOG are not suppressed. We can also use MOG to model the acoustic power spectrum of the cosmic microwave background. A suitably adapted semi-analytical model offers a first indication that MOG may pass this test, and correctly model the peak of the acoustic spectrum.

We explore the cosmological consequences of Modified Gravity (MOG), and find that it provides, using a minimal number of parameters, good fits to data, including CMB temperature anisotropy, galaxy power spectrum, and supernova luminosity-distance observations without exotic dark matter. MOG predicts a bouncing cosmology with a vacuum energy term that yields accelerating expansion and an age of ~13 billion years.

We comment on arXiv:1112.1320 and point out that baryonic oscillations of the matter power spectrum, while predicted by theories that do not incorporate collisionless cold dark matter, are strongly suppressed by the statistical window function that is used to process finite-sized galaxy samples. We assert that with present-day data sets, the slope of the matter power spectrum is a much stronger indicator of a theory's validity. We also argue that MOND should not be used as a strawman theory as it is not in general representative of modified gravity theories; some theories, notably our scalar-vector-tensor MOdified Gravity (MOG), offer much more successful predictions of cosmological observations.

Our modified gravity theory (MOG) was used successfully in the past to explain a range of astronomical and cosmological observations, including galaxy rotation curves, the CMB acoustic peaks, and the galaxy mass power spectrum. MOG was also used successfully to explain the unusual features of the Bullet Cluster 1E0657-558 without exotic dark matter. In the present work, we derive the relativistic equations of motion in the spherically symmetric field of a point source in MOG and, in particular, we derive equations for light bending and lensing. Our results also have broader applications in the case of extended distributions of matter, and they can be used to validate the Bullet Cluster results and provide a possible explanation for the merging clusters in Abell 520.

We compare the MOG and MOND predictions of the Tully-Fisher relationship to observations using gas-rich galaxies for which the baryonic mass can be estimated accurately. We find that both theories are in good agreement with the data. Future observations of large, gas-rich galaxies may offer a means to distinguish between these two theories.

As an alternative to the LCDM concordance model, Scalar-Tensor-Vector Modified Gravity (MOG) theory reproduces key cosmological observations without postulating the presence of an exotic dark matter component. MOG is a field theory based on an action principle, with a variable gravitational constant and a repulsive vector field with variable range. MOG yields a phenomenological acceleration law that includes strong tensorial gravity partially canceled by a repulsive massive vector force. This acceleration law can be used to model the CMB acoustic spectrum and the matter power spectrum yielding good agreement with observation. A key prediction of MOG is the presence of strong baryonic oscillations, which will be detectable by future surveys. MOG is also consistent with Type Ia supernova data. We also describe on-going research of the coupling between MOG and continuous matter, consistent with the weak equivalence principle and solar system observations.

The Sloan Digital Sky Survey (SDSS) provides data on several hundred thousand galaxies. Precise location of these galaxies in the sky, along with information about their luminosities and line-of-sight (Doppler) velocities allows one to construct a three-dimensional map of their location and estimate their line-of-sight velocity dispersion. This information, in principle, allows one to test dynamical gravity models, specifically models of satellite galaxy velocity dispersions near massive hosts. A key difficulty is the separation of true satellites from interlopers. We sidestep this problem by not attempting to derive satellite galaxy velocity dispersions from the data, but instead incorporate an interloper background into the mathematical models and compare the result to the actual data. We find that due to the presence of interlopers, it is not possible to exclude several gravitational theories on the basis of the SDSS data.

(abridged) We perform a detailed investigation of sources from the Cold Cores Catalogue of Planck Objects (C3PO). Our goal is to probe the reliability of the detections, validate the separation between warm and cold dust emission components, provide the first glimpse at the nature, internal morphology and physical characterictics of the Planck-detected sources. We focus on a sub-sample of ten sources from the C3PO list, selected to sample different environments, from high latitude cirrus to nearby (150pc) and remote (2kpc) molecular complexes. We present Planck surface brightness maps and derive the dust temperature, emissivity spectral index, and column densities of the fields. With the help of higher resolution Herschel and AKARI continuum observations and molecular line data, we investigate the morphology of the sources and the properties of the substructures at scales below the Planck beam size.

We present the statistical properties of the first version of the Cold Core Catalogue of Planck Objects (C3PO), in terms of their spatial distribution, temperature, distance, mass, and morphology. We also describe the statistics of the Early Cold Core Catalogue (ECC, delivered with the Early Release Compact Source Catalogue, ERCSC) that is the subset of the 915 most reliable detections of the complete catalogue. We have used the CoCoCoDeT algorithm to extract 10783 cold sources. Temperature and dust emission spectral index {\beta} values are derived using the fluxes in the IRAS 100 \mum band and the three highest frequency Planck bands. Temperature spans from 7K to 17K, and peaks around 13K. Data are not consistent with a constant value of {\beta} over the all temperature range. {\beta} ranges from 1.4 to 2.8 with a mean value around 2.1, and several possible scenarios are possible, including {\beta}(T) and the effect of multiple T components folded into the measurements. For one third of the objects the distances are obtained. Most of the detections are within 2 kpc in the Solar neighbourhood, but a few are at distances greater than 4 kpc. The cores are distributed from the deep Galactic plane, despite the confusion, to high latitudes (>30$^{\circle}$). The associated mass estimates range from 1 to $10^5$ solar masses. Using their physical properties these cold sources are shown to be cold clumps, defined as the intermediate cold sub-structures between clouds and cores. These cold clumps are not isolated but mostly organized in filaments associated with molecular clouds. The Cold Core Catalogue of Planck Objects (C3PO) is the first unbiased all-sky catalogue of cold objects. It gives an unprecedented statistical view to the properties of these potential pre-stellar clumps and offers a unique possibility for their classification in terms of their intrinsic properties and environment.

We apply our scalar-tensor-vector (STVG) modified gravity theory (MOG) to calculate the infall velocities of the two clusters constituting the Bullet Cluster 1E0657-06. In the absence of an applicable two-body solution to the MOG field equations, we adopt an approximate acceleration formula based on the spherically symmetric, static, vacuum solution of the theory in the presence of a point source. We find that this formula predicts an infall velocity of the two clusters that is consistent with estimates based on hydrodynamic simulations.

MOG is a fully relativistic modified theory of gravity based on an action principle. The MOG field equations are exactly solvable numerically in two important cases. In the spherically symmetric, static case of a gravitating mass, the equations also admit an approximate solution that closely resembles the Reissner-Nordstrom metric. Furthermore, for weak gravitational fields, a Yukawa-type modification to the Newtonian acceleration law can be obtained, which can be used to model a range of astronomical observations. Without nonbaryonic dark matter, MOG provides good agreement with the data for galaxy rotation curves, galaxy cluster masses, and gravitational lensing, while predicting no appreciable deviation from Einstein's predictions on the scale of the solar system. Another solution of the field equations is obtained for the case of a a spatially homogeneous, isotropic cosmology. MOG predicts an accelerating universe without introducing Einstein's cosmological constant; it also predicts a CMB acoustic power spectrum and a mass power spectrum that are consistent with observations without relying on non-baryonic dark matter. Increased sensitivity in future observations or space-based experiments may be sufficient to distinguish MOG from other theories, notably the LCDM "standard model" of cosmology.

Modified gravity theory is known to violate Birkhoff's theorem. We explore a key consequence of this violation, the effect of distant matter in the Universe on the motion of test particles. We find that when a particle is accelerated, a force is experienced that is proportional to the particle's mass and acceleration and acts in the direction opposite to that of the acceleration. We identify this force with inertia. At very low accelerations, our inertial law deviates slightly from that of Newton, yielding a testable prediction that may be verified with relatively simple experiments. Our conclusions apply to all gravity theories that reduce to a Yukawa-like force in the weak field approximation.

Modified Gravity (MOG) has been used successfully to explain the rotation curves of galaxies, the motion of galaxy clusters, the Bullet Cluster, and cosmological observations without the use of dark matter or Einstein's cosmological constant. We now have the ability to demonstrate how these solutions can be obtained directly from the action principle, without resorting to the use of fitted parameters or empirical formulae. We obtain numerical solutions to the theory's field equations that are exact in the sense that no terms are omitted, in two important cases: the spherically symmetric, static vacuum solution and the cosmological case of an homogeneous, isotropic universe. We compare these results to selected astrophysical and cosmological observations.

A modified gravity (MOG) theory, that has been successfully fitted to galaxy rotational velocity data, cluster data, the Bullet Cluster 1E0657-56 and cosmological observations, is shown to be in good agreement with the motion of satellite galaxies around host galaxies at distances 50-400 kpc.

Globular clusters (GCs) in the Milky Way have characteristic velocity dispersions that are consistent with the predictions of Newtonian gravity, and may be at odds with Modified Newtonian Dynamics (MOND). We discuss a modified gravity (MOG) theory that successfully predicts galaxy rotation curves, galaxy cluster masses and velocity dispersions, lensing, and cosmological observations, yet produces predictions consistent with Newtonian theory for smaller systems, such as GCs. MOG produces velocity dispersion predictions for GCs that are independent of the distance from the galactic center, which may not be the case for MOND. New observations of distant GCs may produce strong criteria that can be used to distinguish between competing gravitational theories.