We consider a new form of theories of gravity in which the action is written in terms of the Ricci scalar and its first and second derivatives. Despite the higher derivative nature of the action, the theory is free from ghost under an appropriate choice of the functional form of the Lagrangian. This model possesses $2+2$ physical degrees of freedom, namely $2$ scalar degrees and $2$ tensor degrees. We exhaust all such theories with the Lagrangian of the form $f(R, (\nabla R)^2, \Box R)$, where $R$ is the Ricci scalar, and then show some examples beyond this ansatz. In course of analysis, we prove the equivalence between these examples and generalized bi-Galileon theories.

The presence of a light vector field coupled to a scalar field during inflation makes a distinct prediction: the observed correlation functions of the cosmic microwave background (CMB) become statistically anisotropic. We study the implications of the current bound on statistical anisotropy derived from the Planck 2013 CMB temperature data for such a model. The previous calculations based on the attractor solution indicate that the magnitude of anisotropy in the power spectrum is proportional to $N^2$, where $N$ is the number of $e$-folds of inflation counted from the end of inflation. In this paper, we show that the attractor solution is not compatible with the current bound, and derive new predictions using another branch of anisotropic inflation. In addition, we improve upon the calculation of the mode function of perturbations by including the leading-order slow-roll corrections. We find that the anisotropy is roughly proportional to $[2(\varepsilon_H+4\eta_H)/3-4(c-1)]^{-2}$, where $\varepsilon_H$ and $\eta_H$ are the usual slow-roll parameters and $c$ is the parameter in the model, regardless of the form of potential of an inflaton field. The bound from Planck implies that breaking of rotational invariance during inflation (characterized by the background homogeneous shear divided by the Hubble rate) is limited to be less than ${\cal O}(10^{-9})$. This bound is many orders of magnitude smaller than the amplitude of breaking of time translation invariance, which is observed to be ${\cal O}(10^{-2})$.

We show that there is a class of spacetime curvature singularities which can be resolved with metric and matter field transformations. As an example, we consider an anisotropic power-law inflation model with gauge and scalar fields in which a space-like curvature singularity exists at the beginning of time. First, we provide a transformation of the metric to the flat geometry. The transformation is regular in the whole region of spacetime except for the singularity. In general, matter fields are still singular after such a metric transformation. However, we explicitly show that there is a case in which the singular behavior of the matter fields can be completely removed by a field re-definition. Since the action is invariant under any metric and matter field transformations, the regularity of the action at the original singularity is a necessary condition for the complete removal of a singularity.

We study cosmological applications of extended vector-tensor theories, whose Lagrangians contain up to two derivatives with respect to metric and vector field. We derive background equations under the assumption of homogeneous and isotropic universe and study the nature of cosmological perturbations on top of the background. We found an example of healthy cosmological solutions in non-trivial degenerate theory of gravity for the first time, where those perturbations do not suffer from any instabilities, that is, ghost and gradient instabilities.

Recently, several extensions of massive vector theory in curved space-time have been proposed in many literatures. In this paper, we consider the most general vector-tensor theories that contain up to two derivatives with respect to metric and vector field. By imposing a degeneracy condition of the Lagrangian in the context of ADM decomposition of space-time to eliminate an unwanted mode, we construct a new class of massive vector theories where five degrees of freedom can propagate, corresponding to three for massive vector modes and two for massless tensor modes. We find that the generalized Proca and the beyond generalized Proca theories up to the quartic Lagrangian, which should be included in this formulation, are degenerate theories even in curved space-time. Finally, introducing new metric and vector field transformations, we investigate the properties of thus obtained theories under such transformations.

The invariance of physical observables under disformal transformations is considered. It is known that conformal transformations leave physical observables invariant. However, whether it is true for disformal transformations is still an open question. In this paper, it is shown that a pure disformal transformation without any conformal factor is equivalent to rescaling the time coordinate. Since this rescaling applies equally to all the physical quantities, physics must be invariant under a disformal transformation, that is, neither causal structure, propagation speed nor any other property of the fields are affected by a disformal transformation itself. This fact is presented at the action level for gravitational and matter fields and it is illustrated with some examples of observable quantities. We also find the physical invariance for cosmological perturbations at linear and high orders in perturbation, extending previous studies. Finally, a comparison with Horndeski and beyond Horndeski theories under a disformal transformation is made.

The Square Kilometre Array will revolutionize pulsar studies with its wide field-of-view, wide-band observation and high sensitivity, increasing the number of observable pulsars by more than an order of magnitude. Pulsars are of interest not only for the study of neutron stars themselves but for their usage as tools for probing fundamental physics such as general relativity, gravitational waves and nuclear interaction. In this article, we summarize the activity and interests of SKA-Japan Pulsar Science Working Group, focusing on an investigation of modified gravity theory with the supermassive black hole in the Galactic Centre, gravitational-wave detection from cosmic strings and binary supermassive black holes, a study of the physical state of plasma close to pulsars using giant radio pulses and determination of magnetic field structure of Galaxy with pulsar pairs.

We study metric transformations which depend on a scalar field $\phi$ and its first derivatives and confirm that the number of physical degrees of freedom does not change under such transformations, as long as they are not singular. We perform a Hamiltonian analysis of a simple model in the gauge $\phi = t$. In addition, we explicitly show that the transformation and the gauge fixing do commute in transforming the action. We then extend the analysis to more general gravitational theories and transformations in general gauges. We verify that the set of all constraints and the constraint algebra are left unchanged by such transformations and conclude that the number of degrees of freedom is not modified by a regular and invertible generic transformation among two metrics. We also discuss the implications on the recently called "hidden" constraints and on the case of a singular transformation, a.k.a. mimetic gravity.

We study disformal transformations of the metric in the cosmological context. We first consider the disformal transformation generated by a scalar field $\phi$ and show that the curvature and tensor perturbations on the uniform $\phi$ slicing, on which the scalar field is homogeneous, are non-linearly invariant under the disformal transformation. Then we discuss the transformation properties of the evolution equations for the curvature and tensor perturbations at full non-linear order in the context of spatial gradient expansion as well as at linear order. In particular, we show that the transformation can be described in two typically different ways: one that clearly shows the physical invariance and the other that shows an apparent change of the causal structure. Finally we consider a new type of disformal transformation in which a multi-component scalar field comes into play, which we call a "multi-disformal transformation". We show that the curvature and tensor perturbations are invariant at linear order, and also at non-linear order provided that the system has reached the adiabatic limit.

We give new constraints on small-scale non-Gaussianity of primordial curvature perturbations by the use of anisotropies in acoustic reheating. Mixing of local thermal or local kinetic equilibrium systems with different temperatures yields a locally averaged temperature rise, which is proportional to the square of temperature perturbations damping in the photon diffusion scale. Such secondary temperature perturbations are indistinguishable from the standard temperature perturbations linearly coming from primordial curvature perturbations and hence should be subdominant compared to the standard ones. We show that small-scale higher order correlation functions (connected non-Gaussian and disconnected Gaussian parts) of primordial curvature perturbations can be probed by investigating auto power spectrum of the generated secondary perturbations and the cross power spectrum with the standard perturbations. This is simply because these power spectra come from higher order correlation functions of primordial curvature perturbations with non-linear parameters such as $f_{\rm NL}$ and $\tau_{\rm NL}$ since secondary temperature perturbations are second order effects. Thus, the observational results $l(l+1)C^{TT}_l\simeq 6\times 10^{-10}$ at large scales give a robust and universal upper bound on small-scale non-Gaussianities of primordial curvature perturbations.

In this paper, we introduce a new approach to a treatment of the gravitational effects (redshift, time delay and lensing) on the observed cosmic microwave background (CMB) anisotropies based on the Boltzmann equation. From the Liouville's theorem in curved spacetime, the intensity of photons is conserved along a photon geodesic when non-gravitational scatterings are absent. Motivated by this fact, we derive a second-order line-of-sight formula by integrating the Boltzmann equation along a perturbed geodesic (curve) instead of a background geodesic (line). In this approach, the separation of the gravitational and intrinsic effects are manifest. This approach can be considered as a generalization of the remapping approach of CMB lensing, where all the gravitational effects can be treated on the same footing.

We investigate the cosmology of SO(3)-invariant massive gravity with 5 degrees of freedom. In contrast with previous studies, we allow for a non-trivial fiducial metric, which can be justified by invoking, for example, a dilaton-like global symmetry. We write the homogeneous and isotropic equations of motion in this more general setup and identify, in particular, de Sitter solutions. We then study the linear perturbations around the homogeneous cosmological solutions, by deriving the quadratic Lagrangian governing the dynamics of scalar, vector and tensor modes. We thus obtain the conditions for the perturbations to be well-behaved. We show that it is possible to find de Sitter solutions whose perturbations are weakly coupled and stable, i.e. without ghost-like or gradient instabilities.

In the framework of massive gravity with a de Sitter reference metric, we study homogeneous and isotropic solutions with positive spatial curvature. Remarkably, we find that bounces can occur when cosmological matter satisfies the strong energy condition, in contrast to what happens in classical general relativity. This is due to the presence in the Friedmann equations of additional terms, which depend on the scale factor and its derivatives and can be interpreted as an effective fluid. We present a detailed study of the system using a phase space analysis. After having identified the fixed points of the system and investigated their stability properties, we discuss the cosmological evolution in the global physical phase space. We find that bouncing solutions

In the context of cosmological perturbation theory, we derive the second order Boltzmann equation describing the evolution of the distribution function of radiation without a specific gauge choice. The essential steps in deriving the Boltzmann equation are revisited and extended given this more general framework: i) the polarisation of light is incorporated in this formalism by using a tensor-valued distribution function; ii) the importance of a choice of the tetrad field to define the local inertial frame in the description of the distribution function is emphasized; iii) we perform a separation between temperature and spectral distortion, both for the intensity and for polarisation for the first time; iv) the gauge dependence of all perturbed quantities that enter the Boltzmann equation is derived, and this enables us to check the correctness of the perturbed Boltzmann equation by explicitly showing its gauge-invariance for both intensity and polarization. We finally discuss several implications of the gauge dependence for the observed temperature.

We develop a theory of nonlinear cosmological perturbations on superhorizon scales for a multi-component scalar field with a general kinetic term and a general form of the potential in the context of inflationary cosmology. We employ the ADM formalism and the spatial gradient expansion approach, characterised by O(\epsilon^2), where \epsilon=1/(HL) is a small parameter representing the ratio of the Hubble radius to the characteristic length scale L of perturbations. We provide a formalism to obtain the solution in the multi-field case. This formalism can be applied to the superhorizon evolution of a primordial non-Gaussianity beyond the so-called \delta N formalism which is equivalent to O(\epsilon^0) of the gradient expansion. In doing so, we also derive fully nonlinear gauge transformation rules valid through O(\epsilon^2). These fully nonlinear gauge transformation rules can be used to derive the solution in a desired gauge from the one in a gauge where computations are much simpler. As a demonstration, we consider an analytically solvable model and construct the solution explicitly.

In the framework of the recently proposed models of massive gravity, defined with respect to a de Sitter reference metric, we obtain new homogeneous and isotropic solutions for arbitrary cosmological matter and arbitrary spatial curvature. These solutions can be classified into three branches. In the first two, the massive gravity terms behave like a cosmological constant. In the third branch, the massive gravity effects can be described by a time evolving effective fluid with rather remarkable features, including the property to behave as a cosmological constant at late time.

Recently, the equivalence between the \delta N and covariant formalisms has been shown (Suyama et al. 2012), but they essentially assumed Einstein gravity in their proof. They showed that the evolution equation of the curvature covector in the covariant formalism on uniform energy density slicings coincides with that of the curvature perturbation in the \delta N formalism assuming the coincidence of uniform energy and uniform expansion (Hubble) slicings, which is the case on superhorizon scales in Einstein gravity. In this short note, we explicitly show the equivalence between the \delta N and covariant formalisms without specifying the slicing condition and the associated slicing coincidence, in other words, regardless of the gravity theory.

Open inflation scenario is attracting a renewed interest in the context of string landscape. Since there are a large number of metastable de Sitter vacua in string landscape, tunneling transitions to lower metastable vacua through the bubble nucleation occur quite naturally. Although the deviation of Omega_0 from unity is small by the observational bound, we argue that the effect of this small deviation on the large angle CMB anisotropies can be significant for tensor-type perturbation in open inflation scenario. We consider the situation in which there is a large hierarchy between the energy scale of the quantum tunneling and that of the slow-roll inflation in the nucleated bubble. If the potential just after tunneling is steep enough, a rapid-roll phase appears before the slow-roll inflation. In this case the power spectrum is basically determined by the Hubble rate during the slow-roll inflation. If such rapid-roll phase is absent, the power spectrum keeps the memory of the high energy density there in the large angular components. The amplitude of large angular components can be enhanced due to the effects of the wall fluctuation mode if the bubble wall tension is small. Therefore, one can construct some models in which the deviation of Omega_0 from unity is large enough to produce measurable effects. We also consider a more general class of models, where the false vacuum decay may occur due to Hawking-Moss tunneling, as well as the models involving more than one scalar field. We discuss scalar perturbations in these models and point out that a large set of such models is already ruled out by observational data, unless there was a very long stage of slow-roll inflation after the tunneling. These results show that observational data allow us to test various assumptions concerning the structure of the string theory potentials and the duration of the last stage of inflation.

It is known that the curvature perturbation on uniform energy density (or comoving or uniform Hubble) slices on superhorizon scales is conserved to full nonlinear order if the pressure is only a function of the energy density (ie, if the perturbation is purely adiabatic), independent of the gravitational theory. Here we explicitly show that the same conservation holds for a universe dominated by a single scalar field provided that the field is in an attractor regime, for a very general class of scalar field theories. However, we also show that if the scalar field equation contains a second time derivative of the metric, as in the case of the Galileon theory, one has to invoke the gravitational field equations to show the conservation.

We compute the one-point probability distribution function of small-angle cosmic microwave background temperature fluctuations due to curved cosmic (super-)strings with a simple model of string network by performing Monte Carlo simulations. Taking into account of the correlation between the curvature and the velocity of string segments, there appear non-Gaussian features, specifically non-Gaussian tails and a skewness, in the one-point pdf. The obtained sample skewness for the conventional field-theoretic cosmic strings is $g_1\approx -0.14$, which is consistent with the result reported by Fraisse et al.. We also discuss the dependence of the pdf on the intercommuting probability. We find that the standard deviation of the Gaussian part increases and non-Gaussian features are suppressed as the intercommuting probability decreases. For sufficiently small intercommuting probability, the skewness is given by $\lesssim$ (a\ few) $\times 10^{-2}$.

We compute analytically the small-scale temperature fluctuations of the cosmic microwave background from cosmic (super-)strings and study the dependence on the string intercommuting probability $P$. We develop an analytical model which describes the evolution of a string network and calculate the numbers of string segments and kinks in a horizon volume. Then we derive the probability distribution function (pdf) which takes account of finite angular resolution of observation. The resultant pdf consists of a Gaussian part due to frequent scatterings by long string segments and a non-Gaussian tail due to close encounters with kinks. The dispersion of the Gaussian part is reasonably consistent with that obtained by numerical simulations by Fraisse et al.. On the other hand, the non-Gaussian tail contains two phenomenological parameters which are determined by comparison with the numerical results for P=1. Extrapolating the pdf to the cases with $P<1$, we predict that the non-Gaussian feature is suppressed for small $P$.

A model of multi-component hybrid inflation, dubbed multi-brid inflation, in which various observable quantities including the non-Gaussianity parameter f_{NL} can be analytically calculated was proposed recently. In particular, for a two-brid inflation model with an exponential potential and the condition that the end of inflation is an ellipse in the field space, it was found that, while keeping the other observational quantities within the range consistent with observations, large non-Gaussianity is possible for certain inflationary trajectories, provided that the ratio of the two masses is large. One might question whether the resulting large non-Gaussianity is specific to this particular form of the potential and the condition for the end of inflation. In this paper, we consider a model of multi-brid inflation with a potential given by an exponential function of terms quadratic in the scalar field components. We also consider a more general class of ellipses for the end of inflation than those studied previously. Then, focusing on the case of two-brid inflation, we find that large non-Gaussianity is possible in the present model even for the equal-mass case. Then by tuning the model parameters, we find that there exist models for which both the non-Gaussianity and the tensor-to-scalar ratio are large enough to be detected in the very near future.