Infinity, in various guises, has been invoked recently in order to `explain' a number of important questions regarding observable phenomena in science, and in particular in cosmology. Such explanations are by their nature speculative. Here we introduce the notions of relative infinity, closure, and economy of explanation and ask: to what extent explanations involving relative or real constructed infinities can be treated as reasonable?

We study the effects of inhomogeneities on the evolution of the Universe, by considering a range of cosmological models with discretized matter content. This is done using exact and fully relativistic methods that exploit the symmetries in and about submanifolds of spacetimes that themselves possess no continuous global symmetries. These methods allow us to follow the evolution of our models throughout their entire history, far beyond what has previously been possible. We find that while some space-like curves collapse to anisotropic singularities in finite time, others remain non-singular forever. The resulting picture is of a cosmological spacetime in which some behaviour remains close to Friedmann-like, while other behaviours deviate radically. In particular, we find that large-scale acceleration is possible without any violation of the energy conditions.

An important open question in cosmology is the degree to which the Friedmann-Lemaitre-Robertson-Walker (FLRW) solutions of Einstein's equations are able to model the large-scale behaviour of the locally inhomogeneous observable universe. We investigate this problem by considering a range of exact n-body solutions of Einstein's constraint equations. These solutions contain discrete masses, and so allow arbitrarily large density contrasts to be modelled. We restrict our study to regularly arranged distributions of masses in topological 3-spheres. This has the benefit of allowing straightforward comparisons to be made with FLRW solutions, as both spacetimes admit a discrete group of symmetries. It also provides a time-symmetric hypersurface at the moment of maximum expansion that allows the constraint equations to be solved exactly. We find that when all the mass in the universe is condensed into a small number of objects (<10) then the amount of backreaction in dust models can be large, with O(1) deviations from the predictions of the corresponding FLRW solutions. When the number of masses is large (>100), however, then our measures of backreaction become small (<1%). This result does not rely on any averaging procedures, which are notoriously hard to define uniquely in general relativity, and so provides (to the best of our knowledge) the first exact and unambiguous demonstration of backreaction in general relativistic cosmological modelling. Discrete models such as these can therefore be used as laboratories to test ideas about backreaction that could be applied in more complicated and realistic settings.

Planck data has not found the 'smoking gun' of non-Gaussianity that would have necessitated consideration of inflationary models beyond the simplest canonical single field scenarios. This raises the important question of what these results do imply for more general models, and in particular, multi-field inflation. In this paper we revisit four ways in which two-field scenarios can behave differently from single field models; two-field slow-roll dynamics, curvaton-type behaviour, inflation ending on an inhomogeneous hypersurface and modulated reheating. We study the constraints that Planck data puts on these classes of behaviour, focusing on the latter two which have been least studied in the recent literature. We show that these latter classes are almost equivalent, and extend their previous analyses by accounting for arbitrary evolution of the isocurvature mode which, in particular, places important limits on the Gaussian curvature of the reheating hypersurface. In general, however, we find that Planck bispectrum results only constrain certain regions of parameter space, leading us to conclude that inflation sourced by more than one scalar field remains an important possibility.

We propose a thermodynamically motivated measure of gravitational entropy based on the Bel-Robinson tensor, which has a natural interpretation as the effective super-energy-momentum tensor of free gravitational fields. The specific form of this measure differs depending on whether the gravitational field is Coulomb-like or wave-like, and reduces to the Bekenstein-Hawking value when integrated over the interior of a Schwarzschild black hole. For scalar perturbations of a Robertson-Walker geometry we find that the entropy goes like the Hubble weighted anisotropy of the gravitational field, and therefore increases as structure formation occurs. This is in keeping with our expectations for the behaviour of gravitational entropy in cosmology, and provides a thermodynamically motivated arrow of time for cosmological solutions of Einstein's field equations. It is also in keeping with Penrose's Weyl curvature hypothesis.

We study inflationary perturbations in multiple-field models, for which zeta typically evolves until all isocurvature modes decay--the "adiabatic limit". We use numerical methods to explore the sensitivity of the nonlinear parameter fNL to the process by which this limit is achieved, finding an appreciable dependence on model-specific data such as the time at which slow-roll breaks down or the timescale of reheating. In models with a sum-separable potential where the isocurvature modes decay before the end of the slow-roll phase we give an analytic criterion for the asymptotic value of fNL to be large. Other examples can be constructed using a waterfall field to terminate inflation while fNL is transiently large, caused by descent from a ridge or convergence into a valley. We show that these two types of evolution are distinguished by the sign of the bispectrum, and give approximate expressions for the peak fNL.

We compute the covariant three-point function near horizon-crossing for a system of slowly-rolling scalar fields during an inflationary epoch, allowing for an arbitrary field-space metric. We show explicitly how to compute its subsequent evolution using a covariantized version of the separate universe or "delta-N" expansion, which must be augmented by terms measuring curvature of the field-space manifold, and give the nonlinear gauge transformation to the comoving curvature perturbation. Nonlinearities induced by the field-space curvature terms are a new and potentially significant source of non-Gaussianity. We show how inflationary models with non-minimal coupling to the spacetime Ricci scalar can be accommodated within this framework. This yields a simple toolkit allowing the bispectrum to be computed in models with non-negligible field-space curvature.

We calculate the conditions required to produce a large local trispectrum during two-field slow-roll inflation. This is done by extending and simplifying the 'heatmap' approach developed by Byrnes et al. The conditions required to generate a large trispectrum are broadly the same as those that can produce a large bispectrum. We derive a simple relation between tauNL and fNL for models with separable potentials, and furthermore show that gNL and tauNL can be related in specific circumstances. Additionally, we interpret the heatmaps dynamically, showing how they can be used as qualitative tools to understand the evolution of non-Gaussianity during inflation. We also show how fNL, tauNL and gNL are sourced by generic shapes in the inflationary potential, namely ridges, valleys and inflection points.

In a Universe with a detectable nontrivial spatial topology the last scattering surface contains pairs of matching circles with the same distribution of temperature fluctuations --- the so-called circles-in-the-sky. Searches undertaken for nearly antipodal pairs of such circles in cosmic microwave background maps have so far been unsuccessful. Previously we had shown that the negative outcome of such searches, if confirmed, should in principle be sufficient to exclude a detectable non-trivial spatial topology for most observers in very nearly flat ($0<\mid\Omega_{\text{tot}}-1\mid \lesssim10^{-5}$) (curved) universes. More recently, however, we have shown that this picture is fundamentally changed if the universe turns out to be {\it exactly} flat. In this case there are many potential pairs of circles with large deviations from antipodicity that have not yet been probed by existing searches. Here we study under what conditions the detection of a single pair of circles-in-the-sky can be used to uniquely specify the topology and the geometry of the spatial section of the Universe. We show that from the detection of a \emph{single} pair of matching circles one can infer whether the spatial geometry is flat or not, and if so we show how to determine the topology (apart from one case) of the Universe using this information. An important additional outcome of our results is that the dimensionality of the circles-in-the-sky parameter space that needs to be spanned in searches for matching pair of circles is reduced from six to five degrees of freedom, with a significant reduction in the necessary computational time.

We study the evolution of non-Gaussianity in multiple-field inflationary models, focusing on three fundamental questions: (a) How is the sign and peak magnitude of the non-linearity parameter fNL related to generic features in the inflationary potential? (b) How sensitive is fNL to the process by which an adiabatic limit is reached, where the curvature perturbation becomes conserved? (c) For a given model, what is the appropriate tool -- analytic or numerical -- to calculate fNL at the adiabatic limit? We summarise recent results obtained by the authors and further elucidate them by considering an inflection point model.

We study chaotic inflation in the context of modified gravitational theories. Our analysis covers models based on (i) a field coupling $\omega(\phi)$ with the kinetic energy $X$ and a nonmimimal coupling $\zeta \phi^{2} R/2$ with a Ricci scalar $R$, (ii) Brans-Dicke (BD) theories, (iii) Gauss-Bonnet (GB) gravity, and (iv) gravity with a Galileon correction. Dilatonic coupling with the kinetic energy and/or negative nonminimal coupling are shown to lead to compatibility with observations of the Cosmic Microwave Background (CMB) temperature anisotropies for the self-coupling inflaton potential $V(\phi)=\lambda \phi^{4}/4$. BD theory with a quadratic inflaton potential, which covers Starobinsky's $f(R)$ model $f(R)=R+R^{2}/(6M^{2})$ with the BD parameter $\omega_{BD}=0$, gives rise to a smaller tensor-to-scalar ratio for decreasing $\omega_{BD}$. In the presence of a GB term coupled to the field $\phi$, we express the scalar/tensor spectral indices $n_{s}$ and $n_{t}$ as well as the tensor-to-scalar ratio $r$ in terms of two slow-roll parameters and place bounds on the strength of the GB coupling from the joint data analysis of WMAP 7yr combined with other observations. We also study the Galileon-like self-interaction $\Phi(\phi) X \square\phi$ with exponential coupling $\Phi(\phi) \propto e^{\mu\phi}$. Using a CMB likelihood analysis we put bounds on the strength of the Galileon coupling and show that the self coupling potential can in fact be made compatible with observations in the presence of the exponential coupling with $\mu>0$.

We study the ability of magnetic helicity expulsion to alleviate catastrophic $\alpha$-quenching in mean field dynamos in two--dimensional spherical wedge domains. Motivated by the physical state of the outer regions of the Sun, we consider $\alpha^2\Omega$ mean field models with a dynamical $\alpha$ quenching. We include two mechanisms which have the potential to facilitate helicity expulsion, namely advection by a mean flow ("solar wind") and meridional circulation. We find that a wind alone can prevent catastrophic quenching, with the field saturating at finite amplitude. In certain parameter ranges, the presence of a large-scale meridional circulation can reinforce this alleviation. However, the saturated field strengths are typically below the equipartition field strength. We discuss possible mechanisms that might increase the saturated field.

We study the dynamics of states perturbatively expanded about a harmonic system of loop quantum cosmology, exhibiting a bounce. In particular, the evolution equations for the first and second order moments of the system are analyzed. These moments back-react on the trajectories of the expectation values of the state and hence alter the energy density at the bounce. This analysis is performed for isotropic loop quantum cosmology coupled to a scalar field with a small but non-zero constant potential, hence in a regime in which the kinetic energy of matter dominates. Analytic restrictions on the existence of dynamical coherent states and the meaning of semi-classicality within these systems are discussed. A numerical investigation of the trajectories of states that remain semi-classical across the bounce demonstrates that, at least for such states, the bounce persists and that its properties are similar to the standard case, in which the moments of the states are entirely neglected. However the bounce density does change, implying that a quantum bounce may not be guaranteed to happen when the potential is no longer negligible.

This paper has been withdrawn by the authors due to numerical problems to get viable results.

An important, and potentially detectable, signature of a non-trivial topology for the universe is the presence of so called circles-in-the-sky in the cosmic microwave background (CMB). Recent searches, confined to antipodal and nearly antipodal circles, have however failed to detect any. This outcome, coupled with recent theoretical results concerning the detectability of very nearly flat universes, is sufficient to exclude a detectable non-trivial cosmic topology for most observers in the inflationary limit ($0< |\Omega_{tot}-1| \lesssim 10^{-5}$). In a recent paper we have studied the consequences of these searches for circles if the Universe turns out to be exactly flat ($\Omega_{tot} = 1 $) as is often assumed. More specifically, we have derived the maximum angles of deviation possible from antipodicity of pairs of matching circles associated with the shortest closed geodesic for all multiply-connected flat orientable $3$-manifolds. These upper bounds on the deviation from antipodicity demonstrate that in a flat universe for some classes of topology there remains a substantial fraction of observers for whom the deviation from antipodicity of the matching circles is considerably larger than zero, which implies that the searches for circles-in-the-sky undertaken so far are not enough to exclude the possibility of a detectable non-trivial flat topology. Here we briefly review these results and discuss their consequences in the search for circles-in-the-sky in a flat universes.

We present direct numerical simulations of the equations of compressible magnetohydrodynamics in a wedge-shaped spherical shell, without shear, but with random helical forcing which has negative (positive) helicity in the northern (southern) hemisphere. We find a large-scale magnetic field that is nearly uniform in the azimuthal direction and approximately antisymmetric about the equator. Furthermore, the large-scale field in each hemisphere oscillates on nearly dynamical time scales with reversals of polarity and equatorward migration. Corresponding mean-field models also show similar migratory oscillations with a frequency that is nearly independent of the magnetic Reynolds number. This mechanism may be relevant for understanding equatorward migration seen in the solar dynamo.

Self-consistent convective dynamo simulations in wedge-shaped spherical shells are presented. Differential rotation is generated by the interaction of convection with rotation. Equatorward acceleration and dynamo action are obtained only for sufficiently rapid rotation. The angular velocity tends to be constant along cylinders. Oscillatory large-scale fields are found to migrate in the poleward direction. Comparison with earlier simulations in full spherical shells and Cartesian domains is made.

[Abridged] In a Universe with a detectable nontrivial spatial topology the last scattering surface contains pairs of matching circles with the same distribution of temperature fluctuations - the so-called circles-in-the-sky. Searches for nearly antipodal circles in maps of cosmic microwave background have so far been unsuccessful. This negative outcome along with recent theoretical results concerning the detectability of nearly flat compact topologies is sufficient to exclude a detectable nontrivial topology for most observers in very nearly flat positively and negatively curved Universes ($0<|\Omega_{tot}-1| \lesssim 10^{-5}$). Here we investigate the consequences of these searches for observable nontrivial topologies if the Universe turns out to be exactly flat ($\Omega_{tot}=1$). We demonstrate that in this case the conclusions deduced from such searches can be radically different. We show that for all multiply-connected orientable flat manifolds it is possible to directly study the action of the holonomies in order to obtain a general upper bound on the angle that characterizes the deviation from antipodicity of pairs of matching circles associated with the shortest closed geodesic. This bound is valid for all observers and all possible values of the compactification length parameters. We also show that in a flat Universe there are observers for whom the circles-in-the-sky searches already undertaken are insufficient to exclude the possibility of a detectable nontrivial spatial topology. It is remarkable how such small variations in the spatial curvature of the Universe, which are effectively indistinguishable geometrically, can have such a drastic effect on the detectability of cosmic topology.

We use direct numerical simulations of forced MHD turbulence with a forcing function that produces two different signs of kinetic helicity in the upper and lower parts of the domain. We show that the mean flux of magnetic helicity from the small-scale field between the two parts of the domain can be described by a Fickian diffusion law with a diffusion coefficient that is approximately independent of the magnetic Reynolds number and about one third of the estimated turbulent magnetic diffusivity. The data suggest that the turbulent diffusive magnetic helicity flux can only be expected to alleviate catastrophic quenching at Reynolds numbers of more than several thousands. We further calculate the magnetic helicity density and its flux in the domain for three different gauges. We consider the Weyl gauge, in which the electrostatic potential vanishes, the pseudo-Lorenz gauge, where the speed of light is replaced by the sound speed, and the `resistive gauge' in which the Laplacian of the magnetic vector potential acts as resistive term. We find that, in the statistically steady state, the time-averaged magnetic helicity density and the magnetic helicity flux are the same in all three gauges.

The Einstein static universe has played a central role in a number of emergent scenarios recently put forward to deal with the singular origin of the standard cosmological model. Here we study the existence and stability of the Einstein static solution in presence of vacuum energy corresponding to conformally-invariant fields. We show that the presence of vacuum energy stabilizes this solution by changing it to a centre equilibrium point, which is cyclically stable. This allows non-singular emergent cosmological models to be constructed in which initially the universe oscillates indefinitely about an initial Einstein static solution and is thus past eternal.

Standard cosmology poses a number of important questions. Apart from its singular origin, it possesses early and late accelerating phases required to account for observations. The vacuum energy has been considered as a possible way to resolve some of these questions. The vacuum energy density induced by free fields in an early de Sitter phase has earlier been estimated to be proportional to $H^4$, while more recently it has been suggested that the QCD condensate induces a term proportional to H at late times. These results have been employed in models which are non-singular and inflationary at early times and accelerating at late times. Here we cast these models in terms of scalar fields and study the corresponding spectrum of primordial perturbations. At early times the spectrum is found to be not scale-invariant, thus implying that slow-roll inflation is still required after the phase transition induced by the vacuum. At late times the corresponding scalar-field potential is harmonic, with a mass of the order of the Hubble scale, a result that may be understood in the light of the holographic conjecture.

We use three-dimensional direct numerical simulations of the helically forced magnetohydrodynamic equations in spherical shell segments in order to study the effects of changes in the geometrical shape and size of the domain on the growth and saturation of large-scale magnetic fields. We inject kinetic energy along with kinetic helicity in spherical domains via helical forcing using Chandrasekhar-Kendall functions. We take perfect conductor boundary conditions for the magnetic field to ensure that no magnetic helicity escapes the domain boundaries. We find dynamo action giving rise to magnetic fields at scales larger than the characteristic scale of the forcing. The magnetic energy exceeds the kinetic energy over dissipative time scales, similar to that seen earlier in Cartesian simulations in periodic boxes. As we increase the size of the domain in the azimuthal direction we find that the nonlinearly saturated magnetic field organizes itself in long-lived cellular structures with aspect ratios close to unity. These structures tile the domain along the azimuthal direction, thus resulting in very small longitudinally averaged magnetic fields for large domain sizes. The scales of these structures are determined by the smallest scales of the domain, which in our simulations is usually the radial scale. We also find that increasing the meridional extent of the domains produces little qualitative change, except a marginal increase in the large-scale field. We obtain qualitatively similar results in Cartesian domains with similar aspect ratios.

The conditions for the existence and stability of cosmological power-law scaling solutions are established when the Einstein-Hilbert action is modified by the inclusion of a function of the Gauss-Bonnet curvature invariant. The general form of the action that leads to such solutions is determined for the case where the universe is sourced by a barotropic perfect fluid. It is shown by employing an equivalence between the Gauss-Bonnet action and a scalar-tensor theory of gravity that the cosmological field equations can be written as a plane autonomous system. It is found that stable scaling solutions exist when the parameters of the model take appropriate values.

We study the field profile of a scalar field $\phi$ that couples to a matter fluid (dubbed a chameleon field) in the relativistic gravitational background of a spherically symmetric spacetime. Employing a linear expansion in terms of the gravitational potential $\Phi_c$ at the surface of a compact object with a constant density, we derive the thin-shell field profile both inside and outside the object, as well as the resulting effective coupling with matter, analytically. We also carry out numerical simulations for the class of inverse power-law potentials $V(\phi)=M^{4+n} \phi^{-n}$ by employing the information provided by our analytical solutions to set the boundary conditions around the centre of the object and show that thin-shell solutions in fact exist if the gravitational potential $\Phi_c$ is smaller than 0.3, which marginally covers the case of neutron stars. Thus the chameleon mechanism is present in the relativistic gravitational backgrounds, capable of reducing the effective coupling. Since thin-shell solutions are sensitive to the choice of boundary conditions, our analytic field profile is very helpful to provide appropriate boundary conditions for $\Phi_c \lesssim O(0.1)$.

We study the dependence of turbulent transport coefficients, such as the components of the $\alpha$ tensor ($\alpha_{ij}$) and the turbulent magnetic diffusivity tensor ($\eta_{ij}$), on shear and magnetic Reynolds number in the presence of helical forcing. We use three-dimensional direct numerical simulations with periodic boundary conditions and measure the turbulent transport coefficients using the kinematic test field method. In all cases the magnetic Prandtl number is taken as unity. We find that with increasing shear the diagonal components of $\alpha_{ij}$ quench, whereas those of $\eta_{ij}$ increase. The antisymmetric parts of both tensors increase with increasing shear. We also propose a simple expression for the turbulent pumping velocity (or $\gamma$ effect). This pumping velocity is proportional to the kinetic helicity of the turbulence and the vorticity of the mean flow. For negative helicity, i.e. for a positive trace of $\alpha_{ij}$, it points in the direction of the mean vorticity, i.e. perpendicular to the plane of the shear flow. Our simulations support this expression for low shear and magnetic Reynolds number. The transport coefficients depend on the wavenumber of the mean flow in a Lorentzian fashion, just as for non-shearing turbulence.

While the topology of the Universe is at present not specified by any known fundamental theory, it may in principle be determined through observations. In particular, a non-trivial topology will generate pairs of matching circles of temperature fluctuations in maps of the cosmic microwave background, the so-called circles-in-the-sky. A general search for such pairs of circles would be extremely costly and would therefore need to be confined to restricted parameter ranges. To draw quantitative conclusions from the negative results of such partial searches for the existence of circles we need a concrete theoretical framework. Here we provide such a framework by obtaining constraints on the angular parameters of these circles as a function of cosmological density parameters and the observer's position. As an example of the application of our results, we consider the recent search restricted to pairs of nearly back-to-back circles with negative results. We show that assuming the Universe to be very nearly flat, with its total matter-energy density satisfying the bounds $ 0 <|\Omega_0 - 1| \lesssim 10^{-5}$, compatible with the predictions of typical inflationary models, this search, if confirmed, could in principle be sufficient to exclude a detectable non-trivial cosmic topology for most observers. We further relate explicitly the fraction of observers for which this result holds to the cosmological density parameters.

Recollapsing homogeneous and isotropic models present one of the key ingredients for cyclic scenarios. This is considered here within a quantum cosmological framework in presence of a free scalar field with, in turn, a negative cosmological constant and spatial curvature. Effective equations shed light on the quantum dynamics around a recollapsing phase and the evolution of state parameters such as fluctuations and correlations through such a turn around. In the models considered here, the squeezing of an initial state is found to be strictly monotonic in time during the expansion, turn around and contraction phases. The presence of such monotonicity is of potential importance in relation to a long standing intensive debate concerning the (a)symmetry between the expanding and contracting phases in a recollapsing universe. Furthermore, together with recent analogous results concerning a bounce one can extend this monotonicity throughout an entire cycle. This provides a strong motivation for employing the degree of squeezing as an alternative measure of (quantum) entropy. It may also serve as a new concept of emergent time described by a variable without classical analog. The evolution of the squeezing in emergent oscillating scenarios can in principle provide constraints on the viability of such models.

We construct a family of viable scalar-tensor models of dark energy (DE) which possess a phase of late-time acceleration preceded by a standard matter era, while at the same time satisfying the local gravity constraints (LGC). The coupling Q between the scalar field and the non-relativistic matter in the Einstein frame is assumed to be constant in our scenario, which is a generalization of f(R) gravity theories corresponding to the coupling Q=-1/sqrt{6}. We find that these models can be made compatible with local gravity constraints even when |Q| is of the order of unity through a chameleon mechanism, if the scalar-field potential is chosen to have a sufficiently large mass in the high-curvature regions. We show that these models generally lead to the divergence of the equation of state of DE, which occur at smaller redshifts as the deviation from the LambdaCDM model become more significant. We also study the evolution of matter density perturbations and employ them to place bounds on the coupling |Q| as well as model parameters of the field potential from observations of the matter power spectrum and the CMB anisotropies. We find that, as long as |Q| is smaller than the order of unity, there exist allowed parameter regions that are consistent with both observational and local gravity constraints.

We make a detailed study of matter density perturbations in both metric and Palatini formalisms in theories whose Lagrangian density is a general function, f(R), of the Ricci scalar. We derive these equations in a number of gauges. We show that for viable models that satisfy cosmological and local gravity constraints (LGC), matter perturbation equations derived under a sub-horizon approximation are valid even for super-Hubble scales provided the oscillating mode (scalaron) does not dominate over the matter-induced mode. Such approximate equations are especially reliable in the Palatini formalism because of the absence of scalarons. Using these equations we make a comparative study of the behaviour of density perturbations as well as gravitational potentials for a number of classes of theories. In the metric formalism the parameter m=Rf_{,RR}/f_{,R} characterising the deviation from the Lambda CDM model is constrained to be very small during the matter era in order to ensure compatibility with LGC, but the models in which m grows to the order of 10^{-1} around the present epoch can be allowed. These models also suffer from an additional fine tuning due to the presence of scalaron modes which are absent in the Palatini case. In Palatini formalism LGC and background cosmological constraints provide only weak bounds on |m| by constraining it to be smaller than ~ 0.1. This is in contrast to matter density perturbations which, on galactic scales, place far more stringent constraints on the present deviation parameter m of the order of |m| < 10^{-5} - 10^{-4}. This is due to the peculiar evolution of matter perturbations in the Palatini case which exhibits a rapid growth or a damped oscillation depending on the sign of m.

The motivation for considering distributed large scale dynamos in the solar context is reviewed in connection with the magnetic helicity constraint. Preliminary accounts of 3-dimensional direct numerical simulations (in spherical shell segments) and simulations of 2-dimensional mean field models (in spherical shells) are presented. Interesting similarities as well as some differences are noted.

Two approaches to the study of cosmological density perturbations in modified theories of Palatini gravity have recently been discussed. These utilise, respectively, a generalisation of Birkhoff's theorem and a direct linearization of the gravitational field equations. In this paper these approaches are compared and contrasted. The general form of the gravitational lagrangian for which the two frameworks yield identical results in the long-wavelength limit is derived. This class of models includes the case where the lagrangian is a power-law of the Ricci curvature scalar. The evolution of density perturbations in theories of the type $f(R)=R-c /R^ b$ is investigated numerically. It is found that the results obtained by the two methods are in good agreement on sufficiently large scales when the values of the parameters (b,c) are consistent with current observational constraints. However, this agreement becomes progressively poorer for models that differ significantly from the standard concordance model and as smaller scales are considered.

Recent analyses of the WMAP data seem to indicate the possible presence of large-angle anisotropy in the Universe. If confirmed, these can have important consequences for our understanding of the Universe. A number of attempts have recently been made to establish the reality and nature of such anisotropies in the CMB data. Among these is a directional indicator recently proposed by the authors. A distinctive feature of this indicator is that it can be used to generate a sky map of the large-scale anisotropies of the CMB maps. Applying this indicator to full-sky temperature maps we found a statistically significant preferred direction. The full-sky maps used in these analyses are known to have residual foreground contamination as well as complicated noise properties. Thus, here we performed the same analysis for a map where regions with high foreground contamination were removed. We find that the main feature of the full-sky analysis, namely the presence of a significant axis of asymmetry, is robust with respect to this masking procedure. Other subtler anomalies of the full-sky are on the other hand no longer present.

A number of scalar field models proposed as alternatives to the standard inflationary scenario involve contracting phases which precede the universe's present phase of expansion. An important question concerning such models is whether there are effects which could potentially distinguish them from purely expanding cosmologies. Vector perturbations have recently been considered in this context. At first order such perturbations are not supported by a scalar field. In this paper, therefore, we consider second order vector perturbations. We show that such perturbations are generated by first order scalar mode-mode couplings, and give an explicit expression for them. We compare the magnitude of vector perturbations produced in collapsing models with the corresponding amplitudes produced during inflation, using a number of suitable power-law solutions to model the inflationary and collapsing scenarios. We conclude that the ratios of the magnitudes of these perturbations depend on the details of the collapsing scenario as well as on how the hot big bang is recovered, but for certain cases could be large, growing with the duration of the collapse.

We make a systematic study of the cosmological dynamics for a number of f(R) gravity theories in Palatini formalism. We find a number of interesting results: (i) models based on theories of the type (a) f(R)=R-beta R^n and (b) f(R)=R+alpha ln R -beta, unlike the metric formalism, are capable of producing the sequence of radiation-dominated, matter-dominated and de-Sitter periods, and (ii) models based on theories of the type (c) f(R)=R+alpha R^m -beta R^n can produce early as well as late accelerating phases. However for the classes of models considered here, we have been unable to find the sequence of all four dynamical epochs required to account for the complete cosmological dynamics, even though three out of four phases are possible. We also place observational constraints on these models using the recently released supernovae data (SNLS) as well as the baryon acoustic oscillation peak and the CMB shift parameter. The best-fit values are found to be n=0.027, alpha=4.63 for the models based on (a) and alpha=0.11, beta=4.62 for the models based on (b), neither of which are significantly preferred over the LCDM model. The models based on (c) are also consistent with the data with suitable choices of their parameters.

Analyses of recent cosmic microwave background (CMB) observations have provided increasing indications for the existence of large scale anisotropy in the universe. Given the far reaching consequences of such an anisotropy for our understanding of the universe, it is important to employ alternative indicators in order to determine whether the reported anisotropy is cosmological in origin, and if so extract further information that may be helpful for identifying its causes. Here we propose a new directional indicator, based on angular-separation histograms of pairs of pixels with similar temperatures in the CMB map, which provides a measure of departure from statistical isotropy. The main advantage of this indicator is that it can be used to generate a sky map of large-scale anisotropies in the CMB temperature map, thus allowing a possible additional window into their causes. Using this indicator, we find a statistically significant (at 95% CL) preferred direction in the CMB data and discuss how it compares with other such axes recently reported. We also show that our findings are robust with respect to both the details of the method used, and the choice of the WMAP CMB maps employed, including the three-year CMB data released recently.

We put forward a new model-independent reconstruction scheme for dark energy which utilises the expected geometrical features of the luminosity-distance relation. The important advantage of this scheme is that it does not assume explicit ansatzes for cosmological parameters but only some very general cosmological properties via the geometrical features of the reconstructed luminosity-distance relation. Using the recently released supernovae data by the Supernova Legacy Survey together with a phase space representation, we show that the reconstructed luminosity-distance curves best fitting the data correspond to a slightly varying dark energy density with the Universe expanding slightly slower than the Lambda CDM model. However, the Lambda CDM model fits the data at 1 sigma significance level and the fact that our best fitting luminosity-distance curve is lower than that of the corresponding Lambda CDM model could be due to systematics. The transition from an accelerating to a decelerating expansion occurs at a redshift larger than z=0.35. Interpreting the dark energy as a minimally coupled scalar field we also reconstruct the scalar field and its potential. We constrain $\Omega_{m_0}$ using the baryon acoustic oscillation peak in the SDSS luminous red galaxy sample and find that the best fit is obtained with $\Omega_{m_0}=0.27$, in agreement with the CMB data.

Recent high precision data by WMAP and SDSS have provided strong evidence to suggest that the universe is nearly flat. They are also making it possible to probe the topology of the universe. Motivated by these results, we have recently studied the consequences of taking the inflationary limit, i.e. $|\Omega_0 - 1| << 1$. We have shown that in this limit a generic detectable spherical or hyperbolic topology is locally indistinguishable from either $R^2 x S^1$ or $R x T^2$, irrespective of its global shape. Here we briefly present these results and further discuss their observational implications.

Closed, singularity-free, inflationary cosmological models have recently been studied in the context of general relativity. Despite their appeal, these so called emergent models suffer from a number of limitations. These include the fact that they rely on an initial Einstein static state to describe the past eternal phase of the universe. Given the instability of such a state within the context of general relativity, this amounts to a very severe fine tuning. Also in order to be able to study the dynamics of the universe within the context of general relativity, they set the initial conditions for the universe in the classical phase. Here we study the existence and stability of such models in the context of Loop Quantum Cosmology and show that both these limitations can be partially remedied, once semi-classical effects are taken into account. An important consequence of these effects is to give rise to a static solution (not present in GR), which dynamically is a centre equilibrium point and located in the more natural semi-classical regime. This allows the construction of emergent models in which the universe oscillates indefinitely about such an initial static state. We construct an explicit emergent model of this type, in which a non-singular past eternal oscillating universe enters a phase where the symmetry of the oscillations is broken, leading to an emergent inflationary epoch, while satisfying all observational and semi-classical constraints. We also discuss emergent models in which the universe possesses both early- and late-time accelerating phases.

Increasing evidence is becoming available about not only the surface differential rotation of rapidly rotating cool stars but, in a small number of cases, also about temporal variations, which possibly are analogous to the solar torsional oscillations. Given the present difficulties in resolving the precise nature of such variations, due to both the short length and poor resolution of the available data, theoretical input is vital to help assess the modes of behaviour that might be expected, and will facilitate interpretation of the observations. Here we take a first step in this direction by studying the variations in the convection zones of such stars, using a two dimensional axisymmetric mean field dynamo model operating in a spherical shell in which the only nonlinearity is the action of the azimuthal component of the Lorentz force of the dynamo generated magnetic field on the stellar angular velocity. We consider three families of models with different depths of dynamo-active regions. For moderately supercritical dynamo numbers we find torsional oscillations that penetrate all the way down to the bottom of the convection zones, similar to the case of the Sun. For larger dynamo numbers we find fragmentation in some cases and sometimes there are other dynamical modes of behaviour, including quasi-periodicity and chaos. We find that the largest deviations in the angular velocity distribution caused by the Lorentz force are of the order of few percent, implying that the original assumed `background' rotation field is not strongly distorted.

Positively-curved, oscillatory universes are studied within the context of Loop Quantum Cosmology subject to a consistent semi-classical treatment. The semi-classical effects are reformulated in terms of an effective phantom fluid with a variable equation of state. In cosmologies sourced by a massless scalar field, these effects lead to a universe that undergoes ever-repeating cycles of expansion and contraction. The presence of a self-interaction potential for the field breaks the symmetry of the cycles and can enable the oscillations to establish the initial conditions for successful slow-roll inflation, even when the field is initially at the minimum of its potential with a small kinetic energy. The displacement of the field from its minimum is enhanced for lower and more natural values of the parameter that sets the effective quantum gravity scale. For sufficiently small values of this parameter, the universe can enter a stage of eternal self-reproduction.

In loop quantum gravity, modifications to the geometrical density cause a self-interacting scalar field to accelerate away from a minimum of its potential. In principle, this mechanism can generate the conditions that subsequently lead to slow-roll inflation. The consequences for this mechanism of various quantization ambiguities arising within loop quantum cosmology are considered. For the case of a quadratic potential, it is found that some quantization procedures are more likely to generate a phase of slow--roll inflation. In general, however, loop quantum cosmology is robust to ambiguities in the quantization and extends the range of initial conditions for inflation.

The observed values of density parameters inevitably involve uncertainties. We study the conditions for detectability and undetectability of cosmic topology in presence of such uncertainties. We present closed analytical forms of (un)detectability conditions for infinite redshift, which are important because: (i) they allow the examine of the detectability of cosmic topology not only for individual manifolds (topologies), but also in whole classes of manifolds; (ii) they are, to a very good approximation, (un)detectibility conditions for z=1100.

Analyses of recent helioseismic data have produced ample evidence for substantial dynamical variation of the differential rotation within the solar convection zone. Given the inevitable difficulties in resolving the precise nature of variations at deeper layers, much effort has recently gone into determining theoretically the expected modes of behaviour, using nonlinear dynamo models. Two important limitations of these models are that they have so far included only one form of nonlinearity, and as yet they have not taken into account the density stratification in the solar convection zone. Here we address both of these issues by studying the effects of including density stratification, as well as including an alpha--quenching nonlinearity in addition to the previously studied effects of the Lorentz force on the differential rotation. We find that observationally important features found in the earlier uniform density models remain qualitatively unchanged, although there are quantitative differences. This is important as it provides more realistic theoretical predictions to be compared with and guide observations, especially in the deeper regions where the uncertainties in the inversions are larger. However the presence of an effective alpha-quenching nonlinearity significantly reduces the amplitudes of the oscillations.

Recent observational results suggest that our universe is nearly flat and well modelled within a $\Lambda$CDM framework. The observed values of $\Omega_{m}$ and $\Omega_{\Lambda}$ inevitably involve uncertainties. Motivated by this, we make a systematic study of the necessary and sufficient conditions for undetectability as well as detectability (in principle) of cosmic topology (using pattern repetition) in presence of such uncertainties. We do this by developing two complementary methods to determine detectability for nearly flat universes. Using the first method we derive analytical conditions for undetectability for infinite redshift, the accuracy of which is then confirmed by the second method. Estimates based on WMAP data together with other measurements of the density parameters are used to illustrate both methods, which are shown to provide very similar results for high redshifts.

Recent studies of the detectability of cosmic topology of nearly flat universes have often concentrated on the range of values of $\Omega_{0}$ given by current observations. Here we study the consequences of taking the bounds on $\Omega_{0}$ given by inflationary models, i.e. $|\Omega_0 - 1| \ll 1$. We show that in this limit, a generic detectable non-flat manifold is locally indistinguishable from either a cylindrical ($R^2 \times S^1$) or toroidal ($R \times T^2$) manifold, irrespective of its global shape, with the former being more likely. Importantly this is compatible with some recent indications based on the analysis of high resolution CMB data. It also implies that in this limit an observer would not be able to distinguish topologically whether the universe is spherical, hyperbolic or flat. By severely restricting the expected topological signatures of detectable isometries, our results provide an effective theoretical framework for interpreting cosmological observations, and can be used to confine any parameter space which realistic search strategies, such as the `circles in the sky' method, need to concentrate on. This is particularly important in the inflationary limit, where the precise nature of cosmic topology becomes undecidable.

Some of the consequences for inflationary cosmology of a scale dependence (running) in the tilt of the scalar perturbation spectrum are considered. In the limit where the running is itself approximately scale-invariant, a relationship is found between the scalar and tensor perturbation amplitudes, the scalar spectral index and its running. This relationship is independent of the functional form of the inflaton potential. More general settings, including that of braneworld cosmological models, are also considered. It is found that for the Randall-Sundrum single braneworld scenario, the corresponding relation between the observables takes precisely the same form as that arising in the standard cosmology. Some implications of the observations failing to satisfy such a relationship are discussed.

We reexamine the possibility of the detection of the cosmic topology in nearly flat hyperbolic Friedmann-Lemaitre-Robertson-Walker (FLRW) universes by using patterns repetition. We update and extend our recent results in two important ways: by employing recent observational constraints on the cosmological density parameters as well as the recent mathematical results concerning small hyperbolic 3-manifolds. This produces new bounds with consequences for the detectability of the cosmic topology. In addition to obtaining new bounds, we also give a concrete example of the sensitive dependence of detectability of cosmic topology on the uncertainties in the observational values of the density parameters.

We summarize recent results concerning the evolution of second order perturbations in flat dust irrotational FLRW models with $\Lambda\ne 0$. We show that asymptotically these perturbations tend to constants in time, in agreement with the cosmic no-hair conjecture. We solve numerically the second order scalar perturbation equation, and very briefly discuss its all time behaviour and some possible implications for the structure formation.

Recent analyses of the helioseismic data have produced evidence for a variety of interesting dynamical behaviour associated with torsional oscillations. What is not so far clear is whether these oscillations extend all the way to the bottom of the convection zone and, if so, whether the oscillatory behaviour at the top and the bottom of the convection zone is different. Attempts have been made to understand such modes of behaviour within the framework of nonlinear dynamo models which include the nonlinear action of the Lorentz force of the dynamo generated magnetic field on the solar angular velocity. One aspect of these models that remains uncertain is the nature of the boundary conditions on the magnetic field. Here by employing a range of physically plausible boundary conditions, we show that for near-critical and moderately supercritical dynamo regimes, the oscillations extend all the way down to the bottom of the convection zone. Thus, such penetration is an extremely robust feature of the models considered. We also find parameter ranges for which the supercritical models show spatiotemporal fragmentation for a range of choices of boundary conditions. Given their observational importance, we also make a comparative study of the amplitude of torsional oscillations as a function of the boundary conditions.

We report our first results on comparing the variations of the solar internal rotation with solar activity, as predicted by non-linear solar dynamo modelling, with helioseismic measurements using the SOHO MDI data.

Using recent observational constraints on cosmological density parameters, together with recent mathematical results concerning small volume hyperbolic manifolds, we argue that, by employing pattern repetitions, the topology of nearly flat small hyperbolic universes can be observationally undetectable. This is important in view of the facts that quantum cosmology may favour hyperbolic universes with small volumes, and from the expectation coming from inflationary scenarios, that $\Omega_0$ is likely to be very close to one.

Recent observations suggest that the ratio of the total density to the critical density of the universe, $\Omega_0$, is likely to be very close to one, with a significant proportion of this energy being in the form of a dark component with negative pressure. Motivated by this result, we study the question of observational detection of possible non-trivial topologies in universes with $\Omega_0 \sim 1$, which include a cosmological constant. Using a number of indicators we find that as $\Omega_0 \to 1$, increasing families of possible manifolds (topologies) become either undetectable or can be excluded observationally. Furthermore, given a non-zero lower bound on $|\Omega_0 - 1|$, we can rule out families of topologies (manifolds) as possible candidates for the shape of our universe. We demonstrate these findings concretely by considering families of topologies and employing bounds on cosmological parameters from recent observations. We find that given the present bounds on cosmological parameters, there are families of both hyperbolic and spherical manifolds that remain undetectable and families that can be excluded as the shape of our universe. These results are of importance in future search strategies for the detection of the shape of our universe, given that there are an infinite number of theoretically possible topologies and that the future observations are expected to put a non-zero lower bound on $|\Omega_0 - 1|$ which is more accurate and closer to zero.

We study the robustness of the quintessence tracking scenario in the context of more general cosmological models that derive from high-energy physics. We consider the effects of inclusion of multiple scalar fields, corrections to the Hubble expansion law (such as those that arise in brane cosmological models), and potentials that decay with expansion of the Universe. We find that in a successful tracking quintessence model the average equation of state must remain nearly constant. Overall, the conditions for successful tracking become more complex in these more general settings. Tracking can become more fragile in presence of multiple scalar fields, and more stable when temperature dependent potentials are present. Interestingly though, most of the cases where tracking is disrupted are those in which the cosmological model is itself non-viable due to other constraints. In this sense tracking remains robust in models that are cosmologically viable.

We study the asymptotic stability of de Sitter spacetime with respect to non-linear perturbations, by considering second order perturbations of a flat Robertson-Walker universe with dust and a positive cosmological constant. Using the synchronous comoving gauge we find that, as in the case of linear perturbations, the non-linear perturbations also tend to constants, asymptotically in time. Analysing curvature and other spacetime invariants we show, however, that these quantities asymptotically tend to their de Sitter values, thus demonstrating that the geometry is indeed locally asymptotically de Sitter, despite the fact that matter inhomogeneities tend to constants in time. Our results support the inflationary picture of frozen amplitude matter perturbations that are stretched outside the horizon, and demonstrate the validity of the cosmic no-hair conjecture in the nonlinear inhomogeneous settings considered here.

Analyses of recent helioseismic data indicate that the dynamical regimes at the base of the convection zone can be different from those observed at the top, having either significantly shorter periods or non--periodic behaviour. Recently spatiotemporal fragmentation/bifurcation has been proposed as a dynamical mechanism to account for the multi-mode behaviour that is possibly observed in the solar convection zone, without requiring separate physical mechanisms with different time scales at different depths. Here we study the robustness of this mechanism with respect to changes to the zero order rotation profile, motivated by the uncertainties of and differences between the various reductions of the helioseimic data. We find that spatiotemporal fragmentation is a common feature of the reductions investigated

Observations of the Sun and solar-type stars continue to reveal phenomena whose understanding is very likely to require a nonlinear framework. Here we shall concentrate on two such phenomena, namely the grand minima type behaviour observed in the Sun and solar-type stars and the recent dynamical variations of the differential rotation in the solar convection zone, deduced from the helioseismic observations, and discuss how their explanations have recently motivated the development/employment of novel ideas from nonlinear dynamics.

We find concrete evidence for a recently discovered form of intermittency, referred to as in--out intermittency, in both PDE and ODE models of mean field dynamos. This type of intermittency (introduced in Ashwin et al 1999) occurs in systems with invariant submanifolds and, as opposed to on--off intermittency which can also occur in skew product systems, it requires an absence of skew product structure. By this we mean that the dynamics on the attractor intermittent to the invariant manifold cannot be expressed simply as the dynamics on the invariant subspace forcing the transverse dynamics; the transverse dynamics will alter that tangential to the invariant subspace when one is far enough away from the invariant manifold. Since general systems with invariant submanifolds are not likely to have skew product structure, this type of behaviour may be of physical relevance in a variety of dynamical settings. The models employed here to demonstrate in--out intermittency are axisymmetric mean--field dynamo models which are often used to study the observed large scale magnetic variability in the Sun and solar-type stars. The occurrence of this type of intermittency in such models may be of interest in understanding some aspects of such variabilities.

Recent analyses of helioseismological observations seem to suggest the presence of two new phenomena connected with the dynamics of the solar convective zone. Firstly, there are present torsional oscillations with periods of about 11 years, which penetrate significantly into the solar convection zone and secondly, oscillatory regimes exist near the base of the convection which are markedly different from those observed near the top, having either significantly reduced periods or being non-periodic. Recently spatiotemporal fragmentation/bifurcation has been proposed as a possible dynamical mechanism to account for such observed multi-mode behaviours in different parts of the solar convection zone. Evidence for this scenario was produced in the context of an axisymmetric mean field dynamo model operating in a spherical shell, with a semi-open outer boundary condition and a zero order angular velocity obtained by the inversion of the MDI data, in which the only nonlinearity was the action of the Lorentz force of the dynamo generated magnetic field on the solar angular velocity. Here we make a detailed study of the robustness of this model with respect to plausible changes to its main ingredients, including changes to the alpha and eta profiles as well as the inclusion of a nonlinear alpha quenching. We find that spatiotemporal fragmentation is present in this model for different choices of the rotation data and as the details of the model are varied. Taken together, these results give strong support to the idea that spatiotemporal fragmentation is likely to occur in general dynamo settings.

Recent analyses of the helioseismic observations indicate that the previously observed surface torsional oscillations with periods of about 11 years extend significantly downwards into the solar convective zone. Furthermore, there are indications that the dynamical regimes at the base of the convection zone are different from those observed at the top, having either significantly shorter periods or non-periodic behaviour. We propose that this behaviour can be explained by the occurrence of {\it spatiotemporal fragmentation}, a crucial feature of which is that such behaviour can be explained solely through nonlinear spatiotemporal dynamics, without requiring separate mechanisms with different time scales at different depths. We find evidence for this mechanism in the context of a two dimensional axisymmetric mean field dynamo model operating in a spherical shell, with a semi-open outer boundary condition, in which the only nonlinearity is the action of the azimuthal component of the Lorentz force of the dynamo generated magnetic field on the solar angular velocity.

Recent analysis of the helioseismic observations indicate that the previously observed surface torsional oscillations extend significantly downwards into the solar convection zone. In an attempt to understand these oscillations, we study the nonlinear coupling between the magnetic field and the solar differential rotation in the context of a mean field dynamo model, in which the nonlinearity is due to the action of the azimuthal component of the Lorentz force of the dynamo generated magnetic field on the solar angular velocity. The underlying zero order angular velocity is chosen to be consistent with the most recent helioseismic data. The model produces butterfly diagrams which are in qualitative agreement with the observations. It displays torsional oscillations that penetrate into the convection zone, and which with time migrate towards the equator. The period of these oscillations is found to be half that of the period of the global magnetic fields. This is compatible with the observed period of the surface torsional oscillations. Inside the convection zone, this is a testable prediction that is not ruled out by the observations so far available.

We consider a recent generalisation by Bousso of an earlier holography proposal by Fischler and Susskind. We demonstrate that in general inhomogeneous universes such a proposal would involve extremely complicated - possibly fractal - light sheets. Furthermore, in general such a light sheet cannot be known a priori on the basis of theory and moreover, the evolution of the universe makes it clear that in general such bounds cannot remain invariant under time reversal and will change with epoch. We propose a modified version of this proposal in which the light sheets end on the boundary of the past, and hence avoid contact with the caustics. In this way the resulting light sheets and projections can be made much simpler. We discuss the question of operational definability of these sheets within the context of both proposals and conclude that in both cases the theoretical existence of such sheets must be clearly distinguished from their complexity and the difficulty of their construction in practice. This puts into perspective the likely practical difficulties one would face in applying the holographic principle to the real cosmos. These issues may also be of relevance in debates regarding the applications of the holographic principle to other settings such as string theory.

We find concrete evidence for the presence of crisis-induced and Pomeau-Manneville Type-I intermittencies in an axisymmetric PDE mean-field dynamo model. These findings are of potential importance for two different reasons. Firstly, as far as we are aware, this is the first time detailed evidence has been produced for the occurrence of these types of intermittency for such deterministic PDE models. And secondly, despite the rather idealised nature of these models, the concrete evidence for the occurrence of more than one type of intermittency in such models makes it in principle possible that different types of intermittency may occur in different solar-type stars or even in the same star over different epochs. In this way a multiple intermittency framework may turn out to be of importance in understanding the mechanisms responsible for grand-minima type behaviour in the Sun and solar-type stars and in particular in the interpretation of the corresponding observational and proxy evidence.

We briefly discuss the status of the intermittency hypothesis, according to which the grand minima type variability in solar-type stars may be understood in terms of dynamical intermittency. We review concrete examples which establish this hypothesis in the mean-field setting. We discuss some difficulties and open problems regarding the establishment of this hypothesis in more realistic settings as well as its operationally decidability.

Global symmetries of the string effective action are employed to generate tilted, homogeneous Bianchi type VI_h string cosmologies from a previously known stiff perfect fluid solution to Einstein gravity. The dilaton field is not constant on the surfaces of homogeneity. The future asymptotic state of the models is interpreted as a plane wave and is itself an exact solution to the string equations of motion to all orders in the inverse string tension. An inhomogeneous generalization of the Bianchi type III model is also found.

We study the changes in the dynamical behaviour of axisymmetric spherical mean-field dynamo models produced by changes in their geometry and topology, by considering a two parameter family of models, ranging from a full sphere to spherical shell, torus and disc-like configurations, within a unified framework. We find that the two parameter space of the family of models considered here separates into at least three different regions with distinct characteristics for the onset of dynamo action. In two of these regions, the most easily excited fields are oscillatory, in one case with dipolar symmetry, and in the other with quadrupolar, whereas in the third region the most easily excited field is steady and quadrupolar. In the nonlinear regime, we find that topological changes can alter significantly the dynamical behaviour, whilst modest changes in geometry can produce qualitative changes, particularly for thin disc-like configurations. This is of potential importance, since the exact shapes of astrophysical bodies, especially accretion discs and galaxies, are usually not precisely known.

Calculations for mean field dynamo models (in both full spheres and spherical shells), with both algebraic and dynamic $\alpha$--quenchings, show qualitative as well as quantitative differences and similarities in the dynamical behaviour of these models. We summarise and enhance recent results with extra examples. Overall, the effect of using a dynamic $\alpha$ appears to be complicated and is affected by the region of parameter space examined.

An important question regarding the study of mean field dynamo models is how to make precise the nature of their underlying dynamics. This is difficult both because relatively little is known about the dynamical behaviour of infinite dimensional systems and also due to the numerical cost of studying the related partial differential equations. As a first step towards their understanding, it is useful to consider the corresponding truncated models. Here we summarise some recent results of the study of a class of truncated axisymmetric mean field dynamo models. We find conclusive evidence in these models for various types of intermittency as well as multiple attractors and final state sensitivity. We also find that the understanding of the underlying dynamics of such dynamo models requires the study of a new class of dynamical systems, referred to as the {\it non-normal} systems. Current work demonstrates that these types of systems are capable of a novel type of intermittency and also of relevance for the understanding of the full axisymmetric PDE dynamo models.

We introduce the notion of topological fragility and briefly discuss some examples from the literature. An important example of this type of fragility is the way globally anisotropic Bianchi V generalisations of the FLRW $k=-1$ model result in a radical restriction on the allowed topology of spatial sections, thereby excluding compact cosmological models with negatively curved three-sections with anisotropy. An outcome of this is to exclude chaotic mixing in such models, which may be relevant, given the many recent attempts at employing compact FLRW $k=-1$ models to produce chaotic mixing in the cosmic microwave background radiation, if the Universe turns out to be globally anisotropic.

We study axisymmetric mean field spherical and spherical shell dynamo models, with both dynamic and algebraic $\alpha$--quenchings. Our results show that there are qualitative as well as quantitative differences and similarities between these models. Regarding similarities, both groups of models exhibit symmetric, antisymmetric and mixed modes of behaviour. As regards differences, the important feature in the full sphere models is the occurrence of chaotic behaviour in the algebraic $\alpha$--quenching models. For the spherical shell models with dynamic $\alpha$ the main features include the possibility of multi-attractor regimes with final state sensitivity with respect to small changes in the magnitude of $\alpha$ and the initial parity. We find the effect of introducing a dynamic $\alpha$ is likely to be complicated and depend on the region of the parameter space considered, rather than a uniform change towards simplicity or complexity.

We investigate the detailed dynamics of a truncated $\alpha\omega$ dynamo model with a dynamic $\alpha$ effect. We find the presence of multiple attractors, including two chaotic attractors with a fractal basin boundary which merge to form a single attractor as the control parameter is increased. By considering phase portraits and the scaling of averaged times of transitions between the two attractors, we demonstrate that this merging is accompanied by a crisis-induced intermittency. We also find a range of parameter values over which the system has a fractal parameter dependence for fixed initial conditions. This is the first time this type of intermittency has been observed in a dynamo model and it could be of potential importance in accounting for some forms of intermittency in the solar and stellar output.

We investigate the behaviour of $\alpha\Omega$ dynamos with a dynamic $\alpha$, whose evolution is governed by the imbalance between a driving and a damping term. We focus on truncated versions of such dynamo models which are often studied in connection with solar and stellar variability. Given the approximate nature of such models, it is important to study how robust they are with respect to reasonable changes in the formulation of the driving and damping terms. For each case, we also study the effects of changes of the dynamo number and its sign, the truncation order and initial conditions. Our results show that changes in the formulation of the driving term have important consequences for the dynamical behaviour of such systems, with the detailed nature of these effects depending crucially on the form of the driving term assumed, the value and the sign of the dynamo number and the initial conditions. On the other hand, the change in the damping term considered here seems to produce little qualitative effect.

In a recent work (Covas et al. 1997), the behaviour and the robustness of truncated $\alpha\Omega$ dynamos with a dynamic $\alpha$ were studied with respect to a number of changes in the driving term of the dynamic $\alpha$ equation, which was considered previously by Schmalz and Stix (1991) to be of the form $\sim A_{\phi}B_{\phi}$. Here we review and extend our previous work and consider the effect of adding a quadratic quenching term of the form $\alpha|B|^2$. We find that, as before, such a change can have significant effects on the dynamics of the related truncated systems. We also find intervals of (negative) dynamo numbers, in the system considered by Schmalz and Stix (1991), for which there is sensitivity with respect to small changes in the dynamo number and the initial conditions, similar to what was found in our previous work. This latter behaviour may be of importance in producing the intermittent type of behaviour observed in the sun.

We consider the domain of applicability of general relativity (GR), as a classical theory of gravity, by considering its applications to a variety of settings of physical interest as well as its relationship with real observations. We argue that, as it stands, GR is deficient whether it is treated as a microscopic or a macroscopic theory of gravity. We briefly discuss some recent attempts at removing this shortcoming through the construction of a macroscopic theory of gravity. We point out that such macroscopic extensions of GR are likely to be non-unique and involve non-Riemannian geometrical frameworks.