The population of Milky Way satellite galaxies is of great interest for cosmology, fundamental physics, and astrophysics. They represent the faint end of the galaxy luminosity function, are the most dark-matter dominated objects in the local Universe, and contain the oldest and most metal-poor stellar populations. Recent surveys have revealed around 60 satellites, but this could represent less than half of the total. Characterization of these systems remains a challenge due to their low luminosity. We consider the gravitational wave observatory LISA as a potential tool for studying these satellites through observations of their short-period double white dwarf populations. LISA will observe the entire sky without selection effects due to dust extinction, complementing optical surveys, and could potentially discover massive satellites hidden behind the disk of the galaxy.

Milky Way dwarf satellites are unique objects that encode the early structure formation and therefore represent a window into the high redshift Universe. So far, their study was conducted using electromagnetic waves only. The future Laser Interpreter Space Antenna (LISA) has the potential to reveal Milky Way satellites in gravitational waves emitted by double white dwarf (DWD) binaries. We investigate gravitational wave (GW) signals detectable by LISA as a possible tool for the identification and characterisation of the Milky Way satellites. We use the binary population synthesis technique to model the population of DWDs in dwarf satellites and we assess the impact on the number of LISA detections when making changes to the total stellar mass, distance, star formation history and metallicity of satellites. We calibrate predictions for the known Milky Way satellites on their observed properties. We find that DWDs emitting at frequencies $\gtrsim 3\,$mHz can be detected in Milky Way satellites at large galactocentric distances. The number of these high frequency DWDs per satellite primarily depends on its mass, distance, age and star formation history, and only mildly depends on the other assumptions regarding their evolution such as metallicity. We find that dwarf galaxies with $M_\star>10^6\,$M$_{\odot}$ can host detectable LISA sources with a number of detections that scales linearly with the satellite's mass. We forecast that out of the known satellites, Sagittarius, Fornax, Sculptor and the Magellanic Clouds can be detected with LISA. As an all-sky survey that does not suffer from contamination and dust extinction, LISA will provide observations of the Milky Way and dwarf satellites galaxies valuable for Galactic archaeology and near-field cosmology.

The label switching problem arises in the Bayesian analysis of models containing multiple indistinguishable parameters with arbitrary ordering. Any permutation of these parameters is equivalent, therefore models with many such parameters have extremely multi-modal posterior distributions. It is difficult to sample efficiently from such posteriors. This paper discusses a solution to this problem which involves carefully mapping the input parameter space to a high dimensional hypertriangle. It is demonstrated that this solution is efficient even for large numbers of parameters and can be easily applied alongside any stochastic sampling algorithm. This method is illustrated using two example problems from the field of gravitational wave astronomy.

Any unambiguous detection of a stochastic gravitational wave background by a pulsar timing array will rest on the measurement of a characteristic angular correlation between pulsars. The ability to measure this correlation will depend on the geometry of the array. However, spatially correlated sources of noise, such as errors in the planetary ephemeris or clock errors, can produce false-positive correlations. The severity of this contamination will also depend on the geometry of the array. This paper quantifies these geometric effects with a spherical harmonic analysis of the pulsar timing residuals. At least 9 well-spaced pulsars are needed to simultaneously measure a gravitational wave background and separate it from ephemeris and clock errors. Uniform distributions of pulsars can eliminate the contamination for arrays with large numbers of pulsars, but pulsars following the galactic distribution of known millisecond pulsars will always be affected. We quantitatively demonstrate the need for arrays to include many pulsars and for the pulsars to be distributed as uniformly as possible. Finally, we suggest a technique to cleanly separate the effect of ephemeris and clock errors from the gravitational wave signal.

We introduce ScannerBit, the statistics and sampling module of the public, open-source global fitting framework GAMBIT. ScannerBit provides a standardised interface to different sampling algorithms, enabling the use and comparison of multiple computational methods for inferring profile likelihoods, Bayesian posteriors, and other statistical quantities. The current version offers random, grid, raster, nested sampling, differential evolution, Markov Chain Monte Carlo (MCMC) and ensemble Monte Carlo samplers. We also announce the release of a new standalone differential evolution sampler, Diver, and describe its design, usage and interface to ScannerBit. We subject Diver and three other samplers (the nested sampler MultiNest, the MCMC GreAT, and the native ScannerBit implementation of the ensemble Monte Carlo algorithm T-Walk) to a battery of statistical tests. For this we use a realistic physical likelihood function, based on the scalar singlet model of dark matter. We examine the performance of each sampler as a function of its adjustable settings, and the dimensionality of the sampling problem. We evaluate performance on four metrics: optimality of the best fit found, completeness in exploring the best-fit region, number of likelihood evaluations, and total runtime. For Bayesian posterior estimation at high resolution, T-Walk provides the most accurate and timely mapping of the full parameter space. For profile likelihood analysis in less than about ten dimensions, we find that Diver and MultiNest score similarly in terms of best fit and speed, outperforming GreAT and T-Walk; in ten or more dimensions, Diver substantially outperforms the other three samplers on all metrics.

In this paper, we propose a new framework for treating the angular information in the pulsar timing array response to a gravitational wave background based on standard cosmic microwave background techniques. We calculate the angular power spectrum of the all-sky gravitational redshift pattern induced at the earth for both a single bright source of gravitational radiation and a statistically isotropic, unpolarized Gaussian random gravitational wave background. The angular power spectrum is the harmonic transform of the Hellings & Downs curve. We use the power spectrum to examine the expected variance in the Hellings & Downs curve in both cases. Finally, we discuss the extent to which pulsar timing arrays are sensitive to the angular power spectrum and find that the power spectrum sensitivity is dominated by the quadrupole anisotropy of the gravitational redshift map.

We use large N-body simulations and empirical scaling relations between dark matter halos, galaxies, and supermassive black holes to estimate the formation rates of supermassive black hole binaries and the resulting low-frequency stochastic gravitational wave background (GWB). We find this GWB to be relatively insensitive ($\lesssim10\%$) to cosmological parameters, with only slight variation between WMAP5 and Planck cosmologies. We find that uncertainty in the astrophysical scaling relations changes the amplitude of the GWB by a factor of $\sim 2$. Current observational limits are already constraining this predicted range of models. We investigate the Poisson variance in the amplitude of the GWB for randomly-generated populations of supermassive black holes, finding a scatter of order unity per frequency bin below 10 nHz, and increasing to a factor of $\sim 10$ near 100 nHz. This variance is a result of the rarity of the most massive binaries, which dominate the signal, and acts as a fundamental uncertainty on the amplitude of the underlying power law spectrum. This Poisson uncertainty dominates at $\gtrsim 20$ nHz, while at lower frequencies the dominant uncertainty is related to our poor understanding of the astrophysical scaling relations, although very low frequencies may be dominated by uncertainties related to the final parsec problem and the processes which drive binaries to the gravitational wave dominated regime. Cosmological effects are negligible at all frequencies.

A large peculiar velocity of the intergalactic medium produces a Doppler shift of the cosmic microwave background with a frequency-dependent quadrupole term. This quadrupole will act as a source for polarization of the cosmic microwave background, creating a large-scale polarization anisotropy if the bulk flow is local and coherent on large scales. In the case where we are near the center of the moving region, the polarization signal is a pure quadrupole. We show that the signal is small, but detectable with future experiments for bulk flows as large as some recent reports.

We perform a detailed and systematic investigation of the possible impacts of dark stars upon the reionization history of the Universe, and its signatures in the cosmic microwave background (CMB). We compute hydrogen reionization histories, CMB optical depths and anisotropy power spectra for a range of stellar populations including dark stars. If dark stars capture large amounts of dark matter via nuclear scattering, reionization can be substantially delayed, leading to decreases in the integrated optical depth to last scattering and large-scale power in the EE polarization power spectrum. Using the integrated optical depth observed by WMAP7, in our canonical reionization model we rule out the section of parameter space where dark stars with high scattering-induced capture rates tie up more than ~90% of all the first star-forming baryons, and live for over ~250 Myr. When nuclear scattering delivers only moderate amounts of dark matter, reionization can instead be sped up slightly, modestly increasing the CMB optical depth. If dark stars do not obtain any dark matter via nuclear scattering, effects upon reionization and the CMB are negligible. The effects of dark stars upon reionization and its CMB markers can be largely mimicked or compensated for by changes in the existing parameters of reionization models, making dark stars difficult to disentangle from astrophysical uncertainties, but also widening the range of standard parameters in reionization models that can be made consistent with observations.