In a class of generalized Proca theories, we study the existence of neutron star solutions with a nonvanishing temporal component of the vector field $A_\mu$ approaching 0 toward spatial infinity, as they may be the endpoints of tachyonic instabilities of neutron star solutions in general relativity with $A_{\mu}=0$. Such a phenomenon is called spontaneous vectorization, which is analogous to spontaneous scalarization in scalar-tensor theories with nonminimal couplings to the curvature or matter. For the nonminimal coupling $\beta X R$, where $\beta$ is a coupling constant and $X=-A_{\mu}A^{\mu}/2$, we show that there exist both 0-node and 1-node vector-field solutions, irrespective of the choice of the equations of state of nuclear matter. The 0-node solution, which is present only for $\beta=-{\cal O}(0.1)$, may be induced by some nonlinear effects such as the selected choice of initial conditions. The 1-node solution exists for $\beta=-{\cal O}(1)$, which suddenly emerges above a critical central density of star and approaches the general relativistic branch with the increasing central density. We compute the mass $M$ and radius $r_s$ of neutron stars for some realistic equations of state and show that the $M$-$r_s$ relations of 0-node and 1-node solutions exhibit notable difference from those of scalarized solutions in scalar-tensor theories. Finally, we discuss the possible endpoints of tachyonic instabilities.

The canonical scalar-tensor theory model which exhibits spontaneous scalarization in the strong-gravity regime of neutron stars has long been known to predict a cosmological evolution for the scalar field which generically results in severe violations of present-day Solar System constraints on deviations from general relativity. We study if this tension can be alleviated by generalizing this model to include a disformal coupling between the scalar field $\varphi$ and matter, where the Jordan frame metric ${\tilde g}_{\mu\nu}$ is related to the Einstein frame one $g_{\mu\nu}$ by ${\tilde g}_{\mu\nu}=A(\varphi)^2 (g_{\mu\nu}+\Lambda\, \partial_\mu \varphi \, \partial_\nu\varphi)$. We find that this broader theory admits a late-time attractor mechanism towards general relativity. However, the existence of this attractor requires a value of disformal scale of the order $\Lambda\gtrsim H_0^{-2}$, where $H_0$ is the Hubble parameter of today, which is much larger than the scale relevant for spontaneous scalarization of neutron stars $\Lambda \sim R_s^{2}$ with $R_s (\sim 10^{-22} H_0^{-1})$ being the typical radius of these stars. The large values of $\Lambda$ necessary for the attractor mechanism (i) suppress spontaneous scalarization altogether inside neutron stars and (ii) induce ghost instabilities on scalar field fluctuations, thus preventing a resolution of the tension. We argue that the problem arises because our disformal coupling involves a dimensionful parameter.

We find various exact black hole solutions in the shift-symmetric subclass of the quadratic degenerate higher-order scalar-tensor (DHOST) theories with linearly time dependent scalar field whose kinetic term is constant. The exact solutions are the Schwarzschild and Schwarzschild-(anti-)de Sitter solutions, and the Schwarzschild-type solution with a deficit solid angle, which are accompanied by nontrivial scalar field regular at the black hole event horizon. We derive the conditions for the coupling functions in the DHOST Lagrangian that allow the exact solutions, clarify their compatibility with the degeneracy conditions, and provide general form of coupling functions as well as simple models that satisfy the conditions.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

The generalized Proca theories with second-order equations of motion can be healthily extended to a more general framework in which the number of propagating degrees of freedom remains unchanged. In the presence of a quartic-order nonminimal coupling to gravity arising in beyond-generalized Proca theories, the speed of gravitational waves $c_t$ on the Friedmann-Lemaitre-Robertson-Walker (FLRW) cosmological background can be equal to that of light $c$ under a certain condition. By using this condition alone, we show that the speed of gravitational waves in the vicinity of static and spherically symmetric black holes is also equivalent to $c$ for the propagation of odd-parity perturbations along both radial and angular directions. As a by-product, the black holes arising in our beyond-generalized Proca theories are plagued by neither ghost nor Laplacian instabilities against odd-parity perturbations. We show the existence of both exact and numerical black hole solutions endowed with vector hairs induced by the quartic-order coupling.

Recent gravitational wave observations of binary black hole mergers and a binary neutron star merger by LIGO and Virgo Collaborations associated with its optical counterpart constrain deviation from General Relativity (GR) both on strong-field regime and cosmological scales with high accuracy, and further strong constraints are expected by near-future observations. Thus, it is important to identify theories of modified gravity that intrinsically possess the same solutions as in GR among a huge number of theories. We clarify the three conditions for theories of modified gravity to allow GR solutions, i.e., solutions with the metric satisfying the Einstein equations in GR and the constant profile of the scalar fields. Our analysis is quite general, as it applies a wide class of single-/multi-field scalar-tensor theories of modified gravity in the presence of matter component, and any spacetime geometry including cosmological background as well as spacetime around black hole and neutron star, for the latter of which these conditions provide a necessary condition for no-hair theorem. The three conditions will be useful for further constraints on modified gravity theories as they classify general theories of modified gravity into three classes, each of which possesses i) unique GR solutions (i.e., no-hair cases), ii) only hairy solutions (except the cases that GR solutions are realized by cancellation between singular coupling functions in the Euler-Lagrange equations), and iii) both GR and hairy solutions, for the last of which one of the two solutions may be selected dynamically.

We study relativistic star solutions in second-order generalized Proca theories characterized by a $U(1)$-breaking vector field with derivative couplings. In the models with cubic and quartic derivative coupling, the mass and radius of stars become larger than those in general relativity for negative derivative coupling constants. This phenomenon is mostly attributed to the increase of star radius induced by a slower decrease of the matter pressure compared to general relativity. There is a tendency that the relativistic star with a smaller mass is not gravitationally bound for a low central density and hence dynamically unstable, but that with a larger mass is gravitationally bound. On the other hand, we show that the intrinsic vector-mode couplings give rise to general relativistic solutions with a trivial field profile, so the mass and radius are not modified from those in general relativity.

In generalized Proca theories with vector-field derivative couplings, a bunch of hairy black hole solutions have been derived on a static and spherically symmetric background. In this paper, we formulate the odd-parity black hole perturbations in generalized Proca theories by expanding the corresponding action up to second order and investigate whether or not black holes with vector hair suffer ghost or Laplacian instabilities. We show that the models with cubic couplings $G_3(X)$, where $X=-A_{\mu}A^{\mu}/2$ with a vector field $A_{\mu}$, do not provide any additional stability condition as in General Relativity. On the other hand, the exact charged stealth Schwarzschild solution with a nonvanishing longitudinal vector component $A_1$, which originates from the coupling to the Einstein tensor $G^{\mu\nu}A_\mu A_\nu$ equivalent to the quartic coupling $G_4(X)$ containing a linear function of $X$, is unstable in the vicinity of the event horizon. The same instability problem also persists for hairy black holes arising from general quartic power-law couplings $G_4(X) \supset \beta_4 X^n$ with the nonvanishing $A_1$, while the other branch with $A_1=0$ can be consistent with conditions for the absence of ghost and Laplacian instabilities. We also discuss the case of other exact and numerical black hole solutions associated with intrinsic vector-field derivative couplings and show that there exists a wide range of parameter spaces in which the solutions suffer neither ghost nor Laplacian instabilities against odd-parity perturbations.

We present a family of exact black-hole solutions on a static spherically symmetric background in second-order generalized Proca theories with derivative vector-field interactions coupled to gravity. We also derive non-exact solutions in power-law coupling models including vector Galileons and numerically show the existence of regular black holes with a primary hair associated with the longitudinal propagation. The intrinsic vector-field derivative interactions generally give rise to a secondary hair induced by non-trivial field profiles. The deviation from General Relativity is most significant around the horizon and hence there is a golden opportunity for probing the Proca hair by the measurements of gravitational waves (GWs) in the regime of strong gravity.

We study static and spherically symmetric black hole (BH) solutions in second-order generalized Proca theories with nonminimal vector field derivative couplings to the Ricci scalar, the Einstein tensor, and the double dual Riemann tensor. We find concrete Lagrangians which give rise to exact BH solutions by imposing two conditions of the two identical metric components and the constant norm of the vector field. These exact solutions are described by either Reissner-Nordstr\"{o}m (RN), stealth Schwarzschild, or extremal RN solutions with a non-trivial longitudinal mode of the vector field. We then numerically construct BH solutions without imposing these conditions. For cubic and quartic Lagrangians with power-law couplings which encompass vector Galileons as the specific cases, we show the existence of BH solutions with the difference between two non-trivial metric components. The quintic-order power-law couplings do not give rise to non-trivial BH solutions regular throughout the horizon exterior. The sixth-order and intrinsic vector-mode couplings can lead to BH solutions with a secondary hair. For all the solutions, the vector field is regular at least at the future or past horizon. The deviation from General Relativity induced by the Proca hair can be potentially tested by future measurements of gravitational waves in the nonlinear regime of gravity.

Horndeski gravity holds a special position as the most general extension of Einstein's theory of general relativity with a single scalar degree of freedom and second-order field equations. Because of these features, Horndeski gravity is an attractive phenomenological playground to investigate the consequences of modifications of general relativity in cosmology and astrophysics. We present a review of the progress made so far in the study of compact objects (black holes and neutron stars) within Horndeski gravity. In particular, we review our recent work on slowly rotating black holes and present some new results on slowly rotating neutron stars.

We present a general formulation to analyze the structure of slowly rotating relativistic stars in a broad class of scalar-tensor theories with disformal coupling to matter. Our approach includes theories with generalized kinetic terms, generic scalar field potentials and contains theories with conformal coupling as particular limits. In order to investigate how the disformal coupling affects the structure of relativistic stars, we propose a minimal model of a massless scalar-tensor theory and investigate in detail how the disformal coupling affects the spontaneous scalarization of slowly rotating neutron stars. We show that for negative values of the disformal coupling parameter between the scalar field and matter, scalarization can be suppressed, while for large positive values of the disformal coupling parameter stellar models cannot be obtained. This allows us to put a mild upper bound on this parameter. We also show that these properties can be qualitatively understood by linearizing the scalar field equation of motion in the background of a general-relativistic incompressible star. To address the intrinsic degeneracy between uncertainties in the equation-of-state of neutron stars and gravitational theory, we also show the existence of universal equation-of-state-independent relations between the moment of inertia and compactness of neutron stars in this theory. We show that in a certain range of the theory's parameter space the universal relation largely deviates from that of general relativity, allowing, in principle, to probe the existence of spontaneous scalarization with future observations.

Horndeski's theory of gravity is the most general scalar-tensor theory with a single scalar whose equations of motion contain at most second-order derivatives. A subsector of Horndeski's theory known as "Fab Four" gravity allows for dynamical self-tuning of the quantum vacuum energy, and therefore it has received particular attention in cosmology as a possible alternative to the $\Lambda$CDM model. Here we study compact stars in Fab Four gravity, which includes as special cases general relativity ("George"), Einstein-dilaton-Gauss-Bonnet gravity ("Ringo"), theories with a nonminimal coupling with the Einstein tensor ("John") and theories involving the double-dual of the Riemann tensor ("Paul"). We generalize and extend previous results in theories of the John class and were not able to find realistic compact stars in theories involving the Paul class.

We study black hole solutions at first order in the Hartle-Thorne slow-rotation approximation in Horndeski gravity theories. We derive the equations of motion including also cases where the scalar depends linearly on time. In the Hartle-Thorne formalism, all first-order rotational corrections are described by a single frame-dragging function. We show that the frame-dragging function is exactly the same as in general relativity for all known black hole solutions in shift symmetric Horndeski theories, with the exception of theories with a linear coupling to the Gauss-Bonnet invariant. Our results extend previous no-hair theorems for a broad class of Horndeski gravity theories.

We investigate the gauge-invariant cosmological perturbations in the gravity and matter frames in the general scalar-tensor theory where two frames are related by the disformal transformation. The gravity and matter frames are the extensions of the Einstein and Jordan frames in the scalar-tensor theory where two frames are related by the conformal transformation, respectively. First, it is shown that the curvature perturbation in the comoving gauge to the scalar field is disformally invariant as well as conformally invariant, which gives the predictions from the cosmological model where the scalar field is responsible both for inflation and cosmological perturbations. Second, in case that the disformally coupled matter sector also contributes to curvature perturbations, we derive the evolution equations of the curvature perturbation in the uniform matter energy density gauge from the energy (non)conservation in the matter sector, which are independent of the choice of the gravity sector. While in the matter frame the curvature perturbation in the uniform matter energy density gauge is conserved on superhorizon scales for the vanishing nonadiabatic pressure, in the gravity frame it is not conserved even if the nonadiabatic pressure vanishes. The formula relating two frames gives the amplitude of the curvature perturbation in the matter frame, once it is evaluated in the gravity frame.

Using the \delta N formalism we consider the non-linear curvature perturbation in multi-field models of inflation with non-minimal coupling. In particular, we focus on the relation between the \delta N formalism as applied in the conformally related Jordan and Einstein frames. Exploiting results already known in the Einstein frame, we give expressions for the power spectrum, spectral tilt and non-gaussianity associated with the Jordan frame curvature perturbation. In the case that an adiabatic limit has not been reached, we find that in general these quantities differ from those associated with the Einstein frame curvature perturbation, and also confirm their equivalence in the absence of isocurvature modes. We then proceed to consider two analytically soluble examples, the first involving a non-minimally coupled `spectator' field and the second being a non-minimally coupled extension of the multi-brid inflation model. In the first model we find that predictions can easily be brought into agreement with the recent Planck results, as the tensor-to-scalar ratio is generally small, the spectral tilt tuneable and the non-gaussianity suppressed. In the second model we find that predictions for all three parameters can differ substantially from those predicted in the minimally coupled case, and that the recent Planck results for the spectral tilt can be used to constrain the non-minimal coupling parameters.

We revisit the issue on signatures of pre-inflationary background anisotropy by considering the quantization of a massless and minimally coupled scalar field in an axially symmetric Kasner background, mimicking cosmological perturbations. We show that the power spectrum of the scalar field fluctuation has a negligible difference from the standard inflation in the non-planar directions, but it has a sharp peak around the symmetry plane. For the non-planar high-momentum modes, we use the WKB approximation for the first period and the asymptotic approximation based on the de Sitter solution for the next period. At the boundary, two mode functions have the same accuracy with error of $O({H_{i}}/k)$. We calculate the approximation up to the order of $({H_{i}}/k)^6$ and show that the power spectrum of the scalar field fails to get corrections until we execute the approximation up to $6^{\rm th}$ order.

In this paper we discuss a multi-field model of inflation in which generally all fields are non-minimally coupled to the Ricci scalar and have non-canonical kinetic terms. The background evolution and first-order perturbations for the model are evaluated in both the Jordan and Einstein frames, and the respective curvature perturbations compared. We confirm that they are indeed not the same - unlike in the single-field case - and also that the difference is a direct consequence of the isocurvature perturbations inherent to multi-field models. This result leads us to conclude that the notion of adiabaticity is not invariant under conformal transformations. Using a two-field example we show that even if in one frame the evolution is adiabatic, meaning that the curvature perturbation is conserved on super-horizon scales, in general in the other frame isocurvature perturbations continue to source the curvature perturbation. We also find that it is possible to realise a particular model in which curvature perturbations in both frames are conserved but with each being of different magnitude. These examples highlight that the curvature perturbation itself, despite being gauge-invariant, does not correspond directly to an observable. The non-equivalence of the two curvature perturbations would also be important when considering the addition of Standard Model matter into the system.

In this paper, we present the higher order spectra of a scalar field produced through the higher derivative interactions in the initially anisotropic universe. Although we ignore the backreaction of the scalar field on the geometry, our analysis should have much overlap with the quantum fluctuations of the inflaton field in the anisotropic universe. We also include the planar modes whose momenta are along the plane which is perpendicular to the primordial preferred direction, for which effects of the initial anisotropy are not suppressed. The presence of a negative frequency mode produces features distinguishable from the case of the de Sitter inflation. We also show that richer features appear in the trispectra due to the primordial anisotropy.

We present the analytic forms for the spectra of the cosmological perturbations from an initially anisotropic universe for the high momentum modes in the context of WKB approximations, as the continuation of the work [29]. We consider the Einstein gravity coupled to a light scalar field. We then assume that the scalar field has the zero velocity initially and then slowly rolls down on the potential toward the origin. In the slow-roll approximations, the Kasner-de Sitter universe with a planar symmetry is a good approximation as the background evolution. Quantization of the perturbations in the adiabatic vacuum, which we call the anisotropic vacuum, is carried out. For non-planar high momentum modes whose comoving momentum component orthogonal to the plane is bigger than the Hubble parameter at the inflationary phase, the WKB approximation is valid for the whole stage of the isotropization. On the other hand, the planar modes whose comoving momentum component orthogonal to the plane is comparable to the Hubble parameter, is amplified during the process of the anisotropic expansion. In the final gravitational wave spectra, we find that there is an asymmetry between the two polarizations of the gravitational wave because the initial mode mixing does not vanish.

We discuss the hybrid inflation model where the inflaton field is nonminimally coupled to gravity. In the Jordan frame, the potential contains $\phi^4$ term as well as terms in the original hybrid inflation model. In our model, inflation can be classified into the type (I) and the type (II). In the type (I), inflation is terminated by the tachyonic instability of the waterfall field, while in the type (II) by the violation of slow-roll conditions. In our model, the reheating takes place only at the true minimum and even in the case (II) finally the tachyonic instability occurs after the termination of inflation. For a negative nonminimal coupling, inflation takes place in the vacuum-dominated region, in the large field region, or near the local minimum/maximum. Inflation in the vacuum dominated region becomes either the type (I) or (II), resulting in blue or red spectrum of the curvature perturbations, respectively. Inflation around the local maximum can be either the type (I) or the type (II), which results in the red spectrum of the curvature perturbations, while it around the local minimum must be the type (I), which results in the blue spectrum. In the large field region, to terminate inflation, potential in the Einstein frame must be positively tilted, always resulting in the red spectrum. We then numerically solve the equations of motion to investigate the whole dynamics of inflaton and confirm that the spectrum of curvature perturbations changes from red to blue ones as scales become smaller.

We discuss the primordial spectrum of a massless and minimally coupled scalar field, produced during the initial anisotropic epoch before the onset of inflation. We consider two models of the anisotropic cosmology, the (planar) Kasner de Sitter solution (Bianchi I) and the Taub-NUT de Sitter solution (Bianchi IX), where the 3-space geometry is initially anisotropic, followed by the de Sitter phase due to the presence of a positive cosmological constant. We discuss the behavior of a quantized, massless and minimally coupled scalar field in the anisotropic stage. We focus on the quantization procedure and evolution in the pre-inflationary anisotropic background. The initial condition is set by the requirement that the scalar field is initially in an adiabatic state. In our case, such a vacuum state exists only for limited solutions of the anisotropic Universe, whose spacetime structure is regular in the initial times. In this paper, we call our adiabatic vacuum state the anisotropic vacuum. In the Kasner de Sitter model, for one branch of planar solutions there is an anisotropic vacuum unless $k_3\neq 0$, where $k_3$ is the comoving momentum along the third direction, while in the other branch there is no anisotropic vacuum state. In the first branch, for the moderate modes, $k_3\sim k$, where $k$ is the total comoving momentum, the scalar power spectrum has an oscillatory behavior and its direction dependence is suppressed. For the planar modes, $k_3\ll k$, in contrast, the direction dependence becomes more important, because of the amplification of the scalar amplitude during this interval of the violation of WKB approximation in the initial anisotropic stage. The qualitative behaviors in the Taub-NUT de Sitter models are very similar to the case of the first branch of the planar Kasner de Sitter model.

In this short note, we are going to classify the cosmological evolutions in the recently proposed, nonrelativistic gravitational theory, the Ho\v{r}ava-Lifshitz (HL) theory. We consider the original HL theory (theory I), and the modified version obtained by an analytic continuation of parameters (theory II). We discuss the possiblilities of cosmological evolutions with arbitrary matter.

We present thick de Sitter brane solutions which are supported by two interacting {\it phantom} scalar fields in five-, six- and seven-dimensional spacetime. It is shown that for all cases regular solutions with anti-de Sitter asymptotic (5D problem) and a flat asymptotic far from the brane (6D and 7D cases) exist. We also discuss the stability of our solutions.

We discuss the comsological evolution of a brane in the $D(>6)$-dimensional black brane spacetime in the context of the Kaluza-Klein (KK) braneworld scheme, i.e., to consider KK compactification on the brane. The bulk spacetime is composed of two copies of a patch of $D$-dimensional black three-brane solution. The near-horizon geometry is given by $AdS_{5}\times S^{(D-5)}$ while in the asymptotic infinity the spacetime approaches $D$-dimensional Minkowski. We consider the brane motion from the near-horizon region toward the spatial infinity, which induces cosmology on the brane. As is expected, in the early times, namely when the brane is located in the near-horizon region, the effective cosmology on the brane coincides with that in the second Randall-Sundrum (RS II) model. Then, the brane cosmology starts to deviate from the RS type one since the dynamics of KK compactified dimensions becomes significant. We find that the brane Universe cannot reach the asymptotic infinity, irrespectively of the components of matter on the brane.

We discuss the stability of the higher-dimensional de Sitter (dS) brane solutions with two-dimensional internal space in the Einstein-Maxwel theory. We show that an instability appears in the scalar-type perturbations with respect to the dS spacetime. We derive a differential relation which has the very similar structure to the ordinary laws of thermodynamics as an extension of the work for the six-dimensional model [20]. In this relation, the area of dS horizon (integrated over the two internal dimensions) exactly behaves as the thermodynamical entropy. The dynamically unstable solutions are in the thermodynamically unstable branch. An unstable dS compactification either evolves toward a stable configuration or two-dimensional internal space is decompactified. These dS brane solutions are equivalent to the accelerating cosmological solutions in the six-dimensional Einstein-Maxwell-dilaton theory via dimensional reduction. Thus, if the seed higher-dimensional solution is unstable, the corresponding six-dimensional solution is also unstable. From the effective four-dimensional point of view, a cosmological evolution from an unstable cosmological solution in higher dimensions may be seen as a process of the transition from the initial cosmological inflation to the current dark energy dominated Universe.

We study the cosmological evolution of extended branes in 6D warped flux compactification models. The branes are endowed with the three ordinary spatial dimensions, which are assumed to be homogeneous and isotropic, as well as an internal extra dimension compactified on a circle. We embed these codimension 1 branes in a static bulk 6D spacetime, whose geometry is a solution of 6D Einstein-Maxwell or Einstein-Maxwell-dilaton theories, corresponding to a warped flux compactification. The brane matter consists of a complex scalar field which is coupled to the bulk U(1) gauge field. In both models, we show that there is critical point which the brane cannot cross as it moves in the bulk. We study the cosmological behaviour, especially when the brane approaches this critical point or one of the two conical singularities. In the present setup where the bulk geometry is fixed, we find that the brane cosmology does not coincide with the standard one in the low energy limit.

We study cosmology on a conical brane in the six-dimensional Einstein-Maxwell-dilaton system, where the extra dimensions are compactified by a magnetic flux. We systematically construct exact cosmological solutions using the fact that the system is equivalently described by (6+n)-dimensional pure Einstein-Maxwell theory via dimensional reduction. In particular, we find a power-law inflationary solution for a general dilatonic coupling. When the dilatonic coupling is given by that of Nishino-Sezgin chiral supergravity, this reduces to the known solution which is not inflating. The power-law solution is shown to be the late-time attractor. We also investigate cosmological tensor perturbations in this model using the (6+n)-dimensional description. We obtain the separable equation of motion and find that there always exist a zero mode, while tachyonic modes are absent in the spectrum. The mass spectrum of Kaluza-Klein modes is obtained numerically.

We discuss Casimir effect of a massless, minimally coupled scalar field in a 6D warped flux compactification model and its implications for the hierarchy and cosmological constant problems, which are longstanding puzzles in phenomenology and cosmology. Due to the scale invariance of the background theory, the 4D effective theory contains a volume modulus. To stabilize the modulus, we consider one-loop corrections to the effective potential by the Casimir effect. The one-loop effective potential for the volume modulus has a form which is very similar to Coleman-Weinberg potential. We evaluate coefficients appearing in the effective potential by employing zeta function regularization and heat kernel analysis. The volume modulus is stabilized for smaller degrees of warping, below a critical value, which depends on deficit angle of the reference brane. After stabilizing the modulus, it is possible to obtain observed values of the mass ratio between the fundamental energy scales and a tiny effective cosmological constant (though its sign is negative). The degree of warping should be tuned to be close to the critical value, not as severely as the original fine-tuning.

As a continuation of the work in \cite{mns}, we discuss the Casimir effect for a massless bulk scalar field in a 4D toy model of a 6D warped flux compactification model,to stabilize the volume modulus. The one-loop effective potential for the volume modulus has a form similar to the Coleman-Weinberg potential. The stability of the volume modulus against quantum corrections is related to an appropriate heat kernel coefficient. However, to make any physical predictions after volume stabilization, knowledge of the derivative of the zeta function, $\zeta'(0)$ (in a conformally related spacetime) is also required. By adding up the exact mass spectrum using zeta function regularization, we present a revised analysis of the effective potential. Finally, we discuss some physical implications, especially concerning the degree of the hierarchy between the fundamental energy scales on the branes. For a larger degree of warping our new results are very similar to the previous ones \cite{mns} and imply a larger hierarchy. In the non-warped (rugby-ball) limit the ratio tends to converge to the same value, independently of the bulk dilaton coupling.

We study scalar cosmological perturbations in a braneworld model with a bulk Gauss-Bonnet term. For an anti-de Sitter bulk, the five-dimensional perturbation equations share the same form as in the Randall-Sundrum model, which allows us to obtain metric perturbations in terms of a master variable. We derive the boundary conditions for the master variable from the generalized junction conditions on the brane. We then investigate several limiting cases in which the junction equations are reduced to a feasible level. In the low energy limit, we confirm that the standard result of four-dimensional Einstein gravity is reproduced on large scales, whereas on small scales we find that the perturbation dynamics is described by the four-dimensional Brans-Dicke theory. In the high energy limit, all the non-local contributions drop off from the junction equations, leaving a closed system of equations on the brane. We show that, for inflation models driven by a scalar field on the brane, the Sasaki-Mukhanov equation holds on the high energy brane in its original four-dimensional form.

We investigate quantum fluctuations on a de Sitter (dS) brane, which has its own thickness, in order to examine whether or not the finite thickness of the brane can act as a natural cut-off for the Kaluza-Klein (KK) spectrum. We calculate the amplitude of the KK modes and the bound state by using the zeta function method after a dimensional reduction.We show that the KK amplitude is finite for a given brane thickness and in the thin wall limit the standard surface divergent behavior is recovered. The strength of the divergence in the thin wall limit depends on the number of dimensions, e.g., logarithmic on a two dimensional brane and quadratic on a four dimensional brane. We also find that the amplitude of the bound state mode and KK modes depends on the choice of renormalization scale; and for fixed renormalization scales the bound state mode is insensitive to the brane thickness both for two and four-dimensional dS branes.

We discuss the backreaction of a massless, minimally coupled, quantized scalar field on a thick, two-dimensional de Sitter (dS) brane as an extension of our previous work. We show that a finite brane thickness naturally regularizes the backreaction on the brane. The quantum backreaction exhibits a quadratic divergence in the thin wall limit. We also give a theoretical bound on the brane thickness, in terms of the brane self-consistency of the quantum corrected Einstein equation, namely the requirement that the size of the backreaction should be smaller than that of the background stress-energy at the center of the brane. Finally, we discuss the brane self-consistency for the case of a four-dimensional dS brane.

We discuss the effect of Kaluza-Klein (KK) modes of bulk metric perturbations on the second Randall-Sundrum (RS II) type brane cosmology, taking the possible backreaction in the bulk and on the brane into account. KK gravitons may be produced via quantum fluctuations during a de Sitter (dS) inflating phase of our brane universe. In an effective 4-dimensional theory in which one integrates out the extra-dimensional dependence in the action, KK gravitons are equivalent to massive gravitons on the brane with masses $m>3H/2$, where $H$ represents the expansion rate of a dS brane. Thus production of even a tiny amount of KK gravitons may eventually have a significant impact on the late-time brane cosmology. As a first step to quantify the effect of KK gravitons on the brane, we calculate the effective energy density and pressure for a single KK mode. Surprisingly, we find that a KK mode behaves as cosmic dust with a negative energy density on the brane. We note that the bulk energy density of a KK mode is positive definite and there occurs no singular phenomenon in the bulk.

We investigate the linearized gravity on a single de Sitter brane in the anti-de Sitter (AdS) bulk in the Einstein Gauss-Bonnet (EGB) theory. We find that the Einstein gravity is recovered for a high energy brane, i.e., in the limit of the large expansion rate, $H\ell\gg 1$, where $H$ is the de Sitter expansion rate and $\ell$ is the curvature radius of the AdS bulk. We also show that, in the short distance limit $r\ll \min\{\ell,H^{-1}\}$, the Brans-Dicke gravity is obtained, whereas in the large distance limit $r\gg\max\{\ell,H^{-1}\}$, the Brans-Dicke gravity is obtained for $H\ell=O(1)$, and the Einstein gravity is recovered both for $H\ell\gg1$ and $H\ell\ll1$. In the limit $H\ell\to0$, these results smoothly match to the results known for the Minkowski brane.

In the context of the Randall-Sundrum (RS) single-brane scenario, we discuss the bulk geometry and dynamics of a cosmological brane in terms of the local energy conservation law which exists for the bulk that allows slicing with a maximally symmetric 3-space. This conservation law enables us to define a local mass in the bulk. We show that there is a unique generalization of the dark radiation on the brane, which is given by the local mass. We find there also exists a conserved current associated with the Weyl tensor, and the corresponding local charge, which we call the Weyl charge, is given by the sum of the local mass and a certain linear combination of the components of the bulk energy-momentum tensor. This expression of the Weyl charge relates the local mass with the projected Weyl tensor, $E_{\mu\nu}$, which plays a central role in the geometrical formalism of the RS braneworld. On the brane, in particular, this gives a decomposition of the projected Weyl tensor into the local mass and the bulk energy-momentum tensor. Then, as an application of these results, we consider a null dust model for the bulk energy-momentum tensor and discuss the black hole formation in the bulk. We investigate the causal structure by identifying the locus of the apparent horizon and clarify possible brane trajectories in the bulk. We find that the brane stays always outside the black hole as long as it is expanding. We also find an upper bound on the value of the Hubble parameter in terms of the matter energy density on the brane, irrespective of the energy flux emitted from the brane.

We consider a braneworld inflation model driven by the dynamics of a scalar field living in the 5-dimensional bulk, the so-called ``bulk inflaton model'', and investigate the geometry in the bulk and large scale cosmological perturbations on the brane. The bulk gravitational effects on the brane are described by a projection of the 5-dimensional Weyl tensor, which we denote by $E_{\mu\nu}$. Focusing on a tachionic potential model, we take a perturbative approach in the anti-de Sitter (AdS$_5$) background with a single de Sitter brane. We first formulate the evolution equations for $E_{\mu\nu}$ in the bulk. Next, applying them to the case of a spatially homogeneous brane, we obtain two different integral expressions for $E_{\mu\nu}$. One of them reduces to the expression obtained previously when evaluated on the brane. The other is a new expression that may be useful for analyzing the bulk geometry. Then we consider superhorizon scale cosmological perturbations and evaluate the bulk effects onto the brane. In the limit $H^2\ell^2\ll1$, where $H$ is the Hubble parameter on the brane and $\ell$ is the bulk curvature radius, we find that the effective theory on the brane is identical to the 4-dimensional Einstein-scalar theory with a simple rescaling of the potential even under the presence of inhomogeneities. % atleast on super-Hubble horizon scales. In particular, it is found that the anticipated non-trivial bulk effect due to the spatially anisotropic part of $E_{\mu\nu}$ may appear only at %second order in the low energy expansion, i.e., at $O(H^4\ell^4)$.