The space-borne missions CoRoT and Kepler have revealed numerous mixed modes in red-giant stars. These modes carry a wealth of information about red-giant cores, but are of limited use when constraining rapid structural variations in their envelopes. This limitation can be circumvented if we have access to the frequencies of the pure acoustic dipolar modes in red giants, i.e. the dipole modes that would exist in the absence of coupling between gravity and acoustic waves. We present a pilot study aimed at evaluating the implications of using these pure acoustic mode frequencies in seismic studies of the helium structural variation in red giants. The study is based on artificial seismic data for a red-giant-branch stellar model, bracketing seven acoustic dipole radial orders around vmax. The pure acoustic dipole-mode frequencies are derived from a fit to the mixed-mode period spacings and then used to compute the pure acoustic dipole-mode second differences. The pure acoustic dipole-mode second differences inferred through this procedure follow the same oscillatory function as the radial modes second differences. The additional constraints brought by the dipolar modes allow us to adopt a more complete description of the glitch signature when performing the fit to the second differences. The amplitude of the glitch retrieved from this fit is 15% smaller than that from the fit based on the radial modes alone. Also, we find that thanks to the additional constraints, a bias in the inferred glitch location, found when adopting the simpler description of the glitch, is avoided.

In this paper we demonstrate that Boltzmann's $H$-theorem does not necessarily hold in the context of theories of gravity with nonminimally coupled matter fields. We find sufficient conditions for the violation of Boltzmann's $H$-theorem and derive an expression for the evolution of Boltzmann's $H$ in terms of the nonminimal coupling function, valid in the case of a collisionless gas in a homogeneous and isotropic Friedmann-Lema\^itre-Robertson-Walker universe. We highlight the implications of this result for the evolution of the entropy of the matter fields, briefly discussing the role played by collisions between particles whenever they are relevant. We also suggest a possible link between the high entropy of the Universe and the weakness of gravity in the context of these theories.

In this brief reply we respond to the note of Bertolami and Gomes (arXiv:2005.03968) on our recent paper (arXiv:2003.10154).

We investigate the appropriateness of the use of different Lagrangians to describe various components of the cosmic energy budget, discussing the degeneracies between them, in the absence of non-minimal couplings to gravity or other fields, and clarifying some misconceptions in the literature. We then show that models with the same on-shell Lagrangian may have different proper energy densities and use this result to map dark energy models into unified dark energy models in which dark matter and dark energy are described by the same perfect fluid. We determine the correspondence between their equation of state parameters and sound speeds, briefly discussing the linear sound speed problem of unified dark energy models as well as a possible way out associated to the non-linear dynamics.

We perform a detailed comparison between a recently proposed parameter-free velocity-dependent one-scale model and the standard parametric model for the cosmological evolution of domain wall networks. We find that the latter overestimates the damping of the wall motion due to the Hubble expansion and neglects the direct impact of wall decay on the evolution of the root-mean-square velocity of the network. We show that these effects are significant but may be absorbed into a redefinition of the momentum parameter. We also discuss the implications of these findings for cosmic strings. We compute the energy loss and momentum parameters of the standard parametric model for cosmological domain wall evolution using our non-parametric velocity-dependent one-scale model in the context of cosmological models having a power law evolution of the scale factor $a$ with the cosmic time $t$ ($a \propto t^\lambda$, $0 < \lambda < 1$), and compare with the results obtained from numerical field theory simulations. We further provide simple linear functions which roughly approximate the dependence of the energy loss and momentum parameters on $\lambda$.

In the present work we show that the second law of thermodynamics does not generally hold if the matter and gravitational fields are nonminimally coupled. We demonstrate this result by explicitly computing the evolution of the entropy of the matter fields in the case of a closed homogeneous and isotropic universe filled with dust and radiation, showing that, in this case, the sign of the entropy variation is determined by the evolution of the universe. The preservation of the second law of thermodynamics in these modified theories would require its generalization to account for a gravitational entropy contribution.

We derive a full analytical approximation to the stochastic gravitational wave background generated by the loops that are produced throughout the cosmological evolution of cosmic string networks. We show that this approximation not only predicts the amplitude of the radiation-era plateau exactly, but also provides a good fit to the high-frequency cut-off and to the low-frequency peak generated by the loops that decay during the matter era, irrespective of cosmic string tension and of the length of loops created. We then find that it provides a good quantitative description of the full stochastic gravitational wave background across the relevant frequency range.

We consider the impact of stochastic perturbations on otherwise coherent oscillations of classical pulsators. The resulting dynamics are modelled by a driven damped harmonic oscillator subject to either an external or an internal forcing and white noise velocity fluctuations. We characterize the phase and relative amplitude variations using analytical and numerical tools. When the forcing is internal the phase variation displays a random walk behaviour and a red noise power spectrum with a ragged erratic appearance. We determine the dependence of the root mean square phase and relative amplitude variations ($\sigma_{\Delta \varphi}$ and $\sigma_{\Delta A/A}$, respectively) on the amplitude of the stochastic perturbations, the damping constant $\eta$, and the total observation time $t_{\rm obs}$ for this case, under the assumption that the relative amplitude variations remain small, showing that $\sigma_{\Delta \varphi}$ increases with $t_{\rm obs}^{1/2}$ becoming much larger than $\sigma_{\Delta A/A}$ for $t_{\rm obs} \gg \eta^{-1}$. In the case of an external forcing the phase and relative amplitude variations remain of the same order, independent of the observing time. In the case of an internal forcing, we find that $\sigma_{\Delta \varphi}$ does not depend on $\eta$. Hence, the damping time cannot be inferred from fitting the power of the signal, as done for solar-like pulsators, but the amplitude of the stochastic perturbations may be constrained from the observations. Our results imply that, given sufficient time, the variation of the phase associated to the stochastic perturbations in internally driven classical pulsators will become sufficiently large to be probed observationally.

We develop a parameter-free velocity-dependent one-scale model for the evolution of the characteristic length $L$ and root-mean-square velocity $\sigma_v$ of standard domain wall networks in homogeneous and isotropic cosmologies. We compare the frictionless scaling solutions predicted by our model, in the context of cosmological models having a power law evolution of the scale factor $a$ as a function of the cosmic time $t$ ($a \propto t^\lambda$, $0< \lambda < 1$), with the corresponding results obtained using field theory numerical simulations. We show that they agree well (within a few $\%$) for root-mean-square velocities $\sigma_v$ smaller than $0.2 \, c$ ($\lambda \ge 0.9$), where $c$ is the speed of light in vacuum, but significant discrepancies occur for larger values of $\sigma_v$ (smaller values of $\lambda$). We identify problems with the determination of $L$ and $\sigma_v$ from numerical field theory simulations which might potentially be responsible for these discrepancies.

The characterisation of stellar cores may be accomplished through the modelling of asteroseismic data from stars exhibiting either gravity-mode or mixed-mode pulsations, potentially shedding light on the physical processes responsible for the production, mixing, and segregation of chemical elements. In this work we validate against model data an analytical expression for the period spacing that will facilitate the inference of the properties of stellar cores, including the detection and characterisation of buoyancy glitches (strong chemical gradients). This asymptotically-based analytical expression is tested both in models with and without buoyancy glitches. It does not assume that glitches are small and, consequently, predicts non-sinusoidal glitch-induced period-spacing variations, as often seen in model and real data. We show that the glitch position and width inferred from the fitting of the analytical expression to model data consisting of pure gravity modes are in close agreement (typically better than 7$\%$ relative difference) with the properties measured directly from the stellar models. In the case of fitting mixed-mode model data, the same expression is shown to reproduce well the numerical results, when the glitch properties are known a priori. In addition, the fits performed to mixed-mode model data reveal a frequency dependence of the coupling coefficient, $q$, for a moderate-luminosity red-giant-branch model star. Finally, we find that fitting the analytical expression to the mixed-mode period spacings may provide a way to infer the frequencies of the pure acoustic dipole modes that would exist if no coupling took place between acoustic and gravity waves.

Recently, a measurement of the pressure distribution experienced by the quarks inside the proton has found a strong repulsive (positive) pressure at distances up to 0.6 femtometers from its center and a (negative) confining pressure at larger distances. In this paper we show that this measurement puts significant constraints on modified theories of gravity in which the strength of the gravitational interaction on microscopic scales is enhanced with respect to general relativity. We consider the particular case of Eddington-inspired Born-Infeld gravity, showing that strong limits on $\kappa$, the only additional parameter of the theory with respect to general relativity, may be derived from the quark pressure measurement ($|\kappa| \lsim 10^{-1} \, {\rm m^5 \, kg^{-1} \, s^{-2}}$). Furthermore, we show how these limits may be significantly improved with precise measurements of the first and second moments of the pressure distribution inside the proton.

In extended models of gravity a nonminimal coupling to matter has been assumed to lead to irreversible particle creation. In this paper we challenge this assumption. We argue that a non-minimal coupling of the matter and gravitational sectors results in a change in particle-momentum on a cosmological timescale, irrespective of particle creation or decay. We further argue that particle creation or decay associated with a non-minimal coupling to gravity could only happen as a result of significant deviations from a homogeneous Friedmann-Robertson-Walker geometry on microscopic scales, and provide a phenomenological description of the impact of particle creation or decay on the cosmological evolution of the density of the matter fields.

A cosmic string network created during an inflationary stage in the early Universe - here defined as i-string network - is expected to enter a transient stretching regime during inflation, in which its characteristic length is stretched to scales much larger than the Hubble radius, before attaining a standard evolution once the network re-enters the Hubble volume after inflation. During the stretching regime, the production of cosmic string loops and consequent emission of gravitational radiation are significantly suppressed. Here, we compute the power spectrum of the stochastic gravitational wave background generated by i-string networks using the velocity-dependent one scale model to describe the network dynamics, and demonstrate that this regime introduces a high-frequency signature on an otherwise standard spectrum of the stochastic gravitational wave background generated by cosmic strings. We argue that, if observed by current or forthcoming experiments, this signature would provide strong evidence for i-strings and, therefore, for (primordial) inflation. We also develop a simple single-parameter algorithm for the computation of the stochastic gravitational wave background generated by i-strings from that of a standard cosmic string network, which may be useful in the determination of the observational constraints to be obtained by current and forthcoming gravitational wave experiments.

The generalized Chaplygin gas is usually defined as a barotropic perfect fluid with an equation of state $p=-A \rho^{-\alpha}$, where $\rho$ and $p$ are the proper energy density and pressure, respectively, and $A$ and $\alpha$ are positive real parameters. It has been extensively studied in the literature as a quartessence prototype unifying dark matter and dark energy. Here, we consider an extended family of generalized Chaplygin gas models parameterized by three positive real parameters $A$, $\alpha$ and $\beta$, which, for two specific choices of $\beta$ [$\beta=1$ and $\beta=\left(1+\alpha\right)/(2\alpha)$], is described by two different Lagrangians previously identified in the literature with the generalized Chaplygin gas. We show that, for $\beta > 1/2$, the linear stability conditions and the maximum value of the sound speed $c_s$ are regulated solely by $\beta$, with $0 \le c_s \le 1$ if $\beta \ge 1$. We further demonstrate that in the non-relativistic regime the standard equation of state $p=-A \rho^{-\alpha}$ of the generalized Chaplygin gas is always recovered, while in the relativistic regime this is true only if $\beta=\left(1+\alpha\right)/(2\alpha)$. We present a regularization of the ($\alpha\rightarrow 0$, $A \rightarrow \infty$) limit of the generalized Chaplygin gas, showing that it leads to a logarithmic Chaplygin gas model with an equation of state of the form $p = {\mathcal A} \ln\left(\rho/\rho_{*}\right)$, where ${\mathcal A}$ is a real parameter and $\rho_*>0$ is an arbitrary energy density. We finally derive its Lagrangian formulation.

In this paper we show that the baryon-to-photon ratio $\eta$ is in general not conserved before decoupling in non-minimally coupled theories of gravity. We use big-bang nucleosynthesis and cosmic microwave background limits on the baryon-to-photon ratio $\eta$ to derive new constraints on modified gravity theories with a universal non-minimal coupling between matter and curvature, showing that they rule out a specific class of models previously considered in the literature as a substitute for the dark matter. We also compare these new constraints with the ones obtained from the COBE-FIRAS limits on CMB spectral distortions, highlighting the complementarity between them.

We study the proposal to solve the coincidence problem in the non-local version of the vacuum energy sequestering mechanism by means of a scalar field in a linear potential. We show that there is no solution in the theory compatible with observations if one requires the scalar field to drive the present period of acceleration and the collapse.

In this paper we show that the on-shell Lagrangian of a perfect fluid depends on microscopic properties of the fluid, giving specific examples of perfect fluids with different on-shell Lagrangians but with the same energy-momentum tensor. We demonstrate that if the fluid is constituted by localized concentrations of energy with fixed rest mass and structure (solitons) then the average on-shell Lagrangian of a perfect fluid is given by ${\mathcal L}_m=T$, where $T$ is the trace of the energy-momentum tensor. We show that our results have profound implications for theories of gravity where the matter Lagrangian appears explicitly in the equations of motion of the gravitational and matter fields, potentially leading to observable deviations from a nearly perfect cosmic microwave background black body spectrum: $n$-type spectral distortions, affecting the normalization of the spectral energy density. Finally, we put stringent constraints on $f(R,{\mathcal L}_m)$ theories of gravity using the COBE-FIRAS measurement of the spectral radiance of the cosmic microwave background.

We consider a model where particles are described as localized concentrations of energy, with fixed rest mass and structure, which are not significantly affected by their self-induced gravitational field. We show that the volume average of the on-shell matter Lagrangian ${\mathcal L_m}$ describing such particles, in the proper frame, is equal to the volume average of the trace $T$ of the energy-momentum tensor in the same frame, independently of the particle's structure and constitution. Since both ${\mathcal L_m}$ and $T$ are scalars, and thus independent of the reference frame, this result is also applicable to collections of moving particles and, in particular, to those which can be described by a perfect fluid. Our results are expected to be particularly relevant in the case of modified theories of gravity with nonminimal coupling to matter where the matter Lagrangian appears explicitly in the equations of motion of the gravitational and matter fields, such as $f(R,{\mathcal L_m})$ and $f(R,T)$ gravity. In particular, they indicate that, in this context, $f(R,{\mathcal L_m})$ theories may be regarded as a subclass of $f(R,T)$ gravity.

Observations of Sun-like stars over the last half-century have improved our understanding of how magnetic dynamos, like that responsible for the 11-year solar cycle, change with rotation, mass and age. Here we show for the first time how metallicity can affect a stellar dynamo. Using the most complete set of observations of a stellar cycle ever obtained for a Sun-like star, we show how the solar analog HD 173701 exhibits solar-like differential rotation and a 7.4-year activity cycle. While the duration of the cycle is comparable to that generated by the solar dynamo, the amplitude of the brightness variability is substantially stronger. The only significant difference between HD 173701 and the Sun is its metallicity, which is twice the solar value. Therefore, this provides a unique opportunity to study the effect of the higher metallicity on the dynamo acting in this star and to obtain a comprehensive understanding of the physical mechanisms responsible for the observed photometric variability. The observations can be explained by the higher metallicity of the star, which is predicted to foster a deeper outer convection zone and a higher facular contrast, resulting in stronger variability.

We revisit the physical properties of global and local monopoles and discuss their implications in the dynamics of monopole networks. In particular, we review the Velocity-dependent One-Scale (VOS) model for global and local monopoles and propose physically motivated changes to its equations. We suggest a new form for the acceleration term of the evolution equation of the root-mean-squared velocity and show that, with this change, the VOS model is able to describe the results of radiation and matter era numerical simulations of global monopole networks with a single value of the acceleration parameter $k$, thus resolving the tension previously found in the literature. We also show that the fact that the energy of global monopoles is not localized within their cores affects their dynamics and, thus, the Hubble damping terms in the VOS equations. We study the ultra-relativistic linear scaling regime predicted by the VOS equations and demonstrate that it cannot be attained either on radiation or matter eras and, thus, cannot arise from the cosmological evolution of a global monopole network. We also briefly discuss the implications of our findings for the VOS model for local monopoles.

Deciphering the assembly history of the Milky Way is a formidable task, which becomes possible only if one can produce high-resolution chrono-chemo-kinematical maps of the Galaxy. Data from large-scale astrometric and spectroscopic surveys will soon provide us with a well-defined view of the current chemo-kinematical structure of the Milky Way, but will only enable a blurred view on the temporal sequence that led to the present-day Galaxy. As demonstrated by the (ongoing) exploitation of data from the pioneering photometric missions CoRoT, Kepler, and K2, asteroseismology provides the way forward: solar-like oscillating giants are excellent evolutionary clocks thanks to the availability of seismic constraints on their mass and to the tight age-initial-mass relation they adhere to. In this paper we identify five key outstanding questions relating to the formation and evolution of the Milky Way that will need precise and accurate ages for large samples of stars to be addressed, and we identify the requirements in terms of number of targets and the precision on the stellar properties that are needed to tackle such questions. By quantifying the asteroseismic yields expected from PLATO for red-giant stars, we demonstrate that these requirements are within the capabilities of the current instrument design, provided that observations are sufficiently long to identify the evolutionary state and allow robust and precise determination of acoustic-mode frequencies. This will allow us to harvest data of sufficient quality to reach a 10% precision in age. This is a fundamental pre-requisite to then reach the more ambitious goal of a similar level of accuracy, which will only be possible if we have to hand a careful appraisal of systematic uncertainties on age deriving from our limited understanding of stellar physics, a goal which conveniently falls within the main aims of PLATO's core science.

We consider a model with two real scalar fields which admits phantom domain wall solutions. We investigate the structure and evolution of these phantom domain walls in an expanding homogeneous and isotropic universe. In particular, we show that the increase of the tension of the domain walls with cosmic time, associated to the evolution of the phantom scalar field, is responsible for an additional damping term in their equations of motion. We describe the macroscopic dynamics of phantom domain walls, showing that extended phantom defects whose tension varies on a cosmological timescale cannot be the dark energy.

A unification of dark matter and dark energy in terms of a logotropic perfect dark fluid has recently been proposed, where deviations with respect to the standard $\Lambda {\rm CDM}$ model are dependent on a single parameter $B$. In this paper we show that the requirement that the linear growth of cosmic structures on comoving scales larger than $8 h^{-1} \, {\rm Mpc}$ is not significantly affected with respect to the standard $\Lambda {\rm CDM}$ result provides the strongest constraint to date on the model ($B <6 \times 10^{-7}$), an improvement of more than three orders of magnitude over previous constraints on the value of $B$. We further show that this constraint rules out the logotropic Unified Dark Energy model as a possible solution to the small scale problems of the $\Lambda$CDM model, including the cusp problem of Dark Matter halos or the missing satellite problem, as well as the original version of the model where the Planck energy density was taken as one of the two parameters characterizing the logotropic dark fluid.

Quasi-periodic modulations of the stellar light curve may result from dark spots crossing the visible stellar disc. Due to differential rotation, spots at different latitudes generally have different rotation periods. Hence, by studying spot-induced modulations, one can learn about stellar surface (differential) rotation and magnetic activity. Recently, Reinhold & Arlt (2015) proposed a method based on the Lomb-Scargle periodogram of light curves to identify the sign of the differential rotation at the stellar surface. Our goal is to understand how the modulation of the stellar light curve due to the presence of spots and the corresponding periodogram are affected by both the stellar and spot properties. We generate synthetic light curves of stars with different properties (inclination angle, limb darkening, and rotation rate) and spot configurations (number of spots, latitude, intensity contrast, and size). By analysing their Lomb-Scargle periodograms, we compute the ratio between the heights of the second and first harmonics of the rotation period (peak-height-ratio). We find that the peak-height ratios are essentially a function of a single parameter, the fraction of time the spot is visible, which is related to the sinusoidality of the spot modulation. We identify the conditions under which the periodogram analysis can actually provide an estimate of the spot latitudes and/or the stellar inclination angle. We also identify possible sources of error in the identification of the sign of the differential rotation.

The activity-related variations in the solar acoustic frequencies have been known for 30 years. However, the importance of the different contributions is still not well established. With this in mind, we developed an empirical model to estimate the spot-induced frequency shifts, which takes into account the sunspot properties, such as area and latitude. The comparison between the model frequency shifts obtained from the daily sunspot records and those observed suggests that the contribution from a stochastic component to the total frequency shifts is about 30%. The remaining 70% is related to a global, long-term variation. We also propose a new observable to investigate the short- and mid-term variations of the frequency shifts, which is insensitive to the long-term variations contained in the data. On the shortest time scales the variations in the frequency shifts are strongly correlated with the variations in the total area covered by sunspots. However, a significant loss of correlation is still found, which cannot be fully explained by ignoring the invisible side of the Sun when accounting for the total sunspot area. We also verify that the times when the frequency shifts and the sunspot areas do not vary in a similar way tend to coincide with the times of the maximum amplitude of the quasi-biennial variations found in the seismic data.

Starspots are cooler and darker than the stellar surface. Therefore, the emitted flux of a star changes when spots are visible on its surface. The presence of spots together with the stellar rotation leads to a periodic modulation on the light curve. By studying that modulation one can then learn about the stellar rotation and also magnetic activity. Recently, Reinhold & Arlt (2015) proposed a method based on the analysis of the Lomb Scargle periodogram of the light curve to identify the sign of the differential rotation, i.e. whether the equator rotates faster than the poles or the opposite. We have been studying in detail the spots' impact on the light curve and on the resulting periodogram. We find that, under some conditions, the periodogram can actually provide an estimate of the true spot latitudes and/or the stellar inclination angle. Moreover, we find that the impact of the spot on the ratio between the heights of the second and first harmonics of the main peaks in the periodogram can be described by a single parameter, the visibility time of the spot. Finally, we also identify possible sources of false positives/negatives for the sign of the differential rotation.

The frequencies of the solar acoustic oscillations vary over the activity cycle. The variations in other activity proxies are found to be well correlated with the variations in the acoustic frequencies. However, each proxy has a slightly different time behaviour. Our goal is to characterize the differences between the time behaviour of the frequency shifts and of two other activity proxies, namely, the area covered by sunspots and the 10.7cm flux. We define a new observable that is particularly sensitive to the short-term frequency variations. We then compare the observable when computed from model frequency shifts and from observed frequency shifts obtained with the Global Oscillation Network Group (GONG) for cycle 23. Our analysis shows that on the shortest time-scales the variations in the frequency shifts seen in the GONG observations are strongly correlated with the variations in the area covered by sunspots. However, a significant loss of correlation is still found. We verify that the times when the frequency shifts and the sunspot area do not vary in a similar way tend to coincide with the times of the maxima of the quasi-biennial variations seen in the solar seismic data. A similar analysis of the relation between the 10.7cm flux and the frequency shifts reveals that the short-time variations in the frequency shifts follow even more closely those of the 10.7cm flux than those of the sunspot area. However, a loss of correlation between frequency shifts and 10.7cm flux variations is still found around the same times.

We explore the dynamics and evolution of the Universe at early and late times, focusing on both dark energy and extended gravity models and their astrophysical and cosmological consequences. Modified theories of gravity not only provide an alternative explanation for the recent expansion history of the universe, but they also offer a paradigm fundamentally distinct from the simplest dark energy models of cosmic acceleration. In this review, we perform a detailed theoretical and phenomenological analysis of different modified gravity models and investigate their consistency. We also consider the cosmological implications of well motivated physical models of the early universe with a particular emphasis on inflation and topological defects. Astrophysical and cosmological tests over a wide range of scales, from the solar system to the observable horizon, severely restrict the allowed models of the Universe. Here, we review several observational probes -- including gravitational lensing, galaxy clusters, cosmic microwave background temperature and polarization, supernova and baryon acoustic oscillations measurements -- and their relevance in constraining our cosmological description of the Universe.

We compute the stochastic gravitational wave background generated by cosmic superstrings using a semi-analytical velocity-dependent model to describe their dynamics. We show that heavier string types may leave distinctive signatures on the stochastic gravitational wave background spectrum within the reach of present and upcoming gravitational wave detectors. We examine the physically motivated scenario in which the physical size of loops is determined by the gravitational backreaction scale and use NANOGRAV data to derive a conservative constraint of $G\mu_F<3.2 \times 10^{-9}$ on the tension of fundamental strings. We demonstrate that approximating the gravitational wave spectrum generated by cosmic superstring networks using the spectrum generated by ordinary cosmic strings with reduced intercommuting probability (which is often done in the literature) leads, in general, to weaker observational constraints on $G\mu_F$. We show that the inclusion of heavier string types is required for a more accurate characterization of the region of the $(g_s,G\mu_F)$ parameter space that may be probed using direct gravitational wave detectors. In particular, we consider the observational constraints that result from NANOGRAV data and show that heavier strings generate a secondary exclusion region of parameter space.

Activity-related variations in the solar oscillation properties have been known for 30 years. However, the relative importance of the different contributions to the observed variations is not yet fully understood. Our goal is to estimate the relative contribution from sunspots to the observed activity-related variations in the frequencies of the acoustic modes. We use a variational principle to relate the phase differences induced by sunspots on the acoustic waves to the corresponding changes in the frequencies of the global acoustic oscillations. From the sunspot properties (area and latitude as a function of time), we are able to estimate the spot-induced frequency shifts. These are then combined with a smooth frequency shift component, associated with long-term solar-cycle variations, and the results compared with the frequency shifts derived from the Global Oscillation Network Group (GONG) data. The result of this comparison is consistent with a sunspot contribution to the observed frequency shifts of roughly 30 per cent, with the remaining 70 per cent resulting mostly from a global, non-stochastic variation, possibly related to the changes in the overall magnetic field. Moreover, analysis of the residuals obtained after the subtraction of the model frequency shifts from the observations indicates the presence of a 1.5-yr periodicity in the data in phase with the quasi-biennial variations reported in the literature.

In the letter by Stadnik and Flambaum [Phys. Rev. Lett. 113, 151301 (2014)] it is claimed that topological defects passing through pulsars could be responsible for the observed pulsar glitches. Here, we show that, independently of the detailed network dynamics and defect dimensionality, such proposal is faced with serious difficulties.

We study the inner dynamics of accreting Eddington-inspired Born-Infeld black holes using the homogeneous approximation and taking charge as a surrogate for angular momentum. We show that there is a minimum of the accretion rate below which mass inflation does not occur, and we derive an analytical expression for this threshold as a function of the fundamental scale of the theory, the accretion rate, the mass, and the charge of the black hole. Our result explicitly demonstrates that, no matter how close Eddington-inspired Born-Infeld gravity is to general relativity, there is always a minimum accretion rate below which there is no mass inflation. For larger accretion rates, mass inflation takes place inside the black hole as in general relativity until the extremely rapid density variations bring it to an abrupt end. We derive analytical scaling solutions for the value of the energy density and of the Misner-Sharp mass attained at the end of mass inflation as a function of fundamental scale of the theory, the accretion rate, the mass, and the charge of the black hole, and compare these with the corresponding numerical solutions. We find that, except for unreasonably high accretion rates, our analytical results appear to provide an accurate description of homogeneous mass inflation inside accreting Eddington-inspired Born-Infeld black holes.

We critically assess recent claims suggesting that upper limits on the time variation of the fine-structure constant tightly constrain the coupling of a dark energy scalar field to the electromagnetic sector, and, indirectly, the violation of the weak equivalence principle. We show that such constraints depend crucially on the assumed priors, even if the dark energy was described by a dynamical scalar field with a constant equation of state parameter $w$ linearly coupled to the electromagnetic sector through a dimensionless coupling $\zeta$. We find that, although local atomic clock tests, as well as other terrestrial, astrophysical and cosmological data, put stringent bounds on $|\zeta| {\sqrt {|w+1|}}$, the time variation of the fine-structure constant cannot be used to set or to improve upper limits on $|\zeta|$ or $|w+1|$ without specifying priors, consistent but not favoured by current data, which strongly disfavour low values of $|w+1|$ or $|\zeta|$, respectively. We briefly discuss how this might change with a new generation of high-resolution ultra-stable spectrographs, such as ESPRESSO and ELT-HIRES, in combination with forthcoming missions to map the geometry of the Universe, such as Euclid, or to test the equivalence principle, such as MICROSCOPE or STEP.

We investigate the interior dynamics of accreting black holes in Eddington-inspired Born-Infeld gravity using the homogeneous approximation and taking charge as a surrogate for angular momentum, showing that accretion can have an enormous impact on their inner structure. We find that, unlike in general relativity, there is a minimum accretion rate bellow which the mass inflation instability, which drives the centre-of-mass streaming density to exponentially high values in an extremely short interval of time, does not occur. We further show that, above this threshold, mass inflation takes place inside black holes very much in the same way as in general relativity, but is brought to a halt at a maximum energy density which is, in general, much smaller than the fundamental energy density of the theory. We conjecture that some of these results may be a common feature of modified gravity theories in which significant deviations from general relativity manifest themselves at very high densities.

We show that simple thermodynamic conditions determine, to a great extent, the equation of state and dynamics of cosmic defects of arbitrary dimensionality. We use these conditions to provide a more direct derivation of the Velocity-dependent One-Scale (VOS) model for the macroscopic dynamics of topological defects of arbitrary dimensionality in a $N+1$-dimensional homogeneous and isotropic universe. We parameterize the modifications to the VOS model associated to the interaction of the topological defects with other fields, including, in particular, a new dynamical degree of freedom associated to the variation of the mass per unit $p$-area of the defects, and compute the corresponding scaling solutions. The observational impact of this new dynamical degree of freedom is also briefly discussed.

The increasing interest in understanding stellar magnetic activity cycles is a strong motivation for the development of parameterized starspot models which can be constrained observationally. In this work we develop an empirical tool for the stochastic reconstruction of sunspot cycles, using the average solar properties as a reference. The synthetic sunspot cycle is compared with the sunspot data extracted from the National Geophysical Data Center, in particular using the Kolmogorov-Smirnov test. This tool yields synthetic spot group records, including date, area, latitude, longitude, rotation rate of the solar surface at the group's latitude, and an identification number. Comparison of the stochastic reconstructions with the daily sunspot records confirms that our empirical model is able to successfully reproduce the main properties of the solar sunspot cycle. As a by-product of this work, we show that the Gnevyshev-Waldmeier rule, which describes the spots' area-lifetime relation, is not adequate for small groups and we propose an effective correction to that relation which leads to a closer agreement between the synthetic sunspot cycle and the observations.

We develop a numerical tool for the fast computation of the temperature and polarization power spectra generated by domain wall networks, by extending the publicly available CMBACT code --- that calculates the CMB signatures generated by active sources --- to also describe domain wall networks. In order to achieve this, we adapt the Unconnected Segment model for cosmic strings to also describe domain wall networks, and use it to model the energy-momentum of domain wall networks throughout their cosmological history. We use this new tool to compute and study the TT, EE, TE and BB power spectra generated by standard domain wall networks, and derive a conservative constraint on the energy scale of the domain wall-forming phase transition of $\upeta <0.92\,\,{\rm MeV}$ (which is a slight improvement over the original Zel'dovich bound of $1\,\,{\rm MeV}$).

We consider the possibility that current hints of spatial variations of the fine structure constant at high redshift, based on VLT/UVES and Keck/HIRES observations, could be caused by a biased domain wall network described by a scalar field non-minimally coupled to the electromagnetic field. We show that in order to be responsible for the reported spatial variations of the fine structure constant, the fractional contribution of the domain wall network to the energy density of the Universe should be tightly constrained within the range $10^{-10} < \Omega_{w0} < 10^{-5}$. We also show that the domain wall dynamics should be essentially frictionless, so that its characteristic scale is in the order of the Hubble radius at the present time.

With recent advances in asteroseismology it is now possible to peer into the cores of red giants, potentially providing a way to study processes such as nuclear burning and mixing through their imprint as sharp structural variations -- glitches -- in the stellar cores. Here we show how such core glitches can affect the oscillations we observe in red giants. We derive an analytical expression describing the expected frequency pattern in the presence of a glitch. This formulation also accounts for the coupling between acoustic and gravity waves. From an extensive set of canonical stellar models we find glitch-induced variation in the period spacing and inertia of non-radial modes during several phases of red-giant evolution. Significant changes are seen in the appearance of mode amplitude and frequency patterns in asteroseismic diagrams such as the power spectrum and the \'echelle diagram. Interestingly, along the red-giant branch glitch-induced variation occurs only at the luminosity bump, potentially providing a direct seismic indicator of stars in that particular evolution stage. Similarly, we find the variation at only certain post-helium-ignition evolution stages, namely, in the early phases of helium-core burning and at the beginning of helium-shell burning signifying the asymptotic-giant-branch bump. Based on our results, we note that assuming stars to be glitch-free, while they are not, can result in an incorrect estimate of the period spacing. We further note that including diffusion and mixing beyond classical Schwarzschild, could affect the characteristics of the glitches, potentially providing a way to study these physical processes.

We investigate the potential impact of mass inflation inside black holes on the dynamics of the Universe, considering a recent reformulation of general relativity, proposed in [1], which prevents the vacuum energy from acting as a gravitational source. The interior dynamics of accreting black holes is studied, at the classical level, using the homogeneous approximation and taking charge as a surrogate for angular momentum. We show that, depending on the accreting fluid properties, mass inflation inside black holes could influence the value of the cosmological constant and thus the dynamics of the Universe. A full assessment of the cosmological role played by black holes will require a deeper understanding of the extremely energetic regimes expected inside real astrophysical black holes, including their relation with the physics of the very early Universe, and may eventually lead to an entirely new paradigm for the origin and evolution of the Universe.

We investigate the impact on the classical dynamics of dark matter particles and dark energy of a non-minimal coupling in the dark sector, assuming that the mass of the dark matter particles is coupled to a dark energy scalar field. We show that standard results can only be recovered if the space-time variation of the dark energy scalar field is sufficiently smooth on the characteristic length scale of the dark matter particles, and we determine the associated constraint dependent on both the mass and radius of the dark matter particles and the coupling to the dark energy scalar field. We further show, using field theory numerical simulations, that a violation of such constraint results in a microscopic feedback effect strongly affecting the dynamics of dark matter particles, with a potential impact on structure formation and on the space-time evolution of the dark energy equation of state.

We show that the observed rotation curves of spiral galaxies constrain the sound speed of the dark matter to be $c_s < 10^{-4} c$, where $c$ is the speed of light in vacuum. Using the Modified Chaplygin Gas as a representative example of a class of unified dark energy models incorporating an effective dark matter component with a non-zero sound speed, we determine the most stringent constraint to date on the value of the constant contribution to the equation of state parameter in this class of models. Finally, we explain the reason why previous constraints using the Cosmic Microwave Background and Baryonic Acoustic Oscillations were not as competitive as the one presented in this paper and discuss the limitations of the recently proposed Extended Chaplygin Gas.

We explicitly compute the dynamics of closed homogeneous and isotropic universes permeated by a single perfect fluid with a constant equation of state parameter $w$ in the context of a recent reformulation of general relativity, proposed in [1], which prevents the vacuum energy from acting as a gravitational source. This is done using an iterative algorithm, taking as an initial guess the background cosmological evolution obtained using standard general relativity in the absence of a cosmological constant. We show that, in general, the impact of the vacuum energy sequestering mechanism on the dynamics of the universe is significant, except for the $w=1/3$ case where the results are identical to those obtained in the context of general relativity with a null cosmological constant. We also show that there are well behaved models in general relativity that do not have a well behaved counterpart in the vacuum energy sequestering paradigm studied in this paper, highlighting the specific case of a quintessence scalar field with a linear potential.

We study the nonlinear regime of Unified Dark Energy models, using Generalized Chaplygin Gas cosmologies as a representative example, and introduce a new parameter characterizing the level of small scale clustering in these scenarios. We show that viable Generalized Chaplygin Gas cosmologies, consistent with the most recent observational constraints, may be constructed for any value of the Generalized Chaplygin Gas parameter by considering models with a sufficiently high level of nonlinear clustering.

We consider an alternative approach for the computation of the stochastic gravitational wave background generated by small loops produced throughout the cosmological evolution of cosmic string networks and use it to derive an analytical approximation to the corresponding power spectrum. We show that this approximation produces an excellent fit to more elaborate results obtained using the Velocity-dependent One-Scale model to describe cosmic string network dynamics, over a wide frequency range, in the small-loop regime.

In this letter we give specific examples of Z_N Lotka-Volterra competition models leading to the formation of string networks. We show that, in order to promote coexistence, the species may arrange themselves around regions with a high number density of empty sites generated by predator-prey interactions between competing species. These configurations extend into the third dimension giving rise to string networks. We investigate the corresponding dynamics using both stochastic and mean field theory simulations, showing that the coarsening of these string networks follows a scaling law which is analogous to that found in other physical systems in condensed matter and cosmology.

We derive, using the spherical collapse model, a generalized Layzer-Irvine equation which can be used to describe the gravitational collapse of cold dark matter in a dark energy background. We show that the usual Layzer-Irvine equation is valid if the dark matter and the dark energy are minimally coupled to each other and the dark energy distribution is homogeneous, independently of its equation of state. We compute the corrections to the standard Layzer-Irvine equation which arise in the presence of dark energy inhomogeneities. We show that, in the case of a dark energy component with a constant equation of state parameter consistent with the latest observational constraints, these corrections are expected to be small, even if the dark energy has a negligible sound speed. However, we find that, in more general models, the impact of dark energy perturbations on the dynamics of clusters of galaxies, which will be constrained by ESA's Euclid mission with unprecedented precision, might be significant.

The redshift dependence of the abundance of galaxy clusters is very sensitive to the statistical properties of primordial density perturbations. It can thus be used to probe small deviations from Gaussian initial conditions. Such deviations constitute a very important signature of many inflationary scenarios, and are thus expected to provide crucial information on physical processes which took place in the very early Universe. We have determined the biases which may be introduced in the estimation of cosmological parameters by wrongly assuming the absence of primordial non-Gaussianities. Although we find that the estimation of the present-day dark energy density using cluster counts is not very sensitive to the non-Gaussian properties of the density field, we show that the biases can be considerably larger in the estimation of the dark energy equation of state parameter $w$ and of the amplitude of the primordial density perturbations. Our results suggest that a significant level of non-Gaussianity at cluster scales may be able to reconcile the constraint on the amplitude of the primordial perturbations obtained using galaxy cluster number counts from the Planck Sunyaev-Zeldovich Catalog with that obtained from the primary Cosmic Microwave Background anisotropies measured by the Planck satellite.

We compute the power spectrum of the stochastic gravitational wave background generated by cosmic string networks described by the Velocity-Dependent One-Scale (VOS) model, for a wide range of macroscopic and microscopic parameters. The VOS model --- which has been shown to provide an accurate macroscopic description of the evolution of cosmic string networks --- is used to demonstrate that cosmic string networks are unable to rapidly attain scale-invariant evolution after the transition between the radiation and matter eras. However, in computations of the stochastic gravitational wave background, it is often assumed that the networks experience scale-invariant evolution throughout cosmological history. We demonstrate that this assumption leads to an underestimation of the amplitude and broadness of the peak of the spectrum, that may consequently lead to inaccurate observational constraints on the cosmic string tension.

We assess the viability of successful reconstruction of the evolution of the dark energy equation of state using varying fundamental couplings, such as the fine structure constant or the proton-to-electron mass ratio. We show that the same evolution of the dark energy equation of state parameter with cosmic time may be associated with arbitrary variations of the fundamental couplings. Various examples of models with the same (different) background evolution and different (the same) time variation of fundamental couplings are studied in the letter. Although we demonstrate that, for a broad family of models, it is possible to redefine the scalar field in such a way that its dynamics is that of a standard quintessence scalar field, in general such redefinition leads to the breakdown of the linear relation between the scalar field and the variation of fundamental couplings. This implies that the assumption of a linear coupling is not sufficient to guarantee a successful reconstruction of the dark energy dynamics and consequently additional model dependent assumptions about the scalar field responsible for the dark energy need to be made.

In this paper we investigate the classical non-relativistic limit of the Eddington-inspired Born-Infeld theory of gravity. We show that strong bounds on the value of the only additional parameter of the theory \kappa, with respect to general relativity, may be obtained by requiring that gravity plays a subdominant role compared to electromagnetic interactions inside atomic nuclei. We also discuss the validity of the continuous fluid approximation used in this and other astrophysical and cosmological studies. We argue that although the continuous fluid approximation is expected to be valid in the case of sufficiently smooth density distributions, its use should eventually be validated at a quantum level.

We study the dynamics of a homogeneous and isotropic Friedmann-Robertson-Walker universe in the context of the Eddington-inspired Born-Infeld theory of gravity. We generalize earlier results, obtained in the context of a radiation dominated universe, to account for the evolution of a universe permeated by a perfect fluid with an arbitrary equation of state parameter $w$. We show that a bounce may occur for $\kappa >0$, if $w$ is time-dependent, and we demonstrate that it is free from tensor singularities. We argue that Eddington-inspired Born-Infeld cosmologies may be a viable alternative to the inflationary paradigm as a solution to fundamental problems of the standard cosmological model.

The evolution with time of the abundance of galaxy clusters is very sensitive to the statistical properties of the primordial density perturbations. It can thus be used to probe small deviations from Gaussianity in the initial conditions. The characterization of such deviations would help distinguish between different inflationary scenarios, and provide us with information on physical processes which took place in the early Universe. We have found that when the information contained in the galaxy cluster counts is used to reconstruct the dark energy equation of state as a function of redshift, assuming erroneously that no primordial non-Gaussianities exist, an apparent evolution with time in the effective dark energy equation of state arises,characterized by the appearance of a clear discontinuity.

In this paper we find new scaling laws for the evolution of $p$-brane networks in $N+1$-dimensional Friedmann-Robertson-Walker universes in the weakly-interacting limit, giving particular emphasis to the case of cosmic superstrings ($p=1$) living in a universe with three spatial dimensions (N=3). In particular, we show that, during the radiation era, the root-mean-square velocity is ${\bar v} =1/{\sqrt 2}$ and the characteristic length of non-interacting cosmic string networks scales as $L \propto a^{3/2}$ ($a$ is the scale factor), thus leading to string domination even when gravitational backreaction is taken into account. We demonstrate, however, that a small non-vanishing constant loop chopping efficiency parameter $\tilde c$ leads to a linear scaling solution with constant $L H \ll 1$ ($H$ is the Hubble parameter) and ${\bar v} \sim 1/{\sqrt 2}$ in the radiation era, which may allow for a cosmologically relevant cosmic string role even in the case of light strings. We also determine the impact that the radiation-matter transition has on the dynamics of weakly interacting cosmic superstring networks.

In this letter we compute stringent astrophysical and cosmological constraints on a recently proposed Eddington-inspired Born-Infeld theory of gravity. We find, using a generalized version of the Zel'dovich approximation, that in this theory a pressureless cold dark matter fluid has a non-zero effective sound speed. We compute the corresponding effective Jeans length and show that it is approximately equal to the fundamental length of the theory $R_*=\kappa^{1/2} G^{-1/2}$, where $\kappa$ is the only additional parameter of theory with respect to general relativity and $G$ is the gravitational constant. This scale determines the minimum size of compact objects which are held together by gravity. We also estimate the critical mass above which pressureless compact objects are unstable to colapse into a black hole, showing that it is approximately equal to the fundamental mass $M_* = \kappa^{1/2} c^2 G^{-3/2}$, and we show that the maximum density attainable inside stable compact stars is roughly equal to the fundamental density $\rho_*=\kappa^{-1} c^2$, where $c$ is the speed of light in vacuum. We find that the mere existence of astrophysical objects of size $R$ which are held together by their own gravity leads to the constraint $\kappa < G R^2$. In the case of neutron stars this implies that $\kappa < 10^{-2} \, {\rm m^5 \, kg^{-1} \, s^{-2}}$, a limit which is stronger by about 10 orders of magnitude than big bang nucleosynthesis constraints and by more than 7 orders of magnitude than solar constraints.

In models where dark matter and dark energy interact non-minimally, the total amount of matter in a fixed comoving volume may vary from the time of recombination to the present time due to energy transfer between the two components. This implies that, in interacting dark energy models, the fractional matter density estimated using the cosmic microwave background assuming no interaction between dark matter and dark energy will in general be shifted with respect to its true value. This may result in an incorrect determination of the equation of state of dark energy if the interaction between dark matter and dark energy is not properly accounted for, even if the evolution of the Hubble parameter as a function of redshift is known with arbitrary precision. In this paper we find an exact expression, as well as a simple analytical approximation, for the evolution of the effective equation of state of dark energy, assuming that the energy transfer rate between dark matter and dark energy is described by a simple two-parameter model. We also provide analytical examples where non-phantom interacting dark energy models mimic the background evolution and primary cosmic microwave background anisotropies of phantom dark energy models.

We derive the Layzer-Irvine equation in the presence of a homogeneous (or quasi-homogeneous) dark energy component with an arbitrary equation of state. We extend the Layzer-Irvine equation to homogeneous and isotropic universes with an arbitrary number of dimensions and obtain the corresponding virial relation for sufficiently relaxed objects. We find analogous equations describing the dynamics of cosmic string loops and other p-branes of arbitrary dimensionality, discussing the corresponding relativistic and non-relativistic limits. Finally, we generalize the Layzer-Irvine equation to account for a non-minimal interaction between dark matter and dark energy, discussing its practical use as a signature of such an interaction.

In this paper we derive, directly from the Nambu-Goto action, the relevant components of the acceleration of cosmological featureless $p$-branes, extending previous analysis based on the field theory equations in the thin-brane limit. The component of the acceleration parallel to the velocity is at the core of the velocity-dependent one-scale model for the evolution of $p$-brane networks. We use this model to show that, in a decelerating expanding universe in which the $p$-branes are relevant cosmologically, interactions cannot lead to frustration, except for fine-tuned non-relativistic networks with a dimensionless curvature parameter $k \ll 1$. We discuss the implications of our findings for the cosmological evolution of $p$-brane networks.

In this paper we compare the performances of the chi-square and median likelihood analysis in the determination of cosmological constraints using type Ia supernovae data. We perform a statistical analysis using the 307 supernovae of the Union 2 compilation of the Supernova Cosmology Project and find that the chi-square statistical analysis yields tighter cosmological constraints than the median statistic if only supernovae data is taken into account. We also show that when additional measurements from the Cosmic Microwave Background and Baryonic Acoustic Oscillations are considered, the combined cosmological constraints are not strongly dependent on whether one applies the chi-square statistic or the median statistic to the supernovae data. This indicates that, when complementary information from other cosmological probes is taken into account, the performances of the chi-square and median statistics are very similar, demonstrating the robustness of the statistical analysis.

We explicitly demonstrate the existence of static global defect solutions of arbitrary dimensionality whose energy does not diverge at spatial infinity, by considering maximally symmetric solutions described by an action with non-standard kinetic terms in a D+1 dimensional Minkowski space-time. We analytically determine the defect profile both at small and large distances from the defect centre. We verify the stability of such solutions and discuss possible implications of our findings, in particular for dark matter and charge fractionalization in graphene.

In this work we determine the correspondence between quintessence and tachyon dark energy models with a constant dark energy equation of state parameter, $w_e$. Although the evolution of both the Hubble parameter and the scalar field potential with redshift is the same, we show that the evolution of quintessence/tachyon scalar fields with redshift is, in general, very different. We explicity demonstrate that if $w_e \neq -1$ the potentials need to be very fine-tuned for the relative perturbation on the equation of state parameter, $\Delta w_e/(1+w_e) \ll 1$, to be very small around the present time. We also discuss possible implications of our results for the reconstruction of the evolution of $w_e$ with redshift using varying couplings.

In this paper we investigate possible solutions to the coincidence problem in flat phantom dark energy models with a constant dark energy equation of state and quintessence models with a linear scalar field potential. These models are representative of a broader class of cosmological scenarios in which the universe has a finite lifetime. We show that, in the absence of anthropic constraints, including a prior probability for the models inversely proportional to the total lifetime of the universe excludes models very close to the $\Lambda {\rm CDM}$ model. This relates a cosmological solution to the coincidence problem with a dynamical dark energy component having an equation of state parameter not too close to -1 at the present time. We further show, that anthropic constraints, if they are sufficiently stringent, may solve the coincidence problem without the need for dynamical dark energy.

We develop a velocity-dependent one-scale model describing p-brane dynamics in flat homogeneous and isotropic backgrounds in a unified framework. We find the corresponding scaling laws in frictionless and friction dominated regimes considering both expanding and collapsing phases.

The Press-Ryden-Spergel (PRS) algorithm is a modification to the field theory equations of motion, parametrized by two parameters ($\alpha$ and $\beta$), implemented in numerical simulations of cosmological domain wall networks, in order to ensure a fixed comoving resolution. In this paper we explicitly demonstrate that the PRS algorithm provides the correct domain wall dynamics in $N+1$-dimensional Friedmann-Robertson-Walker (FRW) universes if $\alpha+\beta/2=N$, fully validating its use in numerical studies of cosmic domain evolution. We further show that this result is valid for generic thin featureless domain walls, independently of the Lagrangian of the model.

We develop a velocity-dependent one-scale model for the evolution of domain wall networks in flat expanding or collapsing homogeneous and isotropic universes with an arbitrary number of spatial dimensions, finding the corresponding scaling laws in frictionless and friction dominated regimes. We also determine the allowed range of values of the curvature parameter and the expansion exponent for which a linear scaling solution is possible in the frictionless regime.

In this paper we develop a common theoretical framework for the dynamics of thin featureless interfaces. We explicitly demonstrate that the same phase field and velocity dependent one-scale models characterizing the dynamics of relativistic domain walls, in a cosmological context, can also successfully describe, in a friction dominated regime, the dynamics of nonrelativistic interfaces in a wide variety of material systems. We further show that a statistical version of von Neumann's law applies in the case of scaling relativistic interface networks, implying that, although relativistic and nonrelativistic interfaces have very different dynamics, a single simulation snapshot is not able to clearly distinguish the two regimes. We highlight that crucial information is contained in the probability distribution function for the number of edges of domains bounded by the interface network and explain why laboratory tests with nonrelativistic interfaces can be used to rule out cosmological domain walls as a significant dark energy source.

If you fall into a real astronomical black hole (choosing a supermassive black hole, to make sure that the tidal forces don't get you first), then you will probably meet your fate not at a central singularity, but rather in the exponentially growing, relativistic counter-streaming instability at the inner horizon first pointed out by Poisson & Israel (1990), who called it mass inflation. The chief purpose of this paper is to present a clear exposition of the physical cause and consequence of inflation in spherical, charged black holes. Inflation acts like a particle accelerator in that it accelerates cold ingoing and outgoing streams through each other to prodigiously high energies. Inflation feeds on itself: the acceleration is powered by the gravity produced by the streaming energy. The paper: (1) uses physical arguments to develop simple approximations that follow the evolution of inflation from ignition, through inflation itself, to collapse; (2) confirms that the simple approximations capture accurately the results of fully nonlinear one- and two-fluid self-similar models; (3) demonstrates that, counter-intuitively, the smaller the accretion rate, the more rapidly inflation exponentiates; (4) shows that in single perfect-fluid models, inflation occurs only if the sound speed equals the speed of light, supporting the physical idea that inflation in single fluids is driven by relativistic counter-streaming of waves; (5) shows that what happens during inflation up to the Planck curvature depends not on the distant past or future, but rather on events happening only a few hundred black hole crossing times into the past or future; (6) shows that, if quantum gravity does not intervene, then the generic end result of inflation is not a general relativistic null singularity, but rather a spacelike singularity at zero radius.

In this work we investigate the duality linking standard and tachyon scalar field cosmologies. We determine the transformation between standard and tachyon scalar fields and between their associated potentials, corresponding to the same background evolution. We show that, in general, the duality is broken at a perturbative level, when deviations from a homogeneous and isotropic background are taken into account. However, we find that for slow-rolling fields the duality is still preserved at a linear level. We illustrate our results with specific examples of cosmological relevance, where the correspondence between scalar and tachyon scalar field models can be calculated explicitly.

We investigate the potential role of string and monopole-type junctions in the frustration of domain wall networks using a velocity-dependent one-scale model for the characteristic velocity, $v$, and the characteristic length, $L$, of the network. We show that, except for very special network configurations, $v^2 \lsim (HL)^2 \lsim (\rho_\sigma + \rho_\mu)/\rho_m$ where $H$ is the Hubble parameter and $\rho_\sigma$, $\rho_\mu$ and $\rho_m$ are the average density of domain walls, strings and monopole-type junctions. We further show that if domain walls are to provide a significant contribution to the dark energy without generating exceedingly large CMB temperature fluctuations then, at the present time, the network must have a characteristic length $ L_0 \lsim 10 \Omega_{\sigma 0}^{-2/3} {\rm kpc}$ and a characteristic velocity $v_0 \lsim 10^{-5} \Omega_{\sigma 0}^{-2/3}$ where $\Omega_{\sigma 0}=\rho_{\sigma 0}/\rho_{c 0}$ and $\rho_c$ is the critical density. In order to satisfy these constraints with $\Omega_{\sigma 0} \sim 1$, $\rho_{m 0}$ would have to be at least 10 orders of magnitude larger than $\rho_{\sigma 0}$, which would be in complete disagreement with observations. This result provides very strong additional support for the conjecture that no natural frustration mechanism, which could lead to a significant contribution of domain walls to the dark energy budget, exists.

The role of domain wall junctions in Carter's pentahedral model is investigated both analytically and numerically. We perform, for the first time, field theory simulations of such model with various initial conditions. We confirm that there are very specific realizations of Carter's model corresponding to square lattice configurations with X-type junctions which could be stable. However, we show that more realistic realizations, consistent with causality constraints, do lead to a scaling domain wall network with Y-type junctions. We determine the network properties and discuss the corresponding cosmological implications, in particular for dark energy.

Current constraints on the dark energy equation of state parameter, $w$, are expected to be improved by more than one order of magnitude in the next decade. If $|w-1| \gsim 0.01$ around the present time, but the dark energy dynamics is sufficiently slow, it is possible that future constraints will rule out a cosmological constant while being consistent with a time-independent equation of state parameter. In this paper, we show that although models with such behavior can be constructed, they do require significant fine-tuning. Therefore, if the observed acceleration of the Universe is induced by a dark energy component, then finding $w \neq -1$ would, on its own, constitute very strong evidence for a dynamical dark energy equation of state.

In this paper we study the impact of the fractional matter density uncertainty in the reconstruction of the equation of state of dark energy. We consider both standard reconstruction methods, based on the dynamical effect that dark energy has on the expansion of the Universe, as well as non-standard methods, in which the evolution of the dark energy equation of state with redshift is inferred through the variation of fundamental couplings such as the fine structure constant, $\alpha$, or the proton-to-electron mass ratio, $\mu$. We show that the negative impact of the matter density uncertainty in the dark energy reconstruction using varying couplings may be very small compared to standard reconstruction methods. We also briefly discuss other fundamental questions which need to be answered before varying couplings can be successfully used to probe the nature of the dark energy.

This work deals with bifurcation and pattern changing in models described by two real scalar fields. We consider generic models with quartic potentials and show that the number of independent polynomial coefficients affecting the ratios between the various domain wall tensions can be reduced to 4 if the model has a superpotential. We then study specific one-parameter families of models and compute the wall tensions associated with both BPS and non-BPS sectors. We show how bifurcation can be associated to modification of the patterns of domain wall networks and illustrate this with some examples which may be relevant to describe realistic situations of current interest in high energy physics. In particular, we discuss a simple solution to the cosmological domain wall problem.

A detailed non-linear analysis of the internal structure of spherical, charged black holes that are accreting scalar matter is performed in the framework of the Brans-Dicke theory of gravity. We choose the lowest value of the Brans-Dicke parameter that is compatible with observational constraints. First, the homogeneous approximation is used. It indicates that mass inflation occurs and that the variations of the Brans-Dicke scalar inside the black hole, which could in principle be large in the absence of mass inflation, become small when mass inflation does occur. Then, a full non-linear numerical study of the black hole interior perturbed by a self-gravitating massless uncharged scalar-field is performed. We use an algorithm with adaptive mesh refinement capabilities. In this way, the changes in the internal structure of the black hole caused by mass inflation are determined, as well as the induced variations of the Brans-Dicke scalar, confirming, qualitatively, the indications given by the homogeneous approximation.

This work deals with the dynamics of inflation in the context of a scalar-vector-tensor theory of gravity exhibiting spontaneous Lorentz violation at early times. We describe a first-order formalism which we use to obtain new exact Lorentz violating inflationary solutions for a broad family of models, some in the absence of a potential for the inflaton field. Our results show that different conditions are required to solve the horizon and flatness problems. In particular, we find a necessary condition for inflation to provide a solution to both problems and we show that in inflationary models with no inflaton potential a period of superinflation might be necessary to solve the flatness problem.

We study the evolution of maximally symmetric $p$-branes with a $S_{p-i}\otimes \mathbbm{R}^i$ topology in flat expanding or collapsing homogeneous and isotropic universes with $N+1$ dimensions (with $N \ge 3$, $p < N$, $0 \le i < p$). We find the corresponding equations of motion and compute new analytical solutions for the trajectories in phase space. For a constant Hubble parameter, $H$, and $i=0$ we show that all initially static solutions with a physical radius below a certain critical value, $r_c^0$, are periodic while those with a larger initial radius become frozen in comoving coordinates at late times. We find a stationary solution with constant velocity and physical radius, $r_c$, and compute the root mean square velocity of the periodic $p$-brane solutions and the corresponding (average) equation of state of the $p$-brane gas. We also investigate the $p$-brane dynamics for $H \neq {\rm constant}$ in models where the evolution of the universe is driven by a perfect fluid with constant equation of state parameter, $w={\cal P}_p/\rho_p$, and show that a critical radius, $r_c$, can still be defined for $ -1 \le w < w_c$ with $w_c=(2-N)/N$. We further show that for $w \sim w_c$ the critical radius is given approximately by $r_c H \propto (w_c-w)^{\gamma_c}$ with $\gamma_c=-1/2$ ($r_c H \to \infty$ when $w \to w_c$). Finally, we discuss the impact that the large scale dynamics of the universe can have on the macroscopic evolution of very small loops.

We discuss the possible dynamical role of extended cosmic defects on galactic scales, specifically focusing on the possibility that they may provide the dark matter suggested by the classical problem of galactic rotation curves. We emphasize that the more standard defects (such as Goto-Nambu strings) are unsuitable for this task, but show that more general models (such as transonic wiggly strings) could in principle have a better chance. In any case, we show that observational data severely restricts any such scenarios.

We use a combination of analytic tools and an extensive set of the largest and most accurate three-dimensional field theory numerical simulations to study the dynamics of domain wall networks with junctions. We build upon our previous work and consider a class of models which, in the limit of large number $N$ of coupled scalar fields, approaches the so-called `ideal' model (in terms of its potential to lead to network frustration). We consider values of $N$ between N=2 and N=20, and a range of cosmological epochs, and we also compare this class of models with other toy models used in the past. In all cases we find compelling evidence for a gradual approach to scaling, strongly supporting our no-frustration conjecture. We also discuss the various possible types of junctions (including cases where there is a hierarchy of them) and their roles in the dynamics of the network. Finally, we revise the Zel'dovich bound and provide an updated cosmological bound on the energy scale of this type of defect network: it must be lower than $10 {\rm keV}$.

We study the evolution of biased domain walls in the early universe. We explicitly discuss the roles played by the surface tension and volume pressure in the evolution of the walls, and quantify their effects by looking at the collapse of spherical wall solutions. We then apply our results to a particular mechanism, known as the devaluation scenario, in which the dynamics of biased domain walls was suggested as a possible solution to the cosmological constant problem. Our results indicate that devaluation will in general lead to values of the cosmological constant that differ by several orders of magnitude from the observationally inferred value, $\rho^{1/4}_{vac}\sim10^{-3} \rm eV$. We also argue that the reasons behind this are not specific to a particular realization, and are expected to persist in any scenario of this kind, except if a low energy cut-off on the spectra of vacuum energy densities, of the order of the critical density at the present time, is postulated. This implies that any such scenario will require a fine-tuning similar to the usual one.

We study the cosmological evolution of the fine structure constant, $\alpha$, and the proton-to-electron mass ratio, $\mu=m_p/m_e$, in the context of a generic class of models where the gauge kinetic function is a linear function of a quintessence-type real scalar field, $\phi$, described by a Lagrangian with a standard kinetic term and a scalar field potential, $V(\phi)$. We further assume that the scalar field potential is a monotonic function of $\phi$ and that the scalar field is always rolling down the potential. We show that, for this class of models, low-redshift constrains on the evolution of $\alpha$ and $\mu$ can provide very stringent limits on the corresponding variations at high-redshift. We also demonstrate that these limits may be relaxed by considering more general models for the dynamics of $\alpha$ and $\mu$. However, in this case, the ability to reconstruct the evolution of the dark energy equation of state using varying couplings could be seriously compromised.

We revisit the paradigm of unified dark energy discussing in detail the averaging problem in this type of scenarios, highlighting the need for a full non-linear treatment. We also address the question of if and how models with one or several dark fluids can be observationally distinguished. Simpler and physically clearer derivations of some key results, most notably on the relation between the generalized Chaplygin gas and the standard ($\Lambda$CDM) `concordance' model and on a Jeans-type small-scale instability of some coupled dark energy/dark matter models are presented.

In this paper we show that a \emph{one-to-one} correspondence exists between any dark energy model and an equivalent (from a cosmological point of view, in the absence of perturbations) quartessence model in which dark matter and dark energy are described by a single perfect fluid. We further show that if the density fluctuations are small, the evolution of the sound speed squared, $c_s^2$, is fully coupled to the evolution of the scale factor and that the transition from the dark matter to the dark energy dominated epoch is faster (slower) than in a standard $\Lambda$CDM model if $c_s^2 > 0$ ($c_s^2 < 0$). In particular, we show that the mapping of the simplest quintessence scenario with constant $w_Q \equiv p_Q/ \rho_Q$ into a unified dark energy model requires $c_s^2 < 0$ at the present time (if $w_Q > -1$) contrasting to the Chaplygin gas scenario where one has $c_s^2 > 0$. However, we show that non-linear effects severely complicate the analysis, in particular rendering linear results invalid even on large cosmological scales. Although a detailed analysis of non-linear effects requires solving the full Einstein field equations, some general properties can be understood in simple terms. In particular, we find that in the context of Chaplygin gas models the transition from the dark matter to the dark energy dominated era may be anticipated with respect to linear expectations leading to a background evolution similar to that of standard $\Lambda$CDM models. On the other hand, in models with $c_s^2 > 0$ the expected transition from the decelerating to the accelerating phase may never happen.

We provide a generic but physically clear discussion of the clustering properties of dark energy models. We explicitly show that in quintessence-type models the dark energy fluctuations, on scales smaller than the Hubble radius, are of the order of the perturbations to the Newtonian gravitational potential, hence necessarily small on cosmological scales. Moreover, comparable fluctuations are associated with different gauge choices. We also demonstrate that the often used homogeneous approximation is unrealistic, and that the so-called dark energy mutation is a trivial artifact of an effective, single fluid description. Finally, we discuss the particular case where the dark energy fluid is coupled to dark matter.

We study the evolution of simple cosmic string loop solutions in an inflationary universe. We show, for the particular case of circular loops, that periodic solutions do exist in a de Sitter universe, below a critical loop radius $R_c H=1/2$. On the other hand, larger loops freeze in comoving coordinates, and we explicitly show that they can survive more $e$-foldings of inflation than point-like objects. We discuss the implications of these findings for the survival of realistic cosmic string loops during inflation, and for the general characteristics of post-inflationary cosmic string networks. We also consider the analogous solutions for domain walls, in which case the critical radius is $R_c H=2/3$.

We describe the results of the largest and most accurate three-dimensional field theory simulations of domain wall networks with junctions. We consider a previously introduced class of models which, in the limit of large number $N$ of coupled scalar fields, approaches the so-called `ideal' model (in terms of its potential to lead to network frustration). We consider values of $N$ between N=2 and N=20. In all cases we find compelling evidence for a gradual approach to scaling, with the quantitative scaling parameters having only a mild dependence on $N$. These results strongly support our no-frustration conjecture.

We revisit the idea of using varying couplings to probe the nature of dark energy, in particular by reconstructing its equation of state. We show that this method can be far superior to the standard methods (using type Ia supernovae or weak lensing). We also show that the simultaneous use of measurements of the fine-structure constant $\alpha$ and the electron-to-proton mass ratio $\mu$ allows a direct probe of grand unification scenarios. We present forecasts for the sensitivity of this method, both for the near future and for the next generation of spectrographs -- for the latter we focus on the planned CODEX instrument for ESO's Extremely Large Telescope (formerly known as OWL). A high-accuracy reconstruction of the equation of state may be possible all the way up to redshift $z\sim4$.

We study the reionization history of the Universe in cosmological models with non-Gaussian density fluctuations, taking them to have a renormalized $\chi^2$ probability distribution function parametrized by the number of degrees of freedom, $\nu$. We compute the ionization history using a simple semi-analytical model, considering various possibilities for the astrophysics of reionization. In all our models we require that reionization is completed prior to $z=6$, as required by the measurement of the Gunn--Peterson optical depth from the spectra of high-redshift quasars. We confirm previous results demonstrating that such a non-Gaussian distribution leads to a slower reionization as compared to the Gaussian case. We further show that the recent WMAP three-year measurement of the optical depth due to electron scattering, $\tau=0.09 \pm 0.03$, weakly constrains the allowed deviations from Gaussianity on the small scales relevant to reionization if a constant spectral index is assumed. We also confirm the need for a significant suppression of star formation in mini-halos, which increases dramatically as we decrease $\nu$.

We study a number of domain wall forming models where various types of defect junctions can exist. These illustrate some of the mechanisms that will determine the evolution of defect networks with junctions. Understanding these mechanisms is vital for a proper assessment of a number of cosmological scenarios: we will focus on the issue of whether or not cosmological frustrated domain wall networks can exist at all, but our results are also relevant for the dynamics of cosmic (super)strings, where junctions are expected to be ubiquitous. We also define and discuss the properties that would make up the ideal model in terms of hypothetical frustrated wall networks, and provide an explicit construction for such a model. We carry out a number of numerical simulations of the evolution of these networks, analyze and contrast their results, and discuss their implications for our no-frustration conjecture.

We discuss necessary conditions for a network of cosmic domain walls to have a chance of providing the dark energy that might explain the recent acceleration of the universe. We derive a strong bound on the curvature of the walls, which shows that viable candidate networks must be fine-tuned and non-standard. We also discuss various requirements that any stable lattice of frustrated walls must obey. We conjecture that, even though one can build (by hand) lattices that would be stable, no such lattices will ever come out of realistic domain wall forming cosmological phase transitions. We provide some simple numerical simulations that illustrate our results and correct some misconceptions in the published literature, but a detailed numerical analysis is left for a companion paper.

We study spatial variations of the fine-structure constant in the presence of static straight cosmic strings in the weak-field approximation in Einstein gravity. We work in the context of a generic Bekenstein-type model and consider a gauge kinetic function linear in the scalar field. We determine an analytical form for the scalar field and the string metric at large distances from the core. We show that the gravitational effects of $\alpha$-varying strings can be seen as a combination of the gravitational effects of global and local strings. We also verify that at large distances to the core the space-time metric is similar to that of a global string. We study the motion of test particles approaching from infinity and show that photons are scattered to infinity while massive particles are trapped in bounded trajectories. We also calculate an overall limit on the magnitude of the variation of $\alpha$ for a GUT string, by considering suitable cosmological constraints coming from the Equivalence Principle.

We study the evolution of the fine-structure constant, $\alpha$, induced by non-linear density perturbations in the context of the simplest class of quintessence models with a non-minimal coupling to the electromagnetic field, in which the two available free functions (potential and gauge kinetic function) are Taylor-expanded up to linear order. We show that the results obtained using the spherical infall model for an infinite wavelength inhomogeneity are inconsistent with the results of a local linearized gravity study and we argue in favour of the second approach. We also discuss recent claims that the value of $\alpha$ inside virialised regions could be significantly different from the background one on the basis of these findings.

We determine the abundance of primordial black holes (PBHs) formed in the context of non-gaussian models with primordial density perturbations. We consider models with a renormalized $\chi^2$ probability distribution function parametrized by the number, $\nu$, of degrees of freedom. We show that if $\nu$ is not too large then the PBH abundance will be altered by several orders of magnitude with respect to the standard gaussian result obtained in the $\nu \to \infty$ limit. We also study the dependence of the spectral index constraints on the nature of the cosmological perturbations for a power-law primordial power spectrum.

We study static magnetic monopoles in the context of varying alpha theories and show that there is a group of models for which the t'Hooft-Polyakov solution is still valid. Nevertheless, in general static magnetic monopole solutions in varying alpha theories depart from the classical t'Hooft-Polyakov solution with the electromagnetic energy concentrated inside the core seeding spatial variations of the fine structure constant. We show that Equivalence Principle constraints impose tight limits on the allowed variations of alpha induced by magnetic monopoles which confirms the difficulty to generate significant large-scale spatial variation of the fine structure constant found in previous works. This is true even in the most favorable case where magnetic monopoles are the source for these variations.

We introduce a new phenomenological one-scale model for the evolution of domain wall networks, and test it against high-resolution field theory numerical simulations. We argue that previous numerical estimates of wall velocities are inaccurate, and suggest a more accurate method of measurement. We show that the model provides an adequate approximation to the evolution of key parameters characterizing the evolution of the network. We use the model to study possible scaling solutions for domain wall networks, and discuss some of their cosmological consequences.

We study the cosmological evolution of domain wall networks in two and three spatial dimensions in the radiation and matter eras using a large number of high-resolution field theory simulations with a large dynamical range. We investigate the dependence of the uncertainty in key parameters characterising the evolution of the network on the size, dynamical range and number of spatial dimensions of the simulations and show that the analytic prediction compares well with the simulation results. We find that there is ample evidence from the simulations of a slow approach of domain wall networks towards a linear scaling solution. However, while at early times the uncertainty in the value of the scaling exponent is small enough for deviations from the scaling solution to be measured, at late times the error bars are much larger and no strong deviations from the scaling solution are found.

We study static cosmic string solutions in the context of Bekenstein-type models. We show that there is a class of models of this type for which the classical Nielsen-Olesen vortex is still a valid solution. However, in general static string solutions in Bekenstein-type models strongly depart from the standard Nielsen-Olesen solution with the electromagnetic energy concentrated along the string core seeding spatial variations of the fine structure constant, $\alpha$. We consider models with a generic gauge kinetic function and show that equivalence principle constraints impose tight limits on the allowed variations of $\alpha$ induced by string networks on cosmological scales.

We show that a solution to the the coincidence problem can be found in the context of a generic class of dark energy models with a scalar field, $\phi$, with a linear effective potential $V(\phi)$. We determine the fraction, $f$, of the total lifetime of the universe, $t_U$, which lies within the interval $[t_0-\Delta t_A,t_0+ \Delta t_A]$, where $t_0$ is the age of the universe at the present time, $\Delta t_A \equiv t_0-t_A$ and $t_A$ is the age of the universe when it starts to accelerate. We find that if we require $f$ to be larger than 0.1 (0.01) then $1+\omega_{\phi0} \gapp 2 \times 10^{-2}$ ($1 \times 10^{-3}$), where $\omega_{\phi} \equiv p_\phi/\rho_\phi$. These results depend mainly on the linearity of the scalar field potential for $-V(\phi_0) \lapp V(\phi) \lapp V(\phi_0)$ and are weakly dependent on the specific form of $V(\phi)$ outside this range. We also show that if $\omega_{\phi0}$ is close to -1 then $\omega_{\phi0}+1 \sim 1.6 ({\tilde \omega}_\phi +1)$, where ${\tilde \omega}_\phi$ is the weighted average value of $\omega_{\phi}$ in the time interval $[0,t_0]$. We independently confirm current observational constraints on this class of models which give $\omega_{\phi0} \lapp -0.6$ and $t_U \gapp 2.4 t_0$ at the $2 \sigma$ level.

We have studied the cosmological evolution of domain wall networks in two, three and four spatial dimensions using high-resolution field theory simulations. The dynamical range and number of our simulations is larger than in previous works, but does not allow us to exclude previous hints of deviations to the naively expected scale-invariant evolution. These results therefore suggest that the approach of domain wall networks to linear scaling is a much slower process than that of cosmic strings, which has been previously characterized in detail.

We consider a model for the reionization history of the Universe in which a significant fraction of the observed optical depth is a result of direct reionization by the decay products of a scaling cosmic defect network. We show that such network can make a significant contribution to the reionization history of the Universe even if its energy density is very small (the defect energy density has to be greater than about $10^{-11}$ of the background density). We compute the Cosmic Microwave Background temperature, polarization and temperature-polarization cross power spectrum and show that a contribution to the observed optical depth due to the decay products of a scaling defect network may help to reconcile a high optical depth with a low redshift of complete reionization suggested by quasar data. However, if the energy density of defects is approximately a constant fraction of the background density then these models do not explain the large scale bump in the temperature-polarization cross power spectrum observed by WMAP.

We study the simplest class of Bekenstein-type, varying $\alpha$ models, in which the two available free functions (potential and gauge kinetic function) are Taylor-expanded up to linear order. Any realistic model of this type reduces to a model in this class for a certain time interval around the present day. Nevertheless, we show that no such model is consistent with all existing observational results. We discuss possible implications of these findings, and in particular clarify the ambiguous statement (often found in the literature) that ``the Webb results are inconsistent with Oklo''.

We present up-to-date constraints on the value of the fine-structure constant at the epoch of decoupling from the recent observations made by the Wilkinson Microwave Anisotropy Probe (WMAP) satellite. In the framework of models we considered, a positive (negative) variation of the value of $\alpha$ at decoupling with respect to the present-day value is now bounded to be smaller than 2% (6%) at 95% confidence level. We point out that the existence of an early reionization epoch as suggested by the above measurements will, when more accurate cosmic microwave background polarization data is available, lead to considerably tighter constraints. The so-called `reionization bump', in particular, will be extremely useful for this purpose. We find that constraints on $\alpha$ can in principle only be as small as $0.1 %$ using CMB data alone--tighter constraints will require further (non-CMB) priors.

We present a detailed analysis of present and future Cosmic Microwave Background constraints of the value of the fine-structure constant, $\alpha$. We carry out a more detailed analysis of the WMAP first-year data, deriving state-of-the-art constraints on $\alpha$ and discussing various other issues, such as the possible hints for the running of the spectral index. We find, at 95% C.L. that $0.95 < \alpha_{\text{dec}} / \alpha_0 < 1.02$. Setting $dn_S /dlnk=0$, yields $0.94< \alpha_{\text{dec}} / \alpha_0 < 1.01$ as previously reported. We find that a lower value of $\alpha / \alpha_0$ makes a value of $d n_S /dlnk = 0$ more compatible with the data. We also perform a thorough Fisher Matrix Analysis (including both temperature and polarization, as well as $\alpha$ and the optical depth $\tau$), in order to estimate how future CMB experiments will be able to constrain $\alpha$ and other cosmological parameters. We find that Planck data alone can constrain $\tau$ with a accuracy of the order 4% and that this constraint can be as small as 1.7% for an ideal cosmic variance limited experiment. Constraints on $\alpha$ are of the order 0.3% for Planck and can in principle be as small as 0.1% using CMB data alone - tighter constraints will require further (non-CMB) priors.

We present a summary of recent constraints on the value of the fine-structure constant at the epoch of decoupling from the recent observations made by the Wilkinson Microwave Anisotropy Probe (WMAP) satellite. Within the set of models considered, a variation of the value of $\alpha$ at decoupling with respect to the present-day value is now bounded to be smaller than 2% (6%) at 95% confidence level. We point out that the existence of an early reionization epoch as suggested by the above measurements will, when more accurate cosmic microwave background polarization data is available, lead to considerably tighter constraints. We find that the tightest possible constraint on $\alpha$ is about $0.1 %$ using CMB data alone--tighter constraints will require further (non-CMB) priors.

We study the behaviour of cosmic string networks in contracting universes, and discuss some of their possible consequences. We note that there is a fundamental time asymmetry between defect network evolution for an expanding universe and a contracting universe. A string network with negligible loop production and small-scale structure will asymptotically behave during the collapse phase as a radiation fluid. In realistic networks these two effects are important, making this solution only approximate. We derive new scaling solutions describing this effect, and test them against high-resolution numerical simulations. A string network in a contracting universe, together with the gravitational radiation background it has generated, can significantly affect the dynamics of the universe both locally and globally. The network can be an important source of radiation, entropy and inhomogeneity. We discuss the possible implications of these findings for bouncing and cyclic cosmological models.

We explicitly demonstrate that, contrary to recent claims, the dynamics of a generalized Chaplygin gas model with an equation of state $p=-C$ (where $C$ is a positive constant) is equivalent to that of a standard $\Lambda$CDM model to first order in the metric perturbations. We further argue that the analogy between the two models goes well beyond linear theory and conclude that they cannot be distinguished based on gravity alone.

We discuss the onset of the non-linear regime in the context of unified dark matter models involving a generalised Chaplygin gas. We show that the transition from dark matter-like to dark energy-like behaviour will never be smooth. In some regions of space the transition will never take place while in others it may happen sooner or later than naively expected. As a result the linear theory used in previous studies may break down late in the matter dominated era even on large cosmological scales. We study the importance of this effect showing that its magnitude depends on the exact form of the equation of state in the low density regime. We expect that our results will be relevant for other unified dark matter scenarios particularly those where the quartessence candidate is a perfect fluid.

We derive the equations governing the dynamics of cosmic strings in a flat anisotropic universe of Bianchi type I and study the evolution of simple cosmic string loop solutions. We show that the anisotropy of the background can have a characteristic effect in the loop motion. We discuss some cosmological consequences of these findings and, by extrapolating our results to cosmic string networks, we comment on their ability to survive an inflationary epoch, and hence be a possible fossil remnant (still visible today) of an anisotropic phase in the very early universe.

We investigate the dependence of the epoch of reionization on the properties of cosmological perturbations, in the context of cosmologies permitted by WMAP. We compute the redshift of reionization using a simple model based on the Press-Schechter approximation. For a power-law initial spectrum we estimate that reionization is likely to occur at a redshift $z_{reion} = 17^{+10}_{-7}$, consistent with the WMAP determination based on the temperature-polarization cross power spectrum. We estimate the delay in reionization if there is a negative running of the spectral index, as weakly indicated by WMAP. We then investigate the dependence of the reionization redshift on the nature of the initial perturbations. We consider chi-squared probability distribution functions with various degrees of freedom, motivated both by non-standard inflationary scenarios and by defect models. We find that in these models reionization is likely occur much earlier, and to be a slower process, than in the case of initial gaussian fluctuations. We also consider a hybrid model in which cosmic strings make an important contribution to the seed fluctuations on scales relevant for reionization. We find that in order for that model to agree with the latest WMAP results, the string contribution to the matter power spectrum on the standard $8 h^{-1} Mpc$ scale is likely to be at most at the level of one percent, which imposes tight constraints on the value of the string mass per unit length.

We discuss the importance of including baryons in analyses of unified dark matter scenarios, focusing on toy models involving a generalized Chaplygin gas. We determine observational constraints on this unified dark matter scenario coming from large scale structure, type Ia Supernovae and CMB data showing how this component can bring about a different behaviour from classical $\Lambda$CDM and thus motivate further studies of this type of models. We also speculate on interesting new features which are likely to be important on non-linear scales in this context.

We discuss the issue of toy model building for the dark energy component of the universe. Specifically, we consider two generic toy models recently proposed as alternatives to quintessence models, known as Cardassian expansion and the Chaplygin gas. We show that the former is enteriely equivalent to a class of quintessence models. We determine the observational constraints on the latter, coming from recent supernovae results and from the shape of the matter power spectrum. As expected, these restrict the model to a behaviour that closely matches that of a standard cosmological constant $\Lambda$.

We introduce the notion of the cosmic numbers of a cosmological model, and discuss how they can be used to naturally classify models according to their ability to solve some of the problems of the standard cosmological model.

We study the behaviour and consequences of cosmic string networks in contracting universes. They approximately behave during the collapse phase as a radiation fluids. Scaling solutions describing this are derived and tested against high-resolution numerical simulations. A string network in a contracting universe, together with the gravitational radiation it generates, can affect the dynamics of the universe both locally and globally, and be an important source of radiation, entropy and inhomogeneity. We discuss possible implications for bouncing and cyclic models.

We extend our recent work on the effects of a time-varying fine-structure constant $\alpha$ in the cosmic microwave background, by providing a thorough analysis of the degeneracies between $\alpha$ and the other cosmological parameters, and discussing ways to break these with both existing and/or forthcoming data. In particular, we present the state-of-the-art CMB constraints on $\alpha$, through a combined analysis of the BOOMERanG, MAXIMA and DASI datasets. We also present a novel discussion of the constraints on $\alpha$ coming from large-scale structure observations, focusing in particular on the power spectrum from the 2dF survey. Our results are consistent with no variation in $\alpha$ from the epoch of recombination to the present day, and restrict any such (relative) variation to be less than about 4%. We show that the forthcoming MAP and (particularly) Planck experiments will be able to break most of the currently existing degeneracies between $\alpha$ and other parameters, and measure $\alpha$ to better than percent accuracy.

The median statistic has recently been discussed by Gott \textit{et al.} as a more reliable alternative to the standard $\chi^2$ likelihood analysis, in the sense of requiring fewer assumptions about the data and being almost as constraining. We apply this statistic to the currently available combined dataset of 92 distant type Ia supernovae, and also to a mock SNAP-class dataset. We find that the performances of the modified median and $\chi^2$ statistics are comparable, particularly in the latter case. We further extend the work of Gott \textit{et al.} by modifying the median statistic to account for the number and size of sequences of consecutive points above or below the median. We also comment on how the performance of the statistic depends on the choice of free parameters that one is estimating.

We consider a `brane-world scenario' recently introduced by Dvali, Gabadadze and Porrati, and subsequently proposed as an alternative to a cosmological constant in explaining the current acceleration of the universe. We show that, contrary to these claims, this proposal is already strongly disfavoured by the available Type Ia Supernovae, Cosmic Microwave Background and cluster data.

We show that our Universe may be inhomogeneous on large sub-horizon scales without us being able to realise it. We assume that a network of domain walls permeates the universe dividing it in domains with slightly different vacuum energy densities. We require that the energy scale of the phase transition which produced the domain walls is sufficiently low so that the walls have a negligible effect on structure formation. Nevertheless, the different vacuum densities of different domains will lead to different values of the cosmological parameters $\Omega_\Lambda^0$, $\Omega_m^0$ and $h$, in each patch thus affecting the growth of cosmological perturbations at recent times. Hence, if our local patch of the universe (with uniform vacuum density) is big enough -- which is likely to happen given that we should have on average about one domain per horizon volume -- we might not notice these large-scale inhomogeneities. This happens because in order to see a patch with a different vacuum density one may have to look back at a time when the universe was still very homogeneous.

We extend our recent work on ``two-metric'' theories of gravity by showing how in such models cosmic defects can produce a spectrum of primordial Gaussian density perturbations. This will happen when the speed characterising the decay products of the defect network is much larger than the speed characterising gravity and all standard model particles. This model will exactly mimic all the standard predictions of inflationary models, and the only way of distinguishing the two will be via the detection of the decay products of the network.

In a recent article, a simple `spherical bubble' toy model for a spatially varying vacuum energy density was introduced, and type Ia supernovae data was used to constrain it. Here we generalize the model to allow for the fact that we may not necessarily be at the centre of a region with a given set of cosmological parameters, and discuss the constraints on these models coming from Cosmic Microwave Background Radiation data. We find tight constraints on possible spatial variations of the vacuum energy density for any significant deviations from the centre of the bubble and we comment on the relevance of our results.

Higher-dimensional theories have the remarkable feature of predicting a time (and hence redshift) dependence of the `fundamental' four dimensional constants on cosmological timescales. In this paper we update the bounds on a possible variation of the fine structure constant alpha at the time of BBN (z =10^10) and CMB (z=10^3). Using the recently-released high-resolution CMB anisotropy data and the latest estimates of primordial abundances of 4He and D, we do not find evidence for a varying alpha at more than one-sigma level at either epoch.

We consider a very simple toy model for a spatially varying `cosmological constant', where we are inside a spherical bubble (with a given set of cosmological parameters) that is surrounded by a larger region where these parameters are different. This model includes essential features of more realistic scenarios with a minimum number of parameters. We calculate the luminosity distance in the presence of spatial variations of the vacuum energy density using linear perturbation theory and discuss the use of type Ia supernovae to impose constraints on this type of models. We find that presently available observations are only constraining at very low redshifts, but also provide independent confirmation that the high red-shift supernovae data does prefer a relatively large positive cosmological constant.

We describe the detailed study and results of high-resolution numerical simulations of string-induced structure formation in open universes and those with a non-zero cosmological constant. The effect from small loops generated from the string network has also been investigated. We provide a semi-analytical model which can reproduce these simulation results. A detailed study of cosmic string network properties regarding structure formation is also given, including the correlation time, the topological analysis of the source spectrum, the correlation between long strings and loops, and the evolution of long-string and loop energy densities. For models with $\Gamma=\Omega h=0.1--0.2 and a cold dark matter background, we show that the linear density fluctuation power spectrum induced by cosmic strings has both an amplitude at $8 h^{-1}$Mpc, $\sigma_8$, and an overall shape which are consistent within uncertainties with those currently inferred from galaxy surveys. The cosmic string scenario with hot dark matter requires a strongly scale-dependent bias in order to agree with observations.

We re-analyze the question of the use of cosmological observations to infer the present state and future evolution of our patch of the universe. In particular, we discuss under which conditions one might be able to infer that our patch will enter an inflationary stage, as a prima facie interpretation of the Type Ia supernovae and CMB data would suggest. We then establish a `physical' criterion for the existence of inflation, to be contrasted with the more `mathematical' one recently proposed by Starkman et al. \cite{STV}.

We consider the evolution of domain walls produced during an anisotropic phase in the very early universe, showing that the resulting network can be very anisotropic. If the domain walls are produced during an inflationary era, the network will soon freeze out in comoving coordinates retaining the imprints of the anisotropic regime, even though inflation makes the universe isotropic. Only at late times, when the typical size of the major axis of the domain walls becomes smaller than the Hubble radius, does the network evolve rapidly towards isotropy. Hence, we may hope to see imprints of the anisotropic era if by today the typical size of the major axis of the domain walls is of the order of the Hubble radius, or if the walls re-entered it only very recently. Depending on the mass scale of the domain walls, there is also the possibility that they re-enter at earlier times, but their evolution remains friction-dominated until recently, in which case the signatures of the anisotropic era will be much better preserved. These effects are expected to occur in generic topological defect models.

We perform a likelihood analysis of the recently released BOOMERanG and MAXIMA data, allowing for the possibility of a time-varying fine-structure constant. We find that in general this data prefers a value of $\alpha$ that was smaller in the past (which is in agreement with measurements of $\alpha$ from quasar observations). However, there are some interesting degeneracies in the problem which imply that strong statements about $\alpha$ can not be made using this method until independent accurate determinations of $\Omega_b h^2$ and $H_0$ are available. We also show that a preferred lower value of $\alpha$ comes mainly from the data points around the first Doppler peak, whereas the main effect of the high-$\ell$ data points is to increase the preferred value for $\Omega_b h^2$ (while also tightening the constraints on $\Omega_0$ and $H_0$). We comment on some implications of our results.

We point out that in the context of ``two-metric'' theories of gravity there is the possibility that cosmic defects will produce a spectrum of primordial adiabatic density perturbations. This will happen when the speed characterising the defect-producing scalar field is much larger than the speed characterising gravity and all standard model particles. This model will exactly mimic the standard predictions of inflationary models, with the exception of a small non-Gaussian signal which could be detected by future experiments. We briefly discuss defect evolution in these scenarios and analyze their cosmological consequences.

We discuss constraints on so-called `varying speed of light theories' coming from the physics of the early universe. Specifically, we consider the position of the first acoustic peak of the CMB angular power spectrum, showing that the recent determination of its position by various CMB experiments, in particular BOOMERanG/NA, can be used to constrain the value of the speed of light $c$ after the epoch of last scattering. Specifically, we find that for a flat universe and a standard density of baryonic matter a variation in $c$ of up to 4% is consistent with the current observations. The position of the Doppler peak is fairly sensitive to changes in $c$, and future observations should dramatically improve this bound. On the other hand, we also find that the maximum baryonic density allowed in VSL theories is about $\Omega_B h^2\sim 0.06\Omega_0$. We comment on the significance of these results.

The cloud-in-cloud problem is studied in the context of the extension to non-Gaussian density fields of the Press-Schechter approach for the calculation of the mass function. As an example of a non-Gaussian probability distribution functions (PDFs) we consider the Chi-square, with various degrees of freedom. We generate density fields in cubic boxes with periodic boundary conditions and then determine the number of points considered collapsed at each scale through an hierarchy of smoothing windows. We find that the mass function we obtain differs from that predicted using the Extended Press-Schechter formalism, particularly for low values of $\sigma$ and for those PDFs most distinct from a Gaussian.

We use the present observed number density of large X-ray clusters to constrain the amplitude of matter density perturbations induced by cosmic strings on the scale of $8 h^{-1}$Mpc ($\sigma_8$), in both open cosmologies and flat models with a non-zero cosmological constant. We find a slightly lower value of $\sigma_8$ than that obtained in the context of primordial Gaussian fluctuations generated during inflation. This lower normalization of $\sigma_8$ results from the mild non-Gaussianity on cluster scales, where the one point probability distribution function is well approximated by a $\chi^2$ distribution. We use our estimate of $\sigma_8$ to constrain the string linear energy density $\mu$ and show that it is consistent with the COBE normalization.

We propose a new generalisation of general relativity which incorporates a variation in both the speed of light in vacuum (c) and the gravitational constant (G) and which is both covariant and Lorentz invariant. We solve the generalised Einstein equations for Friedmann universes and show that arbitrary time-variations of c and G never lead to a solution to the flatness, horizon or $\Lambda$ problems for a theory satisfying the strong energy condition. In order to do so, one needs to construct a theory which does not reduce to the standard one for any choice of time, length and energy units. This can be achieved by breaking a number of invariance principles such as covariance and Lorentz invariance.

We investigate the contribution made by small loops from a cosmic string network as seeds for large-scale structure formation. We show that cosmic string loops are highly correlated with the long-string network on large scales and therefore contribute significantly to the power spectrum of density perturbations if the average loop lifetime is comparable to or above one Hubble time. This effect further improves the large-scale bias problem previously identified in earlier studies of cosmic string models.

We study the growth of linear perturbations induced by a generic causal scaling source as a function of the cosmological parameters $h$, $\Omega^m_0$ and $\Omega^\Lambda_0$. We show that for wavenumbers $k \gsim 0.01 h/Mpc$ the spectrum of density and velocity perturbations scale in a similar way to that found in inflationary models with primordial perturbations. We show that this result is independent of the more or less incoherent nature of the source, the small scale power spectrum of the source and of deviations from scaling which naturally occur at late times if $\Omega^m_0 \neq 1$.

We present a study of open inflation cosmological scenarios in which cosmic strings form betwen the two inflationary epochs. It is shown that in these models strings are stretched outside the horizon due to the inflationary expansion but must necessarily re-enter the horizon before the epoch of equal matter and radiation densities. We determine the power spectrum of cold dark matter perturbations in these hybrid models, finding good agreement with observations for values of $\Gamma=\Omega_0h\sim0.3$ and comparable contributions from the active and passive sources to the CMB. Finally, we briefly discuss other cosmological consequences of these models.

We describe the results of high-resolution numerical simulations of string-induced structure formation in open universes and those with a non-zero cosmological constant. For models with $\Gamma = \Omega h = 0.1 - 0.2$ and a cold dark matter background, we show that the linear density fluctuation power spectrum has both an amplitude at $8 h^{-1} Mpc, \sigma_8$, and an overall shape which are consistent within uncertainties with those currently inferred from galaxy surveys. The cosmic string scenario with hot dark matter requires a strongly scale-dependent bias in order to agree with observations.

We investigate the non-gaussian properties of cosmic-string-seeded linear density perturbations with cold and hot dark matter backgrounds, using high-resolution numerical simulations. We compute the one-point probability density function of the resulting density field, its skewness, kurtosis, and genus curves for different smoothing scales. A semi-analytic model is then invoked to provide a physical interpretation of our results. We conclude that on scales smaller than $1.5{(\Omega h^2)}^{-1}$Mpc, perturbations seeded by cosmic strings are very non-gaussian. These scales may still be in a linear or mildly non-linear regime in an open or $\Lambda$-universe with $\Gamma=\Omega h \lsim 0.2$.

We compute the linear power spectrum of cosmic string sedeed fluctuations in the context of neutrinos with a strong self-interaction and show that it is very similar to that obtained in the context of `normal' neutrinos. We compare our results with observational data and show that for any value of the cosmological parameters $h$ and $\Omega_0$ the interacting hot dark matter power spectrum requires a large scale dependent biasing parameter.

The cosmic string scenario in an open universe is developed -- including the equations of motion, a model of network evolution, the large angular scale CMB temperature anisotropy, and the power spectrum of density fluctuations produced by cosmic strings with dark matter. First we derive the equations of motion for cosmic string in an open FRW space-time and construct a quantitative model of the evolution of the gross features of a cosmic string network. Second, we apply this model of network evolution to estimate the rms CMB temperature anisotropy induced by cosmic strings, obtaining the normalization for the mass per unit length $\mu$ as a function of $\Omega$. Third, we consider the effects of the network evolution and normalization in an open universe on the large scale structure formation scenarios with either cold or hot dark matter.

We examine variations in the equation of state of the cosmic string portion of the cosmological fluid which lead to perturbations of the background matter density. These fluctuations in the equation of state are due to variations in the local density of cosmic string loops and gravitational radiation. Constructing a crude model of the distribution of entropy perturbations, we obtain the resulting fluctuation spectrum using a gauge-invariant formalism. We compute the resulting cosmic microwave background anisotropy, and estimate the effect of these perturbations on the cosmic string structure formation scenario.

We study the efficiency of genus statistics in differentiating between different models of structure formation. Simple models which reproduce the salient features of the structure seeded by topological defects are examined. We consider accretion onto static point masses, modeling slow-moving cosmic string loops or other primordial point-like sources. Filamentary structures and wakes are considered as models of the structures seeded by slow and fast moving string, respectively. Comparison is made with predictions of genus statistics for Gaussian fluctuations and with genus curves obtained by the CfA survey. A generic class of density models with wakes and filaments is found to provide results comparable or better than Gaussian models for this suite of tests.