No Arabic abstract
We perform numerical simulations of gravitational waves (GWs) induced by hydrodynamic and hydromagnetic turbulent sources that might have been present at cosmological quantum chromodynamic (QCD) phase transitions. For turbulent energies of about 4% of the radiation energy density, the typical scale of such motions may have been a sizable fraction of the Hubble scale at that time. The resulting GWs are found to have an energy fraction of about $10^{-9}$ of the critical energy density in the nHz range today and may already have been observed by the NANOGrav collaboration. This is further made possible by our findings of shallower spectra proportional to the square root of the frequency for nonhelical hydromagnetic turbulence. This implies more power at low frequencies than for the steeper spectra previously anticipated. The behavior toward higher frequencies depends strongly on the nature of the turbulence. For vortical hydrodynamic and hydromagnetic turbulence, there is a sharp drop of spectral GW energy by up to five orders of magnitude in the presence of helicity, and somewhat less in the absence of helicity. For acoustic hydrodynamic turbulence, the sharp drop is replaced by a power law decay, albeit with a rather steep slope. Our study supports earlier findings of a quadratic scaling of the GW energy with the magnetic energy of the turbulence and inverse quadratic scaling with the peak frequency, which leads to larger GW energies under QCD conditions.
We discuss the theory of pulsar-timing and astrometry probes of a stochastic gravitational-wave background with a recently developed total-angular-momentum (TAM) formalism for cosmological perturbations. We review the formalism, emphasizing in particular the features relevant for this work and describe the observables we consider (i.e. the pulsar redshift and stellar angular displacement). Using the TAM approach, we calculate the angular power spectra for the observables and from them derive angular auto- and cross-correlation functions. We provide the full set of power spectra and correlation functions not only for the standard transverse-traceless propagating degrees of freedom in general relativity, but also for the four additional non-Einsteinian polarizations that may arise in alternative-gravity theories. We discuss how pulsar-timing and astrometry surveys can complement and serve as cross checks to one another and comment on the importance of testing the chirality of the gravitational-wave background as a tool to understand the nature of its sources. A simple rederivation of the power spectra from the plane-wave formalism is provided in an Appendix.
By regularly monitoring the most stable millisecond pulsars over many years, pulsar timing arrays (PTAs) are positioned to detect and study correlations in the timing behaviour of those pulsars. Gravitational waves (GWs) from supermassive black hole binaries (SMBHBs) are an exciting potentially detectable source of such correlations. We describe a straight-forward technique by which a PTA can be phased-up to form time series of the two polarisation modes of GWs coming from a particular direction of the sky. Our technique requires no assumptions regarding the time-domain behaviour of a GW signal. This method has already been used to place stringent bounds on GWs from individual SMBHBs in circular orbits. Here, we describe the methodology and demonstrate the versatility of the technique in searches for a wide variety of GW signals including bursts with unmodeled waveforms. Using the first six years of data from the Parkes Pulsar Timing Array, we conduct an all-sky search for a detectable excess of GW power from any direction. For the lines of sight to several nearby massive galaxy clusters, we carry out a more detailed search for GW bursts with memory, which are distinct signatures of SMBHB mergers. In all cases, we find that the data are consistent with noise.
The standard model of particle physics is known to be intriguingly successful. However their rich phenomena represented by the phase transitions (PTs) have not been completely understood yet, including the possibility of the existence of unknown dark sectors. In this Letter, we investigate the measurement of the equation of state parameter $w$ and the sound speed $c_{rm s}$ of the PT plasma with use of the gravitational waves (GWs) of the universe. Though the propagation of GW is insensitive to $c_{rm s}$ in itself, the sound speed value affects the dynamics of primordial density (or scalar curvature) perturbations and the induced GW by their horizon reentry can then be an indirect probe both $w$ and $c_{rm s}$. We numerically reveal the concrete spectrum of the predicted induced GW with two simple examples of the scalar perturbation spectrum: the monochromatic and scale-invariant spectra. In the monochromatic case, we see that the resonant amplification and cancellation scales of the induced GW depend on the $c_{rm s}$ values at different time respectively. The scale-invariant case gives a more realistic spectrum and its specific shape will be compared with observations. In particular, the QCD phase transition corresponds with the frequency range of the pulsar timing array (PTA) observations. If the amplitude of primordial scalar power is in the range of $10^{-4}lesssim A_zetalesssim10^{-2}$, the induced GW is consistent with current observational constraints and detectable in the future observation in Square Kilometer Array. Futhermore the recent possible detection of stochastic GWs by NANOGrav 12.5 yr analysis~[1] can be explained by the induced GW if $A_zetasimsqrt{7}times10^{-3}$.
The detection of a stochastic background of low-frequency gravitational waves by pulsar-timing and astrometric surveys will enable tests of gravitational theories beyond general relativity. These theories generally permit gravitational waves with non-Einsteinian polarization modes, which may propagate slower than the speed of light. We use the total-angular-momentum wave formalism to derive the angular correlation patterns of observables relevant for pulsar timing arrays and astrometry that arise from a background of subluminal gravitational waves with scalar, vector, or tensor polarizations. We find that the pulsar timing observables for the scalar longitudinal mode, which diverge with source distance in the luminal limit, are finite in the subluminal case. Furthermore, we apply our results to $f(R)$ gravity, which contains a massive scalar degree of freedom in addition to the standard transverse-traceless modes. The scalar mode in this $f(R)$ theory is a linear combination of the scalar-longitudinal and scalar-transverse modes, exciting only the monopole and dipole for pulsar timing arrays and only the dipole for astrometric surveys.
Models of Dark Matter (DM) can leave unique imprints on the Universes small scale structure by boosting density perturbations on small scales. We study the capability of Pulsar Timing Arrays to search for, and constrain, subhalos from such models. The models of DM we consider are ordinary adiabatic perturbations in $Lambda$CDM, QCD axion miniclusters, models with early matter domination, and vector DM produced during inflation. We show that $Lambda$CDM, largely due to tidal stripping effects in the Milky Way, is out of reach for PTAs (as well as every other probe proposed to detect DM small scale structure). Axion miniclusters may be within reach, although this depends crucially on whether the axion relic density is dominated by the misalignment or string contribution. Models where there is matter domination with a reheat temperature below 1 GeV may be observed with future PTAs. Lastly, vector DM produced during inflation can be detected if it is lighter than $10^{-16} ,{rm GeV}$. We also make publicly available a Python Monte Carlo tool for generating the PTA time delay signal from any model of DM substructure.