Showing 1–50 of 74 results for author: Futamase, T

New highly precise weak gravitational lensing flexions measurement method based on ERA method

Authors: Yuki Okura, Toshifumi Futamase

Abstract: Weak gravitational lensing flexions are a kind of weak lensing distortion which are defined as the spin 1 and spin 3 combinations of the third order derivatives of gravitational lensing potential. Since the shear has spin 2 combination of the second order derivative, the flexion signal gives a partly independent information from shear signal and is more sensitive to the local mass distribution tha… ▽ More Weak gravitational lensing flexions are a kind of weak lensing distortion which are defined as the spin 1 and spin 3 combinations of the third order derivatives of gravitational lensing potential. Since the shear has spin 2 combination of the second order derivative, the flexion signal gives a partly independent information from shear signal and is more sensitive to the local mass distribution than shear signal. Thus its measurement is expected to play important roles in observational cosmology. However, since the weakness of the flexion signal as well as the complicatedness of its intrinsic noise made its accurate observation very difficult. We propose a new method of measuring the flexion signal using ERA method which is a method to measure weak lensing shear without any approximation. We find two particular combinations of the flexions which provide the quantities with only lensing information and free of intrinsic noise when taken average. It is confirmed by simple numerical simulation that the statistical average of these combinations do not in fact depend on the strength of the intrinsic distortion. △ Less

Submitted 22 November, 2021; v1 submitted 31 August, 2021; originally announced September 2021.

Detection of Cosmic Magnification via Galaxy Shear -- Galaxy Number Density Correlation from HSC Survey Data

Authors: Xiangkun Liu, Dezi Liu, Zucheng Gao, Chengliang Wei, Guoliang Li, Liping Fu, Toshifumi Futamase, Zuhui Fan

Abstract: We propose a novel method to detect cosmic magnification signals by cross-correlating foreground convergence fields constructed from galaxy shear measurements with background galaxy positional distributions, namely shear-number density correlation. We apply it to the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) survey data. With 27 non-independent data points and their full covariance,… ▽ More We propose a novel method to detect cosmic magnification signals by cross-correlating foreground convergence fields constructed from galaxy shear measurements with background galaxy positional distributions, namely shear-number density correlation. We apply it to the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) survey data. With 27 non-independent data points and their full covariance, $χ_0^2\approx 34.1$ and $χ_T^2\approx 24.0$ with respect to the null and the cosmological model with the parameters from HSC shear correlation analyses in Hamana et al. 2020 (arXiv:1906.06041), respectively. The Bayes factor of the two is $\log_{10}B_{T0}\approx 2.2$ assuming equal model probabilities of null and HSC cosmology, showing a clear detection of the magnification signals. Theoretically, the ratio of the shear-number density and shear-shear correlations can provide a constraint on the effective multiplicative shear bias $\bar m$ using internal data themselves. We demonstrate the idea with the signals from our HSC-SSP mock simulations and rescaling the statistical uncertainties to a survey of $15000°^2$. For two-bin analyses with background galaxies brighter than $m_{lim}=23$, the combined analyses lead to a forecasted constraint of $σ(\bar m) \sim 0.032$, $2.3$ times tighter than that of using the shear-shear correlation alone. Correspondingly, $σ(S_8)$ with $S_8=σ_8(Ω_\mathrm{m}/0.3)^{0.5}$ is tightened by $\sim 2.1$ times. Importantly, the joint constraint on $\bar m$ is nearly independent of cosmological parameters. Our studies therefore point to the importance of including the shear-number density correlation in weak lensing analyses, which can provide valuable consistency tests of observational data, and thus to solidify the derived cosmological constraints. △ Less

Submitted 28 April, 2021; originally announced April 2021.

Comments: 9 pages, 4 figures; Accepted for publication in PRD

Journal ref: Published in Physical Review D, Volume 103, Issue 12, p.123504, 2021

Standard candles and sirens rescue $H_0$

Authors: Aniket Agrawal, Teppei Okumura, Toshifumi Futamase

Abstract: We show that future observations of binary neutron star systems with electromagnetic counterparts together with the traditional probes of low- and high-redshift Type Ia supernovae (SNe Ia) can help resolve the Hubble tension. The luminosity distance inferred from these probes and its scatter depend on the underlying cosmology. By using the gravitational lensing of light or gravitational waves emit… ▽ More We show that future observations of binary neutron star systems with electromagnetic counterparts together with the traditional probes of low- and high-redshift Type Ia supernovae (SNe Ia) can help resolve the Hubble tension. The luminosity distance inferred from these probes and its scatter depend on the underlying cosmology. By using the gravitational lensing of light or gravitational waves emitted by, and peculiar motion of, these systems we derive constraints on the sum of neutrino masses, the equation of state of dark energy parametrized in the form $w_0 + w_a (1-a)$, along with the Hubble constant and cold dark matter density in the universe. We show that even after marginalizing over poorly constrained physical quantities, such as the sum of neutrino masses and the nature of dark energy, low-redshift gravitational-wave observations, in combination with SNe Ia, have the potential to rule out new physics as the underlying cause of the Hubble tension at $\gtrsim 5.5σ$. △ Less

Submitted 1 December, 2020; v1 submitted 3 August, 2020; originally announced August 2020.

Comments: 6 pages, 3 figures, comments welcome; matches version published in ApJ

Journal ref: ApJ, 904, 169 (2020)

Survey of Gravitationally-lensed Objects in HSC Imaging (SuGOHI). V. Group-to-cluster scale lens search from the HSC-SSP Survey

Authors: Anton T. Jaelani, Anupreeta More, Masamune Oguri, Alessandro Sonnenfeld, Sherry H. Suyu, Cristian E. Rusu, Kenneth C. Wong, James H. H. Chan, Issha Kayo, Chien-Hsiu Lee, Dani C. -Y. Chao, Jean Coupon, Kaiki T. Inoue, Toshifumi Futamase

Abstract: We report the largest sample of candidate strong gravitational lenses belonging to the Survey of Gravitationally-lensed Objects in HSC Imaging for group-to-cluster scale (SuGOHI-c) systems. These candidates are compiled from the S18A data release of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) Survey. We visually inspect $\sim39,500$ galaxy clusters, selected from several catalogs, ove… ▽ More We report the largest sample of candidate strong gravitational lenses belonging to the Survey of Gravitationally-lensed Objects in HSC Imaging for group-to-cluster scale (SuGOHI-c) systems. These candidates are compiled from the S18A data release of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) Survey. We visually inspect $\sim39,500$ galaxy clusters, selected from several catalogs, overlapping with the Wide, Deep, and UltraDeep fields, spanning the cluster redshift range $0.05<z_{cl}<1.38$. We discover 641 candidate lens systems, of which 536 are new. From the full sample, 47 are almost certainly bonafide lenses, 181 of them are highly probable lenses and 413 are possible lens systems. Additionally, we present 131 lens candidates at galaxy-scale serendipitously discovered during the inspection. We obtained spectroscopic follow-up of 10 candidates using the X-shooter. With this follow-up, we confirm 8 systems as strong gravitational lenses. Of the remaining two, one of the sources is too faint to detect any emission, and the other has a tentative redshift close to the lens redshift, but additional arcs in this system are yet to be observed spectroscopically. Since the HSC-SSP is an ongoing survey, we expect to find $\sim600$ definite or probable lenses using this procedure and even more if combined with other lens finding methods. △ Less

Submitted 16 April, 2020; v1 submitted 4 February, 2020; originally announced February 2020.

Comments: Accepted for publication in MNRAS, 19 pages, 11 figures

Constraining neutrino mass and dark energy with peculiar velocities and lensing dispersions of Type Ia supernovae

Authors: Aniket Agrawal, Teppei Okumura, Toshifumi Futamase

Abstract: We show that peculiar velocities of Type Ia supernovae can be used to derive constraints on the sum of neutrino masses, $Σm_ν$, and dark energy equation of state, $w = w_0+w_a(1-a)$, from measurements of the magnitude-redshift relation, complementary to galaxy redshift and weak lensing surveys. Light from a supernova propagates through a perturbed Universe so the luminosity distance is modified fr… ▽ More We show that peculiar velocities of Type Ia supernovae can be used to derive constraints on the sum of neutrino masses, $Σm_ν$, and dark energy equation of state, $w = w_0+w_a(1-a)$, from measurements of the magnitude-redshift relation, complementary to galaxy redshift and weak lensing surveys. Light from a supernova propagates through a perturbed Universe so the luminosity distance is modified from its homogeneous prediction. This modification is proportional to the matter density fluctuation and its time derivative due to gravitational lensing and peculiar velocity respectively. At low redshifts, the peculiar velocity signal dominates while at high redshifts lensing does. We show that using lensing and peculiar velocity of supernovae from the upcoming surveys WFIRST and ZTF, without other observations, we can constrain $Σm_ν \lesssim 0.31$ eV, $σ(w_0) \lesssim 0.02$, and ${σ(w_a)} \lesssim 0.18$ ($1-σ$ CL) in the $Σm_ν$-$w_0$-$w_a$ parameter space, where all the other cosmological parameters are fixed. We find that adding peculiar velocity information from low redshifts shrinks the volume of the parameter ellipsoid in this space by $\sim 33$%. We also allow $Ω_{\text{CDM}}$ to vary as well as $Σm_ν$, $w_0$ and $w_a$, and demonstrate how these constraints degrade as a consequence. △ Less

Submitted 25 September, 2019; v1 submitted 4 July, 2019; originally announced July 2019.

Comments: 7+1 pages, 6 figures; comments welcome; matches published version

Journal ref: Phys. Rev. D 100, 063534 (2019)

Forecasts for warm dark matter from weakly lensed Type Ia supernovae

Authors: Mei Takeda, Ryuichiro Hada, Ken-ichi Oohara, Toshifumi Futamase

Abstract: We investigate the possibility to have a constraint on the mass of thermal warm dark matter (WDM) particle from the expected data of the Wide Field Infrared Survey Telescope (WFIRST) survey if all the dark matter is warm. For this purpose we consider the lensing effect of large scale structure based on the warm dark matter scenario on the apparent magnitude of SNe Ia. We use HALOFIT as non-linear… ▽ More We investigate the possibility to have a constraint on the mass of thermal warm dark matter (WDM) particle from the expected data of the Wide Field Infrared Survey Telescope (WFIRST) survey if all the dark matter is warm. For this purpose we consider the lensing effect of large scale structure based on the warm dark matter scenario on the apparent magnitude of SNe Ia. We use HALOFIT as non-linear matter power spectrum and the variance of PDF. We preform a Fisher matrix analysis and obtain the lower bound of $m_{\rm WDM}>0.167$keV. △ Less

Submitted 18 May, 2019; originally announced May 2019.

Comments: 6 pages, 4 figures

A possible solution to the Hubble constant discrepancy -- Cosmology where the local volume expansion is driven by the domain average density

Authors: Masumi Kasai, Toshifumi Futamase

Abstract: The Hubble constant problem is the discrepancy between different measurements of the Hubble constant in different scales. We show that this problem can be resolved within the general relativistic framework of the perturbation theory in the inhomogeneous universe, with the help of spatial averaging procedure over a finite local domain in the $t=\mbox{const.}$ hypersurface. The idea presented in thi… ▽ More The Hubble constant problem is the discrepancy between different measurements of the Hubble constant in different scales. We show that this problem can be resolved within the general relativistic framework of the perturbation theory in the inhomogeneous universe, with the help of spatial averaging procedure over a finite local domain in the $t=\mbox{const.}$ hypersurface. The idea presented in this paper is unique in the sense that it has all of the following properties. a) It is based on the general relativistic perturbation theory, with ordinary dust matter only. No strange matter nor energy components are required. b) The employment of the spatially invariant averaging procedure on the finite domain is essential. c) The key is the first-order effect of the inhomogeneities in the linear perturbation theory. No non-linear effects are required. △ Less

Submitted 4 June, 2019; v1 submitted 21 April, 2019; originally announced April 2019.

Journal ref: Progress of Theoretical and Experimental Physics, Volume 2019, Issue 7, July 2019, 073E01

Forecasts of cosmological constraints from Type Ia supernovae including the weak-lensing convergence

Authors: Ryuichiro Hada, Toshifumi Futamase

Abstract: We investigate how the cosmological constraints from SNe Ia are improved by including the effects of weak-lensing convergence. To do so, we introduce the lognormal function as the convergence PDF modeling the lensing scatter of SN Ia magnitude, and apply a sample selection for SNe Ia to avoid strongly lensed samples. Comparing with the contribution of other uncertainties (e.g., the intrinsic magni… ▽ More We investigate how the cosmological constraints from SNe Ia are improved by including the effects of weak-lensing convergence. To do so, we introduce the lognormal function as the convergence PDF modeling the lensing scatter of SN Ia magnitude, and apply a sample selection for SNe Ia to avoid strongly lensed samples. Comparing with the contribution of other uncertainties (e.g., the intrinsic magnitude scatter), we find that the lensing effect is dominant at $z > 1$. Then forecasting the parameter constraints for the Wide-Field InfraRed Survey Telescope survey, we show that considering the weak-lensing effect, the constraints on the density parameters $Ω_m$ or $Ω_Λ$, and the dark energy equation of state $w$ are improved, especially for SNe Ia samples at higher redshift $z > 1$. Furthermore, we see that the degeneracy between the total mass of neutrino $Σm_ν$ and the (cold) dark matter density parameter $Ω_{c}$ can be resolved. △ Less

Submitted 25 June, 2019; v1 submitted 11 October, 2018; originally announced October 2018.

Comments: 15 pages, 4 figures, accepted for publication in JCAP; some changes to match published version

Analytical Noise Bias Correction for Weak Lensing Shear Analysis with ERA

Authors: Yuki Okura, Toshifumi Futamase

Abstract: Highly precise weak lensing shear measurement is required for statistical weak gravitational lensing analysis such as cosmic shear measurement to achieve severe constraint on the cosmological parameters. For this purpose, the accurate shape measurement of background galaxies is absolutely important in which any systematic error in the measurement should be carefully corrected. One of the main syst… ▽ More Highly precise weak lensing shear measurement is required for statistical weak gravitational lensing analysis such as cosmic shear measurement to achieve severe constraint on the cosmological parameters. For this purpose, the accurate shape measurement of background galaxies is absolutely important in which any systematic error in the measurement should be carefully corrected. One of the main systematic error comes from photon noise which is Poisson noise of flux from the atmosphere(noise bias). We investigate how the photon noise makes a systematic error in shear measurement within the framework of ERA method we developed in earlier papers and gives a practical correction method. The method is tested by simulations with real galaxy images with various conditions and it is confirmed that it can correct $80 \sim 90\%$ of the noise bias except for galaxies with very low signal to noise ratio. △ Less

Submitted 6 February, 2018; originally announced February 2018.

Comments: submitted to MNRAS. arXiv admin note: substantial text overlap with arXiv:1706.03903

Pixel Noise Effect and its Correction method by ERA Method for Precise Weak Gravitational Lensing Shear Measurement

Authors: Yuki Okura, Toshifumi Futamase

Abstract: Highly precise weak lensing shear measurement is required for statistical weak gravitational lensing analysis such as cosmic shear measurement to achieve severe constrain on the cosmological parameters. For this purpose any systematic error in the measurement should be corrected. One of the main systematic error comes from Pixel noise which is Poisson noise of flux from atmosphere. We investigate… ▽ More Highly precise weak lensing shear measurement is required for statistical weak gravitational lensing analysis such as cosmic shear measurement to achieve severe constrain on the cosmological parameters. For this purpose any systematic error in the measurement should be corrected. One of the main systematic error comes from Pixel noise which is Poisson noise of flux from atmosphere. We investigate how the pixel noise makes systematic error in shear measurement based on ERA method and develop the correction method. This method is tested by simulations with various conditions and it is confirmed that the correction method can correct $80 \sim 90\%$ of the systematic error except very low signal to noise ratio galaxies. △ Less

Submitted 13 June, 2017; originally announced June 2017.

Comments: 15 figures

Constraints on neutrino masses from the lensing dispersion of Type Ia supernovae

Authors: Ryuichiro Hada, Toshifumi Futamase

Abstract: We investigate how accurately the total mass of neutrinos is constrained from the magnitude dispersion of SNe Ia due to the effects of gravitational lensing. For this purpose, we use the propagation equation of light bundles in a realistic inhomogeneous universe and propose a sample selection for supernovae to avoid difficulties associated with small-scale effects such as strong lensing or shear e… ▽ More We investigate how accurately the total mass of neutrinos is constrained from the magnitude dispersion of SNe Ia due to the effects of gravitational lensing. For this purpose, we use the propagation equation of light bundles in a realistic inhomogeneous universe and propose a sample selection for supernovae to avoid difficulties associated with small-scale effects such as strong lensing or shear effects. With a fitting formula for the nonlinear matter power spectrum taking account of the effects of massive neutrinos, we find that in our model it is possible to obtain the upper limit $Σm_ν \simeq 1.0[{\rm eV}]$ for future optical imaging surveys with the Wide-Field InfraRed Survey Telescope and Large Synoptic Survey Telescope. Furthermore, we discuss how far we need to observe SNe Ia and to what extent we have to reduce the magnitude error except for lensing in order to realize the current tightest limit $Σm_ν < 0.2[{\rm eV}]$. △ Less

Submitted 28 September, 2016; v1 submitted 29 June, 2016; originally announced June 2016.

Comments: 19 pages, 5 figures, accepted for publication in ApJ; minor changes to match published version

Journal ref: Astrophys. J. 828: 112, 2016

The systematic error test for PSF correction in weak gravitational lensing shear measurement by The ERA Method by Idealizing PSF

Authors: Yuki Okura, Toshifumi Futamase

Abstract: We improve the ERA(Ellipticity of Re-smeared Artificial image) method of PSF(Point Spread Function) correction in weak lensing shear analysis in order to treat realistic shape of galaxies and PSF. This is done by re-smearing PSF and the observed galaxy image smeared by a RSF(Re-Smearing Function), and allows us to use a new PSF with a simple shape and to correct PSF effect without any approximat… ▽ More We improve the ERA(Ellipticity of Re-smeared Artificial image) method of PSF(Point Spread Function) correction in weak lensing shear analysis in order to treat realistic shape of galaxies and PSF. This is done by re-smearing PSF and the observed galaxy image smeared by a RSF(Re-Smearing Function), and allows us to use a new PSF with a simple shape and to correct PSF effect without any approximations and assumptions. We perform numerical test to show that the method applied for galaxies and PSF with some complicated shapes can correct PSF effect with systematic error less than 0.1%. We also apply ERA method for real data of Abell 1689 cluster to confirm that it is able to detect the systematic weak lensing shear pattern. The ERA method requires less than 0.1 or 1 second to correct PSF for each object in numerical test and real data analysis, respectively. △ Less

Submitted 24 June, 2016; originally announced June 2016.

Comments: Accepted by Astrophysical Journal

LoCuSS: Testing hydrostatic equilibrium in galaxy clusters

Authors: G. P. Smith, P. Mazzotta, N. Okabe, F. Ziparo, S. L. Mulroy, A. Babul, A. Finoguenov, I. G. McCarthy, M. Lieu, Y. Bahe, H. Bourdin, A. E. Evrard, T. Futamase, C. P. Haines, M. Jauzac, D. P. Marrone, R. Martino, P. E. May, J. E. Taylor, K. Umetsu

Abstract: We test the assumption of hydrostatic equilibrium in an X-ray luminosity selected sample of 50 galaxy clusters at $0.15<z<0.3$ from the Local Cluster Substructure Survey (LoCuSS). Our weak-lensing measurements of $M_{500}$ control systematic biases to sub-4 per cent, and our hydrostatic measurements of the same achieve excellent agreement between XMM-Newton and Chandra. The mean ratio of X-ray to… ▽ More We test the assumption of hydrostatic equilibrium in an X-ray luminosity selected sample of 50 galaxy clusters at $0.15<z<0.3$ from the Local Cluster Substructure Survey (LoCuSS). Our weak-lensing measurements of $M_{500}$ control systematic biases to sub-4 per cent, and our hydrostatic measurements of the same achieve excellent agreement between XMM-Newton and Chandra. The mean ratio of X-ray to lensing mass for these 50 clusters is $β_{\rm X}=0.95\pm0.05$, and for the 44 clusters also detected by Planck, the mean ratio of Planck mass estimate to LoCuSS lensing mass is $β_{\rm P}=0.95\pm0.04$. Based on a careful like-for-like analysis, we find that LoCuSS, the Canadian Cluster Comparison Project (CCCP), and Weighing the Giants (WtG) agree on $β_{\rm P}\simeq0.9-0.95$ at $0.15<z<0.3$. This small level of hydrostatic bias disagrees at $\sim5σ$ with the level required to reconcile Planck cosmology results from the cosmic microwave background and galaxy cluster counts. △ Less

Submitted 5 November, 2015; originally announced November 2015.

Comments: 5 pages, 1 figure, accepted by MNRAS

The ERA Method with Idealizing PSF for Precise Weak Gravitational Lensing Shear Analysis

Authors: Yuki Okura, Toshifumi Futamase

Abstract: We generalize ERA method of PSF correction for more realistic situations. The method re-smears the observed galaxy image(galaxy image smeared by PSF) and PSF image by an appropriate function called Re-Smearing Function(RSF) to make new images which have the same ellipticity with the lensed (before smeared by PSF) galaxy image. It has been shown that the method avoids a systematic error arising f… ▽ More We generalize ERA method of PSF correction for more realistic situations. The method re-smears the observed galaxy image(galaxy image smeared by PSF) and PSF image by an appropriate function called Re-Smearing Function(RSF) to make new images which have the same ellipticity with the lensed (before smeared by PSF) galaxy image. It has been shown that the method avoids a systematic error arising from an approximation in the usual PSF correction in moment method such as KSB for simple PSF shape. By adopting an idealized PSF we generalize ERA method applicable for arbitrary PSF. This is confirmed with simulated complex PSF shapes. We also consider the effect of pixel noise and found that the effect causes systematic overestimation. △ Less

Submitted 30 June, 2015; originally announced June 2015.

Comparison of Strong Gravitational Lens Model Software III. Direct and indirect semi-independent lens model comparisons of COSMOS J095930+023427, SDSS J1320+1644, SDSSJ1430+4105 and J1000+0021

Authors: Alan T. Lefor, Toshifumi Futamase

Abstract: Analysis of strong gravitational lensing data is important in this era of precision cosmology. The objective of the present study is to directly compare the analysis of strong gravitational lens systems using different lens model software and similarly parameterized models to understand the differences and limitations of the resulting models. The software lens model translation tool, HydraLens, wa… ▽ More Analysis of strong gravitational lensing data is important in this era of precision cosmology. The objective of the present study is to directly compare the analysis of strong gravitational lens systems using different lens model software and similarly parameterized models to understand the differences and limitations of the resulting models. The software lens model translation tool, HydraLens, was used to generate multiple models for four strong lens systems including COSMOS J095930+023427, SDSS J1320+1644, SDSSJ1430+4105 and J1000+0021. All four lens systems were modeled with PixeLens, Lenstool, glafic, and Lensmodel. The input data and parameterization of each lens model was similar for the four model programs used to highlight differences in the output results. The calculation of the Einstein radius and enclosed mass for each lens model was comparable. The results were more dissimilar if the masses of more than one lens potential were free-parameters. The image tracing algorithms of the software are different, resulting in different output image positions and differences in time delay and magnification calculations, as well as ellipticity and position angle of the resulting lens model. In a comparison of different software versions using identical model input files, results differed significantly when using two versions of the same software. These results further support the need for future lensing studies to include multiple lens models, use of open software, availability of lens model files use in studies, and computer challenges to develop new approaches. Future studies need a standard nomenclature and specification of the software used to allow improved interpretation, reproducibility and transparency of results. △ Less

Submitted 3 May, 2015; originally announced May 2015.

Comments: 15 pages, 4 figures

Central Mass Profiles of the Nearby Cool-core Galaxy Clusters Hydra A and A478

Authors: N. Okabe, K. Umetsu, T. Tamura, Y. Fujita, M. Takizawa, K. Matsushita, Y. Fukazawa, T. Futamase, M. Kawaharada, S. Miyazaki, Y. Mochizuki, K. Nakazawa, T. Ohashi, N. Ota, T. Sasaki, K. Sato, S. I. Tam

Abstract: We perform a weak-lensing study of the nearby cool-core galaxy clusters, Hydra A ($z=0.0538$) and A478 ($z=0.0881$), of which brightest cluster galaxies (BCGs) host powerful activities of active galactic nuclei (AGNs). For each cluster, the observed tangential shear profile is well described either by a single Navarro--Frenk--White model or a two-component model including the BCG as an unresolved… ▽ More We perform a weak-lensing study of the nearby cool-core galaxy clusters, Hydra A ($z=0.0538$) and A478 ($z=0.0881$), of which brightest cluster galaxies (BCGs) host powerful activities of active galactic nuclei (AGNs). For each cluster, the observed tangential shear profile is well described either by a single Navarro--Frenk--White model or a two-component model including the BCG as an unresolved point mass. For A478, we determine the BCG and its host-halo masses from a joint fit to weak-lensing and stellar photometry measurements. We find that the choice of initial mass functions (IMFs) can introduce a factor of two uncertainty in the BCG mass, whereas the BCG host halo mass is well constrained by data. We perform a joint analysis of weak-lensing and stellar kinematics data available for the Hydra A cluster, which allows us to constrain the central mass profile without assuming specific IMFs. We find that the central mass profile ($r<300$kpc) determined from the joint analysis is in excellent agreement with those from independent measurements, including dynamical masses estimated from the cold gas disk component, X-ray hydrostatic total mass estimates, and the central stellar mass estimated based on the Salpeter IMF. The observed dark-matter fraction around the BCG for Hydra A is found to be smaller than those predicted by adiabatic contraction models, suggesting the importance of other physical processes, such as the the AGN feedback and/or dissipationless mergers. △ Less

Submitted 10 December, 2015; v1 submitted 15 March, 2015; originally announced March 2015.

Comments: accepted for the publication in MNRAS, 15 pages 10 figures

The magnitude-redshift relation in a realistic inhomogeneous universe

Authors: Ryuichiro Hada, Toshifumi Futamase

Abstract: The light rays from a source are subject to a local inhomogeneous geometry generated by inhomogeneous matter distribution as well as the existence of collapsed objects. In this paper we investigate the effect of inhomogeneities and the existence of collapsed objects on the propagation of light rays and evaluate changes in the magnitude-redshift relation from the standard relationship found in a ho… ▽ More The light rays from a source are subject to a local inhomogeneous geometry generated by inhomogeneous matter distribution as well as the existence of collapsed objects. In this paper we investigate the effect of inhomogeneities and the existence of collapsed objects on the propagation of light rays and evaluate changes in the magnitude-redshift relation from the standard relationship found in a homogeneous FRW universe. We give the expression of the correlation function and the variance for the perturbation of apparent magnitude, and calculate it numerically by using the non-linear matter power spectrum. We use the lognormal probability distribution function for the density contrast and spherical collapse model to truncate the power spectrum in order to estimate the blocking effect by collapsed objects. We find that the uncertainties in $Ω_m$ is $\sim 0.02$, and that of $w$ is $\sim 0.04$. We also discuss a possible method to extract these effects from real data which contains intrinsic ambiguities associated with the absolute magnitude. △ Less

Submitted 23 December, 2014; v1 submitted 9 October, 2014; originally announced October 2014.

Comments: 19 pages, 4 figures, accepted to JCAP

Journal ref: JCAP12(2014)042

Universal Profiles of the Intracluster Medium from Suzaku X-Ray and Subaru Weak Lensing Obesrvations

Authors: N. Okabe, K. Umetsu, T. Tamura, Y. Fujita, M. Takizawa, Y. -Y. Zhang, K. Matsushita, T. Hamana, Y. Fukazawa, T. Futamase, M. Kawaharada, S. Miyazaki, Y. Mochizuki, K. Nakazawa, T. Ohashi, N. Ota, T. Sasaki, K. Sato, S. I. Tam

Abstract: We conduct a joint X-ray and weak-lensing study of four relaxed galaxy clusters (Hydra A, A478, A1689 and A1835) observed by both Suzaku and Subaru out to virial radii, with an aim to understand recently-discovered unexpected feature of the ICM in cluster outskirts. We show that the average hydrostatic-to-lensing total mass ratio for the four clusters decreases from \sim 70% to \sim 40% as the ove… ▽ More We conduct a joint X-ray and weak-lensing study of four relaxed galaxy clusters (Hydra A, A478, A1689 and A1835) observed by both Suzaku and Subaru out to virial radii, with an aim to understand recently-discovered unexpected feature of the ICM in cluster outskirts. We show that the average hydrostatic-to-lensing total mass ratio for the four clusters decreases from \sim 70% to \sim 40% as the overdensity contrast decreases from 500 to the virial value.The average gas mass fraction from lensing total mass estimates increases with cluster radius and agrees with the cosmic mean baryon fraction within the virial radius, whereas the X-ray-based gas fraction considerably exceeds the cosmic values due to underestimation of the hydrostatic mass. We also develop a new advanced method for determining normalized cluster radial profiles for multiple X-ray observables by simultaneously taking into account both their radial dependence and multivariate scaling relations with weak-lensing masses. Although the four clusters span a range of halo mass, concentration, X-ray luminosity and redshift, we find that the gas entropy, pressure, temperature and density profiles are all remarkably self-similar when scaled with the lensing M_200 mass and r_200 radius.The entropy monotonically increases out to \sim 0.5r_200 following the accretion shock heating model K(r)\propto r^1.1, and flattens at \simgt 0.5r_200.The universality of the scaled entropy profiles indicates that the thermalization mechanism over the entire cluster region (>0.1r_200) is controlled by gravitation in a common to all clusters, although the heating efficiency in the outskirts needs to be modified from the standard law.The bivariate scaling functions of the gas density and temperature reveal that the flattening of the outskirts entropy profile is caused by the steepening of the temperature, rather than the flattening of the gas density. △ Less

Submitted 18 July, 2014; v1 submitted 13 June, 2014; originally announced June 2014.

Comments: accepted for publication in PASJ : 13 pages 6 figures and 4 tables

A new weak lensing shear analysis method using ellipticity defined by 0th order moments

Authors: Yuki Okura, Toshifumi Futamase

Abstract: We developed a new method that uses ellipticity defined by 0th order moments (0th-ellipticity) for weak gravitational lensing shear analysis. Although there is a strong correlation between the ellipticity calculated using this approach and the usual ellipticity defined by the 2nd order moment, the ellipticity calculated here has a higher signal-to-noise ratio because it is weighted to the central… ▽ More We developed a new method that uses ellipticity defined by 0th order moments (0th-ellipticity) for weak gravitational lensing shear analysis. Although there is a strong correlation between the ellipticity calculated using this approach and the usual ellipticity defined by the 2nd order moment, the ellipticity calculated here has a higher signal-to-noise ratio because it is weighted to the central region of the image. These results were confirmed using data for Abell 1689 from the Subaru telescope. For shear analysis, we adopted the ellipticity of re-smeared artificial image (ERA) method for point spread function (PSF) correction, and we tested the precision of this 0th-ellipticity with simple simulation, then we obtained the same level of precision with the results of ellipticity defined by quadrupole moments. Thus, we can expect that weak lensing analysis using 0 shear will be improved in proportion to the statistical error. △ Less

Submitted 23 February, 2015; v1 submitted 7 May, 2014; originally announced May 2014.

Comments: 25pages, 8figures; published by A&A 2015

Journal ref: A&A 576, A63 (2015)

A New Method for PSF correction using the Ellipticity of Re-smeared Artificial Images in Weak Gravitational Lensing Shear Analysis

Authors: Yuki Okura, Toshifumi Futamase

Abstract: We propose a new method for Point Spread Function (PSF) correction in weak gravitational lensing shear analysis using an artificial image with the same ellipticity as the lensed image. This avoids the systematic error associated with the approximation in PSF correction used in previous approaches. We test the new method with simulated objects which have Gaussian or Cersic profiles smeared by a Gau… ▽ More We propose a new method for Point Spread Function (PSF) correction in weak gravitational lensing shear analysis using an artificial image with the same ellipticity as the lensed image. This avoids the systematic error associated with the approximation in PSF correction used in previous approaches. We test the new method with simulated objects which have Gaussian or Cersic profiles smeared by a Gaussian PSF, and confirm that there is no systematic error. △ Less

Submitted 28 March, 2014; originally announced March 2014.

Comparison of Strong Gravitational Lens Model Software I. Time delay and mass calculations are sensitive to changes in redshift and are model dependent

Authors: Alan T. Lefor, Toshifumi Futamase

Abstract: Analysis of strong gravitational lensing depends on software analysis of observational data. The purpose of this study was to evaluate the behavior of strong gravitational lens modeling software with changes in redshift. Four different strong gravitational lens software modeling codes were directly compared (Lenstool / glafic, two light traces mass codes, and GRALE / PixeLens, two non-light traces… ▽ More Analysis of strong gravitational lensing depends on software analysis of observational data. The purpose of this study was to evaluate the behavior of strong gravitational lens modeling software with changes in redshift. Four different strong gravitational lens software modeling codes were directly compared (Lenstool / glafic, two light traces mass codes, and GRALE / PixeLens, two non-light traces mass codes) in the analysis of model data as well as analysis of the giant gravitational quasar SDSSJ1004+4112. A generalized model for time delay calculation shows that calculated time delay is proportional to $D_{d}D_{s}/D_{ds}$.The percent change in time delays calculated for each system at each redshift tested were compared with percent change in the value of $D_{d}D_{s}/D_{ds}$. A simple point mass model was tested with each code. Five models were used with a constant $z_{lens}$ and a varying $z_{source}$, and five models with a constant $z_{source}$ and a varying $z_{lens}$. The effects of changing geometry were similarly investigated for SDSSJ1004+4112. In general, the changes in time delay were of a similar magnitude and direction although, some calculated time delays varied by as much as 30 percent from changes in $D_{d}D_{s}/D_{ds}$. Changes in calculated mass for the point mass model with a constant $z_{source}$ were almost identical to changes in $D_{d}D_{s}/D_{ds}$ for three of the four codes tested. These data demonstrate the effect of changes in redshift on parameters calculated by each of the codes as compared to changes in $D_{d}D_{s}/D_{ds}$. The paucity of existing direct comparison studies of strong gravitational lensing supports the need for more studies of this kind. These results show that even small changes in redshift affects the calculation of time delay and mass, and that the effect on the calculations is dependent on the particular software used. △ Less

Submitted 17 July, 2013; originally announced July 2013.

Toward understanding the anisotropic point spread function of Suprime-Cam and its impact on cosmic shear measurement

Authors: Takashi Hamana, Satoshi Miyazaki, Yuki Okura, Tomohiro Okamura, Toshifumi Futamase

Abstract: We examined the anisotropic point spread function (PSF) of Suprime-Cam data utilizing dense star field data. We decomposed the PSF ellipticities into three components, the optical aberration, atmospheric turbulence, and chip-misalignment in an empirical manner, and evaluated the amplitude of each component. We found that, for long-exposure data, the optical aberration has the largest contribution… ▽ More We examined the anisotropic point spread function (PSF) of Suprime-Cam data utilizing dense star field data. We decomposed the PSF ellipticities into three components, the optical aberration, atmospheric turbulence, and chip-misalignment in an empirical manner, and evaluated the amplitude of each component. We found that, for long-exposure data, the optical aberration has the largest contribution to the PSF ellipticities, which could be modeled well by a simple analytic function based on the lowest-order aberration theory. The statistical properties of PSF ellipticities resulting from the atmospheric turbulence were investigated by using the numerical simulations. The simulation results are in a reasonable agreement with the observed data. It is also found that the optical PSF can be well corrected by the standard correction method with a polynomial fitting function. However, for the atmospheric PSF, its correction is affected by the common limitation caused by sparse sampling of PSFs due to a limited number of stars. We also examined the effects of the residual PSF anisotropies on Suprime-Cam cosmic shear data. We found that the shape and amplitude of the B-mode shear variance are broadly consistent with those of the residual PSF ellipticities measured from the dense star field data. This indicates that most of the sources of residual systematic are understood, which is an important step for cosmic shear statistics to be a practical tool of the precision cosmology. △ Less

Submitted 20 November, 2013; v1 submitted 17 April, 2013; originally announced April 2013.

Comments: 17 pages, 21 figures, accepted for publication in PASJ

Journal ref: Publ. Astron. Soc. Japan 65, 104, 2013

Subaru Weak-Lensing Survey of Dark Matter Subhalos in the Coma Cluster : Subhalo Mass Function and Statistical Properties

Authors: Nobuhiro Okabe, Toshifumi Futamase, Masaru Kajisawa, Risa Kuroshima

Abstract: We present a 4 deg^2 weak gravitational lensing survey of subhalos in the very nearby Coma cluster using the Subaru/Suprime-Cam. The large apparent size of cluster subhalos allows us to measure the mass of 32 subhalos detected in a model-independent manner, down to the order of 10^-3 of the virial mass of the cluster. Weak-lensing mass measurements of these shear-selected subhalos enable us to inv… ▽ More We present a 4 deg^2 weak gravitational lensing survey of subhalos in the very nearby Coma cluster using the Subaru/Suprime-Cam. The large apparent size of cluster subhalos allows us to measure the mass of 32 subhalos detected in a model-independent manner, down to the order of 10^-3 of the virial mass of the cluster. Weak-lensing mass measurements of these shear-selected subhalos enable us to investigate subhalo properties and the correlation between subhalo masses and galaxy luminosities for the first time. The mean distortion profiles stacked over subhalos show a sharply truncated feature which is well-fitted by a Navarro-Frenk-White (NFW) mass model with the truncation radius, as expected due to tidal destruction by the main cluster. We also found that subhalo masses, truncation radii, and mass-to-light ratios decrease toward the cluster center. The subhalo mass function, dn/dln M_sub, in the range of 2 orders of magnitude in mass, is well described by a single power law or a Schechter function. Best-fit power indices of 1.09_-0.32^+0.42 for the former model and 0.99_-0.23^+0.34 for the latter, are in remarkable agreement with slopes of ~0.9-1.0 predicted by the cold dark matter paradigm. The tangential distortion signals in the radial range of 0.02-2Mpc/h from the cluster center show a complex structure which is well described by a composition of three mass components of subhalos, the NFW mass distribution as a smooth component of the main cluster, and a lensing model from a large scale structure behind the cluster. Although the lensing signals are 1 order of magnitude lower than those for clusters at z~0.2, the total signal-to-noise ratio, S/N=13.3, is comparable to, or higher, because the enormous number of background source galaxies compensates for the low lensing efficiency of the low lensing efficiency of the nearby cluster. △ Less

Submitted 3 March, 2014; v1 submitted 8 April, 2013; originally announced April 2013.

Comments: 30 pages, 18 figures, 9 tables, ApJ in press. Full resolution version is available at http://www.asiaa.sinica.edu.tw/~okabe/files/coma_survey.pdf

Relation between standard perturbation theory and regularized multi-point propagator method

Authors: Naonori S. Sugiyama, Toshifumi Futamase

Abstract: We investigate the relation between the regularized multi-propagator method, called "Reg PT", and the standard perturbation theory. Reg PT is one of the most successful models to describe nonlinear evolution of dark matter fluctuations. However, Reg PT is a mathematically unproven interpolation formula between the large-scale solution calculated by the standard perturbation theory and the limiting… ▽ More We investigate the relation between the regularized multi-propagator method, called "Reg PT", and the standard perturbation theory. Reg PT is one of the most successful models to describe nonlinear evolution of dark matter fluctuations. However, Reg PT is a mathematically unproven interpolation formula between the large-scale solution calculated by the standard perturbation theory and the limiting solution in the small scale calculated by the multi-point propagator method. In this paper, we give an alternative explanation for Reg PT in the context of the standard perturbation theory, showing that Reg PT does not ever have more effective information on nonlinear matter evolution than the standard perturbation theory. In other words, the solutions of the standard perturbation theory reproduce the results of $N$-body simulations better than those of Reg PT, especially at the high-$k$ region. This fact means that the standard perturbation theory at the two-loop level is still one of the best predictions of the nonlinear power spectrum to date. Nevertheless, the standard perturbation theory has not been preferred because of the divergent behavior of the solution at small scales. To solve this problem, we also propose a modified standard perturbation theory which avoids the divergence. △ Less

Submitted 16 May, 2013; v1 submitted 11 March, 2013; originally announced March 2013.

Comments: 15pages and 2 figures

Journal ref: ApJ, 769, 106 (2013)

LoCuSS: The Mass Density Profile of Massive Galaxy Clusters at z=0.2

Authors: Nobuhiro Okabe, Graham P. Smith, Keiichi Umetsu, Masahiro Takada, Toshifumi Futamase

Abstract: We present a stacked weak-lensing analysis of an approximately mass-selected sample of 50 galaxy clusters at 0.15<z<0.3, based on observations with Suprime-Cam on the Subaru Telescope. We develop a new method for selecting lensed background galaxies from which we estimate that our sample of red background galaxies suffers just 1% contamination. We detect the stacked tangential shear signal from th… ▽ More We present a stacked weak-lensing analysis of an approximately mass-selected sample of 50 galaxy clusters at 0.15<z<0.3, based on observations with Suprime-Cam on the Subaru Telescope. We develop a new method for selecting lensed background galaxies from which we estimate that our sample of red background galaxies suffers just 1% contamination. We detect the stacked tangential shear signal from the full sample of 50 clusters, based on this red sample of background galaxies, at a total signal-to-noise ratio of S/N=32.7. The Navarro-Frenk-White model is an excellent fit to the data, yielding sub-10% statistical precision on mass and concentration: Mvir=7.19^{+0.53}_{-0.50}\times10^{14}h^{-1}Msol, cvir=5.41^{+0.49}_{-0.45} (c_{200}=4.22^{+0.40}_{-0.36}). Tests of a range of possible systematic errors, including shear calibration and stacking-related issues, indicate that they are sub-dominant to the statistical errors. The concentration parameter obtained from stacking our approximately mass-selected cluster sample is broadly in line with theoretical predictions. Moreover, the uncertainty on our measurement is comparable with the differences between the different predictions in the literature. Overall our results highlight the potential for stacked weak-lensing methods to probe the mean mass density profile of cluster-scale dark matter halos with upcoming surveys, including Hyper-Suprime-Cam, Dark Energy Survey, and KIDS. △ Less

Submitted 6 May, 2013; v1 submitted 12 February, 2013; originally announced February 2013.

Comments: 7 pages, 3 figures, and 2 tables. accepted to ApJL

Journal ref: ApJ Letters 769, 35 (2013)

Toward an accurate Dark Matter power spectrum beyond BAO scales

Authors: Naonori S. Sugiyama, Toshifumi Futamase

Abstract: We propose a new approximated expression for non-linear Dark Matter power spectrum much beyond BAO scales. The proposed expression agrees with the result of N-body simulation with the accuracy better than 2 % up to k=1.0 [h/Mpc] and k=0.7 [h/Mpc] at z=3.0 and z=1.0, respectively. Even at z=0.35, the accuracy remains within 10 % up to k=0.8 [h/Mpc]. In doing so, we used an approximation for the ker… ▽ More We propose a new approximated expression for non-linear Dark Matter power spectrum much beyond BAO scales. The proposed expression agrees with the result of N-body simulation with the accuracy better than 2 % up to k=1.0 [h/Mpc] and k=0.7 [h/Mpc] at z=3.0 and z=1.0, respectively. Even at z=0.35, the accuracy remains within 10 % up to k=0.8 [h/Mpc]. In doing so, we used an approximation for the kernel functions used in the Standard Perturbation Theory (SPT) which is also used to prove the Reg PT proposed by Bernardeau et al. (2011). △ Less

Submitted 30 November, 2012; v1 submitted 28 October, 2012; originally announced October 2012.

Comments: 4 pages, 1 figure

An Application of Wiener Hermite Expansion to Non-linear Evolution of Dark Matter

Authors: Naonori S. Sugiyama, Toshifumi Futamase

Abstract: We apply the Wiener Hermite (WH) expansion to the non-linear evolution of Large-Scale Structure, and obtain an approximate expression for the matter power spectrum in full order of the expansion. This method allows us to expand any random function in terms of an orthonormal bases in space of random functions in such a way that the first order of the expansion expresses the Gaussian distribution, a… ▽ More We apply the Wiener Hermite (WH) expansion to the non-linear evolution of Large-Scale Structure, and obtain an approximate expression for the matter power spectrum in full order of the expansion. This method allows us to expand any random function in terms of an orthonormal bases in space of random functions in such a way that the first order of the expansion expresses the Gaussian distribution, and others are the deviation from the Gaussianity. It is proved that the Wiener Hermite expansion is mathematically equivalent with the $Γ$-expansion approach in the Renormalized Perturbation Theory (RPT). While an exponential behavior in the high-$k$ limit has been proved for the mass density and velocity fluctuations of Dark Matter in the RPT, we reprove the behavior in the context of the Wiener Hermite expansion using the result of Standard Perturbation Theory (SPT). We propose a new approximate expression for the matter power spectrum which interpolates the low-$k$ expression corresponding to 1-loop level in SPT and the high-$k$ expression obtained by taking a high-$k$ limit of the WH expansion. The validity of our prescription is specifically verified by comparing with the 2-loop solutions of the SPT. The proposed power spectrum agrees with the result of $N$-body simulation with the accuracy better than 1% or 2% in a range of the BAO scales, where the wave number is about $k= 0.2 \sim 0.4$ $h{\rm Mpc^{-1}}$ at $z=0.5-3.0$. This accuracy is comparable to or slightly less than the ones in the Closure Theory whose the fractional difference from N body result is within 1%. A merit of our method is that the computational time is very short because only single and double integrals are involved in our solution. △ Less

Submitted 30 November, 2012; v1 submitted 5 October, 2012; originally announced October 2012.

Comments: 30 pages

Journal ref: ApJ, 760, 114 (2012)

Elliptical Weighted HOLICs for Weak Lensing Shear Measurement part3:Random Count Noise Effect for Image's Moments in Weak Lensing Analysis

Authors: Yuki Okura, Toshifumi Futamase

Abstract: This is the third paper on the improvements of systematic errors in our weak lensing analysis using an elliptical weight function, called E-HOLICs. In the previous papers we have succeeded in avoiding error which depends on ellipticity of background image. In this paper, we investigate the systematic error which depends on signal to noise ratio of background image. We find that the origin of the e… ▽ More This is the third paper on the improvements of systematic errors in our weak lensing analysis using an elliptical weight function, called E-HOLICs. In the previous papers we have succeeded in avoiding error which depends on ellipticity of background image. In this paper, we investigate the systematic error which depends on signal to noise ratio of background image. We find that the origin of the error is the random count noise which comes from Poisson noise of sky counts. Random count noise makes additional moments and centroid shift error, and those 1st orders are canceled in averaging, but 2nd orders are not canceled. We derived the equations which corrects these effects in measuring moments and ellipticity of the image and test their validity using simulation image. We find that the systematic error becomes less than 1% in the measured ellipticity for objects with $S/N>3$. △ Less

Submitted 17 August, 2012; originally announced August 2012.

Comments: 44pages, 20figures, submitted for ApJ

δN formalism

Authors: Naonori S. Sugiyama, Eiichiro Komatsu, Toshifumi Futamase

Abstract: Precise understanding of nonlinear evolution of cosmological perturbations during inflation is necessary for the correct interpretation of measurements of non-Gaussian correlations in the cosmic microwave background and the large-scale structure of the universe. The "δN formalism" is a popular and powerful technique for computing non-linear evolution of cosmological perturbations on large scales.… ▽ More Precise understanding of nonlinear evolution of cosmological perturbations during inflation is necessary for the correct interpretation of measurements of non-Gaussian correlations in the cosmic microwave background and the large-scale structure of the universe. The "δN formalism" is a popular and powerful technique for computing non-linear evolution of cosmological perturbations on large scales. In particular, it enables us to compute the curvature perturbation, ζ, on large scales without actually solving perturbed field equations. However, people often wonder why this is the case. In order for this approach to be valid, the perturbed Hamiltonian constraint and matter-field equations on large scales must, with a suitable choice of coordinates, take on the same forms as the corresponding unperturbed equations. We find that this is possible when (1) the unperturbed metric is given by a homogeneous and isotropic Friedmann-Lemaître-Robertson-Walker metric; and (2) on large scales and with a suitable choice of coordinates, one can ignore the shift vector (g0i) as well as time-dependence of tensor perturbations to gij/a2(t) of the perturbed metric. While the first condition has to be assumed a priori, the second condition can be met when (3) the anisotropic stress becomes negligible on large scales. However, in order to explicitly show that the second condition follows from the third condition, one has to use gravitational field equations, and thus this statement may depend on the details of theory of gravitation. Finally, as the δN formalism uses only the Hamiltonian constraint and matter-field equations, it does not a priori respect the momentum constraint. We show that the violation of the momentum constraint only yields a decaying mode solution for ζ, and the violation vanishes when the slow-roll conditions are satisfied. △ Less

Submitted 22 February, 2013; v1 submitted 5 August, 2012; originally announced August 2012.

Comments: 10 pages

Journal ref: Phys. Rev. D 87, 023530 (2013)

A Systematic Review of Strong Gravitational Lens Modeling Software

Authors: Alan T. Lefor, Toshifumi Futamase, Mohammad Akhlaghi

Abstract: Despite expanding research activity in gravitational lens modeling, there is no particular software which is considered a standard. Much of the gravitational lens modeling software is written by individual investigators for their own use. Some gravitational lens modeling software is freely available for download but is widely variable with regard to ease of use and quality of documentation. This r… ▽ More Despite expanding research activity in gravitational lens modeling, there is no particular software which is considered a standard. Much of the gravitational lens modeling software is written by individual investigators for their own use. Some gravitational lens modeling software is freely available for download but is widely variable with regard to ease of use and quality of documentation. This review of 13 software packages was undertaken to provide a single source of information. Gravitational lens models are classified as parametric models or non-parametric models, and can be further divided into research and educational software. Software used in research includes the GRAVLENS package (with both gravlens and lensmodel), Lenstool, LensPerfect, glafic, PixeLens, SimpLens, Lensview, and GRALE. In this review, GravLensHD, G-Lens, Gravitational Lensing, lens and MOWGLI are categorized as educational programs that are useful for demonstrating various aspects of lensing. Each of the 13 software packages is reviewed with regard to software features (installation, documentation, files provided, etc.) and lensing features (type of model, input data, output data, etc.) as well as a brief review of studies where they have been used. Recent studies have demonstrated the utility of strong gravitational lensing data for mass mapping, and suggest increased use of these techniques in the future. Coupled with the advent of greatly improved imaging, new approaches to modeling of strong gravitational lens systems are needed. This is the first systematic review of strong gravitational lens modeling software, providing investigators with a starting point for future software development to further advance gravitational lens modeling research. △ Less

Submitted 20 May, 2013; v1 submitted 20 June, 2012; originally announced June 2012.

Elliptical Weighted HOLICs for Weak Lensing Shear Measurement part2:PSF correction and application to Abell 370

Authors: Yuki Okura, Toshifumi Futamase

Abstract: We have developed a new method(E-HOLICs) of estimating gravitational shear by adopting an elliptical weight function to measure background galaxy images in our previous paper. Following to the previous paper where isotropic Point Spread Function(PSF) correction is calculated, we consider an anisotropic PSF correction in this paper in order to apply E-HOLICs for real data. A an example, E-HOLICs is… ▽ More We have developed a new method(E-HOLICs) of estimating gravitational shear by adopting an elliptical weight function to measure background galaxy images in our previous paper. Following to the previous paper where isotropic Point Spread Function(PSF) correction is calculated, we consider an anisotropic PSF correction in this paper in order to apply E-HOLICs for real data. A an example, E-HOLICs is applied to Subaru data of massive and compact galaxy clusters A370, and is able to detect double peaks in the central region of the cluster consistent with the analysis of strong lensing. We also study the systematic error in E-HOLICs using STEP II simulation. In particular we consider the dependences of signal to noise ratio "S/N" of background galaxies in the shear estimation. Although E-HOLICs does improve systematic error due to the ellipticity dependence as shown in part 1, a systematic error due to the S/N dependence remains, namely E-HOLICs underestimates shear when background galaxies with low S/N objects are used. We discuss a possible improvement of the S/N dependence. △ Less

Submitted 28 July, 2011; originally announced July 2011.

Comments: 28 pages, 15 figures, submitted ApJ

Non-Gaussianity Consistency Relation for Multi-field Inflation

Authors: Naonori S. Sugiyama, Eiichiro Komatsu, Toshifumi Futamase

Abstract: While detection of the "local form" bispectrum of primordial perturbations would rule out all single-field inflation models, multi-field models would still be allowed. We show that multi-field models described by the $δN$ formalism obey an inequality between $f_{\rm NL}$ and one of the local-form {\it trispectrum} amplitudes, $τ_{\rm NL}$, such that $τ_{\rm NL}>\frac12(\frac65f_{\rm NL})^2$ with a… ▽ More While detection of the "local form" bispectrum of primordial perturbations would rule out all single-field inflation models, multi-field models would still be allowed. We show that multi-field models described by the $δN$ formalism obey an inequality between $f_{\rm NL}$ and one of the local-form {\it trispectrum} amplitudes, $τ_{\rm NL}$, such that $τ_{\rm NL}>\frac12(\frac65f_{\rm NL})^2$ with a possible logarithmic scale dependence, provided that 2-loop terms are small. Detection of a violation of this inequality would rule out most of multi-field models, challenging inflation as a mechanism for generating the primoridal perturbations. △ Less

Submitted 24 May, 2011; v1 submitted 19 January, 2011; originally announced January 2011.

Comments: 5 pages. Accepted for publication in Physical Review Letters

Report number: TCC-040-10

Journal ref: Phys.Rev.Lett.106:251301,2011

Elliptical Weighted HOLICs for Weak Lensing Shear Measurement. part1:Definitions and isotropic PSF correction

Authors: Yuki Okura, Toshifumi Futamase

Abstract: We develop a new method to estimate gravitational shear by adopting an elliptical weight function to measure background galaxy images. In doing so, we introduce a new concept of "zero plane" which is an imaginal source plane where shapes of all sources are perfect circles, and regard the intrinsic shear as the result of an imaginal lensing distortion. This makes the relation between the observed s… ▽ More We develop a new method to estimate gravitational shear by adopting an elliptical weight function to measure background galaxy images. In doing so, we introduce a new concept of "zero plane" which is an imaginal source plane where shapes of all sources are perfect circles, and regard the intrinsic shear as the result of an imaginal lensing distortion. This makes the relation between the observed shear, the intrinsic shear and lensing distortion more simple and thus higher-order calculation more easy. The elliptical weight function allows us to measure the mutiplemoment of shape of background galaxies more precisely by weighting highly to brighter parts of image and moreover to reduce systematic error due to insufficient expansion of the weight function in the original approach of KSB. Point Spread Function(PSF) correction in E-HOLICs methods becomes more complicated than those in KSB methods. In this paper we studied isotropic PSF correction in detail. By adopting the lensing distortion as the ellipticity of the weight function, we are able to show that the shear estimation in E-HOLICs method reduces to solve a polynomial in the absolute magnitude of the distortion. We compare the systematic errors between our approach and KSB using STEP2 simulation. It is confirmed that KSB method overestimate the input shear for images with large ellipticities, and E-HOLICs correctly estimate the input shear even for such images. Anisotropic PSF correction and analysis of real data will be presented in forthcoming paper. △ Less

Submitted 3 October, 2010; originally announced October 2010.

Comments: 30 pages, 8 figures, submitted to ApJ

LoCuSS: Calibrating Mass-Observable Scaling Relations for Cluster Cosmology with Subaru Weak Lensing Observations

Authors: Nobuhiro Okabe, Yu-Ying Zhang, Alexis Finoguenov, Masahiro Takada, Graham P. Smith, Keiichi Umetsu, Toshifumi Futamase

Abstract: We present a joint weak-lensing/X-ray study of galaxy cluster mass-observable scaling relations, motivated by the critical importance of accurate calibration of mass proxies for future X-ray missions, including eROSITA. We use a sample of 12 clusters at z\simeq0.2 that we have observed with Subaru and XMM-Newton to construct relationships between the weak-lensing mass (M), and three X-ray observab… ▽ More We present a joint weak-lensing/X-ray study of galaxy cluster mass-observable scaling relations, motivated by the critical importance of accurate calibration of mass proxies for future X-ray missions, including eROSITA. We use a sample of 12 clusters at z\simeq0.2 that we have observed with Subaru and XMM-Newton to construct relationships between the weak-lensing mass (M), and three X-ray observables: gas temperature (T), gas mass (Mgas), and quasi-integrated gas pressure (Yx) at overdensities of Δ=2500, 1000, and 500 with respect to the critical density. We find that Mgas at Δ\le1000 appears to be the most promising mass proxy of the three, because it has the lowest intrinsic scatter in mass at fixed observable: σ_lnM\simeq0.1, independent of cluster dynamical state. The scatter in mass at fixed T and Yx is a factor of \sim2-3 larger than at fixed Mgas, which are indicative of the structural segregation that we find in the M-T and M-Yx relationships. Undisturbed clusters are found to be \sim40% and \sim20% more massive than disturbed clusters at fixed T and Yx respectively at \sim2σsignificance. In particular, A1914 - a well-known merging cluster - significantly increases the scatter and lowers the the normalization of the relation for disturbed clusters. We also investigated the covariance between intrinsic scatter in M-Mgas and M-T relations, finding that they are positively correlated. This contradicts the adaptive mesh refinement simulations that motivated the idea that Yx may be a low scatter mass proxy, and agrees with more recent smoothed particle hydrodynamic simulations based on the Millennium Simulation. We also propose a method to identify a robust mass proxy based on principal component analysis. The statistical precision of our results are limited by the small sample size and the presence of the extreme merging cluster in our sample. △ Less

Submitted 20 August, 2010; v1 submitted 22 July, 2010; originally announced July 2010.

Comments: 13 pages, 6 figures : ApJ in press : proof ver

Weak Lensing Mass Measurements of Substructures in COMA Cluster with Subaru/Suprime-Cam

Authors: Nobuhiro Okabe, Yuki Okura, Toshifumi Futamase

Abstract: We obtain the projected mass distributions for two Subaru/Suprime-Cam fields in the southwest region (r\simlt 60') of the Coma cluster (z=0.0236) by weak lensing analysis and detect eight subclump candidates. We quantify the contribution of background large-scale structure (LSS) on the projected mass distributions using SDSS multi-bands and photometric data, under the assumption of mass-to-light… ▽ More We obtain the projected mass distributions for two Subaru/Suprime-Cam fields in the southwest region (r\simlt 60') of the Coma cluster (z=0.0236) by weak lensing analysis and detect eight subclump candidates. We quantify the contribution of background large-scale structure (LSS) on the projected mass distributions using SDSS multi-bands and photometric data, under the assumption of mass-to-light ratio for field galaxies. We find that one of eight subclump candidates, which is not associated with any member galaxies, is significantly affected by LSS lensing. The mean projected mass for seven subclumps extracted from the main cluster potential is <M_2D^(corr)> = (5.06\pm1.30)10^12h^-1 M_sun after a LSS correction. A tangential distortion profile over an ensemble of subclumps is well described by a truncated singular-isothermal sphere model and a truncated NFW model. A typical truncated radius of subclumps, r_t\simeq 35 h^-1 kpc, is derived without assuming any relations between mass and light for member galaxies. The radius coincides well with the tidal radius, \sim42 h^-1 kpc, of the gravitational force of the main cluster. Taking into account the incompleteness of data area, a projection effect and spurious lensing peaks, it is expected that mass of cluster substructures account for 19 percent of the virial mass, with 13 percent statistical error. The mass fraction of cluster substructures is in rough agreement with numerical simulations. △ Less

Submitted 27 February, 2010; v1 submitted 14 January, 2010; originally announced January 2010.

Comments: ApJ, accepted, 16 pages, 10 figures and 4 tables. High-resolution pictures available at http://www.asiaa.sinica.edu.tw/~okabe/files/comaWL.pdf

Distance-Redshift Relation in a Realistic Inhomogeneous Universe

Authors: Tomohiro Okamura, Toshifumi Futamase

Abstract: We investigate the distance-redshift relation in a realistic inhomogeneous universe where the mass distribution is described by the mass function of Sheth and Tormen. It is found that the derived distance deviates systematically from the standard distance up to 10% depending on the choice of the lowest halo mass in which baryonic matter condensed to form luminous object such as galaxies. Remarka… ▽ More We investigate the distance-redshift relation in a realistic inhomogeneous universe where the mass distribution is described by the mass function of Sheth and Tormen. It is found that the derived distance deviates systematically from the standard distance up to 10% depending on the choice of the lowest halo mass in which baryonic matter condensed to form luminous object such as galaxies. Remarkably the derived distance is well approximated by the Dyer-Roeder distance if we choose the clumpiness parameter αcalculated by our model.We also discuss the effect of inhomogeneities in the determination of dark energy parameter in the supernovae observation, and find that this effect must be taken into account for the future high redshift supernovae observation. △ Less

Submitted 21 July, 2009; v1 submitted 7 May, 2009; originally announced May 2009.

Comments: accepted for publication in PTP, 9 pages including 2 tables and 3 figures

Journal ref: Prog.Theor.Phys.122:511-520,2009

LoCuSS: Subaru Weak Lensing Study of 30 Galaxy Clusters

Authors: Nobuhiro Okabe, Masahiro Takada, Keiichi Umetsu, Toshifumi Futamase, Graham P. Smith

Abstract: (Abridged) We use Subaru data to conduct a detailed weak-lensing study of the dark matter distribution in a sample of 30 X-ray luminous galaxy clusters at 0.15<z<0.3. A weak-lensing signal is detected at high statistical significance in each cluster, the total detection S/N ranging from 5 to 13. In this paper we concentrate on fitting spherical models to the tangential distortion profiles of the c… ▽ More (Abridged) We use Subaru data to conduct a detailed weak-lensing study of the dark matter distribution in a sample of 30 X-ray luminous galaxy clusters at 0.15<z<0.3. A weak-lensing signal is detected at high statistical significance in each cluster, the total detection S/N ranging from 5 to 13. In this paper we concentrate on fitting spherical models to the tangential distortion profiles of the clusters. When the models are fitted to the clusters individually, we are unable to discriminate statistically between SIS and NFW models. However when the tangential distortion profiles of the individual clusters are combined, and models fitted to the stacked profile, the SIS model is rejected at 6- and 11-sigma, respectively, for low- and high-mass bins. We also use the individual cluster NFW model fits to investigate the relationship between cluster mass (M_vir) and concentration (c_vir), finding an anti-correlation of c_vir and M_vir. The best-fit c_vir-M_vir relation is: c_vir(M_vir) propto M_vir^{-alpha} with alpha=0.41+/-0.19 -- i.e. a non-zero slope is detected at 2sigma significance. We then investigate the optimal radius within which to measure cluster mass, finding that the typical fractional errors are improved to sigma(M_Delta)/M_Delta ~ 0.1-0.2 for cluster masses at higher over-densities Delta=500-2000, from 0.2-0.3 for the virial over-density (~110). Further comparisons between mass measurements based on spherical model fitting and the model-independent aperture mass method reveal that the 2D aperture mass enclosed within a cylinder of a given aperture radius is systematically greater than the 3D spherical mass obtained from NFW model fitting: M_2D/M_3D= 1.34 and 1.40 for Delta=500 and 110, respectively. The amplitude of this effect agrees well with that predicted by integrating the NFW model along the line-of-sight. △ Less

Submitted 22 March, 2011; v1 submitted 5 March, 2009; originally announced March 2009.

Comments: 30 pages, 15 figures, revised version accepted for publication in PASJ. The version including mass maps of each cluster (the submitted version) is available from http://www.astr.tohoku.ac.jp/~okabe/files/locuss_astroph.pdf

Journal ref: Publ. Astron. Soc. Jap. 62:811-870,2010

Gravitational Lensing Effect on the Two-Point Correlation Function of Hot Spots in 21 cm fluctuations

Authors: H. Tashiro, T. Futamase

Abstract: We investigate the weak gravitational lensing effect on the two-point correlation function of local maxima (hot spots) in the cosmic 21 cm fluctuation map. The intrinsic two-point function has a pronounced depression feature around the angular scale of $θ\sim 40$ arcmin, which depends on the observed frequency and corresponds to the scale of the acoustic oscillation of cosmic plasma before the r… ▽ More We investigate the weak gravitational lensing effect on the two-point correlation function of local maxima (hot spots) in the cosmic 21 cm fluctuation map. The intrinsic two-point function has a pronounced depression feature around the angular scale of $θ\sim 40$ arcmin, which depends on the observed frequency and corresponds to the scale of the acoustic oscillation of cosmic plasma before the recombination. It is found that the weak lensing induces a large w-dependent smoothing at that scale where w is the equation of state parameter of dark energy and thus provides a useful constraints on the dark energy property combined with the depression angular scales on the two-point correlation function. △ Less

Submitted 15 January, 2009; originally announced January 2009.

Comments: 7 pages, submitted to MNRAS

A New Method for measuring Weak Gravitational Lensing Shear using Higher Order Spin-2 HOLICs

Authors: Yuki Okura, Toshifumi Futamase

Abstract: We investigate the possibility to use higher order moments of gravitational lensed images in the weak lensing analysis. For this purpose we employ spin-2 component of HOLICs(Higher Order Lensing Characteristics) developed by us. We test the weak lensing analysis with spin-2 HOLICs using actual, ground based Subaru observations of the massive galaxy cluster A1689(z=0.183). It turns out that spin-… ▽ More We investigate the possibility to use higher order moments of gravitational lensed images in the weak lensing analysis. For this purpose we employ spin-2 component of HOLICs(Higher Order Lensing Characteristics) developed by us. We test the weak lensing analysis with spin-2 HOLICs using actual, ground based Subaru observations of the massive galaxy cluster A1689(z=0.183). It turns out that spin-2 HOLICs of order up to 8 are sufficiently applicable for weak lensing analysis after correcting PSF anisotropy as well as isotropic PSF smearing. △ Less

Submitted 14 April, 2009; v1 submitted 29 May, 2008; originally announced May 2008.

Comments: Accepted for publication in ApJ; 10 pages, 7 figures, uses aastex.cls

Journal ref: Astrophys.J.699:143-149,2009

Validity of strong lensing statistics for constraints on the galaxy evolution model

Authors: Akiko Matsumoto, Toshifumi Futamase

Abstract: We examine the usefulness of the strong lensing statistics to constrain the evolution of the number density of lensing galaxies by adopting the values of the cosmological parameters determined by recent WMAP observation. For this purpose, we employ the lens-redshift test proposed by Kochanek (1992) and constrain the parameters in two evolution models, simple power-law model characterized by the… ▽ More We examine the usefulness of the strong lensing statistics to constrain the evolution of the number density of lensing galaxies by adopting the values of the cosmological parameters determined by recent WMAP observation. For this purpose, we employ the lens-redshift test proposed by Kochanek (1992) and constrain the parameters in two evolution models, simple power-law model characterized by the power law indexes $nu_{n}$ and $ν_{v}$ and the evolution model by Mitchell et al. (2005) based on CDM structure formation scenario. We use the well-defined lens sample from the Sloan Digital Sky Survey (SDSS) and this is similarly sized samples used in the previous studies. Furthermore, we adopt the velocity dispersion function of early-type galaxies based on SDSS DR1 and DR5. It turns out that the indexes of power-law model are consistent with the previous studies, thus our results indicate the mild evolution in the number and velocity dispersion of early-type galaxies out to z = 1. However we found that the values for p and q used by Mitchell et al. are inconsistent with the presently available observational data. More complete sample is necessary to withdraw more realistic determination on these parameters. △ Less

Submitted 3 December, 2007; originally announced December 2007.

Comments: Accepted to MNRAS, 7 pages, 3 figures

Journal ref: Mon.Not.Roy.Astron.Soc.384:843-848,2008

A Method for Weak Lensing Flexion Analysis by the HOLICs Moment Approach

Authors: Yuki Okura, Keiichi Umetsu, Toshifumi Futamase

Abstract: We have developed a method for measuring higher-order weak lensing distortions of faint background galaxies, namely the weak gravitational flexion, by fully extending the Kaiser, Squires & Broadhurst method to include higher-order lensing image characteristics (HOLICs) introduced by Okura, Umetsu, & Futamase. We take into account explicitly the weight function in calculations of noisy shape mome… ▽ More We have developed a method for measuring higher-order weak lensing distortions of faint background galaxies, namely the weak gravitational flexion, by fully extending the Kaiser, Squires & Broadhurst method to include higher-order lensing image characteristics (HOLICs) introduced by Okura, Umetsu, & Futamase. We take into account explicitly the weight function in calculations of noisy shape moments and the effect of higher-order PSF anisotropy, as well as isotropic PSF smearing. Our HOLICs formalism allows accurate measurements of flexion from practical observational data in the presence of non-circular, anisotropic PSF. We test our method using mock observations of simulated galaxy images and actual, ground-based Subaru observations of the massive galaxy cluster A1689 ($z=0.183$). From the high-precision measurements of spin-1 first flexion, we obtain a high-resolution mass map in the central region of A1689. The reconstructed mass map shows a bimodal feature in the central $4'\times 4'$ region of the cluster. The major, pronounced peak is associated with the brightest cluster galaxy and central cluster members, while the secondary mass peak is associated with a local concentration of bright galaxies. The refined, high-resolution mass map of A1689 demonstrates the power of the generalized weak lensing analysis techniques for quantitative and accurate measurements of the weak gravitational lensing signal. △ Less

Submitted 20 February, 2008; v1 submitted 11 October, 2007; originally announced October 2007.

Comments: Accepted for publication in ApJ, 36 pages, including 9 figures and 1 table. A version with higher resolution figures and Full Appendix can be found at http://www.asiaa.sinica.edu.tw/keiichi/OUF2/ouf2008.pdf

A Phantom does not Result from a Backreaction

Authors: Hisako Tanaka, Toshifumi Futamase

Abstract: The backreaction of nonlinear inhomogeneities to the cosmic expansion is analyzed in the framework of general relativity with a cosmological constant. By defining the spatially averaged matter energy density, we find that the cosmological constant induces a new type of backreaction whose equation of state is $p=-4/3 ρ$, where $ρ$ and $p$ are the effective energy density and effective pressure of… ▽ More The backreaction of nonlinear inhomogeneities to the cosmic expansion is analyzed in the framework of general relativity with a cosmological constant. By defining the spatially averaged matter energy density, we find that the cosmological constant induces a new type of backreaction whose equation of state is $p=-4/3 ρ$, where $ρ$ and $p$ are the effective energy density and effective pressure of the backreaction in the averaged Friedmann universe. However, the effective density is negative, and thus it decreases the acceleration caused by the cosmological constant. △ Less

Submitted 4 January, 2007; v1 submitted 6 December, 2006; originally announced December 2006.

Comments: 6pages, accepted by Progress of Theoretical Physics

Journal ref: Prog.Theor.Phys.117:183-188,2007

A New Measure for Weak Lensing Flexion

Authors: Yuki Okura, Keiichi Umetsu, Toshifumi Futamase

Abstract: We study a possibility to use the octopole moment of gravitationally lensed images as a direct measure of the third-order weak gravitational lensing effect, or the gravitational flexion. It turns out that there is a natural relation between flexion and certain combinations of octopole/higher-multipole moments which we call the Higher Order Lensing Image's Characteristics (HOLICs). This will allo… ▽ More We study a possibility to use the octopole moment of gravitationally lensed images as a direct measure of the third-order weak gravitational lensing effect, or the gravitational flexion. It turns out that there is a natural relation between flexion and certain combinations of octopole/higher-multipole moments which we call the Higher Order Lensing Image's Characteristics (HOLICs). This will allow one to measure directly flexion from observable octopole and higher-multipole moments of background images. We show based on simulated observations how the use of HOLICs can improve the accuracy and resolution of a reconstructed mass map, in which we assume Gaussian uncertainties in the shape measurements estimated using deep i'-band data of blank fields observed with Suprime-Cam on the Subaru telescope. △ Less

Submitted 24 January, 2007; v1 submitted 13 July, 2006; originally announced July 2006.

Comments: 18 pages, 2 figures, accepted for publication in ApJ

Journal ref: Astrophys.J.660:995-1002,2007

On the Equation of Motion for a Fast Moving Small Object in the Strong Field Point Particle Limit

Authors: Takashi Fukumoto, Toshifumi Futamase, Yousuke Itoh

Abstract: We give a straightforward and divergence free derivation of the equation of motion for a small but finite object in an arbitrary background using strong field point particle limit. The resulting equation becomes a generalized geodesic for a non-rotating spherical object which is consisitent with previous studies. We give a straightforward and divergence free derivation of the equation of motion for a small but finite object in an arbitrary background using strong field point particle limit. The resulting equation becomes a generalized geodesic for a non-rotating spherical object which is consisitent with previous studies. △ Less

Submitted 4 January, 2007; v1 submitted 27 June, 2006; originally announced June 2006.

Comments: accepted by Progress of Theoretical Physics

Journal ref: Prog.Theor.Phys. 116 (2006) 423-428

Toward a No-Go Theorem for an Accelerating Universe through a Nonlinear Backreaction

Authors: Masumi Kasai, Hideki Asada, Toshifumi Futamase

Abstract: The backreaction of nonlinear inhomogeneities to the cosmic expansion is re-analyzed in the framework of general relativity. Apparent discrepancies regarding the effect of the nonlinear backreaction, which exist among the results of previous works in different gauges, are resolved. By defining the spatially averaged matter energy density as a conserved quantity in the large comoving volume, it i… ▽ More The backreaction of nonlinear inhomogeneities to the cosmic expansion is re-analyzed in the framework of general relativity. Apparent discrepancies regarding the effect of the nonlinear backreaction, which exist among the results of previous works in different gauges, are resolved. By defining the spatially averaged matter energy density as a conserved quantity in the large comoving volume, it is shown that the nonlinear backreaction neither accelerates nor decelerates the cosmic expansion in a matter-dominated universe. The present result in the Newtonian gauge is consistent with the previous results obtained in the comoving synchronous gauge. Although our work does not give a complete proof, it strongly suggests the following no-go theorem: No cosmic acceleration occurs as a result of the nonlinear backreaction via averaging. △ Less

Submitted 29 March, 2006; v1 submitted 23 February, 2006; originally announced February 2006.

Comments: 6 pages (PTPTeX); accepted for publication in Prog. Theor. Phys; minor changes to match published version

Journal ref: Prog.Theor.Phys. 115 (2006) 827-832

Cosmology with High-redshift Galaxy Survey: Neutrino Mass and Inflation

Authors: Masahiro Takada, Eiichiro Komatsu, Toshifumi Futamase

Abstract: (abridged) High-z galaxy redshift surveys open up exciting possibilities for precision determinations of neutrino masses and inflationary models. The high-z surveys are more useful for cosmology than low-z ones owing to much weaker non-linearities in matter clustering, redshift-space distortion and galaxy bias. We can then utilize the two-dimensional information of the linear power spectrum in a… ▽ More (abridged) High-z galaxy redshift surveys open up exciting possibilities for precision determinations of neutrino masses and inflationary models. The high-z surveys are more useful for cosmology than low-z ones owing to much weaker non-linearities in matter clustering, redshift-space distortion and galaxy bias. We can then utilize the two-dimensional information of the linear power spectrum in angular and redshift space to measure the scale-dependent suppression of matter clustering due to neutrino free-streaming as well as the shape of the primordial power spectrum. To illustrate capabilities of high-z surveys for constraining neutrino masses and the primordial power spectrum, we compare three future redshift surveys covering 300 square degrees at 0.5<z<2, 2<z<4, and 3.5<z<6.5. We find that, combined with the cosmic microwave background data expected from the Planck satellite, these surveys allow precision determination of the total neutrino mass with the projected errors of sigma(m_nu)=0.059, 0.043, and 0.025 eV, respectively, thus yielding a positive detection of the neutrino mass rather than an upper limit, as sigma(m_nu) is smaller than the lower limits to the neutrino masses implied from the neutrino oscillation experiments. The accuracies of constraining the tilt and running index of the primordial power spectrum, sigma(n_s)=(3.8, 3.7, 3.0)x10^-3, and sigma(alpha_s)=(5.9, 5.7, 2.4)x10^-3, respectively, are smaller than the current uncertainties by more than an order of magnitude, which will allow us to discriminate between candidate inflationary models. In particular, the error on alpha_s from the highest redshift survey is not very far away from the prediction of a class of simple inflationary models driven by a massive scalar field with self-coupling, alpha_s=-(0.8-1.2)x10^-3. △ Less

Submitted 5 June, 2006; v1 submitted 15 December, 2005; originally announced December 2005.

Comments: 23 pages, 8 figures, to appear in PRD. Added Figures 5 and 8 and Table 7

Journal ref: Phys.Rev.D73:083520,2006

Discovery of a Strongly Lensed Galaxy at z=3.9 behind a z=0.83 Galaxy Cluster

Authors: Keiichi Umetsu, Masayuki Tanaka, Tadayuki Kodama, Ichi Tanaka, Toshifumi Futamase, Nobunari Kashikawa, Takako Hoshi

Abstract: We report the discovery and spectroscopic confirmation of three gravitationally-lensed images of a galaxy at z=3.9 in the background of a distant, rich cluster of galaxies at z=0.83, on the basis of observations with Faint Object Camera And Spectrograph (FOCAS) on the Subaru telescope. We construct a simple lens model of the cluster mass distribution based on Jee et al.'s weak lensing mass estim… ▽ More We report the discovery and spectroscopic confirmation of three gravitationally-lensed images of a galaxy at z=3.9 in the background of a distant, rich cluster of galaxies at z=0.83, on the basis of observations with Faint Object Camera And Spectrograph (FOCAS) on the Subaru telescope. We construct a simple lens model of the cluster mass distribution based on Jee et al.'s weak lensing mass estimates from deep, high-resolution images by Advanced Camera for Surveys (ACS) on the Hubble Space Telescope. The lens model can account simultaneously for the observed image configuration and the flux ratio of the closer pair located close to the critical curve. The parities of the three images are also consistent with the lensing hypothesis. Since this galaxy is apparently bright (i'_{AB} \sim 23.7) for its redshift due to the magnification, it serves as a good high redshift target on which we can make extensive and detailed studies based on multi-wavelength imaging and spectroscopy. △ Less

Submitted 30 September, 2005; v1 submitted 29 June, 2005; originally announced June 2005.

Comments: 4 pages, 3 figures, 1 table, accepted for publication in PASJ

Journal ref: PASJ, 57, 877 (2005)

The Surprisingly Steep Mass Profile of Abell 1689, from a Lensing Analysis of Subaru Images

Authors: Tom Broadhurst, Masahiro Takada, Keiichi Umetsu, Xu Kong, Nobuo Arimoto, Masashi Chiba, Toshifumi Futamase

Abstract: Subaru observations of A1689 (z=0.183) are used to derive an accurate, model-independent mass profile for the entire cluster, r<2 Mpc/h, by combining magnification bias and distortion measurements. The projected mass profile steepens quickly with increasing radius, falling away to zero at r~1.0 Mpc/h, well short of the anticipated virial radius. Our profile accurately matches onto the inner prof… ▽ More Subaru observations of A1689 (z=0.183) are used to derive an accurate, model-independent mass profile for the entire cluster, r<2 Mpc/h, by combining magnification bias and distortion measurements. The projected mass profile steepens quickly with increasing radius, falling away to zero at r~1.0 Mpc/h, well short of the anticipated virial radius. Our profile accurately matches onto the inner profile, r<200 kpc/h, derived from deep HST/ACS images. The combined ACS and Subaru information is well fitted by an NFW profile with virial mass, (1.93 \pm 0.20)10^15 M_sun, and surprisingly high concentration, c_vir=13.7^{+1.4}_{-1.1}, significantly larger than theoretically expected (c_vir~4), corresponding to a relatively steep overall profile. A slightly better fit is achieved with a steep power-law model that has its 2D logarithmic slope -3 and core radius theta_c~1.7' (r_c~210 kpc/h), whereas an isothermal profile is strongly rejected. These results are based on a reliable sample of background galaxies selected to be redder than the cluster E/S0 sequence. By including the faint blue galaxy population a much smaller distortion signal is found, demonstrating that blue cluster members significantly dilute the true signal for r~400 kpc/h. This contamination is likely to affect most weak lensing results to date. △ Less

Submitted 8 December, 2004; originally announced December 2004.

Comments: 5 pages, 3 figures, to appear in ApJL

Journal ref: Astrophys.J.619:L143,2005

The topology of weak lensing field in the neighborhood of MS1054-03

Authors: Jun'ichi Sato, Keiichi Umetsu, Toshifumi Futamase, Toru Yamada

Abstract: We report the first measurement of genus curves for the two-dimensional mass map in the neighborhood of rich, X-ray luminous galaxy cluster MS1054-03 at z=0.83, reconstructed from weak lensing data obtained by Suprime-Cam on the prime focus of 8.2m Subaru telescope. We find that the genus curve measured in the whole survey field deviates from that expected from a random Gaussian field. We show t… ▽ More We report the first measurement of genus curves for the two-dimensional mass map in the neighborhood of rich, X-ray luminous galaxy cluster MS1054-03 at z=0.83, reconstructed from weak lensing data obtained by Suprime-Cam on the prime focus of 8.2m Subaru telescope. We find that the genus curve measured in the whole survey field deviates from that expected from a random Gaussian field. We show that this non-Gaussianity is induced by the rich cluster in this region, and that the genus curve for the region without the cluster is consistent with the prediction for a random Gaussian field. The results demonstrate the feasibility of the genus statistics to examine the non-Gaussianity due to the large scale structures and to probe the statistical properties of the large scale clustering. △ Less

Submitted 8 November, 2002; v1 submitted 6 November, 2002; originally announced November 2002.

Comments: 13 pages, 4 figures, submitted to The Astrophysical Journal Letters

Journal ref: Astrophys.J.582:L67-L70,2003

Subaru Spectroscopy of the Gravitational Lens HST 14176+5226: Implications for a Large Cosmological Constan

Authors: Youichi Ohyama, Takashi Hamana, Nobunari Kashikawa, Masashi Chiba, Toshifumi Futamase, Masanori Iye, Koji S. Kawabata, Kentaro Aoki, Toshiyuki Sasaki, George Kosugi, Tadafumi Takata

Abstract: We present new optical spectroscopy of the lens elliptical galax in the ``Einstein Cross'' lens system HST 14176+5226, using the Faint Object Camera and Spectrograph (FOCAS) of the Subaru t Our spectroscopic observations are aimed at measuring the stella dispersion of the lens galaxy, located at high redshift of z_L= as an important component to lens models. We have measured this 230 +- 14 km s^… ▽ More We present new optical spectroscopy of the lens elliptical galax in the ``Einstein Cross'' lens system HST 14176+5226, using the Faint Object Camera and Spectrograph (FOCAS) of the Subaru t Our spectroscopic observations are aimed at measuring the stella dispersion of the lens galaxy, located at high redshift of z_L= as an important component to lens models. We have measured this 230 +- 14 km s^{-1} (1 sigma) inside 0.35 effective radi based on the comparison between the observed galaxy spectrum and templates of three G-K giants by means of the Fourier cross-corr To extract the significance of this information on the geometry universe which also affects the lensing of the background image, to fit three different lens models to the available data of the Provided that the lens galaxy has the structural and dynamical p (i.e., its radial density profile, core radius, and velocity ani similar to those of local elliptical galaxies, we calculate the function for the simultaneous reproduction of both the observed and newly measured velocity dispersion of the lens. Although the interval depends rather sensitively on the adopted lens models o parameters, our experiments suggest the larger likelihood for a cosmological constant, Omega_Lambda: formal 1 sigma lower Omega_Lambda in the flat universe ranges 0.73 to 0.97, where lower limit is basically unavailable. This method for determinin model is thus dependent on lens models but is insensitive to oth ambiguities, such as the dust absorption or the evolutionary eff galaxies. Exploring spectroscopic observations of more lens gala redshift may minimize the model uncertainties and thus place a m constraint on Omega_Lambda. △ Less

Submitted 19 March, 2002; originally announced March 2002.

Comments: 23 pages, 5 figures, accepted for publication in The Astronomical Journal