We study the frame dependence/independence of cosmological observables under disformal transformations, extending the previous results regarding conformal transformations, and provide the correspondence between Jordan-frame and Einstein-frame variables. We consider quantities such as the gravitational constant in the Newtonian limit, redshift, luminosity and angular diameter distances, as well as the distance-duality relation. Also, the Boltzmann equation, the observed specific intensity, and the adiabaticity condition are discussed. Since the electromagnetic action changes under disformal transformations, photons in the Einstein frame no longer propagate along null geodesics. As a result, several quantities of cosmological interest are modified. Nevertheless, we show that the redshift is invariant and the distance-duality relation (the relation between the luminosity distance and the angular diameter distance) still holds in general spacetimes even though the reciprocity relation (the relation between two geometrical distances) is modified.

Redshifts of an astronomical body measured at multiple epochs (e.g., separated by 10 years) are different due to the cosmic expansion. This so-called Sandage-Loeb test offers a direct measurement of the expansion rate of the Universe. However, acceleration in the motion of Solar System with respect to the cosmic microwave background also changes redshifts measured at multiple epochs. If not accounted for, it yields a biased cosmological inference. To address this, we calculate the acceleration of Solar System with respect to the Local Group of galaxies to quantify the change in the measured redshift due to local motion. Our study is motivated by the recent determination of the mass of Large Magellanic Cloud (LMC), which indicates a significant fraction of the Milky Way mass. We find that the acceleration towards the Galactic Center dominates, which gives a redshift change of 7 cm/s in 10 years, while the accelerations due to LMC and M31 cannot be ignored depending on lines of sight. We create all-sky maps of the expected change in redshift and the corresponding uncertainty, which can be used to correct for this effect.

Stochastic gravitational wave backgrounds, predicted in many models of the early universe and also generated by various astrophysical processes, are a powerful probe of the Universe. The spectral shape is key information to distinguish the origin of the background since different production mechanisms predict different shapes of the spectrum. In this paper, we investigate how precisely future gravitational wave detectors can determine the spectral shape using single and broken power-law templates. We consider the detector network of Advanced-LIGO, Advanced-Virgo and KAGRA and the space-based gravitational-wave detector DECIGO, and estimate the parameter space which could be explored by these detectors. We find that, when the spectrum changes its slope in the frequency range of the sensitivity, the broken power-law templates dramatically improve the $\chi^2$ fit compared with the single power-law templates and help to measure the shape with a good precision.

Recent cosmological observations are in good agreement with the scalar spectral index $n_s$ with $n_s-1\simeq -2/N$, where $N$ is the number of e-foldings. In the previous work, the reconstruction of the inflaton potential for a given $n_s$ was studied, and it was found that for $n_s-1=-2/N$, the potential takes the form of either $\alpha$-attractor model or chaotic inflation model with $\phi^2$ to the leading order in the slow-roll approximation. Here we consider the reconstruction of $f(R)$ gravity model for a given $n_s$ both in the Einstein frame and in the Jordan frame. We find that for $n_s-1=-2/N$ (or more general $n_s-1=-p/N$), $f(R)$ is given parametrically and is found to asymptote to $R^2$ for large $R$. This behavior is generic as long as the scalar potential is of slow-roll type.

We estimate the spin distribution of primordial black holes based on the recent study of the critical phenomena in the gravitational collapse of a rotating radiation fluid. We find that primordial black holes are mostly slowly rotating.

We present cosmological constraints on the scalar-tensor theory of gravity by analyzing the angular power spectrum data of the cosmic microwave background obtained from the Planck 2015 results together with the baryon acoustic oscillations (BAO) data. We find that the inclusion of the BAO data improves the constraints on the time variation of the effective gravitational constant by more than $10\%$, that is, the time variation of the effective gravitational constant between the recombination and the present epochs is constrained as $G_{\rm rec}/G_0-1 <1.9\times 10^{-3}\ (95.45\%\ {\rm C.L.})$ and $G_{\rm rec}/G_0-1 <5.5\times 10^{-3}\ (99.99 \%\ {\rm C.L.})$. We also discuss the dependence of the constraints on the choice of the prior.

Cosmological constraints on the scalar-tensor theory of gravity by analyzing the angular power spectrum data of the cosmic microwave background (CMB) obtained from the Planck 2015 results are presented. We consider the harmonic attractor model, in which the scalar field has a harmonic potential with curvature ($\beta$) in the Einstein frame and the theory relaxes toward the Einstein gravity with time. Analyzing the {\it TT}, {\it EE}, {\it TE} and lensing CMB data from Planck by the Markov chain Monte Carlo method, we find that the present-day deviation from the Einstein gravity (${\alpha_0}^2$) is constrained as ${\alpha_0}^2<2.5\times10^{-4-4.5\beta^2}\ (95.45\% {\rm\ C.L.})$ and ${\alpha_0}^2<6.3\times10^{-4-4.5\beta^2}\ (99.99\%\ {\rm C.L.})$ for $0<\beta<0.4$. The time variation of the effective gravitational constant between the recombination and the present epochs is constrained as $G_{\rm rec}/G_0<1.0056\ (95.45\% {\rm\ C.L.})$ and $G_{\rm rec}/G_0<1.0115\ (99.99 \%{\rm\ C.L.})$. We also find that the constraints are little affected by extending to nonflat cosmological models because the diffusion damping effect revealed by Planck breaks the degeneracy of the projection effect.

Recent cosmological observations are in good agreement with the scalar spectral index $n_s$ with $n_s-1\sim -2/N$, where $N$ is the number of e-foldings. Quadratic chaotic model, Starobinsky model and Higgs inflation or $\alpha$-attractors connecting them are typical examples predicting such a relation. We consider the problem in the opposite: given $n_s$ as a function of $N$, what is the inflaton potential $V(\phi)$. We find that for $n_s-1=-2/N$, $V(\phi)$ is either $\tanh^2(\gamma\phi/2)$ ("T-model") or $\phi^2$ (chaotic inflation) to the leading order in the slow-roll approximation. $\gamma$ is the ratio of $1/V$ at $N\rightarrow \infty$ to the slope of $1/V$ at a finite $N$ and is related to "$\alpha$" in the $\alpha$-attractors by $\gamma^2=2/3\alpha$. The tensor-to-scalar ratio $r$ is $r=8/N(\gamma^2 N +1) $. The implications for the reheating temperature are also discussed. We also derive formulas for $n_s-1=-p/N$. We find that if the potential is bounded from above, only $p>1$ is allowed. Although $r$ depends on a parameter, the running of the spectral index is independent of it, which can be used as a consistency check of the assumed relation of $n_s(N)$.

We derive the consistency relations for a chaotic inflation model with a non-minimal coupling to gravity. For a quadratic potential in the limit of a small non-minimal coupling parameter $\xi$ and for a quartic potential without assuming small $\xi$, we give the consistency relations among the spectral index $n_s$, the tensor-to-scalar ratio $r$ and the running of the spectral index $\alpha$. We find that unlike $r$, $\alpha$ is less sensitive to $\xi$. If $r<0.1$, then $\xi$ is constrained to $\xi<0$ and $\alpha$ is predicted to be $\alpha\simeq -8\times 10^{-4}$ for a quartic potential. For a general monomial potential, $\alpha$ is constrained in the range $-2.7\times 10^{-3}<\alpha< -6\times 10^{-4}$ as long as $|\xi|\leq 10^{-3}$ if $r<0.1$.

We study the motion of a charged particle around a weakly magnetized rotating black hole. We classify the fate of a charged particle kicked out from the innermost stable circular orbit. We find that the final fate of the charged particle depends mostly on the energy of the particle and the radius of the orbit. The energy and the radius in turn depend on the initial velocity, the black hole spin, and the magnitude of the magnetic field. We also find possible evidence for the existence of bound motion in the vicinity of the equatorial plane.

The prospects for direct measurements of inflationary gravitational waves by next generation interferometric detectors inferred from the possible detection of B-mode polarization of the cosmic microwave background are studied. We compute the spectra of the gravitational wave background and the signal-to-noise ratios by two interferometric detectors (DECIGO and BBO) for large-field inflationary models in which the tensor-to-scalar ratio is greater than the order of 0.01. If the reheating temperature $T_{\rm RH}$ of chaotic inflation with the quadratic potential is high ($T_{\rm RH}>7.9\times10^6$ GeV for upgraded DECIGO and $T_{\rm RH}> 1.8\times 10^{6}$ GeV for BBO), it will be possible to reach the sensitivity of the gravitational background in future experiments at $3\sigma$ confidence level. The direct detection is also possible for natural inflation with the potential $V(\phi)=\Lambda^4 [1-\cos(\phi/f)]$, provided that $f>4.2 M_{\rm pl}$ (upgraded DECIGO) and $f>3.6 M_{\rm pl}$ (BBO) with $T_{\rm RH}$ higher than $10^8$ GeV. The quartic potential $V(\phi)=\lambda \phi^4/4$ with a non-minimal coupling $\xi$ between the inflaton field $\phi$ and the Ricci scalar $R$ gives rise to a detectable level of gravitational waves for $|\xi|$ smaller than the order of 0.01, irrespective of the reheating temperature.

Consistency relations for chaotic inflation with a monomial potential and natural inflation and hilltop inflation are given which involve the scalar spectral index $n_s$, the tensor-to-scalar ratio $r$ and the running of the spectral index $\alpha$. The measurement of $\alpha$ with $O(10^{-3})$ and the improvement in the measurement of $n_s$ could discriminate monomial model from natural/hilltop inflation models. A consistency region for general large field models is also presented.

The general k-essence Lagrangian for the existence of cosmological scaling solutions is derived in the presence of multiple scalar fields coupled to a barotropic perfect fluid. In addition to the scaling fixed point associated with the dynamics during the radiation and matter eras, we also obtain a scalar-field dominated solution relevant to dark energy and discuss the stability of them in the two-field set-up. We apply our general results to a model of two canonical fields with coupled exponential potentials arising in string theory. Depending on model parameters and initial conditions, we show that the scaling matter-dominated epochs followed by an attractor with cosmic acceleration can be realized with/without the couplings to scalar fields. The different types of scaling solutions can be distinguished from each other by the evolution of the dark energy equation of state from high-redshifts to today.

We provide the correspondence between the variables in the Jordan frame and those in the Einstein frame in scalar-tensor gravity and consider the frame-(in)dependence of the cosmological observables. In particular, we show that the cosmological observables/relations (redshift, luminosity distance, temperature anisotropies) are frame-independent. We also study the frame-dependence of curvature perturbations and find that the curvature perturbations are conformal invariant if the perturbation is adiabatic and the entropy perturbation between matter and the Brans-Dicke scalar is vanishing. The relation among various definitions of curvature perturbations in the both frames is also discussed, and the condition for the equivalence is clarified.

For two types of quintessence models having thawing and tracking properties, there exist analytic solutions for the dark energy equation of state w expressed in terms of several free parameters. We put observational bounds on the parameters in such scenarios by using the recent data of Supernovae type Ia, Cosmic Microwave Background, and Baryon Acoustic Oscillations. The observational constraints are quite different depending on whether or not the recent BAO data from BOSS are taken into account. With the BOSS data the upper bounds of today's values of w (=w_0) in thawing models is very close to -1, whereas without this data the values of w_0 away from -1 can be still allowed. The tracker equation of state w_{(0)} during the matter era is constrained to be w_{(0)}<-0.949 at 95 % confidence level even without the BOSS data, so that the tracker models with w away from -1 are severely disfavored. We also study observational constraints on scaling models in which w starts to evolve from 0 in the deep matter era and show that the transition to the equation of state close to w=-1 needs to occur at an early cosmological epoch. In the three classes of quintessence models studied in this paper, the past evolution of the Hubble parameters in the best-fit models shows only less than the 2.5 % difference compared to the LCDM.

Shadows of multi-black holes have structures distinct from the mere superposition of the shadow of a single black hole: the eyebrow-like structures outside the main shadows and the deformation of the shadows. We present analytic estimates of these structures using the static multi-black hole solution (Majumdar-Papapetrou solution). We show that the width of the eyebrow is related with the distance between the black holes and that the shadows are deformed into ellipses due to the presence of the second black holes. These results are helpful to understand qualitatively the features of the shadows of colliding black holes. We also present the shadows of colliding/coalescing black holes in the Kastor-Traschen solution.

The current observational and experimental bounds on the time variation of the constants of nature (the fine structure constant $\alpha$, the gravitational constant $G$ and the proton-electron mass ratio $\mu=m_p/m_e$) are reviewed.

We compute the shadows of colliding black holes using the Kastor-Traschen cosmological multiblack hole solution that is an exact solution describing the collision of maximally charged black holes with a positive cosmological constant. We find that in addition to the shadow of each black hole, an eyebrowlike structure appears as the black holes come close to each other. These features can be used as probes to find the multiblack hole system at the final stage of its merger process.

We present HST and UKIRT spectra and images of the 2 kpc binary quasar LBQS 0103-2753 (z=0.858). The HST images (V- and I-band) show tidal features demonstrating that this system is a major galaxy merger in progress. A two-color composite image brings out knots of star formation along the tidal arc and elsewhere. The infrared spectrum shows that both objects are at the same redshift, and that the discrepant redshift of C IV in component A is a consequence of the BAL absorption in the spectrum of this component. LBQS 0103-2753 is one of the most closely spaced binary QSOs known, and is one of relatively few dual AGN showing confirmed broad emission lines from both components. While statistical studies of binary QSOs suggest that simultaneous fueling of both black holes during a merger may be relatively rare, LBQS 0103-2753 demonstrates that such fueling can occur at high luminosity at a late stage in the merger at nuclear spacing of only a few kpc, without severe obscuration of the nuclei.

Recently spatial as well as temporal variations of the fine structure constant alpha have been reported. We show that a "runaway domain wall", which arises for the scalar field potential without minima, can account for such variations simultaneously. The time variation is induced by a runaway potential and the spatial variation is induced by the formation of a domain wall. The model is consistent with the current cosmological data and can be tested by the future experiments to test the equivalence principle.

We study the potential impact of detecting the inflationary gravitational wave background by the future space-based gravitational wave detectors, such as DECIGO and BBO. The signal-to-noise ratio of each experiment is calculated for chaotic/natural/hybrid inflation models by using the precise predictions of the gravitational wave spectrum based on numerical calculations. We investigate the dependence of each inflation model on the reheating temperature which influences the amplitude and shape of the spectrum, and find that the gravitational waves could be detected for chaotic/natural inflation models with high reheating temperature. From the detection of the gravitational waves, a lower bound on the reheating temperature could be obtained. The implications of this lower bound on the reheating temperature for particle physics are also discussed.

We study the detectability of the gravitational waves (GWs) from the Q-ball formation associated with the Affleck-Dine (AD) mechanism, taking into account both of the dilution effect due to Q-ball domination and of finite temperature effects. The AD mechanism predicts the formation of non-topological solitons, Q-balls, from which GWs are generated. Q-balls with large conserved charge $Q$ can produce a large amount of GWs. On the other hand, the decay rate of such Q-balls is so small that they may dominate the energy density of the universe, which implies that GWs are significantly diluted and that their frequencies are redshifted during Q-ball dominated era. Thus, the detectability of the GWs associated with the formation of Q-balls is determined by these two competing effects. We find that there is a finite but small parameter region where such GWs may be detected by future detectors such as DECIGO or BBO, only in the case when the thermal logarithmic potential dominates the potential of the AD field. Otherwise GWs from Q-balls would not be detectable even by these futuristic detectors: $\Omega_{\rm GW}^0<10^{-21}$. Unfortunately, for such parameter region the present baryon asymmetry of the universe can hardly be explained unless one fine-tunes A-terms in the potential. However the detection of such a GW background may give us an information about the early universe, for example, it may suggest that the flat directions with $B-L=0$ are favored.

We study time evolution of the $Q$ ball in thermal logarithmic potential using lattice simulations. As the temperature decreases due to the cosmic expansion, the thermal logarithmic term in the potential is eventually overcome by a mass term, and we confirm that the $Q$ ball transforms from the thick-wall type to the thin-wall type for a positive coefficient of radiative corrections to the mass term, as recently suggested. Moreover, we find that the $Q$ ball finally "melts down" when the $Q$-ball solution disappears. We also discuss the effects of this phenomenon on the detectability of gravitational waves from the $Q$-ball formation.

We derive the slow-roll conditions for a non-minimally coupled scalar field (extended quintessence) during the radiation/matter dominated era extending our previous results for thawing quintessence. We find that the ratio $\ddot\phi/3H\dot\phi$ becomes constant but negative, in sharp contrast to the ratio for the minimally coupled scalar field. We also find that the functional form of the equation of state of the scalar field asymptotically approaches that of the minimally coupled thawing quintessence.

We derive the equation of state of tracker fields, which are typical examples of freezing quintessence (quintessence with the equation of state approaching toward -1), taking into account of the late-time departure from the tracker solution due to the nonzero density parameter of dark energy $\Omp$. We calculate the equation of state as a function of $\Omp$ for constant $\Gamma=VV"/(V')^2$ (during matter era) models. The derived equation of state contains a single parameter, $w_{(0)}$, which parametrizes the equation of state during the matter-dominated epoch. We derive observational constraints on $w_{(0)}$ and find that observational data are consistent with the cosmological constant: $-1.11< \wzero< -0.96 (1 \sigma)$.

Important observables to reveal the nature of dark energy are the equation of state $w$ and its time derivative in units of the Hubble time $w'$. Recently, it is shown that the simplest scalar field models of dark energy (quintessence) occupy rather narrow regions in the $w-w'$ plane. We extend the $w-w'$ plane to $w<-1$ and derive bounds on $w'$ as a function of $w$ for tracker phantom dark energy. We also derive bounds on tracker k-essence.

We derive slow-roll conditions for thawing k-essence with a separable Lagrangian $p(X,\phi)=F(X)V(\phi)$. We examine the evolution of the equation of state parameter, $w$, as a function of the scale factor $a$, for the case where $w$ is close to -1. We find two distinct cases, corresponding to $X \approx 0$ and $F_X \approx 0$, respectively. For the case where $X\approx0$ the evolution of $\phi$ and hence $w$ is described by only two parameters, and $w(a)$ is model-independent and coincides with similar behavior seen in thawing quintessence models. This result also extends to non-separable Lagrangians where $X\approx0$. For the case $F_X \approx 0$, an expression is derived for $w(a)$, but this expression depends on the potential $V(\phi)$, so there is no model-independent limiting behavior. For the $X \approx 0$ case, we derive observational constraints on the two parameters of the model, $w_0$ (the present-day value of $w$), and the $K$, which parametrizes the curvature of the potential. We find that the observations sharply constrain $w_0$ to be close to -1, but provide very poor constraints on $K$.

We derive slow-roll conditions for thawing quintessence. We solve the equation of motion of $\phi$ for a Taylor expanded potential (up to the quadratic order) in the limit where the equation of state $w$ is close to -1 to derive the equation of state as a function of the scale factor. We find that the evolution of $\phi$ and hence $w$ are described by only two parameters. The expression for $w(a)$, which can be applied to general thawing models, coincides precisely with that derived recently by Dutta and Scherrer for hilltop quintessence. The consistency conditions of $|w+1|\ll 1$ are derived. The slow-roll conditions for freezing quintessence are also derived.

The spectrum of the gravitational wave background originating from quantum fluctuations during inflation is calculated numerically for various inflation models over a wide range of frequencies. We take into account four ingredients : the scalar field dynamics during inflation making no use of the slow-roll approximation, the fermionic decay of the scalar field with a small coupling constant during the reheating process, the change of the effective number of degrees of freedom g_* in the radiation-dominated era, and the anisotropic stress of free-streaming neutrinos. By numerically solving the evolution of gravitational waves during and after inflation up to the present, all of these effects can be examined comprehensively and accurately over a broad spectrum, even at very high frequencies. We find that the spectrum shows (i) a large deviation from the spectrum less accurate obtained by Taylor expanding around the CMB scale using the slow-roll approximation (ii) a characteristic frequency dependence due to the reheating effect, and (iii) damping due to the g_* changes and the neutrino anisotropic stress. We suggest that future analysis of the gravitational wave background should take into consideration the fact that analytical estimates using the Taylor expansion overestimate the amplitude of the spectrum.

As an extension of our previous study, we derive slow-roll conditions for multiple scalar fields which are non-minimally coupled with gravity and for generalized gravity theories of the form $f(\phi,R)$. We provide simple formulae of the spectral indices of scalar/tensor perturbations in terms of the slow-roll parameters.

We derive slow-roll conditions for a scalar field which is non-minimally coupled with gravity in a consistent manner and express spectral indices of scalar/tensor perturbations in terms of the slow-roll parameters. The conformal invariance of the curvature perturbation is proved without linear approximations. Rapid-roll conditions are also derived, and the relation with the slow-roll conditions is discussed.

Recently, a model of inflation using non-minimally coupled massive vector fields has been proposed. For a particular choice of non-minimal coupling parameter and for a flat FRW model, the model is reduced to the model of chaotic inflation with massive scalar field. We study the effect of non-zero curvature of the universe on the onset of vector inflation. We find that in a curved universe the dynamics of vector inflation can be different from chaotic inflation, and the fraction of the initial conditions leading to inflationary solutions is reduced compared with the chaotic inflation case.

We study weak gravitational lensing of galaxy clusters in terms of the MOND (MOdified Newtonian Dynamics) theory. We calculate shears and convergences of background galaxies for three clusters (A1689, CL0024+1654, CL1358+6245) and the mean profile of 42 SDSS (Sloan Digital Sky Survey) clusters and compare them with observational data. The mass profile is modeled as a sum of X-ray gas, galaxies and dark halo. For the shear as a function of the angular radius, MOND predicts a shallower slope than the data irrespective of the critical acceleration parameter $g_0$. The dark halo is necessary to explain the data for any $g_0$ and for three interpolation functions. If the dark halo is composed of massive neutrinos, its mass should be heavier than 2 eV. However the constraint still depends on the dark halo model and there are systematic uncertainties, and hence the more careful study is necessary to put a stringent constraint.

We provide a consistency relation between cosmological observables in general relativity with the cosmological constant. Breaking of this relation at any redshift would imply the breakdown of the hypothesis of the cosmological constant as an explanation of the current acceleration of the universe.

We calculate the signal-to-noise ratio (SNR) of the stochastic gravitational-wave background in an extreme case that its spectrum has a sharp falloff with its amplitude close to the detection threshold. Such a spectral feature is a characteristic imprint of the change in the number of relativistic degrees of freedom on the stochastic background generated during inflation in the early Universe. We find that, although SNR is maximal with the correct template which is proportional to the assumed real spectrum, its sensitivity to the shape of template is fairly weak indicating that a simple power-law template is sufficient to detect the signature.

We provide a consistency relation between cosmological observables in general relativity without relying on the equation of state of dark energy. The consistency relation should be satisfied if general relativity is the correct theory of gravity and dark energy clustering is negligible. As an extension, we also provide the DGP counterpart of the relation.

It has been proposed that cosmic acceleration or inflation can be driven by replacing the Einstein-Hilbert action of general relativity with a function f(R) of the Ricci scalar R. Such f(R) gravity theories have been shown to be equivalent to scalar-tensor theories of gravity that are incompatible with Solar System tests of general relativity, as long as the scalar field propagates over Solar System scales. Specifically, the PPN parameter in the equivalent scalar-tensor theory is gamma=1/2, which is far outside the range allowed by observations. In response to a flurry of papers that questioned the equivalence of f(R) theory to scalar-tensor theories, it was recently shown explicitly, without resorting to the scalar-tensor equivalence, that the vacuum field equations for 1/R gravity around a spherically symmetric mass also yield gamma= 1/2. Here we generalize this analysis to f(R) gravity and enumerate the conditions that, when satisfied by the function f(R), lead to the prediction that gamma=1/2.

We revisit the three-body problem in the framework of general relativity. The Newtonian N-body problem admits choreographic solutions, where a solution is called choreographic if every massive particles move periodically in a single closed orbit. One is a stable figure-eight orbit for a three-body system, which was found first by Moore (1993) and re-discovered with its existence proof by Chenciner and Montgomery (2000). In general relativity, however, the periastron shift prohibits a binary system from orbiting in a single closed curve. Therefore, it is unclear whether general relativistic effects admit a choreographic solution such as the figure eight. We carefully examine general relativistic corrections to initial conditions so that an orbit for a three-body system can be closed and a figure eight. This solution is still choreographic. This illustration suggests that the general relativistic N-body problem also may admit a certain class of choreographic solutions.

Recent astrophysical observations indicate that the proton-electron mass ratio and the fine structure constant have gone through nontrivial time evolution. We discuss their time variation in the context of a dilaton runaway scenario with gauge coupling unification at the string scale $M_{\rm s}$. We show that the choice of adjustable parameters allows them to fit the same order magnitude of both variations and their (opposite) signs in such a scenario.

None of N-body gravitating systems have been considered to emit periodic gravitational waves because of their chaotic orbits when N=3 (or more). We employ a figure-eight orbit as a specific model for a 3-body system in order to illustrate that some of triple stars are capable of generating periodic waves. This illustration would imply that a certain class of N-body gravitating systems may be relevant to the gravitational waves generation. We show also that the total angular momentum of this 3-body system is not carried away by gravitational waves. A waveform generated by this system is volcano-shaped and thus different from that of a binary system. Finally, by evaluating the radiation reaction time scale, we give an order-of-magnitude estimate of merging event rates. The estimate suggests that figure-eight sources, which require carefully prepared initial states, may be too rare to detect.

We classify the future of the universe for general cosmological models including matter and dark energy. If the equation of state of dark energy is less then -1, the age of the universe becomes finite. We compute the rest of the age of the universe for such universe models. The behaviour of the future growth of matter density perturbation is also studied. We find that the collapse of spherical overdensity region is greatly changed if the equation of state of dark energy is less than -1.

We find numerical solutions of Einstein equations and scalar field equation for a global defect in higher dimensional spacetimes ($\geq 6$). We examine in detail the relation among the expansion rate $H$ and the symmetry-breaking scale $\eta$ and the number of extra dimensions $n$ for these solutions. We find that even if the extra dimensions do not have a cigar geometry, the expansion rate $H$ grows as $\eta$ increases, which is opposite to what is needed for the recently proposed mechanism for solving the cosmological constant problem. We also find that the expansion rate $H$ decreases as $n$ increases.

We show that generalized gravity theories involving the curvature invariants of the Ricci tensor and the Riemann tensor as well as the Ricci scalar are equivalent to multi- scalar-tensor gravities with four derivatives terms. By expanding the action around a vacuum spacetime, the action is reduced to that of the Einstein gravity with four derivative terms, and consequently there appears a massive spin-2 ghost in such generalized gravity theories in addition to a massive spin-0 field.

Positronium (the bound state of electron and positron) has been thought to be formed after proton decay ($>10^{34}$yr) through collisional recombination and then decays by pair annihilation, thereby changing the matter content of the universe. We revisit the issue of the formation of positronium in the long-term future of the universe in light of recent indication that the universe is dominated by dark energy and dark matter. We find that if the equation of state of dark energy $w$ is less than -1/3 (including the cosmological constant $w=-1$), then the formation of positronium would not be possible, while it is possible through bound-bound transitions for $-1/3\siml w\siml-0.2$, or through collisional recombination for $w\simg-0.2$. The radiation from $\epm$ pair annihilation cannot dominate over $\epm$, while that from proton decay will dominate over baryon and $\epm$ for a while but not over dark matter.

We investigate the dependence of the gravitational wave spectrum from quintessential inflation on the reheating process. We consider two extreme reheating processes. One is the gravitational reheating by particle creation in the expanding universe in which the beginning of the radiation dominated epoch is delayed due to the presence of the epoch of domination of the kinetic energy of the inflaton (kination). The other is the instant preheating considered by Felder et al. in which the Universe becomes radiation dominated soon after the end of inflation. We find that the spectrum of the gravitational waves at $\sim 100$ MHz is quite sensitive to the reheating process. This result is not limited to quintessential inflation but applicable to various inflation models. Conversely, the detection or non-detection of primordial gravitational waves at $\sim$100 MHz would provide useful information regarding the reheating process in inflation.

We present a new mechanism of spontaneous baryogenesis. Usually such mechanisms require a derivative coupling between a scalar field and baryon current. In our model, the scalar field instead derivatively couples to a current associated with some global symmetry $U(1)_Q$. In this case, baryogenesis is still possible provided that an interaction exists, which violates the baryon and $U(1)_Q$ symmetries simultaneously. As a concrete example, we discuss baryogenesis in a flat direction with neither baryon nor lepton charge.

We present observational constraints on a scalar-tensor gravity theory by $\chi^2$ test for CMB anisotropy spectrum. We compare the WMAP temperature power spectrum with the harmonic attractor model, in which the scalar field has its harmonic effective potential with curvature $\beta$ in the Einstein conformal frame and the theory relaxes toward Einstein gravity with time. We found that the present value of the scalar coupling, i.e. the present level of deviation from Einstein gravity $(\alpha_0^2)$, is bounded to be smaller than $5\times 10^{-4-7\beta}$ ($2\sigma$), and $10^{-2-7\beta}$ ($4\sigma$) for $0< \beta<0.45$. This constraint is much stronger than the bound from the solar system experiments for large $\beta$ models, i.e., $\beta> 0.2$ and 0.3 in $2\sigma$ and $4\sigma$ limits, respectively. Furthermore, within the framework of this model, the variation of the gravitational constant at the recombination epoch is constrained as $|G(z=z_{rec})-G_0|/G_0 < 0.05(2\sigma)$, and $0.23(4\sigma)$.

We study the dependence of the peak luminosity of Type Ia supernovae on the fine structure constant $\alpha$. We find that decreasing (increasing) $\alpha$ enhances (reduces) the luminosity. Future experiments like SNAP could determine the variation of $\alpha$ to a precision of $10^{-2}$.

We point out that extended gravity theories, the Lagrangian of which is an arbitrary function of scalar curvature $R$, are equivalent to a class of the scalar tensor theories of gravity. The corresponding gravity theory is $\omega=0$ Brans-Dicke gravity with a potential for the Brans-Dicke scalar field, which is not compatible with solar system experiments if the field is very light: the case when such modifications are important recently.

We estimate the sensitivity of lensing observables to the parameters in the lens model (isothermal sphere/Navarro-Frenk-White profile) parameters and to cosmological parameters. We find that the observables are primarily dependent on the lens model parameters, while the dependence on cosmological parameters is minor (especially so for the dark energy parameters). We demonstrate the lens model degeneracy by deriving both the projected mass density profile and the circular velocity profile. We also identify a possible source of the problem of fitting the averaged mass profile of CL 0024+1654 with the Navarro-Frenk-White profile.

Evolution of primordial fluctuations in a Brans-Dicke type scalar-tensor gravity theory is comprehensively investigated. The harmonic attractor model, in which the scalar field has its harmonic effective potential in the Einstein conformal frame and the theory relaxes toward Einstein gravity with time, is considered. The evolution of adiabatic initial perturbations in flat SCDM models is examined from the radiation-dominated epoch up to the present. We discuss how the scalar-tensor gravity affects the evolution of metric and matter perturbations, mainly focusing on the observational consequences, i.e., the matter power spectrum and the power spectrum of cosmic microwave background temperature. We find that the early time deviation is characterized only by the large static gravitational constant while the late time behavior is qualitatively different from that in Einstein gravity because the time variation of the gravitational constant and its fluctuation have non-negligible effects. The attracting scalar-tensor gravity affects only small scale modes due to its attracting nature, the degree of which is far beyond the post-Newtonian deviation at the present epoch.

We derive a condition for converging a common evolutionary track for k-essence (a scalar field dark energy with non-canonical kinetic terms). For the Lagrangian density V(phi)W(X) with X=dot{phi}^2/2, we find tracker solutions with w_{phi} < w_B exist if Gamma equivalent V''V/(V')^2 > 3/2. Here w_{phi}(w_B) is the equation-of-state of the scalar field (background radiation/matter). Our condition may be useful for examining the existence of the attractor-like behavior in cosmology with k-essence.

The current observational and experimental bounds on time variation of the constants of Nature are briefly reviewed.

If the reported measurements of the time variation of the fine structure constant from observations of distant QSOs are correct, combined with the Oklo limit they would strongly constrain the class of the quintessence potential. If these results prove valid, future satellite experiment (STEP) should measure the induced violation of the weak equivalence principle. Future cosmological observations of nearby $(z \siml 0.5)$ absorption systems would make it clear whether the variation is significant or not.

We investigate the influence of the inner profile of lens objects on gravitational lens statistics taking into account of the effect of magnification bias and both the evolution and the scatter of halo profiles. We take the dark halos as the lens objects and consider the following three models for the density profile of dark halos; SIS (singular isothermal sphere), the NFW (Navarro Frenk White) profile, and the generalized NFW profile which has a different slope at smaller radii. The mass function of dark halos is assumed to be given by the Press-Schechter function. We find that magnification bias for the NFW profile is order of magnitude larger than that for SIS. We estimate the sensitivity of the lensing probability of distant sources to the inner profile of lenses and to the cosmological parameters. It turns out that the lensing probability is strongly dependent on the inner density profile as well as on the cosmological constant. We compare the predictions with the largest observational sample, the Cosmic Lens All-Sky Survey. The absence or presence of large splitting events in larger surveys currently underway such as the 2dF and SDSS could set constraints on the inner density profile of dark halos.

Various astronomical observations point towards the evidence for dark energy. One of the most mysterious problem is the coincidence problem: why dark energy becomes dominant only recently. We present a scenario based on extended quintessence models to explain the late-time domination of dark energy without severe fine-tuning of initial conditions and model parameters.

We investigate the feasibility of the method for reconstructing the equation of state and the effective potential of the quintessence field from SNIa data. We introduce a useful functional form to fit the luminosity distance with good accuracy (the relative error is less than 0.1%). We assess the ambiguity in reconstructing the equation of state and the effective potential which originates from the uncertainty in $\Omega_M$. We find that the equation of state is sensitive to the assumed $\Omega_M$, while the shape of the effective potential is not. We also demonstrate the actual reconstruction procedure using the data created by Monte-Carlo simulation. Future high precision measurements of distances to thousands of SNIa could reveal the shape of the quintessential potential.

To determine the equation of state of the universe, we propose to use a new independent variable $R\equiv (H_0/c)(d_L(z)/(1+z))$, where $H_0$ and $d_L(z)$ are the present Hubble parameter and the luminosity distance, respectively. For the flat universe suggested from the observation of the anisotropy of cosmic microwave background, the density and the pressure are expressed as $\rho/\rho_0=4(df/dR)^2/f^6$ and $p/\rho_0=-4/3(d^2f/dR^2)/f^5$ where $\rho_0$ is the present density and $f(R)=1/\sqrt{1+z(R)}$. In $(R, f)$ plane the sign as well as the strength of the pressure is in proportion to the curvature of the curve $f(R)$. We propose to adopt a Pade-like expression of $f(R)=1/\sqrt{u}$ with $u\equiv 1+\sum\limits_{n=1}^{N}u_nR^n$. For flat $\Lambda$ model the expansion up to N=7 has at most an error $< 0.2%$ for $z < 1.7$ and any value of $\Lambda$. We also propose a general method to determine the equation of state of the universe which has $N-1$ free parameters. If the number of parameters are smaller than $N-1$, there is a consistency check of the equation of state so that we may confirm or refute each model.

Recently, a novel class of models for inflation has been found in which the inflationary dynamics is driven solely by (non-canonical) kinetic terms rather than by potential terms. As an obvious extension, we show that a scalar field with non-canonical kinetic terms alone behaves like an energy component which is time-varying and has negative pressure presently, i.e. quintessence. We present a model which has a constant equation of state, that is, a ``kinetic'' counterpart of the Ratra-Peebles model of a quintessence field with a potential term. We make clear the structure of the phase plane and show that the quintessential solution is a late-time attractor. We also give a model for the ``phantom'' component which has an equation of state with $w=p/\rho <-1$.

We present a new type of one-field models for open inflation utilizing a nonminimally coupled scalar field with polynomial potentials in which the Coleman-de Luccia instanton does exist and slow-roll inflation after the bubble nucleation is realized.

Dynamical vacuum energy or quintessence, a slowly varying and spatially inhomogeneous component of the energy density with negative pressure, is currently consistent with the observational data. One potential difficulty with the idea of quintessence is that couplings to ordinary matter should be strongly suppressed so as not to lead to observable time variations of the constants of nature. We further explore the possibility of an explicit coupling between the quintessence field and the curvature. Since such a scalar field gives rise to another gravity force of long range ($\simg H^{-1}_0$), the solar system experiments put a constraint on the non-minimal coupling: $|\xi| \siml 10^{-2}$.

Even if the luminosity distance as a function of redshift is obtained accurately using, for example, Type Ia supernovae, the equation of state of the Universe cannot be determined uniquely but depends on one free parameter $\Omega_{k0} ={k}/(a_0^2H_0^2)$, where $a_0$ and $H_0$ are the present scale factor and the Hubble parameter, respectively. This degeneracy might be resolved if, for example, the time variations of the redshift of quasars are measured as proposed recently by Loeb. Therefore the equation of state of the Universe (or the metric of the universe) might be determined without any theoretical assumption on the matter content of the Universe in future.

We study numerically the inhomogeneous pre-big-bang inflation in a spherically symmetric space-time. We find that a large initial inhomogeneity suppresses the onset of the pre-big-bang inflation. We also find that even if the pre-big-bang inflationary stage is realized, the initial inhomogeneities are not homogenized. Namely, during the pre-big-bang inflation ``hairs''(irregularities) do not fall, in sharp contrast to the usual (potential energy dominated) inflation where initial inhomogeneity and anisotropy are damped and thus the resulting universe is less sensitive to initial conditions.

We point out that the singular instanton of Hawking-Turok type, in which the singularity occurs due to the divergence of a massless scalar field, can be generated by Euclideanized regular $p$-brane solutions in string (or M-) theory upon compactification to four dimensions.

The second derivative of the luminosity distance with respect to the redshift is written in terms of the deceleration parameter $q_0$. We point out that the third derivative contains the information regarding the sound speed of cosmic matter as well as the curvature of the universe. We restrict physically possible parameter ranges of the coefficients. It is found that there is a relation between the coefficients in a flat universe model with matter such that $c_{s0}(1+w_{\rm x0})=0$ ($c_{s0}$ is the total sound speed of the matter component and $w_{\rm x0}$=$p_{\rm x0}/\rho_{\rm x0}$).

We study gravitational lensing statistics, matter power spectra and the angular power spectra of the cosmic microwave background (CMB) radiation in x-matter models. We adopt an equation-of-state of x-matter which can express a wide range of matter from the pressureless dust to the cosmological constant. A new ingredient in this model is the sound speed of the x-component in addition to the equation-of-state $w_0=p_{x0}/\rho_{x0}$. Except for the cosmological constant case, the perturbations of x-matter itself are considered. Our primary interest is in the effect of non-zero sound speed on the structure formation and the CMB spectra. It is found that there exist parameter ranges where x-matter models are consistent with all current observations. The x-matter generally leaves imprints in the CMB anisotropy and the matter power spectrum, which should be detectable in future observations.

We study black hole MACHO binary formation through three-body interactions in the early universe at $t\sim 10^{-5}$s. The probability distribution functions of the eccentricity and the semimajor axis of binaries as well as of the coalescence time are obtained assuming that the black holes are randomly formed in space. We confirm that the previous order-of-magnitude estimate for the binary parameters is valid within $\sim 50%$ error. We find that the coalescence rate of the black hole MACHO binaries is $\sim 5 \times 10^{-2} \times 2^{\pm 1}$ events/year/galaxy taking into consideration several possible factors which may affect this estimate. This suggests that the event rate of coalescing binary black holes will be at least several events per year within 15 Mpc. The first LIGO/VIRGO interferometers in 2001 will be able to verify whether the MACHOs are black holes or not.

Spectra of density perturbations produced during chaotic inflation are calculated, taking both adiabatic and isocurvature modes into account in a class of scalar-tensor theories of gravity in which the dilaton is metrically coupled. Comparing the predicted spectrum of the cosmic microwave background radiation anisotropies with the one observed by the COBE-DMR we calculate constraints on the parameters of these theories, which turn out to be stronger by an order-of-magnitude than those obtained from post-Newtonian experiments.

We consider the horizon problem in a homogeneous but anisotropic universe (Bianchi type I). We show that the problem cannot be solved if (1) the matter obeys the strong energy condition with the positive energy density and (2) the Einstein equations hold. The strong energy condition is violated during cosmological inflation.

Motivated by the possibility of $H_0t_0 > 1$ where $H_0$ and $t_0$ are the Hubble parameter and the age of the universe, respectively, we investigate the cosmology including x-matter. x-matter is expressed by the equation of state $p_x=w_0\rho_{x0} + c_s^2(\rho_x-\rho_{x0})$, where $p_x$, $\rho_x$ and $\rho_{x0}$ are the pressure, the density of x-matter and the density at present, respectively. $w_0$ and $c_s^2$ are functions of $\rho_x$ in general. x-matter has the most general form of the equation of state which is characterized by 1)violation of strong energy condition at present for $w_0 < -1/3$ ; 2) locally stable (i.e. $c_s^2\geq 0$); 3)causality is guaranteed ($c_s\leq 1$). Considering the case that $w_0$ and $c_s^2$ are constants, we find that there is a large parameter space of $(w_0,c_s^2,\Omega_{x0})$ in which the model universe is consistent with the age of the universe and the observations of distant Type I supernovae.

Unlike general relativity, scalar-tensor theories of gravity predict scalar gravitational waves even from a spherically symmetric gravitational collapse. We solve numerically the generation and propagation of the scalar gravitational wave from a spherically symmetric and homogeneous dust collapse under the approximation that we can neglect the back reaction of the scalar wave on the space-time, and examine how the amplitude, characteristic frequency and wave form of the observed scalar gravitational wave depend on the initial radius and mass of the dust and parameters contained in the theory. In the Brans-Dicke theory, through the observation of the scalar gravitational wave, it is possible to determine the initial radius and mass and a parameter contained in the theory. In the scalar-tensor theories, it would be possible to get the information of the first derivative of the coupling function contained in the theory because the wave form of the scalar gravitational wave greatly depends on it.

The lower limit of the mass of the first stars suggested recently may imply the formation of massive stars of mass greater than 8 solar mass irrespective of the details of the initial mass function. The production of heavy metals from the first stars will ensure a requisite for the existence of life without the anthropic principle.

Under four assumptions such that 1) Einstein's theory of gravity is correct, 2) existence of foliation by geodesic slicing, 3) the trace of the extrinsic curvature, $K$, is negative at the present time in observed region, that is, the observed universe is now expanding, 4) the strong energy condition is satisfied, we show that $H_{v0} t_{0} \leq 1$, where $H_v \equiv - K/3$ agrees with the Hubble parameter in the case of a homogeneous and isotropic universe, and $t_0$ is the age of the Universe. If $H_{v0} t_0 > 1$ is confirmed observationally, at least one of the four assumptions is incorrect.

We investigate the distance-redshift relation in the simple void model. As discussed by Moffat and Tatarski, if the observer stays at the center of the void, the observed Hubble parameter is not so different from the background Hubble parameter. However, if the position of observer is off center of the void, we must consider the peculiar velocity correction which is measured by the observed dipole anisotropy of cosmic microwave background. This peculiar velocity correction for the redshift is crucial to determine the Hubble parameter and we shall discuss this effect. Further the results of Turner et al by the N-body simulation will be also considered.

We have performed the most comprehensive predictions of the temperature fluctuations \dtt in the primeval isocurvature baryon models to see whether or not the models are consistent with the recent data on the cosmic microwave background anisotropies. More specifically, we computed the \dtt corresponding to the experimental set-up by the South-Pole and the Owens Valley experiments as well as the COBE satellite. The amplitudes of the predicted \dtt are normalized by means of the COBE 10$^\circ$ data. The resulting constraints on the models are presented on $n - \Omega_b$ plane in the case of $\lambda_0=1-\Omega_b$ (flat models) and $\lambda_0=0$ (open models), where $n$ is the primordial spectral index of entropy fluctuations and $\Omega_b$ is the present baryon density parameter. Our results imply that the PIB models cannot be reconciled with the current observations for any reasonable set of cosmological parameters.