[Highly abridged, from executive summary] As much as NewSpace presents opportunities, there are significant challenges that must be overcome, requiring engagement with policy makers to influence domestic and international space governance. Failure to do so could result in a range of long-lasting negative outcomes for science and space stewardship. How will the Canadian astronomical community engage with NewSpace? What are the implications for NewSpace on the astro-environment, including Earth orbits, lunar and cis-lunar orbits, and surfaces of celestial bodies? This white paper analyzes the rapid changes in space use and what those changes could mean for Canadian astronomers. Our recommendations are as follows: Greater cooperation between the astronomical and the Space Situational Awareness communities is needed. Build closer ties between the astronomical community and Global Affairs Canada (GAC). Establish a committee for evaluating the astro-environmental impacts of human space use, including on and around the Moon and other bodies. CASCA and the Tri-Council should coordinate to identify programs that would enable Canadian astronomers to participate in pay-for-use services at appropriate funding levels. CASCA should continue to foster a relationship with CSA, but also build close ties to the private space industry. Canadian-led deep space missions are within Canada's capabilities, and should be pursued.

We model the fastest moving (v_tot > 300 km/s) local (D < 3 kpc) halo stars using cosmological simulations and 6-dimensional Gaia data. Our approach is to use our knowledge of the assembly history and phase-space distribution of halo stars to constrain the form of the high velocity tail of the stellar halo. Using simple analytical models and cosmological simulations, we find that the shape of the high velocity tail is strongly dependent on the velocity anisotropy and number density profile of the halo stars --- highly eccentric orbits and/or shallow density profiles have more extended high velocity tails. The halo stars in the solar vicinity are known to have a strongly radial velocity anisotropy, and it has recently been shown the origin of these highly eccentric orbits is the early accretion of a massive (M_star ~ 10^9 M_Sun) dwarf satellite. We use this knowledge to construct a prior on the shape of the high velocity tail. Moreover, we use the simulations to define an appropriate outer boundary of 2r_200, beyond which stars can escape. After applying our methodology to the Gaia data, we find a local (r_0=8.3 kpc) escape speed of v_esc(r_0) = 528(+24,-25) km/s. We use our measurement of the escape velocity to estimate the total Milky Way mass, and dark halo concentration: M_200,tot = 1.00(+0.31,-0.24) x 10^12 M_Sun, c_200 = 10.9(+4.4,-3.3). Our estimated mass agrees with recent results in the literature that seem to be converging on a Milky Way mass of M_200,tot ~ 10^12 M_Sun.

In the rest frame of the Local Group (LG), the total momentum of the Milky Way (MW) and Andromeda (M31) should balance to zero. We use this fact to constrain new solutions for the solar motion with respect to the LG centre-of-mass, the total mass of the LG, and the individual masses of M31 and the MW. Using the set of remote LG galaxies at $>350$ kpc from the MW and M31, we find that the solar motion has amplitude $V_{\odot}=299\pm 15 {\rm ~km~s^{-1}}$ in a direction pointing toward galactic longitude $l_{\odot}=98.4^{\circ}\pm 3.6^{\circ}$ and galactic latitude $b_{\odot}=-5.9^{\circ}\pm 3.0^{\circ}$. The velocities of M31 and the MW in this rest frame give a direct measurement of their mass ratio, for which we find $\log_{10} (M_{\rm M31}/M_{\rm MW})=0.36 \pm 0.29$. We combine these measurements with the virial theorem to estimate the total mass within the LG as $M_{\rm LG}=(2.5\pm 0.4)\times 10^{12}~{\rm M}_{\odot}$. Our value for $M_{\rm LG}$ is consistent with the sum of literature values for $M_{\rm MW}$ and $M_{\rm M31}$. This suggests that the mass of the LG is almost entirely located within the two largest galaxies rather than being dispersed on larger scales or in a background medium. The outskirts of the LG are seemingly rather empty. Combining our measurement for $M_{\rm LG}$ and the mass ratio, we estimate the individual masses of the MW and M31 to be $M_{\rm MW}=(0.8\pm 0.5)\times 10^{12}~{\rm M}_{\odot}$ and $M_{\rm M31}=(1.7\pm 0.3)\times 10^{12}~{\rm M}_{\odot}$, respectively. Our analysis favours M31 being more massive than the MW by a factor of $\sim$2.3, and the uncertainties allow only a small probability (9.8%) that the MW is more massive. This is consistent with other properties such as the maximum rotational velocities, total stellar content, and numbers of globular clusters and dwarf satellites, which all suggest that $M_{\rm M31}/M_{\rm MW}>1$.

The Sagittarius dwarf galaxy is the archetype of a tidally disrupting system. Both leading and trailing tails can be observed across at least 180 degrees of the sky and measurements of their luminosity density profiles have recently become available. Using numerical simulations, we explore the factors that control the appearance of such profiles. We use two possible models for the Sgr progenitor. The first is a one-component Plummer model, which may represent either a dark matter free progenitor, or one in which pre-existing dark matter has already been largely stripped. The second is a two-component model in which the stars are represented by a Hernquist sphere embedded in a cosmologically modish Navarro-Frenk-White dark halo. Disruption of the models in the Milky Way galaxy provides us with two tellings of the tale of the formation of the Sgr stream. The initial disintegration of the baryons proceeds more slowly for the two-component models because of the protective cocoon of dark matter. Once this has been stripped, though, matters proceed apace. In both cases, the profiles after ~6 pericentric passages provide good matches to the observational data, but the tails are more extended for the two-component models. The leading and trailing tails are symmetric at apocentre or pericentre. At other orbital phases, asymmetries are present, as tails are compressed as they approach apocentre and stretched out as they approach pericentre. There may exist density enhancements corresponding to such pile-ups which may be observable in current survey data. We re-visit the calculation of Niederste-Ostholt et al. (2010) and slightly revise upwards the luminosity of the Sgr progenitor to 9.9-14.4 x 10^7 solar luminosities based on insights from the simulations.

Some of the most massive globular clusters of our Milky Way, such as for example omega-Centauri, show a mixture of stellar populations spanning a few Gyr in age and 1.5 dex in metallicities. In contrast, standard formation scenarios predict that globular and open clusters form in one single star-burst event of duration less than about 10 Myr and therefore should exhibit only one age and one metallicity in its stars. Here, we investigate the possibility that a massive stellar super-cluster may trap older galactic field stars during its formation process that are later detectable in the cluster as an apparent population of stars with a very different age and metallicity. With a set of numerical N-body simulations, we are able to show that, if the mass of the stellar super-cluster is high enough and the stellar velocity dispersion in the cluster is comparable to the dispersion of the surrounding disc stars in the host galaxy, then up to about 40 per cent of its initial mass can be additionally gained from trapped disc stars. We also suggest that a super-cluster may capture in excess of its own mass under certain conditions.

We increase the number of remote halo tracers by using blue horizontal branch (BHB) stars out to Galactocentric distances of 130 kpc. We use SDSS EDR photometry and the VLT to detect 16 BHB stars at Galactocentric distances 70 < r < 130, and to measure their radial velocities. We find the mass of the Milky Way is M =1.7^{+3.0}_{-0.6} x 10^{12} solar masses. When completed this survey will: (i) substantially reduce the errors in the total mass and extent of the Milky Way halo, (ii) map the velocity space in a hitherto unexplored region of the halo.

Ultra-high energy cosmic rays (UHECRs) may originate from the decay of massive relic particles in the dark halo of the Galaxy, or they may be produced by supermassive black holes in the nuclei of nearby galaxies. The anisotropy in the arrival directions is studied in four dark halo models (cusped, isothermal, triaxial and tilted) and in four galaxy samples (galaxies intrinsically brighter than Centaurus A within 50 and 100 Mpc, and brighter than M32 within 50 and 100 Mpc). In decaying dark matter models, the amplitude of the anisotropy is controlled by the size of the Galactic halo, while the phase is controlled by the shape. In the northern hemisphere, the amplitude is about 0.5 for cusped haloes, but falls to roughly 0.3 for isothermal haloes. The phase points in the direction of the Galactic Centre, with deviations of up to 30 degrees possible for triaxial and tilted haloes. The effect of the halo of M31 is too weak to provide conclusive evidence for the decaying dark matter origin of UHECRs. In extragalactic models, samples of galaxies brighter than Centaurus A produce substantial anisotropies (roughly 1.8), much larger than the limits set by the available data. If all galaxies brighter than M32 contribute, then the anisotropy is more modest (roughly 0.5) and is directed toward mass concentrations in the supergalactic plane, like the Virgo cluster. Predictions are made for the southern hemisphere station of the Pierre Auger Observatory. If the UHECRs have a Galactic origin, then the phase points towards the Galactic Centre. If they have an extragalactic origin, then it points in the rough direction of the Fornax cluster. This provides an unambiguous discriminant and requires about 350-500 events at South Auger.

This paper studies multistep methods for the integration of reversible dynamical systems, with particular emphasis on the planar Kepler problem. It has previously been shown by Cano & Sanz-Serna that reversible linear multisteps for first-order differential equations are generally unstable. Here, we report on a subset of these methods -- the zero-growth methods -- that evade these instabilities. We provide an algorithm for identifying these rare methods. We find and study all zero-growth, reversible multisteps with six or fewer steps. This select group includes two well-known second-order multisteps (the trapezoidal and explicit midpoint methods), as well as three new fourth-order multisteps -- one of which is explicit. Variable timesteps can be readily implemented without spoiling the reversibility. Tests on Keplerian orbits show that these new reversible multisteps work well on orbits with low or moderate eccentricity, although at least 100 steps/radian are required for stability.

The recent years have witnessed a carnival of discoveries in the Outer Solar System. Here we provide evidence from numerical simulations of orbital stability to suggest the possible existence of two long lived belts of primordial planetesimals in the Solar System. The first is the domain of the Vulcanoids (between 0.09 - 0.21 AU) between the Sun and Mercury, where remnant planetesimals may survive on dynamically stable orbits provided they possess a characteristic radius greater than 0.1 km. The second is a belt between the Earth and Mars (1.08 -1.28 AU) on which an initial population of particles on circular orbits may survive for the age of the Solar System. A search through the catalogues of Near-Earth Objects reveals an excess of asteroids with low eccentricities and inclinations occupying this belt, such as the recently discovered objects 1996 XB27, 1998 HG49 and 1998 KG3.

Scale-free disks have no preferred length or time scale. The question has been raised whether such disks have a continuum of unstable linear modes or perhaps no unstable modes at all. We resolve this paradox by analysing the particular case of a gaseous, isentropic disk with a completely flat rotation curve (the Mestel disk) exactly. The heart of the matter is this: what are the correct boundary conditions to impose at the origin or central cusp? We argue that the linear stability problem is ill-posed. From any finite radius, waves reach the origin after finite time but with logarithmically divergent phase. Instabilities exist, but their pattern speeds depend upon an undetermined phase with which waves are reflected from the origin. For any definite choice of this phase, there is an infinite but discrete set of growing modes. Similar ambiguities may afflict general disk models with power-law central cusps. The ratio of growth rate to pattern speed, however, is independent of the central phase. This ratio is derived in closed form for non self-gravitating normal modes. The ratio for self-gravitating normal modes is found numerically by solving recurrence relations in Mellin-transform space.

Lens models appropriate for representing cusped galaxies and clusters are developed. The analogue of the odd number theorem for cusped density distributions is given. Density cusps are classified into strong, isothermal or weak, according to their lensing properties. Strong cusps cause multiple imaging for any source position, whereas isothermal and weak cusps give rise to only one image for distant sources. Isothermal cusps always possess a pseudo-caustic. When the source crosses the pseudo-caustic, the number of images changes by unity.Two families of cusped galaxy and cluster models are examined in detail. The double power-law family has an inner cusp, followed by a transition region and an outer envelope. One member of this family -- the isothermal double power-law model -- possesses an exceedingly scarce property, namely, the lens equation is exactly solvable for any source position. This means that the magnifications, the time delay and the lensing cross-sections are all readily available. The model has a three dimensional density that is cusped like $r^{-2}$ at small radii and falls off like $r^{-4}$ asymptotically. Thus, it provides a very useful representation of the lensing properties of a galaxy or cluster of finite total mass with a flat rotation curve. The second set of models studied is the single power-law family. These are single density cusps of infinite extent. The properties of the critical curves and caustics and the behaviour of the lenses in the presence of external shear are all discussed in some detail.

Microlensing has established itself as a powerful new method for the detection of baryonic dark matter in the Galaxy. The theory of microlensing is sketched and its similarity with the optical effect of twinkling is explained. The bulk of the article presents a new analysis of the data-set on microlensing towards the Large Magellanic Cloud. The extent, flattening and velocity anisotropy of the Galactic halo are unknown. So, it is vital to analyse the microlensing data-set with families of models that span the viable ranges of these structural parameters. Also crucial is proper modelling of the dark matter halo of the Large Magellanic Cloud. Despite all the unknowns, a robust conclusion is that the Galactic and LMC haloes cannot be primarily built from objects in the mass range 10^{-7} - 0.1 solar masses. The baryonic fraction of the Galactic halo probably lies in the range 0.4 - 0.5. The mass of the deflectors probably lies between 0.05 - 1.0 solar masses. Stronger claims concerning the masses of the dark objects are unwarranted because the estimates are sensitive to the uncertain velocity anisotropy.

In a recent Letter (astro-ph/9605001), Cowsik et al. claim that the best fit value of the velocity dispersion of the Galactic dark matter is 600 kms. This result is obtained by a numerical method of solving the coupled self-gravitation equations, assuming Maxwellian velocity distributions. The rotation curve presented in their original paper is not the whole story, as it does not match onto the spherical self-gravitating halo with a Maxwellian velocity distribution -- as it must. This was pointed out in the original version of this Comment. In their Reply (astro-ph/9608035), Cowsik et al. provide a rotation curve that does behave with the correct asymptotics. Unfortunately, it violates mass estimates of the Local Group by over an order of magnitude. It is the pressure required to balance the weight of the overlying layers of this phenomenal mass at large radii that is the physical origin of Cowsik et al's anomalously high dark matter velocity dispersion.