Researcher profile

Reza Tavakol

Reza Tavakol contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
18works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

18 published item(s)

preprint2013arXiv

An exact quantification of backreaction in relativistic cosmology

An important open question in cosmology is the degree to which the Friedmann-Lemaitre-Robertson-Walker (FLRW) solutions of Einstein&#39;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&#39;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.

preprint2013arXiv

What Planck does not tell us about inflation

Planck data has not found the &#39;smoking gun&#39; 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.

preprint2012arXiv

Evolution of fNL to the adiabatic limit

We study inflationary perturbations in multiple-field models, for which zeta typically evolves until all isocurvature modes decay--the &#34;adiabatic limit&#34;. 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.

preprint2012arXiv

Large trispectrum in two-field slow-roll inflation

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 &#39;heatmap&#39; 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.

preprint2012arXiv

The inflationary bispectrum with curved field-space

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 &#34;delta-N&#34; 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.

preprint2011arXiv

Alleviating alpha quenching by solar wind and meridional flow

We study the ability of magnetic helicity expulsion to alleviate catastrophic $α$-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 $α^2Ω$ mean field models with a dynamical $α$ quenching. We include two mechanisms which have the potential to facilitate helicity expulsion, namely advection by a mean flow (&#34;solar wind&#34;) 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.

preprint2011arXiv

Evolution of non-Gaussianity in multi-scalar field models

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.

preprint2010arXiv

The high-density regime of kinetic-dominated loop quantum cosmology

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.

preprint2009arXiv

Chameleon scalar fields in relativistic gravitational backgrounds

We study the field profile of a scalar field $ϕ$ 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 $Φ_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(ϕ)=M^{4+n} ϕ^{-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 $Φ_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 $Φ_c \lesssim O(0.1)$.

preprint2009arXiv

Cosmological scaling solutions in generalised Gauss-Bonnet gravity theories

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.

preprint2009arXiv

Equatorial magnetic helicity flux in simulations with different gauges

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&#39; 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.

preprint2009arXiv

On vacuum density, the initial singularity and dark energy

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.

preprint2009arXiv

Turbulent dynamos in spherical shell segments of varying geometrical extent

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.

preprint2007arXiv

Cosmological perturbations in Palatini modified gravity

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&#39;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.

preprint2005arXiv

Geodesics, the Equivalence Principle and Singularities in Higher-dimensional General Relativity and Braneworlds

The geodesics of a spacetime seldom coincide with those of an embedded submanifold of codimension one. We investigate this issue for higher-dimensional general relativity-like models, firstly in the simpler case without branes to isolate which features are already present, and then in the more complicated case with branes. The framework in which we consider branes is general enough to include asymmetric braneworlds but not thick branes. We apply our results on geodesics to study both the equivalence principle and cosmological singularities. Among the models we study these considerations favour $Z_2$ symmetric braneworlds with a negative bulk cosmological constant.

preprint2001arXiv

Cosmic no-hair: non-linear asymptotic stability of de Sitter universe

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.

preprint1994arXiv

Quantum Cosmology and Higher-Order Lagrangian Theories

In this paper the quantum cosmological consequences of introducing a term cubic in the Ricci curvature scalar $R$ into the Einstein--Hilbert action are investigated. It is argued that this term represents a more generic perturbation to the action than the quadratic correction usually considered. A qualitative argument suggests that there exists a region of parameter space in which neither the tunneling nor the no-boundary boundary conditions predict an epoch of inflation that can solve the horizon and flatness problems of the big bang model. This is in contrast to the $R^2$--theory.