Source author record

Toshifumi Futamase

Toshifumi Futamase appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

40works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

40 published item(s)

preprint2022arXiv

On the gauge invariance of the locally averaged Friedmann Universe and the Hubble tension

The Hubble tension cast a blight on the standard cosmology. As a possible attitude to the problem, the local variation of the expansion rate in an inhomogeneous cosmology has been proposed where the spatial averaging over a finite domain was introduced in order to construct local Friedmann spacetime. However, it is not clear that the concept of the spatial averaging itself is gauge invariant or not. Namely the conclusion obtained by the averaging in a particular gauge is physically meaningful or not. In this paper we address this question, namely the gauge invariance of the spatial averaging in cosmology . We show that the answer is positive by studying the spatial average in the gauge-invariant cosmological perturbation theory.

preprint2021arXiv

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

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.

preprint2020arXiv

Current status of space gravitational wave antenna DECIGO and B-DECIGO

Deci-hertz Interferometer Gravitational Wave Observatory (DECIGO) is the future Japanese space mission with a frequency band of 0.1 Hz to 10 Hz. DECIGO aims at the detection of primordial gravitational waves, which could be produced during the inflationary period right after the birth of the universe. There are many other scientific objectives of DECIGO, including the direct measurement of the acceleration of the expansion of the universe, and reliable and accurate predictions of the timing and locations of neutron star/black hole binary coalescences. DECIGO consists of four clusters of observatories placed in the heliocentric orbit. Each cluster consists of three spacecraft, which form three Fabry-Perot Michelson interferometers with an arm length of 1,000 km. Three clusters of DECIGO will be placed far from each other, and the fourth cluster will be placed in the same position as one of the three clusters to obtain the correlation signals for the detection of the primordial gravitational waves. We plan to launch B-DECIGO, which is a scientific pathfinder of DECIGO, before DECIGO in the 2030s to demonstrate the technologies required for DECIGO, as well as to obtain fruitful scientific results to further expand the multi-messenger astronomy.

preprint2020arXiv

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

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.

preprint2016arXiv

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

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}]$.

preprint2016arXiv

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

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.

preprint2015arXiv

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

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.

preprint2015arXiv

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

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.

preprint2015arXiv

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

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.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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.

preprint2014arXiv

The magnitude-redshift relation in a realistic inhomogeneous universe

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.

preprint2013arXiv

A Systematic Review of Strong Gravitational Lens Modeling Software

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.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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.

preprint2013arXiv

δN formalism

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.

preprint2012arXiv

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

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.

preprint2012arXiv

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

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$.

preprint2012arXiv

New solutions of charged regular black holes and their stability

We construct new regular black hole solutions by matching the de Sitter solution and the Reissner-Nordstrom solution with a timelike thin shell. The thin shell is assumed to have mass but no pressure and obeys an equation of motion derived from Israel's junction conditions. By investigating the equation of motion for the shell, we obtain stationary solutions of charged regular black holes and examine stability of the solutions. Stationary solutions are found in limited ranges of 0.87L < m < 1.99L, and they are stable against small radial displacement of the shell with fixed values of m, M, and Q if M>0, where L is the de Sitter horizon radius, m the black hole mass, M the proper mass of the shell and Q the black hole charge. All the solutions obtained are highly charged in the sense of Q/m >0.866. By taking the massless limit of the shell in the present regular black hole solutions, we obtain the charged regular black hole with a massless shell obtained by Lemos and Zanchin and investigate stability of the solutions. It is found that Lemos and Zanchin's regular black hole solutions given by the massless limit of the present regular black hole solutions permit stable solutions, which are obtained by the limit of M -> 0.

preprint2012arXiv

Toward an accurate Dark Matter power spectrum beyond BAO scales

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).

preprint2011arXiv

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

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.

preprint2011arXiv

LoCuSS: Subaru Weak Lensing Study of 30 Galaxy Clusters

(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.

preprint2011arXiv

Non-Gaussianity Consistency Relation for Multi-field Inflation

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.

preprint2010arXiv

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

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.

preprint2010arXiv

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

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.

preprint2009arXiv

A method to measure a relative transverse velocity of source-lens-observer system using gravitational lensing of gravitational waves

Gravitational waves propagate along null geodesics like light rays in the geometrical optics approximation, and they may have a chance to suffer from gravitational lensing by intervening objects, as is the case for electromagnetic waves. Long wavelength of gravitational waves and compactness of possible sources may enable us to extract information in the interference among the lensed images. We point out that the interference term contains information of relative transverse velocity of the source-lens-observer system, which may be obtained by possible future space-borne gravitational wave detectors such as BBO/DECIGO.

preprint2007arXiv

A New Measure for Weak Lensing Flexion

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.

preprint2005arXiv

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

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.

preprint2004arXiv

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

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.

preprint2002arXiv

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

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.

preprint2001arXiv

Equation of motion for relativistic compact binaries with the strong field point particle limit : the second and half post-Newtonian order

We study the equation of motion appropriate to an inspiralling binary star system whose constituent stars have strong internal gravity. We use the post-Newtonian approximation with the strong field point particle limit by which we can introduce into general relativity a notion of a point-like particle with strong internal gravity without using Dirac delta distribution. Besides this limit, to deal with strong internal gravity we express the equation of motion in surface integral forms and calculate these integrals explicitly. As a result we obtain the equation of motion for a binary of compact bodies accurate through the second and half post-Newtonian (2.5 PN) order. This equation is derived in the harmonic coordinate. Our resulting equation perfectly agrees with Damour and Deruelle 2.5 PN equation of motion. Hence it is found that the 2.5 PN equation of motion is applicable to a relativistic compact binary.

preprint2001arXiv

On the abundance of collapsed objects

The redshift dependence of the abundance of collapsed objects places strong constraints on cosmological models of structure formation. We apply a recently proposed model describing the anisotropic collapse of inhomogeneous spatial domains. Compared with the spherical top-hat model, this generic model leads to significantly more collapsed objects at high redshifts: at redshift one and on the scale of rich clusters a factor of 65. Furthermore, for a fixed normalization of the initial fluctuation spectrum (σ_8=1), we predict four times as much presently collapsed objects on the mass-scale of rich clusters within the standard CDM cosmogony, compared to the spherical collapse.

preprint2001arXiv

R-mode Instability of Slowly Rotating Non-isentropic Relativistic Stars

We investigate properties of $r$-mode instability in slowly rotating relativistic polytropes. Inside the star slow rotation and low frequency formalism that was mainly developed by Kojima is employed to study axial oscillations restored by Coriolis force. At the stellar surface, in order to take account of gravitational radiation reaction effect, we use a near-zone boundary condition instead of the usually imposed boundary condition for asymptotically flat spacetime. Due to the boundary condition, complex frequencies whose imaginary part represents secular instability are obtained for discrete $r$-mode oscillations in some polytropic models. It is found that such discrete $r$-mode solutions can be obtained only for some restricted polytropic models. Basic properties of the solutions are similar to those obtained by imposing the boundary condition for asymptotically flat spacetime. Our results suggest that existence of a continuous part of spectrum cannot be avoided even when its frequency becomes complex due to the emission of gravitational radiation.

preprint2000arXiv

Cluster Mass Reconstruction by a Weak Shear Field

The tidal gravitational field of galaxy clusters causes a coherent distortion of the images of background sources. Since the distribution of image distortions, namely the shear field, traces the local gravitational potential of a deflector, it can be used to reconstruct the two-dimensional mass distribution of clusters of galaxies. Moreover, the shear field can provide unique information on the redshift distribution of high-redshift galaxies. In this review we summarize recently-developed parameter-free methods of cluster-mass reconstruction based on the shear field, and we apply a mass-reconstruction method to the cluster Abell 370 at redshift 0.375.

preprint2000arXiv

Equation of motion for relativistic compact binaries with the strong field point particle limit : Formulation, the first post-Newtonian and multipole terms

We derive the equation of motion for the relativistic compact binaries in the post-Newtonian approximation taking explicitly their strong internal gravity into account. For this purpose we adopt the method of the point particle limit where the equation of motion is expressed in terms of the surface integrals. We examine carefully the behavior of the surface integrals in the derivation. As a result, we obtain the Einstein-Infeld-Hoffman equation of motion at the first post-Newtonian (1PN) order, and a part of the 2PN order which depends on the quadrupole moments and the spins of component stars. Hence, it is found that the equation of motion in the post-Newtonian approximation is valid for the compact binaries by a suitable definition of the mass, spin and quadrupole moment.

preprint2000arXiv

Gravitational Lensing Effect on the Two-point Correlation of Hotspots in the Cosmic Microwave Background

We investigate the weak gravitational lensing effect due to the large-scale structure of the universe on two-point correlations of local maxima ({\em hotspots}) in the 2D sky map of the cosmic microwave background (CMB) anisotropy. According to the Gaussian random statistics as most inflationary scenarios predict, the hotspots are discretely distributed with some {\em characteristic} angular separations on the last scattering surface owing to oscillations of the CMB angular power spectrum. The weak lensing then causes pairs of hotspots which are separated with the characteristic scale to be observed with various separations. We found that the lensing fairly smoothes the oscillatory features of the two-point correlation function of hotspots. This indicates that the hotspots correlations can be a new statistical tool for measuring shape and normalization of the power spectrum of matter fluctuations from the lensing signatures.

preprint1999arXiv

How is the local-scale gravitational instability influenced by the surrounding large-scale structure formation?

We develop the formalism to investigate the relation between the evolution of the large-scale (quasi) linear structure and that of the small-scale nonlinear structure in Newtonian cosmology within the Lagrangian framework. In doing so, we first derive the standard Friedmann expansion law using the averaging procedure over the present horizon scale. Then the large-scale (quasi) linear flow is defined by averaging the full trajectory field over a large-scale domain, but much smaller than the horizon scale. The rest of the full trajectory field is supposed to describe small-scale nonlinear dynamics. We obtain the evolution equations for the large-scale and small-scale parts of the trajectory field. These are coupled to each other in most general situations. It is shown that if the shear deformation of fluid elements is ignored in the averaged large-scale dynamics, the small-scale dynamics is described by Newtonian dynamics in an effective Friedmann-Robertson-Walker (FRW) background with a local scale factor. The local scale factor is defined by the sum of the global scale factor and the expansion deformation of the averaged large-scale displacement field. This means that the evolution of small-scale fluctuations is influenced by the surrounding large-scale structure through the modification of FRW scale factor. The effect might play an important role in the structure formation scenario. Furthermore, it is argued that the so-called {\it optimized} or {\it truncated} Lagrangian perturbation theory is a good approximation in investigating the large-scale structure formation up to the quasi nonlinear regime, even when the small-scale fluctuations are in the non-linear regime.

preprint1997arXiv

Constraints on the chaotic inflationary scenario with a nonminimally coupled "inflaton" field from the cosmic microwave background radiation anisotropy

We investigate a possibility to restrict the chaotic inflationary scenario with "large" nonminimally coupled inflaton("phi") field considered by Fakir and Unruh by means of the observed CMBR anisotropy. We calculate the contributions of the long wavelength gravitational waves generated in the inflationary period to CMBR anisotropy quadrupole moment. We obtain the constraint lambda/(xi^2) = 1.8 X 10^(-8)/psi_i, where "lambda" is the self coupling, "xi" is the nonminimal coupling constant and "psi_i" means the initial value of psi = 8piG(phi^2). Combining this with previously obtained constraint that 2.1 times squre root of lambda/(xi^2) = dT/T = 1.1 X 10^(-5), we conclude that the initial value has to be psi_i = 1.6 X 10^3. If the self-coupling has a reasonably values of order 10^(-2), then xi = 10^4 and phi_i = 10^(-1) times Planck mass.