Showing 1–50 of 58 results for author: Yoshikawa, K

High-energy extension of the gamma-ray band observable with an electron-tracking Compton camera

Authors: Tomohiko Oka, Shingo Ogio, Mitsuru Abe, Kenji Hamaguchi, Tomonori Ikeda, Hidetoshi Kubo, Shunsuke Kurosawa, Kentaro Miuchi, Yoshitaka Mizumura, Yuta Nakamura, Tatsuya Sawano, Atsushi Takada, Taito Takemura, Toru Tanimori, Kei Yoshikawa

Abstract: Although the MeV gamma-ray band is a promising energy-band window in astrophysics, the current situation of MeV gamma-ray astronomy significantly lags behind those of the other energy bands in angular resolution and sensitivity. An electron-tracking Compton camera (ETCC), a next-generation MeV detector, is expected to revolutionize the situation. An ETCC tracks each Compton-recoil electron with a… ▽ More Although the MeV gamma-ray band is a promising energy-band window in astrophysics, the current situation of MeV gamma-ray astronomy significantly lags behind those of the other energy bands in angular resolution and sensitivity. An electron-tracking Compton camera (ETCC), a next-generation MeV detector, is expected to revolutionize the situation. An ETCC tracks each Compton-recoil electron with a gaseous electron tracker and determines the incoming direction of each gamma-ray photon; thus, it has a strong background rejection power and yields a better angular resolution than classical Compton cameras. Here, we study ETCC events in which the Compton-recoil electrons do not deposit all energies to the electron tracker but escape and hit the surrounding pixel scintillator array (PSA). We developed an analysis method for this untapped class of events and applied it to laboratory and simulation data. We found that the energy spectrum obtained from the simulation agreed with that of the actual data within a factor of 1.2. We then evaluated the detector performance using the simulation data. The angular resolution for the new-class events was found to be twice as good as in the previous study at the energy range 1.0--2.0~MeV, where both analyses overlap. We also found that the total effective area is dominated by the contribution of the double-hit events above an energy of 1.5~MeV. Notably, applying this new method extends the sensitive energy range with the ETCC from 0.2--2.1 MeV in the previous studies to up to 3.5~MeV. Adjusting the PSA dynamic range should improve the sensitivity in even higher energy gamma-rays. The development of this new analysis method would pave the way for future observations by ETCC to fill the MeV-band sensitivity gap in astronomy. △ Less

Submitted 8 March, 2024; originally announced March 2024.

Comments: 11 pages, 8 figures, Accepted for publication in NIM A

Background contributions in the electron-tracking Compton camera onboard SMILE-2+

Authors: Tomonori Ikeda, Atsushi Takada, Taito Takemura, Kei Yoshikawa, Yuta nakamura, Ken Onozaka, Mitsuru Abe, Toru Tanimori

Abstract: The Mega electron volt (MeV) gamma-ray observation is a promising diagnostic tool for observing the universe. However, the sensitivity of MeV gamma-ray telescopes is limited due to peculiar backgrounds, restricting the application of MeV gamma rays for observation. Identification of backgrounds is crucial for designing next-generation telescopes. Therefore, herein, we assessed the background contr… ▽ More The Mega electron volt (MeV) gamma-ray observation is a promising diagnostic tool for observing the universe. However, the sensitivity of MeV gamma-ray telescopes is limited due to peculiar backgrounds, restricting the application of MeV gamma rays for observation. Identification of backgrounds is crucial for designing next-generation telescopes. Therefore, herein, we assessed the background contribution in the electron-tracking Compton camera (ETCC) on board the SMILE- 2+ balloon experiment. This assessment was performed using the Monte Carlo simulation. The results revealed that the background below 400 keV existed due to the atmospheric gamma-ray background, the cosmic-ray/secondary-particle background, and the accidental background. On the other hand, the unresolved background component, which was not likely to be relevant to direct Compton-scattering events in the ETCC, was confirmed above 400 keV. Overall, this study demonstrated that the Compton-kinematics test provides a powerful tool to remove the background and principally improves the signal-to-noise ratio at 400 keV by an order of magnitude. △ Less

Submitted 28 September, 2023; v1 submitted 5 June, 2023; originally announced June 2023.

Comments: 11 pages, 18 figures

The spatial distribution of impact craters on Ryugu

Authors: Naoyuki Hirata, Tomokatsu Morota, Yuichiro Cho, Masanori Kanamaru, Sei-ichiro Watanabe, Seiji Sugita, Naru Hirata, Yukio Yamamoto, Rina Noguchi, Yuri Shimaki, Eri Tatsumi, Kazuo Yoshioka, Hirotaka Sawada, Yasuhiro Yokota, Naoya Sakatani, Masahiko Hayakawa, Moe Matsuoka, Rie Honda, Shingo Kameda, Mamabu Yamada, Toru Kouyama, Hidehiko Suzuki, Chikatoshi Honda, Kazunori Ogawa, Yuichi Tsuda , et al. (18 additional authors not shown)

Abstract: Asteroid 162173 Ryugu has numerous craters. The initial measurement of impact craters on Ryugu, by Sugita et al. (2019), is based on Hayabusa2 ONC images obtained during the first month after the arrival of Hayabusa2 in June 2018. Utilizing new images taken until February 2019, we constructed a global impact crater catalogue of Ryugu, which includes all craters larger than 20 m in diameter on the… ▽ More Asteroid 162173 Ryugu has numerous craters. The initial measurement of impact craters on Ryugu, by Sugita et al. (2019), is based on Hayabusa2 ONC images obtained during the first month after the arrival of Hayabusa2 in June 2018. Utilizing new images taken until February 2019, we constructed a global impact crater catalogue of Ryugu, which includes all craters larger than 20 m in diameter on the surface of Ryugu. As a result, we identified 77 craters on the surface of Ryugu. Ryugu shows variation in crater density which cannot be explained by the randomness of cratering; there are more craters at lower latitudes and fewer at higher latitudes, and fewer craters in the western bulge (160 E - 290 E) than in the region around the meridian (300 E - 30 E). This variation implies a complicated geologic history for Ryugu. It seems that the longitudinal variation in crater density simply indicates variation in the crater ages; the cratered terrain around the meridian seems to be geologically old while the western bulge is relatively young. The latitudinal variation in crater density suggests that the equatorial ridge of Ryugu is a geologically old structure; however, this could be alternatively explained by a collision with many fission fragments during a short rotational period of Ryugu in the past. △ Less

Submitted 11 May, 2022; originally announced May 2022.

Comments: 18 pages 3 figures

Journal ref: Icarus Volume 338, 1 March 2020, 113527

A 400 Trillion-Grid Vlasov Simulation on Fugaku Supercomputer: Large-Scale Distribution of Cosmic Relic Neutrinos in a Six-dimensional Phase Space

Authors: Kohji Yoshikawa, Satoshi Tanaka, Naoki Yoshida

Abstract: We report a Vlasov simulation of cosmic relic neutrinos combined with N-body simulation of cold dark matter in the context of large-scale structure formation in the Universe performed on Fugaku supercomputer. Gravitational dynamics of the neutrinos is followed, for the first time, by directly integrating the Vlasov equation in a six-dimensional phase space. Our largest simulation combines the Vlas… ▽ More We report a Vlasov simulation of cosmic relic neutrinos combined with N-body simulation of cold dark matter in the context of large-scale structure formation in the Universe performed on Fugaku supercomputer. Gravitational dynamics of the neutrinos is followed, for the first time, by directly integrating the Vlasov equation in a six-dimensional phase space. Our largest simulation combines the Vlasov simulation on 400 trillion grids and 330 billion-body calculations in a self-consistent manner, and reproduces accurately the nonlinear dynamics of neutrinos in the Universe. The novel high-order Vlasov solver is optimized by combining an array of state-of-the-art numerical schemes and fully utilizing the SIMD instructions on the A64FX processors. Time-To-Solution of our simulation is an order of magnitude shorter than the largest N-body simulations. The performance scales excellently with up to 147,456 nodes (7 million CPU cores) on Fugaku; the weak and strong scaling efficiencies are 82% - 96% and 82% - 93%, respectively. △ Less

Submitted 29 October, 2021; originally announced October 2021.

Comments: 11 pages, 8 figures. 2021 ACM Gordon-Bell Prize Finalist. A simulation movie is available at https://youtu.be/p2fHt8f7chQ

Journal ref: SC '21: Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis, November 2021, Article No.5, pages 1-11

First observation of MeV gamma-ray universe with bijective imaging spectroscopy using the Electron-Tracking Compton Telescope aboard SMILE-2+

Authors: Atsushi Takada, Taito Takemura, Kei Yoshikawa, Yoshitaka Mizumura, Tomonori Ikeda, Yuta Nakamura, Ken Onozaka, Mitsuru Abe, Kenji Hamaguchi, Hidetoshi Kubo, Shunsuke Kurosawa, Kentaro Miuchi, Kaname Saito, Tatsuya Sawano, Toru Tanimori

Abstract: MeV gamma-rays provide a unique window for the direct measurement of line emissions from radioisotopes, but observations have made little significant progress after COMPTEL/{\it CGRO}. To observe celestial objects in this band, we are developing an electron-tracking Compton camera (ETCC), which realizes both bijective imaging spectroscopy and efficient background reduction gleaned from the recoil… ▽ More MeV gamma-rays provide a unique window for the direct measurement of line emissions from radioisotopes, but observations have made little significant progress after COMPTEL/{\it CGRO}. To observe celestial objects in this band, we are developing an electron-tracking Compton camera (ETCC), which realizes both bijective imaging spectroscopy and efficient background reduction gleaned from the recoil electron track information. The energy spectrum of the observation target can then be obtained by a simple ON-OFF method using a correctly defined point spread function on the celestial sphere. The performance of celestial object observations was validated on the second balloon SMILE-2+ installed with an ETCC having a gaseous electron tracker with a volume of 30$\times$30$\times$30 cm$^3$. Gamma-rays from the Crab nebula were detected with a significance of 4.0$σ$ in the energy range 0.15--2.1 MeV with a live time of 5.1 h, as expected before launching. Additionally, the light curve clarified an enhancement of gamma-ray events generated in the Galactic center region, indicating that a significant proportion of the final remaining events are cosmic gamma rays. Independently, the observed intensity and time variation were consistent with the pre-launch estimates except in the Galactic center region. The estimates were based on the total background of extragalactic diffuse, atmospheric, and instrumental gamma-rays after accounting for the variations in the atmospheric depth and rigidity during the level flight. The Crab results and light curve strongly support our understanding of both the detection sensitivity and the background in real observations. This work promises significant advances in MeV gamma-ray astronomy. △ Less

Submitted 24 March, 2022; v1 submitted 30 June, 2021; originally announced July 2021.

Journal ref: The Astrophyssical Journal, Volume 930 (2022), 6

Development of Convolutional Neural Networks for an Electron-Tracking Compton Camera

Authors: Tomonori Ikeda, Atsushi Takada, Mitsuru Abe, Kei Yoshikawa, Masaya Tsuda, Shingo Ogio, Shinya Sonoda, Yoshitaka Mizumura, Yura Yoshida, Toru Tanimori

Abstract: Electron-tracking Compton camera, which is a complete Compton camera with tracking Compton scattering electron by a gas micro time projection chamber, is expected to open up MeV gamma-ray astronomy. The technical challenge for achieving several degrees of the point spread function is the precise determination of the electron-recoil direction and the scattering position from track images. We attemp… ▽ More Electron-tracking Compton camera, which is a complete Compton camera with tracking Compton scattering electron by a gas micro time projection chamber, is expected to open up MeV gamma-ray astronomy. The technical challenge for achieving several degrees of the point spread function is the precise determination of the electron-recoil direction and the scattering position from track images. We attempted to reconstruct these parameters using convolutional neural networks. Two network models were designed to predict the recoil direction and the scattering position. These models marked 41$~$degrees of the angular resolution and 2.1$~$mm of the position resolution for 75$~$keV electron simulation data in Argon-based gas at 2$~$atm pressure. In addition, the point spread function of ETCC was improved to 15$~$degrees from 22$~$degrees for experimental data of 662$~$keV gamma-ray source. These performances greatly surpassed that using the traditional analysis. △ Less

Submitted 1 July, 2021; v1 submitted 6 May, 2021; originally announced May 2021.

Journal ref: PTEP (2021), 083F01

Cosmological Vlasov-Poisson Simulations of Structure Formation with Relic Neutrinos: Nonlinear Clustering and the Neutrino Mass

Authors: Kohji Yoshikawa, Satoshi Tanaka, Naoki Yoshida, Shun Saito

Abstract: We present the results of cosmological simulations of large-scale structure formation with massive neutrinos. The phase-space distribution of the cosmic relic neutrinos is followed, for the first time, by directly integrating the six-dimensional Vlasov-Poisson equations. Our novel approach allows us to represent free streaming and clustering of neutrinos, and their gravitational interaction with c… ▽ More We present the results of cosmological simulations of large-scale structure formation with massive neutrinos. The phase-space distribution of the cosmic relic neutrinos is followed, for the first time, by directly integrating the six-dimensional Vlasov-Poisson equations. Our novel approach allows us to represent free streaming and clustering of neutrinos, and their gravitational interaction with cold dark matter accurately. We thus obtain solutions for the collisionless dynamics independent of conventional N-body methods. We perform a suite of hybrid N-body/Vlasov simulations with varying the neutrino mass, and systematically examine the dynamical effects of massive neutrinos on large-scale structure formation. Our simulations show characteristic large-scale clustering of the neutrinos and their coherent streaming motions relative to dark matter. The effective local neutrino "temperature" around massive galaxy clusters varies by several percent with respect to the cosmic mean; the neutrinos in clusters can be hotter or colder depending on the neutrino mass. We study a number of statistics of the large-scale structure and of dark matter halos in comparison with those obtained by N-body simulations and/or by perturbation theory. Our simulations mark an important milestone in numerical cosmology, and pave a new way to study cosmic structure formation with massive neutrinos. △ Less

Submitted 1 October, 2020; v1 submitted 1 October, 2020; originally announced October 2020.

Comments: 21 pages, accepted for publication in ApJ

Deeply cooled core of the Phoenix galaxy cluster imaged by ALMA with the Sunyaev-Zel'dovich effect

Authors: T. Kitayama, S. Ueda, T. Akahori, E. Komatsu, R. Kawabe, K. Kohno, S. Takakuwa, T. Tsutsumi, M. Takizawa, K. Yoshikawa

Abstract: We present measurements of the thermal Sunyaev-Zel'dovich effect (SZE) toward SPT-CL J2334-4243 (the Phoenix galaxy cluster) at z=0.597 by the Atacama Large Millimeter/submillimeter Array (ALMA) in Band 3. The SZE is imaged at 5" resolution (corresponding to the physical scale of 23kpc/h) within 200kpc/h from the central galaxy with the peak signal-to-noise ratio exceeding 11. Combined with the Ch… ▽ More We present measurements of the thermal Sunyaev-Zel'dovich effect (SZE) toward SPT-CL J2334-4243 (the Phoenix galaxy cluster) at z=0.597 by the Atacama Large Millimeter/submillimeter Array (ALMA) in Band 3. The SZE is imaged at 5" resolution (corresponding to the physical scale of 23kpc/h) within 200kpc/h from the central galaxy with the peak signal-to-noise ratio exceeding 11. Combined with the Chandra X-ray image, the ALMA SZE data further allow for non-parametric deprojection of electron temperature, density, and entropy. Our method can minimize contamination by the central AGN and the X-ray absorbing gas within the cluster, both of which largely affect the X-ray spectrum. We find no significant asymmetry or disturbance in the SZE image within the current measurement errors. The detected SZE signal agrees well with the average pressure profile of local cool-core clusters. Unlike any other known clusters, however, gas temperature drops by at least a factor of 5 toward the center. We identify ~6x10^{11} M_sun cool gas with temperature ~3keV in the inner 20kpc/h. Taken together, our results imply that the gas is indeed cooling efficiently and nearly isobarically down to this radius in the Phoenix cluster. △ Less

Submitted 30 January, 2020; originally announced January 2020.

Comments: 18 pages, 16 figures. Accepted for publication in PASJ

Cool core disturbed: Observational evidence for coexistence of sub-sonic sloshing gas and stripped shock-heated gas around the core of RX J1347.5-1145

Authors: Shutaro Ueda, Tetsu Kitayama, Masamune Oguri, Eiichiro Komatsu, Takuya Akahori, Daisuke Iono, Takumi Izumi, Ryohei Kawabe, Kotaro Kohno, Hiroshi Matsuo, Naomi Ota, Yasushi Suto, Shigehisa Takakuwa, Motokazu Takizawa, Takahiro Tsutsumi, Kohji Yoshikawa

Abstract: RXJ1347.5-1145 (z = 0.451) is one of the most luminous X-ray galaxy clusters, which hosts a prominent cool core and exhibits a signature of a major merger. We present the first direct observational evidence for sub-sonic nature of sloshing motion of the cool core. We find that a residual X-ray image from the Chandra X-ray Observatory after removing the global emission shows a clear dipolar pattern… ▽ More RXJ1347.5-1145 (z = 0.451) is one of the most luminous X-ray galaxy clusters, which hosts a prominent cool core and exhibits a signature of a major merger. We present the first direct observational evidence for sub-sonic nature of sloshing motion of the cool core. We find that a residual X-ray image from the Chandra X-ray Observatory after removing the global emission shows a clear dipolar pattern characteristic of gas sloshing, whereas we find no significant residual in the Sunyaev-Zel'dovich effect (SZE) image from the Atacama Large Millimeter/submillimeter Array (ALMA). We estimate the equation of state of perturbations in the gas from the X-ray and SZE residual images. The inferred velocity is 420 +310 -420 km s-1, which is much lower than the adiabatic sound speed of the intracluster medium in the core. We thus conclude that the perturbation is nearly isobaric, and gas sloshing motion is consistent with being in pressure equilibrium. Next, we report evidence for gas stripping of an infalling subcluster, which likely shock-heats gas to high temperature well in excess of 20 keV. Using mass distribution inferred from strong lensing images of the Hubble Space Telescope (HST), we find that the mass peak is located away from the peak position of stripped gas with statistical significance of > 5σ. Unlike for the gas sloshing, the velocity inferred from the equation of state of the excess hot gas is comparable to the adiabatic sound speed expected for the 20 keV intracluster medium. All of the results support that the southeast substructure is created by a merger. On the other hand, the positional offset between the mass and the gas limits the self-interaction cross section of dark matter to be less than 3.7 h-1 cm2 g-1 (95% CL). △ Less

Submitted 28 August, 2018; originally announced August 2018.

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

Properties of the cosmological filament between two clusters: A possible detection of a large-scale accretion shock by $Suzaku$

Authors: H. Akamatsu, Y. Fujita, T. Akahori, Y. Ishisaki, K. Hayashida, A. Hoshino, F. Mernier, K. Yoshikawa, K. Sato, J. S. Kaastra

Abstract: We report on the results of a $Suzaku$ observation of the plasma in the filament located between the two massive clusters of galaxies Abell 399 and Abell 401. Abell 399 ($z$=0.0724) and Abell 401 ($z$=0.0737) are expected to be in the initial phase of a cluster merger. In the region between the two clusters, we find a clear enhancement in the temperature of the filament plasma from 4 keV (expected… ▽ More We report on the results of a $Suzaku$ observation of the plasma in the filament located between the two massive clusters of galaxies Abell 399 and Abell 401. Abell 399 ($z$=0.0724) and Abell 401 ($z$=0.0737) are expected to be in the initial phase of a cluster merger. In the region between the two clusters, we find a clear enhancement in the temperature of the filament plasma from 4 keV (expected value from a typical cluster temperature profile) to $kT\sim$6.5 keV. Our analysis also shows that filament plasma is present out to a radial distance of 15' (1.3 Mpc) from a line connecting the two clusters. The temperature profile is characterized by an almost flat radial shape with $kT\sim$6-7 keV within 10' or $\sim$0.8 Mpc. Across $r$=8'~from the axis, the temperature of the filament plasma shows a drop from 6.3 keV to 5.1 keV, indicating the presence of a shock front. The Mach number based on the temperature drop is estimated to be ${\cal M}\sim$1.3. We also successfully determined the abundance profile up to 15' (1.3 Mpc), showing an almost constant value ($Z$=0.3 solar) at the cluster outskirt. We estimated the Compton $y$-parameter to be $\sim$14.5$\pm1.3\times10^{-6}$, which is in agreement with $Planck$'s results (14-17$\times10^{-6}$ on the filament). The line of sight depth of the filament is $l\sim$1.1 Mpc, indicating that the geometry of filament is likely a pancake shape rather than cylindrical. The total mass of the filamentary structure is $\sim$7.7$\times10^{13}~\rm M_{\odot}$. We discuss a possible interpretation of the drop of X-ray emission at the rim of the filament, which was pushed out by the merging activity and formed by the accretion flow induced by the gravitational force of the filament. △ Less

Submitted 19 April, 2017; originally announced April 2017.

Comments: 8 pages, 8 figures, accepted for publication in A&A

Imaging Polarimeter for a Sub-MeV Gamma-Ray All-Sky Survey using an Electron-Tracking Compton Camera

Authors: Shotaro Komura, Atsushi Takada, Yoshitaka Mizumura, Shohei Miyamoto, Taito Takemura, Tetsuro Kishimoto, Hidetoshi Kubo, Shunsuke Kurosawa, Yoshihiro Matsuoka, Kentaro Miuchi, Tetsuya Mizumoto, Yuma Nakamasu, Kiseki Nakamura, Makoto Oda, Joseph D. Parker, Tatsuya Sawano, Shinya Sonoda, Toru Tanimori, Dai Tomono, Kei Yoshikawa

Abstract: X-ray and gamma-ray polarimetry is a promising tool to study the geometry and the magnetic configuration of various celestial objects, such as binary black holes or gamma-ray bursts (GRBs). However, statistically significant polarizations have been detected in few of the brightest objects. Even though future polarimeters using X-ray telescopes are expected to observe weak persistent sources, there… ▽ More X-ray and gamma-ray polarimetry is a promising tool to study the geometry and the magnetic configuration of various celestial objects, such as binary black holes or gamma-ray bursts (GRBs). However, statistically significant polarizations have been detected in few of the brightest objects. Even though future polarimeters using X-ray telescopes are expected to observe weak persistent sources, there are no effective approaches to survey transient and serendipitous sources with a wide field of view (FoV). Here we present an electron-tracking Compton camera (ETCC) as a highly-sensitive gamma-ray imaging polarimeter. The ETCC provides powerful background rejection and a high modulation factor over a FoV of up to 2$π$ sr thanks to its excellent imaging based on a well-defined point spread function. Importantly, we demonstrated for the first time the stability of the modulation factor under realistic conditions of off-axis incidence and huge backgrounds using the SPring-8 polarized X-ray beam. The measured modulation factor of the ETCC was 0.65 $\pm$ 0.01 at 150 keV for an off-axis incidence with an oblique angle of 30$^\circ$ and was not degraded compared to the 0.58 $\pm$ 0.02 at 130 keV for on-axis incidence. These measured results are consistent with the simulation results. Consequently, we found that the satellite-ETCC proposed in Tanimori et al. (2015) would provide all-sky surveys of weak persistent sources of 13 mCrab with 10% polarization for a 10$^{7}$ s exposure and over 20 GRBs down to a $6\times10^{-6}$ erg cm$^{-2}$ fluence and 10% polarization during a one-year observation. △ Less

Submitted 22 March, 2017; originally announced March 2017.

Comments: 22 pages, 19 figures, 1 table, accepted for publication in ApJ

Multidimensional VlasovPoisson Simulations with High-order Monotonicity- and Positivity-preserving Schemes

Authors: Satoshi Tanaka, Kohji Yoshikawa, Takashi Minoshima, Naoki Yoshida

Abstract: We develop new numerical schemes for Vlasov--Poisson equations with high-order accuracy. Our methods are based on a spatially monotonicity-preserving (MP) scheme and are modified suitably so that positivity of the distribution function is also preserved. We adopt an efficient semi-Lagrangian time integration scheme that is more accurate and computationally less expensive than the three-stage TVD R… ▽ More We develop new numerical schemes for Vlasov--Poisson equations with high-order accuracy. Our methods are based on a spatially monotonicity-preserving (MP) scheme and are modified suitably so that positivity of the distribution function is also preserved. We adopt an efficient semi-Lagrangian time integration scheme that is more accurate and computationally less expensive than the three-stage TVD Runge-Kutta integration. We apply our spatially fifth- and seventh-order schemes to a suite of simulations of collisionless self-gravitating systems and electrostatic plasma simulations, including linear and nonlinear Landau damping in one dimension and Vlasov--Poisson simulations in a six-dimensional phase space. The high-order schemes achieve a significantly improved accuracy in comparison with the third-order positive-flux-conserved scheme adopted in our previous study. With the semi-Lagrangian time integration, the computational cost of our high-order schemes does not significantly increase, but remains roughly the same as that of the third-order scheme. Vlasov--Poisson simulations on $128^3 \times 128^3$ mesh grids have been successfully performed on a massively parallel computer. △ Less

Submitted 31 October, 2017; v1 submitted 23 February, 2017; originally announced February 2017.

Comments: 27 pages, 21 figures. Accepted for publication in Astrophysical Journal

Establishment of Imaging Spectroscopy of Nuclear Gamma-Rays based on Geometrical Optics

Authors: Toru Tanimori, Yoshitaka Mizumura, Atsushi Takada, Shohei Miyamoto, Taito Takemura, Tetsuro Kishimoto, Shotaro Komura, Hidetoshi Kubo, Shunsuke Kurosawa, Yoshihiro Matsuoka, Kentaro Miuchi, Tetsuya Mizumoto, Yuma Nakamasu, Kiseki Nakamura, Joseph D. Parker, Tatsuya Sawano, Shinya Sonoda, Dai Tomono, Kei Yoshikawa

Abstract: Since the discovery of nuclear gamma-rays, its imaging has been limited to pseudo imaging, such as Compton Camera (CC) and coded mask. Pseudo imaging does not keep physical information (intensity, or brightness in Optics) along a ray, and thus is capable of no more than qualitative imaging of bright objects. To attain quantitative imaging, cameras that realize geometrical optics is essential, whic… ▽ More Since the discovery of nuclear gamma-rays, its imaging has been limited to pseudo imaging, such as Compton Camera (CC) and coded mask. Pseudo imaging does not keep physical information (intensity, or brightness in Optics) along a ray, and thus is capable of no more than qualitative imaging of bright objects. To attain quantitative imaging, cameras that realize geometrical optics is essential, which would be, for nuclear MeV gammas, only possible via complete reconstruction of the Compton process. Recently we have revealed that "Electron Tracking Compton Camera" (ETCC) provides a well-defined Point Spread Function (PSF). The information of an incoming gamma is kept along a ray with the PSF and that is equivalent to geometrical optics. Here we present an imaging-spectroscopic measurement with the ETCC. Our results highlight the intrinsic difficulty with CCs in performing accurate imaging, and show that the ETCC surmounts this problem. The imaging capability also helps the ETCC suppress the noise level dramatically by ~3 orders of magnitude without a shielding structure. Furthermore, full reconstruction of Compton process with the ETCC provides spectra free of Compton edges. These results mark the first proper imaging of nuclear gammas based on the genuine geometrical optics. △ Less

Submitted 5 February, 2017; originally announced February 2017.

Comments: 22 pages, 8 figures

Journal ref: Scientific Reports 7, 41511 (2017)

The Sunyaev-Zel'dovich Effect at Five Arc-seconds: RXJ1347.5-1145 Imaged by ALMA

Authors: T. Kitayama, S. Ueda, S. Takakuwa, T. Tsutsumi, E. Komatsu, T. Akahori, D. Iono, T. Izumi, R. Kawabe, K. Kohno, H. Matsuo, N. Ota, Y. Suto, M. Takizawa, K. Yoshikawa

Abstract: We present the first image of the thermal Sunyaev-Zel'dovich effect (SZE) obtained by the Atacama Large Millimeter/submillimeter Array (ALMA). Combining 7-m and 12-m arrays in Band 3, we create an SZE map toward a galaxy cluster RXJ1347.5-1145 with 5 arc-second resolution (corresponding to the physical size of 20 kpc/h), the highest angular and physical spatial resolutions achieved to date for ima… ▽ More We present the first image of the thermal Sunyaev-Zel'dovich effect (SZE) obtained by the Atacama Large Millimeter/submillimeter Array (ALMA). Combining 7-m and 12-m arrays in Band 3, we create an SZE map toward a galaxy cluster RXJ1347.5-1145 with 5 arc-second resolution (corresponding to the physical size of 20 kpc/h), the highest angular and physical spatial resolutions achieved to date for imaging the SZE, while retaining extended signals out to 40 arc-seconds. The 1-sigma statistical sensitivity of the image is 0.017 mJy/beam or 0.12 mK_CMB at the 5 arc-second full width at half maximum. The SZE image shows a good agreement with an electron pressure map reconstructed independently from the X-ray data and offers a new probe of the small-scale structure of the intracluster medium. Our results demonstrate that ALMA is a powerful instrument for imaging the SZE in compact galaxy clusters with unprecedented angular resolution and sensitivity. As the first report on the detection of the SZE by ALMA, we present detailed analysis procedures including corrections for the missing flux, to provide guiding methods for analyzing and interpreting future SZE images by ALMA. △ Less

Submitted 29 July, 2016; originally announced July 2016.

Comments: 20 pages, 13 figures. Accepted for publication in PASJ

Japanese Cosmic Dawn/Epoch of Reionization Science with the Square Kilometre Array

Authors: Kenji Hasegawa, Shinsuke Asaba, Kiyotomo Ichiki, Akio K. Inoue, Susumu Inoue, Tomoaki Ishiyama, Hayato Shimabukuro, Keitaro Takahashi, Hiroyuki Tashiro, Hidenobu Yajima, Shu-ichiro Yokoyama, Kohji Yoshikawa, Shintaro Yoshiura

Abstract: Cosmic reionization is known to be a major phase transition of the gas in the Universe. Since astronomical objects formed in the early Universe, such as the first stars, galaxies and black holes, are expected to have caused cosmic reionization, the formation history and properties of such objects are closely related to the reionization process. In spite of the importance of exploring reionization,… ▽ More Cosmic reionization is known to be a major phase transition of the gas in the Universe. Since astronomical objects formed in the early Universe, such as the first stars, galaxies and black holes, are expected to have caused cosmic reionization, the formation history and properties of such objects are closely related to the reionization process. In spite of the importance of exploring reionization, our understandings regarding reionization is not sufficient yet. Square Kilometre Array (SKA) is a next-generation large telescope that will be operated in the next decade. Although several programs of next-generation telescopes are currently scheduled, the SKA will be the unique telescope with a potential to directly observe neutral hydrogen up to z~30, and provide us with valuable information on the Cosmic Dawn (CD) and the Epoch of Reionization (EoR). The early science with the SKA will start in a few years; it is thus the time for us to elaborate a strategy for CD/EoR Science with the SKA. The purpose of this document is to introduce Japanese scientific interests in the SKA project and to report results of our investigation. △ Less

Submitted 7 March, 2016; originally announced March 2016.

Comments: 14 pages, 12 figures, Japan SKA consortium (SKA-JP) Cosmic Dawn/Epoch of Reionization Science book

Cosmology with the Square Kilometre Array by SKA-Japan

Authors: Daisuke Yamauchi, Kiyotomo Ichiki, Kazunori Kohri, Toshiya Namikawa, Yoshihiko Oyama, Toyokazu Sekiguchi, Hayato Shimabukuro, Keitaro Takahashi, Tomo Takahashi, Shuichiro Yokoyama, Kohji Yoshikawa

Abstract: In the past several decades, the standard cosmological model has been established and its parameters have been measured to a high precision, while there are still many of the fundamental questions in cosmology; such as the physics in the very early Universe, the origin of the cosmic acceleration and the nature of the dark matter. The future world's largest radio telescope, Square Kilometre Array (… ▽ More In the past several decades, the standard cosmological model has been established and its parameters have been measured to a high precision, while there are still many of the fundamental questions in cosmology; such as the physics in the very early Universe, the origin of the cosmic acceleration and the nature of the dark matter. The future world's largest radio telescope, Square Kilometre Array (SKA), will be able to open the new frontier of cosmology and will be one of the most powerful tools for cosmology in the next decade. The cosmological surveys conducted by the SKA would have the potential not only to answer these fundamental questions but also deliver the precision cosmology. In this article we briefly review the role of the SKA from the view point of the modern cosmology. The cosmology science led by the SKA-Japan Consortium (SKA-JP) Cosmology Science Working Group is also discussed. △ Less

Submitted 12 December, 2016; v1 submitted 7 March, 2016; originally announced March 2016.

Comments: 18 pages, 11 figures; v2: reference added; v3: accepted for publication in PASJ

Report number: RESCEU-2/16, KEK-TH-1893, KEK-Cosmo-192

Journal ref: Publ. Astron. Soc. Jpn. 68 (2016) 6

A new ray-tracing scheme for 3D diffuse radiation transfer on highly parallel architectures

Authors: Satoshi Tanaka, Kohji Yoshikawa, Takashi Okamoto, Kenji Hasegawa

Abstract: We present a new numerical scheme to solve the transfer of diffuse radiation on three-dimensional mesh grids which is efficient on processors with highly parallel architecture such as recently popular GPUs and CPUs with multi- and many-core architectures. The scheme is based on the ray-tracing method and the computational cost is proportional to $N_{\rm m}^{5/3}$ where $N_{\rm m}$ is the number of… ▽ More We present a new numerical scheme to solve the transfer of diffuse radiation on three-dimensional mesh grids which is efficient on processors with highly parallel architecture such as recently popular GPUs and CPUs with multi- and many-core architectures. The scheme is based on the ray-tracing method and the computational cost is proportional to $N_{\rm m}^{5/3}$ where $N_{\rm m}$ is the number of mesh grids, and is devised to compute the radiation transfer along each light-ray completely in parallel with appropriate grouping of the light-rays. We find that the performance of our scheme scales well with the number of adopted CPU cores and GPUs, and also that our scheme is nicely parallelized on a multi-node system by adopting the multiple wave front scheme, and the performance scales well with the amount of the computational resources. As numerical tests to validate our scheme and to give a physical criterion for the angular resolution of our ray-tracing scheme, we perform several numerical simulations of the photo-ionization of neutral hydrogen gas by ionizing radiation sources without the "on-the-spot" approximation, in which the transfer of diffuse radiation by radiative recombination is incorporated in a self-consistent manner. △ Less

Submitted 16 April, 2015; v1 submitted 3 October, 2014; originally announced October 2014.

Comments: 26 pages, 15 figures, accepted for publication in PASJ

Direct Integration of the Collisionless Boltzmann Equation in Six-dimensional Phase Space: Self-gravitating Systems

Authors: Kohji Yoshikawa, Naoki Yoshida, Masayuki Umemura

Abstract: We present a scheme for numerical simulations of collisionless self-gravitating systems which directly integrates the Vlasov--Poisson equations in six-dimensional phase space. By the results from a suite of large-scale numerical simulations, we demonstrate that the present scheme can simulate collisionless self-gravitating systems properly. The integration scheme is based on the positive flux cons… ▽ More We present a scheme for numerical simulations of collisionless self-gravitating systems which directly integrates the Vlasov--Poisson equations in six-dimensional phase space. By the results from a suite of large-scale numerical simulations, we demonstrate that the present scheme can simulate collisionless self-gravitating systems properly. The integration scheme is based on the positive flux conservation method recently developed in plasma physics. We test the accuracy of our code by performing several test calculations including the stability of King spheres, the gravitational instability and the Landau damping. We show that the mass and the energy are accurately conserved for all the test cases we study. The results are in good agreement with linear theory predictions and/or analytic solutions. The distribution function keeps the property of positivity and remains non-oscillatory. The largest simulations are run on 64^6 grids. The computation speed scales well with the number of processors, and thus our code performs efficiently on massively parallel supercomputers. △ Less

Submitted 26 June, 2012; originally announced June 2012.

Comments: 35 pages, 19 figures. Submitted to the Astrophysical Journal

Imaging Simulations of the Sunyaev-Zel'dovich Effect for ALMA

Authors: K. Yamada, T. Kitayama, S. Takakuwa, D. Iono, T. Tsutsumi, K. Kohno, M. Takizawa, K. Yoshikawa, T. Akahori, E. Komatsu, Y. Suto, H. Matsuo, R. Kawabe

Abstract: We present imaging simulations of the Sunyaev-Zel'dovich effect of galaxy clusters for the Atacama Large Millimeter/submillimeter Array (ALMA) including the Atacama Compact Array (ACA). In its most compact configuration at 90GHz, ALMA will resolve the intracluster medium with an effective angular resolution of 5 arcsec. It will provide a unique probe of shock fronts and relativistic electrons prod… ▽ More We present imaging simulations of the Sunyaev-Zel'dovich effect of galaxy clusters for the Atacama Large Millimeter/submillimeter Array (ALMA) including the Atacama Compact Array (ACA). In its most compact configuration at 90GHz, ALMA will resolve the intracluster medium with an effective angular resolution of 5 arcsec. It will provide a unique probe of shock fronts and relativistic electrons produced during cluster mergers at high redshifts, that are hard to spatially resolve by current and near-future X-ray detectors. Quality of image reconstruction is poor with the 12m array alone but improved significantly by adding ACA; expected sensitivity of the 12m array based on the thermal noise is not valid for the Sunyaev-Zel'dovich effect mapping unless accompanied by an ACA observation of at least equal duration. The observations above 100 GHz will become excessively time-consuming owing to the narrower beam size and the higher system temperature. On the other hand, significant improvement of the observing efficiency is expected once Band 1 is implemented in the future. △ Less

Submitted 30 June, 2012; v1 submitted 27 March, 2012; originally announced March 2012.

Comments: 16 pages, 12 figures. Accepted for publication in PASJ. Note added in proof is included

Report number: TCC-003-12

Phantom-GRAPE: numerical software library to accelerate collisionless $N$-body simulation with SIMD instruction set on x86 architecture

Authors: Ataru Tanikawa, Kohji Yoshikawa, Keigo Nitadori, Takashi Okamoto

Abstract: (Abridged) We have developed a numerical software library for collisionless N-body simulations named "Phantom-GRAPE" which highly accelerates force calculations among particles by use of a new SIMD instruction set extension to the x86 architecture, AVX, an enhanced version of SSE. In our library, not only the Newton's forces, but also central forces with an arbitrary shape f(r), which has a finite… ▽ More (Abridged) We have developed a numerical software library for collisionless N-body simulations named "Phantom-GRAPE" which highly accelerates force calculations among particles by use of a new SIMD instruction set extension to the x86 architecture, AVX, an enhanced version of SSE. In our library, not only the Newton's forces, but also central forces with an arbitrary shape f(r), which has a finite cutoff radius r_cut (i.e. f(r)=0 at r>r_cut), can be quickly computed. Using an Intel Core i7--2600 processor, we measure the performance of our library for both the forces. In the case of Newton's forces, we achieve 2 x 10^9 interactions per second with 1 processor core, which is 20 times higher than the performance of an implementation without any explicit use of SIMD instructions, and 2 times than that with the SSE instructions. With 4 processor cores, we obtain the performance of 8 x 10^9 interactions per second. In the case of the arbitrarily shaped forces, we can calculate 1 x 10^9 and 4 x 10^9 interactions per second with 1 and 4 processor cores, respectively. The performance with 1 processor core is 6 times and 2 times higher than those of the implementations without any use of SIMD instructions and with the SSE instructions. These performances depend weakly on the number of particles. It is good contrast with the fact that the performance of force calculations accelerated by GPUs depends strongly on the number of particles. Substantially weak dependence of the performance on the number of particles is suitable to collisionless N-body simulations, since these simulations are usually performed with sophisticated N-body solvers such as Tree- and TreePM-methods combined with an individual timestep scheme. Collisionless N-body simulations accelerated with our library have significant advantage over those accelerated by GPUs, especially on massively parallel environments. △ Less

Submitted 9 October, 2012; v1 submitted 19 March, 2012; originally announced March 2012.

Comments: 19 pages, 11 figures, 4tables, accepted for publication in New Astronomy

Non-Equilibrium Ionization State and Two-Temperature Structure in the Bullet Cluster 1E0657-56

Authors: Takuya Akahori, Kohji Yoshikawa

Abstract: We investigate a non-equilibrium ionization state and an electron-ion two-temperature structure of the intracluster medium in the merging galaxy cluster, 1E0657-56 (the Bullet cluster), using a series of N-body and hydrodynamic simulations. We find that the electron temperature at the shock layer associated with the X-ray sub peak (bullet) is quite different depending on the thermal relaxation mod… ▽ More We investigate a non-equilibrium ionization state and an electron-ion two-temperature structure of the intracluster medium in the merging galaxy cluster, 1E0657-56 (the Bullet cluster), using a series of N-body and hydrodynamic simulations. We find that the electron temperature at the shock layer associated with the X-ray sub peak (bullet) is quite different depending on the thermal relaxation model between electrons and ions; ~25 keV for the Coulomb thermal relaxation model and ~45 keV for the instantaneous thermal relaxation model in the simulations which reproduce the observed X-ray morphology. Furthermore, both of Fe xxv and Fe xxvi are overabundant compared with the ionization equilibrium state around the shock layer, and thus, the intensity ratio between Fe xxv and Fe xxvi K alpha lines are significantly altered from that in the ionization equilibrium state. We also carry out the simulations with various sets of merger parameters, and discuss a possible range of the non-equilibrium effects in this system. Our results could be tested with future X-ray observations such as Astro-H with better sensitivity in high energy band. △ Less

Submitted 5 September, 2011; originally announced September 2011.

Comments: 11 pages, 9 figures. To appear in PASJ

A high-frequency study of the Sunyaev-Zel'dovich effect morphology in galaxy clusters

Authors: D. A. Prokhorov, S. Colafrancesco, T. Akahori, E. T. Million, S. Nagataki, K. Yoshikawa

Abstract: High-frequency, high-resolution imaging of the Sunyaev-Zel'dovich (SZ) effect is an important technique to study the complex structures of the atmospheres of merging galaxy clusters. Such observations are sensitive to the details of the electron spectrum. We show that the morphology of the SZ intensity maps in simulated galaxy clusters observed at 345 GHz, 600 GHz, and 857 GHz are significantly di… ▽ More High-frequency, high-resolution imaging of the Sunyaev-Zel'dovich (SZ) effect is an important technique to study the complex structures of the atmospheres of merging galaxy clusters. Such observations are sensitive to the details of the electron spectrum. We show that the morphology of the SZ intensity maps in simulated galaxy clusters observed at 345 GHz, 600 GHz, and 857 GHz are significantly different because of SZ relativistic corrections. These differences can be revealed by high-resolution imaging instruments. We calculate relativistically corrected SZ intensity maps of a simulated, massive, merging galaxy cluster and of the massive, merging clusters 1E0657-558 (the Bullet Cluster) and Abell 2219. The morphologies of the SZ intensity maps are remarkably different between 345 GHz and 857 GHz for each merging cluster. We show that high-resolution imaging observations of the SZ intensity maps at these frequencies, obtainable with the LABOCA and HERSCHEL-SPIRE instruments, allow to fully exploit the astrophysical relevance of the predicted SZ morphological effect. △ Less

Submitted 21 May, 2011; originally announced May 2011.

Comments: 9 pages, 14 figures, accepted for publication in MNRAS

Can electron distribution functions be derived through the Sunyaev-Zel'dovich effect?

Authors: D. A. Prokhorov, S. Colafrancesco, T. Akahori, K. Yoshikawa, S. Nagataki, K. -I. Seon

Abstract: Measurements of the Sunyaev-Zel'dovich (hereafter SZ) effect distortion of the cosmic microwave background provide methods to derive the gas pressure and temperature of galaxy clusters. Here we study the ability of SZ effect observations to derive the electron distribution function (DF) in massive galaxy clusters. Our calculations of the SZ effect include relativistic corrections considered with… ▽ More Measurements of the Sunyaev-Zel'dovich (hereafter SZ) effect distortion of the cosmic microwave background provide methods to derive the gas pressure and temperature of galaxy clusters. Here we study the ability of SZ effect observations to derive the electron distribution function (DF) in massive galaxy clusters. Our calculations of the SZ effect include relativistic corrections considered within the framework of the Wright formalism and use a decomposition technique of electron DFs into Fourier series. Using multi-frequency measurements of the SZ effect, we find the solution of a linear system of equations that is used to derive the Fourier coefficients; we further analyze different frequency samples to decrease uncertainties in Fourier coefficient estimations. We propose a method to derive DFs of electrons using SZ multi-frequency observations of massive galaxy clusters. We found that the best frequency sample to derive an electron DF includes high frequencies $ν$=375, 600, 700, 857 GHz. We show that it is possible to distinguish a Juttner DF from a Maxwell-Bolzman DF as well as from a Juttner DF with the second electron population by means of SZ observations for the best frequency sample if the precision of SZ intensity measurements is less than 0.1%. We demonstrate by means of 3D hydrodynamic numerical simulations of a hot merging galaxy cluster that the morphologies of SZ intensity maps are different for frequencies $ν$=375, 600, 700, 857 GHz. We stress that measurements of SZ intensities at these frequencies are a promising tool for studying electron distribution functions in galaxy clusters. △ Less

Submitted 17 May, 2011; originally announced May 2011.

Comments: 11 pages, 12 figures, published in Astronomy and Astrophysics

ARGOT: Accelerated radiative transfer on grids using oct-tree

Authors: Takashi Okamoto, Kohji Yoshikawa, Masayuki Umemura

Abstract: We present two types of numerical prescriptions that accelerate the radiative transfer calculation around point sources within a three-dimensional Cartesian grid by using the oct-tree structure for the distribution of radiation sources. In one prescription, distant radiation sources are grouped as a bright extended source when the group's angular size, $θ_{\rm s}$, is smaller than a critical value… ▽ More We present two types of numerical prescriptions that accelerate the radiative transfer calculation around point sources within a three-dimensional Cartesian grid by using the oct-tree structure for the distribution of radiation sources. In one prescription, distant radiation sources are grouped as a bright extended source when the group's angular size, $θ_{\rm s}$, is smaller than a critical value, $θ_{\rm crit}$, and radiative transfer is solved on supermeshes whose angular sizes are similar to that of the group of sources. The supermesh structure is constructed by coarse-graining the mesh structure. With this method, the computational time scales with $N_{\rm m} \log(N_{\rm m}) \log(N_{\rm s})$ where $N_{\rm m}$ and $N_{\rm s}$ are the number of meshes and that of radiation sources, respectively. While this method is very efficient, it inevitably overestimates the optical depth when a group of sources acts as an extended powerful radiation source and affects distant meshes. In the other prescription, a distant group of sources is treated as a bright point source ignoring the spatial extent of the group and the radiative transfer is solved on the meshes rather than the supermeshes. This prescription is simply a grid-based version of {\scriptsize START} by Hasegawa & Umemura and yields better results in general with slightly more computational cost ($\propto N_{\rm m}^{4/3} \log(N_{\rm s})$) than the supermesh prescription. Our methods can easily be implemented to any grid-based hydrodynamic codes and are well-suited to the adaptive mesh refinement methods. △ Less

Submitted 7 September, 2011; v1 submitted 16 May, 2011; originally announced May 2011.

Comments: 13 pages, 12 figures, submitted to MNRAS. Revised according to referee's comments

Journal ref: MNRAS, 419, 2855-2866, 2012

N-body simulation for self-gravitating collisional systems with a new SIMD instruction set extension to the x86 architecture, Advanced Vector eXtensions

Authors: Ataru Tanikawa, Kohji Yoshikawa, Takashi Okamoto, Keigo Nitadori

Abstract: We present a high-performance N-body code for self-gravitating collisional systems accelerated with the aid of a new SIMD instruction set extension of the x86 architecture: Advanced Vector eXtensions (AVX), an enhanced version of the Streaming SIMD Extensions (SSE). With one processor core of Intel Core i7-2600 processor (8 MB cache and 3.40 GHz) based on Sandy Bridge micro-architecture, we implem… ▽ More We present a high-performance N-body code for self-gravitating collisional systems accelerated with the aid of a new SIMD instruction set extension of the x86 architecture: Advanced Vector eXtensions (AVX), an enhanced version of the Streaming SIMD Extensions (SSE). With one processor core of Intel Core i7-2600 processor (8 MB cache and 3.40 GHz) based on Sandy Bridge micro-architecture, we implemented a fourth-order Hermite scheme with individual timestep scheme (Makino and Aarseth, 1992), and achieved the performance of 20 giga floating point number operations per second (GFLOPS) for double-precision accuracy, which is two times and five times higher than that of the previously developed code implemented with the SSE instructions (Nitadori et al., 2006b), and that of a code implemented without any explicit use of SIMD instructions with the same processor core, respectively. We have parallelized the code by using so-called NINJA scheme (Nitadori et al., 2006a), and achieved 90 GFLOPS for a system containing more than N = 8192 particles with 8 MPI processes on four cores. We expect to achieve about 10 tera FLOPS (TFLOPS) for a self-gravitating collisional system with N 105 on massively parallel systems with at most 800 cores with Sandy Bridge micro-architecture. This performance will be comparable to that of Graphic Processing Unit (GPU) cluster systems, such as the one with about 200 Tesla C1070 GPUs (Spurzem et al., 2010). This paper offers an alternative to collisional N-body simulations with GRAPEs and GPUs. △ Less

Submitted 5 September, 2011; v1 submitted 14 April, 2011; originally announced April 2011.

Comments: 14 pages, 9 figures, 3 tables, accepted for publication in New Astronomy. The code is publicly available at http://code.google.com/p/phantom-grape/

Unveiling the 3D temperature structure of galaxy clusters by means of the thermal SZ effect

Authors: D. A. Prokhorov, Y. Dubois, S. Nagataki, T. Akahori, K. Yoshikawa

Abstract: The Sunyaev-Zel'dovich (hereafter SZ) effect is a promising tool to derive the gas temperature of galaxy clusters. The approximation of a spherically symmetric gas distribution is usually used to determine the temperature structure of galaxy clusters, but this approximation cannot properly describe merging galaxy clusters. The methods used so far, which do not assume the spherically symmetric dist… ▽ More The Sunyaev-Zel'dovich (hereafter SZ) effect is a promising tool to derive the gas temperature of galaxy clusters. The approximation of a spherically symmetric gas distribution is usually used to determine the temperature structure of galaxy clusters, but this approximation cannot properly describe merging galaxy clusters. The methods used so far, which do not assume the spherically symmetric distribution, permit us to derive 2D temperature maps of merging galaxy clusters. In this paper, we propose a method to derive the standard temperature deviation and temperature variance along the line-of-sight, which permits us to analyze the 3D temperature structure of galaxy clusters by means of the thermal SZ effect. We also propose a method to reveal merger shock waves in galaxy clusters by analyzing the presence of temperature inhomogeneities along the line-of-sight. △ Less

Submitted 12 April, 2011; originally announced April 2011.

Comments: 8 pages, 3 figures, accepted for publication in MNRAS

ORIGIN: Metal Creation and Evolution from the Cosmic Dawn

Authors: J. W. den Herder, L. Piro, T. Ohashi, C. Kouveliotou, D. H. Hartmann, J. S. Kaastra, L. Amati, M. I. Andersen, M. Arnaud, J-L. Attéia, S. Bandler, M. Barbera, X. Barcons, S. Barthelmy, S. Basa, S. Basso, M. Boer, E. Branchini, G. Branduardi-Raymont, S. Borgani, A. Boyarsky, G. Brunetti, C. Budtz-Jorgensen, D. Burrows, N. Butler , et al. (134 additional authors not shown)

Abstract: ORIGIN is a proposal for the M3 mission call of ESA aimed at the study of metal creation from the epoch of cosmic dawn. Using high-spectral resolution in the soft X-ray band, ORIGIN will be able to identify the physical conditions of all abundant elements between C and Ni to red-shifts of z=10, and beyond. The mission will answer questions such as: When were the first metals created? How does the… ▽ More ORIGIN is a proposal for the M3 mission call of ESA aimed at the study of metal creation from the epoch of cosmic dawn. Using high-spectral resolution in the soft X-ray band, ORIGIN will be able to identify the physical conditions of all abundant elements between C and Ni to red-shifts of z=10, and beyond. The mission will answer questions such as: When were the first metals created? How does the cosmic metal content evolve? Where do most of the metals reside in the Universe? What is the role of metals in structure formation and evolution? To reach out to the early Universe ORIGIN will use Gamma-Ray Bursts (GRBs) to study their local environments in their host galaxies. This requires the capability to slew the satellite in less than a minute to the GRB location. By studying the chemical composition and properties of clusters of galaxies we can extend the range of exploration to lower redshifts (z ~ 0.2). For this task we need a high-resolution spectral imaging instrument with a large field of view. Using the same instrument, we can also study the so far only partially detected baryons in the Warm-Hot Intergalactic Medium (WHIM). The less dense part of the WHIM will be studied using absorption lines at low redshift in the spectra for GRBs. △ Less

Submitted 15 April, 2011; v1 submitted 11 April, 2011; originally announced April 2011.

Comments: 34 pages, 13 figures. ESA Cosmic Vision medium-class mission (M3) proposal. Accepted for publication in Experimental Astronomy. Including minor corrections in the author list

Journal ref: Exp Astron (2012) 34, 519

Hydrodynamic Simulations of Merging Galaxy Clusters: Non-Equilibrium Ionization State and Two-Temperature Structure

Authors: Takuya Akahori, Kohji Yoshikawa

Abstract: We investigate a non-equilibrium ionization state and an electron-ion two-temperature structure of the intracluster medium (ICM) in merging galaxy clusters using a series of N-body and hydrodynamic simulations. Mergers with various sets of mass ratios and impact parameters are systematically investigated, and it is found that, in most cases, ICM significantly departs from the ionization equilibr… ▽ More We investigate a non-equilibrium ionization state and an electron-ion two-temperature structure of the intracluster medium (ICM) in merging galaxy clusters using a series of N-body and hydrodynamic simulations. Mergers with various sets of mass ratios and impact parameters are systematically investigated, and it is found that, in most cases, ICM significantly departs from the ionization equilibrium state at the shock layers with a Mach number of ~1.5-2.0 in the outskirts of the clusters, and the shock layers with a Mach number of ~2-4 in front of the ICM cores. Accordingly, the intensity ratio between Fe xxv and Fe xxvi K alpha line emissions is significantly altered from that in the ionization equilibrium state. If the effect of the two-temperature structure of ICM is incorporated, the electron temperature is ~10-20 % and ~30-50 % lower than the mean temperature of ICM at the shock layers in the outskirts and in front of the ICM cores, respectively, and the deviation from the ionization equilibrium state becomes larger. We also address the dependence of the intensity ratio on the viewing angle with respect to the merging plane. △ Less

Submitted 24 January, 2010; v1 submitted 28 September, 2009; originally announced September 2009.

Comments: 11 pages, 10 figures. Submitted to PASJ; Accepted for publication in PASJ

Modeling Nonlinear Evolution of Baryon Acoustic Oscillations: Convergence Regime of N-body Simulations and Analytic Models

Authors: Takahiro Nishimichi, Akihito Shirata, Atsushi Taruya, Kazuhiro Yahata, Shun Saito, Yasushi Suto, Ryuichi Takahashi, Naoki Yoshida, Takahiko Matsubara, Naoshi Sugiyama, Issha Kayo, Yipeng Jing, Kohji Yoshikawa

Abstract: We use a series of cosmological N-body simulations and various analytic models to study the evolution of the matter power spectrum in real space in a ΛCold Dark Matter universe. We compare the results of N-body simulations against three analytical model predictions; standard perturbation theory, renormalized perturbation theory, and the closure approximation. We include the effects from finite s… ▽ More We use a series of cosmological N-body simulations and various analytic models to study the evolution of the matter power spectrum in real space in a ΛCold Dark Matter universe. We compare the results of N-body simulations against three analytical model predictions; standard perturbation theory, renormalized perturbation theory, and the closure approximation. We include the effects from finite simulation box size in the comparison. We determine the values of the maximum wavenumbers, k^{lim}_{1%} and k^{lim}_{3%}, below which the analytic models and the simulation results agree to within 1 and 3 percent, respectively. We then provide a simple empirical function which describes the convergence regime determined by comparison between our simulations and the analytical models. We find that if we use the Fourier modes within the convergence regime alone, the characteristic scale of baryon acoustic oscillations can be determined within 1% accuracy from future surveys with a volume of a few h^{-3}Gpc^3 at z\sim1 or z\sim3 in the absence of any systematic distortion of the power spectrum. △ Less

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

Comments: 23 pages, 14 figures, accepted by PASJ

Journal ref: Publ. Astron. Soc. Japan 61, pp.321-332 (2009)

Non-Equilibrium Ionization State and Two-Temperature Structure in the Linked Region of Abell 399/401

Authors: Takuya Akahori, Kohji Yoshikawa

Abstract: We investigate a non-equilibrium ionization state and two-temperature structure of the intracluster medium in the linked region of Abell 399/401, using a series of N-body + SPH simulations, and find that there exist significant shock layers at the edge of the linked region, and that the ionization state of iron departs from the ionization equilibrium state at the shock layers and around the cent… ▽ More We investigate a non-equilibrium ionization state and two-temperature structure of the intracluster medium in the linked region of Abell 399/401, using a series of N-body + SPH simulations, and find that there exist significant shock layers at the edge of the linked region, and that the ionization state of iron departs from the ionization equilibrium state at the shock layers and around the center of the linked region. As for the two-temperature structure, an obvious difference of temperature between electrons and ions is found in the edge of the linked regions. K alpha line emissions of Fe xxiv and Fe xxv are not severely affected by the deviation from the ionization equilibrium state around the center of the linked region, suggesting that the detection of relatively high metallicity in this area cannot be ascribed to the non-equilibrium ionization state of the intracluster medium. On the other hand, the K alpha emissions are significantly deviated from the equilibrium values at the shock layers, and the intensity ratio of K alpha lines between Fe xxiv-xxv and Fe xxvi is found to be significantly altered from that in the ionization equilibrium state. △ Less

Submitted 5 June, 2008; v1 submitted 16 May, 2008; originally announced May 2008.

Comments: 4 pages, 2 figures, submitted to PASJ Letter; accepted June 2

Journal ref: PASJ vol.60, No.4, L19-22 (2008)

Suzaku broad-band spectroscopy of RX J1347.5-1145: constraints on the extremely hot gas and non-thermal emission

Authors: N. Ota, K. Murase, T. Kitayama, E. Komatsu, M. Hattori, H. Matsuo, T. Oshima, Y. Suto, K. Yoshikawa

Abstract: We present the results from the analysis of long Suzaku observations of the most X-ray luminous galaxy cluster RX J1347.5-1145 at z=0.451. Aims: We study physical properties of the hot (~20 keV) gas clump in the south-east (SE) region discovered by the Sunyaev-Zel'dovich (SZ) effect observations, to understand the gas physics of a violent cluster merger. We also explore a signature of non-therma… ▽ More We present the results from the analysis of long Suzaku observations of the most X-ray luminous galaxy cluster RX J1347.5-1145 at z=0.451. Aims: We study physical properties of the hot (~20 keV) gas clump in the south-east (SE) region discovered by the Sunyaev-Zel'dovich (SZ) effect observations, to understand the gas physics of a violent cluster merger. We also explore a signature of non-thermal emission using the hard X-ray data. Results: We find that the single-temperature model fails to reproduce the continuum emission and Fe-K lines measured by XIS simultaneously. The two-temperature model with a very hot component improves the fit, although the XIS data can only give a lower bound on its temperature. We detect the hard X-ray emission in the 12-40 keV band at the 7 sigma level; however, the significance becomes marginal when the systematic error in the background estimation is included. With the Suzaku + Chandra joint analysis, we determine the temperature of the SE excess component to be 25.3^{+6.1}_{-4.5} ^{+6.9}_{-9.5} keV (90% statistical and systematic errors), which is in an excellent agreement with the previous SZ + X-ray analysis. This is the first time that the X-ray spectroscopy alone gives a good measurement of the temperature of the hot component in the SE region, which is made possible by Suzaku's unprecedented sensitivity to the wide X-ray band. These results strongly indicate that the cluster has undergone a recent, violent merger. The spectral analysis shows that the SE component is consistent with being thermal. We find the 3 sigma upper limit on the non-thermal flux, F < 8e-12 erg s^{-1} cm^{-2} in the 12-60 keV band. Combining this limit with a recent discovery of the radio mini halo at 1.4 GHz, we find a lower limit on the strength of the intracluster magnetic field, B > 0.007 micro G. △ Less

Submitted 1 September, 2008; v1 submitted 5 May, 2008; originally announced May 2008.

Comments: 15 pages, 13 figures. Accepted for publication in A&A

DIOS: the Diffuse Intergalactic Oxygen Surveyor

Authors: T. Ohashi, M. Ishida, S. Sasaki, Y. Ishisaki, K. Mitsuda, N. Y. Yamasaki, R. Fujimoto, Y. Takei, Y. Tawara, A. Furuzawa, Y. Suto, K. Yoshikawa, H. Kawahara, N. Kawai, T. G. Tsuru, K. Matsushita, T. Kitayama

Abstract: We present our proposal for a small X-ray mission DIOS (Diffuse Intergalactic Oxygen Surveyor), consisting of a 4-stage X-ray telescope and an array of TES microcalorimeters, cooled with mechanical coolers, with a total weight of about 400 kg. The mission will perform survey observations of warm-hot intergalactic medium using OVII and OVIII emission lines, with the energy coverage up to 1.5 keV.… ▽ More We present our proposal for a small X-ray mission DIOS (Diffuse Intergalactic Oxygen Surveyor), consisting of a 4-stage X-ray telescope and an array of TES microcalorimeters, cooled with mechanical coolers, with a total weight of about 400 kg. The mission will perform survey observations of warm-hot intergalactic medium using OVII and OVIII emission lines, with the energy coverage up to 1.5 keV. The wide field of view of about 50' diameter, superior energy resolution close to 2 eV FWHM, and very low background will together enable us a wide range of science for diffuse X-ray sources. We briefly describe the design of the satellite, performance of the subsystems and the expected results. △ Less

Submitted 7 July, 2006; originally announced July 2006.

Comments: 9 pages, 11 figures, a proceedings of SPIE "Astronomical Telescopes and Instrumentation" 2006

Proto-clusters in the Lambda CDM Universe

Authors: Tamon Suwa, Asao Habe, Kohji Yoshikawa

Abstract: We compare the highly clustered populations of very high redshift galaxies with proto-clusters identified numerically in a standard $Λ$CDM universe ($Ω_0=0.3, λ_0=0.7$) simulation. We evolve 256^3 dark matter particles in a comoving box of side 150h^{-1}Mpc. By the present day there are 63 cluster sized objects of mass in excess of 10^{14}h^{-1}Mo in this box. We trace these clusters back to hig… ▽ More We compare the highly clustered populations of very high redshift galaxies with proto-clusters identified numerically in a standard $Λ$CDM universe ($Ω_0=0.3, λ_0=0.7$) simulation. We evolve 256^3 dark matter particles in a comoving box of side 150h^{-1}Mpc. By the present day there are 63 cluster sized objects of mass in excess of 10^{14}h^{-1}Mo in this box. We trace these clusters back to higher redshift finding that their progenitors at z=4--5 are extended regions of typically 20--40 Mpc (comoving) in size, with dark halos of mass in excess of 10^{12}h^{-1}Mo and are overdense by typically 1.3--13 times the cosmological mean density. Comparison with the observation of Lyman alpha emitting (LAEs) galaxies at z=4.86 and at z=4.1 indicates that the observed excess clustering is consistent with that expected for a proto-cluster region if LAEs typically correspond to massive dark halos of more than 10^{12}h^{-1}Mo. We give a brief discussion on the relation between high redshift concentration of massive dark halos and present day rich clusters of galaxies. △ Less

Submitted 15 June, 2006; v1 submitted 13 June, 2006; originally announced June 2006.

Comments: 4 pages, 5 figures, Accepted for publication in ApJ Letters

Journal ref: Astrophys.J.646:L5-L8,2006

Non-equilibrium Ionization State of Warm-Hot Intergalactic Medium

Authors: Kohji Yoshikawa, Shin Sasaki

Abstract: Time evolution of the ionization state of metals in the cosmic baryons is investigated in a cosmological context without the assumption of ionization equilibrium. We find that a significant fraction of ionized oxygen ions (OVII and OVIII) in the warm-hot intergalactic medium (WHIM) is not in the ionization equilibrium state at a redshift of z~0. We also investigate the effect on the detectabilit… ▽ More Time evolution of the ionization state of metals in the cosmic baryons is investigated in a cosmological context without the assumption of ionization equilibrium. We find that a significant fraction of ionized oxygen ions (OVII and OVIII) in the warm-hot intergalactic medium (WHIM) is not in the ionization equilibrium state at a redshift of z~0. We also investigate the effect on the detectability and observables of WHIM as a consequence of such deviation from ionization equilibrium. It is found that the detectability of WHIM is not altered very much both through its emission and absorption signatures, but line ratios between OVII and OVIII are significantly different from those in the ionization equilibrium state. △ Less

Submitted 13 June, 2006; v1 submitted 2 March, 2006; originally announced March 2006.

Comments: 25 pages, 18 figures, accepted for publication in PASJ

Soft X-ray Transmission Spectroscopy of Warm/Hot Intergalactic Medium with XEUS

Authors: Hajime Kawahara, Kohji Yoshikawa, Shin Sasaki, Yasushi Suto, Nobuyuki Kawai, Kazuhisa Mitsuda, Takaya Ohashi, Noriko Yamasaki

Abstract: We discuss the detectability of Warm/Hot Intergalactic medium (WHIM) via the absorption lines toward bright point sources with a future X-ray satellite mission, XEUS. While we consider bright QSOs as specific examples, the methodology can be applied to bright GRB afterglows. We create mock absorption spectra for bright QSOs (more than 20 QSOs over the all sky) using a light-cone output of a cosm… ▽ More We discuss the detectability of Warm/Hot Intergalactic medium (WHIM) via the absorption lines toward bright point sources with a future X-ray satellite mission, XEUS. While we consider bright QSOs as specific examples, the methodology can be applied to bright GRB afterglows. We create mock absorption spectra for bright QSOs (more than 20 QSOs over the all sky) using a light-cone output of a cosmological hydrodynamic simulation. We assume that WHIM is under collisional and photo-ionization equilibrium. If WHIM has a constant metallicity of $Z=0.1Z_\odot$, approximately 2 O{\sc vii} absorption line system with $>3σ$ will be detected on average along a random line-of-sight toward bright QSOs up to $z=0.3$ for 30 ksec exposure. △ Less

Submitted 18 July, 2006; v1 submitted 27 April, 2005; originally announced April 2005.

Comments: 34 pages, 17 figures. Substantial revision. PASJ, in press (2006, vol.58, August 25 issue)

PPPM and TreePM Methods on GRAPE Systems for Cosmological N-body Simulations

Authors: Kohji Yoshikawa, Toshiyuki Fukushige

Abstract: We present Particle-Particle-Particle-Mesh (PPPM) and Tree Particle-Mesh (TreePM) implementations on GRAPE-5 and GRAPE-6A systems, special-purpose hardware accelerators for gravitational many-body simulations. In our PPPM and TreePM implementations on GRAPE, the computational time is significantly reduced compared with the conventional implementations without GRAPE, especially under the strong p… ▽ More We present Particle-Particle-Particle-Mesh (PPPM) and Tree Particle-Mesh (TreePM) implementations on GRAPE-5 and GRAPE-6A systems, special-purpose hardware accelerators for gravitational many-body simulations. In our PPPM and TreePM implementations on GRAPE, the computational time is significantly reduced compared with the conventional implementations without GRAPE, especially under the strong particle clustering, and almost constant irrespective of the degree of particle clustering. We carry out the survey of two simulation parameters, the PM grid spacing and the opening parameter for the most optimal combination of force accuracy and computational speed. We also describe the parallelization of these implementations on a PC-GRAPE cluster, in which each node has one GRAPE board, and present the optimal configuration of simulation parameters for good parallel scalability. △ Less

Submitted 12 September, 2005; v1 submitted 4 April, 2005; originally announced April 2005.

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

Journal ref: Publ.Astron.Soc.Jap. 57 (2005) 849-860

Locating the Warm-Hot Intergalactic Medium in the Simulated Local Universe

Authors: K. Yoshikawa, K. Dolag, Y. Suto, S. Sasaki, N. Y. Yamasaki, T. Ohashi, K. Mitsuda, Y. Tawara, R. Fujimoto, T. Furusho, A. Furuzawa, M. Ishida, Y. Ishisaki, Y. Takei

Abstract: We present an analysis of mock spectral observation of warm-hot intergalactic medium (WHIM) using a constrained simulation of the local universe. The simulated map of oxygen emission lines from local WHIM reproduces well the observed structures traced by galaxies in the real local universe. We further attempt to perform mock observations of outer parts of simulated Coma cluster and A3627 adoptin… ▽ More We present an analysis of mock spectral observation of warm-hot intergalactic medium (WHIM) using a constrained simulation of the local universe. The simulated map of oxygen emission lines from local WHIM reproduces well the observed structures traced by galaxies in the real local universe. We further attempt to perform mock observations of outer parts of simulated Coma cluster and A3627 adopting the expected performance of DIOS (Diffuse Intergalactic Oxygen Surveyor), which is proposed as a dedicated soft X-ray mission to search for cosmic missing baryons. We find that WHIMs surrounding nearby clusters are detectable with a typical exposure time of a day, and thus constitute realistic and promising targets for DIOS. We also find that an X-ray emitting clump in front of Coma cluster, recently reported in the XMM-Newton observation, has a counterpart in the simulated local universe, and its observed spectrum can be well reproduced in the simulated local universe if the gas temperature is set to the observationally estimated value. △ Less

Submitted 27 September, 2004; v1 submitted 7 August, 2004; originally announced August 2004.

Comments: 25 pages, 3 tables, 16 figures. Accepted for publication in PASJ. High resolution PS/PDF files are available at http://www-utap.phys.s.u-tokyo.ac.jp/~kohji/research/x-ray/index.html

Journal ref: Publ.Astron.Soc.Jap. 56 (2004) 939-957

Future Japanese missions for the study of warm-hot intergalactic medium

Authors: T. Ohashi, M. Ishida, S. Sasaki, Y. Ishisaki, K. Mitsuda, N. Y. Yamasaki, R. Fujimoto, T. Furusho, H. Kunieda, Y. Tawara, A. Furuzawa, Y. Suto, K. Yoshikawa

Abstract: We present our proposal for a small X-ray mission DIOS (Diffuse Intergalactic Oxygen Surveyor) to perform survey observations of warm-hot intergalactic medium using OVII and OVIII emission lines. This will be proposed to a small satellite program planned by ISAS/JAXA in Japan for a launch in 2008. The instrument consists of an array of TES microcalorimeters with an energy resolution 2 eV, cooled… ▽ More We present our proposal for a small X-ray mission DIOS (Diffuse Intergalactic Oxygen Surveyor) to perform survey observations of warm-hot intergalactic medium using OVII and OVIII emission lines. This will be proposed to a small satellite program planned by ISAS/JAXA in Japan for a launch in 2008. The instrument consists of an array of TES microcalorimeters with an energy resolution 2 eV, cooled by mechanical coolers. The X-ray telescope will employ 4-stage reflection mirrors with a focal length as short as 70 cm and an angular resolution 2'. In addition to DIOS, we briefly describe the NeXT (New X-ray Telescope) mission, which is a larger Japanese X-ray observatory to be launched in 2010 and plans to explore non-thermal processes in the universe. △ Less

Submitted 23 February, 2004; originally announced February 2004.

Comments: 6 pages, 9 figures. To appear in the proceedings of "Modelling the Intergalactic and Intracluster Media", Vulcano Island, October 1-4, 2003

Searching for cosmic missing baryons with DIOS -- Diffuse Intergalactic Oxygen Surveyor --

Authors: Yasushi Suto, Kohji Yoshikawa, Noriko Y. Yamasaki, Kazuhisa Mitsuda, Ryuichi Fujimoto, Tae Furusho, Takaya Ohashi, Manabu Ishida, Shin Sasaki, Yoshitaka Ishisaki, Yuzuru Tawara, Akihiro Furuzawa

Abstract: Approximately 30 to 50 percent of the total baryons in the present universe is supposed to take a form of warm/hot intergalactic medium (WHIM) whose X-ray continuum emission is very weak. In order to carry out a direct and homogeneous survey of elusive cosmic missing baryons, we propose a dedicated soft-X-ray mission, {\it DIOS} (Diffuse Intergalactic Oxygen Surveyor). The unprecedented energy r… ▽ More Approximately 30 to 50 percent of the total baryons in the present universe is supposed to take a form of warm/hot intergalactic medium (WHIM) whose X-ray continuum emission is very weak. In order to carry out a direct and homogeneous survey of elusive cosmic missing baryons, we propose a dedicated soft-X-ray mission, {\it DIOS} (Diffuse Intergalactic Oxygen Surveyor). The unprecedented energy resolution ($\sim 2$eV) of the XSA (X-ray Spectrometer Array) on-board {\it DIOS} enables us to identify WHIM with gas temperature $T=10^{6-7}$K and overdensity $δ=10-100$ located at $z<0.3$ through emission lines of O{\sc vii} and O{\sc viii}. {\it DIOS}, hopefully launched in several years, promises to open a new window of detection and characterization of cosmic missing baryons, and to provide yet another important and complementary tool to trace the large-scale structure of the universe. △ Less

Submitted 16 February, 2004; originally announced February 2004.

Comments: 8 pages, 10 figures, invited talk at the VI International Conference on Gravitation and Astrophysics of Asian-Pacfic Countries, to be published in the proceedings

Exploring Cluster Physics with High-Resolution Sunyaev-Zel'dovich Effect Images and X-Ray Data: The Case of the Most X-Ray-Luminous Galaxy Cluster RX J1347-1145

Authors: T. Kitayama, E. Komatsu, N. Ota, T. Kuwabara, Y. Suto, K. Yoshikawa, M. Hattori, H. Matsuo

Abstract: Foreseeing the era of high spatial resolution measurements of the Sunyaev-Zel'dovich effect (SZE) in clusters of galaxies, we present a prototype analysis of this sort combined with Chandra X-ray data. It is applied specifically to RX J1347-1145 at z=0.451, the most X-ray-luminous galaxy cluster known, for which the highest resolution SZE and X-ray images are currently available. We demonstrate… ▽ More Foreseeing the era of high spatial resolution measurements of the Sunyaev-Zel'dovich effect (SZE) in clusters of galaxies, we present a prototype analysis of this sort combined with Chandra X-ray data. It is applied specifically to RX J1347-1145 at z=0.451, the most X-ray-luminous galaxy cluster known, for which the highest resolution SZE and X-ray images are currently available. We demonstrate that the combined analysis yields a unique probe of complex structures in the intracluster medium, offering determinations of their temperature, density, and line-of-sight extent. For a subclump in RX J1347-1145, previously discovered in our SZE map, the temperature inferred after removing the foreground and background components is well in excess of 20 keV, indicating that the cluster has recently undergone a violent merger. Excluding the region around this subclump, the SZE signals in submillimeter to centimeter bands (350, 150, and 21 GHz) are all consistent with those expected from Chandra X-ray observations. We further present a temperature deprojection technique based on the SZE and X-ray images, without any knowledge of spatially resolved X-ray spectroscopy. The methodology presented here will be applicable to a statistical sample of clusters available in the future SZE surveys. △ Less

Submitted 18 February, 2004; v1 submitted 26 November, 2003; originally announced November 2003.

Comments: 13 pages, 8 figures; minor corrections, final version to appear in PASJ

Journal ref: Publ.Astron.Soc.Jap. 56 (2004) 17-28

Angular momentum distribution of hot gas and implications for disk galaxy formation

Authors: D. N. Chen, Y. P. Jing, Kohji Yoshikawa

Abstract: We study the angular momentum profiles both for dark matter and for gas within virialized halos, using a statistical sample of halos drawn from cosmological hydrodynamics simulations. Three simulations have been analyzed, one is the ``non-radiative'' simulation, and the other two have radiative cooling. We find that the gas component on average has a larger spin and contains a smaller fraction o… ▽ More We study the angular momentum profiles both for dark matter and for gas within virialized halos, using a statistical sample of halos drawn from cosmological hydrodynamics simulations. Three simulations have been analyzed, one is the ``non-radiative'' simulation, and the other two have radiative cooling. We find that the gas component on average has a larger spin and contains a smaller fraction of mass with negative angular momentum than its dark matter counterpart in the non-radiative model. As to the cooling models, the gas component shares approximately the same spin parameter as its dark matter counterpart, but the hot gas has a higher spin and is more aligned in angular momentum than dark matter, while the opposite holds for the cold gas. After the mass of negative angular momentum is excluded, the angular momentum profile of the hot gas component approximately follows the universal function originally proposed by Bullock et al. for dark matter, though the shape parameter $μ$ is much larger for hot gas and is comfortably in the range required by observations of disk galaxies. Since disk formation is related to the distribution of hot gas that will cool, our study may explain the fact that the disk component of observed galaxies contains a smaller fraction of low angular momentum material than dark matter in halos. △ Less

Submitted 13 July, 2003; originally announced July 2003.

Comments: 30 pages, 12 figures, 4 tables, accepted for publication in ApJ

Journal ref: Astrophys.J.597:35-47,2003

Detectability of the Warm/Hot Intergalactic Medium Through Emission Lines of OVII and OVIII

Authors: Kohji Yoshikawa, Noriko Y. Yamasaki, Yasushi Suto, Takaya Ohashi, Kazuhisa Mitsuda, Yuzuru Tawara, Akihiro Furuzawa

Abstract: Most of cosmic baryons predicted by the big-bang nucleosynthesis has evaded the direct detection. Recent numerical simulations indicate that approximately 30 to 50 percent of the total baryons in the present universe is supposed to take a form of warm/hot intergalactic medium (WHIM) whose X-ray continuum emission is very weak. To identify those missing baryons, we consider in detail the detectab… ▽ More Most of cosmic baryons predicted by the big-bang nucleosynthesis has evaded the direct detection. Recent numerical simulations indicate that approximately 30 to 50 percent of the total baryons in the present universe is supposed to take a form of warm/hot intergalactic medium (WHIM) whose X-ray continuum emission is very weak. To identify those missing baryons, we consider in detail the detectability of WHIM directly through emission lines of OVII (561, 568, 574, 665eV) and OVIII (653eV). For this purpose, we create mock spectra of the emission lines of WHIM using a light-cone output of the cosmological hydrodynamic simulations. Since the predicted fluxes are generally below the current detection limit, the unambiguous detection requires a dedicated X-ray satellite mission that we also discuss in detail. We find that our proposed mission is sensitive to the WHIM with gas temperature $T=10^{6-7}$K and overdensity $δ=10-100$ up to a redshift of 0.3 without being significantly contaminated by the cosmic X-ray background and the Galactic emissions. Thus such a mission provides a unique and important tool to identify a large fraction of otherwise elusive baryons in the universe. △ Less

Submitted 1 July, 2003; v1 submitted 13 March, 2003; originally announced March 2003.

Comments: 22 pages, 10 figures. To appear in PASJ, v55, Oct. 25 (2003) issue

Journal ref: Publ.Astron.Soc.Jap. 55 (2003) 879-890

The Luminosity Function of IRAS PSCz Galaxies

Authors: Tsutomu T. Takeuchi, Kohji Yoshikawa, Takako T. Ishii

Abstract: We estimated the luminosity function (LF) of IRAS galaxies in the PSCz catalogue. The faint end of the PSCz LF is slightly steeper than that of the LF derived by Saunders et al. (1990; S90). Using an analytical form for the LF used by S90, we obtain the following parameters: α= 1.23 \pm 0.04, L_*=(8.85 \pm 1.75) \times 10^8 h^{-2} L_\odot, σ=0.724 \pm 0.010, and φ_* = (2.34 \pm 0.30) \times 10^{… ▽ More We estimated the luminosity function (LF) of IRAS galaxies in the PSCz catalogue. The faint end of the PSCz LF is slightly steeper than that of the LF derived by Saunders et al. (1990; S90). Using an analytical form for the LF used by S90, we obtain the following parameters: α= 1.23 \pm 0.04, L_*=(8.85 \pm 1.75) \times 10^8 h^{-2} L_\odot, σ=0.724 \pm 0.010, and φ_* = (2.34 \pm 0.30) \times 10^{-2} h^3 Mpc^{-3}. We also examined the evolution in the sample by a simple assumption φ_*(z) \propto (1+z)^P, and found P=3.40 \pm 0.70. It does not affect the three parameters, α, L_*, and σ, but φ_*(z=0) is overestimated up to \sim 15% if we ignore evolution. We estimated the temperature dependence of the LF. The LFs of warm and cool galaxies are quite different: the LF of warm galaxies has a very steep faint end with α=1.37. We also discuss a lump found at the brightest end of the LF. △ Less

Submitted 8 March, 2003; originally announced March 2003.

Comments: ApJ Letters (part II) in press, 4 pages, 3 figures, using aastex.cls, emulateapj5.sty and apjfonts.sty

Journal ref: Astrophys.J. 587 (2003) L89-L92; Erratum-ibid. 606 (2004) L171-L172

Spatial and Dynamical Biases in Velocity Statistics of Galaxies

Authors: Kohji Yoshikawa, Y. P. Jing, Gerhard Börner

Abstract: We present velocity statistics of galaxies and their biases inferred from the statistics of the underlying dark matter using a cosmological hydrodynamic simulation of galaxy formation in low-density and spatially flat cold dark matter cosmogony. We find that the pairwise velocity dispersion (PVD) of all galaxies is significantly lower than that of the dark matter particles, and that the PVD of t… ▽ More We present velocity statistics of galaxies and their biases inferred from the statistics of the underlying dark matter using a cosmological hydrodynamic simulation of galaxy formation in low-density and spatially flat cold dark matter cosmogony. We find that the pairwise velocity dispersion (PVD) of all galaxies is significantly lower than that of the dark matter particles, and that the PVD of the young galaxies is lower than that of the old types, and even of all galaxies together, especially at small separations. These results are in reasonable agreement with the recent measurements of PVDs in the Las Campanas redshift survey, the PSCz catalogue and the SDSS data. We also find that the low PVD of young galaxies is due to the effects of dynamical friction as well as the different spatial distribution. We also consider the mean infall velocity and the POTENT density reconstruction that are often used to measure the cosmological parameters, and investigate the effects of spatial bias and dynamical friction. In our simulation, the mean infall velocity of young galaxies is significantly lower than that of all the galaxies or of the old galaxies, and the dynamical bias becomes important on scales less than 3Mpc/h. The mass density field reconstructed from the velocity field of young galaxies using the POTENT-style method suffers in accuracy both from the spatial bias and the dynamical friction on the smoothing scale of R_s=8Mpc/h. On the other hand, in the case of R_s=12Mpc/h, which is typically adopted in the actual POTENT analysis, the density reconstruction based on various tracers of galaxies is reasonably accurate. △ Less

Submitted 3 March, 2003; originally announced March 2003.

Comments: 29 pages, 11 figures, accepted for publication in the ApJ

Journal ref: Astrophys.J. 590 (2003) 654-663

The UV (GALEX) and FIR (ASTRO-F) All Sky Surveys: the measure of the dust extinction in the local universe

Authors: V. Buat, T. T. Takeuchi, A. Boselli, D. Burgarella, H. Hirashita, A. Tomita, H. Shibai, B. Milliard, J. Donas, K. Yoshikawa, A. K. Inoue, Y. Y. Tajiri

Abstract: Before the end of 2002 will be launched the GALEX satellite (a NASA/SMEX project) which will observe all the sky in Ultraviolet (UV) through filters at 1500 and 2300 A down to m(AB) 21. In 2004 will be launched the ASTRO-F satellite which will perform an all sky survey at Far-Infrared (FIR) wavelengths. The cross-correlation of both surveys will lead to very large samples of galaxies for which… ▽ More Before the end of 2002 will be launched the GALEX satellite (a NASA/SMEX project) which will observe all the sky in Ultraviolet (UV) through filters at 1500 and 2300 A down to m(AB) 21. In 2004 will be launched the ASTRO-F satellite which will perform an all sky survey at Far-Infrared (FIR) wavelengths. The cross-correlation of both surveys will lead to very large samples of galaxies for which FIR and UV fluxes will be available. Using the FIR to UV flux ratio as a quantitative tracer of the dust extinction we will be able to measure the extinction in the nearby universe (z<0.2) and to perform a statistically significant analysis of the extinction as a function of galactic properties. Of particular interest is the construction of pure FIR and UV selected samples for which the extinction will be measured as templates for the observation of high redshift galaxies. △ Less

Submitted 12 September, 2002; originally announced September 2002.

Journal ref: Proceedings of the 8th IAU Asian-Pacific Regional meeting, july 2002

Cluster Morphologies as a Test of Different Cosmological Models

Authors: Tamon Suwa, Asao Habe, Kohji Yoshikawa, Takashi Okamoto

Abstract: We investigate how cluster morphology is affected by the cosmological constant in low-density universes. Using high-resolution cosmological N-body/SPH simulations of flat (Ω_0 = 0.3, λ_0 = 0.7, ΛCDM) and open (Ω_0 = 0.3, λ_0 = 0, OCDM) cold dark matter universes, we calculate statistical indicators to quantify the irregularity of the cluster morphologies. We study axial ratios, center shifts, cl… ▽ More We investigate how cluster morphology is affected by the cosmological constant in low-density universes. Using high-resolution cosmological N-body/SPH simulations of flat (Ω_0 = 0.3, λ_0 = 0.7, ΛCDM) and open (Ω_0 = 0.3, λ_0 = 0, OCDM) cold dark matter universes, we calculate statistical indicators to quantify the irregularity of the cluster morphologies. We study axial ratios, center shifts, cluster clumpiness, and multipole moment power ratios as indicators for the simulated clusters at z=0 and 0.5. Some of these indicators are calculated for both the X-ray surface brightness and projected mass distributions. In ΛCDM all these indicators tend to be larger than those in OCDM at z=0. This result is consistent with the analytical prediction of Richstone, Loeb, & Turner, that is, clusters in ΛCDM are formed later than in OCDM, and have more substructure at z=0. We make a Kolmogorov-Smirnov test on each indicator for these two models. We then find that the results for the multipole moment power ratios and the center shifts for the X-ray surface brightness are under the significance level (5%). We results also show that these two cosmological models can be distinguished more clearly at z=0 than z = 0.5 by these indicators. △ Less

Submitted 9 January, 2003; v1 submitted 19 August, 2001; originally announced August 2001.

Comments: 30pages, 6figures, Accepted for publication in ApJ

Journal ref: Astrophys.J. 588 (2003) 7-17

Theoretical Models of Multi-waveband QSO Luminosity Functions

Authors: T. Hosokawa, S. Mineshige, T. Kawaguchi, K. Yoshikawa, M. Umemura

Abstract: Cosmological evolution of the QSO luminosity functions (LFs) at NIR/optical/X-ray bands for 1.3 < z < 3.5 is investigated based on the realistic QSO spectra. The accretion-disk theory predicts that although QSO luminosities only depend on mass-accretion rate, \Mdot, QSO spectra have a dependence on black-hole mass, M_{BH}, as well. The smaller M_{BH} is and/or the larger \Mdot is, the harder bec… ▽ More Cosmological evolution of the QSO luminosity functions (LFs) at NIR/optical/X-ray bands for 1.3 < z < 3.5 is investigated based on the realistic QSO spectra. The accretion-disk theory predicts that although QSO luminosities only depend on mass-accretion rate, \Mdot, QSO spectra have a dependence on black-hole mass, M_{BH}, as well. The smaller M_{BH} is and/or the larger \Mdot is, the harder becomes the QSO NIR/optical/UV spectrum. We model disk spectra which can reproduce these features and calculated LFs for redshift z ~ 3 with the assumption of new-born QSOs being shining at the Eddington luminosity. The main results are: (i) the observed LFs at optical and X-rays can be simultaneously reproduced. (ii) LFs at optical and X-ray bands are not sensitive to M_{BH}, while LFs at NIR bands are; about one order of magnitude difference is expected in volume number densities at L_{I, J} ~ 10^{46} erg s^{-1} between the case that all QSOs would have the same spectral shape as that of M_{BH} = 10^{9} M_sun and the case with M_{BH} = 10^{11} M_sun. (iii) The resultant LFs at NIR are dominated by 10^{7} M_sun black-holes at L_{I, J} ~ 10^{44} erg s^{-1}, and by 10^{11} M_sun black-holes at L_{I, J} \~ 10^{46} erg s^{-1}. Future infrared observations from space(e.g.NGST) will probe cosmological evolution of black hole masses. For redshift z < 3, on the other hand, the observed optical/X-ray LFs can be fitted, if the initial QSO luminosity L_0 is below the Eddington luminosity. Interestingly, the best fitting values of l = L_0/L_{Edd} are different in B- and X-ray bands; l_B ~ 2.5 l_X. The reason for this discrepancy is briefly discussed. △ Less

Submitted 17 July, 2001; v1 submitted 14 July, 2001; originally announced July 2001.

Comments: 10 pages,7 Figures,to be published in Publications of the Astronomical Society of Japan

Nonlinear Stochastic Biasing of Galaxies and Dark Halos in Cosmological Hydrodynamic Simulations

Authors: Kohji Yoshikawa, Atsushi Taruya, Y. P. Jing, Yasushi Suto

Abstract: We perform an extensive analysis of nonlinear and stochastic biasing of galaxies and dark halos in spatially flat low-density CDM universe using cosmological hydrodynamic simulations. We compare their biasing properties with the predictions of an analytic halo biasing model. Dark halos in our simulations exhibit reasonable agreement with the predictions only on scales larger than 10h^{-1}Mpc, an… ▽ More We perform an extensive analysis of nonlinear and stochastic biasing of galaxies and dark halos in spatially flat low-density CDM universe using cosmological hydrodynamic simulations. We compare their biasing properties with the predictions of an analytic halo biasing model. Dark halos in our simulations exhibit reasonable agreement with the predictions only on scales larger than 10h^{-1}Mpc, and on smaller scales the volume exclusion effect of halos due to their finite size becomes substantial. Interestingly the biasing properties of galaxies are better described by extrapolating the halo biasing model predictions. We also find the clear dependence of galaxy biasing on their formation epoch; the distribution of old populations of galaxies tightly correlates with the underlying mass density field, while that of young populations is slightly more stochastic and anti-biased relative to dark matter. The amplitude of two-point correlation function of old populations is about 3 times larger than that of the young populations. Furthermore, the old population of galaxies reside within massive dark halos while the young galaxies are preferentially formed in smaller dark halos. Assuming that the observed early and late-type galaxies correspond to the simulated old and young populations of galaxies, respectively, all of these segregations of galaxies are consistent with observational ones for the early and late-type of galaxies such as the morphology--density relation of galaxies. △ Less

Submitted 23 April, 2001; originally announced April 2001.

Comments: 28 pages, 14 figures, accepted for publication in ApJ, Abstract abridged. For preprint with higher-resolution PS files, see ftp://www.kusastro.kyoto-u.ac.jp/pub/kohji/ytjs2001/

Journal ref: Astrophys.J.558:520-534,2001

Impact of Future Submillimeter and Millimeter Large Facilities on the Studies of Galaxy Formation and Evolution

Authors: Tsutomu T. Takeuchi, Ryohei Kawabe, Kotaro Kohno, Koichiro Nakanishi, Takako T. Ishii, Hiroyuki Hirashita, Kohji Yoshikawa

Abstract: We investigate what we can learn about galaxy formation and evolution from the data which will be obtained by the forthcoming large submillimeter/millimeter facilities, mainly by the Atacama Submillimeter Telescope Experiment (ASTE) and the Atacama Large Millimeter Array/Large Millimeter and Submillimeter Array (ALMA/LMSA). We first calculate the source counts from 350 um to 3 mm using the empir… ▽ More We investigate what we can learn about galaxy formation and evolution from the data which will be obtained by the forthcoming large submillimeter/millimeter facilities, mainly by the Atacama Submillimeter Telescope Experiment (ASTE) and the Atacama Large Millimeter Array/Large Millimeter and Submillimeter Array (ALMA/LMSA). We first calculate the source counts from 350 um to 3 mm using the empirical infrared galaxy number count model of Takeuchi et al. (2001). Based on the number counts, we evaluate the source confusion and determine the confusion limit at various wavebands as a function of the characteristic beam size. At submillimeter wavelengths, source confusion with the 10-15-m class facilities becomes severe at 0.1 to 1 mJy level, and astrometry and flux measurements will be difficult. However, we show that very a large-area survey of submillimeter sources brighter than 10-50 mJy can provide a unique constraint on infrared galaxy evolution at z = 1-2. We also find that the 5σ-confusion limit of the LMSA is fainter than 1 μJy. The source counts at such faint flux levels give important information of the epoch of galaxy formation. We then show that multiband photometry from the infrared (by ASTRO-F) to the millimeter can be utilized as a redshift estimator. In addition, we compare the observed 1.4, 5, and 8-GHz source counts with our model counts to examine the contribution of star forming galaxies to the faint radio galaxies. △ Less

Submitted 31 January, 2001; originally announced January 2001.

Comments: 17 pages LaTeX, PASP accepted. Abstract abridged

Exploring Galaxy Evolution from Infrared Number Counts and Cosmic Infrared Background

Authors: Tsutomu T. Takeuchi, Takako T. Ishii, Hiroyuki Hirashita, Kohji Yoshikawa, Hideo Matsuhara, Kimiaki Kawara, Haruyuki Okuda

Abstract: Recently reported infrared (IR) galaxy number counts and cosmic infrared background (CIRB) all suggest that galaxies have experienced a strong evolution sometime in their lifetime. We statistically estimate the galaxy evolution history from these data. We find that an order of magnitude increase of the far-infrared (FIR) luminosity at redshift z = 0.5 - 1.0 is necessary to reproduce the very hig… ▽ More Recently reported infrared (IR) galaxy number counts and cosmic infrared background (CIRB) all suggest that galaxies have experienced a strong evolution sometime in their lifetime. We statistically estimate the galaxy evolution history from these data. We find that an order of magnitude increase of the far-infrared (FIR) luminosity at redshift z = 0.5 - 1.0 is necessary to reproduce the very high CIRB intensity at 140 um reported by Hauser et al. (1998). z \sim 0.75 and decreases to, even at most, a factor of 10 toward z \sim 5, though many variants are allowed within these constraints. This evolution history also satisfies the constraints from the galaxy number counts obtained by IRAS, ISO and, roughly, SCUBA. The rapid evolution of the comoving IR luminosity density required from the CIRB well reproduces the very steep slope of galaxy number counts obtained by ISO. We also estimate the cosmic star formation history (SFH) from the obtained FIR luminosity density, considering the effect of the metal enrichment in galaxies. The derived SFH increases steeply with redshift in 0 < z < 0.75, and becomes flat or even declines at z > 0.75. This is consistent with the SFH estimated from the reported ultraviolet luminosity density. In addition, we present the performance of the Japanese ASTRO-F FIR galaxy survey. We show the expected number counts in the survey. We also evaluate how large a sky area is necessary to derive a secure information of galaxy evolution up to z \sim 1 from the survey, and find that at least 50 - 300 deg^2 is required. △ Less

Submitted 28 September, 2000; originally announced September 2000.

Comments: 18 pages LaTeX, PASJ in press. Abstract abridged