We present the first astrometry catalog from the Japanese VLBI (very long baseline interferometer) project VERA (VLBI Exploration of Radio Astrometry). We have compiled all the astrometry results from VERA, providing accurate trigonometric annual parallax and proper motion measurements. In total, 99 maser sources are listed in the VERA catalog. Among them, 21 maser sources are newly reported while the rest of 78 sources are referred to previously published results or those in preparation for forthcoming papers. The accuracy in the VERA astrometry are revisited and compared with those from the other VLBI astrometry projects such as BeSSeL (The Bar and Spiral Structure Legacy) Survey and GOBELINS (the Gould's Belt Distances Survey) with the VLBA (Very Long Baseline Array). We have confirmed that most of the astrometry results are consistent with each other, and the largest error sources are due to source structure of the maser features and their rapid variation, along with the systematic calibration errors and different analysis methods. Combined with the BeSSeL results, we estimate the up-to-date fundamental Galactic parameter of $R_{0}=7.92\pm0.16_{\rm{stat.}}\pm0.3_{\rm{sys.}}$~kpc and $\Omega_{\odot}=30.17\pm0.27_{\rm{stat.}}\pm0.3_{\rm{sys.}}$~km~s$^{-1}$~kpc$^{-1}$, where $R_{0}$ and $\Omega_{\odot}$ are the distance from the Sun to the Galactic center and the Sun's angular velocity of the Galactic circular rotation, respectively.

We have investigated inflation models that predict a very small value of the tensor-to-scalar ratio, $r$. The spectral index $n_s$, and the tensor-to-scalar ratio $r$, are strictly constrained by the Planck data. $n_s$ and $r$ are sensitive to the shape and magnitude of the inflaton potential, respectively.The constraints by the Planck 2018 data combined with other cosmological observations are compared with the predictions from the inflation models regarding $n_s$ and $r$. Furthermore, we discuss the comparison of future tensor-to-scalar ratio data with predictions from the inflation models with a focus on part of the quantum fluctuation origin.

Many studies have shown that there are clear sequences in the period-luminosity relationship (PLR) for Mira variables and semiregular variables (SRVs) in the Large Magellanic Cloud (LMC). To investigate the PLR for SRVs in our galaxy, we examined the annual parallax measurement and conducted K'-band photometric monitoring of an SRV star SV Pegasus (SV Peg). We measured the position change of the associating H$_2$O maser spots by phase-referencing VLBI observations with VERA at 22 GHz, spanning approximately 3 yr, and detected an annual parallax of $\pi = 3.00 \pm 0.06$ mas, corresponding to a distance of $D=333 \pm 7$ pc. This result is in good agreement with the Hipparcos parallax and improves the accuracy of the distance from 35 $\%$ to 2 $\%$. However, the GAIA DR2 catalog gave a parallax of $\pi=1.12\pm0.28$ mas for SV Peg. This indicates that the GAIA result might be blurred by the effect of the stellar size because the estimated stellar radius was $\sim 5$ mas, which is comparable to the parallax. We obtained a K'-band mean magnitude of $m_{K'} = -0.48$ mag and a period of $P=177$ days from our photometric monitoring with a 1-m telescope. Using the trigonometric distance, we derived an absolute magnitude of $M_{K'}=-8.09 \pm 0.05$ mag. This result shows that the position of SV Peg in the PLR falls on the C' sequence found in the PLR in the LMC, which is similar to other SRVs in our galaxy.

A probability of distinguishing between $\Lambda$CDM model and modified gravity is studied by using the future observations for the growth rate of cosmic structure (Euclid redshift survey). Adopting extended DGP model, Kinetic Gravity Braiding model, and Galileon model as modified gravity, we compare predicted cosmic growth rate by models to the mock observational data. The growth rate $f\sigma_8$ in the original DGP model is suppressed compared with the $\Lambda$CDM case, for the same value of the current density parameter of matter $\Omega_{m,0}$, because of the suppression of effective gravitational constant. In case of the kinetic gravity braiding model and the Galileon model, the growth rate $f\sigma_8$ is enhanced compared with the $\Lambda$CDM case, for the same value of $\Omega_{m,0}$, because of the enhancement of effective gravitational constant. For future observational data of the cosmic growth rate (Euclid), compatible value of $\Omega_{m,0}$ are different by models, furthermore $\Omega_{m,0}$ can be stringently constrained. Thus, we find the $\Lambda$CDM model is distinguishable from modified gravity by combining the growth rate data of the Euclid with other observations.

Constraints on neutrino masses are estimated based on future observations of the cosmic microwave background (CMB), which includes the B-mode polarization produced by CMB lensing from the Planck satellite, and the growth rate of cosmic structure from the Euclid redshift survey by using the Markov-Chain Monte-Carlo (MCMC) method. The error in the bound on the total neutrino mass is estimated to be $\Delta\sum m_{\nu} = 0.075$ eV with a 68\% confidence level. By using the growth rate rather than the galaxy power spectrum, accurate constraints are obtained, since the growth rate is less influenced by the uncertainty regarding galaxy bias than by the galaxy power spectrum.

Constraints on neutrino masses are estimated based on future observations of the cosmic microwave background (CMB) including the B-mode polarization produced by CMB lensing using the Planck satellite, and baryon acoustic oscillations distance scale and the galaxy power spectrum from all-sky galaxy redshift survey in the BigBOSS experiment. We estimate the error in the bound on the total neutrino mass to be $\Delta\sum m_{\nu} = 0.012$ eV with a 68% confidence level. If the fiducial value of the total neutrino mass is $\sum m_{\nu}= 0.06$ eV, this result implies that the neutrino mass hierarchy must be normal.

We studied the cosmological constraints on the Galileon gravity obtained from observational data of the growth rate of matter density perturbations, the supernovae Ia (SN Ia), the cosmic microwave background (CMB), and baryon acoustic oscillations (BAO). For the same value of the energy density parameter of matter $\Omega_{m,0}$, the growth rate $f$ in Galileon models is enhanced, relative to the $\Lambda$CDM case, because of an increase in Newton's constant. The smaller $\Omega_{m,0}$ is, the more growth rate is suppressed. Therefore, the best fit value of $\Omega_{m,0}$ in the Galileon model, based only the growth rate data, is quite small. This is incompatible with the value of $\Omega_{m,0}$ obtained from the combination of SN Ia, CMB, and BAO data. On the other hand, in the $\Lambda$CDM model, the values of $\Omega_{m,0}$ obtained from different observational data sets are consistent. In the analysis of this paper, we found that the Galileon model is less compatible with observations than the $\Lambda$CDM model. This result seems to be qualitatively the same in most of the generalized Galileon models in which Newton's constant is enhanced.

We have succeeded in establishing a cosmological model with a non-minimally coupled scalar field $\phi$ that can account not only for the spatial periodicity or the picket-fence structure exhibited by the galaxy $N$-$z$ relation of the 2dF survey, but also for the spatial power spectrum of the cosmic microwave background radiation (CMB) temperature anisotropy observed by the WMAP satellite. The scalar field of our model universe starts from an extremely small value at around the nucleosynthesis epoch, remains in that state for sufficiently long periods, allowing sufficient time for the CMB temperature anisotropy to form, and then starts to grow in magnitude at the redshift $z$ of $\sim 1$, followed by a damping oscillation which is required to reproduce the observed picket-fence structure of the $N$-$z$ relation. To realize such behavior of the scalar field, we have found it necessary to introduce a new form of potential $V(\phi)\propto \phi^2\exp(-q\phi^2)$, with $q$ being a constant. Through this parameter $q$, we can control the epoch at which the scalar field starts growing.

We compare observational data of growth rate with the prediction by Galileon theory. For the same value of the energy density parameter $\Omega_{m,0}$, the growth rate in Galileon models is enhanced compared with the $\Lambda$CDM case, due to the enhancement of Newton's constant. The smaller $\Omega_{m,0}$ is, the more suppressed growth rate is. Hence the best fit value of $\Omega_{m,0}$ in the Galileon model is 0.16 from only the growth rate data, which is considerably smaller than such value obtained from observations of supernovae Ia, the cosmic microwave background and baryon acoustic oscillations. This result seems to be qualitatively the same in most of the generalized Galileon models. We also find the upper limit of the Brans-Dicke parameter to be $\omega < -40$, from the growth rate data. More and better growth rate data are required to distinguish between dark energy and modified gravity.

We study the observational constraints on the Phantom Crossing DGP model. We demonstrate that the crossing of the phantom divide does not occur within the framework of the original Dvali-Gabadadze-Porrati (DGP) model or the DGP model developed by Dvali and Turner. By extending their model in the framework of an extra dimension scenario, we study a model that realizes crossing of the phantom divide. We investigate the cosmological constraints obtained from the recent observational data of Type Ia Supernovae, Cosmic Microwave Background anisotropies, and Baryon Acoustic Oscillations. The best fit values of the parameters with 1$\sigma$ (68%) errors for the Phantom Crossing DGP model are $\Omega_{m,0}=0.27^{+0.02}_{-0.02}$, $\beta=0.54^{+0.24}_{-0.30}$. We find that the Phantom Crossing DGP model is more compatible with the observations than the original DGP model or the DGP model developed by Dvali and Turner. Our model can realize late-time acceleration of the universe, similar to that of $\Lambda$CDM model, without dark energy due to the effect of DGP gravity. In our model, crossing of the phantom divide occurs at a redshift of $z \sim 0.2$.

We investigate the observational constraints on the oscillating scalar field model using data from type Ia supernovae, cosmic microwave background anisotropies, and baryon acoustic oscillations. According to a Fourier analysis, the galaxy number count $N$ from redshift $z$ data indicates that galaxies have preferred periodic redshift spacings. We fix the mass of the scalar field as $m_\phi=3.2\times 10^{-31}h$ ${\rm eV}$ such that the scalar field model can account for the redshift spacings, and we constrain the other basic parameters by comparing the model with accurate observational data. We obtain the following constraints: $\Omega_{m,0}=0.28\pm 0.03$ (95% C.L.), $\Omega_{\phi,0} < 0.035$ (95% C.L.), $\xi > -158$ (95% C.L.) (in the range $\xi \le 0$). The best fit values of the energy density parameter of the scalar field and the coupling constant are $\Omega_{\phi,0}= 0.01$ and $\xi= -25$, respectively. The value of $\Omega_{\phi,0}$ is close to but not equal to $0$. Hence, in the scalar field model, the amplitude of the galaxy number count cannot be large. However, because the best fit values of $\Omega_{\phi,0}$ and $\xi$ are not $0$, the scalar field model has the possibility of accounting for the periodic structure in the $N$--$z$ relation of galaxies. The variation of the effective gravitational constant in the scalar field model is not inconsistent with the bound from observation.

We investigate whether or nor it is possible to find a scalar field model or models that are capable of explaining simultaneously the observed $N$-$z$ relation given by the 2dF Galaxy Redshift Survey, which still seems to exhibit a spatial periodicity of the galaxy distribution(the 'picket-fence structure'), and the CMB spectrum obtained by the WMAP experiments. It is found that both the observed size of the spatial periodicity and the amplitude of the 2dF $N$-$z$ relation can be fairly well fitted by the theoretical computations based on the scalar field models with $-20\le \xi\le -10$, and $140\le m_{\rm s} \le 160$, where $\xi$ is the gravitational coupling parameter, and $m_{\rm s}$ the normalized mass of the scalar field, respectively. To reproduce the CMB spectrum of the WMAP, we find that it is very crucial to have a null state of the scalar field in the early epochs of evolution of the universe.

We propose a phantom crossing Dvali--Gabadadze--Porrati (DGP) model. In our model, the effective equation of state of the DGP gravity crosses the phantom divide line. We demonstrate crossing of the phantom divide does not occur within the framework of the original DGP model or the DGP model developed by Dvali and Turner. By extending their model, we construct a model that realizes crossing of the phantom divide. We find that the smaller the value of the new introduced parameter $\beta$ is, the older epoch crossing of the phantom divide occurs in. Our model can account for late-time acceleration of the universe without dark energy. We investigate and show the property of Phantom crossing DGP model.

We propose a phantom crossing Dvali--Gabadadze--Porrati (DGP) model. In our model, the effective equation of state of the DGP gravity crosses the phantom divide line. We demonstrate crossing of the phantom divide does not occur within the framework of the original DGP model or the DGP model developed by Dvali and Turner. By extending their model, we construct a model that realizes crossing of the phantom divide. DGP models can account for late-time acceleration of the universe without dark energy. Phantom Crossing DGP model is more compatible with recent observational data from Type Ia Supernovae (SNIa), Cosmic Microwave Background (CMB) anisotropies, and Baryon Acoustic Oscillations (BAO) than the original DGP model or the DGP model developed by Dvali and Turner.

A Fourier analysis on galaxy number counts from redshift data of both the Sloan Digital Sky Survey and the 2dF Galaxy Redshift Survey indicates that galaxies have preferred periodic redshift spacings of $\Delta z =$ 0.0102, 0.0246, and 0.0448 in the SDSS and strong agreement with the results from the 2dF GRS. The redshift spacings are confirmed by the mass density fluctuations, the power spectrum $P(z)$ and $N_{pairs}$ calculations. Application of the Hubble law results in galaxies preferentially located on co-moving concentric shells with periodic spacings. The combined results from both surveys indicate regular co-moving radial distance spacings of $31.7 \pm 1.8$ $h^{-1} Mpc$, $73.4 \pm 5.8$ $h^{-1} Mpc$ and $127 \pm 21$ $h^{-1} Mpc$. The results are consistent with oscillations in the expansion rate of the universe over past epochs.

We have succeeded in establishing a cosmological model with a non-minimally coupled scalar field $\phi$ that can account not only for the spatial periodicity or the {\it picket-fence structure} exhibited by the galaxy $N$-$z$ relation of the 2dF survey but also for the spatial power spectrum of the cosmic microwave background radiation (CMB) temperature anisotropy observed by the WMAP satellite. The Hubble diagram of our model also compares well with the observation of Type Ia supernovae. The scalar field of our model universe starts from an extremely small value at around the nucleosynthesis epoch, remains in that state for sufficiently long periods, allowing sufficient time for the CMB temperature anisotropy to form, and then starts to grow in magnitude at the redshift $z$ of $\sim 1$, followed by a damping oscillation which is required to reproduce the observed picket-fence structure of the $N$-$z$ relation. To realize such behavior of the scalar field, we have found it necessary to introduce a new form of potential $V(\phi)\propto \phi^2\exp(-q\phi^2)$, with $q$ being a constant. Through this parameter $q$, we can control the epoch at which the scalar field starts growing.