Researcher profile

Masahiro Takada

Masahiro Takada contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
26works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

26 published item(s)

preprint2024arXiv

The Milky Way Tomography with Subaru Hyper Suprime-Cam. I. Halo substructures

We analyze the photometric data in the Wide layer of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) over $\sim 1,200$ deg$^{2}$ to uncover new halo substructures beyond the distance, $D_{\odot}\sim$ 30 kpc, from the Sun. For this purpose, we develop an isochrone filter for an old, metal-poor stellar system to extract the faint main-sequence stars at a range of distances. With this method, we detect, not only the previously discovered substructures such as the Orphan Stream, but also the new overdensity toward Boötes at about $D_{\odot}\sim$ 60 kpc and the new stream-like feature toward Pisces at around $D_{\odot}\sim$ 60 kpc. It has been suggested that a small-scale overdensity exists in this direction of Pisces (the so-called Pisces Overdensity), but our results show that the overdensity is widely spread with a tidally elongated feature. Combining our results with the ongoing Hyper Suprime-Cam narrow-band survey and the near-future spectroscopic survey with Prime Focus Spectrograph (PFS) will allow us to place strong constraints on the origin of these halo substructures.

preprint2022arXiv

Analysis method for 3D power spectrum of projected tensor field with fast estimator and window convolution modelling: an application to intrinsic alignments

Rank-2 tensor fields of large-scale structure, e.g. a tensor field inferred from shapes of galaxies, open up a window to directly access 2-scalar, 2-vector and 2-tensor modes, where the scalar fields can be measured independently from the standard density field that is traced by distribution of galaxies. Here we develop an estimator of the multipole moments of power spectra for the three-dimensional tensor field, taking into account the projection onto plane perpendicular to the line-of-sight direction. To do this, we find that a convenient representation of the power spectrum multipoles can be obtained by the use of the associated Legendre polynomials in the form which allows for the fast Fourier transform estimations under the local plane-parallel (LPP) approximation. The formulation also allows us to obtain the Hankel transforms to connect the two-point statistics in Fourier and configuration space, which are needed to derive theoretical templates of the power spectrum including convolution of a survey window. To validate our estimators, we use the simulation data of the projected tidal field assuming a survey window that mimics the BOSS-like survey footprint. We show that the LPP estimators fairly well recover the multipole moments that are inferred from the global plane-parallel approximation. We find that the survey window causes a more significant change in the multipole moments of projected tensor power spectrum at $k\lesssim 0.1\,h{\rm Mpc}^{-1}$ from the input power spectrum, than in the density power spectrum. Nevertheless, our method to compute the theory template including the survey window effects successfully reproduces the window-convolved multipole moments measured from the simulations. The analysis method presented here paves the way for a cosmological analysis using three-dimensional tensor-type tracers of large-scale structure for current and future surveys.

preprint2022arXiv

Constraining Dark matter annihilation with Dark Energy Survey Y3 LSBG sample

To reveal natures of the dark matter (DM) particles, a gamma-ray signal produced in annihilation processes of DM into the standard model particles has been one of the major probes. The cross-correlation between highly DM dominated structures, such as local dwarf galaxies, and observed photons in the direction of the structures, has been explored and provided stringent constraints on the annihilation rate. In our previous work, we have shown that it is sufficient to know the distance distribution of the galaxy sample and individual distance measurement is not required for constraining the annihilation rate. In this work, we apply the method to low-surface brightness galaxies (LSBGs) with unknown individual redshifts provided from the Dark Energy Survey (DES) Year 3 data. With all DES-LSBGs of about 24,000 objects, we find that the upper limits of the cross section for bb channel with 95% C.L. is 3*10^-25 cm^3/s at DM mass of 100 GeV. To be more conservative, we remove about 7000 LSBGs within 1 degree from resolved gamma-ray point sources and then the constraint becomes 30 % weaker than the one with all samples in all DM mass ranges.

preprint2022arXiv

Cosmological constraints without fingers of God

Non-linear redshift-space distortions ("fingers of God") are challenging to model analytically, a fact that limits the applicability of perturbation theory in redshift space as compared to real space. We show how this problem can be mitigated using a new observable, $Q_0$, which can be easily estimated from the redshift space clustering data and is approximately equal to the real space power spectrum. The new statistic does not suffer from fingers of God and can be accurately described with perturbation theory down to $k_{\rm max}\simeq 0.4~h~\text{Mpc}^{-1}$. It can be straightforwardly included in the likelihood at negligible additional computational cost, and yields noticeable improvements on cosmological parameters compared to standard power spectrum multipole analyses. Using both simulations and observational data from the Baryon Oscillation Spectroscopic Survey, we show that improvements vary from $10\%$ to $100\%$ depending on the cosmological parameter considered, the galaxy sample and the survey volume.

preprint2022arXiv

Dark matter constraints with stacked gamma rays scales with the number of galaxies

Low-surface-brightness galaxies (LSBGs) are interesting targets for searches of dark matter emission due to their low baryonic content. However, predicting their expected dark matter emissivities is difficult because of observational challenges in their distance measurements. Here we present a stacking method that makes use of catalogs of LSBGs and maps of unresolved gamma-ray emission measured by the Fermi Gamma-Ray Space Telescope. We show that, for relatively large number of LSBGs, individual distance measurements to the LSBGs are not necessary, instead the overall distance distribution of the population is sufficient in order to impose dark matter constraints. Further, we demonstrate that the effect of the covariance between two galaxies located closely -- at an angular distance comparable to the size of the Fermi point spread function -- is negligibly small. As a case in point, we apply our pipeline to a sample of 800 faint LSBGs discovered by Hyper Suprime-Cam and find that, the 95 per cent confidence level upper limits on the dark matter annihilation cross-section scales with inverse of the number LSBGs. In light of this linear dependence with the number of objects, we argue this methodology could provide extremely powerful limits if it is applied to the more than 105 LSBGs readily available with the Legacy Survey of Space and Time.

preprint2022arXiv

Full-shape cosmology analysis of SDSS-III BOSS galaxy power spectrum using emulator-based halo model: a $5\%$ determination of $σ_8$

We present the results obtained from the full-shape cosmology analysis of the redshift-space power spectra for 4 galaxy samples of the SDSS-III BOSS DR12 galaxy catalog over $0.2 < z < 0.75$. For the theoretical template, we use an emulator that was built from an ensemble set of $N$-body simulations, which enables fast and accurate computation of the redshift-space power spectrum of halos. Combining with the halo occupation distribution to model the halo-galaxy connection, we can compute the redshift-space power spectrum of BOSS-like galaxies in the flat $Λ$CDM cosmology. In our cosmology inference, we use the power spectrum monopole, quadrupole and hexadecapole and include 7 nuisance parameters to model uncertainties in the halo-galaxy connection for each galaxy sample, but do not use any information on the abundance of galaxies. We demonstrate a validation of our analysis pipeline using the mock catalogs of BOSS-like galaxies, generated using different recipes of the halo-galaxy connection and including the assembly bias effect. Assuming weak priors on cosmological parameters, except for $Ω_{\rm b}h^2$ and $n_{\rm s}$, we show that our model well reproduces the BOSS power spectra. Including the power spectrum information up to $k_{\rm max}=0.25\,h{\rm Mpc}^{-1}$, we find $Ω_\mathrm{m}=0.301^{+0.012}_{-0.011}$, $H_0=68.2 \pm 1.4~\mathrm{km\,s}^{-1}\mathrm{Mpc}^{-1}$, and $σ_8=0.786^{+0.036}_{-0.037}$, for the mode and 68\% credible interval, after marginalization over nuisance parameters. We find little improvement in the cosmological parameters beyond a maximum wavelength $k_{\rm max}\simeq 0.2\,h\,{\rm Mpc}^{-1}$ due to the shot noise domination and marginalization of the halo-galaxy connection parameters. Our results are consistent with the Planck CMB results within $1σ$ statistical uncertainties.

preprint2022arXiv

Mock catalogues of emission line galaxies based on the local mass density in dark-matter only simulations

The high-precision measurement of spatial clustering of emission line galaxies (ELGs) is a primary objective for upcoming cosmological spectroscopic surveys. The source of strong emission of ELGs is nebular emission from surrounding ionized gas irradiated by massive short-lived stars in star-forming galaxies. As a result, ELGs are more likely to reside in newly-formed halos and this leads to a nonlinear relation between ELG number density and matter density fields. In order to estimate the covariance matrix of cosmological observables, it is essential to produce many independent realisations to simulate ELG distributions for large survey volumes. To this end, we present a novel and fast scheme to populate ELGs in dark-matter only $N$-body simulations based on local density field. This method enables fast production of mock ELG catalogues suitable for verifying analysis methods and quantifying observational systematics in upcoming spectroscopic surveys and can populate ELGs in moderately high-density regions even though the halo structure cannot be resolved due to low resolution. The power spectrum of simulated ELGs is consistent with results of hydrodynamical simulations up to fairly small scales ($\lesssim 1 h \, \mathrm{Mpc}^{-1}$), and the simulated ELGs are more likely to be found in filamentary structures, which is consistent with results of semi-analytic and hydrodynamical simulations. Furthermore, we address the redshift-space power spectrum of simulated ELGs. The measured multipole moments of simulated ELGs clearly exhibit a weaker Finger-of-God effect than those of matter due to infalling motions towards halo centre, rather than random virial motions inside halos.

preprint2022arXiv

The Subaru HSC weak lensing mass-observable scaling relations of spectroscopic galaxy groups from the GAMA survey

We utilize the galaxy shape catalogue from the first-year data release of the Subaru Hyper Suprime-cam Survey (HSC) to study the dark matter content of galaxy groups in the Universe using weak lensing. We use galaxy groups from the Galaxy Mass and Assembly galaxy survey in approximately $100$ sq. degrees of the sky that overlap with the HSC survey as lenses. We restrict our analysis to the $1587$ groups with at least five members. We divide these groups into six bins each of group luminosity and group member velocity dispersion and measure the lensing signal with a signal-to-noise ratio of $55$ and $51$ for these two different selections, respectively. We use a Bayesian halo model framework to infer the halo mass distribution of our groups binned in the two different observable properties and constrain the power-law scaling relation, and the scatter between mean halo masses and the two group observable properties. We obtain a 5 percent constraint on the amplitude of the scaling relation between halo mass and group luminosity with $\avg{M} = (0.81\pm 0.04)\times10^{14}\hinvMsun$ for $L_{\rm grp}=10^{11.5}\hinvsqLsun$, and a power-law index of $α=1.01\pm 0.07$. We constrain the amplitude of the scaling relation between halo mass and velocity dispersion to be $\avg{M}=(0.93\pm 0.05)\times10^{14}\hinvMsun$ for $σ=500 \kms$ and a power-law index to be $α=1.52\pm0.10$. However, these scaling relations are sensitive to the exact cuts applied to the number of group members. Comparisons with similar scaling relations from the literature show that our results are consistent and have significantly reduced errors.

preprint2022arXiv

The three-year shear catalog of the Subaru Hyper Suprime-Cam SSP Survey

We present the galaxy shear catalog that will be used for the three-year cosmological weak gravitational lensing analyses using data from the Wide layer of the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) Survey. The galaxy shapes are measured from the $i$-band imaging data acquired from 2014 to 2019 and calibrated with image simulations that resemble the observing conditions of the survey based on training galaxy images from the Hubble Space Telescope in the COSMOS region. The catalog covers an area of 433.48 deg$^2$ of the northern sky, split into six fields. The mean $i$-band seeing is 0.59 arcsec. With conservative galaxy selection criteria (e.g., $i$-band magnitude brighter than 24.5), the observed raw galaxy number density is 22.9 arcmin$^{-2}$, and the effective galaxy number density is 19.9 arcmin$^{-2}$. The calibration removes the galaxy property-dependent shear estimation bias to a level: $|δm|<9\times 10^{-3}$. The bias residual $δm$ shows no dependence on redshift in the range $0<z\leq 3$. We define the requirements for cosmological weak lensing science for this shear catalog, and quantify potential systematics in the catalog using a series of internal null tests for systematics related to point-spread function modelling and shear estimation. A variety of the null tests are statistically consistent with zero or within requirements, but (i) there is evidence for PSF model shape residual correlations; and (ii) star-galaxy shape correlations reveal additive systematics. Both effects become significant on $>1$ degree scales and will require mitigation during the inference of cosmological parameters using cosmic shear measurements.

preprint2021arXiv

HSC Year 1 cosmology results with the minimal bias method: HSC$\times$BOSS galaxy-galaxy weak lensing and BOSS galaxy clustering

We present cosmological parameter constraints from a blinded joint analysis of galaxy-galaxy weak lensing, $Δ\!Σ(R)$, and projected correlation function, $w_\mathrm{p}(R)$, measured from the first-year HSC (HSC-Y1) data and SDSS spectroscopic galaxies over $0.15<z<0.7$. We use luminosity-limited samples as lens samples for $Δ\!Σ$ and as large-scale structure tracers for $w_\mathrm{p}$ in three redshift bins, and use the HSC-Y1 galaxy catalog to define a secure sample of source galaxies at $z_\mathrm{ph}>0.75$ for the $Δ\!Σ$ measurements, selected based on their photometric redshifts. For theoretical template, we use the &#34;minimal bias&#34; model for the cosmological clustering observables for the flat $Λ$CDM cosmological model. We compare the model predictions with the measurements in each redshift bin on large scales, $R>12$ and $8~h^{-1}\mathrm{Mpc}$ for $Δ\!Σ(R)$ and $w_\mathrm{p}(R)$, respectively, where the perturbation theory-inspired model is valid. When we employ weak priors on cosmological parameters, without CMB information, we find $S_8=0.936^{+0.092}_{-0.086}$, $σ_8=0.85^{+0.16}_{-0.11}$, and $Ω_\mathrm{m}=0.283^{+0.12}_{-0.035}$ for the flat $Λ$CDM model. Although the central value of $S_8$ appears to be larger than those inferred from other cosmological experiments, we find that the difference is consistent with expected differences due to sample variance, and our results are consistent with the other results to within the statistical uncertainties. (abriged)

preprint2021arXiv

Precision analysis of the redshift-space galaxy bispectrum

We study the information content of the angle-averaged (monopole) redshift space galaxy bispectrum. The main novelty of our approach is the use of a systematic tree-level perturbation theory model that includes galaxy bias, IR resummation, and also accounts for nonlinear redshift space distortions, binning, and projection effects. We analyze data from the PT challenge simulations, whose cumulative volume of 566 $h^{-3}$Gpc$^3$ allows for a precise comparison to theoretical predictions. Fitting the power spectrum and bispectrum of our simulated data, and varying all necessary cosmological and nuisance parameters in a consistent Markov chain Monte Carlo analysis, we find that our tree-level bispectrum model is valid up to $k_{\max}=0.08~h{\rm Mpc}^{-1}$ (at $z=0.61$). We also find that inclusion of the bispectrum monopole improves constraints on cosmological parameters by $(5-15)\%$ relative to the power spectrum. The improvement is more significant for the quadratic bias parameters of our simulated galaxies, which we also show to deviate from biases of the host dark matter halos at the $\sim 3σ$ level. Finally, we adjust the covariance and scale cuts to match the volume of the BOSS survey, and estimate that within the minimal $Λ$CDM model the bispectrum data can tighten the constraint on the mass fluctuation amplitude $σ_8$ by roughly $10\%$.

preprint2021arXiv

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). XVI. 69 New Quasars at 5.8 < z < 7.0

We present the spectroscopic discovery of 69 quasars at 5.8 < z < 7.0, drawn from the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) imaging survey data. This is the 16th publication from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, and completes identification of all but the faintest candidates (i.e., i-band dropouts with zAB < 24 and y-band detections, and z-band dropouts with yAB < 24) with Bayesian quasar probability Pq > 0.1 in the HSC-SSP third public data release (PDR3). The sample reported here also includes three quasars with Pq < 0.1 at z ~ 6.6, which we selected in an effort to completely cover the reddest point sources with simple color cuts. The number of high-z quasars discovered in SHELLQs has now grown to 162, including 23 type-II quasar candidates. This paper also presents identification of seven galaxies at 5.6 < z < 6.7, an [O III] emitter at z = 0.954, and 31 Galactic cool stars and brown dwarfs. High-z quasars and galaxies comprise 75 % and 16 % respectively of all the spectroscopic SHELLQs objects that pass our latest selection algorithm with the PDR3 photometry. That is, a total of 91 % of the objects lie at z > 5.6. This demonstrates that the algorithm has very high efficiency, even though we are probing an unprecedentedly low-luminosity population down to M1450 ~ -21 mag.

preprint2021arXiv

Super-sample covariance of the thermal Sunyaev-Zel&#39;dovich effect

The thermal Sunyaev-Zel&#39;dovich (tSZ) effect is a powerful probe of cosmology. The statistical errors in the tSZ power spectrum measurements are dominated by the presence of massive clusters in a survey volume that are easy to identify on individual cluster basis. First, we study the impact of super sample covariance (SSC) on the tSZ power spectrum measurements, and find that the sample variance is dominated by the connected non-Gaussian (cNG) covariance arising mainly from Poisson number fluctuations of massive clusters in the survey volume. Second, we find that removing such individually-detected, massive clusters from the analysis significantly reduces the cNG contribution, thereby leading the SSC to be a leading source of the sample variance. We then show, based on Fisher analysis, that the power spectrum measured from the remaining diffuse tSZ effects can be used to obtain tight constraints on cosmological parameters as well as the hydrostatic mass bias parameter. Our method offers complementary use of individual tSZ cluster counts and the power spectrum measurements of diffuse tSZ signals for cosmology and intracluster gas physics.

preprint2021arXiv

Testing Stochastic Gravitational Wave Signals from Primordial Black Holes with Optical Telescopes

Primordial black holes (PBHs) can constitute the predominant fraction of dark matter (DM) if PBHs reside in the currently unconstrained &#34;sublunar&#34; mass range. PBHs originating from scalar perturbations generated during inflation can naturally appear with a broad spectrum in a class of models. The resulting stochastic gravitational wave (GW) background generated from such PBH production can account for the recently reported North American Nanohertz Observatory for Gravitational Waves (NANOGrav) pulsar timing array data signal, and will be testable in future GW observations by interferometer-type experiments such as Laser Interferometer Space Antenna (LISA). We show that the broad mass function of such PBH DM is already being probed by Subaru Hyper Suprime-Cam (HSC) microlensing data and is consistent with a detected candidate event. Upcoming observations of HSC will be able to provide an independent definitive test of the stochastic GW signals originating from such PBH DM production scenarios.

preprint2021arXiv

Third Data Release of the Hyper Suprime-Cam Subaru Strategic Program

The paper presents the third data release of Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP), a wide-field multi-band imaging survey with the Subaru 8.2m telescope. HSC-SSP has three survey layers (Wide, Deep, and UltraDeep) with different area coverages and depths, designed to address a wide array of astrophysical questions. This third release from HSC-SSP includes data from 278 nights of observing time and covers about 670 square degrees in all five broad-band filters at the full depth ($\sim26$~mag at $5σ$) in the Wide layer. If we include partially observed area, the release covers 1,470 square degrees. The Deep and UltraDeep layers have $\sim80\%$ of the originally planned integration times, and are considered done, as we have slightly changed the observing strategy in order to compensate for various time losses. There are a number of updates in the image processing pipeline. Of particular importance is the change in the sky subtraction algorithm; we subtract the sky on small scales before the detection and measurement stages, which has significantly reduced false detections. Thanks to this and other updates, the overall quality of the processed data has improved since the previous release. However, there are limitations in the data (for example, the pipeline is not optimized for crowded fields), and we encourage the user to check the quality assurance plots as well as a list of known issues before exploiting the data. The data release website is https://hsc-release.mtk.nao.ac.jp/.

preprint2020arXiv

Anisotropic separate universe simulations

The long-wavelength coherent overdensity and tidal force, which are not direct observables for a finite-volume survey, affect time evolution of cosmic structure formation and therefore clustering observables through the mode coupling. In this paper we develop an &#34;anisotropic&#34; separate universe (SU) simulation technique to simulate large-scale structure formation taking into account the effect of large-scale tidal force into the anisotropic expansion of local background. We modify the TreePM N-body simulation code to implement the anisotropic SU simulations, and then study the &#34;response&#34; function of matter power spectrum that describes how the matter power spectrum responds to the large-scale tidal effect as a function of wavenumber and redshift for a given global cosmology. We test and validate the SU simulation results from the comparison with the perturbation theory predictions and the results from high-resolution PM simulations. We find that the response function displays characteristic scale dependences over the range of scales down to nonlinear scales, up to k ~ 6 h/Mpc.

preprint2020arXiv

Blinded challenge for precision cosmology with large-scale structure: results from effective field theory for the redshift-space galaxy power spectrum

An accurate theoretical template for the galaxy power spectrum is a key for the success of ongoing and future spectroscopic surveys. We examine to what extent the Effective Field Theory of Large Scale Structure is able to provide such a template and correctly estimate cosmological parameters. To that end, we initiate a blinded challenge to infer cosmological parameters from the redshift-space power spectrum of high-resolution mock catalogs mimicking the BOSS galaxy sample but covering a hundred times larger cumulative volume. This gigantic simulation volume allows us to separate systematic bias due to theoretical modeling from the statistical error due to sample variance. The challenge task was to measure three unknown input parameters used in the simulation: the Hubble constant, the matter density fraction, and the clustering amplitude. We present analyses done by two independent teams, who have fitted the mock simulation data generated by yet another independent group. This allows us to avoid any confirmation bias by analyzers and pin down possible tuning of the specific EFT implementations. Both independent teams have recovered the true values of the input parameters within sub-percent statistical errors corresponding to the total simulation volume.

preprint2020arXiv

Constraining dark matter annihilation with HSC Low Surface Brightness Galaxies

Searches for dark matter annihilation signals have been carried out in a number of target regions such as the Galactic Center and Milky Way dwarf spheroidal galaxies (dSphs), among a few others. Here we propose low surface brightness galaxies (LSBGs) asnovel targets for the indirect detection of dark matter emission. In particular, LSBGs are known to have very large dark matter contents and be less contaminated by extragalactic gamma-ray sources (e.g., blazars) compared to star forming galaxies. We report on an analysis that uses eight LSBGs (detected by Subaru Hyper Suprime-Cam survey data) with known redshifts to conduct a search for gamma-ray emission at the positions of these new objects in Fermi Large Area Telescope data. We found no excesses of gamma-ray emission and set constraints on the dark matter annihilation cross-section. We exclude (at the 95% C.L.) dark matter scenarios predicting a cross-section higher than 10^-23[cm^3/s] for dark matter particles of mass 10 GeV self-annihilating in the b_b channel. Although this constraint is weaker than the ones reported in recent studies using other targets, we note that in the near future, the number of detections of new LSBGs will increase by a few orders of magnitude. We forecast that with the use of the full catalog of soon-to-be-detected LSBGs the constraint will reach cross-section sensitivities of ~ 3*10^-25 [cm^3/s] for dark matter particles with masses less than 10 GeV.

preprint2020arXiv

Cosmological information content in redshift-space power spectrum of SDSS-like galaxies in the quasi-nonlinear regime up to $k=0.3\,h\,$Mpc$^{-1}$

Clustering properties and peculiar velocities of halos in large-scale structure carry a wealth of cosmological information over a wide range of scales from linear to nonlinear scales. We use halo catalogs in a suite of high-resolution $N$-body simulations to construct mock catalogs of galaxies that resemble the SDSS-like luminous early-type galaxies at three redshift bins in the range $0.15\le z\le 0.7$. To do this we include 10 nuisance parameters to model variations in halo-galaxy connections for each redshift bin; the halo occupation distribution, and the spatial and velocity distributions of galaxies in the host halos. We evaluate the Fisher information matrix for the redshift-space power spectrum of SDSS-like galaxies using different sets of the mock catalogs that are generated from changes in each of model parameters; cosmological parameters ($σ_8$ and $Ω_{\rm m}$), the halo-galaxy connection parameters, and the cosmological distances ($D_{\rm A}$ and $H$ parameters at each redshift bin) for modeling an apparent geometrical distortion of the redshift-space power spectrum (the Alcock-Paczynski effect). We show that combining the monopole and quadrupole power spectra of galaxies allows for accurate estimations of the cosmological parameters and the cosmological distances, even after marginalization over the halo-galaxy parameters, by lifting the parameter degeneracies that are otherwise inevitable if either of the two spectra alone is used. When including the galaxy power spectrum information up to $k=0.3\,h\,$Mpc$^{-1}$, we find about factor of 6 gain in the cosmological information content of ($σ_8,Ω_{\rm m},D_{\rm A}$&#39;s and $H$&#39;s) compared to $k=0.2\,h\,$Mpc$^{-1}$.

preprint2020arXiv

Hunting gravitational wave black holes with microlensing

Gravitational microlensing is a powerful tool to search for a population of invisible black holes (BHs) in the Milky Way (MW), including isolated BHs and binary BHs at wide orbits that are complementary to gravitational wave observations. By monitoring highly populated regions of source stars like the MW bulge region, one can pursue microlensing events due to these BHs. We find that if BHs have a Salpeter-like mass function extended beyond $30M_\odot$ and a similar velocity and spatial structure to stars in the Galactic bulge and disk regions, the BH population is a dominant source of the microlensing events at long timescales of the microlensing light curve $\gtrsim 100~$days. This is due to a boosted sensitivity of the microlensing event rate to lens mass, given as $M^2$, for such long-timescale events. A monitoring observation of $2 \times 10^{10}$ stars in the bulge region over 10 years with the Rubin Observatory Legacy Survey of Space and Time (LSST) would enable one to find about $6\times 10^5$ BH microlensing events. We evaluate the efficiency of potential LSST cadences for characterizing the light curves of BH microlensing and find that nearly all events of long timescales can be detected.

preprint2020arXiv

On the wave optics effect on primordial black hole constraints from optical microlensing search

Microlensing of stars, e.g. in the Galactic bulge and Andromeda galaxy (M31), is among the most robust, powerful method to constrain primordial black holes (PBHs) that are a viable candidate of dark matter. If PBHs are in the mass range $M_{\rm PBH} \lower.5ex\hbox{$\; \buildrel < \over \sim \;$} 10^{-10}M_\odot$, its Schwarzschild radius ($r_{\rm Sch}$) becomes comparable with or shorter than optical wavelength ($λ)$ used in a microlensing search, and in this regime the wave optics effect on microlensing needs to be taken into account. For a lensing PBH with mass satisfying $r_{\rm Sch}\sim λ$, it causes a characteristic oscillatory feature in the microlensing light curve, and it will give a smoking gun evidence of PBH if detected, because any astrophysical object cannot have such a tiny Schwarzschild radius. Even in a statistical study, e.g. constraining the abundance of PBHs from a systematic search of microlensing events for a sample of many source stars, the wave effect needs to be taken into account. We examine the impact of wave effect on the PBH constraints obtained from the $r$-band (6210Å) monitoring observation of M31 stars in Niikura et al. (2019), and find that a finite source size effect is dominant over the wave effect for PBHs in the mass range $M_{\rm PBH}\simeq[10^{-11},10^{-10}]M_\odot$. We also discuss that, if a denser-cadence (10~sec), $g$-band monitoring observation for a sample of white dwarfs over a year timescale is available, it would allow one to explore the wave optics effect on microlensing light curve, if it occurs, or improve the PBH constraints in $M_{\rm PBH}\lower.5ex\hbox{$\; \buildrel < \over \sim \;$} 10^{-11}M_\odot$ even from a null detection.

preprint2020arXiv

Proper motion measurements for stars up to $100$ kpc with Subaru HSC and SDSS Stripe 82

We present proper motion measurements for more than $0.55$ million main-sequence stars, by comparing astrometric positions of matched stars between the multi-band imaging datasets from the Hyper Suprime-Cam (HSC) Survey and the SDSS Stripe 82. In doing this we use $3$ million galaxies to recalibrate the astrometry and set up a common reference frame between the two catalogues. The exquisite depth and the nearly $12$ years of time baseline between HSC and SDSS enable high-precision measurements of statistical proper motions for stars down to $i\simeq 24$. A validation of our method is demonstrated by the agreement with the $Gaia$ proper motions, to the precision better than $0.1$ mas yr$^{-1}$. To retain the precision, we make a correction of the subtle effects due to the differential chromatic refraction in the SDSS images based on the comparison with the $Gaia$ proper motions against colour of stars, which is validated using the SDSS spectroscopic quasars. Combining with the photometric distance estimates for individual stars based on the precise HSC photometry, we show a significant detection of the net proper motions for stars in each bin of distance out to $100$ kpc. The two-component tangential velocities after subtracting the apparent motions due to our own motion display rich phase-space structures including a clear signature of the Sagittarius stream in the halo region of distance range $[10,\ 35]$ kpc. We also measure the tangential velocity dispersion in the distance range $5-20$ kpc and find that the data are consistent with a constant isotropic dispersion of $80\pm 10\ {\rm km/s}$. More distant stars appear to have random motions with respect to the Galactic centre on average.

preprint2020arXiv

The impact of projection effects on cluster observables: stacked lensing and projected clustering

An optical cluster finder inevitably suffers from projection effects, where it misidentifies a superposition of galaxies in multiple halos along the line-of-sight as a single cluster. Using mock cluster catalogs built from cosmological N-body simulations, we quantify the impact of these projection effects with a particular focus on the observables of interest for cluster cosmology, namely the cluster lensing and the cluster clustering signals. We find that &#34;observed&#34; clusters, i.e. clusters identified by our cluster finder algorithm, exhibit lensing and clustering signals that deviate from expectations based on a statistically isotropic halo model -- while both signals agree with halo model expectations on small scales, they show unexpected boosts on large scales, by up to a factor of 1.2 or 1.4 respectively. We identify the origin of these boosts as the inherent selection bias of optical cluster finders for clusters embedded within filaments aligned with the line-of-sight, and show that a minority ($\sim 30\%$) of such clusters within the entire sample is responsible for this observed boost. We discuss the implications of our results on previous studies of optical cluster, as well as prospects for identifying and mitigating projection effects in future cluster cosmology analyses.

preprint2020arXiv

Towards a non-Gaussian model of redshift space distortions

To understand the nature of the accelerated expansion of the Universe, we need to combine constraints on the expansion rate and growth of structure. The growth rate is usually extracted from three dimensional galaxy maps by exploiting the effects of peculiar motions on galaxy clustering. However, theoretical models of the probability distribution function (PDF) of galaxy pairwise peculiar velocities are not accurate enough on small scales to reduce the error on theoretical predictions to the level required to match the precision expected for measurements from future surveys. Here, we improve the modelling of the pairwise velocity distribution by using the Skew-T PDF, which has nonzero skewness and kurtosis. Our model accurately reproduces the redshift-space multipoles (monopole, quadrupole and hexadecapole) predicted by N-body simulations, above scales of about $10\,h^{-1}{\rm Mpc}$. We illustrate how a Taylor expansion of the streaming model can reveal the contributions of the different moments to the clustering multipoles, which are independent of the shape of the velocity PDF. The Taylor expansion explains why the Gaussian streaming model works well in predicting the first two redshift-space multipoles, although the velocity PDF is non-Gaussian even on large scales. Indeed, any PDF with the correct first two moments would produce precise results for the monopole down to scales of about $10\,h^{-1}{\rm Mpc}$, and for the quadrupole down to about $30\,h^{-1}{\rm Mpc}$. An accurate model for the hexadecapole needs to include higher-order moments.

preprint2019arXiv

A Search for a Contribution from Axion-Like Particles to the X-Ray Diffuse Background Utilizing the Earth&#39;s Magnetic Field

The Axion Like Particle (ALP) is a hypothetical pseudo-scalar particle beyond the Standard Model, with a compelling possible connection to dark matter and early universe physics. ALPs can be converted into photons via interactions with magnetic fields in the universe, i.e., the so-called inverse Primakoff effect. In this paper, we propose a novel method to explore ALP-induced photons from X-ray data obtained from the {\it Suzaku} satellite, arising from a possible interaction of ALPs with the direction-dependent Earth&#39;s magnetic field viewed from the satellite. {\it Suzaku} data is suitable for this purpose because its low-altitude Earth orbit result in intrinsically low cosmic-ray background radiation. We study whether the X-ray diffuse background (XDB) spectra estimated from the four deep fields collected over eight years, vary with the integrated Earth&#39;s magnetic strength in the direction of each target field at each observation epoch, which amounts to $10^2$ Tm-a value greater than that achieved by terrestrial experiments due to the large coherent length. From the detailed analysis, we did not find evidence of the XDB confidence level spectra having dependence on the Earth&#39;s magnetic strength. We obtained 99 % confidence level upper limit on a possible residual contribution to the cosmic X-ray background (CXB) surface brightness to be $1.6\times 10^{-9}~{\rm ergs~s}^{-1}{\rm cm}^{-2}{\rm sr}^{-1}$ normalized at $10^4$ T${}^2$ m${}^2$ in the 2-6 keV range, which corresponds to 6-15 % of the observed CXB brightness, depending on which model of unresolved point sources are used in the interpretation. It is consistent with 80-90 % of the CXB now being resolved into point sources.

preprint2019arXiv

Mitigating the impact of fiber assignment on clustering measurements from deep galaxy redshift surveys

We examine the impact of fiber assignment on clustering measurements from fiber-fed spectroscopic galaxy surveys. We identify new effects which were absent in previous, relatively shallow galaxy surveys such as Baryon Oscillation Spectroscopic Survey . Specifically, we consider deep surveys covering a wide redshift range from z=0.6 to z=2.4, as in the Subaru Prime Focus Spectrograph survey. Such surveys will have more target galaxies than we can place fibers on. This leads to two effects. First, it eliminates fluctuations with wavelengths longer than the size of the field of view, as the number of observed galaxies per field is nearly fixed to the number of available fibers. We find that we can recover the long-wavelength fluctuation by weighting galaxies in each field by the number of target galaxies. Second, it makes the preferential selection of galaxies in under-dense regions. We mitigate this effect by weighting galaxies using the so-called individual inverse probability. Correcting these two effects, we recover the underlying correlation function at better than 1 percent accuracy on scales greater than 10 Mpc/h.