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.

We present a significant update of the constraints on the Abelian Higgs cosmic string tension by cosmic microwave background (CMB) data, enabled both by the use of new high-resolution CMB data from suborbital experiments as well as the latest results of the WMAP satellite, and by improved predictions for the impact of Abelian Higgs cosmic strings on the CMB power spectra. The new cosmic string spectra (presented in a previous work) were improved especially for small angular scales, through the use of larger Abelian Higgs string simulations and careful extrapolation. If Abelian Higgs strings are present then we find improved bounds on their contribution to the CMB anisotropies, f10< 0.095, and on their tension, G\mu< 0.57 x 10^-6, both at 95% confidence level using WMAP7 data; and f10 < 0.048 and G\mu < 0.42 x 10^-6 using all the CMB data. We also find that using all the CMB data, a scale invariant initial perturbation spectrum, ns=1, is now disfavoured at 2.4\sigma\ even if strings are present. A Bayesian model selection analysis no longer indicates a preference for strings.

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{\sigma} level. We also confirm that a tensor contribution of r = 0.0018 can be detected at over 3{\sigma}, in agreement with the CMBpol mission concept study. These results are supported by a model selection analysis.

In previous work it was shown that MOND theories predict anomalously strong tidal stresses near the saddle points of the Newtonian gravitational potential. An analytical examination of the saddle between two bodies revealed a linear and a non-linear solution, valid for the outer and inner regions. Here we present a numerical algorithm for solving the MOND equations. We check the code against the two-body analytical solutions and explore the region transitioning between them. We then develop a a realistic model for the MONDian effects on the saddles of the Sun-Earth-Moon system (including further sources is straightforward). For the Sun-Earth saddle we find that the two-body results are almost unchanged, with corrections increasing from full to new Moon. In contrast, the Moon saddle is an intrinsically three-body problem, but we numerically find a recipe for adapting the two-body solution to this case, by means of a suitable re-scaling and axis re-orientation. We explore possible experimental scenarios for LISA Pathfinder, and the prospect of a visit to the saddle(s) at the end of the mission. Given the chaotic nature of the orbits, awareness of the full range of the possibilities is crucial for a realistic prediction. We conclude that even with very conservative assumptions on the impact parameter, the accelerometers are abundantly sensitive to vindicate or rule out the theory.

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.

We suggest that LISA Pathfinder could be used to subject TEVES, and in particular the non-relativistic MOND phenomenology it incorporates, to a direct, controlled experimental test, in just a few years' time. The basic concept is to fly LISA Pathfinder through the region around the Sun-Earth saddle point, following its nominal mission, in order to look for anomalous gravity gradients. We examine various strategies to reach the saddle point, and conclude that the preferred strategy, resulting in relatively short transfer times of order one year, probably involves a lunar fly-by. We present robust estimates of the MOND gravity gradients that LISA Pathfinder should be exposed to, and conclude that if the gradiometer on-board the spacecraft achieves its nominal performance, these gradients will not just be detected, but measured and characterised in some detail, should they exist. Conversely, given the large predicted signal based on standard assumptions, a null result would most likely spell the end of TEVES/MOND.

We present the first ever calculation of cosmic microwave background CMB anisotropy power spectra from semilocal cosmic strings, obtained via simulations of a classical field theory. Semilocal strings are a type of non-topological defect arising in some models of inflation motivated by fundamental physics, and are thought to relax the constraints on the symmetry breaking scale as compared to models with (topological) cosmic strings. We derive constraints on the model parameters, including the string tension parameter mu, from fits to cosmological data, and find that in this regard BPS semilocal strings resemble global textures more than topological strings. The observed microwave anisotropy at l = 10 is reproduced if Gmu = 5.3x10^{-6} (G is Newton's constant). However as with other defects the spectral shape does not match observations, and in models with inflationary perturbations plus semilocal strings the 95% confidence level upper bound is Gmu<2.0x10^{-6} when CMB data, Hubble Key Project and Big Bang Nucleosynthesis data are used (c.f. Gmu<0.9x10^{-6} for cosmic strings). We additionally carry out a Bayesian model comparison of several models with and without defects, showing models with defects are neither conclusively favoured nor disfavoured at present.

While observations indicate that the predominant source of cosmic inhomogeneities are adiabatic perturbations, there are a variety of candidates to provide auxiliary trace effects, including inflation-generated primordial tensors and cosmic defects which both produce B-mode cosmic microwave background (CMB) polarization. We investigate whether future experiments may suffer confusion as to the true origin of such effects, focusing on the ability of Planck to distinguish tensors from cosmic strings, and show that there is no significant degeneracy.

We perform a multiparameter likelihood analysis to compare measurements of the cosmic microwave background (CMB) power spectra with predictions from models involving cosmic strings. Adding strings to the standard case of a primordial spectrum with power-law tilt n, we find a 2-sigma detection of strings: f_10 = 0.11 +/- 0.05, where f_10 is the fractional contribution made by strings in the temperature power spectrum (at multipole l = 10). CMB data give moderate preference to the model n = 1 with cosmic strings over the standard zero-strings model with variable tilt. When additional non-CMB data are incorporated, the two models become on a par. With variable n and these extra data, we find that f_10 < 0.11, which corresponds to G mu < 0.7x10^-6 (where mu is the string tension and G is the gravitational constant).

We present the first calculation of the possible (local) cosmic string contribution to the cosmic microwave background polarization spectra from simulations of a string network (rather than a stochastic collection of unconnected string segments). We use field theory simulations of the Abelian Higgs model to represent local U(1) strings, including their radiative decay and microphysics. Relative to previous estimates, our calculations show a shift in power to larger angular scales, making the chance of a future cosmic string detection from the B-mode polarization slightly greater. We explore a future ground-based polarization detector, taking the CLOVER project as our example. In the null hypothesis (that cosmic strings make a zero contribution) we find that CLOVER should limit the string tension \mu to G\mu<0.12e-6 (where G is the gravitational constant), above which it is likely that a detection would be possible.

We present the first field-theoretic calculations of the contribution made by cosmic strings to the temperature power spectrum of the cosmic microwave background (CMB). Unlike previous work, in which strings were modeled as idealized one-dimensional objects, we evolve the simplest example of an underlying field theory containing local U(1) strings, the Abelian Higgs model. Limitations imposed by finite computational volumes are overcome using the scaling property of string networks and a further extrapolation related to the lessening of the string width in comoving coordinates. The strings and their decay products, which are automatically included in the field theory approach, source metric perturbations via their energy-momentum tensor, the unequal-time correlation functions of which are used as input into the CMB calculation phase. These calculations involve the use of a modified version of CMBEASY, with results provided over the full range of relevant scales. We find that the string tension $\mu$ required to normalize to the WMAP 3-year data at multipole $\ell = 10$ is $G\mu = [2.04\pm0.06\textrm{(stat.)}\pm0.12\textrm{(sys.)}] \times 10^{-6}$, where we have quoted statistical and systematic errors separately, and $G$ is Newton's constant. This is a factor 2-3 higher than values in current circulation.

We present predictions for the angular correlation function and redshift distribution for SHADES, the SCUBA HAlf-Degree Extragalactic Survey, which will yield a sample of around 300 sub-mm sources in the 850 micron waveband in two separate fields. Complete and unbiased photometric redshift information on these sub-mm sources will be derived by combining the SCUBA data with i) deep radio imaging already obtained with the VLA, ii) guaranteed-time Spitzer data at mid-infrared wavelengths, and iii) far-infrared maps to be produced by BLAST, the Balloon-borne Large-Aperture Sub-millimeter Telescope. Predictions for the redshift distribution and clustering properties of the final anticipated SHADES sample have been computed for a wide variety of models, each constrained to fit the observed number counts. Since we are dealing with around 150 sources per field, we use the sky-averaged angular correlation function to produce a more robust fit of a power-law shape w(theta)=(theta/A)^{-delta} to the model data. Comparing the predicted distributions of redshift and of the clustering amplitude A and slope delta, we find that models can be constrained from the combined SHADES data with the expected photometric redshift information.

We use the cosmic microwave background angular power spectra to place upper limits on the degree to which global defects may have aided cosmic structure formation. We explore this under the inflationary paradigm, but with the addition of textures resulting from the breaking of a global O(4) symmetry during the early stages of the Universe. As a measure of their contribution, we use the fraction of the temperature power spectrum that is attributed to the defects at a multipole of 10. However, we find a parameter degeneracy enabling a fit to the first-year WMAP data to be made even with a significant defect fraction. This degeneracy involves the baryon fraction and the Hubble constant, plus the normalization and tilt of the primordial power spectrum. Hence, constraints on these cosmological parameters are weakened. Combining the WMAP data with a constraint on the physical baryon fraction from big bang nucleosynthesis calculations and high-redshift deuterium abundance, limits the extent of the degeneracy and gives an upper bound on the defect fraction of 0.13 (95% confidence).