Showing 1–50 of 100 results for author: Hindmarsh, M

Recovering a phase transition signal in simulated LISA data with a modulated galactic foreground

Authors: Mark Hindmarsh, Deanna C. Hooper, Tiina Minkkinen, David J. Weir

Abstract: Stochastic backgrounds of gravitational waves from primordial first-order phase transitions are a key probe of physics beyond the Standard Model. They represent one of the best prospects for observing or constraining new physics with the LISA gravitational wave observatory. However, the large foreground population of galactic binaries in the same frequency range represents a challenge, and will hi… ▽ More Stochastic backgrounds of gravitational waves from primordial first-order phase transitions are a key probe of physics beyond the Standard Model. They represent one of the best prospects for observing or constraining new physics with the LISA gravitational wave observatory. However, the large foreground population of galactic binaries in the same frequency range represents a challenge, and will hinder the recovery of a stochastic background. To test the recoverability of a stochastic gravitational wave background, we use the LISA Simulation Suite to generate data incorporating both a stochastic background and an annually modulated foreground modelling the galactic binary population, and the Bayesian analysis code Cobaya to attempt to recover the model parameters. By applying the Deviance Information Criterion to compare models with and without a stochastic background we place bounds on the detectability of gravitational waves from first-order phase transitions. By further comparing models with and without the annual modulation, we show that exploiting the modulation improves the goodness-of-fit and gives a modest improvement to the bounds on detectable models. △ Less

Submitted 7 June, 2024; originally announced June 2024.

Comments: 29 pages, 9 figures

Report number: HIP-2024-8/TH

Revised bounds on local cosmic strings from NANOGrav observations

Authors: Jun'ya Kume, Mark Hindmarsh

Abstract: In a recent paper, the NANOGrav collaboration studied new physics explanations of the observed pulsar timing residuals consistent with a stochastic gravitational wave background (SGWB), including cosmic strings in the Nambu-Goto (NG) approximation. Analysing one of current models for the loop distribution, it was found that the cosmic string model is disfavored compared to other sources, for examp… ▽ More In a recent paper, the NANOGrav collaboration studied new physics explanations of the observed pulsar timing residuals consistent with a stochastic gravitational wave background (SGWB), including cosmic strings in the Nambu-Goto (NG) approximation. Analysing one of current models for the loop distribution, it was found that the cosmic string model is disfavored compared to other sources, for example, super massive black hole binaries (SMBHBs). When both SMBHB and cosmic string models are included in the analysis, an upper bound on a string tension $Gμ\lesssim 10^{-10}$ was derived. However, the analysis did not accommodate results from cosmic string simulations in an underlying field theory, which indicate that at most a small fraction of string loops survive long enough to emit GW. Following and extending our previous study, we suppose that a fraction $f_{\rm NG}$ of string loops follow NG dynamics and emit only GWs, and study the three different models of the loop distribution discussed in the LIGO-Virgo-KAGRA (LVK) collaboration analyses. We re-analyze the NANOGrav 15yrs data with our signal models by using the NANOGrav $\texttt{ENTERPRISE}$ analysis code via the wrapper $\texttt{PTArcade}$. We find that loop distributions similar to LVK Model B and C yield higher Bayes factor than Model A analyzed in the NANOGrav paper, as they can more easily accommodate a blue-tilted spectrum of the observed amplitude. Furthermore, because of the degeneracy of $Gμ$ and $f_{\rm NG}$ in determining the signal amplitude, our posterior distribution extends to higher values of $Gμ$, and in some cases the uppermost value of credible intervals is close to the Cosmic Microwave Background limit $Gμ\lesssim 10^{-7}$. Hence, in addition to the pulsar timing array data, further information about the fraction of long-lived loops in a cosmic string network is required to constrain the string tension. △ Less

Submitted 15 April, 2024; v1 submitted 3 April, 2024; originally announced April 2024.

Comments: 22 pages, 9 figures; v2: references are added

Report number: RESCEU-4/24, HIP-2024-11/TH

A-B transition in superfluid $^3$He and cosmological phase transitions

Authors: Mark Hindmarsh, J. A. Sauls, Kuang Zhang, S. Autti, Richard P. Haley, Petri J. Heikkinen, Stephan J. Huber, Lev V. Levitin, Asier Lopez-Eiguren, Adam J. Mayer, Kari Rummukainen, John Saunders, Dmitry Zmeev

Abstract: First order phase transitions in the very early universe are a prediction of many extensions of the Standard Model of particle physics and could provide the departure from equilibrium needed for a dynamical explanation of the baryon asymmetry of the Universe. They could also produce gravitational waves of a frequency observable by future space-based detectors such as the Laser Interferometer Space… ▽ More First order phase transitions in the very early universe are a prediction of many extensions of the Standard Model of particle physics and could provide the departure from equilibrium needed for a dynamical explanation of the baryon asymmetry of the Universe. They could also produce gravitational waves of a frequency observable by future space-based detectors such as the Laser Interferometer Space Antenna (LISA). All calculations of the gravitational wave power spectrum rely on a relativistic version of the classical nucleation theory of Cahn-Hilliard and Langer, due to Coleman and Linde. The high purity and precise control of pressure and temperature achievable in the laboratory made the first-order A to B transition of superfluid $^3$He an ideal for test of classical nucleation theory. As Leggett and others have noted the theory fails dramatically. The lifetime of the metastable A phase is measurable, typically of order minutes to hours, far faster than classical nucleation theory predicts. If the nucleation of B phase from the supercooled A phase is due to a new, rapid intrinsic mechanism that would have implications for first-order cosmological phase transitions as well as predictions for gravitational wave (GW) production in the early universe. Here we discuss studies of the AB phase transition dynamics in $^3$He, both experimental and theoretical, and show how the computational technology for cosmological phase transition can be used to simulate the dynamics of the A-B transition, support the experimental investigations of the A-B transition in the QUEST-DMC collaboration with the goal of identifying and quantifying the mechanism(s) responsible for nucleation of stable phases in ultra-pure metastable quantum phases. △ Less

Submitted 15 January, 2024; originally announced January 2024.

Comments: 34 pp, 4 figures. Based on a talk given at Quantum Fluids and Solids 2023, Manchester, U.K, on behalf of the QUEST-DMC collaboration

Report number: HIP-2024-3/TH

Nanofluidic platform for studying the first-order phase transitions in superfluid helium-3

Authors: Petri J. Heikkinen, Nathan Eng, Lev V. Levitin, Xavier Rojas, Angadjit Singh, Samuli Autti, Richard P. Haley, Mark Hindmarsh, Dmitry E. Zmeev, Jeevak M. Parpia, Andrew Casey, John Saunders

Abstract: The symmetry-breaking first-order phase transition between superfluid phases $^3$He-A and $^3$He-B can be triggered extrinsically by ionising radiation or heterogeneous nucleation arising from the details of the sample cell construction. However, the role of potential homogeneous intrinsic nucleation mechanisms remains elusive. Discovering and resolving the intrinsic processes may have cosmologica… ▽ More The symmetry-breaking first-order phase transition between superfluid phases $^3$He-A and $^3$He-B can be triggered extrinsically by ionising radiation or heterogeneous nucleation arising from the details of the sample cell construction. However, the role of potential homogeneous intrinsic nucleation mechanisms remains elusive. Discovering and resolving the intrinsic processes may have cosmological consequences, since an analogous first-order phase transition, and the production of gravitational waves, has been predicted for the very early stages of the expanding Universe in many extensions of the Standard Model of particle physics. Here we introduce a new approach for probing the phase transition in superfluid $^3$He. The setup consists of a novel stepped-height nanofluidic sample container with close to atomically smooth walls. The $^3$He is confined in five tiny nanofabricated volumes and assayed non-invasively by NMR. Tuning of the state of $^3$He by confinement is used to isolate each of these five volumes so that the phase transitions in them can occur independently and free from any obvious sources of heterogeneous nucleation. The small volumes also ensure that the transitions triggered by ionising radiation are strongly suppressed. Here we present the preliminary measurements using this setup, showing both strong supercooling of $^3$He-A and superheating of $^3$He-B, with stochastic processes dominating the phase transitions between the two. The objective is to study the nucleation as a function of temperature and pressure over the full phase diagram, to both better test the proposed extrinsic mechanisms and seek potential parallel intrinsic mechanisms. △ Less

Submitted 29 May, 2024; v1 submitted 11 January, 2024; originally announced January 2024.

Comments: 18 pages, 6 figures, part of the Proceedings of the International Conference on Quantum Fluids and Solids 2023. This is a copy of the version of record of this article, first published in Journal of Low Temperature Physics, available online at Publisher's website: https://doi.org/10.1007/s10909-024-03146-6

Shallow relic gravitational wave spectrum with acoustic peak

Authors: Ramkishor Sharma, Jani Dahl, Axel Brandenburg, Mark Hindmarsh

Abstract: We study the gravitational wave (GW) spectrum produced by acoustic waves in the early universe, such as would be produced by a first order phase transition, focusing on the low-frequency side of the peak. We confirm with numerical simulations the Sound Shell model prediction of a steep rise with wave number $k$ of $k^9$ to a peak whose magnitude grows at a rate $(H/k_\text{p})H$, where $H$ is the… ▽ More We study the gravitational wave (GW) spectrum produced by acoustic waves in the early universe, such as would be produced by a first order phase transition, focusing on the low-frequency side of the peak. We confirm with numerical simulations the Sound Shell model prediction of a steep rise with wave number $k$ of $k^9$ to a peak whose magnitude grows at a rate $(H/k_\text{p})H$, where $H$ is the Hubble rate and $k_\text{p}$ the peak wave number, set by the peak wave number of the fluid velocity power spectrum. We also show that hitherto neglected terms give a shallower part with amplitude $(H/k_\text{p})^2$ in the range $H \lesssim k \lesssim k_\text{p}$, which in the limit of small $H/k$ rises as $k$. This linear rise has been seen in other modelling and also in direct numerical simulations. The relative amplitude between the linearly rising part and the peak therefore depends on the peak wave number of the velocity spectrum and the lifetime of the source, which in an expanding background is bounded above by the Hubble time $H^{-1}$. For slow phase transitions, which have the lowest peak wave number and the loudest signals, the acoustic GW peak appears as a localized enhancement of the spectrum, with a rise to the peak less steep than $k^9$. The shape of the peak, absent in vortical turbulence, may help to lift degeneracies in phase transition parameter estimation at future GW observatories. △ Less

Submitted 10 November, 2023; v1 submitted 24 August, 2023; originally announced August 2023.

Comments: 21 pages, 10 figures, Accepted for publication in JCAP

Report number: NORDITA-2023-051, HIP-2023-13/TH

Journal ref: J. Cosmol. Astropart. Phys. 12 (2023) 042

General relativistic bubble growth in cosmological phase transitions

Authors: Lorenzo Giombi, Mark Hindmarsh

Abstract: We use a full general relativistic framework to study the self-similar expansion of bubbles of the stable phase into a flat Friedmann-Lemaître-Robertson-Walker Universe in a first order phase transition in the early Universe. With a simple linear barotropic equation of state in both phases, and in the limit of a phase boundary of negligible width, we find that self-similar solutions exist, which a… ▽ More We use a full general relativistic framework to study the self-similar expansion of bubbles of the stable phase into a flat Friedmann-Lemaître-Robertson-Walker Universe in a first order phase transition in the early Universe. With a simple linear barotropic equation of state in both phases, and in the limit of a phase boundary of negligible width, we find that self-similar solutions exist, which are qualitatively similar to the analogous solutions in Minkowski space, but with distinguishing features. Rarefaction waves extend to the centre of the bubble, while spatial sections near the centre of the bubble have negative curvature. Gravitational effects redistribute the kinetic energy of the fluid around the bubble, and can change the kinetic energy fraction significantly. The kinetic energy fraction of the gravitating solution can be enhanced over the analogous Minkowski solution by as much as $\mathcal{O}(1)$, and suppressed by a factor as larger as $\mathcal{O}(10)$ in case of fast detonations. The amount of negative spatial curvature at the centre of the bubble is of the same order of magnitude of the naive expectation based on considerations of the energy density perturbation in Minkowski solutions, with gravitating deflagrations less negatively curved, and detonations more. We infer that general relativistic effects might have a significant impact on accurate calculations of the gravitational wave power spectrum when the bubble size becomes comparable to the cosmological Hubble radius, affecting the primary generation from the fluid shear stress, and inducing secondary generation by scalar perturbations. △ Less

Submitted 29 August, 2023; v1 submitted 22 July, 2023; originally announced July 2023.

Comments: 42 pages, 11 figures

Report number: HIP-2023-12/TH

Multi-messenger constraints on Abelian-Higgs cosmic string networks

Authors: Mark Hindmarsh, Jun'ya Kume

Abstract: Nielsen-Olesen vortices in the Abelian-Higgs (AH) model are the simplest realisations of cosmic strings in a gauge field theory. Large-scale numerical solutions show that the dominant decay channel of a network of AH strings produced from random initial conditions is classical field radiation. However, they also show that with special initial conditions, loops of string can be created for which cl… ▽ More Nielsen-Olesen vortices in the Abelian-Higgs (AH) model are the simplest realisations of cosmic strings in a gauge field theory. Large-scale numerical solutions show that the dominant decay channel of a network of AH strings produced from random initial conditions is classical field radiation. However, they also show that with special initial conditions, loops of string can be created for which classical field radiation is suppressed, and which behave like Nambu-Goto (NG) strings with a dominant decay channel into gravitational radiation. This indicates that cosmic strings are generically sources of both high-energy particles and gravitational waves. Here we adopt a simple parametrisation of the AH string network allowing for both particle and gravitational wave production. With a reference to a specific model for NG-like loop distribution, this sets the basis for a ``multi-messenger'' investigation of this model. We find that, in order to explain the NANOGrav detection of a possible gravitational wave background, while satisfying the constraint on NG-like loop production from simulations and bounds from the cosmic microwave background, the tension of the AH string in Planck units $Gμ$ and the fraction of the NG-like loops $f_{\rm NG}$ should satisfy $Gμf_{\rm NG}^{2.6} \gtrsim 3.2\times 10^{-13}$ at 95$\%$ confidence. On the other hand, for such string tensions, constraints from the diffuse gamma-ray background (DGRB) indicate that more than 97$\%$ of the total network energy should be converted to dark matter (DM) or dark radiation. We also consider joint constraints on the annihilation cross-section, the mass, and the relic abundance of DM produced by decays of strings. For example, for a DM mass of 500 GeV, the observed relic abundance can be explained by decaying AH strings that also account for the NANOGrav signal. △ Less

Submitted 13 March, 2023; v1 submitted 12 October, 2022; originally announced October 2022.

Comments: 20 pages, 4 figures; v3: references are added, assumptions are clarified

Report number: RESCEU-15/22, HIP-2022-25/TH

Reconstructing physical parameters from template gravitational wave spectra at LISA: first order phase transitions

Authors: Chloe Gowling, Mark Hindmarsh, Deanna C. Hooper, Jesús Torrado

Abstract: A gravitational wave background from a first order phase transition in the early universe may be observable at millihertz gravitational wave (GW) detectors such as the Laser Interferometer Space Antenna (LISA). In this paper we introduce and test a method for investigating LISA's sensitivity to gravitational waves from a first order phase transition using parametrised templates as an approximation… ▽ More A gravitational wave background from a first order phase transition in the early universe may be observable at millihertz gravitational wave (GW) detectors such as the Laser Interferometer Space Antenna (LISA). In this paper we introduce and test a method for investigating LISA's sensitivity to gravitational waves from a first order phase transition using parametrised templates as an approximation to a more complete physical model. The motivation for developing the method is to provide a less computationally intensive way to perform Markov Chain Monte Carlo (MCMC) inference on the thermodynamic parameters of a first order phase transition, or on generally computationally intensive models. Starting from a map between the physical parameters and the parameters of an empirical template, we first construct a prior on the empirical parameters that contains the necessary information about the physical parameters; we then use the inverse mapping to reconstruct approximate posteriors on the physical parameters from a fast MCMC on the empirical template. We test the method on a double broken power law approximation to spectra in the sound shell model. The reconstruction method substantially reduces the proposal evaluation time, and despite requiring some precomputing of the mapping, this method is still cost-effective overall. In two test cases, with signal-to-noise $\sim 40$, the method recovers the physical parameters and the spectrum of the injected gravitational wave power spectrum to $95\%$ confidence. In previous Fisher matrix analysis we found the phase boundary speed $v_{\rm w}$ was expected to be the best constrained of the thermodynamic parameters. In this work, for an injected phase transition GW power spectrum with $v_{\rm w} = 0.55$, with a direct sample on the thermodynamic parameters we recover $0.630^{+0.17}_{-0.059}$ and for our reconstructed sample $0.646^{+0.098}_{-0.075}$. △ Less

Submitted 27 September, 2022; originally announced September 2022.

Prospects for LISA to detect a gravitational-wave background from first order phase transitions

Authors: Guillaume Boileau, Nelson Christensen, Chloe Gowling, Mark Hindmarsh, Renate Meyer

Abstract: First order phase transitions in the early universe could produce a gravitational-wave background that might be detectable by the Laser Interferometer Space Antenna (LISA). Such an observation would provide evidence for physics beyond the Standard Model. We study the ability of LISA to observe a gravitational-wave background from phase transitions in the presence of an extragalactic foreground fro… ▽ More First order phase transitions in the early universe could produce a gravitational-wave background that might be detectable by the Laser Interferometer Space Antenna (LISA). Such an observation would provide evidence for physics beyond the Standard Model. We study the ability of LISA to observe a gravitational-wave background from phase transitions in the presence of an extragalactic foreground from binary black hole mergers throughout the universe, a galactic foreground from white dwarf binaries, and LISA noise. Modelling the phase transition gravitational wave background as a double broken power law, we use the deviance information criterion as a detection statistic, and Fisher matrix and Markov Chain Monte Carlo methods to assess the measurement accuracy of the parameters of the power spectrum. While estimating all the parameters associated with the gravitational-wave backgrounds, foregrounds, and LISA noise, we find that LISA could detect a gravitational-wave background from phase transitions with a peak frequency of 1 mHz and normalized energy density amplitude of $Ω_{\text p} \simeq 3 \times 10^{-11}$. With $Ω_{\text p} \simeq 10^{-10}$, the signal is detectable if the peak frequency is in the range $4 \times 10^{-4}$ to $9 \times 10^{-3}$ Hz, and the peak amplitude and frequency can be estimated to an accuracy of 10\% to 1\%. △ Less

Submitted 28 February, 2023; v1 submitted 27 September, 2022; originally announced September 2022.

Journal ref: JCAP02(2023)056

Generation of gravitational waves from freely decaying turbulence

Authors: Pierre Auclair, Chiara Caprini, Daniel Cutting, Mark Hindmarsh, Kari Rummukainen, Danièle A. Steer, David J. Weir

Abstract: We study the stochastic gravitational wave background (SGWB) produced by freely decaying vortical turbulence in the early Universe. We thoroughly investigate the time correlation of the velocity field, and hence of the anisotropic stresses producing the gravitational waves. With hydrodynamical simulations, we show that the unequal time correlation function (UETC) of the Fourier components of the v… ▽ More We study the stochastic gravitational wave background (SGWB) produced by freely decaying vortical turbulence in the early Universe. We thoroughly investigate the time correlation of the velocity field, and hence of the anisotropic stresses producing the gravitational waves. With hydrodynamical simulations, we show that the unequal time correlation function (UETC) of the Fourier components of the velocity field is Gaussian in the time difference, as predicted by the "sweeping" decorrelation model. We introduce a decorrelation model that can be extended to wavelengths around the integral scale of the flow. Supplemented with the evolution laws of the kinetic energy and of the integral scale, this provides a new model UETC of the turbulent velocity field consistent with the simulations. We discuss the UETC as a positive definite kernel, and propose to use the Gibbs kernel for the velocity UETC as a natural way to ensure positive definiteness of the SGWB. The SGWB is given by a 4-dimensional integration of the resulting anisotropic stress UETC with the gravitational wave Green's function. We perform this integration using a Monte Carlo algorithm based on importance sampling, and find that the result matches that of the simulations. Furthermore, the SGWB obtained from the numerical integration and from the simulations show close agreement with a model in which the source is constant in time and abruptly turns off after a few eddy turnover times. Based on this assumption, we provide an approximate analytical form for the SGWB spectrum and its scaling with the initial kinetic energy and integral scale. Finally, we use our model and numerical integration algorithm to show that including an initial growth phase for the turbulent flow heavily influences the spectral shape of the SGWB. This highlights the importance of a complete understanding of the turbulence generation mechanism. △ Less

Submitted 12 September, 2022; v1 submitted 5 May, 2022; originally announced May 2022.

Comments: 54 pages, 30 figures

Report number: HIP-2021-35/TH

Journal ref: JCAP 09 (2022) 029

Cosmology with the Laser Interferometer Space Antenna

Authors: Pierre Auclair, David Bacon, Tessa Baker, Tiago Barreiro, Nicola Bartolo, Enis Belgacem, Nicola Bellomo, Ido Ben-Dayan, Daniele Bertacca, Marc Besancon, Jose J. Blanco-Pillado, Diego Blas, Guillaume Boileau, Gianluca Calcagni, Robert Caldwell, Chiara Caprini, Carmelita Carbone, Chia-Feng Chang, Hsin-Yu Chen, Nelson Christensen, Sebastien Clesse, Denis Comelli, Giuseppe Congedo, Carlo Contaldi, Marco Crisostomi , et al. (155 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) has two scientific objectives of cosmological focus: to probe the expansion rate of the universe, and to understand stochastic gravitational-wave backgrounds and their implications for early universe and particle physics, from the MeV to the Planck scale. However, the range of potential cosmological applications of gravitational wave observations exten… ▽ More The Laser Interferometer Space Antenna (LISA) has two scientific objectives of cosmological focus: to probe the expansion rate of the universe, and to understand stochastic gravitational-wave backgrounds and their implications for early universe and particle physics, from the MeV to the Planck scale. However, the range of potential cosmological applications of gravitational wave observations extends well beyond these two objectives. This publication presents a summary of the state of the art in LISA cosmology, theory and methods, and identifies new opportunities to use gravitational wave observations by LISA to probe the universe. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Report number: LISA CosWG-22-03

Gravitational Waves at Strong Coupling from an Effective Action

Authors: Fëanor Reuben Ares, Oscar Henriksson, Mark Hindmarsh, Carlos Hoyos, Niko Jokela

Abstract: Using a holographic derivation of a quantum effective action for a scalar operator at strong coupling, we compute quasi-equilibrium parameters relevant for the gravitational wave signal from a first order phase transition in a simple dual model. We discuss how the parameters of the phase transition vary with the effective number of degrees of freedom of the dual field theory. Our model can produce… ▽ More Using a holographic derivation of a quantum effective action for a scalar operator at strong coupling, we compute quasi-equilibrium parameters relevant for the gravitational wave signal from a first order phase transition in a simple dual model. We discuss how the parameters of the phase transition vary with the effective number of degrees of freedom of the dual field theory. Our model can produce an observable signal at LISA if the critical temperature is around a TeV, in a parameter region where the field theory has an approximate conformal symmetry. △ Less

Submitted 4 November, 2021; v1 submitted 27 October, 2021; originally announced October 2021.

Comments: 6 pages, 3 figures

Report number: HIP-2021-24/TH

Journal ref: Phys. Rev. Lett. 128, 131101 (2022)

Effective actions and bubble nucleation from holography

Authors: Fëanor Reuben Ares, Oscar Henriksson, Mark Hindmarsh, Carlos Hoyos, Niko Jokela

Abstract: We discuss the computation of the quantum effective action of strongly interacting field theories using holographic duality, and its use to determine quasi-equilibrium parameters of first order phase transitions relevant for gravitational wave production. A particularly simple holographic model is introduced, containing only the metric and a free massive scalar field. Despite the simplicity, the m… ▽ More We discuss the computation of the quantum effective action of strongly interacting field theories using holographic duality, and its use to determine quasi-equilibrium parameters of first order phase transitions relevant for gravitational wave production. A particularly simple holographic model is introduced, containing only the metric and a free massive scalar field. Despite the simplicity, the model contains a rich phase diagram, including first order phase transitions at non-zero temperature, due to various multi-trace deformations. We obtain the leading terms in the effective action from homogeneous black brane solutions in the gravity dual, and linearised perturbations around them. We then employ the effective action to construct bubble and domain wall solutions in the field theory side and study their properties. In particular, we show how the scaling of the effective action with the effective number of degrees of freedom of the quantum field theory determines the corresponding scaling of gravitational wave parameters. △ Less

Submitted 4 November, 2021; v1 submitted 28 September, 2021; originally announced September 2021.

Comments: 40 pages, 17 figures

Report number: HIP-2021-30/TH

Journal ref: Phys. Rev. D 105, 066020 (2022)

Comment on "More Axions from Strings''

Authors: Mark Hindmarsh, Joanes Lizarraga, Asier Lopez-Eiguren, Jon Urrestilla

Abstract: We comment on a claim that axion strings show a long-term logarithmic increase in the number of Hubble lengths per Hubble volume [arXiv:2007.04990], thereby violating the standard "scaling" expectation of an O(1) constant. We demonstrate that the string density data presented in [arXiv:2007.04990] are consistent with standard scaling, at a string density consistent with that obtained by us [arXiv:… ▽ More We comment on a claim that axion strings show a long-term logarithmic increase in the number of Hubble lengths per Hubble volume [arXiv:2007.04990], thereby violating the standard "scaling" expectation of an O(1) constant. We demonstrate that the string density data presented in [arXiv:2007.04990] are consistent with standard scaling, at a string density consistent with that obtained by us [arXiv:1908.03522, arXiv:2102.07723] and other groups. A transient slow growth in Hubble lengths per Hubble volume towards its constant scaling value is explained by standard network modelling [arXiv:2102.07723]. △ Less

Submitted 10 December, 2021; v1 submitted 20 September, 2021; originally announced September 2021.

Comments: 10 pages, 4 figures. Revised version: interpretation error corrected, improved discussion, global fit added

Observational prospects for phase transitions at LISA: Fisher matrix analysis

Authors: Chloe Gowling, Mark Hindmarsh

Abstract: A first order phase transition at the electroweak scale would lead to the production of gravitational waves that may be observable at upcoming space-based gravitational wave (GW) detectors such as LISA (Laser Interferometer Space Antenna). As the Standard Model has no phase transition, LISA can be used to search for new physics by searching for a stochastic gravitational wave background. In this w… ▽ More A first order phase transition at the electroweak scale would lead to the production of gravitational waves that may be observable at upcoming space-based gravitational wave (GW) detectors such as LISA (Laser Interferometer Space Antenna). As the Standard Model has no phase transition, LISA can be used to search for new physics by searching for a stochastic gravitational wave background. In this work we investigate LISA's sensitivity to the thermodynamic parameters encoded in the stochastic background produced by a phase transition, using the sound shell model to characterise the gravitational wave power spectrum, and the Fisher matrix to estimate uncertainties. We explore a parameter space with transition strengths $α< 0.5$ and phase boundary speeds $0.4 < v_\text{w} < 0.9$, for transitions nucleating at $T_{\text{N}} = 100$ GeV, with mean bubble spacings $0.1$ and $0.01$ of the Hubble length, and sound speed $c/\sqrt{3}$. We show that the power spectrum in the sound shell model can be well approximated by a four-parameter double broken power law, and find that the peak power and frequency can be measured to approximately 10% accuracy for signal-to-noise ratios (SNRs) above 20. Determinations of the underlying thermodynamic parameters are complicated by degeneracies, but in all cases the phase boundary speed will be the best constrained parameter. Turning to the principal components of the Fisher matrix, a signal-to-noise ratio above 20 produces a relative uncertainty less than 3% in the two highest-order principal components, indicating good prospects for combinations of parameters. The highest-order principal component is dominated by the wall speed. These estimates of parameter sensitivity provide a preliminary accuracy target for theoretical calculations of thermodynamic parameters. △ Less

Submitted 10 June, 2021; originally announced June 2021.

Comments: 37 pages, 11 figures

Loop decay in Abelian-Higgs string networks

Authors: Mark Hindmarsh, Joanes Lizarraga, Ander Urio, Jon Urrestilla

Abstract: We study the decay of cosmic string loops in the Abelian-Higgs model. We confirm earlier results that loops formed by intersections of infinite strings formed from random-field initial conditions disappear quickly, with lifetime proportional to their initial rest-frame length $\ell_\text{init}$. We study a population with $\ell_\text{init}$ up to $6000$ inverse mass units, and measure the proporti… ▽ More We study the decay of cosmic string loops in the Abelian-Higgs model. We confirm earlier results that loops formed by intersections of infinite strings formed from random-field initial conditions disappear quickly, with lifetime proportional to their initial rest-frame length $\ell_\text{init}$. We study a population with $\ell_\text{init}$ up to $6000$ inverse mass units, and measure the proportionality constant to be $0.14\pm0.04$, independently of the initial lengths. We propose a new method to construct oscillating non-self intersecting loops from initially stationary strings, and show that by contrast these loops have lifetimes scaling approximately as $\ell_\text{init}^2$, in line with previous works on artificially created string configurations. We show that the oscillating strings have mean-square velocity $\bar{v}^2 \simeq 0.500 \pm 0.004$, consistent with the Nambu-Goto value of $1/2$, while the network loops have $\bar{v}^2 \simeq 0.40 \pm 0.04$. We argue that whatever the mechanism behind the network loop decay is, it is non-linear, can only be suppressed by careful tuning of initial conditions, and is much stronger than gravitational radiation. An implication is that one cannot use the Nambu-Goto model to derive robust constraints on the tension of field theory strings. We advocate parametrising the uncertainty as the fraction $f_\text{NG}$ of Nambu-Goto-like loops surviving to radiate gravitationally. None of the 31 large network loops created survived longer than 0.25 of their initial length, so one can estimate that $f_\text{NG}<0.1$ at $95$% confidence level. If the recently reported NANOgrav signal is due to cosmic strings, $f_\text{NG}$ must be greater than $10^{-3}$ in order not to violate bounds from the Cosmic Microwave Background. △ Less

Submitted 30 March, 2021; originally announced March 2021.

Comments: 16 pages, 20 figures. Movies can be found as ancillary files

Report number: HIP-2021-15/TH

Journal ref: Phys. Rev. D 104, 043519 (2021)

Approach to scaling in axion string networks

Authors: Mark Hindmarsh, Joanes Lizarraga, Asier Lopez-Eiguren, Jon Urrestilla

Abstract: We study the approach to scaling in axion string networks in the radiation era, through measuring the root-mean-square velocity $v$ as well as the scaled mean string separation $x$. We find good evidence for a fixed point in the phase-space analysis in the variables $(x,v)$, providing a strong indication that standard scaling is taking place. We show that the approach to scaling can be well descri… ▽ More We study the approach to scaling in axion string networks in the radiation era, through measuring the root-mean-square velocity $v$ as well as the scaled mean string separation $x$. We find good evidence for a fixed point in the phase-space analysis in the variables $(x,v)$, providing a strong indication that standard scaling is taking place. We show that the approach to scaling can be well described by a two parameter velocity-one-scale (VOS) model, and show that the values of the parameters are insensitive to the initial state of the network. The string length has also been commonly expressed in terms of a dimensionless string length density $ζ$, proportional to the number of Hubble lengths of string per Hubble volume. In simulations with initial conditions far from the fixed point $ζ$ is still evolving after half a light-crossing time, which has been interpreted in the literature as a long-term logarithmic growth. We show that all our simulations, even those starting far from the fixed point, are accounted for by a VOS model with an asymptote of $ζ_*=1.20\pm0.09$ (calculated from the string length in the cosmic rest frame) and $v_* = 0.609\pm 0.014$. △ Less

Submitted 27 May, 2021; v1 submitted 15 February, 2021; originally announced February 2021.

Comments: 14 pages, 9 figures. v2: Minor changes, matches published version

Report number: HIP-2021-7/TH

Journal ref: Phys. Rev. D 103, 103534 (2021)

Gravitational waves from a holographic phase transition

Authors: Fëanor Reuben Ares, Mark Hindmarsh, Carlos Hoyos, Niko Jokela

Abstract: We investigate first order phase transitions in a holographic setting of five-dimensional Einstein gravity coupled to a scalar field, constructing phase diagrams of the dual field theory at finite temperature. We scan over the two-dimensional parameter space of a simple bottom-up model and map out important quantities for the phase transition: the region where first order phase transitions take pl… ▽ More We investigate first order phase transitions in a holographic setting of five-dimensional Einstein gravity coupled to a scalar field, constructing phase diagrams of the dual field theory at finite temperature. We scan over the two-dimensional parameter space of a simple bottom-up model and map out important quantities for the phase transition: the region where first order phase transitions take place; the latent heat, the transition strength parameter $α$, and the stiffness. We find that $α$ is generically in the range 0.1 to 0.3, and is strongly correlated with the stiffness (the square of the sound speed in a barotropic fluid). Using the LISA Cosmology Working Group gravitational wave power spectrum model corrected for kinetic energy suppression at large $α$ and non-conformal stiffness, we outline the observational prospects at the future space-based detectors LISA and TianQin. A TeV-scale hidden sector with a phase transition described by the model could be observable at both detectors. △ Less

Submitted 21 December, 2020; v1 submitted 25 November, 2020; originally announced November 2020.

Comments: 36 pages, 17 figures

Report number: HIP-2020-31/TH, Sussex-94886

Journal ref: J. High Energ. Phys. 2021, 100 (2021)

Challenges and Opportunities of Gravitational Wave Searches at MHz to GHz Frequencies

Authors: N. Aggarwal, O. D. Aguiar, A. Bauswein, G. Cella, S. Clesse, A. M. Cruise, V. Domcke, D. G. Figueroa, A. Geraci, M. Goryachev, H. Grote, M. Hindmarsh, F. Muia, N. Mukund, D. Ottaway, M. Peloso, F. Quevedo, A. Ricciardone, J. Steinlechner, S. Steinlechner, S. Sun, M. E. Tobar, F. Torrenti, C. Unal, G. White

Abstract: The first direct measurement of gravitational waves by the LIGO and Virgo collaborations has opened up new avenues to explore our Universe. This white paper outlines the challenges and gains expected in gravitational wave searches at frequencies above the LIGO/Virgo band, with a particular focus on Ultra High-Frequency Gravitational Waves (UHF-GWs), covering the MHz to GHz range. The absence of kn… ▽ More The first direct measurement of gravitational waves by the LIGO and Virgo collaborations has opened up new avenues to explore our Universe. This white paper outlines the challenges and gains expected in gravitational wave searches at frequencies above the LIGO/Virgo band, with a particular focus on Ultra High-Frequency Gravitational Waves (UHF-GWs), covering the MHz to GHz range. The absence of known astrophysical sources in this frequency range provides a unique opportunity to discover physics beyond the Standard Model operating both in the early and late Universe, and we highlight some of the most promising gravitational sources. We review several detector concepts which have been proposed to take up this challenge, and compare their expected sensitivity with the signal strength predicted in various models. This report is the summary of the workshop "Challenges and opportunities of high-frequency gravitational wave detection" held at ICTP Trieste, Italy in October 2019, that set up the stage for the recently launched Ultra-High-Frequency Gravitational Wave (UHF-GW) initiative. △ Less

Submitted 13 December, 2021; v1 submitted 24 November, 2020; originally announced November 2020.

Comments: Published in Living Reviews in Relativity

Report number: CERN-TH-2020-185, HIP-2020-28/TH, DESY 20-195

Journal ref: Living Reviews in Relativity volume 24, Article number: 4 (2021)

Phase transitions in the early universe

Authors: Mark B. Hindmarsh, Marvin Lüben, Johannes Lumma, Martin Pauly

Abstract: These lecture notes are based on a course given by Mark Hindmarsh at the 24th Saalburg Summer School 2018 and written up by Marvin Lüben, Johannes Lumma and Martin Pauly. The aim is to provide the necessary basics to understand first-order phase transitions in the early universe, to outline how they leave imprints in gravitational waves, and advertise how those gravitational waves could be detecte… ▽ More These lecture notes are based on a course given by Mark Hindmarsh at the 24th Saalburg Summer School 2018 and written up by Marvin Lüben, Johannes Lumma and Martin Pauly. The aim is to provide the necessary basics to understand first-order phase transitions in the early universe, to outline how they leave imprints in gravitational waves, and advertise how those gravitational waves could be detected in the future. A first-order phase transition at the electroweak scale is a prediction of many theories beyond the Standard Model, and is also motivated as an ingredient of some theories attempting to provide an explanation for the matter-antimatter asymmetry in our Universe. Starting from bosonic and fermionic statistics, we derive Boltzmann's equation and generalise to a fluid of particles with field dependent mass. We introduce the thermal effective potential for the field in its lowest order approximation, discuss the transition to the Higgs phase in the Standard Model and beyond, and compute the probability for the field to cross a potential barrier. After these preliminaries, we provide a hydrodynamical description of first-order phase transitions as it is appropriate for describing the early Universe. We thereby discuss the key quantities characterising a phase transition, and how they are imprinted in the gravitational wave power spectrum that might be detectable by the space-based gravitational wave detector LISA in the 2030s. △ Less

Submitted 25 May, 2021; v1 submitted 20 August, 2020; originally announced August 2020.

Comments: Introduction to the topic of phase transitions in the early universe, 63 pages

Report number: MPP-2020-163, HIP-2020-27/TH

Journal ref: SciPost Phys. Lect. Notes 24 (2021)

Irreducible background of gravitational waves from a cosmic defect network: update and comparison of numerical techniques

Authors: Daniel G. Figueroa, Mark Hindmarsh, Joanes Lizarraga, Jon Urrestilla

Abstract: Cosmological phase transitions in the early Universe may produce relics in the form of a network of cosmic defects. Independently of the order of a phase transition, topology of the defects, and their global or gauge nature, the defects are expected to emit gravitational waves (GWs) as the network energy-momentum tensor adapts itself to maintaining {scaling}. We show that the evolution of any defe… ▽ More Cosmological phase transitions in the early Universe may produce relics in the form of a network of cosmic defects. Independently of the order of a phase transition, topology of the defects, and their global or gauge nature, the defects are expected to emit gravitational waves (GWs) as the network energy-momentum tensor adapts itself to maintaining {scaling}. We show that the evolution of any defect network (and for that matter any scaling source) emits a GW background with spectrum $Ω_{\rm GW} \propto f^3$ for $f \ll f_0$, $Ω_{\rm GW} \propto 1/f^2$ for $f_0 \lesssim f \lesssim f_{\rm eq}$, and $Ω_{\rm GW} \propto~const$ (i.e.~exactly scale-invariant) for $f \gg f_{\rm eq}$, where $f_0$ and $ f_{\rm eq}$ denote respectively the frequencies corresponding to the present and matter-radiation equality horizons. This background represents an irreducible emission of GWs from any scaling network of cosmic defects, with its amplitude characterized only by the symmetry breaking scale and the nature of the defects. Using classical lattice simulations we calculate the GW signal emitted by defects created after the breaking of a global symmetry $O(N) \rightarrow O(N-1)$. We obtain the GW spectrum for $N$ between 2 and 20 with two different techniques: integrating over unequal time correlators of the energy momentum tensor, updating our previous work on smaller lattices, and for the first time, comparing the result with the real time evolution of the tensor perturbations sourced by the same defects. Our results validate the equivalence of the two techniques. Using CMB upper bounds on the defects' energy scale, we discuss the difficulty of detecting this GW background in the case of global defects. △ Less

Submitted 7 July, 2020; originally announced July 2020.

Comments: 24 pages, 8 figures

Journal ref: Phys. Rev. D 102, 103516 (2020)

Gravitational waves from vacuum first order phase transitions II: from thin to thick walls

Authors: Daniel Cutting, Elba Granados Escartin, Mark Hindmarsh, David J. Weir

Abstract: In a vacuum first-order phase transition, gravitational waves are generated from collision of bubbles of the true vacuum. The spectrum from such collisions takes the form of a broken power law. We consider a toy model for such a phase transition, where the dynamics of the scalar field depends on a single parameter $\overlineλ$, which controls how thin the bubble wall is at nucleation and how close… ▽ More In a vacuum first-order phase transition, gravitational waves are generated from collision of bubbles of the true vacuum. The spectrum from such collisions takes the form of a broken power law. We consider a toy model for such a phase transition, where the dynamics of the scalar field depends on a single parameter $\overlineλ$, which controls how thin the bubble wall is at nucleation and how close to degenerate the vacua are relative to the barrier. We extend on our previous work by performing a series of simulations with a range of $\overlineλ$. The peak of the gravitational-wave power spectrum varies by up to a factor of $1.3$, which is probably an unobservable effect. We find that the ultraviolet (UV) power law in the gravitational-wave spectrum becomes steeper as $\overlineλ \rightarrow 0$, varying between $k^{-1.4}$ and $k^{-2.2}$ for the $\overlineλ$ considered. This provides some evidence that the form of the underlying effective potential of a vacuum first-order phase transition could be determined from the gravitational-wave spectrum it produces. △ Less

Submitted 7 January, 2021; v1 submitted 27 May, 2020; originally announced May 2020.

Comments: 24 pages, 20 figures. Updated to match version accepted for publication. Modification of how lattice uncertainty was incorporated into analysis leads to a moderate quantitative but not qualitative change in the fits reported

Report number: HIP-2020-13/TH

Journal ref: Phys. Rev. D 103, 023531 (2021)

Detecting gravitational waves from cosmological phase transitions with LISA: an update

Authors: Chiara Caprini, Mikael Chala, Glauber C. Dorsch, Mark Hindmarsh, Stephan J. Huber, Thomas Konstandin, Jonathan Kozaczuk, Germano Nardini, Jose Miguel No, Kari Rummukainen, Pedro Schwaller, Geraldine Servant, Anders Tranberg, David J. Weir

Abstract: We investigate the potential for observing gravitational waves from cosmological phase transitions with LISA in light of recent theoretical and experimental developments. Our analysis is based on current state-of-the-art simulations of sound waves in the cosmic fluid after the phase transition completes. We discuss the various sources of gravitational radiation, the underlying parameters describin… ▽ More We investigate the potential for observing gravitational waves from cosmological phase transitions with LISA in light of recent theoretical and experimental developments. Our analysis is based on current state-of-the-art simulations of sound waves in the cosmic fluid after the phase transition completes. We discuss the various sources of gravitational radiation, the underlying parameters describing the phase transition and a variety of viable particle physics models in this context, clarifying common misconceptions that appear in the literature and identifying open questions requiring future study. We also present a web-based tool, PTPlot, that allows users to obtain up-to-date detection prospects for a given set of phase transition parameters at LISA. △ Less

Submitted 12 January, 2021; v1 submitted 29 October, 2019; originally announced October 2019.

Comments: 68 pages, 16 figures, some comments added, typos fixed

Gravitational waves from first order cosmological phase transitions in the Sound Shell Model

Authors: Mark Hindmarsh, Mulham Hijazi

Abstract: We calculate gravitational wave power spectra from first order early Universe phase transitions using the Sound Shell Model. The model predicts that the power spectrum depends on the mean bubble separation, the phase transition strength, the phase boundary speed, with the overall frequency scale set by the nucleation temperature. There is also a dependence on the time evolution of the bubble nucle… ▽ More We calculate gravitational wave power spectra from first order early Universe phase transitions using the Sound Shell Model. The model predicts that the power spectrum depends on the mean bubble separation, the phase transition strength, the phase boundary speed, with the overall frequency scale set by the nucleation temperature. There is also a dependence on the time evolution of the bubble nucleation rate. The gravitational wave peak power and frequency are in good agreement with published numerical simulations, where bubbles are nucleated simultaneously. Agreement is particularly good for detonations, but the total power for deflagrations is predicted higher than numerical simulations show, indicating refinement of the model of the transfer of energy to the fluid is needed for accurate computations. We show how the time-dependence of the bubble nucleation rate affects the shape of the power spectrum: an exponentially rising nucleation rate produces higher amplitude gravitational waves at a longer wavelength than simultaneous nucleation. We present an improved fit for the predicted gravitational wave power spectrum in the form of a double broken power law, where the two breaks in the slope happen at wavenumber corresponding to the mean bubble separation and the thickness of the fluid shell surrounding the expanding bubbles, which in turn is related to the difference of the phase boundary speed from the speed of sound. △ Less

Submitted 22 September, 2019; originally announced September 2019.

Comments: 44 pp

Report number: NORDITA-2019-083, HIP-2019-29/TH

The scaling density of axion strings

Authors: Mark Hindmarsh, Joanes Lizarraga, Asier Lopez-Eiguren, Jon Urrestilla

Abstract: In the QCD axion dark matter scenario with post-inflationary Peccei-Quinn symmetry breaking, the number density of axions, and hence the dark matter density, depends on the length of string per unit volume at cosmic time $t$, by convention written $ζ/t^2$. The expectation has been that the dimensionless parameter $ζ$ tends to a constant $ζ_0$, a feature of a string network known as scaling. It has… ▽ More In the QCD axion dark matter scenario with post-inflationary Peccei-Quinn symmetry breaking, the number density of axions, and hence the dark matter density, depends on the length of string per unit volume at cosmic time $t$, by convention written $ζ/t^2$. The expectation has been that the dimensionless parameter $ζ$ tends to a constant $ζ_0$, a feature of a string network known as scaling. It has recently been claimed that in larger numerical simulations $ζ$ shows a logarithmic increase with time, while theoretical modelling suggests an inverse logarithmic correction. Either case would result in a large enhancement of the string density at the QCD transition, and a substantial revision to the axion mass required for the axion to constitute all of the dark matter. With a set of new simulations of global strings we compare the standard scaling (constant-$ζ$) model to the logarithmic growth and inverse-logarithmic correction models. In the standard scaling model, by fitting to linear growth in the mean string separation $ξ= t/\sqrtζ$, we find $ζ_0 = 1.19 \pm 0.20$. We conclude that the apparent corrections to $ζ$ are artefacts of the initial conditions, rather than a property of the scaling network. The residuals from the constant-$ζ$ (linear $ξ$) fit also show no evidence for logarithmic growth, restoring confidence that numerical simulations can be simply extrapolated from the Peccei-Quinn symmetry-breaking scale to the QCD scale. Re-analysis of previous work on the axion number density suggests that recent estimates of the axion dark matter mass in the post-inflationary symmetry-breaking scenario we study should be increased by about 50%. △ Less

Submitted 5 February, 2020; v1 submitted 9 August, 2019; originally announced August 2019.

Comments: 8 pages including supplemental material; 5 figures. v2 matches published version

Report number: HIP-2020-1/TH

Journal ref: Phys. Rev. Lett. 124, 021301 (2020)

Vorticity, kinetic energy, and suppressed gravitational wave production in strong first order phase transitions

Authors: Daniel Cutting, Mark Hindmarsh, David J. Weir

Abstract: We have performed the first 3-dimensional simulations of strong first-order thermal phase transitions in the early Universe. For deflagrations, we find that the rotational component of the fluid velocity increases as the transition strength is increased. For detonations, however, the rotational velocity component remains constant and small. We also find that the efficiency with which kinetic energ… ▽ More We have performed the first 3-dimensional simulations of strong first-order thermal phase transitions in the early Universe. For deflagrations, we find that the rotational component of the fluid velocity increases as the transition strength is increased. For detonations, however, the rotational velocity component remains constant and small. We also find that the efficiency with which kinetic energy is transferred to the fluid falls below theoretical expectations as we increase the transition strength. The probable origin of the kinetic energy deficit is the formation of reheated droplets of the metastable phase during the collision, slowing the bubble walls. The rate of increase in the gravitational wave energy density for deflagrations in strong transitions is suppressed compared to that predicted in earlier work. This is largely accounted for by the reduction in kinetic energy. Current modelling therefore substantially overestimates the gravitational wave signal for strong transitions with deflagrations, in the most extreme case by a factor of $10^{3}$. Detonations are less affected. △ Less

Submitted 19 July, 2020; v1 submitted 2 June, 2019; originally announced June 2019.

Comments: 6+6 pages, 4+8 figures. Movies illustrating the physics in this paper available at https://vimeo.com/showcase/5968055 Updated to reflect published version

Report number: HIP-2019-15/TH

Journal ref: Phys. Rev. Lett. 125, 021302 (2020)

Type I Abelian Higgs strings: evolution and Cosmic Microwave Background constraints

Authors: Mark Hindmarsh, Joanes Lizarraga, Jon Urrestilla, David Daverio, Martin Kunz

Abstract: We present results from the first simulations of networks of Type I Abelian Higgs cosmic strings to include both matter and radiation eras and Cosmic Microwave Background (CMB) constraints. In Type I strings, the string tension is a slowly decreasing function of the ratio of the scalar and gauge mass-squared, $β$. We find that the mean string separation shows no dependence on $β$, and that the ene… ▽ More We present results from the first simulations of networks of Type I Abelian Higgs cosmic strings to include both matter and radiation eras and Cosmic Microwave Background (CMB) constraints. In Type I strings, the string tension is a slowly decreasing function of the ratio of the scalar and gauge mass-squared, $β$. We find that the mean string separation shows no dependence on $β$, and that the energy-momentum tensor correlators decrease approximately in proportion to the square of the string tension, with additional O(1) correction factors which asymptote to constants below $β\lesssim 0.01$. Strings in models with low self-couplings can therefore satisfy current CMB bounds at higher symmetry-breaking scales. This is particularly relevant for models where the gauge symmetry is broken in a supersymmetric flat direction, for which the effective self-coupling can be extremely small. If our results can be extrapolated to $β\simeq 10^{-15}$, even strings formed at $10^{16}$ GeV (approximately the grand unification scale in supersymmetric extensions of the Standard Model) can be compatible with CMB constraints. △ Less

Submitted 29 April, 2019; v1 submitted 20 December, 2018; originally announced December 2018.

Comments: 11 pages, 9 figures. v2 Discussion of previous Type I cosmic string simulations added (Ref. [36]), main results unchanged. v3 matches published version

Journal ref: Phys. Rev. D 99, 083522 (2019)

Scaling in necklaces of monopoles and semipoles

Authors: Mark Hindmarsh, Anna Kormu, Asier Lopez-Eiguren, David J. Weir

Abstract: Models of symmetry breaking in the early universe can produce networks of cosmic strings threading 't Hooft-Polyakov monopoles. In certain cases there is a larger global symmetry group and the monopoles split into so-called semipoles. These networks are all known as cosmic necklaces. We carry out large-scale field theory simulations of the simplest model containing these objects, confirming that t… ▽ More Models of symmetry breaking in the early universe can produce networks of cosmic strings threading 't Hooft-Polyakov monopoles. In certain cases there is a larger global symmetry group and the monopoles split into so-called semipoles. These networks are all known as cosmic necklaces. We carry out large-scale field theory simulations of the simplest model containing these objects, confirming that the energy density of networks of cosmic necklaces approaches scaling, i.e. that it remains a constant fraction of the background energy density. The number of monopoles per unit comoving string length is constant, meaning that the density fraction of monopoles decreases with time. Where the necklaces carry semipoles rather than monopoles, we perform the first simulations large enough to demonstrate that they also maintain a constant number per unit comoving string length. We also compare our results to a number of analytical models of cosmic necklaces, finding that none explains our results. We put forward evidence that annihilation of poles on the strings is controlled by a diffusive process, a possibility not considered before. The observational constraints derived in our previous work for necklaces with monopoles can now be safely applied to those with semipoles as well. △ Less

Submitted 12 December, 2018; v1 submitted 10 September, 2018; originally announced September 2018.

Comments: 12 pages and 7 figures, updated to match with published version

Report number: HIP-2018-11/TH

Journal ref: Phys. Rev. D 98, 103533 (2018)

Gravitational waves from vacuum first-order phase transitions: from the envelope to the lattice

Authors: Daniel Cutting, Mark Hindmarsh, David J. Weir

Abstract: We conduct large scale numerical simulations of gravitational wave production at a first order vacuum phase transition. We find a power law for the gravitational wave power spectrum at high wavenumber which falls off as $k^{-1.5}$ rather than the $k^{-1}$ produced by the envelope approximation. The peak of the power spectrum is shifted to slightly lower wave numbers from that of the envelope appro… ▽ More We conduct large scale numerical simulations of gravitational wave production at a first order vacuum phase transition. We find a power law for the gravitational wave power spectrum at high wavenumber which falls off as $k^{-1.5}$ rather than the $k^{-1}$ produced by the envelope approximation. The peak of the power spectrum is shifted to slightly lower wave numbers from that of the envelope approximation. The envelope approximation reproduces our results for the peak power less well, agreeing only to within an order of magnitude. After the bubbles finish colliding the scalar field oscillates around the true vacuum. An additional feature is produced in the UV of the gravitational wave power spectrum, and this continues to grow linearly until the end of our simulation. The additional feature peaks at a length scale close to the bubble wall thickness and is shown to have a negligible contribution to the energy in gravitational waves, providing the scalar field mass is much smaller than the Planck mass. △ Less

Submitted 13 June, 2018; v1 submitted 15 February, 2018; originally announced February 2018.

Comments: 15 pages, 15 figures, minor changes, updated to match with published version

Report number: HIP-2018-4/TH

Journal ref: Phys. Rev. D 97, 123513 (2018)

Primordial black holes with an accurate QCD equation of state

Authors: Christian T. Byrnes, Mark Hindmarsh, Sam Young, Michael R. S. Hawkins

Abstract: Making use of definitive new lattice computations of the Standard Model thermodynamics during the quantum chromodynamic (QCD) phase transition, we calculate the enhancement in the mass distribution of primordial black holes (PBHs) due to the softening of the equation of state. We find that the enhancement peaks at approximately $0.7M_\odot$, with the formation rate increasing by at least two order… ▽ More Making use of definitive new lattice computations of the Standard Model thermodynamics during the quantum chromodynamic (QCD) phase transition, we calculate the enhancement in the mass distribution of primordial black holes (PBHs) due to the softening of the equation of state. We find that the enhancement peaks at approximately $0.7M_\odot$, with the formation rate increasing by at least two orders of magnitude due to the softening of the equation of state at this time, with a range of approximately $0.3M_\odot<M<1.4M_\odot$ at full width half-maximum. PBH formation is increased by a smaller amount for PBHs with masses spanning a large range, $10^{-3}M_\odot<M_{\rm PBH}<10^{3}M_\odot$, which includes the masses of the BHs that LIGO detected. The most significant source of uncertainty in the number of PBHs formed is now due to unknowns in the formation process, rather than from the phase transition. A near scale-invariant density power spectrum tuned to generate a population with mass and merger rate consistent with that detected by LIGO should also produce a much larger energy density of PBHs with solar mass. The existence of BHs below the Chandresekhar mass limit would be a smoking gun for a primordial origin and they could arguably constitute a significant fraction of the cold dark matter density. They also pose a challenge to inflationary model building which seek to produce the LIGO BHs without overproducing lighter PBHs. △ Less

Submitted 21 September, 2018; v1 submitted 18 January, 2018; originally announced January 2018.

Comments: 18 pages, 5 figures. v2: Improvements in the presentation and references added, matches the published version

Exploring Cosmic Origins with CORE: Survey requirements and mission design

Authors: J. Delabrouille, P. de Bernardis, F. R. Bouchet, A. Achúcarro, P. A. R. Ade, R. Allison, F. Arroja, E. Artal, M. Ashdown, C. Baccigalupi, M. Ballardini, A. J. Banday, R. Banerji, D. Barbosa, J. Bartlett, N. Bartolo, S. Basak, J. J. A. Baselmans, K. Basu, E. S. Battistelli, R. Battye, D. Baumann, A. Benoît, M. Bersanelli, A. Bideaud , et al. (178 additional authors not shown)

Abstract: Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium… ▽ More Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission. △ Less

Submitted 14 June, 2017; originally announced June 2017.

Comments: 79 pages, 14 figures

Cosmic Microwave Background constraints for global strings and global monopoles

Authors: Asier Lopez-Eiguren, Joanes Lizarraga, Mark Hindmarsh, Jon Urrestilla

Abstract: We present the first CMB power spectra from numerical simulations of the global O(N) linear $σ$-model with N = 2,3, which have global strings and monopoles as topological defects. In order to compute the CMB power spectra we compute the unequal time correlators (UETCs) of the energy-momentum tensor, showing that they fall off at high wave number faster than naive estimates based on the geometry of… ▽ More We present the first CMB power spectra from numerical simulations of the global O(N) linear $σ$-model with N = 2,3, which have global strings and monopoles as topological defects. In order to compute the CMB power spectra we compute the unequal time correlators (UETCs) of the energy-momentum tensor, showing that they fall off at high wave number faster than naive estimates based on the geometry of the defects, indicating non-trivial (anti-)correlations between the defects and the surrounding Goldstone boson field. We obtain source functions for Einstein-Boltzmann solvers from the UETCs, using a recent method that improves the modelling at the radiation- matter transition. We show that the interpolation function that mimics the transition is similar to other defect models, but not identical, confirming the non-universality of the interpolation function. The CMB power spectra for global strings and monopoles have the same overall shape as those obtained using the non-linear $σ$-model approximation, which is well captured by a large-N calculation. However, the amplitudes are larger than the large-N calculation predict, and in the case of global strings much larger: a factor of 20 at the peak. Finally we compare the CMB power spectra with the latest CMB data to put limits on the allowed contribution to the temperature power spectrum at multipole $\ell$ = 10 of 1.7% for global strings and 2.4% for global monopoles. These limits correspond to symmetry-breaking scales of 2.9x1015 GeV (6.3x1014 GeV with the expected logarithmic scaling of the effective string tension between the simulation time and decoupling) and 6.4x1015 GeV respectively. The bound on global strings is a significant one for the ultra-light axion scenario with axion masses ma 10-28 eV. These upper limits indicate that gravitational wave from global topological defects will not be observable at the GW observatory LISA. △ Less

Submitted 27 July, 2017; v1 submitted 11 May, 2017; originally announced May 2017.

Comments: Minor changes, matches published version; 24 pages, 9 figures

Report number: HIP-2017-07/TH

Journal ref: JCAP 1707 (2017) no.07, 026

Shape of the acoustic gravitational wave power spectrum from a first order phase transition

Authors: Mark Hindmarsh, Stephan J. Huber, Kari Rummukainen, David J. Weir

Abstract: We present results from large-scale numerical simulations of a first order thermal phase transition in the early universe, in order to explore the shape of the acoustic gravitational wave and the velocity power spectra. We compare the results with the predictions of the recently proposed sound shell model. For the gravitational wave power spectrum, we find that the predicted $k^{-3}$ behaviour, wh… ▽ More We present results from large-scale numerical simulations of a first order thermal phase transition in the early universe, in order to explore the shape of the acoustic gravitational wave and the velocity power spectra. We compare the results with the predictions of the recently proposed sound shell model. For the gravitational wave power spectrum, we find that the predicted $k^{-3}$ behaviour, where $k$ is the wavenumber, emerges clearly for detonations. The power spectra from deflagrations show similar features, but exhibit a steeper high-$k$ decay and an extra feature not accounted for in the model. There are two independent length scales: the mean bubble separation and the thickness of the sound shell around the expanding bubble of the low temperature phase. It is the sound shell thickness which sets the position of the peak of the power spectrum. The low wavenumber behaviour of the velocity power spectrum is consistent with a causal $k^{3}$, except for the thinnest sound shell, where it is steeper. We present parameters for a simple broken power law fit to the gravitational wave power spectrum for wall speeds well away from the speed of sound where this form can be usefully applied. We examine the prospects for the detection, showing that a LISA-like mission has the sensitivity to detect a gravitational wave signal from sound waves with an RMS fluid velocity of about $0.05c$, produced from bubbles with a mean separation of about $10^{-2}$ of the Hubble radius. The shape of the gravitational wave power spectrum depends on the bubble wall speed, and it may be possible to estimate the wall speed, and constrain other phase transition parameters, with an accurate measurement of a stochastic gravitational wave background. △ Less

Submitted 15 April, 2020; v1 submitted 19 April, 2017; originally announced April 2017.

Comments: 16+1 pages, 9+1 figures; v2: added a few paragraphs about the importance of turbulence, published in PRD; v3: added erratum correcting two equations and the SNR plot, methods and conclusions unchanged

Report number: HIP-2017-02/TH

Journal ref: Phys. Rev. D 96, 103520 (2017); Erratum: Phys. Rev. D 101, 089902 (2020)

Scaling from gauge and scalar radiation in Abelian Higgs string networks

Authors: Mark Hindmarsh, Joanes Lizarraga, Jon Urrestilla, David Daverio, Martin Kunz

Abstract: We investigate cosmic string networks in the Abelian Higgs model using data from a campaign of large-scale numerical simulations on lattices of up to $4096^3$ grid points. We observe scaling or self-similarity of the networks over a wide range of scales, and estimate the asymptotic values of the mean string separation in horizon length units $\dotξ$ and of the mean square string velocity… ▽ More We investigate cosmic string networks in the Abelian Higgs model using data from a campaign of large-scale numerical simulations on lattices of up to $4096^3$ grid points. We observe scaling or self-similarity of the networks over a wide range of scales, and estimate the asymptotic values of the mean string separation in horizon length units $\dotξ$ and of the mean square string velocity $\bar v^2$ in the continuum and large time limits. The scaling occurs because the strings lose energy into classical radiation of the scalar and gauge fields of the Abelian Higgs model. We quantify the energy loss with a dimensionless radiative efficiency parameter, and show that it does not vary significantly with lattice spacing or string separation. This implies that the radiative energy loss underlying the scaling behaviour is not a lattice artefact, and justifies the extrapolation of measured network properties to large times for computations of cosmological perturbations. We also show that the core growth method, which increases the defect core width with time to extend the dynamic range of simulations, does not introduce significant systematic error. We compare $\dotξ$ and $\bar v^2$ to values measured in simulations using the Nambu-Goto approximation, finding that the latter underestimate the mean string separation by about 25%, and overestimate $\bar v^2$ by about 10%. The scaling of the string separation implies that string loops decay by the emission of massive radiation within a Hubble time in field theory simulations, in contrast to the Nambu-Goto scenario which neglects this energy loss mechanism. String loops surviving for only one Hubble time emit much less gravitational radiation than in the Nambu-Goto scenario, and are consequently subject to much weaker gravitational wave constraints on their tension. △ Less

Submitted 21 August, 2017; v1 submitted 20 March, 2017; originally announced March 2017.

Comments: 19 pages, 12 figures. v2: typos fixed and video link added https://vimeo.com/214671412 . v3: minor changes, matches published version

Journal ref: Phys. Rev. D 96, 023525 (2017)

Exploring Cosmic Origins with CORE: Inflation

Authors: CORE Collaboration, Fabio Finelli, Martin Bucher, Ana Achúcarro, Mario Ballardini, Nicola Bartolo, Daniel Baumann, Sébastien Clesse, Josquin Errard, Will Handley, Mark Hindmarsh, Kimmo Kiiveri, Martin Kunz, Anthony Lasenby, Michele Liguori, Daniela Paoletti, Christophe Ringeval, Jussi Väliviita, Bartjan van Tent, Vincent Vennin, Rupert Allison, Frederico Arroja, Marc Ashdown, A. J. Banday, Ranajoy Banerji , et al. (107 additional authors not shown)

Abstract: We forecast the scientific capabilities to improve our understanding of cosmic inflation of CORE, a proposed CMB space satellite submitted in response to the ESA fifth call for a medium-size mission opportunity. The CORE satellite will map the CMB anisotropies in temperature and polarization in 19 frequency channels spanning the range 60-600 GHz. CORE will have an aggregate noise sensitivity of… ▽ More We forecast the scientific capabilities to improve our understanding of cosmic inflation of CORE, a proposed CMB space satellite submitted in response to the ESA fifth call for a medium-size mission opportunity. The CORE satellite will map the CMB anisotropies in temperature and polarization in 19 frequency channels spanning the range 60-600 GHz. CORE will have an aggregate noise sensitivity of $1.7 μ$K$\cdot \,$arcmin and an angular resolution of 5' at 200 GHz. We explore the impact of telescope size and noise sensitivity on the inflation science return by making forecasts for several instrumental configurations. This study assumes that the lower and higher frequency channels suffice to remove foreground contaminations and complements other related studies of component separation and systematic effects, which will be reported in other papers of the series "Exploring Cosmic Origins with CORE." We forecast the capability to determine key inflationary parameters, to lower the detection limit for the tensor-to-scalar ratio down to the $10^{-3}$ level, to chart the landscape of single field slow-roll inflationary models, to constrain the epoch of reheating, thus connecting inflation to the standard radiation-matter dominated Big Bang era, to reconstruct the primordial power spectrum, to constrain the contribution from isocurvature perturbations to the $10^{-3}$ level, to improve constraints on the cosmic string tension to a level below the presumptive GUT scale, and to improve the current measurements of primordial non-Gaussianities down to the $f_{NL}^{\rm local} < 1$ level. For all the models explored, CORE alone will improve significantly on the present constraints on the physics of inflation. Its capabilities will be further enhanced by combining with complementary future cosmological observations. △ Less

Submitted 5 April, 2017; v1 submitted 25 December, 2016; originally announced December 2016.

Comments: Latex 107 pages, revised with updated author list and minor modifications

Breaking symmetry, breaking ground

Authors: Mark Hindmarsh

Abstract: I give a short commentary on a seminal article by T W B Kibble in 1976, "Topology of cosmic domains and strings". I give a short commentary on a seminal article by T W B Kibble in 1976, "Topology of cosmic domains and strings". △ Less

Submitted 23 November, 2016; originally announced December 2016.

Comments: Published in J. Phys. A as part of a series of Viewpoints on 50 influential papers, celebrating the 50th anniversary of the Journal of Physics series. The article is dedicated to Tom, who died while it was being prepared

Journal ref: J. Phys. A: Math. Theor. 49 (2016) 411001

Numerical simulations of necklaces in SU(2) gauge-Higgs field theory

Authors: Mark Hindmarsh, Kari Rummukainen, David J. Weir

Abstract: We perform the first numerical simulations of necklaces in a non-Abelian gauge theory. Necklaces are composite classical solutions which can be interpreted as monopoles trapped on strings, rather generic structures in a Grand Unified Theory. We generate necklaces from random initial conditions, modelling a phase transition in the early Universe, and study the evolution. For all cases, we find that… ▽ More We perform the first numerical simulations of necklaces in a non-Abelian gauge theory. Necklaces are composite classical solutions which can be interpreted as monopoles trapped on strings, rather generic structures in a Grand Unified Theory. We generate necklaces from random initial conditions, modelling a phase transition in the early Universe, and study the evolution. For all cases, we find that the necklace system shows scaling behaviour similar to that of a network of ordinary cosmic strings. Furthermore, our simulations indicate that comoving distance between the monopoles or semipoles along the string asymptotes to a constant value at late times. This means that while the monopole-to-string energy density ratio decreases as the inverse of the scale factor, a horizon-size length of string has a large number of monopoles, significantly affecting the dynamics of string loops. We argue that gravitational wave bounds from millisecond pulsar timing on the string tension in the Nambu-Goto scenario are greatly relaxed. △ Less

Submitted 31 March, 2017; v1 submitted 25 November, 2016; originally announced November 2016.

Comments: 15 pages, 8 figures; v2: Figs. 2 (bottom) and 4 (top) corrected, some discussion added, results unchanged, accepted by PRD

Report number: HIP-2016-28/TH

Journal ref: Phys. Rev. D 95, 063520 (2017)

New CMB constraints for Abelian Higgs cosmic strings

Authors: Joanes Lizarraga, Jon Urrestilla, David Daverio, Mark Hindmarsh, Martin Kunz

Abstract: We present cosmic microwave background (CMB) power spectra from recent numerical simulations of cosmic strings in the Abelian Higgs model and compare them to CMB power spectra measured by Planck. We obtain revised constraints on the cosmic string tension parameter $Gμ$. For example, in the $Λ$CDM model with the addition of strings and no primordial tensor perturbations, we find… ▽ More We present cosmic microwave background (CMB) power spectra from recent numerical simulations of cosmic strings in the Abelian Higgs model and compare them to CMB power spectra measured by Planck. We obtain revised constraints on the cosmic string tension parameter $Gμ$. For example, in the $Λ$CDM model with the addition of strings and no primordial tensor perturbations, we find $Gμ< 2.0 \times 10^{-7}$ at 95% confidence, about 20% lower than the value obtained from previous simulations, which had 1/64 of the spatial volume. We investigate the source of the difference, showing that the main cause is an improved treatment of the string evolution across the radiation-matter transition. The increased computational volume also makes possible to simulate fully the physical equations of motion, in which the string cores shrink in comoving coordinates. This, and the larger dynamic range, changes the amplitude of the power spectra by only about 10%, demonstrating that field theory simulations of cosmic strings have now reached the required dynamic range for CMB calculations. △ Less

Submitted 27 February, 2017; v1 submitted 12 September, 2016; originally announced September 2016.

Comments: 13 pages, 6 figures. Temperature and polarization anisotropy power spectra for Abelian Higgs cosmic strings are included in the ancillary file. v2: Minor changes, matches published version. v3: Ancillary files changed after Erratum PRD 95, 049903 (2027). CMB constraints are not affected by this change

Journal ref: JCAP10(2016)042

Sound shell model for acoustic gravitational wave production at a first-order phase transition in the early Universe

Authors: Mark Hindmarsh

Abstract: A model for the acoustic production of gravitational waves at a first order phase transition is presented. The source of gravitational radiation is the sound waves generated by the explosive growth of bubbles of the stable phase. The model assumes that the sound waves are linear and that their power spectrum is determined by the characteristic form of the sound shell around the expanding bubble. T… ▽ More A model for the acoustic production of gravitational waves at a first order phase transition is presented. The source of gravitational radiation is the sound waves generated by the explosive growth of bubbles of the stable phase. The model assumes that the sound waves are linear and that their power spectrum is determined by the characteristic form of the sound shell around the expanding bubble. The predicted power spectrum has two length scales, the average bubble separation and the sound shell width when the bubbles collide. The peak of the power spectrum is at wavenumbers set by the sound shell width. For higher wavenumber $k$, the power spectrum decreases as $k^{-3}$. At wavenumbers below the inverse bubble separation, the power spectrum goes as $k^5$. For bubble wall speeds near the speed of sound where these two length scales are distinguished, there is an intermediate $k^{1}$ power law. The detailed dependence of the power spectrum on the wall speed and the other parameters of the phase transition raises the possibility of their constraint or measurement at a future space-based gravitational wave observatory such as eLISA. △ Less

Submitted 16 August, 2016; originally announced August 2016.

Comments: 5 pages

Journal ref: Phys. Rev. Lett. 120, 071301 (2018)

New solutions for non-Abelian cosmic strings

Authors: Mark Hindmarsh, Kari Rummukainen, David J. Weir

Abstract: We study the properties of classical vortex solutions in a non-Abelian gauge theory. A system of two adjoint Higgs fields breaks the SU(2) gauge symmetry to $Z_2$, producing 't Hooft-Polyakov monopoles trapped on cosmic strings, termed beads; there are two charges of monopole and two degenerate string solutions. The strings break an accidental discrete $Z_2$ symmetry of the theory, explaining the… ▽ More We study the properties of classical vortex solutions in a non-Abelian gauge theory. A system of two adjoint Higgs fields breaks the SU(2) gauge symmetry to $Z_2$, producing 't Hooft-Polyakov monopoles trapped on cosmic strings, termed beads; there are two charges of monopole and two degenerate string solutions. The strings break an accidental discrete $Z_2$ symmetry of the theory, explaining the degeneracy of the ground state. Further symmetries of the model, not previously appreciated, emerge when the masses of the two adjoint Higgs fields are degenerate. The breaking of the enlarged discrete symmetry gives rise to additional string solutions and splits the monopoles into four types of `semipole': kink solutions that interpolate between the string solutions, classified by a complex gauge invariant magnetic flux and a $Z_4$ charge. At special values of the Higgs self-couplings, the accidental symmetry broken by the string is continuous, giving rise to supercurrents on the strings. The SU(2) theory can be embedded in a wide class of Grand Unified Theories, including SO(10). We argue that semipoles and supercurrents are generic on GUT strings. △ Less

Submitted 14 December, 2016; v1 submitted 4 July, 2016; originally announced July 2016.

Comments: 6 pages, 2 figures. We dedicate this paper to Tom Kibble; v2: some references added, typos fixed, phrasing changed in places; accepted by PRL

Report number: HIP-2016-21/TH

Journal ref: Phys. Rev. Lett. 117, 251601 (2016)

Science with the space-based interferometer eLISA. II: Gravitational waves from cosmological phase transitions

Authors: Chiara Caprini, Mark Hindmarsh, Stephan Huber, Thomas Konstandin, Jonathan Kozaczuk, Germano Nardini, Jose Miguel No, Antoine Petiteau, Pedro Schwaller, Geraldine Servant, David J. Weir

Abstract: We investigate the potential for the eLISA space-based interferometer to detect the stochastic gravitational wave background produced by strong first-order cosmological phase transitions. We discuss the resulting contributions from bubble collisions, magnetohydrodynamic turbulence, and sound waves to the stochastic background, and estimate the total corresponding signal predicted in gravitational… ▽ More We investigate the potential for the eLISA space-based interferometer to detect the stochastic gravitational wave background produced by strong first-order cosmological phase transitions. We discuss the resulting contributions from bubble collisions, magnetohydrodynamic turbulence, and sound waves to the stochastic background, and estimate the total corresponding signal predicted in gravitational waves. The projected sensitivity of eLISA to cosmological phase transitions is computed in a model-independent way for various detector designs and configurations. By applying these results to several specific models, we demonstrate that eLISA is able to probe many well-motivated scenarios beyond the Standard Model of particle physics predicting strong first-order cosmological phase transitions in the early Universe. △ Less

Submitted 1 April, 2016; v1 submitted 19 December, 2015; originally announced December 2015.

Comments: 45 pages, 7 figures. Accepted for publication in JCAP, added discussion on noise from galactic binaries and citations

Gravitational waves from cosmological first order phase transitions

Authors: Mark Hindmarsh, Stephan Huber, Kari Rummukainen, David Weir

Abstract: First order phase transitions in the early Universe generate gravitational waves, which may be observable in future space-based gravitational wave observatiories, e.g. the European eLISA satellite constellation. The gravitational waves provide an unprecedented direct view of the Universe at the time of their creation. We study the generation of the gravitational waves during a first order phase tr… ▽ More First order phase transitions in the early Universe generate gravitational waves, which may be observable in future space-based gravitational wave observatiories, e.g. the European eLISA satellite constellation. The gravitational waves provide an unprecedented direct view of the Universe at the time of their creation. We study the generation of the gravitational waves during a first order phase transition using large-scale simulations of a model consisting of relativistic fluid and an order parameter field. We observe that the dominant source of gravitational waves is the sound generated by the transition, resulting in considerably stronger radiation than earlier calculations have indicated. △ Less

Submitted 14 November, 2015; originally announced November 2015.

Comments: Presented at the 33rd International Symposium on Lattice Field Theory, 14-18 July 2015, Kobe, Japan

Energy-momentum correlations for Abelian Higgs cosmic strings

Authors: David Daverio, Mark Hindmarsh, Martin Kunz, Joanes Lizarraga, Jon Urrestilla

Abstract: We report on the energy-momentum correlators obtained with recent numerical simulations of the Abelian Higgs model, essential for the computation of cosmic microwave background and matter perturbations of cosmic strings. Due to significant improvements both in raw computing power and in our parallel simulation framework, the dynamical range of the simulations has increased four-fold both in space… ▽ More We report on the energy-momentum correlators obtained with recent numerical simulations of the Abelian Higgs model, essential for the computation of cosmic microwave background and matter perturbations of cosmic strings. Due to significant improvements both in raw computing power and in our parallel simulation framework, the dynamical range of the simulations has increased four-fold both in space and time, and for the first time we are able to simulate strings with a constant physical width in both the radiation and matter eras. The new simulations improve the accuracy of the measurements of the correlation functions at the horizon scale and confirm the shape around the peak. The normalization is slightly higher in the high wave-number tails, due to a small increase in the string density. We study for the first time the behaviour of the correlators across cosmological transitions, and discover that the correlation functions evolve adiabatically, ie the network adapts quickly to changes in the expansion rate. We propose a new method for constructing source functions for Einstein-Boltzmann integrators, comparing it with two other methods previously used. The new method is more consistent, easier to implement, and significantly more accurate. △ Less

Submitted 27 February, 2017; v1 submitted 16 October, 2015; originally announced October 2015.

Comments: 22 pages, 18 figures. v2: minor changes, matches published version. v3: added Erratum that corrects the decay of vector UETCs at superhorizon scales; Planck CMB constraints are not affected by this change

Journal ref: Phys. Rev. D 93, 085014 (2016); Erratum: Phys. Rev. D 95, 049903 (2017)

Latfield2: A c++ library for classical lattice field theory

Authors: David Daverio, Mark Hindmarsh, Neil Bevis

Abstract: latfield2 is a C++ library designed to simplify writing parallel codes for solving partial differen- tial equations, developed for application to classical field theories in particle physics and cosmology. It is a significant rewrite of the latfield framework, moving from a slab domain decomposition to a rod decomposition, where the last two dimension of the lattice are scattered into a two dimens… ▽ More latfield2 is a C++ library designed to simplify writing parallel codes for solving partial differen- tial equations, developed for application to classical field theories in particle physics and cosmology. It is a significant rewrite of the latfield framework, moving from a slab domain decomposition to a rod decomposition, where the last two dimension of the lattice are scattered into a two dimensional process grid. Parallelism is implemented using the Message Passing Interface (MPI) standard, and hidden in the basic objects of grid-based simulations: Lattice, Site and Field. It comes with an integrated parallel fast Fourier transform, and I/O server class permitting computation to continue during the writing of large files to disk. latfield2 has been used for production runs on tens of thousands of processor elements, and is expected to be scalable to hundreds of thousands. △ Less

Submitted 19 January, 2016; v1 submitted 23 August, 2015; originally announced August 2015.

Comments: 14 pages, 5 figures

MSC Class: 70S99; 81T80 ACM Class: D.1.3; D.1.5; D.2.2; D.2.13

The bispectrum of cosmic string temperature fluctuations including recombination effects

Authors: Donough Regan, Mark Hindmarsh

Abstract: We calculate the cosmic microwave background temperature bispectrum from cosmic strings, for the first time including the contributions from the last scattering surface, using a well-established Gaussian model for the string energy-momentum correlation functions, and a simplified model for the cosmic fluid. We check our approximation for the integrated Sachs-Wolfe (ISW) contribution against the bi… ▽ More We calculate the cosmic microwave background temperature bispectrum from cosmic strings, for the first time including the contributions from the last scattering surface, using a well-established Gaussian model for the string energy-momentum correlation functions, and a simplified model for the cosmic fluid. We check our approximation for the integrated Sachs-Wolfe (ISW) contribution against the bispectrum obtained from the full sky map of the cosmic string ISW signal used by the Planck team, obtaining good agreement. We validate our model for the last scattering surface contribution by comparing the predicted temperature power spectrum with that obtained from a full Boltzmann code treatment applied to the Unconnected Segment Model of a string network. We find that including the last scattering contribution has only a small impact on the upper limit on the string tension resulting from the bispectrum at Planck resolutions, and argue that the bispectrum is unlikely to be competitive with the power spectrum at any resolution. △ Less

Submitted 10 August, 2015; originally announced August 2015.

Comments: 23 pages, 3 figures

Numerical simulations of acoustically generated gravitational waves at a first order phase transition

Authors: Mark Hindmarsh, Stephan J. Huber, Kari Rummukainen, David J. Weir

Abstract: We present details of numerical simulations of the gravitational radiation produced by a first order thermal phase transition in the early universe. We confirm that the dominant source of gravitational waves is sound waves generated by the expanding bubbles of the low-temperature phase. We demonstrate that the sound waves have a power spectrum with a power-law form between the scales set by the av… ▽ More We present details of numerical simulations of the gravitational radiation produced by a first order thermal phase transition in the early universe. We confirm that the dominant source of gravitational waves is sound waves generated by the expanding bubbles of the low-temperature phase. We demonstrate that the sound waves have a power spectrum with a power-law form between the scales set by the average bubble separation (which sets the length scale of the fluid flow $L_\text{f}$) and the bubble wall width. The sound waves generate gravitational waves whose power spectrum also has a power-law form, at a rate proportional to $L_\text{f}$ and the square of the fluid kinetic energy density. We identify a dimensionless parameter $\tildeΩ_\text{GW}$ characterising the efficiency of this "acoustic" gravitational wave production whose value is $8π\tildeΩ_\text{GW} \simeq 0.8 \pm 0.1$ across all our simulations. We compare the acoustic gravitational waves with the standard prediction from the envelope approximation. Not only is the power spectrum steeper (apart from an initial transient) but the gravitational wave energy density is generically larger by the ratio of the Hubble time to the phase transition duration, which can be 2 orders of magnitude or more in a typical first order electroweak phase transition. △ Less

Submitted 7 January, 2016; v1 submitted 13 April, 2015; originally announced April 2015.

Comments: 22 pages, 13 figures; v2: typos in section IIIE (Source lifetime) corrected, new appendix on GW production in an expanding universe, other minor corrections, version accepted for publication

Report number: HIP-2015-13/TH

Journal ref: Phys. Rev. D 92, 123009 (2015)

Dark Matter with Topological Defects in the Inert Doublet Model

Authors: Mark Hindmarsh, Russell Kirk, Jose Miguel No, Stephen M. West

Abstract: We examine the production of dark matter by decaying topological defects in the high mass region $m_{\mathrm{DM}} \gg m_W$ of the Inert Doublet Model, extended with an extra U(1) gauge symmetry. The density of dark matter states (the neutral Higgs states of the inert doublet) is determined by the interplay of the freeze-out mechanism and the additional production of dark matter states from the dec… ▽ More We examine the production of dark matter by decaying topological defects in the high mass region $m_{\mathrm{DM}} \gg m_W$ of the Inert Doublet Model, extended with an extra U(1) gauge symmetry. The density of dark matter states (the neutral Higgs states of the inert doublet) is determined by the interplay of the freeze-out mechanism and the additional production of dark matter states from the decays of topological defects, in this case cosmic strings. These decays increase the predicted relic abundance compared to the standard freeze-out only case, and as a consequence the viable parameter space of the Inert Doublet Model can be widened substantially. In particular, for a given dark matter annihilation rate lower dark matter masses become viable. We investigate the allowed mass range taking into account constraints on the energy injection rate from the diffuse $γ$-ray background and Big Bang Nucleosynthesis, together with constraints on the dark matter properties coming from direct and indirect detection limits. For the Inert Doublet Model high-mass region, an inert Higgs mass as low as $\sim 200$ GeV is permitted. There is also an upper limit on string mass per unit length, and hence the symmetry breaking scale, from the relic abundance in this scenario. Depending on assumptions made about the string decays, the limits are in the range $10^{12}$ GeV to $10^{13}$ GeV. △ Less

Submitted 29 July, 2015; v1 submitted 15 December, 2014; originally announced December 2014.

Comments: 27 pages, 3 figures. V2: Published version with references added

Fitting BICEP2 with defects, primordial gravitational waves and dust

Authors: Joanes Lizarraga, Jon Urrestilla, David Daverio, Mark Hindmarsh, Martin Kunz, Andrew R. Liddle

Abstract: In this work we discuss the possibility of cosmic defects being responsible for the B-mode signal measured by the BICEP2 collaboration. We also allow for the presence of other cosmological sources of B-modes such as inflationary gravitational waves and polarized dust foregrounds, which might contribute to or dominate the signal. On the one hand, we find that defects alone give a poor fit to the da… ▽ More In this work we discuss the possibility of cosmic defects being responsible for the B-mode signal measured by the BICEP2 collaboration. We also allow for the presence of other cosmological sources of B-modes such as inflationary gravitational waves and polarized dust foregrounds, which might contribute to or dominate the signal. On the one hand, we find that defects alone give a poor fit to the data points. On the other, we find that defects help to improve the fit at higher multipoles when they are considered alongside inflationary gravitational waves or polarized dust. Finally, we derive new defect constraints from models combining defects and dust. This proceeding is based on previous works [1,2]. △ Less

Submitted 1 December, 2014; originally announced December 2014.

Comments: 6 pages and 3 figures. Prepared for the Spanish Relativity Meetings ERE2014

The bispectrum of matter perturbations from cosmic strings

Authors: Donough Regan, Mark Hindmarsh

Abstract: We present the first calculation of the bispectrum of the matter perturbations induced by cosmic strings. The calculation is performed in two different ways: the first uses the unequal time correlators (UETCs) of the string network - computed using a Gaussian model previously employed for cosmic string power spectra. The second approach uses the wake model, where string density perturbations are c… ▽ More We present the first calculation of the bispectrum of the matter perturbations induced by cosmic strings. The calculation is performed in two different ways: the first uses the unequal time correlators (UETCs) of the string network - computed using a Gaussian model previously employed for cosmic string power spectra. The second approach uses the wake model, where string density perturbations are concentrated in sheet-like structures whose surface density grows with time. The qualitative and quantitative agreement of the two gives confidence to the results. An essential ingredient in the UETC approach is the inclusion of compensation factors in the integration with the Green's function of the matter and radiation fluids, and we show that these compensation factors must be included in the wake model also. We also present a comparison of the UETCs computed in the Gaussian model, and those computed in the unconnected segment model (USM) used by the standard cosmic string perturbation package CMBACT. We compare numerical estimates for the bispectrum of cosmic strings to those produced by perturbations from an inflationary era, and discover that, despite the intrinsically non-Gaussian nature of string-induced perturbations, the matter bispectrum is unlikely to produce competitive constraints on a population of cosmic strings. △ Less

Submitted 10 November, 2014; originally announced November 2014.

Comments: 24 pages, 4 figures

Constraining topological defects with temperature and polarization anisotropies

Authors: Joanes Lizarraga, Jon Urrestilla, David Daverio, Mark Hindmarsh, Martin Kunz, Andrew R. Liddle

Abstract: We analyse the possible contribution of topological defects to cosmic microwave anisotropies, both temperature and polarisation. We allow for the presence of both inflationary scalars and tensors, and of polarised dust foregrounds that may contribute to or dominate the B-mode polarisation signal. We confirm and quantify our previous statements that topological defects on their own are a poor fit t… ▽ More We analyse the possible contribution of topological defects to cosmic microwave anisotropies, both temperature and polarisation. We allow for the presence of both inflationary scalars and tensors, and of polarised dust foregrounds that may contribute to or dominate the B-mode polarisation signal. We confirm and quantify our previous statements that topological defects on their own are a poor fit to the B-mode signal. However, adding topological defects to a models with a tensor component or a dust component improves the fit around $\ell=200$. Fitting simultaneously to both temperature and polarisation data, we find that textures fit almost as well as tensors ($Δχ^2 = 2.0$), while Abelian Higgs strings are ruled out as the sole source of the B-mode signal at low $\ell$. The 95% confidence upper limits on models combining defects and dust are $Gμ< 2.7\times 10^{-7}$ (Abelian Higgs strings), $Gμ< 9.8\times 10^{-7}$ (semilocal strings) and $Gμ< 7.3\times 10^{-7}$ (textures), a small reduction on the $Planck$ bounds. The most economical fit overall is obtained by the standard $Λ$CDM model with a polarised dust component. △ Less

Submitted 3 November, 2014; v1 submitted 18 August, 2014; originally announced August 2014.

Comments: 8 pages, 3 figures. Temperature and polarization anisotropy power spectra for Abelian Higgs strings, semilocal strings and textures are included as ancillary files. Minor changes, matches published version

Journal ref: Phys. Rev. D 90, 103504 (2014)