The apparent anisotropies of the galaxy clustering in observable redshift space provide a unique opportunity to simultaneously probe cosmic expansion and gravity on cosmological scales via the Alcock--Paczynski effect and redshift-space distortions. While the improved theoretical models have been proposed and developed to describe the apparent anisotropic clustering at weakly non-linear scales, the applicability of these models is still limited in the presence of the non--perturbative smearing effect caused by the randomness of the relative velocities. Although the cosmological constraint from the anisotropic clustering will be certainly improved with a more elaborate theoretical model, we here consider an alternative way by using the statistical power of both the power spectrum and bispectrum at large scales. Based on the Fisher matrix analysis, we estimate the benefit of combining the power spectra and bispectra, finding that the constraints on the cosmic expansion and growth of structure will be improved by a factor of two. This compensates for the loss of constraining power using the power spectrum alone due to the randomness of the relative velocities.

Redshift-space distortions (RSD) offer an attractive method to measure the growth of cosmic structure on large scales, and combining with the measurement of the cosmic expansion history, it can be used as cosmological tests of gravity. With the advent of future galaxy redshift surveys aiming at precisely measuring the RSD, an accurate modeling of RSD going beyond linear theory is a critical issue in order to detect or disprove small deviations from general relativity (GR). While several improved models of RSD have been recently proposed based on the perturbation theory (PT), the framework of these models heavily relies on GR. Here, we put forward a new PT prescription for RSD in general modified gravity models. As a specific application, we present theoretical predictions of the redshift-space power spectra in f(R) gravity model, and compare them with N-body simulations. Using the PT template that takes into account the effects of both modifications of gravity and RSD properly, we successfully recover the fiducial model parameter in N-body simulations in an unbiased way. On the other hand, we found it difficult to detect the scale dependence of the growth rate in a model-independent way based on GR templates.

We model the apparent clustering anisotropy of Luminous Red Galaxies (LRGs) in the Sloan Digital Sky Survey using subhalos identified in cosmological $N$-body simulations. We first conduct a Markov-chain Monte Carlo analysis on the parameters characterizing subhalos hosting LRGs assuming a specific $\Lambda$CDM cosmology on which we run the simulations. We show that simple models with central and satellite subhalos can explain the observed multipole moments of the power spectrum up to hexadecapole on large scales ($k\lesssim0.3~h\mathrm{Mpc}^{-1}$). A satellite fraction of $20$ to $30$ per cent is favored weakly depending on the detail of the model. The fraction is shown to be robust when we adopt a more refined model based on the halo occupation number from the literature. We then vary cosmological parameters controlling the anisotropy in redshift-space effectively by deforming the simulation box (the Alcock-Paczynski effect) and changing the amplitude of the velocities (the redshift-space distortions). We demonstrate that we can constrain the geometry of the universe, the structure growth rate, and the parameters characterizing LRGs simultaneously. This is a step toward cosmological analysis with realistic bias description beyond empirical bias functions with nuisance parameters.

The anisotropic galaxy clustering on large scales provides us with a unique opportunity to probe into the gravity theory through the redshift-space distortions (RSDs) and the Alcock-Paczynski effect. Using the multipole power spectra up to hexadecapole (ell=4), of the Luminous Red Galaxy (LRG) sample in the data release 7 (DR7) of the Sloan Digital Sky Survey II (SDSS-II), we obtain simultaneous constraints on the linear growth rate f, angular diameter distance D_A, and Hubble parameter H at redshift z = 0.3. For this purpose, we first extensively examine the validity of a theoretical model for the non-linear RSDs using mock subhalo catalogues from N-body simulations, which are constructed to match with the observed multipole power spectra. We show that the input cosmological parameters of the simulations can be recovered well within the error bars by comparing the multipole power spectra of our theoretical model and those of the mock subhalo catalogues. We also carefully examine systematic uncertainties in our analysis by testing the dependence on prior assumption of the theoretical model and the range of wavenumbers to be used in the fitting. These investigations validate that the theoretical model can be safely applied to the real data. Thus, our results from the SDSS DR7 LRG sample are robust including systematics of theoretical modeling; f(z = 0.3) sigma_8(z = 0.3) =0.49+-0.08(stat.)+-0.04(sys.), D_A (z = 0.3) =968+-42(stat.)+-17(sys.)[Mpc], H (z = 0.3) =81.7+-5.0(stat.)+-3.7(sys.)[km/s/Mpc]. We believe that our method to constrain the cosmological parameters using subhaloes catalogues will be useful for more refined samples like CMASS and LOWZ catalogues in the Baryon Oscillation Spectroscopic Survey in SDSS-III.

We investigate the effect of photodesorption on the snow line position at the surface of a protoplanetary disk around a Herbig Ae/Be star, motivated by the detection of water ice particles at the surface of the disk around HD142527 by Honda et al. For this aim, we obtain the density and temperature structure in the disk with a 1+1D radiative transfer and determine the distribution of water ice particles in the disk by the balance between condensation, sublimation, and photodesorption. We find that photodesorption induced by the far-ultraviolet radiation from the central star depresses the ice-condensation front toward the mid-plane and pushes the surface snow line outward significantly when the stellar effective temperature exceeds a certain critical value. This critical effective temperature depends on the stellar luminosity and mass, the water abundance in the disk, and the yield of photodesorption. We present an approximate analytic formula for the critical temperature. We separate Herbig Ae/Be stars into two groups on the HR diagram according to the critical temperature; one is the disks where photodesorption is effective and from which we may not find ice particles at the surface, and the other is the disks where photodesorption is not effective. We estimate the snow line position at the surface of the disk around HD142527 to be 100--300 AU which is consistent with the water ice detection at >140 AU in the disk. All results depend on the dust grain size by a complex way and this point requires more work in future.

Evolution of a snow line in an optically-thick protoplanetary disk is investigated with numerical simulations. The ice-condensing region in the disk is obtained by calculating the temperature and the density with the 1+1D approach. The snow line migrates as the mass accretion rate (\dot{M}) in the disk decreases with time. Calculations are carried out from an early phase with high disk accretion rates (\dot{M} \sim 10^{-7} M_sun/yr) to a later phase with low disk accretion rates (\dot{M} \sim 10^{-12} M_sun/yr) using the same numerical method. It is found that the snow line moves inward for \dot{M} > 10^{-10} M_sun/yr, while it gradually moves outward in the later evolution phase with \dot{M} < 10^{-10} M_sun/yr. In addition to the silicate opacity, the ice opacity is taken into consideration. In the inward migration phase, the additional ice opacity increases the distance of the snow line from the central star by a factor of 1.3 for dust grains < 10 micro-meter in size and 1.6 for > 100 micro-meter. It is inevitable that the snow line comes inside the Earth's orbit in the course of the disk evolution, if the alpha viscosity parameter is in a range 0.001-0.1, the dust-to-gas mass ratio is higher than a tenth of the solar abundance value, and the dust grains are smaller than 1 mm. The formation of water-devoid planetesimals in the terrestrial planet region seems to be difficult throughout the disk evolution, which imposes a new challenge to planet formation theory.

The temperature in the optically thick interior of protoplanetary discs is essential for the interpretation of millimeter observations of the discs, for the vertical structure of the discs, for models of the disc evolution and the planet formation, and for the chemistry in the discs. Since large icy grains have a large albedo even in the infrared, the effect of scattering of the diffuse radiation in the discs on the interior temperature should be examined. We have performed a series of numerical radiation transfer simulations including isotropic scattering by grains with various typical sizes for the diffuse radiation as well as for the incident stellar radiation. We also have developed an analytic model including isotropic scattering to understand the physics concealed in the numerical results. With the analytic model, we have shown that the standard two-layer approach is valid only for grey opacity (i.e. grain size $\ga10$ \micron) even without scattering. A three-layer interpretation is required for grain size $\la10$ \micron. When the grain size is 0.1--10 \micron, the numerical simulations show that isotropic scattering reduces the temperature of the disc interior. This reduction is nicely explained by the analytic three-layer model as a result of the energy loss by scatterings of the incident stellar radiation and of the warm diffuse radiation in the disc atmosphere. For grain size $\ga10$ \micron (i.e. grey scattering), the numerical simulations show that isotropic scattering does not affect the interior temperature. This is nicely explained by the analytic two-layer model; the energy loss by scattering in the disc atmosphere is exactly offset by the "green-house effect" due to scattering of the cold diffuse radiation in the interior.

We discuss how we obtain the spatial distribution of ice on the surface of the circumstellar disk around young stars. Ice in the disks plays a very important role in various issues, for instance, on the disk structure, on the planet formation, on the isotopic anomaly in meteorites, and on the origin of the sea on the Earth. Therefore, the spatially resolved observation of the condensation/sublimation front of ice, so-called ``snow line'' is strongly required. Here, we propose a new method for obtaining the spatially resolved ``snow line'' on the circumstellar disks by observing 3 \micron H$_2$O ice feature in the scattered light. Based on radiative transfer considerations, we show that the feature is clearly imprinted in the spectrum of the scattered light from both optically thick and thin circumstellar disks. We also show that the scattered light and the H$_2$O ice feature from protoplanetary disks are detectable and spatially resolvable with the current instruments through a $H_2O$ narrowband filter around 3 \micron. Finally, we present a diagnostics of disk dust properties on the $K-H_2O$ and $K-L'$ two color diagram.