Source author record

Mark Hindmarsh

Mark Hindmarsh appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

57works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

57 published item(s)

preprint2022arXiv

Generation of gravitational waves from freely decaying turbulence

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.

preprint2022arXiv

Thermal suppression of bubble nucleation at first-order phase transitions in the early Universe

One of the key observables in a gravitational wave power spectrum from a first order phase transition in the early Universe is the mean bubble spacing, which depends on the rate of nucleation of bubbles of the stable phase, as well as the bubble wall speed. When the bubbles expand as deflagrations, it is expected that the heating of the fluid in front of the phase boundary suppresses the nucleation rate. We quantify the effect, showing that it increases the mean bubble separation, and acts to enhance the gravitational wave signal by a factor of up to order 10. The effect is largest for small wall speeds and strong transitions.

preprint2021arXiv

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

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.

preprint2021arXiv

Effective actions and bubble nucleation from holography

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.

preprint2021arXiv

Gravitational Waves at Strong Coupling from an Effective Action

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.

preprint2021arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

The scaling density of axion strings

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%.

preprint2020arXiv

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

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.

preprint2019arXiv

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

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.

preprint2016arXiv

Improving cosmic string network simulations

In real-time lattice simulations of cosmic strings in the Abelian Higgs model, the broken translational invariance introduces lattice artefacts; relativistic strings therefore decelerate and radiate. We introduce two different methods to construct a moving string on the lattice, and study in detail the lattice effects on moving strings. We find that there are two types of lattice artefact: there is an effective maximum speed with which a moving string can be placed on the lattice, and a moving string also slows down, with the deceleration approximately proportional to the exponential of the velocity. To mitigate this, we introduce and study an improved discretisation, based on the tree-level Lüscher-Weisz action, which is found to reduce the deceleration by an order of magnitude, and to increase the string speed limit by an amount equivalent to halving the lattice spacing. The improved algorithm is expected to be very useful for 3D simulations of cosmic strings in the early universe, where one wishes to simulate as large a volume as possible.

preprint2016arXiv

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

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.

preprint2016arXiv

New solutions for non-Abelian cosmic strings

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.

preprint2016arXiv

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

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.

preprint2016arXiv

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

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.

preprint2015arXiv

Dark Matter with Topological Defects in the Inert Doublet Model

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.

preprint2015arXiv

Gravitational waves from cosmological first order phase transitions

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.

preprint2015arXiv

The bispectrum of cosmic string temperature fluctuations including recombination effects

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.

preprint2014arXiv

Can topological defects mimic the BICEP2 B-mode signal?

We show that the B-mode polarization signal detected at low multipoles by BICEP2 cannot be entirely due to topological defects. This would be incompatible with the high-multipole B-mode polarization data and also with existing temperature anisotropy data. Adding cosmic strings to a model with tensors, we find that B-modes on their own provide a comparable limit on the defects to that already coming from Planck satellite temperature data. We note that strings at this limit give a modest improvement to the best-fit of the B-mode data, at a somewhat lower tensor-to-scalar ratio of $r \simeq 0.15$.

preprint2014arXiv

Constraining topological defects with temperature and polarization anisotropies

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.

preprint2014arXiv

Dark Matter from Decaying Topological Defects

We study dark matter production by decaying topological defects, in particular cosmic strings. In topological defect or "top-down" (TD) scenarios, the dark matter injection rate varies as a power law with time with exponent $p-4$. We find a formula in closed form for the yield for all $p < 3/2$, which accurately reproduces the solution of the Boltzmann equation. We investigate two scenarios ($p=1$, $p=7/6$) motivated by cosmic strings which decay into TeV-scale states with a high branching fraction into dark matter particles. For dark matter models annihilating either by s-wave or p-wave, we find the regions of parameter space where the TD model can account for the dark matter relic density as measured by Planck. We find that topological defects can be the principal source of dark matter, even when the standard freeze-out calculation under-predicts the relic density and hence can lead to potentially large "boost factor" enhancements in the dark matter annihilation rate. We examine dark matter model-independent limits on this scenario arising from unitarity and discuss example model-dependent limits coming from indirect dark matter search experiments. In the four cases studied, the upper bound on $Gμ$ for strings with an appreciable channel into TeV-scale states is significantly more stringent than the current Cosmic Microwave Background limits.

preprint2014arXiv

Fitting BICEP2 with defects, primordial gravitational waves and dust

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].

preprint2014arXiv

Gravitational waves from the sound of a first order phase transition

We report on the first 3-dimensional numerical simulations of first-order phase transitions in the early universe to include the cosmic fluid as well as the scalar field order parameter. We calculate the gravitational wave (GW) spectrum resulting from the nucleation, expansion and collision of bubbles of the low-temperature phase, for phase transition strengths and bubble wall velocities covering many cases of interest. We find that the compression waves in the fluid continue to be a source of GWs long after the bubbles have merged, a new effect not taken properly into account in previous modelling of the GW source. For a wide range of models the main source of the GWs produced by a phase transition is therefore the sound the bubbles make.

preprint2014arXiv

The bispectrum of matter perturbations from cosmic strings

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.

preprint2013arXiv

Consistent cosmology with Higgs thermal inflation in a minimal extension of the MSSM

We consider a class of models in which minimal gauged F-term hybrid inflation is coupled renormalisably to the minimal supersymmetric standard model (MSSM), with no extra ingredients; we call this class the "minimal hybrid inflationary supersymmetric standard model" (MHISSM). The singlet inflaton supplies the Higgs mu-term, and allows an exit from inflation to a vacuum characterised by large Higgs vevs. The true ground state is reached after an period of thermal inflation along the Higgs flat direction. The scalar spectral index is reduced from the standard F-term value, to approximately 0.976 in the case where the inflaton potential is dominated by the 1-loop corrections. The reheat temperature following thermal inflation is about 10^9 GeV, solving the gravitino overclosure problem. A Higgs condensate reduces the cosmic string mass per unit length, rendering it compatible with the Cosmic Microwave Background constraints without tuning the inflaton coupling. With the minimal U(1)' gauge symmetry in the inflaton sector, where one of the waterfall fields generates a right-handed neutrino mass, we investigate the Higgs thermal inflation scenario in three popular supersymmetry-breaking schemes: AMSB, GMSB and the CMSSM, focusing on the constraints on the gravitino mass.

preprint2013arXiv

Exact Scale-Invariant Background of Gravitational Waves from Cosmic Defects

We demonstrate that any scaling source in the radiation era produces a background of gravitational waves with an exact scale-invariant power spectrum. Cosmic defects, created after a phase transition in the early Universe, are such a scaling source. We emphasise that the result is independent of the topology of the cosmic defects, the order of phase transition, and the nature of the symmetry broken, global or gauged. As an example, using large-scale numerical simulations, we calculate the scale invariant gravitational wave power spectrum generated by the dynamics of a global O(N) scalar theory. The result approaches the large N theoretical prediction as N^(-2), albeit with a large coefficient. The signal from global cosmic strings is O(100) times larger than the large N prediction.

preprint2012arXiv

f(R) Gravity from the renormalisation group

We explore the cosmological dynamics of an effective f(R) model constructed from a renormalisation group (RG) improvement of the Einstein--Hilbert action, using the non-perturbative beta functions of the exact renormalisation group equation. The resulting f(R) model has some remarkable properties. It naturally exhibits an unstable de Sitter era in the ultraviolet (UV), dynamically connected to a stable de Sitter era in the IR, via a period of radiation and matter domination, thereby describing a non-singular universe. We find that the UV de Sitter point is one of an infinite set, which make the UV RG fixed point inaccessible to classical cosmological evolution. In the vicinity of the fixed point, the model behaves as R^2 gravity, while it correctly recovers General Relativity at solar system scales. In this simplified model, the fluctuations are too large to be the observed ones, and more ingredients in the action are needed.

preprint2012arXiv

Particle motion in weak relativistic gravitational fields

We derive the geodesic equation of motion in the presence of weak gravitational fields produced by relativistic sources such as cosmic strings, decomposed into scalar, vector and tensor parts. We find that the vector (gravito-magnetic) force is an important contributor, and for non-relativistic particles we recover the well-known result for the impulse from a moving straight string. Our results can be straightforwardly incorporated into N-body simulations to allow for the presence of cosmic defects or other sources of weak gravitational fields.

preprint2012arXiv

Radiation and Relaxation of Oscillons

We study oscillons, extremely long-lived localized oscillations of a scalar field, with three different potentials: quartic, sine-Gordon model and in a new class of convex potentials. We use an absorbing boundary at the end of the lattice to remove emitted radiation. The energy and the frequency of an oscillon evolve in time and are well fitted by a constant component and a decaying, radiative part obeying a power law as a function time. The power spectra of the emitted radiation show several distinct frequency peaks where oscillons release energy. In two dimensions, and with suitable initial conditions, oscillons do not decay within the range of the simulations, which in quartic theory reach 10^8 time units. While it is known that oscillons in three-dimensional quartic theory and sine-Gordon model decay relatively quickly, we observe a surprising persistence of the oscillons in the convex potential with no sign of demise up to 10^7 time units. This leads us to speculate that an oscillon in such a potential could actually live infinitely long both in two and three dimensions.

preprint2012arXiv

Renormalisation group improvement of scalar field inflation

We study quantum corrections to Friedmann-Robertson-Walker cosmology with a scalar field under the assumption that the dynamics are subject to renormalisation group improvement. We use the Bianchi identity to relate the renormalisation group scale to the scale factor and obtain the improved cosmological evolution equations. We study the solutions of these equations in the renormalisation group fixed point regime, obtaining the time-dependence of the scalar field strength and the Hubble parameter in specific models with monomial and trinomial quartic scalar field potentials. We find that power-law inflation can be achieved in the renormalisation group fixed point regime with the trinomial potential, but not with the monomial one. We study the transition to the quasi-classical regime, where the quantum corrections to the couplings become small, and find classical dynamics as an attractor solution for late times. We show that the solution found in the renormalisation group fixed point regime is also a cosmological fixed point in the autonomous phase space. We derive the power spectrum of cosmological perturbations and find that the scalar power spectrum is exactly scale-invariant and bounded up to arbitrarily small times, while the tensor perturbations are tilted as appropriate for the background power-law inflation. We specify conditions for the renormalisation group fixed point values of the couplings under which the amplitudes of the cosmological perturbations remain small.

preprint2012arXiv

Strictly Anomaly Mediated Supersymmetry Breaking

We consider an MSSM extension with anomaly mediation as the source of supersymmetry-breaking, and a U(1) symmetry which solves the tachyonic slepton problem, and introduces both the see-saw mechanism for neutrino masses, and the Higgs mu-term. We compare its spectra with those from so-called minimal anomaly mediated supersymmetry breaking. We find a Standard Model-like Higgs of mass 125 GeV with a gravitino mass of 140 TeV and tan(beta)=16. However, the muon anomalous magnetic moment is 3 sigma away from the experimental value. The model naturally produces a period of hybrid inflation, which can exit to a false vacuum characterised by large Higgs vevs, reaching the true ground state after a period of thermal inflation. The scalar spectral index is reduced to approximately 0.975, and the correct abundance of neutralino dark matter can be produced by decays of thermally-produced gravitinos, provided the gravitino mass (and hence the Higgs mass) is high. Naturally light cosmic strings are produced, satisfying bounds from the Cosmic Microwave Background. The complementary pulsar timing and cosmic ray bounds require that strings decay primarily via loops into gravitational waves. Unless the loops are extremely small, the next generation pulsar timing array will rule out or detect the string-derived gravitational radiation background in this model.

preprint2011arXiv

Asymptotically Safe Cosmology

We study quantum modifications to cosmology in a Friedmann-Robertson-Walker universe with and without scalar fields by taking the renormalisation group running of gravitational and matter couplings into account. We exploit the Bianchi identity to relate the renormalisation group scale with scale factor and derive the improved cosmological evolution equations. We find two types of cosmological fixed points where the renormalisation group scale either freezes in, or continues to evolve with scale factor. We discuss the implications of each of these, and classify the different cosmological fixed points with and without gravity displaying an asymptotically safe renormalisation group fixed point. We state conditions of existence for an inflating ultraviolet cosmological fixed point for Einstein gravity coupled to a scalar field. We also discuss other fixed point solutions such as "scaling" solutions, or fixed points with equipartition between kinetic and potential energies.

preprint2011arXiv

Covariant Vertex Operators for Cosmic Strings

We construct complete sets of (open and closed string) covariant coherent state and mass eigenstate vertex operators in bosonic string theory. By minimally extending the standard definition of coherent states so as to include the string theory requirements, we show that the naive construction of the the closed string coherent states requires the existence of a lightlike compactification of spacetime. When the null winding states in the underlying Hilbert space are projected out the resulting vertex operators satisfy the definition of a coherent state and have a classical interpretation. We present explicitly both the covariant and lightcone gauge realization of the resulting states using the DDF map that relates the two. We also identify the corresponding general lightcone gauge classical solutions around which the quantum states are fluctuating. We go on to show that both the covariant gauge coherent vertex operators, the corresponding lightcone gauge coherent states and the classical solutions all share the same mass and angular momenta and conjecture that the covariant and lightcone gauge states are different manifestations of the same state and share identical interactions. This construction can be used to study the evolution of fundamental cosmic strings as predicted by string theory and may also be useful for other applications where massive string vertex operators are of interest.

preprint2011arXiv

String Vertex Operators and Cosmic Strings

We construct complete sets of (open and closed string) covariant coherent state and mass eigenstate vertex operators in bosonic string theory. This construction can be used to study the evolution of fundamental cosmic strings as predicted by string theory, and is expected to serve as a self-contained prototype toy model on which realistic cosmic superstring vertex operators can be based on. It is also expected to be useful for other applications where massive string vertex operators are of interest. We pay particular attention to all the normalization constants, so that these vertices lead directly to unitary S-matrix elements.

preprint2011arXiv

The dynamical equivalence of modified gravity revisited

We revisit the dynamical equivalence between different representations of vacuum modified gravity models in view of Legendre transformations. The equivalence is discussed for both bulk and boundary space, by including in our analysis the relevant Gibbons-Hawking terms. In the f(R) case, the Legendre transformed action coincides with the usual Einstein frame one. We then re-express the R+f(G) action, where G is the Gauss-Bonnet term, as a second order theory with a new set of field variables, four tensor fields and one scalar and study its dynamics. For completeness, we also calculate the conformal transformation of the full Jordan frame R+f(G) action. All the appropriate Gibbons-Hawking terms are calculated explicitly.

preprint2010arXiv

CMB power spectra from cosmic strings: predictions for the Planck satellite and beyond

We present a significant improvement over our previous calculations of the cosmic string contribution to cosmic microwave background (CMB) power spectra, with particular focus on sub-WMAP angular scales. These smaller scales are relevant for the now-operational Planck satellite and additional sub-orbital CMB projects that have even finer resolutions. We employ larger Abelian Higgs string simulations than before and we additionally model and extrapolate the statistical measures from our simulations to smaller length scales. We then use an efficient means of including the extrapolations into our Einstein-Boltzmann calculations in order to yield accurate results over the multipole range 2 < l < 4000. Our results suggest that power-law behaviour cuts in for l > 3000 in the case of the temperature power spectrum, which then allows cautious extrapolation to even smaller scales. We find that a string contribution to the temperature power spectrum making up 10% of power at l=10 would be larger than the Silk-damped primary adiabatic contribution for l > 3500. Astrophysical contributions such as the Sunyaev-Zeldovich effect also become important at these scales and will reduce the sensitivity to strings, but these are potentially distinguishable by their frequency-dependence.

preprint2010arXiv

CMB temperature trispectrum of cosmic strings

We provide an analytical expression for the trispectrum of the Cosmic Microwave Background (CMB) temperature anisotropies induced by cosmic strings. Our result is derived for the small angular scales under the assumption that the temperature anisotropy is induced by the Gott-Kaiser-Stebbins effect. The trispectrum is predicted to decay with a non-integer power-law exponent l^(-r) with 6<r<7, depending on the string microstructure, and thus on the string model. For Nambu-Goto strings, this exponent is related to the string mean square velocity and the loop distribution function. We then explore two classes of wavenumber configuration in Fourier space, the kite and trapezium quadrilaterals. The trispectrum can be of any sign and appears to be strongly enhanced for all squeezed quadrilaterals.

preprint2010arXiv

Covariant Closed String Coherent States

We give the first construction of covariant coherent closed string states, which may be identified with fundamental cosmic strings. We outline the requirements for a string state to describe a cosmic string, and using DDF operators provide an explicit and simple map that relates three different descriptions: classical strings, lightcone gauge quantum states and covariant vertex operators. The naive construction leads to covariant vertex operators whose existence requires a lightlike compactification of spacetime. When the lightlike compactified states in the underlying Hilbert space are projected out the resulting coherent states have a classical interpretation and are in one-to-one correspondence with arbitrary classical closed string loops.

preprint2010arXiv

Detecting and distinguishing topological defects in future data from the CMBPol satellite

The proposed CMBPol mission will be able to detect the imprint of topological defects on the cosmic microwave background (CMB) provided the contribution is sufficiently strong. We quantify the detection threshold for cosmic strings and for textures, and analyse the satellite's ability to distinguish between these different types of defects. We also assess the level of danger of misidentification of a defect signature as from the wrong defect type or as an effect of primordial gravitational waves. A 0.002 fractional contribution of cosmic strings to the CMB temperature spectrum at multipole ten, and similarly a 0.001 fractional contribution of textures, can be detected and correctly identified at the 3σ level. We also confirm that a tensor contribution of r = 0.0018 can be detected at over 3σ, in agreement with the CMBpol mission concept study. These results are supported by a model selection analysis.

preprint2008arXiv

Abelian Higgs Cosmic Strings: Small Scale Structure and Loops

Classical lattice simulations of the Abelian Higgs model are used to investigate small scale structure and loop distributions in cosmic string networks. Use of the field theory ensures that the small-scale physics is captured correctly. The results confirm analytic predictions of Polchinski & Rocha [1] for the two-point correlation function of the string tangent vector, with a power law from length scales of order the string core width up to horizon scale with evidence to suggest that the small scale structure builds up from small scales. An analysis of the size distribution of string loops gives a very low number density, of order 1 per horizon volume, in contrast with Nambu-Goto simulations. Further, our loop distribution function does not support the detailed analytic predictions for loop production derived by Dubath et al. [2]. Better agreement to our data is found with a model based on loop fragmentation [3], coupled with a constant rate of energy loss into massive radiation. Our results show a strong energy loss mechanism which allows the string network to scale without gravitational radiation, but which is not due to the production of string width loops. From evidence of small scale structure we argue a partial explanation for the scale separation problem of how energy in the very low frequency modes of the string network is transformed into the very high frequency modes of gauge and Higgs radiation. We propose a picture of string network evolution which reconciles the apparent differences between Nambu-Goto and field theory simulations.

preprint2006arXiv

Simulations of Cold Electroweak Baryogenesis: Finite time quenches

The electroweak symmetry breaking transition may supply the appropriate out-of-equilibrium conditions for baryogenesis if it is triggered sufficiently fast. This can happen at the end of low-scale inflation, prompting baryogenesis to occur during tachyonic preheating of the Universe, when the potential energy of the inflaton is transfered into Standard Model particles. With the proper amount of CP-violation present, the observed baryon number asymmetry can be reproduced. Within this framework of Cold Electroweak Baryogenesis, we study the dependence of the generated baryon asymmetry on the speed of the quenching transition. We find that there is a separation between ``fast'' and ``slow'' quenches, which can be used to put bounds on the allowed Higgs-inflaton coupling. We also clarify the strong Higgs mass dependence of the asymmetry reported in a companion paper (hep-ph/0604263).

preprint2003arXiv

Kink-boundary collisions in a two dimensional scalar field theory

In a two-dimensional toy model, motivated from five-dimensional heterotic M-theory, we study the collision of scalar field kinks with boundaries. By numerical simulation of the full two-dimensional theory, we find that the kink is always inelastically reflected with a model-independent fraction of its kinetic energy converted into radiation. We show that the reflection can be analytically understood as a fluctuation around the scalar field vacuum. This picture suggests the possibility of spontaneous emission of kinks from the boundary due to small perturbations in the bulk. We verify this picture numerically by showing that the radiation emitted from the collision of an initial single kink eventually leads to a bulk populated by many kinks. Consequently, processes changing the boundary charges are practically unavoidable in this system. We speculate that the system has a universal final state consisting of a stack of kinks, their number being determined by the initial energy.

preprint1998arXiv

Constraints on Axion Models from $K^+\toπ^+ a$

We explore a new class of axion models in which some, but not all, of the left-handed quarks have a Peccei-Quinn symmetry. These models are potentially afflicted by flavour changing neutral currents. We derive the bounds on the Peccei-Quinn symmetry-breaking scale from bounds on the $K^+\to π^+ a$ branching ratio, showing that in some cases they are even stronger than the astrophysical ones, but still not strong enough to kill off the models.

preprint1998arXiv

Cosmic Strings - Dead Again?

I report on recent numerical simulations of the simplest field theory with cosmic string solutions, the Abelian Higgs model. We find that random networks of string quickly converge to a scaling solution in which the network scale length $ξ$ increases linearly with time. There are very few loops with sizes less than $ξ$, and the strings are smooth, showing no signs of ``small scale structure''. We claim that particle production is the dominant energy-loss mechanism, not gravitational radiation as previously thought. For strings in Grand Unified Models, stringent constraints can be placed from cosmic ray observations on the string tension $μ$: we estimate $Gμ< 10^{-9}$, three orders of magnitude lower than the constraint from Cosmic Microwave Background fluctuations.

preprint1998arXiv

Sphalerons with CP-Violating Higgs Potentials

We investigate the effect on the sphaleron in the two Higgs doublet electroweak theory of including CP violation in the Higgs potential. To have better control over the relation between the sphaleron energy and the physical quantities in the theory, we show how to parametrize the Higgs potential in terms of physical masses and mixing angles, one of which causes CP violation. By altering this CP violating angle (and keeping the other physical quantities fixed) the sphaleron energy increases by up to 10%. We also calculate the static minimum energy path between adjacent vacua as a function of Chern-Simons number, using the method of gradient flow. The only effect CP violation has on the barrier is the change in height. As a by-product of our work on parametrization of the potential, we demonstrate that CP violation in the Higgs sector favours nearly degenerate light Higgs masses.

preprint1997arXiv

Constraints on Variant Axion Models

A particular class of variant axion models with two higgs doublets and a singlet is studied. In these models the axion couples either to the $u$-quark or $t$-quark or both, but not to $b$, $c$, $s$, or $d$. When the axion couples to only one quark the models possess the desirable feature of having no domain wall problem, which makes them viable candidates for a cosmological axion string scenario. We calculate the axion couplings to leptons, photons and nucleons, and the astrophysical constraints on the axion decay constant $v_a$ are investigated and compared to the DFSZ axion model. We find that the most restrictive lower bound on $v_a$, that from SN1987a, is lowered by up to a factor of about 30, depending on the model and also the ratio of the vacuum expectation values of the higgs doublets. For scenarios with axionic strings, the allowed window for $v_a$ in the $u$ quark model can be more than two orders of magnitude. For inflationary scenarios, the cosmological upper bound on $v_a/N$, where $N$ is the QCD anomaly factor, is unaffected: however, the variant models have $N$ either 3 or 6 times smaller than the DFSZ model.

preprint1997arXiv

Magnetic Fields from Phase Transitions

The generation of primordial magnetic fields from cosmological phase transitions is discussed, paying particular attention to the electroweak transition and to the various definitions of the `average' field that have been put forward. It is emphasised that only the volume average has dynamical significance as a seed for galactic dynamos. On rather general grounds of causality and energy conservation, it is shown that, in the absence of MHD effects that transfer power in the magnetic field from small to large scales, processes occurring at the electroweak transition cannot generate fields stronger than $10^{-20}$ Gauss on a scale of 0.5 Mpc. However, it is implausible that this upper bound could ever be reached, as it would require all the energy in the Universe to be turned into a magnetic field coherent at the horizon scale. Non-linear MHD effects seem therefore to be necessary if the electroweak transition is to create a primordial seed field.

preprint1997arXiv

On the instability of the one-texture universe

The one-texture universe, introduced by Davis in 1987, is a homogeneous mapping of a scalar field with an $S^3$ vacuum into a closed universe. It has long been known to mathematicians that such solutions, although static, are unstable. We show by explicit construction that there is only one unstable mode, corresponding to collapse of the texture towards a single point, in the case where gravitational backreaction is neglected. We discuss the instability timescale in both static and expanding space-times; in the latter case it is of order of the present age of the universe, suggesting that, though unstable, the one-texture universe could survive to the present. The cosmic microwave background constrains the initial magnitude of this unstable perturbation to be less than of order 10^{-3}.

preprint1997arXiv

The Thermodynamics of Cosmic String densities in U(1) Scalar Field Theory

We present a full characterization of the phase transition in U(1) scalar field theory and of the associated vortex string thermodynamics in 3D. We show that phase transitions in the string densities exist and measure their critical exponents, both for the long string and the short loops. Evidence for a natural separation between these two string populations is presented. In particular our results strongly indicate that an infinite string population will only exist above the critical temperature. Canonical initial conditions for cosmic string evolution are show to correspond to the infinite temperature limit of the theory.

preprint1996arXiv

Analytic scaling solutions for cosmic domain walls

A relativistic generalisation of a well-known method for approximating the dynamics of topological defects in condensed matter is constructed, and applied to the evolution of domain walls in a cosmological context. It is shown that there are self-similar ``scaling'' solutions, for which one can in principle calculate many quantities of interest without recourse to numerical simulations. Here, the area density in the scaling regime is calculated in various backgrounds. Remarkably good agreement with numerical simulations is obtained.

preprint1995arXiv

Where are the Hedgehogs in Nematics?

In experiments which take a liquid crystal rapidly from the isotropic to the nematic phase, a dense tangle of defects is formed. In nematics, there are in principle both line and point defects (``hedgehogs''), but no point defects are observed until the defect network has coarsened appreciably. In this letter the expected density of point defects is shown to be extremely low, approximately $10^{-8}$ per initially correlated domain, as result of the topology (specifically, the homology) of the order parameter space.

preprint1994arXiv

Cosmological Perturbations from Cosmic Strings

Some aspects of the theory of cosmological perturbations from cosmic strings and other topological defects are outlined, with particular reference to a simple example: a spatially flat CDM-dominated universe. The conserved energy-momentum pseudo-tensor is introduced, and the equation for the density perturbation derived from it. It is shown how the scaling hypothesis for defect evolution results in a Harrison-Zel'dovich spectrum for wavelengths well inside the horizon.

preprint1994arXiv

Small Scale Microwave Background Fluctuations from Cosmic Strings

The Cosmic Microwave Background (CMB) fluctuations at very small angular scales (less than $10'$) induced by matter sources are computed in a simplified way. The result corrects a previous formula appearing in the literature. The small scale power spectrum from cosmic strings is then calculated by a new analytic method. The result compares extremely well with the spectrum computed by numerical techniques (when the old, incorrect, formula is used). The upper bound on the string parameters derived from OVRO data is re-examined, taking into account the non-Gaussian nature of stringy perturbations on small scales. Assuming a conventional ionization history, the bound is $γGμ< 11\times 10^{-6}$, where $γ^2$ is the number of horizon lengths of string per horizon volume. Current simulations give $γ^2 = 31\pm 7$.

preprint1992arXiv

Generalized Semilocal Theories and Higher Hopf Maps

\def\mon{S^3\stackrel{S^1}{\rightarrow}S^2} \def\inst{S^7\stackrel{S^3}{\rightarrow}S^4} \def\octo{S^{15}\stackrel{S^7}{\rightarrow}S^8} In semilocal theories, the vacuum manifold is fibered in a non-trivial way by the action of the gauge group. Here we generalize the original semilocal theory (which was based on the Hopf bundle $\mon$) to realize the next Hopf bundle $\inst$, and its extensions $S^{2n+1}\stackrel{S^3}\rightarrow \H P^n$. The semilocal defects in this class of theories are classified by $π_3(S^3)$, and are interpreted as constrained instantons or generalized sphaleron configurations. We fail to find a field theoretic realization of the final Hopf bundle $\octo$, but are able to construct other semilocal spaces realizing Stiefel bundles over Grassmanian spaces.

preprint1992arXiv

Semilocal Topological Defects

Semilocal defects are those formed in field theories with spontaneously broken symmetries, where the vacuum manifold $M$ is fibred by the action of the gauge group in a non-trivial way. Studied in this paper is the simplest such class of theories, in which $M\simeq S^{2N-1}$, fibred by the action of a local $U(1)$ symmetry. Despite $M$ having trivial homotopy groups up to $π_{2N-2}$, this theory exhibits a fascinating variety of defects: vortices, or semilocal strings; monopoles (on which the strings terminate); and (when $N=2$) textures, which may be stabilised by their associated magnetic field to produce a skyrmion.