Statistical Analysis of Random Fields in Cosmology

Asia/Tokyo
Tsukuba Campus, KEK, Tsukuba, Japan

Tsukuba Campus, KEK, Tsukuba, Japan

1-1 Oho, Tsukuba, Ibaraki 305-0801 Japan
Description

For the workshop details, please visit the main page.

 

Venue:

Seminar Hall on the ground floor at Building 4 Seminar Hall, 1F Bldg. 4 (M02 Yon-Go-Kan Bldg.) on the KEK Tsukuba campus.

Come to KEK 

アクセス

Participants
  • Akira Tokwia
  • Asuka Ito
  • Atsushi J. Nishizawa
  • Atsushi Taruya
  • Chiaki Hikage
  • Daisuke Yamauchi
  • Dan Cheng
  • Haruki Takahashi
  • Hayato Shimabukuro
  • Hayato Shimabukuro
  • Hidenaga Watanabe
  • Jessica Cowell
  • Jia Liu
  • Jianing Wang
  • Jingjing Shi
  • Kanmi Nose
  • Kazunori Kohri
  • Kazuyuki Akitsu
  • Keisuke Nagano
  • Ken Osato
  • KRISHANU SENGUPTA
  • Masahiro Takada
  • Masamune Oguri
  • Masanori Iye
  • Matteo Peron
  • Mihoko NOJIRI
  • Mikage Kobayashi
  • Pankaj Saha
  • Riku Yashiki
  • Ritei Shibata
  • Ryo Terasawa
  • Sampurna Bhar
  • Satoshi Kuriki
  • Shintaro Nakano
  • Shiro Ikeda
  • Shogo Masaki
  • Shun ARAI
  • Sora Yamashita
  • Taisei Terawaki
  • Takahiko Matsubara
  • Takahiro Masuda
  • Takahiro Nishimichi
  • Takanori Taniguchi
  • Takashi Hamana
  • Tomo Takahashi
  • Tomoki Katayama
  • Tsutomu Takeuchi
  • Volodymyr Takhistov
  • Yoh-ichi Mototake
  • Yuji CHINONE
  • Yukimura Izawa
  • Yuma Furuta
  • Yuuki Rikimaru
    • 1
      Registration
    • 2
    • Chair: Satoshi Kuriki
      Convener: Satoshi Kuriki
      • 3
        Peak height distributions of Gaussian random fields and their statistical applications

        Motivated by computing p-values for multiple testing of local maxima in signal and change point detections, we study the height distribution of local maxima of smooth isotropic Gaussian random fields parameterized on Euclidean space or spheres. The obtained formulae hold in general in the sense that there are no restrictions on the covariance function of the field except for smoothness and isotropy. The results are based on a characterization of the distribution of the Hessian of the Gaussian field by means of the family of Gaussian orthogonally invariant (GOI) matrices, of which the Gaussian orthogonal ensemble (GOE) is a special case. I will also present related statistical applications such as peak/signal detection in cosmology and change point detection.

        Speaker: Dan Cheng
    • 12:00
      Lunch
    • Chair: Dan Cheng
      Convener: Dan Cheng
      • 4
        Minkowski functionals for isotropic random fields in the Euclidean space and the sphere

        We first introduce the expected Minkowski functional (MF) formulas for the excursion sets of a weakly non-Gaussian smooth isotropic random field. Here, the random field is defined on a bounded index set $T$ in the Euclidean space. The MF formulas contain the boundary correction terms unless the MFs of the index set $T$ vanish. In applications in cosmology, only their leading terms proportional to the volume of the index set $T$ are often used, and the boundary effects are ignored. We numerically evaluate the effects of boundaries in typical settings and demonstrate the importance of boundary corrections.
        Next, we discuss the isotropic (orthogonally invariant) random field on the sphere. We show that modifying the Euclidean case obtains the corresponding MF formulas. The resulting formulas are almost the same as the Euclidean case, and we see that the curvature information cannot be detected from the observed MFs.
        (joint work with T. Matsubara and C. Hikage)

        Speaker: Satoshi Kuriki
      • 5
        Application of Minkowski functional to high energy physics and deep learning
        Speaker: Mihoko Nojiri
    • 14:40
      Tea time
    • Chair: Chiaki Hikage
      • 6
        Minkowski Tensors in Cosmology

        The tensor Minkowski functionals or Minkowski tensors (MTs) are generalizations of the usual Minkowski Functionals, which are scalar quantities. The MTs are a set of statistics defined as integrals over the boundary of an excursion set, with integrands related to symmetric tensor products of position vectors and normals to the boundary. They provide directional or anisotropy information that is not present in the scalar Minkowski functionals. Since 2017, the MTs have been used in cosmology. Efficient and accurate algorithms have been developed for calculation of MTs for analyses of cosmic microwave background radiation anisotropy on the surface of the celestial sphere, and matter and galaxy density fluctuations in 3D space or in 2D slice. Extensive studies have been made to examine the impacts of the nonlinear gravitational evolution, galaxy bias, and redshift space distortion on the MTs of large-scale structures in the universe. The recent developments in the application of the MTs in cosmology are reviewed.

        Speaker: Changbom Park
    • Chair: Yuji Chinone
      • 7
        Precision cosmology with big data
        Speaker: Masahiro Takada
      • 8
        Weak lensing analysis

        I'll describe the basics of weak lensing analysis.

        Speaker: Masamune Oguri
    • Chair: Takahiro Nishimichi
      • 9
        Data analysis methods for the Cosmic Microwave Background (CMB) polarization
        Speaker: Yuji Chinone
      • 10
        Astronomy and Data Science

        The importance of data science is increasing in various fields, and this trend is also evident in the natural sciences. For the past decade, I have been engaged in research applying data science methods to astronomy. One of the highlights of the work is the imaging of a black hole shadow by the Event Horizon Telescope Collaboration. In the presentation, I will provide an overview of my researches in astronomical data analysis, including the results of the EHT collaboration.

        Speaker: Shiro Ikeda
    • Chair: Takahiro Nishimichi
      • 11
        Learning from Topology: Cosmological Parameter Inference from the Large-scale Structure

        A challenge common to different scientific areas is to effectively infer from big, complex, higher-dimensional datasets the underlying theory. Persistent homology is a tool in computational topology developed for recognizing the ``shape” of data. Such topological measures have the advantages that 1) they are stable against experimental noise, 2) they probe multiscale, non-local characteristics of a dataset, 3) they provide interpretable statistics that encode information of all higher-point correlations. In this talk, I will focus on the applications of persistent homology to cosmological inference. I’ll show how the constraints on primordial non-Gaussianities and cosmological parameters derived from persistent homologies are generally tighter than those from the redshift-space power spectrum and bispectrum combined. I’ll also present a recent work in which we proposed a neural network model to map persistence images to cosmological parameters. Through a parameter recovery test, we demonstrate that our model makes accurate and precise estimates, considerably outperforming conventional Bayesian inference approaches. This provides a proof-of-concept that topological data analysis (via persistent homology) and machine learning can be combined for cosmological inference.

        Speaker: Gary Shiu
    • 12
      Photo
    • 12:05
      Lunch
    • Chair: Atsushi J. Nishizawa
      • 13
        Non-Gaussian Statistics for Non-Gaussian Fields

        I will present our recent results analyzing HSC weak lensing data using non-Gaussian statistics, where 30% improvement in cosmological parameter constraints was achieved. I will also summarize our ongoing effort in preparing for future non-Gaussian statistics analysis with Rubin LSST.

        Speaker: Jia Liu
      • 14
        Data Scientific Approach to the Large-Scale Structure in the Universe
        Speaker: Tsutomu Takeuchi
    • 14:40
      Tea time
    • Chair: Hayato Shimabukuro
      • 15
        Exploring the Potential of 3D Scene Reconstruction Techniques in Astronomy

        In recent years, technologies for reconstructing 3D scenes from limited 2D image data, such as Neural Radiance Fields (NeRF) and 3D Gaussian Splatting, have been actively developed. These technologies have gained attention in a wide range of applications, including VR/AR and Simultaneous Localization and Mapping (SLAM). This study explores the potential application of these technologies in the field of astronomy. Specifically, we focus on the possibility of reconstructing 3D dark matter distributions from 2D gravitational lensing data, discussing the feasibility and challenges of this approach. Our research presents new possibilities in the analysis of astronomical data and is expected to contribute to the advancement of related technologies.

        Speaker: Chiaki Hikage
    • Chair: Tsutomu Takeuchi
      • 16
        image based cosmology on the large-scale structure

        After the great success of the perturbation theory of the large-scale structure, we definitely need more accurate and robust tools to expand the analysis to utilize more and more information embedded in the large-scale structure.
        One possible approach which has already got great success is to extent the order of expansion to higher order or to use the interpolation method calibrated based on the set of numerical simulations.
        Another approach recently proposed is not based on the summary statistics but on the field level inference, that is based on the "image".
        In this talk, I will revisit the recent works based on the image of the large-scale structure, and discuss our recent cosmological analysis based on the machine learning methods.

        Speaker: Atsushi J. Nichizawa
      • 17
        21cm cosmology with machine learning
        Speaker: Hayato Shimabukuro
    • 18
      Banquet Andryu

      Andryu

      512-3 Sa, Tsukuba, Ibaraki 300-3251

      https://andryu.owst.jp/

    • Chair: Kazuyuki Akitsu
      • 19
        Consistency relations of the kurtosis in large scale structure and modified gravity

        Skewness and kurtosis parameters are the quantities to probe non-Gaussianities of density fluctuations.
        By introducing some spatial derivatives, we can define multiple skewness and kurtosis parameters.
        We argue that the consistency relations among observables of the skewness and kurtosis parameters can be derived, particularly focusing on the kurtosis parameters in this talk.
        The consistency relations can serve as a test of gravity theories.

        Speaker: Sora Yamashita
      • 20
        Emulators for cosmological inference

        Cosmological large-scale structure is a nonlinear stochastic process governed mainly by gravity. Its statistical properties depend on the initial conditions and the energy components that constitute the universe. Therefore, extracting information on these fundamental cosmological problems from observational data is an inverse problem. The computational cost of numerical simulations is a bottleneck in this problem setting. We discuss the use of emulators as a working example to replace simulations and alleviate this issue.

        Speaker: Takahiro Nishimichi
    • Chair: Satoshi Kuriki
      • 21
        Interpretable reduced modeling of large-scale pattern dynamics - from materials science to astrophysics -
        Speaker: Yoh-ichi Mototake
    • 12:00
      Lunch
    • Chair: Atsushi Taruya
      • 22
        How much information is left in galaxy clustering beyond the power spectrum?
        Speaker: Kazuyuki Akitsu
      • 23
        Field-level inference with differentiable GridSPT

        In galaxy clustering measurements, summary statistics, e.g., 2-pt correlation functions and power spectra, have been widely used to extract cosmological information. These statistics summarise the information of the observed galaxy number density field and theoretical models such as perturbation theory can accurately predict the statistics up to mildly non-linear regime. On the other hand, the statistics do not contain the full information of the observed field and, especially for non-linear fields, a large fraction of the information is not accessible solely from the statistics. In this work, we focus on the field-level approach that directly compares the observed field itself with the non-linear density field calculated by the forward modelling. To this end, we develop differentiable GridSPT, which is a forward modelling method based on standard perturbation theory and utilizes accelerated Fourier transforms and automatic differentiation with GPUs. This functionaltiy makes statistical inference more efficient and convergent. We discuss the validity of the field-level inference through reconstruction of the initial density field and cosmological parameters from non-linear density fields from N-body simulations.

        Speaker: Ken Osato
    • 14:40
      Tea time
    • Chair: Ken Osato
      • 24
        Galaxy shape/alignment statistics as a novel cosmological probe

        Mapping the large-scale structure of the universe with galaxy surveys is one of the main science drivers for cosmology. So far, most of the analysis has been made wth the galaxy positional information, ignoring the individual shapes and orientations. In this talk, we consider the galaxy intrinsic orientation (intrinsic alignment, IA) as a novel probe, and show that the spatial correlation of the IAs is beneficial to tighten cosmological constraints. Precision modeling of the IA statistics is thus very important, and we will present an improved model of the IA power spectrum, which properly takes into account the nonlinear effects of gravity and redshift-space distortions based on perturbation theory.

        Speaker: Atsushi Taruya
      • 25
        Symmetry breaking in the large scale distribution of spin vector of galaxies?

        We are compiling spin catalogs of galaxies from PanStarrs, SDSS, DES and HSC surveys to make statistical studies to search for any symmetry breaking in the large scale distribution of spin vectors of galaxies. The method and some tentative results of analyses will be reported.

        Speaker: Masanori Iye