The Universe is homogeneous and isotropic on large scales, so on those scales it is usually modelled as a Friedmann-Lema\^{i}tre-Robertson-Walker (FLRW) space-time. The non-linearity of the Einstein field equations raises concern over averaging over small-scale deviations form homogeneity and isotropy, with possible implications on the applicability of the FLRW metric to the Universe, even on large scales. Here I present a technique, based on the multiple-scales method of singular perturbation theory, to handle the small-scale inhomogeneities consistently. I obtain a leading order effective Einstein equation for the large-scale space-time metric, which contains a back-reaction term. The derivation of this equation is done in harmonic gauge, and conversion to other gauges is discussed. I estimate the magnitude of the back-reaction term relative to the critical density of the Universe in the example of an NFW halo, and find it to be of the order of a few percent. In this example, the back-reaction term is interpreted as a contribution of the energy-density of gravitational potential energy, averaged over the small-scale, to the total energy-momentum tensor.

We investigate the extent to which resonances between an oscillating background of ultra-light axion and a binary Keplerian system can affect the motion of the latter. These resonances lead to perturbations in the instantaneous time-of-arrivals, and to secular variations in the period of the binary. While the secular changes at exact resonance have recently been explored, the instantaneous effects have been overlooked. In this paper, we examine the latter using N-body simulations including the external oscillatory forcing induced by the axion background. While the secular effects are restricted to a narrow width near the resonance, the instantaneous changes, albeit strongest close to resonances, are apparent for wide range of configurations. We compute the signal-to-noise ratio (SNR) as a function of semi-major axis for a detection of axion oscillations through the R\{o} mer delay. The latter can be extracted from the time-of-arrivals of binary pulsars. The SNR broadly increases with increasing binary eccentricity in agreement with the secular expectation. However, we find that it differs significantly from the scaling a^{5/2} around the lowest orders of resonance. Future observations could probe these effects away from resonances and, therefore, constrain a much broader range of axion masses provided that binary pulsar systems are found near the central region of our Galaxy, and that the time-or-arrival measurement accuracy reaches < 10 ns

Scalar, tensor waves induce oscillatory perturbations in Keplerian systems which can be probed with measurements of pulsar timing residuals. In this paper, we consider the imprint of coherent oscillations produced by ultra-light axion dark matter on the Roemer time delay. We use the angle-action formalism to calculate the time evolution of the observed signal and its dependence on the orbital parameters and the axion phase. We emphasize the similarity of the expected signal-to-noise with the response of an harmonic oscillator to an external oscillatory driving. We validate our theoretical predictions with numerical simulations. Our results furnish a useful benchmark for numerical codes and analysis procedures and, hopefully, will motivate the search for such imprints in real data.

The coalescence of compact binary stars is expected to produce a stochastic background of gravitational waves (GW) observable with future GW detectors. Such backgrounds are usually characterized by their power spectrum as a function of frequency. Here, we present a method to calculate the full 1-point distribution of strain fluctuations. We focus on time series data, but our approach generalizes to the frequency domain. We illustrate how this probability distribution can be evaluated numerically. In addition, we derive accurate analytical asymptotic expressions for the large strain tail, which demonstrate that it is dominated by the nearest source. As an application, we also calculate the distribution of strain fluctuations for the astrophysical GW background produced by binary mergers of compact stars in the Universe, and quantify the extent to which it deviates from a Gaussian distribution. Our approach could be useful for the spectral shape reconstruction of stochastic GW backgrounds.

Detection of gravitational-wave (GW) sources enables the characterisation of binary compact objects and of their in-spiral. However, other dissipative processes can affect the in-spiral. Here we show that the in-spiral of compact objects through a gaseous common-envelope (CE) arising from an evolved stellar companion produces a novel type of GW-sources, whose evolution is dominated by the dissipative gas dynamical friction effects from the CE, rather than the GW-emission itself. The evolution and properties of the GW-signals differ from those of isolated gas-poor mergers significantly. We find characteristic strains of $\sim10^{-23}$-$10^{-21}$ ($10{\rm kpc}/{D}$) for such sources -- observable by next-generation space-based GW-detectors. The evolution of the GW-signal can serve as a probe of the interior regions of the evolved star, and the final stages of CE-evolution, otherwise inaccessible through other observational means. Moreover, such CE-mergers are frequently followed by observable explosive electromagnetic counterparts and/or the formation of exotic stars.

We analyse the scatter in the correlation between super-massive black hole (SMBH) mass and bulge stellar mass of the host galaxy, and infer that it cannot be accounted for by mergers alone. The merger-only scenario, where small galaxies merge to establish a proportionality relation between the SMBH and bulge masses, leads to a scatter around the linear proportionality line that increases with the square root of the SMBH (or bulge) mass. By examining a sample of 103 galaxies we find that the intrinsic scatter increases more rapidly than expected from the merger-only scenario. The correlation between SMBH masses and their host galaxy properties is therefore more likely to be determined by a negative feedback mechanism that is driven by an active galactic nucleus. We find a hint that some galaxies with missing stellar mass reside close to the centre of clusters and speculate that ram-pressure stripping of gas off the young galaxy as it moves near the cluster centre, might explain the missing stellar mass at later times.