We present a novel method for determining the total matter surface density of the Galactic disk by analysing the kinematics of a dynamically cold stellar stream that passes through or close to the Galactic plane. The method relies on the fact that the vertical component of energy for such stream stars is approximately constant, such that their vertical positions and vertical velocities are interrelated via the matter density of the Galactic disk. By testing our method on mock data stellar streams, with realistic phase-space dispersions and Gaia uncertainties, we demonstrate that it is applicable to small streams out to a distance of a few kilo-parsec, and that the surface density of the disk can be determined to a precision of 6 %. This method is complementary to other mass measurements. In particular, it does not rely on any equilibrium assumption for stars in the Galactic disk, and also makes it possible to measure the surface density to good precision at large distances from the Sun. Such measurements would inform us of the matter composition of the Galactic disk and its spatial variation, place stronger constraints on dark disk sub-structure, and even diagnose possible non-equilibrium effects that bias other types of dynamical mass measurements.

The EDGES low-band experiment has measured an absorption feature in the cosmic microwave background radiation (CMB), corresponding to the 21 cm hyperfine transition of hydrogen at redshift $z \simeq 17$, before the era of cosmic reionization. The amplitude of this absorption is connected to the ratio of singlet and triplet hyperfine states in the hydrogen gas, which can be parametrized by a spin temperature. The EDGES result suggests that the spin temperature is lower than the expected temperatures of both the CMB and the hydrogen gas. A variety of mechanisms have been proposed in order to explain this signal, for example by lowering the kinetic temperature of the hydrogen gas via dark matter interactions. We introduce an alternative mechanism, by which a sub-GeV dark matter particle with spin-dependent coupling to nucleons or electrons can cause hyperfine transitions and lower the spin temperature directly, with negligible reduction of the kinetic temperature of the hydrogen gas. We consider a model with an asymmetric dark matter fermion and a light pseudo-vector mediator. Significant reduction of the spin temperature by this simple model is excluded, most strongly by coupling constant bounds coming from stellar cooling. Perhaps an alternative dark sector model, subject to different sets of constraints, can lower the spin temperature by the same mechanism.

We determine the total dynamical matter density in the solar neighbourhood using the second Gaia data release (DR2). The dynamical matter density distribution is inferred in a framework of a Bayesian hierarchical model, which accounts for position and velocity of all individual stars, as well as the full error covariance matrix of astrometric observables, in a joint fit of the vertical velocity distribution and stellar number density distribution. This was done for eight separate data samples, with different cuts in observed absolute magnitude, each containing about 25,000 stars. The model for the total matter density does not rely on any underlying baryonic model, although we assumed that it is symmetrical, smooth, and monotonically decreasing with distance from the mid-plane. We infer a density distribution which is strongly peaked in the region close to the Galactic plane ($\lesssim 60$ pc), for all eight stellar samples. Assuming a baryonic model and a dark matter halo of constant density, this corresponds to a surplus surface density of approximately 5--9 M$_\odot$/pc$^2$. For the Sun's position and vertical velocity with respect to the Galactic plane, we infer $Z_\odot = 4.76\pm 2.27$ pc and $W_\odot = 7.24\pm 0.19$ km/s. These results suggest a surplus of matter close to the Galactic plane, possibly explained by an underestimated density of cold gas. We discuss possible systematic effects that could bias our result, for example unmodelled non-equilibrium effects, and how to account for such effects in future extensions of this work.

We determine the total dynamical density in the solar neighbourhood using the Tycho-Gaia Astrometric Solution (TGAS) catalogue. Astrometric measurements of proper motion and parallax of stars inform us of both the stellar number density distribution and the velocity distribution of stars close to the plane. Assuming equilibrium, these distributions are interrelated through the local dynamical density. For the first time, we do a full joint fit of the velocity and stellar number density distributions while accounting for the astrometric error of individual stars, in the framework of Bayesian Hierarchical Modelling. We use a sample of stars whose distance extends to approximately 160 pc from the Sun. We find a local matter density of $\rho_0 = 0.119^{+0.015}_{-0.012}~M_\odot$pc$^{-3}$, where the result is presented as the median to the posterior distribution, plus/minus its 16th and 84th percentiles. We find the Sun's position above the Galactic plane to be $z_\odot = 15.3^{+2.24}_{-2.16}$ pc, and the Sun's velocity perpendicular to the Galactic plane to be $W_\odot = 7.19^{+0.18}_{-0.18}$ km/s.

The Gaia mission will provide precise astrometry for an unprecedented number of white dwarfs (WDs), encoding information on stellar evolution, Type Ia supernovae progenitor scenarios, and the star formation and dynamical history of the Milky Way. With such a large data set, it is possible to infer properties of the WD population using only astrometric and photometric information. We demonstrate a framework to accomplish this using a mock data set with SDSS ugriz photometry and Gaia astrometric information. Our technique utilises a Bayesian hierarchical model for inferring properties of a WD population while also taking into account all observational errors of individual objects, as well as selection and incompleteness effects. We demonstrate that photometry alone can constrain the WD population's distributions of temperature, surface gravity and phenomenological type, and that astrometric information significantly improves determination of the WD surface gravity distribution. We also discuss the possibility of identifying unresolved binary WDs using only photometric and astrometric information.

Multiple stellar systems are ubiquitous in the Milky Way, but are often unresolved and seen as single objects in spectroscopic, photometric, and astrometric surveys. Yet, modeling them is essential for developing a full understanding of large surveys such as Gaia, and connecting them to stellar and Galactic models. In this paper we address this problem by jointly fitting the Gaia and 2MASS photometric and astrometric data using a data-driven Bayesian hierarchical model that includes populations of binary and trinary systems. This allows us to classify observations into singles, binaries, and trinaries, in a robust and efficient manner, without resorting to external models. We are able to identify multiple systems and, in some cases, make strong predictions for the properties of its unresolved stars. We will be able to compare such predictions with Gaia Data Release 4, which will contain astrometric identification and analysis of binary systems.

I study the process of dark matter capture by the Sun, under the assumption of a Weakly Interacting Massive Particle (WIMP), in the framework of non-relativistic effective field theory. Hypothetically, WIMPs from the galactic halo can scatter against atomic nuclei in the solar interior, settle to thermal equilibrium with the solar core and annihilate to produce an observable flux of neutrinos. In particular, I examine the thermalization process using Monte-Carlo integration of WIMP trajectories. I consider WIMPs in a mass range of 10--1000 GeV and WIMP-nucleon interaction operators with different dependence on spin and transferred momentum. I find that the density profiles of captured WIMPs are in accordance with a thermal profile described by the Sun's gravitational potential and core temperature. Depending on the operator that governs the interaction, the majority of the thermalization time is spent in either the solar interior or exterior. If normalizing the WIMP-nuclei interaction strength to a specific capture rate, I find that the thermalization time differs at most by 3 orders of magnitude between operators. In most cases of interest, the thermalization time is many orders of magnitude shorter than the age of the solar system.

We study the capture of WIMP dark matter by the Sun in the non-relativistic effective theory of dark matter self-interactions. The aim is to assess the impact of self-interactions on the expected neutrino flux from the annihilation of WIMPs trapped in the Sun in a model independent manner. We consider all non-relativistic Galilean invariant self-interaction operators that can arise from the exchange of a heavy particle of spin less than or equal to 1 for WIMPs of spin equal to 0, 1/2 and 1. We show that for interaction operators depending at most linearly on the momentum transfer, the WIMP-induced neutrino flux can be enhanced by several orders of magnitude compared to the same flux in absence of self-interactions. This is true even for standard values of the thermally averaged annihilation cross-section. This conclusion impacts the analysis of present and future observations performed at neutrino telescopes.