Source author record

Richard A. Battye

Richard A. Battye 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

41works
8topics
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

41 published item(s)

preprint2021arXiv

The BINGO Project VII: Cosmological Forecasts from 21cm Intensity Mapping

The 21cm line of neutral hydrogen (HI) opens a new avenue in our exploration of the structure and evolution of the Universe. It provides complementary data to the current large-scale structure observations with different systematics, and thus it will be used to improve our understanding of the $Λ$CDM model. Among several radio cosmological surveys designed to measure this line, BINGO is a single-dish telescope mainly designed to detect baryon acoustic oscillations (BAOs) at low redshifts ($0.127< z<0.449$). Our goal is to assess the fiducial BINGO setup and its capabilities of constraining the cosmological parameters, and to analyze the effect of different instrument configurations. We used the Phase 1 fiducial configuration of the BINGO telescope to perform our cosmological forecasts. In addition, we investigated the impact of several instrumental setups, taking into account some instrumental systematics, and different cosmological models. Combining BINGO with Planck temperature and polarization data, the projected constraint improves from a $13\%$ and $25\%$ precision measurement at the $68\%$ confidence level with Planck only to $1\%$ and $3\%$ for the Hubble constant and the dark energy equation of state (EoS), respectively, within the wCDM model. Assuming a Chevallier-Polarski-Linder parameterization, the EoS parameters have standard deviations given by $σ_{w_0} = 0.30$ and $σ_{w_a} = 1.2$, which are improvements on the order of $30\%$ with respect to Planck alone. Also, we can access information about the HI density and bias, obtaining $\sim 8.5\%$ and $\sim 6\%$ precision, respectively, assuming they vary with redshift at three independent bins. The fiducial BINGO configuration will be able to extract significant cosmological information from the HI distribution and provide constraints competitive with current and future cosmological surveys. (Abridged)

preprint2020arXiv

Dark matter axion detection in the radio/mm-waveband

We discuss axion dark matter detection via two mechanisms: spontaneous decays and resonant conversion in neutron star magnetospheres. For decays, we show that the brightness temperature signal, rather than flux, is a less ambiguous measure for selecting candidate objects. This is owing principally to the finite beam width of telescopes which prevents one from being sensitive to the total flux from the object. With this in mind, we argue that the large surface-mass-density of the galactic centre or the Virgo cluster centre offers the best chance of improving current constraints on the axion-photon coupling via spontaneous decays. For the neutron star case, we first carry out a detailed study of mixing in magnetised plasmas. We derive transport equations for the axion-photon system via a controlled gradient expansion, allowing us to address inhomogeneous mass-shell constraints for arbitrary momenta. We then derive a non-perturbative Landau-Zener formula for the conversion probability valid across the range of relativistic and non-relativistic axions and show that the standard perturbative resonant conversion amplitude is a truncation of this result in the non-adiabatic limit. Our treatment reveals that that infalling dark matter axions typically convert non-adiabatically in magnetospheres. We describe the limitations of one-dimensional mixing equations and explain how three-dimensional effects activate new photon polarisations, including longitudinal modes and illustrate these arguments with numerical simulations in higher dimensions. We find that the bandwidth of the radio signal is dominated by Doppler broadening from the relative motion of the neutron star with respect to the observer. Therefore, we conclude that the radio signal from the resonant decay is weaker than previously thought, which means one relies on local density peaks to probe weaker axion-photon couplings.

preprint2020arXiv

Domain Wall Constraints on Two Higgs Doublet Models with $Z_2$ Symmetry

The Two Higgs Doublet Model (2HDM) with spontaneously broken $Z_2$ symmetry predicts a production of domain walls at the electroweak scale. We derive cosmological constraints on model parameters for both Type-I and Type-II 2HDMs from the requirement that domain walls do not dominate the Universe by the present day. For Type-I 2HDMs, we deduce the lower bound on the key parameter $\tanβ> 10^5$ for a wide range of Higgs-boson masses $\sim$ 100 GeV or greater close to the Standard Model alignment limit. In addition, we perform numerical simulations of the 2HDM with an approximate as well as an exact $Z_2$ symmetry but biased initial conditions. In both cases, we find that domain wall networks are unstable and, hence, do not survive at late times. The domain walls experience an exponential suppression of scaling in these models which can help ameliorate the stringent constraints found in the case of an exact discrete symmetry. For a 2HDM with softly-broken $Z_2$ symmetry, we relate the size of this exponential suppression to the soft-breaking bilinear parameter $m_{12}$ allowing limits to be placed on this parameter of order $μ$eV, such that domain wall domination can be avoided. In particular, for Type-II 2HDMs, we obtain a corresponding lower limit on the CP-odd phase $θ$ generated by QCD instantons, $θ\ \stackrel{>}{{}_\sim}\ 10^{-11}/(\sinβ\cosβ)$, which is in some tension with the upper limit of $θ\ \stackrel{<}{{}_\sim}\ 10^{-11}$--$10^{-10}$, as derived from the non-observation of a non-zero neutron electric dipole moment. For a $Z_2$-symmetric 2HDM with biased initial conditions, we are able to relate the size of the exponential suppression to a biasing parameter $\varepsilon$ so as to avoid domain wall domination.

preprint2020arXiv

Simulations of Domain Walls in Two Higgs Doublet Models

The Two Higgs Doublet Model predicts the emergence of 3 distinct domain wall solutions arising from the breaking of 3 accidental global symmetries, $Z_2$, CP1 and CP2, at the electroweak scale for specific choices of the model parameters. We present numerical kink solutions to the field equations in all three cases along with dynamical simulations of the models in (2+1) and (3+1) dimensions. For each kink solution we define an associated topological current. In all three cases simulations produce a network of domain walls which deviates from power law scaling in Minkowski and FRW simulations. This deviation is attributed to a winding of the electroweak group parameters around the domain walls in our simulations. We observe a local violation of the neutral vacuum condition on the domain walls in our simulations. This violation is attributed to relative electroweak transformations across the domain walls which is a general feature emerging from random initial conditions.

preprint2020arXiv

SuperCLASS -- I. The Super CLuster Assisted Shear Survey: Project overview and Data Release 1

The SuperCLuster Assisted Shear Survey (SuperCLASS) is a legacy programme using the e-MERLIN interferometric array. The aim is to observe the sky at L-band (1.4 GHz) to a r.m.s. of 7 uJy per beam over an area of ~1 square degree centred on the Abell 981 supercluster. The main scientific objectives of the project are: (i) to detect the effects of weak lensing in the radio in preparation for similar measurements with the Square Kilometre Array (SKA); (ii) an extinction free census of star formation and AGN activity out to z~1. In this paper we give an overview of the project including the science goals and multi-wavelength coverage before presenting the first data release. We have analysed around 400 hours of e-MERLIN data allowing us to create a Data Release 1 (DR1) mosaic of ~0.26 square degrees to the full depth. These observations have been supplemented with complementary radio observations from the Karl G. Jansky Very Large Array (VLA) and optical/near infra-red observations taken with the Subaru, Canada-France-Hawaii and Spitzer Telescopes. The main data product is a catalogue of 887 sources detected by the VLA, of which 395 are detected by e-MERLIN and 197 of these are resolved. We have investigated the size, flux and spectral index properties of these sources finding them compatible with previous studies. Preliminary photometric redshifts, and an assessment of galaxy shapes measured in the radio data, combined with a radio-optical cross-correlation technique probing cosmic shear in a supercluster environment, are presented in companion papers.

preprint2020arXiv

SuperCLASS -- III. Weak lensing from radio and optical observations in Data Release 1

We describe the first results on weak gravitational lensing from the SuperCLASS survey: the first survey specifically designed to measure the weak lensing effect in radio-wavelength data, both alone and in cross-correlation with optical data. We analyse 1.53 square degrees of optical data from the Subaru telescope and 0.26 square degrees of radio data from the e-MERLIN and VLA telescopes (the DR1 data set). Using standard methodologies on the optical data only we make a significant (10 sigma) detection of the weak lensing signal (a shear power spectrum) due to the massive supercluster of galaxies in the targeted region. For the radio data we develop a new method to measure the shapes of galaxies from the interferometric data, and we construct a simulation pipeline to validate this method. We then apply this analysis to our radio observations, treating the e-MERLIN and VLA data independently. We achieve source densities of 0.5 per square arcmin in the VLA data and 0.06 per square arcmin in the e-MERLIN data, numbers which prove too small to allow a detection of a weak lensing signal in either the radio data alone or in cross-correlation with the optical data. Finally, we show preliminary results from a visibility-plane combination of the data from e-MERLIN and VLA which will be used for the forthcoming full SuperCLASS data release. This approach to data combination is expected to enhance both the number density of weak lensing sources available and the fidelity with which their shapes can be measured.

preprint2018arXiv

Cosmology with Phase 1 of the Square Kilometre Array; Red Book 2018: Technical specifications and performance forecasts

We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z = 0.35 - 3; and a deep, high-redshift HI intensity mapping survey over 100 sqdeg from z = 3 - 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to z ~ 3 with percent-level precision measurements of the cosmic expansion rate; constraining possible deviations from General Relativity on cosmological scales by measuring the growth rate of structure through multiple independent methods; mapping the structure of the Universe on the largest accessible scales, thus constraining fundamental properties such as isotropy, homogeneity, and non-Gaussianity; and measuring the HI density and bias out to z = 6. These surveys will also provide highly complementary clustering and weak lensing measurements that have independent systematic uncertainties to those of optical surveys like LSST and Euclid, leading to a multitude of synergies that can improve constraints significantly beyond what optical or radio surveys can achieve on their own. This document, the 2018 Red Book, provides reference technical specifications, cosmological parameter forecasts, and an overview of relevant systematic effects for the three key surveys, and will be regularly updated by the Cosmology Science Working Group in the run up to start of operations and the Key Science Programme of SKA1.

preprint2016arXiv

Approximation of the potential in scalar field dark energy models

We study the nature of potentials in scalar field based models for dark energy - with both canonical and noncanonical kinetic terms. We calculate numerically, and using an analytic approximation around $a\approx 1$, potentials for models with constant equation-of-state parameter, $w_ϕ$. We find that for a wide range of models with canonical and noncanonical kinetic terms there is a simple approximation for the potential that holds when the scale factor is in the range $0.6\lesssim a\lesssim 1.4$. We discuss how this form of the potential can also be used to represent models with non-constant $w_ϕ$ and, hence, how it could be used in reconstruction from cosmological data.

preprint2016arXiv

The Equation of State Approach to Cosmological Perturbations in f(R) Gravity

The discovery of apparent cosmological acceleration has spawned a huge number of dark energy and modified gravity theories. The f(R) models of gravity are obtained when one replaces the Ricci scalar in the Einstein-Hilbert action by an arbitrary function of the Ricci scalar. In this work we provide expressions for the equations of state (EoS) for perturbations which completely characterize the linearized perturbations in f(R) gravity, including the scalar, vector, and tensor modes. The EoS formalism is a powerful and elegant parametrization where the modification to General Relativity are treated as a dark-energy-fluid. The perturbed dark-energy-fluid variables such as the anisotropic stress or the entropy perturbation are explicitly given in terms of parameters of the model of interest.

preprint2015arXiv

A demonstration of position angle-only weak lensing shear estimators on the GREAT3 simulations

We develop and apply the position angle-only shear estimator of Whittaker et al. (2014) to realistic galaxy images. This is done by demonstrating the method on the simulations of the GREAT3 challenge (Mandelbaum et al. 2014b), which include contributions from anisotropic PSFs. We measure the position angles of the galaxies using three distinct methods - the integrated light method, quadrupole moments of surface brightness, and using model-based ellipticity measurements provided by IM3SHAPE. A weighting scheme is adopted to address biases in the position angle measurements which arise in the presence of an anisotropic PSF. Biases on the shear estimates, due to measurement errors on the position angles and correlations between the measurement errors and the true position angles, are corrected for using simulated galaxy images and an iterative procedure. The properties of the simulations are estimated using the deep field images provided as part of the challenge. A method is developed to match the distributions of galaxy fluxes and half-light radii from the deep fields to the corresponding distributions in the field of interest. We recover angle-only shear estimates with a performance close to current well-established model and moments-based methods for all three angle measurement techniques. The Q-values for all three methods are found to be Q~400. The code is freely available online at http://www.jb.man.ac.uk/~mbrown/angle_only_shear/.

preprint2015arXiv

Classically isospinning Skyrmion solutions

We investigate how isospin affects the geometrical shape and energy of classical soliton solutions of topological charges $B=1-4,8$ in the Skyrme model. The novel approach in our work is that we study classically isospinning Skyrmions beyond the rigid-body approximation; that is, we explicitly allow the soliton solutions to deform and to break the symmetries of the static configurations. Our fully three-dimensional relaxation calculations reveal that the symmetries of isospinning Skyrme solitons can differ significantly from the ones of the static configurations. In particular, isospinning Skyrmion solutions can break up into lower-charge Skyrmions, can deform into new solution types that do not exist at vanishing angular frequency $ω$ or energy degeneracy can be removed. These types of deformations have been largely ignored in previous work on modeling nuclei by quantized Skyrmion solutions.

preprint2015arXiv

Constraining dark sector perturbations I: cosmic shear and CMB lensing

We present current and future constraints on equations of state for dark sector perturbations. The equations of state considered are those corresponding to a generalized scalar field model and time-diffeomorphism invariant $L(g)$ theories that are equivalent to models of a relativistic elastic medium and also Lorentz violating massive gravity. We develop a theoretical understanding of the observable impact of these models. In order to constrain these models we use CMB temperature data from Planck, BAO measurements, CMB lensing data from Planck and the South Pole Telescope, and weak galaxy lensing data from CFHTLenS. We find non-trivial exclusions on the range of parameters, although the data remains compatible with $w=-1$. We gauge how future experiments will help to constrain the parameters. This is done via a likelihood analysis for CMB experiments such as CoRE and PRISM, and tomographic galaxy weak lensing surveys, focussing in on the potential discriminatory power of Euclid on mildly non-linear scales.

preprint2015arXiv

Constraining dark sector perturbations II: ISW and CMB lensing tomography

Any Dark Energy (DE) or Modified Gravity (MG) model that deviates from a cosmological constant requires a consistent treatment of its perturbations, which can be described in terms of an effective entropy perturbation and an anisotropic stress. We have considered a recently proposed generic parameterisation of DE/MG perturbations and compared it to data from the Planck satellite and six galaxy catalogues, including temperature-galaxy (Tg), CMB lensing-galaxy and galaxy-galaxy (gg) correlations. Combining these observables of structure formation with tests of the background expansion allows us to investigate the properties of DE/MG both at the background and the perturbative level. Our constraints on DE/MG are mostly in agreement with the cosmological constant paradigm, while we also find that the constraint on the equation of state w (assumed to be constant) depends on the model assumed for the perturbation evolution. We obtain $w=-0.92^{+0.20}_{-0.16}$ (95% CL; CMB+gg+Tg) in the entropy perturbation scenario; in the anisotropic stress case the result is $w=-0.86^{+0.17}_{-0.16}$. Including the lensing correlations shifts the results towards higher values of w. If we include a prior on the expansion history from recent Baryon Acoustic Oscillations (BAO) measurements, we find that the constraints tighten closely around $w=-1$, making it impossible to measure any DE/MG perturbation evolution parameters. If, however, upcoming observations from surveys like DES, Euclid or LSST show indications for a deviation from a cosmological constant, our formalism will be a useful tool towards model selection in the dark sector.

preprint2015arXiv

HI intensity mapping with FAST

We discuss the detectability of large-scale HI intensity fluctuations using the FAST telescope. We present forecasts for the accuracy of measuring the Baryonic Acoustic Oscillations and constraining the properties of dark energy. The FAST $19$-beam L-band receivers ($1.05$--$1.45$ GHz) can provide constraints on the matter power spectrum and dark energy equation of state parameters ($w_{0},w_{a}$) that are comparable to the BINGO and CHIME experiments. For one year of integration time we find that the optimal survey area is $6000\,{\rm deg}^2$. However, observing with larger frequency coverage at higher redshift ($0.95$--$1.35$ GHz) improves the projected errorbars on the HI power spectrum by more than $2~σ$ confidence level. The combined constraints from FAST, CHIME, BINGO and Planck CMB observations can provide reliable, stringent constraints on the dark energy equation of state.

preprint2015arXiv

Separating weak lensing and intrinsic alignments using radio observations

We discuss methods for performing weak lensing using radio observations to recover information about the intrinsic structural properties of the source galaxies. Radio surveys provide unique information that can benefit weak lensing studies, such as HI emission, which may be used to construct galaxy velocity maps, and polarized synchrotron radiation; both of which provide information about the unlensed galaxy and can be used to reduce galaxy shape noise and the contribution of intrinsic alignments. Using a proxy for the intrinsic position angle of an observed galaxy, we develop techniques for cleanly separating weak gravitational lensing signals from intrinsic alignment contamination in forthcoming radio surveys. Random errors on the intrinsic orientation estimates introduce biases into the shear and intrinsic alignment estimates. However, we show that these biases can be corrected for if the error distribution is accurately known. We demonstrate our methods using simulations, where we reconstruct the shear and intrinsic alignment auto and cross-power spectra in three overlapping redshift bins. We find that the intrinsic position angle information can be used to successfully reconstruct both the lensing and intrinsic alignment power spectra with negligible residual bias.

preprint2014arXiv

Computing model independent perturbations in dark energy and modified gravity

We present a methodology for computing model independent perturbations in dark energy and modified gravity. This is done from the Lagrangian for perturbations, by showing how field content, symmetries, and physical principles are often sufficient ingredients for closing the set of perturbed fluid equations. The fluid equations close once "equations of state for perturbations" are identified: these are linear combinations of fluid and metric perturbations which construct gauge invariant entropy and anisotropic stress perturbations for broad classes of theories. Our main results are the proof of the equation of state for perturbations presented in a previous paper, and the development of the required calculational tools.

preprint2014arXiv

Evidence for massive neutrinos from CMB and lensing observations

We discuss whether massive neutrinos (either active or sterile) can reconcile some of the tensions within cosmological data that have been brought into focus by the recently released {\it Planck} data. We point out that a discrepancy is present when comparing the primary CMB and lensing measurements both from the CMB and galaxy lensing data using CFHTLenS, similar to that which arises when comparing CMB measurements and SZ cluster counts. A consistent picture emerges and including a prior for the cluster constraints and BAOs we find that: for an active neutrino model with 3 degenerate neutrinos, $\sum m_ν= (0.320 \pm 0.081)\,{\rm eV}$, whereas for a sterile neutrino, in addition to 3 neutrinos with a standard hierarchy and $\sum m_ν= 0.06\,{\rm eV}$, $m_{ν, \, \rm sterile}^{\rm eff}= (0.450 \pm 0.124)\,{\rm eV}$ and $ΔN_{\rm eff} = 0.45 \pm 0.23$. In both cases there is a significant detection of modification to the neutrino sector from the standard model and in the case of the sterile neutrino it is possible to reconcile the BAO and local $H_0$ measurements. However, a caveat to our result is some internal tension between the CMB and lensing/cluster observations, and the masses are in excess of those estimated from the shape of the matter power spectrum from galaxy surveys.

preprint2014arXiv

Impact of baryons on the cluster mass function and cosmological parameter determination

Recent results by the Planck collaboration have shown that cosmological parameters derived from the cosmic microwave background anisotropies and cluster number counts are in tension, with the latter preferring lower values of the matter density parameter, $Ω_\mathrm{m}$, and power spectrum amplitude, $σ_8$. Motivated by this, we investigate the extent to which the tension may be ameliorated once the effect of baryonic depletion on the cluster mass function is taken into account. We use the large-volume Millennium Gas simulations in our study, including one where the gas is pre-heated at high redshift and one where the gas is heated by stars and active galactic nuclei (in the latter, the self-gravity of the baryons and radiative cooling are omitted). In both cases, the cluster baryon fractions are in reasonably good agreement with the data at low redshift, showing significant depletion of baryons with respect to the cosmic mean. As a result, it is found that the cluster abundance in these simulations is around 15 per cent lower than the commonly-adopted fit to dark matter simulations by Tinker et al (2008) for the mass range $10^{14}-10^{14.5}h^{-1} \mathrm{M}_\odot$. Ignoring this effect produces a significant artificial shift in cosmological parameters which can be expressed as $Δ[σ_8(Ω_\mathrm{m}/0.27)^{0.38}]\simeq -0.03$ at $z=0.17$ (the median redshift of the $\mathit{Planck}$ cluster sample) for the feedback model. While this shift is not sufficient to fully explain the $\mathit{Planck}$ discrepancy, it is clear that such an effect cannot be ignored in future precision measurements of cosmological parameters with clusters. Finally, we outline a simple, model-independent procedure that attempts to correct for the effect of baryonic depletion and show that it works if the baryon-dark matter back-reaction is negligible.

preprint2014arXiv

Isospinning baby Skyrmion solutions

We perform full two-dimensional (2D) numerical relaxations of isospinning soliton solutions in the baby Skyrme model in which the global $O(3)$ symmetry is broken by the 2D analogue of the pion mass term in the Skyrme model. In our calculations we explicitely allow the isospinning solitons to deform and to break the symmetries of the static configurations. We find that stable isospinning baby Skyrme solutions can be constructed numerically for all angular frequencies $ω\le \text{min}(μ,1)$, where $μ$ is the mass parameter of the model. Stable, rotationally-symmetric baby Skyrmion solutions for higher angular velocities are simply an artefact of the hedgehog approximation. Isospinning multisoliton solutions of topological charge $B$ turn out to be unstable to break up into their $B$ charge-1 constituents at some critical breakup frequency value. Furthermore, we find that for $μ$ sufficiently large the rotational symmetry of charge-2 baby Skyrmions becomes broken at a critical angular frequency $ω$.

preprint2014arXiv

Parameterizing dark sector perturbations via equations of state

The evolution of perturbations is a crucial part of the phenomenology of the dark sector cosmology. We advocate parameterizing these perturbations using equations of state for the entropy perturbation and the anisotropic stress. For small perturbations, these equations of state will be linear in the density, velocity and metric perturbations, and in principle these can be related back to the field content of the underlying model allowing for confrontation with observations. We illustrate our point by constructing gauge invariant entropy perturbations for theories where the dark sector Lagrangian is a general function of a scalar field, its first and second derivatives, and the metric and its first derivative, ${\cal L}={\cal L}(ϕ,\partial_μϕ,\partial_μ\partial_νϕ,g_{μν},\partial_αg_{μν})$. As an example, we show how to apply this approach to the case of models of Kinetic Gravity Braiding.

preprint2014arXiv

Tension between the power spectrum of density perturbations measured on large and small scales

There is a tension between measurements of the amplitude of the power spectrum of density perturbations inferred using the Cosmic Microwave Background (CMB) and directly measured by Large-Scale Structure (LSS) on smaller scales. We show that this tension exists, and is robust, for a range of LSS indicators including clusters, lensing and redshift space distortions and using CMB data from either $Planck$ or WMAP+SPT/ACT. One obvious way to try to reconcile this is the inclusion of a massive neutrino which could be either active or sterile. Using $Planck$ and a combination of all the LSS data we find that (i) for an active neutrino $\sum m_ν= (0.357\pm0.099)\,{\rm eV}$ and (ii) for a sterile neutrino $m_{\rm sterile}^{\rm eff}= (0.67\pm0.18)\,{\rm eV}$ and $ΔN_{\rm eff}= 0.32\pm0.20$. This is, however, at the expense of a degraded fit to $Planck$ temperature data, and we quantify the residual tension at $2.5σ$ and $1.6 σ$ for massive and sterile neutrinos respectively. We also consider alternative explanations including a lower redshift for reionization that would be in conflict with polarisation measurements made by WMAP and $ad$-$hoc$ modifications to primordial power spectrum.

preprint2013arXiv

Classically Isospinning Hopf Solitons

We perform full three-dimensional numerical relaxations of isospinning Hopf solitons with Hopf charge up to 8 in the Skyrme-Faddeev model with mass terms included. We explicitly allow the soliton solution to deform and to break the symmetries of the static configuration. It turns out that the model with its rich spectrum of soliton solutions, often of similar energy, allows for transmutations, formation of new solution types and the rearrangement of the spectrum of minimal-energy solitons in a given topological sector when isospin is added. We observe that the shape of isospinning Hopf solitons can differ qualitatively from that of the static solution. In particular the solution type of the lowest energy soliton can change. Our numerical results are of relevance for the quantization of the classical soliton solutions.

preprint2013arXiv

Massive gravity, the elasticity of space-time and perturbations in the dark sector

We consider a class of modified gravity models where the terms added to the standard Einstein-Hilbert Lagrangian are just a function of the metric only. For linearized perturbations around an isotropic space-time, this class of models is entirely specified by a rank 4 tensor that encodes possibly time-dependent masses for the gravitons. This tensor has the same symmetries as an elasticity tensor, suggesting an interpretation of massive gravity as an effective rigidity of space-time. If we choose a form for this tensor which is compatible with the symmetries of FRW and enforce full reparameterization invariance, then the only theory possible is a cosmological constant. However, in the case where the theory is only time translation invariant, the ghost-free massive gravity theory is equivalent to the elastic dark energy scenario with the extra Lorentz violating vector giving rise to 2 transverse and 1 longitudinal degrees of freedom, whereas when one demands spatial translation invariance one is left with scalar field theory with a non-standard kinetic term.

preprint2013arXiv

Projected constraints on the cosmic (super)string tension with future gravitational wave detection experiments

We present projected constraints on the cosmic string tension, $Gμ/c^2$, that could be achieved by future gravitational wave detection experiments and express our results as semi-analytic relations of the form $Gμ(Ω_{\rm gw}h^2)/c^2$, to allow for direct computation of the tension constraints for future experiments. These results can be applied to new constraints on $\ogwh$ as they are imposed. Experiments operating in different frequency bands probe different parts of the gravitational wave spectrum of a cosmic string network and are sensitive to different uncertainties in the underlying cosmic string model parameters. We compute the gravitational wave spectra of cosmic string networks based on the one-scale model, covering all the parameter space accessed by each experiment which is strongly dependent on the birth scale of loops relative to the horizon, $α$. The upper limits on the string tension avoid any assumptions on the model parameters. We perform this investigation for Pulsar Timing Array experiments of different durations as well as ground-based and space-borne interferometric detectors.

preprint2012arXiv

Effective action approach to cosmological perturbations in dark energy and modified gravity

In light of upcoming observations modelling perturbations in dark energy and modified gravity models has become an important topic of research. We develop an effective action to construct the components of the perturbed dark energy momentum tensor which appears in the perturbed generalized gravitational field equations, δG_{μν} = 8πGδT_{μν} + δU_{μν} for linearized perturbations. Our method does not require knowledge of the Lagrangian density of the dark sector to be provided, only its field content. The method is based on the fact that it is only necessary to specify the perturbed Lagrangian to quadratic order and couples this with the assumption of global statistical isotropy of spatial sections to show that the model can be specified completely in terms of a finite number of background dependent functions. We present our formalism in a coordinate independent fashion and provide explicit formulae for the perturbed conservation equation and the components of δU_{μν} for two explicit generic examples: (i) the dark sector does not contain extra fields, L = L(g_{μν}) and (ii) the dark sector contains a scalar field and its first derivative L = L(g_{μν}, ϕ, \nabla_μϕ). We discuss how the formalism can be applied to modified gravity models containing derivatives of the metric, curvature tensors, higher derivatives of the scalar fields and vector fields.

preprint2012arXiv

Sunyaev-Zel'dovich clusters in Millennium Gas simulations

We have exploited the large-volume Millennium Gas cosmological N-body hydrodynamics simulations to study the SZ cluster population at low and high redshift, for three models with varying gas physics. We confirm previous results using smaller samples that the intrinsic (spherical) Y_{500}-M_{500} relation has very little scatter (sigma_{log_{10}Y}~0.04), is insensitive to cluster gas physics and evolves to redshift one in accord with self-similar expectations. Our pre-heating and feedback models predict scaling relations that are in excellent agreement with the recent analysis from combined Planck and XMM-Newton data by the Planck Collaboration. This agreement is largely preserved when r_{500} and M_{500} are derived using the hydrostatic mass proxy, Y_{X,500}, albeit with significantly reduced scatter (sigma_{log_{10}Y}~0.02), a result that is due to the tight correlation between Y_{500} and Y_{X,500}. Interestingly, this assumption also hides any bias in the relation due to dynamical activity. We also assess the importance of projection effects from large-scale structure along the line-of-sight, by extracting cluster Y_{500} values from fifty simulated 5x5 square degree sky maps. Once the (model-dependent) mean signal is subtracted from the maps we find that the integrated SZ signal is unbiased with respect to the underlying clusters, although the scatter in the (cylindrical) Y_{500}-M_{500} relation increases in the pre-heating case, where a significant amount of energy was injected into the intergalactic medium at high redshift. Finally, we study the hot gas pressure profiles to investigate the origin of the SZ signal and find that the largest contribution comes from radii close to r_{500} in all cases. The profiles themselves are well described by generalised Navarro, Frenk & White profiles but there is significant cluster-to-cluster scatter.

preprint2012arXiv

X-type and Y-type junction stability in domain wall networks

We develop an analytic formalism that allows one to quantify the stability properties of X-type and Y-type junctions in domain wall networks in two dimensions. A similar approach might be applicable to more general defect systems involving junctions that appear in a range of physical situations, for example, in the context of F- and D-type strings in string theory. We apply this formalism to a particular field theory, Carter's pentavac model, where the strength of the symmetry breaking is governed by the parameter $|ε|< 1$. We find that for low values of the symmetry breaking parameter X-type junctions will be stable, whereas for higher values an X-type junction will separate into two Y-type junctions. The critical angle separating the two regimes is given by α_c = 293^{\circ}\sqrt{|ε|}$ and this is confirmed using simple numerical experiments. We go on to simulate the pentavac model from random initial conditions and we find that the dominant junction is of \ytype for $|ε| \geq 0.02$ and is of \xtype for $|ε| \leq 0.02$. We also find that for small $ε$ the evolution of the number of domain walls $\qsubrm{N}{dw}$ in Minkowski space does not follow the standard $\propto t^{-1}$ scaling law with the deviation from the standard lore being more pronounced as $ε$ is decreased. The presence of dissipation appears to restore the $t^{-1}$ lore.

preprint2011arXiv

Mapping the dark matter with polarized radio surveys

In a recent paper (Brown & Battye 2011), we proposed the use of integrated polarization measurements of background galaxies in radio weak gravitational lensing surveys and investigated the potential impact on the statistical measurement of cosmic shear. Here we extend this idea to reconstruct maps of the projected dark matter distribution, or lensing convergence field. The addition of polarization can, in principle, greatly reduce shape noise due to the intrinsic dispersion in galaxy ellipticities. We show that maps reconstructed using this technique in the radio band can be competitive with those derived using standard lensing techniques which make use of many more galaxies. In addition, since the reconstruction noise is uncorrelated between these standard techniques and the polarization technique, their comparison can serve as a powerful check for systematics and their combination can reduce noise further. We examine the convergence reconstruction which could be achieved with two forthcoming facilities: (i) a deep survey, covering 1.75 square degrees using the e-MERLIN instrument currently being commissioned in the UK and (ii) the high resolution, deep wide field surveys which will eventually be conducted with the Square Kilometre Array.

preprint2011arXiv

Modelling neutral hydrogen in galaxies using cosmological hydrodynamical simulations

The characterisation of the atomic and molecular hydrogen content of high-redshift galaxies is a major observational challenge that will be addressed over the coming years with a new generation of radio telescopes. We investigate this important issue by considering the states of hydrogen across a range of structures within high-resolution cosmological hydrodynamical simulations. Additionally, our simulations allow us to investigate the sensitivity of our results to numerical resolution and to sub-grid baryonic physics (especially feedback from supernovae and active galactic nuclei). We find that the most significant uncertainty in modelling the neutral hydrogen distribution arises from our need to model a self-shielding correction in moderate density regions. Future simulations incorporating radiative transfer schemes will be vital to improve on our empirical self-shielding threshold. Irrespective of the exact nature of the threshold we find that while the atomic hydrogen mass function evolves only mildly from redshift two to zero, the molecular hydrogen mass function increases with increasing redshift, especially at the high-mass end. Interestingly, the weak evolution of the neutral hydrogen mass function is insensitive to the feedback scheme utilised, but the opposite is true for the molecular gas, which is more closely associated with the star formation in the simulations.

preprint2011arXiv

Vacuum Topology of the Two Higgs Doublet Model

We perform a systematic study of generic accidental Higgs-family and CP symmetries that could occur in the two-Higgs-doublet-model potential, based on a Majorana scalar-field formalism which realizes a subgroup of GL(8,C). We derive the general conditions of convexity and stability of the scalar potential and present analytical solutions for two non-zero neutral vacuum expectation values of the Higgs doublets for a typical set of six symmetries, in terms of the gauge-invariant parameters of the theory. By means of a homotopy-group analysis, we identify the topological defects associated with the spontaneous symmetry breaking of each symmetry, as well as the massless Goldstone bosons emerging from the breaking of the continuous symmetries. We find the existence of domain walls from the breaking of Z_2, CP1 and CP2 discrete symmetries, vortices in models with broken U(1)_PQ and CP3 symmetries and a global monopole in the SO(3)_HF-broken model. The spatial profile of the topological defect solutions is studied in detail, as functions of the potential parameters of the two-Higgs doublet model. The application of our Majorana scalar-field formalism in studying more general scalar potentials that are not constrained by the U(1)_Y hypercharge symmetry is discussed. In particular, the same formalism may be used to properly identify seven additional symmetries that may take place in a U(1)_Y-invariant scalar potential.

preprint2010arXiv

Charge, junctions and the scaling dynamics of domain wall networks

It has been shown that superconducting domain walls in a model with U(1) x Z2 symmetry can form long-lived loops called kinky vortons from random initial conditions in the broken field and a uniform charged background in (2+1) dimensions. In this paper we investigate a similar model with a hyper-cubic symmetry coupled to an unbroken U(1) in which the domain walls can form junctions and hence a lattice. We call this model the charge-coupled cubic-anisotropy (CCCA) model. First, we present a detailed parametric study of the U(1) x Z2 model; features which we vary include the nature of the initial conditions and the coupling constants. This allows us to identify interesting parameters to vary in the more complicated, and hence more computationally intensive, CCCA models. In particular we find that the coefficient of the interaction term can be used to engineer three separate regimes: phase mixing, condensation and phase separation with the condensation regime corresponding to a single value of the coupling constant defined by the determinant of the quartic interaction terms being zero. We then identify the condensation regime in the CCCA model and show that, in this regime, the number of domain walls does not scale in the standard way if the initial conditions have a sufficiently high background charge. Instead of forming loops of domain wall, we find that, within the constraints of dynamic range, the network appears to be moving toward a glass-like configuration. We find that the results are independent of the dimension of the hyper-cube.

preprint2010arXiv

Impact of baryon physics on dark matter structures: a detailed simulation study of halo density profiles

The back-reaction of baryons on the dark matter halo density profile is of great interest, not least because it is an important systematic uncertainty when attempting to detect the dark matter. Here, we draw on a large suite of high resolution cosmological hydrodynamical simulations, to systematically investigate this process and its dependence on the baryonic physics associated with galaxy formation. The inclusion of baryons results in significantly more concentrated density profiles if radiative cooling is efficient and feedback is weak. The dark matter halo concentration can in that case increase by as much as 30 (10) per cent on galaxy (cluster) scales. The most significant effects occur in galaxies at high redshift, where there is a strong anti-correlation between the baryon fraction in the halo centre and the inner slope of both the total and the dark matter density profiles. If feedback is weak, isothermal inner profiles form, in agreement with observations of massive, early-type galaxies. However, we find that AGN feedback, or extremely efficient feedback from massive stars, is necessary to match observed stellar fractions in groups and clusters, as well as to keep the maximum circular velocity similar to the virial velocity as observed for disk galaxies. These strong feedback models reduce the baryon fraction in galaxies by a factor of 3 relative to the case with no feedback. The AGN is even capable of reducing the baryon fraction by a factor of 2 in the inner region of group and cluster haloes. This in turn results in inner density profiles which are typically shallower than isothermal and the halo concentrations tend to be lower than in the absence of baryons.

preprint2010arXiv

Polarization as an indicator of intrinsic alignment in radio weak lensing

We propose a new technique for weak gravitational lensing in the radio band making use of polarization information. Since the orientation of a galaxy's polarized emission is both unaffected by lensing and is related to the galaxy's intrinsic orientation, it effectively provides information on the unlensed galaxy position angle. We derive a new weak lensing estimator which exploits this effect and makes full use of both the observed galaxy shapes and the estimates of the intrinsic position angles as provided by polarization. Our method has the potential to both reduce the effects of shot noise, and to reduce to negligible levels, in a model-independent way, all effects of intrinsic galaxy alignments. We test our technique on simulated weak lensing skies, including an intrinsic alignment contaminant consistent with recent observations, in three overlapping redshift bins. Adopting a standard weak lensing analysis and ignoring intrinsic alignments results in biases of 5-10% in the recovered power spectra and cosmological parameters. Applying our new estimator to one tenth the number of galaxies used for the standard case, we recover both power spectra and the input cosmology with similar precision as compared to the standard case and with negligible residual bias, even in the presence of a substantial (astrophysical) scatter in the relationship between the observed orientation of the polarized emission and the intrinsic orientation. Assuming a reasonable polarization fraction for star-forming galaxies, and no cosmological conspiracy in the relationship between polarization direction and intrinsic morphology, our estimator should prove a valuable tool for weak lensing analyses of forthcoming radio surveys, in particular, deep wide field surveys with e-MERLIN, MeerKAT and ASKAP and ultimately, definitive radio lensing surveys with the SKA.

preprint2009arXiv

Statistics of the Sunyaev-Zel'dovich Effect power spectrum

Using large numbers of simulations of the microwave sky, incorporating the Cosmic Microwave Background (CMB) and the Sunyaev-Zel'dovich (SZ) effect due to clusters, we investigate the statistics of the power spectrum at microwave frequencies between spherical multipoles of 1000 and 10000. From these virtual sky maps, we find that the spectrum of the SZ effect has a larger standard deviation by a factor of 3 than would be expected from purely Gaussian realizations, and has a distribution that is significantly skewed towards higher values, especially when small map sizes are used. The standard deviation is also increased by around 10 percent compared to the trispectrum calculation due to the clustering of galaxy clusters. We also consider the effects of including residual point sources and uncertainties in the gas physics. This has implications for the excess power measured in the CMB power spectrum by the Cosmic Background Imager and BIMA experiments. Our results indicate that the observed excess could be explained using a lower value of $σ_8$ than previously suggested, however the effect is not enough to match $σ_8=0.825$. The uncertainties in the gas physics could also play a substantial role. We have made our maps of the SZ effect available online.

preprint2008arXiv

Observations of the Corona Borealis supercluster with the superextended Very Small Array: further constraints on the nature of the non-Gaussian CMB cold spot

We present interferometric imaging at 33 GHz, with the new superextended configuration of the Very Small Array (VSA), of a very deep decrement in the cosmic microwave background (CMB) temperature. This decrement is located in the direction of the Corona Borealis supercluster, at a position with no known galaxy clusters, and was discovered by a previous VSA survey (Genova-Santos et al.). A total area of 3 sq.deg. has now been imaged, with an angular resolution of 7 arcmin and a flux sensitivity of 5 mJy/beam. These observations confirm the presence of this strong and resolved negative spot at -41+/-5 mJy/beam (-258+/-29 muK), with a signal to noise level of 8. This structure is also present in the WMAP 5-year data. The temperature of the W-band (94 GHz) data at the position of the decrement agrees within 1.2-sigma with that observed by the VSA at 33 GHz, and within 0.2-sigma with the Sunyaev-Zel'dovich (SZ) spectrum. Our analyses show that it is a non-Gaussian feature in the CMB at a level of 4.8-sigma. The probability of finding such a deviation or larger in CMB Gaussian simulations is only 0.19 per cent. Therefore, an explanation other than primordial CMB is required. We have considered the possibility of an SZ effect generated in a diffuse, extended warm/hot gas distribution. This hypothesis is especially relevant, as the presence of such structures, if confirmed, could provide the location for a significant fraction of the missing baryons in the Local Universe.

preprint2007arXiv

Do consistent $F(R)$ models mimic General Relativity plus $Λ$?

Modified gravity models are subject to a number of consistency requirements which restrict the form that the function $F(R)$ can take. We study a particular class of $F(R)$ functions which satisfy various constraints that have been found in the literature. These models have a late time accelerating epoch, and an acceptable matter era. We calculate the Friedmann equation for our models, and show that in order to satisfy the constraints we impose, they must mimic General Relativity plus $Λ$ throughout the cosmic history, with exponentially suppressed corrections. We also find that the free parameters in our model must be fine tuned to obtain an acceptable late time accelerating phase. We discuss the generality of this conclusion.

preprint2005arXiv

Constraints on the solid dark universe model

If the dark energy is modelled as a relativistic elastic solid then the standard CDM and $Λ$CDM models, as well as lattice configurations of cosmic strings or domain walls, are points in the two-dimensional parameter space $(w,c_{\rm s}^2)$. We present a detailed analysis of the best fitting cosmological parameters in this model using data from a range of observations. We find that the $χ^2$ is improved by $\sim 10$ by including the two parameters and that the $w=-1$ $Λ$CDM model is only the best fit to the data when a large number of different datasets are included. Using CMB observations alone we find that $w=-0.38\pm 0.16 $ and with the addition of Large-Scale Structure data $w=-0.62\pm 0.15 $ and $\log c_{\rm s}=-0.77\pm 0.28$. We conclude that the models based on topological defects provide a good fit to the current data, although $Λ$CDM cannot be ruled out.

preprint2004arXiv

Multiple-scales analysis of cosmological perturbations in brane-worlds

We present a new approximation method for solving the equations of motion for cosmological tensor perturbations in a Randall-Sundrum brane-world model of the type with one brane in a five-dimensional anti-de Sitter spacetime. This method avoids the problem of coordinate singularities inherent in some methods. At leading order, the zero-mode solution replicates the evolution of perturbations in a four-dimensional Friedmann-Robertson-Walker universe in the absence of any tensor component to the matter perturbation on the brane. At next order, there is a mode-mixing effect, although, importantly, the zero-mode does not source any other modes.

preprint2001arXiv

Einstein equations for an asymmetric brane-world

We consider a brane-world of co-dimension one without the reflection symmetry that is commonly imposed between the two sides of the brane. Using the coordinate-free formalism of the Gauss-Codacci equations, we derive the effective Einstein equations by relating the local curvature to the matter on the brane in the case when its bare tension is much larger than the localized matter, and hence show that Einstein gravity is a natural consequence of such models in the weak field limit. We find agreement with the recently derived cosmological case, which can be solved exactly, and point out that such models can be realized naturally in the case where there is a minimally coupled form field in the bulk.

preprint1999arXiv

Reionization by active sources and its effects on the cosmic microwave background

We investigate the possible effects of reionization by active sources on the cosmic microwave background. We concentrate on the sources themselves as the origin of reionization, rather than early object formation, introducing an extra period of heating motivated by the active character of the perturbations. Using reasonable parameters, this leads to four possibilities depending on the time and duration of the energy input: delayed last scattering, double last scattering, shifted last scattering and total reionization. We show that these possibilities are only very weakly constrained by the limits on spectral distortions from the COBE FIRAS measurements. We illustrate the effects of these reionization possibilities on the angular power spectrum of temperature anisotropies and polarization for simple passive isocurvature models and simple coherent sources, observing the difference between passive and active models. Finally, we comment on the implications of this work for more realistic active sources, such as causal white noise and topological defect models. We show for these models that non-standard ionization histories can shift the peak in the CMB power to larger angular scales.

preprint1997arXiv

The case against scaling defect models of cosmic structure formation

We calculate predictions from defect models of structure formation for both the matter and Cosmic Microwave Background (CMB) over all observable scales. Our results point to a serious problem reconciling the observed large-scale galaxy distribution with the COBE normalization, a result which is robust for a wide range of defect parameters. We conclude that standard scaling defect models are in conflict with the data, and show how attempts to resolve the problem by considering non-scaling defects would require radical departures from the standard scaling picture.