Showing 1–50 of 135 results for author: Fujita, Y

Indications of an offset merger in Abell 3667

Authors: Y. Omiya, K. Nakazawa, T. Tamura, H. Akamatsu, K. Matsushita, N. Okabe, K. Sato, Y. Fujita, L. Gu, A. Simionescu, Y. Ichinohe, C. J. Riseley, T. Akahori, D. Ito, K. Sakai, K. Kurahara

Abstract: Abell 3667 is a nearby merging cluster with a prominent cold front and a pair of two bright radio relics. Assuming a head-on merger, the origin of the cold front is often considered to be a remnant of the cluster core stripped by its surrounding ICM. Some authors have proposed an offset merger scenario in which the subcluster core rotates after the first core crossing. This scenario can reproduce… ▽ More Abell 3667 is a nearby merging cluster with a prominent cold front and a pair of two bright radio relics. Assuming a head-on merger, the origin of the cold front is often considered to be a remnant of the cluster core stripped by its surrounding ICM. Some authors have proposed an offset merger scenario in which the subcluster core rotates after the first core crossing. This scenario can reproduce features such as the cold front and a pair of radio relics. To distinguish between these scenarios, we reanalyzed the ICM distribution and measured the line-of-sight bulk ICM velocity using the XMM-Newton PN data. In the unsharp masked image, we identify several ICM features. The notable feature is a RG1 vortex, which is a clockwise vortex-like enhancement with a radius of about 250 kpc connecting the first BCG to the radio galaxy (RG1). It is particularly enhanced near the north of the 1st BCG, which is named the BCG-N tail. The thermodynamic maps show that the ICM of the RG1 vortex has a relatively high abundance of 0.5-0.6 solar compared to the surrounding regions. The ICM of the BCG-E tail also has a high abundance and low pseudo-entropy and can be interpreted as a remnant of the cluster core's ICM. Including its arc-like shape, the RG1 vortex supports the idea that the ICM around the cluster center is rotating, which is natural for an offset merger scenario. The results of the line-of-sight bulk ICM velocity measurements show that the ICM around the BCG-N tail is redshifted with a velocity difference of 940+/-440 km/s compared to the optical redshift of the first BCG. We obtain other indications of variations in the line-of-sight velocity of the ICM and discuss these in the context of an offset merger. △ Less

Submitted 26 June, 2024; v1 submitted 15 March, 2024; originally announced March 2024.

Comprehensive study of magnetic field evolution in relativistic jets based on 2D simulations

Authors: Amin Esmaeili, Yutaka Fujita

Abstract: We use two-dimensional particle-in-cell simulations to investigate the generation and evolution of the magnetic field associated with the propagation of a jet for various initial conditions. We demonstrate that, in general, the magnetic field is initially grown by the Weibel and Mushroom instabilities. However, the field is saturated by the Alfv'en current limit. For initially non-magnetized plasm… ▽ More We use two-dimensional particle-in-cell simulations to investigate the generation and evolution of the magnetic field associated with the propagation of a jet for various initial conditions. We demonstrate that, in general, the magnetic field is initially grown by the Weibel and Mushroom instabilities. However, the field is saturated by the Alfv'en current limit. For initially non-magnetized plasma, we show that the growth of the magnetic field is delayed when the matter density of the jet environment is lower, which are in agreement with simple analytical predictions. We show that the higher Lorentz factor ($\gtrsim 2$) prevents rapid growth of the magnetic fields. When the initial field is troidal, the position of the magnetic filaments moves away from the jet as the field strength increases. The axial initial field helps the jet maintain its shape more effectively than the troidal initial field. △ Less

Submitted 5 February, 2024; originally announced February 2024.

Report number: Volume 31, Issue 2, 022305

Journal ref: Physics of Plasmas 2024

Broadband non-thermal emission of odd radio circles induced by explosive galactic outflow remnants and their evolution

Authors: Yutaka Fujita, Norita Kawanaka, Susumu Inoue

Abstract: Odd radio circles (ORCs) are mysterious rings of faint, diffuse emission recently discovered in radio surveys, some of which may be associated with galaxies in relatively dense environments. We propose such ORCs to be synchrotron emission from remnants of explosive galactic outflows, calling them OGREs, and discuss their broadband non-thermal emission and evolution. We posit that a large amount of… ▽ More Odd radio circles (ORCs) are mysterious rings of faint, diffuse emission recently discovered in radio surveys, some of which may be associated with galaxies in relatively dense environments. We propose such ORCs to be synchrotron emission from remnants of explosive galactic outflows, calling them OGREs, and discuss their broadband non-thermal emission and evolution. We posit that a large amount of energy was ejected from the central galaxy in the past, creating an outgoing shock that accelerates cosmic rays. Assuming plausible values for the density, temperature and magnetic field of the ambient medium, consistency with the observed spectral index, size and power of the ORCs requires the energy to be as high as ~10^60 erg, suggesting that their sources could be active galactic nuclei. We calculate the spectral energy distributions (SEDs) of the OGREs and their evolution, including synchrotron, inverse Compton (IC) and bremsstrahlung emission from electrons, and pion-decay emission from protons. We find that the SEDs of the younger OGREs are not greatly different from those of older ones currently observable as ORCs if radiative cooling of electrons is effective. As such younger OGREs are expected to be rarer and smaller, they may not be readily observable. However, if radiative cooling of electrons is ineffective, younger OGREs may be detectable in X-rays. △ Less

Submitted 1 May, 2024; v1 submitted 20 November, 2023; originally announced November 2023.

Comments: Accepted for publication in PASJ

The relationships between AGN power and molecular gas mass within 500 pc of the center of elliptical galaxies

Authors: Yutaka Fujita, Takuma Izumi, Hiroshi Nagai, Nozomu Kawakatu, Norita Kawanaka

Abstract: The physical quantity that directly controls the feedback of active galactic nuclei (AGNs) in elliptical galaxies remains to be determined. The discovery of molecular gas around the AGNs suggests that the gas is fueling the AGNs. Therefore, we analyze Atacama Large Millimeter/submillimeter Array (ALMA) data for the CO line (J=1-0, 2-1, 3-2) emission and estimate the mass of molecular gas within 50… ▽ More The physical quantity that directly controls the feedback of active galactic nuclei (AGNs) in elliptical galaxies remains to be determined. The discovery of molecular gas around the AGNs suggests that the gas is fueling the AGNs. Therefore, we analyze Atacama Large Millimeter/submillimeter Array (ALMA) data for the CO line (J=1-0, 2-1, 3-2) emission and estimate the mass of molecular gas within 500 pc of the center of 12 non-central elliptical galaxies (NCEGs) and 10 of the brightest cluster galaxies (BCGs). We find that the mass (M_mol ~ 10^5-10^9 M_sun) is correlated with the jet power of their AGNs, which is represented by P_cav ~ 4.1x10^42 (M_mol/10^7 M_sun)^{1.3} erg s^{-1}, although NCEGs alone do not show the correlation. We also find that M_mol is correlated with the AGN continuum luminosities at ~ 1.4 GHz (L_1.4) and ~ 100-300 GHz (L_con). Since P_cav reflects galactic-scale, long-term AGN activity, while the continuum luminosities reflect local (~< 500 pc), short-term AGN activity, our results suggest that AGN activity depends on the amount of gas, regardless of its time scale. On the other hand, we cannot find a clear correlation between the mass of the black holes in the AGNs (M_BH) and P_cav. This suggests that M_mol, rather than M_BH, is the main factor controlling AGN activity. We confirm that the origin of the continuum emission from the AGNs at ~ 1.4-300 GHz is mostly synchrotron radiation. △ Less

Submitted 11 February, 2024; v1 submitted 5 October, 2023; originally announced October 2023.

Comments: Accepted for publication in ApJ

Supermassive black hole feeding and feedback observed on sub-parsec scales

Authors: Takuma Izumi, Keiichi Wada, Masatoshi Imanishi, Kouichiro Nakanishi, Kotaro Kohno, Yuki Kudoh, Taiki Kawamuro, Shunsuke Baba, Naoki Matsumoto, Yutaka Fujita, Konrad R. W. Tristram

Abstract: Active galaxies contain a supermassive black hole at their center, which grows by accreting matter from the surrounding galaxy. The accretion process in the central ~10 parsecs has not been directly resolved in previous observations, due to the small apparent angular sizes involved. We observed the active nucleus of the Circinus Galaxy using sub-millimeter interferometry. A dense inflow of molecul… ▽ More Active galaxies contain a supermassive black hole at their center, which grows by accreting matter from the surrounding galaxy. The accretion process in the central ~10 parsecs has not been directly resolved in previous observations, due to the small apparent angular sizes involved. We observed the active nucleus of the Circinus Galaxy using sub-millimeter interferometry. A dense inflow of molecular gas is evident on sub-parsec scales. We calculate that less than 3% of this inflow is accreted by the black hole, with the rest being ejected by multiphase outflows, providing feedback to the host galaxy. The observations also reveal a dense gas disk surrounding the inflow; the disk is gravitationally unstable which drives the accretion into the central ~1 parsec. △ Less

Submitted 13 November, 2023; v1 submitted 6 May, 2023; originally announced May 2023.

Comments: First release on Nov 3, 2023 in Science. 32 pages in one column = Main (13 pages, 4 figures) + Supplement (19 pages, 9 figures + 2 tables). This is the accepted version after peer review

Journal ref: Science, Volume 382, Year 2023, Issue 6670, Pages 554-559

The correlation between the 500 pc scale molecular gas masses and AGN powers for massive elliptical galaxies

Authors: Yutaka Fujita, Takuma Izumi, Nozomu Kawakatu, Hiroshi Nagai, Ryo Hirasawa, Yu Ikeda

Abstract: Massive molecular clouds have been discovered in massive elliptical galaxies at the center of galaxy clusters. Some of this cold gas is expected to flow in the central supermassive black holes and activate galactic nucleus (AGN) feedback. In this study, we analyze archival ALMA data of 9 massive elliptical galaxies, focusing on CO line emissions, to explore the circumnuclear gas. We show that the… ▽ More Massive molecular clouds have been discovered in massive elliptical galaxies at the center of galaxy clusters. Some of this cold gas is expected to flow in the central supermassive black holes and activate galactic nucleus (AGN) feedback. In this study, we analyze archival ALMA data of 9 massive elliptical galaxies, focusing on CO line emissions, to explore the circumnuclear gas. We show that the mass of the molecular gas within a fixed radius (500 pc) from the AGNs (M_mol ~ 10^7-10^8 M_sun) is correlated with the jet power estimated from X-ray cavities (P_cav ~ 10^42-10^45 erg/s). The mass accretion rate of the circumnuclear gas \dot{M} also has a correlation with P_cav. On the other hand, the continuum luminosities at ~1.4 GHz and ~100-300 GHz have no correlation with M_mol. These results indicate that the circumnuclear gas is sustaining the long-term AGN activities (~10^7 yr) rather than the current ones. The circumnuclear gas mass is a better indicator of the jet power than the continuum luminosity, which probably changes on a shorter time scale. We also study the origin of the continuum emission from the AGNs at ~100-300 GHz and find that it is mostly synchrotron radiation. For low-luminosity AGNs, however, dust emission appears to contaminate the continuum. △ Less

Submitted 21 July, 2023; v1 submitted 29 March, 2023; originally announced March 2023.

Comments: Accepted for publication in PASJ

A broadband X-ray imaging spectroscopy in the 2030s: the FORCE mission

Authors: Koji Mori, Takeshi G. Tsuru, Kazuhiro Nakazawa, Yoshihiro Ueda, Shin Watanabe, Takaaki Tanaka, Manabu Ishida, Hironori Matsumoto, Hisamitsu Awaki, Hiroshi Murakami, Masayoshi Nobukawa, Ayaki Takeda, Yasushi Fukazawa, Hiroshi Tsunemi, Tadayuki Takahashi, Ann Hornschemeier, Takashi Okajima, William W. Zhang, Brian J. Williams, Tonia Venters, Kristin Madsen, Mihoko Yukita, Hiroki Akamatsu, Aya Bamba, Teruaki Enoto , et al. (27 additional authors not shown)

Abstract: In this multi-messenger astronomy era, all the observational probes are improving their sensitivities and overall performance. The Focusing on Relativistic universe and Cosmic Evolution (FORCE) mission, the product of a JAXA/NASA collaboration, will reach a 10 times higher sensitivity in the hard X-ray band ($E >$ 10~keV) in comparison with any previous hard X-ray missions, and provide simultaneou… ▽ More In this multi-messenger astronomy era, all the observational probes are improving their sensitivities and overall performance. The Focusing on Relativistic universe and Cosmic Evolution (FORCE) mission, the product of a JAXA/NASA collaboration, will reach a 10 times higher sensitivity in the hard X-ray band ($E >$ 10~keV) in comparison with any previous hard X-ray missions, and provide simultaneous soft X-ray coverage. FORCE aims to be launched in the early 2030s, providing a perfect hard X-ray complement to the ESA flagship mission Athena. FORCE will be the most powerful X-ray probe for discovering obscured/hidden black holes and studying high energy particle acceleration in our Universe and will address how relativistic processes in the universe are realized and how these affect cosmic evolution. FORCE, which will operate over 1--79 keV, is equipped with two identical pairs of supermirrors and wideband X-ray imagers. The mirror and imager are connected by a high mechanical stiffness extensible optical bench with alignment monitor systems with a focal length of 12~m. A light-weight silicon mirror with multi-layer coating realizes a high angular resolution of $<15''$ in half-power diameter in the broad bandpass. The imager is a hybrid of a brand-new SOI-CMOS silicon-pixel detector and a CdTe detector responsible for the softer and harder energy bands, respectively. FORCE will play an essential role in the multi-messenger astronomy in the 2030s with its broadband X-ray sensitivity. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: 12 pages, 8 figures. Proceedings of SPIE Astronomical Telescopes and Instrumentation 2022

Journal ref: Proc. of SPIE vol. 12181, 1218122 (2022)

XMM-Newton view of the shock heating in an early merging cluster, CIZA J1358.9$-$4750

Authors: Y. Omiya, K. Nakazawa, K. Matsushita, S. B. Kobayashi, N. Okabe, K. Sato, T. Tamura, Y. Fujita, L. Gu, T. Kitayama, T. Akahori, K. Kurahara, T. Yamaguchi

Abstract: CIZA J1358.9-4750 is a nearby galaxy cluster in the early phase of a major merger. The two-dimensional temperature map using XMM-Newton EPIC-PN observation confirms the existence of a high temperature region, which we call the "hot region", in the "bridge region" connecting the two clusters. The ~ 500 kpc wide region between the southeast and northwest boundaries also has higher pseudo pressure co… ▽ More CIZA J1358.9-4750 is a nearby galaxy cluster in the early phase of a major merger. The two-dimensional temperature map using XMM-Newton EPIC-PN observation confirms the existence of a high temperature region, which we call the "hot region", in the "bridge region" connecting the two clusters. The ~ 500 kpc wide region between the southeast and northwest boundaries also has higher pseudo pressure compared to the unshocked regions, suggesting the existence of two shocks. The southern shock front is clearly visible in the X-ray surface brightness image and has already been reported by Kato et al. (2015). The northern one, on the other hand, is newly discovered. To evaluate their Mach number, we constructed a three-dimensional toy merger model with overlapping shocked and unshocked components in line of sight. The unshocked and preshock ICM conditions are estimated based on those outside the interacting bridge region assuming point symmetry. The hot region spectra are modeled with two-temperature thermal components, assuming that the shocked condition follows the Rankin-Hugoniot relation with the preshock condition. As a result, the shocked region is estimated to have a line-of-sight depth of ~ 1 Mpc with a Mach number of ~ 1.3 in the southeast shock and ~ 1.7 in the northwest shock. The age of the shock waves is estimated to be ~ 260 Myr. This three dimensional merger model is consistent with the Sunyaev-Zeldovich signal obtained using the Planck observation within the CMB fluctuations. The total flow of the kinetic energy of the ICM through the southeast shock was estimated to be ~ 2.2 x $10^{42}$ erg/s. Assuming that 10 % of this energy is converted into ICM turbulence, the line-of-sight velocity dispersion is calculated to be ~ 200 km/s, which is basically resolvable via coming high spectral resolution observations. △ Less

Submitted 23 October, 2022; v1 submitted 5 October, 2022; originally announced October 2022.

Diffuse radio source candidate in CIZA J1358.9-4750

Authors: Kohei Kurahara, Takuya Akahori, Ruta Kale, Hiroki Akamatsu, Yutaka Fujita, Liyi Gu, Huib Intema, Kazuhiro Nakazawa, Nobuhiro Okabe, Yuki Omiya, Viral Parekh, Timothy Shimwell, Motokazu Takizawa, Reinout van Weeren

Abstract: We report on results of our upgraded Giant Metrewave Radio Telescope (uGMRT) observations for an early-stage merging galaxy cluster, CIZA J1358.9-4750 (CIZA1359), in Band-3 (300--500 MHz). We achieved the image dynamic range of $\sim 38,000$ using the direction dependent calibration and found a candidate of diffuse radio emission at 4~$σ_{rms}$ significance. The flux density of the candidate at 40… ▽ More We report on results of our upgraded Giant Metrewave Radio Telescope (uGMRT) observations for an early-stage merging galaxy cluster, CIZA J1358.9-4750 (CIZA1359), in Band-3 (300--500 MHz). We achieved the image dynamic range of $\sim 38,000$ using the direction dependent calibration and found a candidate of diffuse radio emission at 4~$σ_{rms}$ significance. The flux density of the candidate at 400~MHz, $24.04 \pm 2.48$~mJy, is significantly positive compared to noise, where its radio power, $2.40 \times 10^{24}$~W~Hz$^{-1}$, is consistent with those of typical diffuse radio sources of galaxy clusters. The candidate is associated with a part of the X-ray shock front at which the Mach number reaches its maximum value of $\mathcal{M}\sim 1.7$. The spectral index ($F_ν\propto ν^α$) of the candidate, $α= - 1.22 \pm 0.33$, is in agreement with an expected value derived from the standard diffusive shock acceleration (DSA) model. But such a low Mach number with a short acceleration time would require seed cosmic-rays supplied from active galactic nucleus (AGN) activities of member galaxies, as suggested in some other clusters. Indeed, we found seven AGN candidates inside the diffuse source candidate. Assuming the energy equipartition between magnetic fields and cosmic-rays, the magnetic field strength of the candidate was estimated to be $2.1~μ$G. We also find head-tail galaxies and radio phoenixes or fossils near the CIZA1359. △ Less

Submitted 23 November, 2022; v1 submitted 9 August, 2022; originally announced August 2022.

Comments: 17 pages, 9 figures, Accepted to PASJ

Evidence for powerful winds and the associated reverse shock as the origin of the Fermi bubbles

Authors: Yutaka Fujita

Abstract: The Fermi bubbles are large gamma-ray-emitting structures. They are symmetric about the Galactic Centre (GC), and their creation is therefore attributed to intensive energy injection at the GC. In this study, we focus on the non-equilibrium X-ray gas structures associated with the bubbles. We show that a combination of the density, temperature, and shock age profiles of the X-ray gas can be used t… ▽ More The Fermi bubbles are large gamma-ray-emitting structures. They are symmetric about the Galactic Centre (GC), and their creation is therefore attributed to intensive energy injection at the GC. In this study, we focus on the non-equilibrium X-ray gas structures associated with the bubbles. We show that a combination of the density, temperature, and shock age profiles of the X-ray gas can be used to distinguish the energy injection mechanisms. By comparing the results of numerical simulations with observations, we indicate that the bubbles were created by a fast wind from the GC because it generates a strong reverse shock and reproduces the observed temperature peak there. On the other hand, instantaneous energy injection at the GC cannot reproduce the temperature profile. The wind had a speed of ~1000 km/s, and blew for ~10^7 yr. Because the mass flux of the wind is large, the entrainment of interstellar gas by wide-angle outflows from the black hole is required. Thus, the wind may be the same as active galactic nuclei outflows often observed in other galaxies and thought to regulate the growth of galaxies and their central black holes. △ Less

Submitted 14 December, 2022; v1 submitted 2 August, 2022; originally announced August 2022.

Journal ref: 2023, MNRAS, 518, 4551

The Gamma-Ray Emission from the Supernova Remnant RX J1713.7-3946 Interacting with Two-phase Medium

Authors: Yutaka Fujita, Ryo Yamazaki, Yutaka Ohira

Abstract: We study the origin of gamma rays from the supernova remnant (SNR) RX J1713.7-3946. Using an analytical model, we calculate the distribution of cosmic rays (CRs) around the SNRs. Motivated by the results of previous studies, we assume that the SNR is interacting with two-phase interstellar medium (ISM), where dense clumps are surrounded by tenuous interclump medium. We also assume that only higher… ▽ More We study the origin of gamma rays from the supernova remnant (SNR) RX J1713.7-3946. Using an analytical model, we calculate the distribution of cosmic rays (CRs) around the SNRs. Motivated by the results of previous studies, we assume that the SNR is interacting with two-phase interstellar medium (ISM), where dense clumps are surrounded by tenuous interclump medium. We also assume that only higher-energy protons (~> TeV) can penetrate the dense clumps. We find that pi^0-decay gamma rays produced by protons reproduce the observed gamma-ray spectrum peaked at ~TeV. On the other hand, it has recently been indicated that the observed ISM column density (N_p), the X-ray surface brightness (I_X), and the gamma-ray surface brightness (I_g) at grid points across the SNR form a plane in the three-dimensional (3D) space of (N_p, I_X, I_g). We find that the planar configuration is naturally reproduced if the ISM or the CR electron-to-proton ratio is not spherically uniform. We show that the shift of the observed data in the 3D space could be used to identify which of the quantities, the ISM density, the CR electron-to-proton ratio, or the magnetic field, varies in the azimuthal direction of the SNR. △ Less

Submitted 31 July, 2022; v1 submitted 24 May, 2022; originally announced May 2022.

Comments: matched the published version

Suzaku Observations of Fe K-shell Lines in the Supernova Remnant W51C and Hard X-ray Sources in the Proximity

Authors: Aika Shimaguchi, Kumiko K. Nobukawa, Shigeo Yamauchi, Masayoshi Nobukawa, Yutaka Fujita

Abstract: In this paper, we investigated the Fe K-shell lines in the supernova remnant W51C and hard X-ray sources in the proximity. We measured the intensities of Fe I K$α$ and Fe XXV He$α$ lines at 6.40 keV and 6.68 keV, respectively, and found that the intensity of the 6.68 keV line is consistent with the background level expected from previous studies, while that of the 6.40 keV line is higher at the si… ▽ More In this paper, we investigated the Fe K-shell lines in the supernova remnant W51C and hard X-ray sources in the proximity. We measured the intensities of Fe I K$α$ and Fe XXV He$α$ lines at 6.40 keV and 6.68 keV, respectively, and found that the intensity of the 6.68 keV line is consistent with the background level expected from previous studies, while that of the 6.40 keV line is higher at the significance level of 2.0$σ$. Given the presence of gamma-ray emission and high ionization rate point spatially coincident with the remnant, we conclude that the enhanced 6.40 keV line most likely originates from the interaction between low-energy cosmic rays and molecular clouds. Also, we discovered an enhanced 6.68 keV line emission from the compact H II region G49.0-0.3 at the significance level of 3.4$σ$. Spectral analysis revealed that the temperature and abundance of the thermal plasma with the 6.68 keV line is $kT = 3.0^{+0.8}_{-0.7}$ keV and $Z = 0.5 \pm 0.2$ solar, respectively. These values are explained by the thermal plasma generated by the stellar winds of O stars. △ Less

Submitted 7 March, 2022; originally announced March 2022.

Comments: 8 pages, 4 figures, accepted for publication in Publications of the Astronomical Society of Japan (PASJ)

Massive Molecular Gas as a Fuel Tank for Active Galactic Nuclei Feedback in Central Cluster Galaxies

Authors: Yutaka Fujita, Nozomu Kawakatu, Hiroshi Nagai

Abstract: Massive molecular gas has been discovered in giant elliptical galaxies at the centers of galaxy clusters. To reveal its role in active galactic nucleus (AGN) feedback in those galaxies, we construct a semianalytical model of gas circulation. This model especially focuses on the massive molecular gas (interstellar cold gas on a scale of ~ 10 kpc) and the circumnuclear disk (~< 0.5 kpc). We consider… ▽ More Massive molecular gas has been discovered in giant elliptical galaxies at the centers of galaxy clusters. To reveal its role in active galactic nucleus (AGN) feedback in those galaxies, we construct a semianalytical model of gas circulation. This model especially focuses on the massive molecular gas (interstellar cold gas on a scale of ~ 10 kpc) and the circumnuclear disk (~< 0.5 kpc). We consider the destruction of the interstellar cold gas by star formation and the gravitational instability for the circumnuclear disk. Our model can reproduce the basic properties of the interstellar cold gas and the circumnuclear disk, such as their masses. We also find that the circumnuclear disk tends to stay at the boundary between stable and unstable states. This works as an 'adjusting valve' that regulates mass accretion toward the supermassive black hole. On the other hand, the interstellar cold gas serves as a 'fuel tank' in the AGN feedback. Even if the cooling of the galactic hot gas is prevented, the interstellar cold gas can sustain the AGN activity for ~> 0.5 Gyr. We also confirm that the small entropy of the hot gas (~< 30 keV cm^2) or the short cooling time (~< 1 Gyr) is a critical condition for the existence of the massive amounts of molecular gas in the galaxy. The dissipation time of the interstellar cold gas may be related to the critical cooling time. The galaxy behavior is described by a simple relation among the disk stability, the cloud dissipation time, and the gas cooling rate. △ Less

Submitted 6 January, 2022; v1 submitted 3 August, 2021; originally announced August 2021.

Comments: matched the published version

Journal ref: The Astrophysical Journal, 924:24 (14pp), 2022 January 1

X-Ray Emission from the PeVatron-candidate Supernova Remnant G106.3+2.7

Authors: Yutaka Fujita, Aya Bamba, Kumiko K. Nobukawa, Hironori Matsumoto

Abstract: We report a discovery of diffuse X-ray emission around the supernova remnant (SNR) G106.3+2.7, which is associated with VER J2227+608 and HAWC J2227+610 and is known as a candidate for a PeV cosmic ray accelerator (PeVatron). We analyze observational data of Suzaku around the SNR and the adjacent pulsar PSR J2229+6114. We find diffuse X-ray emission that is represented by either thermal or non-the… ▽ More We report a discovery of diffuse X-ray emission around the supernova remnant (SNR) G106.3+2.7, which is associated with VER J2227+608 and HAWC J2227+610 and is known as a candidate for a PeV cosmic ray accelerator (PeVatron). We analyze observational data of Suzaku around the SNR and the adjacent pulsar PSR J2229+6114. We find diffuse X-ray emission that is represented by either thermal or non-thermal one. However, the metal abundance for the thermal emission is <0.13 Z_sun, which may be too small in the Milky Way and suggests that the emission is non-thermal. The intensity of the diffuse emission increases toward PSR J2229+6114 in the same way as radio emission, and it is in contrast with gamma-ray emission concentrated on a molecular cloud. The X-ray photon index does not change with the distance from the pulsar and it indicates that radiative cooling is ineffective and particle diffusion is not extremely slow. The X-ray and radio emissions seem to be of leptonic origin and the parent electrons may originate from the pulsar or its wind nebula. The gamma-ray emission appears to be of hadronic origin because of its spacial distribution. The parent protons may be tightly confined in the cloud separately from the diffusing electrons. △ Less

Submitted 23 May, 2021; v1 submitted 25 January, 2021; originally announced January 2021.

Comments: matched the published version

Journal ref: The Astrophysical Journal, 912:133 (6pp), 2021 May 13

Intrusion of MeV-TeV Cosmic-Rays into Molecular Clouds Studied by Ionization, the Neutral Iron Line, and Gamma-Rays

Authors: Yutaka Fujita, Kumiko K. Nobukawa, Hidetoshi Sano

Abstract: Low-energy (~MeV) cosmic rays (CRs) ionize molecular clouds and create the neutral iron line (Fe I Kα) at 6.4 keV. On the other hand, high-energy (>~ GeV) CRs interact with the dense cloud gas and produce gamma rays. Based on a one-dimensional model, we study the spatial correlation among ionization rates of gas, 6.4 keV line fluxes, and gamma-ray emissions from a molecular cloud illuminated by CR… ▽ More Low-energy (~MeV) cosmic rays (CRs) ionize molecular clouds and create the neutral iron line (Fe I Kα) at 6.4 keV. On the other hand, high-energy (>~ GeV) CRs interact with the dense cloud gas and produce gamma rays. Based on a one-dimensional model, we study the spatial correlation among ionization rates of gas, 6.4 keV line fluxes, and gamma-ray emissions from a molecular cloud illuminated by CRs accelerated at an adjacent supernova remnant. We find that the spatial distributions of these three observables depend on how CRs intrude the cloud and on the internal structure of the cloud. If the intrusion is represented by slow diffusion, the 6.4 keV line should be detected around the cloud edge where ionization rates are high. On the other hand, if CRs freely stream in the cloud, the 6.4 keV line should be observed where gamma rays are emitted. In the former, the cooling time of the CRs responsible for the 6.4 keV line is shorter than their cloud crossing time, and it is opposite in the latter. Although we compare the results with observations, we cannot conclude whether the diffusion or the free-streaming is dominantly realized. Our predictions can be checked in more detail with future X-ray missions such as XRISM and Athena and by observations of ionization rates that cover wider fields. △ Less

Submitted 18 February, 2021; v1 submitted 28 September, 2020; originally announced September 2020.

Journal ref: The Astrophysical Journal, 908:136 (7pp), 2021 February 20

ALMA Observations of PSR B1259-63/LS 2883 in an Inactive Period: Variable Circumstellar Disk?

Authors: Yutaka Fujita, Hiroshi Nagai. Takuya Akahori, Akiko Kawachi, Atsuo T. Okazaki

Abstract: We report Atacama Large Millimeter/submillimeter Array (ALMA) observations of the gamma-ray binary system containing the pulsar PSR B1259-63 orbiting around a massive star LS 2883 in an inactive period between the 2017 and 2021 periastron passages. We detected radio continuum emission from the binary system at 97 GHz (Band 3) and 343 GHz (Band 7). Compared with our previous ALMA observations perfo… ▽ More We report Atacama Large Millimeter/submillimeter Array (ALMA) observations of the gamma-ray binary system containing the pulsar PSR B1259-63 orbiting around a massive star LS 2883 in an inactive period between the 2017 and 2021 periastron passages. We detected radio continuum emission from the binary system at 97 GHz (Band 3) and 343 GHz (Band 7). Compared with our previous ALMA observations performed soon after the 2017 periastron passage, the fluxes have decreased by an factor of six at 97 GHz and two at 343 GHz. The flux at 343 GHz is large relative to that at 97 GHz and appears to be thermal emission from the circumstellar disk around LS 2883. The decrease of the 343 GHz flux may indicate that the disk has expanded and become partially optically-thin since the disk is no longer affected by pulsar winds. The flux at 97 GHz is consistent with that expected from the pulsed emission from the pulsar, which indicates that the unpulsed emission that had been produced through pulsar-disk or pulsar-stellar wind interaction has disappeared. The image of the system is consistent with a point source and shows no sign of ejecta. △ Less

Submitted 15 January, 2021; v1 submitted 30 April, 2020; originally announced May 2020.

Comments: Accepted for publication in PASJ

Journal ref: PASJ, 72, L9 (2020)

XXL Survey groups and clusters in the Hyper Suprime-Cam Survey. Scaling relations between X-ray properties and weak lensing mass

Authors: Mauro Sereno, Keiichi Umetsu, Stefano Ettori, Dominique Eckert, Fabio Gastaldello, Paul Giles, Maggie Lieu, Ben Maughan, Nobuhiro Okabe, Mark Birkinshaw, I-Non Chiu, Yutaka Fujita, Satoshi Miyazaki, David Rapetti, Elias Koulouridis, Marguerite Pierre

Abstract: Scaling relations trace the formation and evolution of galaxy clusters. We exploited multi-wavelength surveys -- the XXL survey at \emph{XMM-Newton} in the X-ray band, and the Hyper Suprime-Cam Subaru Strategic Program for optical weak lensing -- to study an X-ray selected, complete sample of clusters and groups. The scalings of gas mass, temperature, and soft-band X-ray luminosity with the weak l… ▽ More Scaling relations trace the formation and evolution of galaxy clusters. We exploited multi-wavelength surveys -- the XXL survey at \emph{XMM-Newton} in the X-ray band, and the Hyper Suprime-Cam Subaru Strategic Program for optical weak lensing -- to study an X-ray selected, complete sample of clusters and groups. The scalings of gas mass, temperature, and soft-band X-ray luminosity with the weak lensing mass show imprints of radiative cooling and AGN feedback in groups. From the multi-variate analysis, we found some evidence for steeper than self-similar slopes for gas mass ($β_{m_\text{g}|m}=1.73 \pm0.80$) and luminosity ($β_{l|m}=1.91\pm0.94$) and a nearly self-similar slope for the temperature ($β_{t|m}=0.78\pm0.43$). Intrinsic scatters of X-ray properties appear to be positively correlated at a fixed mass (median correlation factor $ρ_{X_1X_2|m}\sim0.34$) due to dynamical state and merger history of the halos. Positive correlations with the weak lensing mass (median correlation factor $ρ_{m_\text{wl}X|m}\sim0.35$) can be connected to triaxiality and orientation. Comparison of weak lensing and hydrostatic masses suggests a small role played by non-thermal pressure support ($9\pm17\%$). △ Less

Submitted 5 December, 2019; originally announced December 2019.

Comments: 19 pages; in press on MNRAS. One of two companion papers on HSC-XXL surveys (Keiichi Umetsu et al., arXiv:1909.10524)

Non-steady heating of cool cores of galaxy clusters by ubiquitous turbulence and AGN

Authors: Yutaka Fujita, Renyue Cen, Irina Zhuravleva

Abstract: Recent cosmological simulations have shown that turbulence should be generally prevailing in clusters because clusters are continuously growing through matter accretion. Using one-dimensional hydrodynamic simulations, we study the heating of cool-core clusters by the ubiquitous turbulence as well as feedback from the central active galactic nuclei (AGNs) for a wide range of cluster and turbulence… ▽ More Recent cosmological simulations have shown that turbulence should be generally prevailing in clusters because clusters are continuously growing through matter accretion. Using one-dimensional hydrodynamic simulations, we study the heating of cool-core clusters by the ubiquitous turbulence as well as feedback from the central active galactic nuclei (AGNs) for a wide range of cluster and turbulence parameters, focusing on the global stability of the core. We find that the AGN shows intermittent activities in the presence of moderate turbulence similar to the one observed with Hitomi. The cluster core maintains a quasi-equilibrium state for most of the time because the heating through turbulent diffusion is nearly balanced with radiative cooling. The balance is gradually lost because of slight dominance of the radiative cooling, and the AGN is ignited by increased gas inflow. Finally, when the AGN bursts, the core is heated almost instantaneously. Thanks to the pre-existing turbulence, the heated gas is distributed throughout the core without becoming globally unstable and causing catastrophic cooling, and the core recovers the quasi-equilibrium state. The AGN bursts can be stronger in lower-mass clusters. Predictions of our model can be easily checked with future X-ray missions like XRISM and Athena. △ Less

Submitted 20 May, 2020; v1 submitted 2 December, 2019; originally announced December 2019.

Journal ref: MNRAS, 494, 5507-5519 (2020)

Weak lensing Analysis of X-Ray-selected XXL Galaxy Groups and Clusters with Subaru HSC Data

Authors: Keiichi Umetsu, Mauro Sereno, Maggie Lieu, Hironao Miyatake, Elinor Medezinski, Atsushi J. Nishizawa, Paul Giles, Fabio Gastaldello, Ian G. McCarthy, Martin Kilbinger, Mark Birkinshaw, Stefano Ettori, Nobuhiro Okabe, I-Non Chiu, Jean Coupon, Dominique Eckert, Yutaka Fujita, Yuichi Higuchi, Elias Koulouridis, Ben Maughan, Satoshi Miyazaki, Masamune Oguri, Florian Pacaud, Marguerite Pierre, David Rapetti , et al. (1 additional authors not shown)

Abstract: We present a weak-lensing analysis of X-ray galaxy groups and clusters selected from the XMM-XXL survey using the first-year data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. Our joint weak-lensing and X-ray analysis focuses on 136 spectroscopically confirmed X-ray-selected systems at 0.031 < z < 1.033 detected in the 25sqdeg XXL-N region. We characterize the mass distributions of in… ▽ More We present a weak-lensing analysis of X-ray galaxy groups and clusters selected from the XMM-XXL survey using the first-year data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. Our joint weak-lensing and X-ray analysis focuses on 136 spectroscopically confirmed X-ray-selected systems at 0.031 < z < 1.033 detected in the 25sqdeg XXL-N region. We characterize the mass distributions of individual clusters and establish the concentration-mass (c-M) relation for the XXL sample, by accounting for selection bias and statistical effects, and marginalizing over the remaining mass calibration uncertainty. We find the mass-trend parameter of the c-M relation to be β= -0.07 \pm 0.28 and the normalization to be c200 = 4.8 \pm 1.0 (stat) \pm 0.8 (syst) at M200=10^{14}Msun/h and z = 0.3. We find no statistical evidence for redshift evolution. Our weak-lensing results are in excellent agreement with dark-matter-only c-M relations calibrated for recent LCDM cosmologies. The level of intrinsic scatter in c200 is constrained as σ(\ln[c200]) < 24% (99.7% CL), which is smaller than predicted for the full population of LCDM halos. This is likely caused in part by the X-ray selection bias in terms of the relaxation state. We determine the temperature-mass (Tx-M500) relation for a subset of 105 XXL clusters that have both measured HSC lensing masses and X-ray temperatures. The resulting Tx-M500 relation is consistent with the self-similar prediction. Our Tx-M500 relation agrees with the XXL DR1 results at group scales, but has a slightly steeper mass trend, implying a smaller mass scale in the cluster regime. The overall offset in the Tx-M500 relation is at the $1.5σ$ level, corresponding to a mean mass offset of (34\pm 20)%. We also provide bias-corrected, weak-lensing-calibrated M200 and M500 mass estimates of individual XXL clusters based on their measured X-ray temperatures. △ Less

Submitted 4 March, 2020; v1 submitted 23 September, 2019; originally announced September 2019.

Comments: Version matching the one published in ApJ. We recommend to use statistically corrected mass estimates (M200MT, M500MT) of Table 2 for a given individual cluster. One of two companion papers presenting initial HSC-XXL results (Mauro Sereno et al., arXiv:1912.02827)

Journal ref: ApJ, 890, 148 (2020)

Neutral iron line in the supernova remnant IC 443 and implications for MeV cosmic rays

Authors: Kumiko K. Nobukawa, Arisa Hirayama, Aika Shimaguchi, Yutaka Fujita, Masayoshi Nobukawa, Shigeo Yamauchi

Abstract: We report a discovery of bright blob-like enhancements of an Fe I K$α$ line in the northwest and the middle of the supernova remnant (SNR) IC 443. The distribution of the line emission is associated with molecular clouds interacting with the shock front, and is totally different from that of the plasma. The Fe I K$α$ line has a large equivalent width. The most plausible scenario for the origin of… ▽ More We report a discovery of bright blob-like enhancements of an Fe I K$α$ line in the northwest and the middle of the supernova remnant (SNR) IC 443. The distribution of the line emission is associated with molecular clouds interacting with the shock front, and is totally different from that of the plasma. The Fe I K$α$ line has a large equivalent width. The most plausible scenario for the origin of the line emission is that the MeV protons accelerated in the shell leak into the molecular clouds and ionize the Fe atoms therein. The observed Fe I K$α$ line intensity is consistent with the prediction of a theoretical model, in which MeV protons are accelerated along with GeV and TeV protons at the SNR. △ Less

Submitted 14 August, 2019; originally announced August 2019.

Comments: 13 pages, 3 figures, accepted for publication in PASJ

The ALMA Discovery of the Rotating Disk and Fast Outflow of Cold Molecular Gas in NGC 1275

Authors: H. Nagai, K. Onishi, N. Kawakatu, Y. Fujita, M. Kino, Y. Fukazawa, J. Lim, W. Forman, J. Vrtilek, K. Nakanishi, H. Noda, K. Asada, K. Wajima, Y. Ohyama, L. David

Abstract: We present ALMA Band 6 observations of the CO(2-1), HCN(3-2), and HCO$^{+}$(3-2) lines in the nearby radio galaxy / brightest cluster galaxy (BCG) of NGC 1275 with the spatial resolution of $\sim20$ pc. In the previous observations, CO(2-1) emission was detected as radial filaments lying in the east-west direction. We resolved the inner filament and found that the filament cannot be represented by… ▽ More We present ALMA Band 6 observations of the CO(2-1), HCN(3-2), and HCO$^{+}$(3-2) lines in the nearby radio galaxy / brightest cluster galaxy (BCG) of NGC 1275 with the spatial resolution of $\sim20$ pc. In the previous observations, CO(2-1) emission was detected as radial filaments lying in the east-west direction. We resolved the inner filament and found that the filament cannot be represented by a simple infalling stream both morphologically and kinematically. The observed complex nature of the filament resembles the cold gas structure predicted by recent numerical simulations of cold chaotic accretion. A crude estimate suggests that the accretion rate of the cold gas can be higher than that of hot gas. Within the central 100 pc, we detected a rotational disk of the molecular gas whose mass is $\sim10^{8} M_{\sun}$. This is the first evidence of the presence of massive cold gas disk on this spatial scale for BCGs. The disk rotation axis is approximately consistent with the axis of the radio jet on subpc scales. This probably suggests that the cold gas disk is physically connected to the innermost accretion disk which is responsible for jet launching. We also detected absorption features in the HCN(3-2) and HCO$^{+}$(3-2) spectra against the radio continuum emission mostly radiated by $\sim1.2$-pc size jet. The absorption features are blue-shifted from the systemic velocity by $\sim$300-600~km~s$^{-1}$, which suggests the presence of outflowing gas from the active galactic nucleus (AGN). We discuss the relation of the AGN feeding with cold accretion, the origin of blue-shifted absorption, and estimate of black hole mass using the molecular gas dynamics. △ Less

Submitted 26 August, 2019; v1 submitted 15 May, 2019; originally announced May 2019.

Comments: Version 2 (accepted version). 18 pages, 16 figures. Accepted for publication in ApJ

First detection of PSR B1259-63/LS 2883 in the Millimeter and Submillimeter Wavelengths with ALMA

Authors: Yutaka Fujita, Akiko Kawachi, Takuya Akahori, Hiroshi Nagai, Masaki Yamaguchi

Abstract: We report Atacama Large Millimeter/submillimeter Array (ALMA) observations of the binary system containing the pulsar PSR B1259-63 orbiting around a Be star LS 2883 after the 2017 periastron passage. We detected radio continuum emission from the binary system in the millimeter/submillimeter wavelengths for the first time. At Band 3 (97 GHz), the flux 84 days after the periastron is almost the same… ▽ More We report Atacama Large Millimeter/submillimeter Array (ALMA) observations of the binary system containing the pulsar PSR B1259-63 orbiting around a Be star LS 2883 after the 2017 periastron passage. We detected radio continuum emission from the binary system in the millimeter/submillimeter wavelengths for the first time. At Band 3 (97 GHz), the flux 84 days after the periastron is almost the same as that 71 days after the periastron. Although the binary system showed intense GeV gamma-ray flares during our observations, the Band 3 flux did not indicate any time correlation with them. The Band 3 fluxes are consistent with an extrapolation of the radio spectrum at lower frequencies. Assuming that it is synchrotron emission, we constrain magnetic fields ($\lesssim 0.6$ G) and the high-energy cutoff of the electrons ($γ\gtrsim 360$). The flux at Band 7 (343 GHz) 69 days after the periastron shows a significant excess from the extrapolation of the radio spectrum at lower frequencies. The flux may be associated with the circumstellar disk around the Be star. We also present the results of Australian Telescope Compact Array (ATCA) observations at 94 GHz for the 2014 periastron passage, which show that the radio spectrum was relatively soft when the pulsar passed the disk. △ Less

Submitted 8 July, 2019; v1 submitted 17 April, 2019; originally announced April 2019.

Comments: Accepted for publication in PASJ Letter

Report number: OU-TAP 440

Unveiling the Galaxy Cluster - Cosmic Web Connection with X-ray observations in the Next Decade

Authors: Stephen A. Walker, Daisuke Nagai, A. Simionescu, M. Markevitch, H. Akamatsu, M. Arnaud, C. Avestruz, M. Bautz, V. Biffi, S. Borgani, E. Bulbul, E. Churazov, K. Dolag, D. Eckert, S. Ettori, Y. Fujita, M. Gaspari, V. Ghirardini, R. Kraft, E. T. Lau, A. Mantz, K. Matsushita, M. McDonald, E. Miller, T. Mroczkowski , et al. (13 additional authors not shown)

Abstract: In recent years, the outskirts of galaxy clusters have emerged as one of the new frontiers and unique laboratories for studying the growth of large scale structure in the universe. Modern cosmological hydrodynamical simulations make firm and testable predictions of the thermodynamic and chemical evolution of the X-ray emitting intracluster medium. However, recent X-ray and Sunyaev-Zeldovich effect… ▽ More In recent years, the outskirts of galaxy clusters have emerged as one of the new frontiers and unique laboratories for studying the growth of large scale structure in the universe. Modern cosmological hydrodynamical simulations make firm and testable predictions of the thermodynamic and chemical evolution of the X-ray emitting intracluster medium. However, recent X-ray and Sunyaev-Zeldovich effect observations have revealed enigmatic disagreements with theoretical predictions, which have motivated deeper investigations of a plethora of astrophysical processes operating in the virialization region in the cluster outskirts. Much of the physics of cluster outskirts is fundamentally different from that of cluster cores, which has been the main focus of X-ray cluster science over the past several decades. A next-generation X-ray telescope, equipped with sub-arcsecond spatial resolution over a large field of view along with a low and stable instrumental background, is required in order to reveal the full story of the growth of galaxy clusters and the cosmic web and their applications for cosmology. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 10 pages, 3 figures, Science white paper submitted to the Astro2020 Decadal Survey

Interaction between Molecular Clouds and MeV-TeV Cosmic-ray Protons Escaped from Supernova Remnants

Authors: Ken Makino, Yutaka Fujita, Kumiko Nobukawa, Hironori Matsumoto, Yutaka Ohira

Abstract: Recent discovery of the X-ray neutral iron line (Fe I Kαat 6.40 keV) around several supernova remnants (SNRs) show that MeV cosmic-ray (CR) protons are distributed around the SNRs and are interacting with neutral gas there. We propose that these MeV CRs are the ones that have been accelerated at the SNRs together with GeV-TeV CRs. In our analytical model, the MeV CRs are still confined in the SNR… ▽ More Recent discovery of the X-ray neutral iron line (Fe I Kαat 6.40 keV) around several supernova remnants (SNRs) show that MeV cosmic-ray (CR) protons are distributed around the SNRs and are interacting with neutral gas there. We propose that these MeV CRs are the ones that have been accelerated at the SNRs together with GeV-TeV CRs. In our analytical model, the MeV CRs are still confined in the SNR when the SNR collides with molecular clouds. After the collision, the MeV CRs leak into the clouds and produce the neutral iron line emissions. On the other hand, GeV-TeV CRs had already escaped from the SNRs and emit gamma-rays through interaction with molecular clouds surrounding the SNRs. We apply this model to the SNRs W28 and W44 and show that it can reproduce the observations of the iron line intensities and the gamma-ray spectra. This can be another support of a hadronic scenario for the gamma-ray emissions from these SNRs. △ Less

Submitted 23 August, 2019; v1 submitted 29 January, 2019; originally announced January 2019.

Comments: Matches the version published on PASJ

Report number: OU-TAP 437

Halo Concentrations and the New Baseline X-ray Luminosity-Temperature and Mass Relations of Galaxy Clusters

Authors: Yutaka Fujita, Han Aung

Abstract: The standard self-similar model of galaxy cluster formation predicts that the X-ray luminosity-temperature ($L_X$-$T_X$) relation of galaxy clusters should have been $L_X\propto T_X^2$ in absence of the baryonic physics, such as radiative cooling and feedback from stars and black holes. However, this baseline relation is predicted without considering the fact that the halo concentration and the ch… ▽ More The standard self-similar model of galaxy cluster formation predicts that the X-ray luminosity-temperature ($L_X$-$T_X$) relation of galaxy clusters should have been $L_X\propto T_X^2$ in absence of the baryonic physics, such as radiative cooling and feedback from stars and black holes. However, this baseline relation is predicted without considering the fact that the halo concentration and the characteristic density of clusters increases as their mass decreases, which is a consequence of hierarchical structure formation of the universe. Here, we show that the actual baseline relation should be $L_X\propto T_X^α$, where $α\sim 1.7$, instead of $α=2$, given the mass dependence of the concentration and the fundamental plane relation of galaxy clusters. Numerical simulations show that $α\sim 1.6$, which is consistent with the prediction. We also show that the baseline luminosity-mass ($L_X$-$M_Δ$) relation should have been $L_X\propto M_Δ^β$, where $β\sim 1.1$-1.2, in contrast with the conventional prediction ($β=4/3$). In addition, some of the scatter in the $L_X$-$M_Δ$ relation can be attributed to the scatter in the concentration-mass ($c$-$M$) relation. The confirmation of the shallow slope could be a proof of hierarchical clustering. As an example, we show that the new baseline relations could be checked by studying the temperature or mass dependence of gas mass fraction of clusters. Moreover, the highest-temperature clusters would follow the shallow baseline relations if the influences of cool cores and cluster mergers are properly removed. △ Less

Submitted 11 April, 2019; v1 submitted 16 January, 2019; originally announced January 2019.

Comments: Matches the version to appear in ApJ

Report number: OU-TAP 436

Journal ref: The Astrophysical Journal, 875, 26 (2019)

Halo Concentrations and the Fundamental Plane of Galaxy Clusters

Authors: Yutaka Fujita, Megan Donahue, Stefano Ettori, Keiichi Umetsu, Elena Rasia, Massimo Meneghetti, Elinor Medezinski, Nobuhiro Okabe, Marc Postman

Abstract: According to the standard cold dark matter (CDM) cosmology, the structure of dark halos including those of galaxy clusters reflects their mass accretion history. Older clusters tend to be more concentrated than younger clusters. Their structure, represented by the characteristic radius $r_s$ and mass $M_s$ of the Navarro--Frenk--White (NFW) density profile, is related to their formation time. In~t… ▽ More According to the standard cold dark matter (CDM) cosmology, the structure of dark halos including those of galaxy clusters reflects their mass accretion history. Older clusters tend to be more concentrated than younger clusters. Their structure, represented by the characteristic radius $r_s$ and mass $M_s$ of the Navarro--Frenk--White (NFW) density profile, is related to their formation time. In~this study, we showed that $r_s$, $M_s$, and the X-ray temperature of the intracluster medium (ICM), $T_X$, form a thin plane in the space of $(\log r_s, \log M_s, \log T_X)$. This tight correlation indicates that the ICM temperature is also determined by the formation time of individual clusters. Numerical simulations showed that clusters move along the fundamental plane as they evolve. The plane and the cluster evolution within the plane could be explained by a similarity solution of structure formation of the universe. The angle of the plane shows that clusters have not achieved "virial equilibrium" in the sense that mass/size growth and pressure at the boundaries cannot be ignored. The distribution of clusters on the plane was related to the intrinsic scatter in the halo concentration--mass relation, which originated from the variety of cluster ages. The well-known mass--temperature relation of clusters ($M_Δ\propto T_X^{3/2}$) can be explained by the fundamental plane and the mass dependence of the halo concentration without the assumption of virial equilibrium. The fundamental plane could also be used for calibration of cluster masses. △ Less

Submitted 31 December, 2018; originally announced January 2019.

Comments: Invited review article, to be published in "From Dark Haloes to Visible Galaxies", special issue of Galaxies

Journal ref: Galaxies 2019, 7(1), 8

Discovery of Kolmogorov-like Magnetic Energy Spectrum in Tycho's Supernova Remnant by Two-point Correlations of Synchrotron Intensity

Authors: Jiro Shimoda, Takuya Akahori, A. Lazarian, Tsuyoshi Inoue, Yutaka Fujita

Abstract: The spectral slope of the magnetic energy in supernova remnants (SNRs) can be obtained by analysis of spatial two-point correlation functions of synchrotron intensities. This method has been originally developed for the analysis of magnetic field structure in diffuse interstellar medium and applied when the geometry of the emission region is simple and known. In this paper, applying this correlati… ▽ More The spectral slope of the magnetic energy in supernova remnants (SNRs) can be obtained by analysis of spatial two-point correlation functions of synchrotron intensities. This method has been originally developed for the analysis of magnetic field structure in diffuse interstellar medium and applied when the geometry of the emission region is simple and known. In this paper, applying this correlation analysis to Tycho's SNR, for which the synchrotron emission region is known to be a spherical shell, we find that the magnetic energy spectrum shows the Kolmogorov-like scaling. Our results can be explained by turbulence developed downstream of the shock front via Richtmyer-Meshkov instability or an amplification of upstream magnetic field induced by cosmic rays. They could be discriminated by future observations with a sub arcsecond resolution such as Square Kilometer Array. △ Less

Submitted 26 July, 2018; originally announced July 2018.

Comments: 6 pages, 3 figures, accepted by MNRAS

A New Interpretation of the Mass-Temperature Relation and Mass Calibration of Galaxy Clusters Based on the Fundamental Plane

Authors: Yutaka Fujita, Keiichi Umetsu, Stefano Ettori, Elena Rasia, Nobuhiro Okabe, Massimo Meneghetti

Abstract: Observations and numerical simulations have shown that the relation between the mass scaled with the critical density of the universe and the X-ray temperature of galaxy clusters is approximately represented by $M_Δ\propto T_X^{3/2}$ (e.g. $Δ=500$). This relation is often interpreted as evidence that clusters are in virial equilibrium. However, the recently discovered fundamental plane (FP) of clu… ▽ More Observations and numerical simulations have shown that the relation between the mass scaled with the critical density of the universe and the X-ray temperature of galaxy clusters is approximately represented by $M_Δ\propto T_X^{3/2}$ (e.g. $Δ=500$). This relation is often interpreted as evidence that clusters are in virial equilibrium. However, the recently discovered fundamental plane (FP) of clusters indicates that the temperature of clusters primarily depends on a combination of the characteristic mass $M_s$ and radius $r_s$ of the Navarro-Frenk-White profile rather than $M_Δ$. Moreover, the angle of the FP revealed that clusters are not in virial equilibrium because of continuous mass accretion from the surrounding matter. By considering both the FP and the mass dependence of the cluster concentration parameter, we show that this paradox can be solved and the relation $M_Δ\propto T_X^{3/2}$ actually reflects the central structure of clusters. We also find that the intrinsic scatter in the halo concentration-mass relation can largely account for the spread of clusters on the FP. We also show that X-ray data alone form the FP and the angle and the position are consistent with those of the FP constructed from gravitational lensing data. We demonstrate that a possible shift between the two FPs can be used to calibrate cluster masses obtained via X-ray observations. △ Less

Submitted 13 August, 2018; v1 submitted 13 April, 2018; originally announced April 2018.

Comments: Published on ApJ. Matched to published version

Report number: OU-TAP 420

Journal ref: The Astrophysical Journal, Volume 863, Issue 1, article id. 37, 9 pp. (2018)

Discovery of a new fundamental plane dictating galaxy cluster evolution from gravitational lensing

Authors: Yutaka Fujita, Keiichi Umetsu, Elena Rasia, Massimo Meneghetti, Megan Donahue, Elinor Medezinski, Nobuhiro Okabe, Marc Postman

Abstract: In cold dark matter (CDM) cosmology, objects in the Universe have grown under the effect of gravity of dark matter. The intracluster gas in a galaxy cluster was heated when the dark-matter halo formed through gravitational collapse. The potential energy of the gas was converted to thermal energy through this process. However, this process and the thermodynamic history of the gas have not been clea… ▽ More In cold dark matter (CDM) cosmology, objects in the Universe have grown under the effect of gravity of dark matter. The intracluster gas in a galaxy cluster was heated when the dark-matter halo formed through gravitational collapse. The potential energy of the gas was converted to thermal energy through this process. However, this process and the thermodynamic history of the gas have not been clearly characterized in connection with with the formation and evolution of the internal structure of dark-matter halos. Here, we show that observational CLASH data of high-mass galaxy clusters lie on a plane in the three-dimensional logarithmic space of their characteristic radius $r_s$, mass $M_s$, and X-ray temperature $T_X$ with a very small orthogonal scatter. The tight correlation indicates that the gas temperature was determined at a specific cluster formation time, which is encoded in $r_s$ and $M_s$. The plane is tilted with respect to $T_X \propto M_s/r_s$, which is the plane expected in case of simplified virial equilibrium. We show that this tilt can be explained by a similarity solution, which indicates that clusters are not isolated but continuously growing through matter accretion from their outer environments. Numerical simulations reproduce the observed plane and its angle. This result holds independently of the gas physics implemented in the code, revealing the fundamental origin of this plane. △ Less

Submitted 25 April, 2018; v1 submitted 11 January, 2018; originally announced January 2018.

Comments: Replaced with a revised version to match the ApJ accepted version

Report number: OU-TAP 413

Journal ref: Astrophys.J. 857 (2018) 118

Science with the Cherenkov Telescope Array

Authors: The Cherenkov Telescope Array Consortium, :, B. S. Acharya, I. Agudo, I. Al Samarai, R. Alfaro, J. Alfaro, C. Alispach, R. Alves Batista, J. -P. Amans, E. Amato, G. Ambrosi, E. Antolini, L. A. Antonelli, C. Aramo, M. Araya, T. Armstrong, F. Arqueros, L. Arrabito, K. Asano, M. Ashley, M. Backes, C. Balazs, M. Balbo, O. Ballester , et al. (558 additional authors not shown)

Abstract: The Cherenkov Telescope Array, CTA, will be the major global observatory for very high energy gamma-ray astronomy over the next decade and beyond. The scientific potential of CTA is extremely broad: from understanding the role of relativistic cosmic particles to the search for dark matter. CTA is an explorer of the extreme universe, probing environments from the immediate neighbourhood of black ho… ▽ More The Cherenkov Telescope Array, CTA, will be the major global observatory for very high energy gamma-ray astronomy over the next decade and beyond. The scientific potential of CTA is extremely broad: from understanding the role of relativistic cosmic particles to the search for dark matter. CTA is an explorer of the extreme universe, probing environments from the immediate neighbourhood of black holes to cosmic voids on the largest scales. Covering a huge range in photon energy from 20 GeV to 300 TeV, CTA will improve on all aspects of performance with respect to current instruments. The observatory will operate arrays on sites in both hemispheres to provide full sky coverage and will hence maximize the potential for the rarest phenomena such as very nearby supernovae, gamma-ray bursts or gravitational wave transients. With 99 telescopes on the southern site and 19 telescopes on the northern site, flexible operation will be possible, with sub-arrays available for specific tasks. CTA will have important synergies with many of the new generation of major astronomical and astroparticle observatories. Multi-wavelength and multi-messenger approaches combining CTA data with those from other instruments will lead to a deeper understanding of the broad-band non-thermal properties of target sources. The CTA Observatory will be operated as an open, proposal-driven observatory, with all data available on a public archive after a pre-defined proprietary period. Scientists from institutions worldwide have combined together to form the CTA Consortium. This Consortium has prepared a proposal for a Core Programme of highly motivated observations. The programme, encompassing approximately 40% of the available observing time over the first ten years of CTA operation, is made up of individual Key Science Projects (KSPs), which are presented in this document. △ Less

Submitted 21 January, 2018; v1 submitted 22 September, 2017; originally announced September 2017.

Comments: 213 pages, including references and glossary. Version 2: credits and references updated, some figures updated, and author list updated

Enhanced Polarized Emission from the One-Parsec-Scale Hotspot of 3C 84 as a Result of the Interaction with Clumpy Ambient Medium

Authors: H. Nagai, Y. Fujita, M. Nakamura, M. Orienti, M. Kino, K. Asada, G. Giovannini

Abstract: We present Very Long Baseline Array polarimetric observations of the innermost jet of 3C$\sim$84 (NGC$\sim$1275) at 43$\sim$GHz. A significant polarized emission is detected at the hotspot of the innermost re-started jet, which is located $\sim$1 pc south from the radio core. While the previous report presented a hotspot at the southern end of the western limb, the hotspot location has been moved… ▽ More We present Very Long Baseline Array polarimetric observations of the innermost jet of 3C$\sim$84 (NGC$\sim$1275) at 43$\sim$GHz. A significant polarized emission is detected at the hotspot of the innermost re-started jet, which is located $\sim$1 pc south from the radio core. While the previous report presented a hotspot at the southern end of the western limb, the hotspot location has been moved to the southern end of the eastern limb. Faraday rotation is detected within an entire bandwidth of the 43-GHz band. The measured rotation measure (RM) is at most (6.3$\pm$1.9)$\times10^{5}$$\sim$rad$\sim$m$^{-2}$ and might be slightly time variable on the timescale of a month by a factor of a few. Our measured RM and the RM previously reported by the CARMA and SMA observations cannot be consistently explained by the spherical accretion flow with a power-law profile. We propose that a clumpy/inhomogeneous ambient medium is responsible for the observed rotation measure. Using equipartition magnetic field, we derive the electron density of $2\times10^{4}$$\sim$cm$^{-3}$. Such an electron density is consistent with the cloud of narrow line emission region around the central engine. We also discuss the magnetic field configuration from black hole scale to pc scale and the origin of low polarization. △ Less

Submitted 19 September, 2017; originally announced September 2017.

Comments: 8 pages, 8 figures, accepted for publication in ApJ

Cherenkov Telescope Array Contributions to the 35th International Cosmic Ray Conference (ICRC2017)

Authors: F. Acero, B. S. Acharya, V. Acín Portella, C. Adams, I. Agudo, F. Aharonian, I. Al Samarai, A. Alberdi, M. Alcubierre, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves Batista, J. -P. Amans, E. Amato, L. Ambrogi, G. Ambrosi, M. Ambrosio, J. Anderson, M. Anduze, E. O. Angüner, E. Antolini, L. A. Antonelli, V. Antonuccio , et al. (1117 additional authors not shown)

Abstract: List of contributions from the Cherenkov Telescope Array Consortium presented at the 35th International Cosmic Ray Conference, July 12-20 2017, Busan, Korea. List of contributions from the Cherenkov Telescope Array Consortium presented at the 35th International Cosmic Ray Conference, July 12-20 2017, Busan, Korea. △ Less

Submitted 24 October, 2017; v1 submitted 11 September, 2017; originally announced September 2017.

Comments: Index of Cherenkov Telescope Array conference proceedings at the ICRC2017, Busan, Korea

Cosmic Magnetism in Centimeter and Meter Wavelength Radio Astronomy

Authors: Takuya Akahori, Hiroyuki Nakanishi, Yoshiaki Sofue, Yutaka Fujita, Kiyotomo Ichiki, Shinsuke Ideguchi, Osamu Kameya, Takahiro Kudoh, Yuki Kudoh, Mami Machida, Yoshimitsu Miyashita, Hiroshi Ohno, Takeaki Ozawa, Keitaro Takahashi, Motokazu Takizawa, Dai G. Yamazaki

Abstract: Magnetic field is ubiquitous in the Universe and it plays essential roles in various astrophysical phenomena, yet its real origin and evolution are poorly known. This article reviews current understanding of magnetic fields in the interstellar medium, the Milky Way Galaxy, external galaxies, active galactic nuclei, clusters of galaxies, and the cosmic web. Particularly, the review concentrates on… ▽ More Magnetic field is ubiquitous in the Universe and it plays essential roles in various astrophysical phenomena, yet its real origin and evolution are poorly known. This article reviews current understanding of magnetic fields in the interstellar medium, the Milky Way Galaxy, external galaxies, active galactic nuclei, clusters of galaxies, and the cosmic web. Particularly, the review concentrates on the achievements that have been provided by centimeter and meter wavelength radio observations. The article also introduces various methods to analyze linear polarization data, including synchrotron radiation, Faraday rotation, depolarization, and Faraday tomography. △ Less

Submitted 7 September, 2017; originally announced September 2017.

Comments: 40 pages, 16 figures, accepted for publication in PASJ

Suzaku Observations of the Outskirts of the Galaxy Cluster Abell 3395 including a Filament toward Abell 3391

Authors: Yuuki Sugawara, Motokazu Takizawa, Madoka Itahana, Hiroki Akamatsu, Yutaka Fujita, Takaya Ohashi, Yoshitaka Ishisaki

Abstract: The results of Suzaku observations of the outskirts of Abell 3395 including a large-scale structure filament toward Abell 3391 are presented. We measured temperature and abundance distributions from the southern outskirt of Abell 3395 to the north at the virial radius, where a filament structure has been found in the former X-ray and Sunyaev-Zel'dovich effect observations between Abell 3391 and 33… ▽ More The results of Suzaku observations of the outskirts of Abell 3395 including a large-scale structure filament toward Abell 3391 are presented. We measured temperature and abundance distributions from the southern outskirt of Abell 3395 to the north at the virial radius, where a filament structure has been found in the former X-ray and Sunyaev-Zel'dovich effect observations between Abell 3391 and 3395. The overall temperature structure is consistent with the universal profile proposed by Okabe et al.(2014) for relaxed clusters except for the filament region. A hint of the ICM heating is found between the two clusters, which might be due to the interaction of them in the early phase of a cluster merger. Although we obtained relatively low metal abundance of $Z=0.169^{+0.164+0.009+0.018 }_{-0.150-0.004-0.015 }$ solar, where the first, second, and third errors are statistical, cosmic X-ray background systematic, and non X-ray background systematic, respectively, at the virial radius in the filament, our results are still consistent with the former results of other clusters ($Z \sim 0.3$ solar) within errors. Therefore, our results are also consistent with the early enrichment scenario. We estimated Compton $y$ parameters only from X-ray results in the region between Abell 3391 and 3395 assuming a simple geometry. They are smaller than the previous SZ results with Planck satellite. The difference could be attributed to a more elaborate geometry such as a filament inclined to the line-of-sight direction, or underestimation of the X-ray temperature because of the unresolved multi-temperature structures or undetected hot X-ray emission of the shock heated gas. △ Less

Submitted 29 August, 2017; originally announced August 2017.

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

Temporal Evolution of the Gamma-ray Burst Afterglow Spectrum for an Observer: GeV--TeV Synchrotron Self-Compton Light Curve

Authors: Takuma Fukushima, Sho To, Katsuaki Asano, Yutaka Fujita

Abstract: We numerically simulate the gamma-ray burst (GRB) afterglow emission with a one-zone time-dependent code. The temporal evolutions of the decelerating shocked shell and energy distributions of electrons and photons are consistently calculated. The photon spectrum and light curves for an observer are obtained taking into account the relativistic propagation of the shocked shell and the curvature of… ▽ More We numerically simulate the gamma-ray burst (GRB) afterglow emission with a one-zone time-dependent code. The temporal evolutions of the decelerating shocked shell and energy distributions of electrons and photons are consistently calculated. The photon spectrum and light curves for an observer are obtained taking into account the relativistic propagation of the shocked shell and the curvature of the emission surface. We find that the onset time of the afterglow is significantly earlier than the previous analytical estimate. The analytical formulae of the shock propagation and light curve for the radiative case are also different from our results. Our results show that even if the emission mechanism is switching from synchrotron to synchrotron self-Compton, the gamma-ray light curves can be a smooth power-law, which agrees with the observed light curve and the late detection of a 32 GeV photon in GRB 130427A. The uncertainty of the model parameters obtained with the analytical formula is discussed, especially in connection with the closure relation between spectral index and decay index. △ Less

Submitted 15 July, 2017; v1 submitted 22 June, 2017; originally announced June 2017.

Comments: 11 pages, 11 figures, Accepted for publication in ApJ

Multiwavelength study of X-ray Luminous Clusters in the Hyper Suprime-Cam Subaru Strategic Program S16A field

Authors: Keita Miyaoka, Nobuhiro Okabe, Takao Kitaguchi, Masamune Oguri, Yasushi Fukazawa, Rachel Mandelbaum, Elinor Medezinski, Yasunori Babazaki, Atsushi J. Nishizawa, Takashi Hamana, Yen-Ting Lin, Hiroki Akamatsu, I-Non Chiu, Yutaka Fujita, Yuto Ichinohe, Yutaka Komiyama, Toru Sasaki, Motokazu Takizawa, Shutaro Ueda, Keiichi Umetsu, Jean Coupon, Chiaki Hikage, Akio Hoshino, Alexie Leauthaud, Kyoko Matsushita , et al. (12 additional authors not shown)

Abstract: We present a joint X-ray, optical and weak-lensing analysis for X-ray luminous galaxy clusters selected from the MCXC (Meta-Catalog of X-Ray Detected Clusters of Galaxies) cluster catalog in the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) survey field with S16A data, As a pilot study of our planned series papers, we measure hydrostatic equilibrium (H.E.) masses using XMM-Newton data for f… ▽ More We present a joint X-ray, optical and weak-lensing analysis for X-ray luminous galaxy clusters selected from the MCXC (Meta-Catalog of X-Ray Detected Clusters of Galaxies) cluster catalog in the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) survey field with S16A data, As a pilot study of our planned series papers, we measure hydrostatic equilibrium (H.E.) masses using XMM-Newton data for four clusters in the current coverage area out of a sample of 22 MCXC clusters. We additionally analyze a non-MCXC cluster associated with one MCXC cluster. We show that H.E. masses for the MCXC clusters are correlated with cluster richness from the CAMIRA catalog (Oguri et al. 2017), while that for the non-MCXC cluster deviates from the scaling relation. The mass normalization of the relationship between the cluster richness and H.E. mass is compatible with one inferred by matching CAMIRA cluster abundance with a theoretical halo mass function. The mean gas mass fraction based on H.E. masses for the MCXC clusters is $\langle f_{\rm gas} \rangle = 0.125\pm0.012$ at spherical overdensity $Δ=500$, which is $\sim80-90$ percent of the cosmic mean baryon fraction, $Ω_b/Ω_m$, measured by cosmic microwave background experiments. We find that the mean baryon fraction estimated from X-ray and HSC-SSP optical data is comparable to $Ω_b/Ω_m$. A weak-lensing shear catalog of background galaxies, combined with photometric redshifts, is currently available only for three clusters in our sample. Hydrostatic equilibrium masses roughly agree with weak-lensing masses, albeit with large uncertainty. This study demonstrates that further multiwavelength study for a large sample of clusters using X-ray, HSC-SSP optical and weak lensing data will enable us to understand cluster physics and utilize cluster-based cosmology. △ Less

Submitted 18 October, 2017; v1 submitted 18 May, 2017; originally announced May 2017.

Comments: 16 pages, 11 figures, accepted for publication in PASJ. Full resolution paper is available from http://home.hiroshima-u.ac.jp/m161855/hscssp1st_mcxc_miyaoka.pdf

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

Prospects for CTA observations of the young SNR RX J1713.7-3946

Authors: The CTA Consortium, :, F. Acero, R. Aloisio, J. Amans, E. Amato, L. A. Antonelli, C. Aramo, T. Armstrong, F. Arqueros, K. Asano, M. Ashley, M. Backes, C. Balazs, A. Balzer, A. Bamba, M. Barkov, J. A. Barrio, W. Benbow, K. Bernlöhr, V. Beshley, C. Bigongiari, A. Biland, A. Bilinsky, E. Bissaldi , et al. (359 additional authors not shown)

Abstract: We perform simulations for future Cherenkov Telescope Array (CTA) observations of RX~J1713.7$-$3946, a young supernova remnant (SNR) and one of the brightest sources ever discovered in very-high-energy (VHE) gamma rays. Special attention is paid to explore possible spatial (anti-)correlations of gamma rays with emission at other wavelengths, in particular X-rays and CO/H{\sc i} emission. We presen… ▽ More We perform simulations for future Cherenkov Telescope Array (CTA) observations of RX~J1713.7$-$3946, a young supernova remnant (SNR) and one of the brightest sources ever discovered in very-high-energy (VHE) gamma rays. Special attention is paid to explore possible spatial (anti-)correlations of gamma rays with emission at other wavelengths, in particular X-rays and CO/H{\sc i} emission. We present a series of simulated images of RX J1713.7$-$3946 for CTA based on a set of observationally motivated models for the gamma-ray emission. In these models, VHE gamma rays produced by high-energy electrons are assumed to trace the non-thermal X-ray emission observed by {\it XMM-Newton}, whereas those originating from relativistic protons delineate the local gas distributions. The local atomic and molecular gas distributions are deduced by the NANTEN team from CO and H{\sc i} observations. Our primary goal is to show how one can distinguish the emission mechanism(s) of the gamma rays (i.e., hadronic vs leptonic, or a mixture of the two) through information provided by their spatial distribution, spectra, and time variation. This work is the first attempt to quantitatively evaluate the capabilities of CTA to achieve various proposed scientific goals by observing this important cosmic particle accelerator. △ Less

Submitted 13 April, 2017; originally announced April 2017.

Comments: 19 pages, 5 figures. Accepted for publication in ApJ

Journal ref: The Astrophysical Journal, Volume 840, Number 2 (2017) p. 74

ATCA 16 cm Observation of CIZA J1358.9-4750: Implication of Merger Stage and Constraint on Non-Thermal Properties

Authors: Takuya Akahori, Yuichi Kato, Kazuhiro Nakazawa, Takeaki Ozawa, Liyi Gu, Motokazu Takizawa, Yutaka Fujita, Hiroyuki Nakanishi, Nobuhiro Okabe, Kazuo Makishima

Abstract: We report the Australia Telescope Compact Array 16 cm observation of CIZA J1358.9-4750. Recent X-ray studies imply that this galaxy cluster is composed of merging, binary clusters. Using the EW367 configuration, we found no significant diffuse radio emission in and around the cluster. An upper limit of the total radio power at 1.4 GHz is $\sim 1.1 \times 10^{22}$ Watt/Hz in $30$ square arcminutes… ▽ More We report the Australia Telescope Compact Array 16 cm observation of CIZA J1358.9-4750. Recent X-ray studies imply that this galaxy cluster is composed of merging, binary clusters. Using the EW367 configuration, we found no significant diffuse radio emission in and around the cluster. An upper limit of the total radio power at 1.4 GHz is $\sim 1.1 \times 10^{22}$ Watt/Hz in $30$ square arcminutes which is a typical size of radio relics. It is known that an empirical relation holds between the total radio power and X-ray luminosity of the host cluster. The upper limit is about one order of magnitude lower than the power expected from the relation. Very young ($\sim 70$ Myr) shocks with low Mach numbers ($\sim 1.3$), which are often seen at an early stage of merger simulations, are suggested by the previous X-ray observation. The shocks may generate cosmic-ray electrons with a steep energy spectrum, which is consistent with non-detection of bright ($>10^{23}$ Watt/Hz) relic in this 16 cm band observation. Based on the assumption of energy equipartition, the upper limit gives the magnetic-field strength below $0.68f(D_{\rm los}/{\rm 1~Mpc})^{-1}(γ_{\rm min}/200)^{-1}$ $μ$G, where $f$ is the cosmic-ray total energy density over the cosmic-ray electron energy density, $D_{\rm los}$ is the depth of the shock wave along the sightline and $γ_{\rm min}$ is the lower cutoff Lorentz factor of the cosmic-ray electron energy spectrum. △ Less

Submitted 27 March, 2018; v1 submitted 20 November, 2016; originally announced November 2016.

Comments: 11 pages, 8 figures, accepted for publication in PASJ

Contributions of the Cherenkov Telescope Array (CTA) to the 6th International Symposium on High-Energy Gamma-Ray Astronomy (Gamma 2016)

Authors: The CTA Consortium, :, A. Abchiche, U. Abeysekara, Ó. Abril, F. Acero, B. S. Acharya, C. Adams, G. Agnetta, F. Aharonian, A. Akhperjanian, A. Albert, M. Alcubierre, J. Alfaro, R. Alfaro, A. J. Allafort, R. Aloisio, J. -P. Amans, E. Amato, L. Ambrogi, G. Ambrosi, M. Ambrosio, J. Anderson, M. Anduze, E. O. Angüner , et al. (1387 additional authors not shown)

Abstract: List of contributions from the Cherenkov Telescope Array (CTA) Consortium presented at the 6th International Symposium on High-Energy Gamma-Ray Astronomy (Gamma 2016), July 11-15, 2016, in Heidelberg, Germany. List of contributions from the Cherenkov Telescope Array (CTA) Consortium presented at the 6th International Symposium on High-Energy Gamma-Ray Astronomy (Gamma 2016), July 11-15, 2016, in Heidelberg, Germany. △ Less

Submitted 17 October, 2016; originally announced October 2016.

Comments: Index of CTA conference proceedings for the Gamma 2016, Heidelberg, Germany

Diffuse gamma-ray emission from the Galactic center and implications of its past activities

Authors: Yutaka Fujita, Shigeo S. Kimura, Kohta Murase

Abstract: It has been indicated that low-luminosity active galactic nuclei (LLAGNs) are accelerating high-energy cosmic-ray (CR) protons in their radiatively inefficient accretion flows (RIAFs). If this is the case, Sagittarius A* (Sgr A*) should also be generating CR protons, because Sgr A* is a LLAGN. Based on this scenario, we calculate a production rate of CR protons in Sgr A* and their diffusion in the… ▽ More It has been indicated that low-luminosity active galactic nuclei (LLAGNs) are accelerating high-energy cosmic-ray (CR) protons in their radiatively inefficient accretion flows (RIAFs). If this is the case, Sagittarius A* (Sgr A*) should also be generating CR protons, because Sgr A* is a LLAGN. Based on this scenario, we calculate a production rate of CR protons in Sgr A* and their diffusion in the central molecular zone (CMZ) around Sgr A*. The CR protons diffusing in the CMZ create gamma-rays through pp interaction. We show that the gamma-ray luminosity and spectrum are consistent with observations if Sgr A* was active in the past. △ Less

Submitted 3 October, 2016; originally announced October 2016.

Comments: Proceedings of IAU Symposium 322 (The Multi-Messenger Astrophysics of the Galactic Centre)

Discovery of a new subparsec counterjet in NGC 1275: the inclination angle and the environment

Authors: Yutaka Fujita, Hiroshi Nagai

Abstract: We report the detection of a new feature at the centre of NGC 1275 in the Perseus cluster, hosting the radio source 3C 84. This feature emerges ~ 2 mas (~ 0.8 pc) north of the central core in recent 15 and 43 GHz VLBA images, and seems to be the counterjet to a known radio jet expanding to the south of the core. Apparently, the two jets were born through an outburst around 2005. From the ratio of… ▽ More We report the detection of a new feature at the centre of NGC 1275 in the Perseus cluster, hosting the radio source 3C 84. This feature emerges ~ 2 mas (~ 0.8 pc) north of the central core in recent 15 and 43 GHz VLBA images, and seems to be the counterjet to a known radio jet expanding to the south of the core. Apparently, the two jets were born through an outburst around 2005. From the ratio of the apparent lengths of the two jets from the core, we found that the jet angle to the line of sight is θ=65^\circ\pm 16^\circ, which is not much different from the angle of the outer jets generated by an activity around 1959 and constrains theories on gamma-ray emission from jets. The new northern jet has a strongly inverted spectrum in contrast with the southern jet. This suggests that the central black hole is surrounded by a subparsec-scale accretion disk with the density of >~ 10^5 cm^-3. The brightness of the counterjet suggests that the disk is highly inhomogeneous. The ambient gas density in the direction of the jet is ~ 8 cm^-3 if the current jet activity is similar to the past average. △ Less

Submitted 25 November, 2016; v1 submitted 13 September, 2016; originally announced September 2016.

Report number: OU-TAP 387

Journal ref: 2017, MNRAS, 465, L94

Probing WHIM around Galaxy Clusters with Fast Radio Bursts and the Sunyaev-Zel'dovich effect

Authors: Yutaka Fujita, Takuya Akahori, Keiichi Umetsu, Craig L. Sarazin, Ka-Wah Wong

Abstract: We propose a new method to probe the Warm Hot Intergalactic Medium (WHIM) beyond the virial radius (R_200) of a cluster of galaxies, where X-ray observations are not easily achievable. In this method, we use dispersion measures (DMs) of Fast Radio Bursts (FRBs) that appear behind the cluster and the Sunyaev-Zel'dovich (SZ) effect towards the cluster. The DMs reflect the density of the intracluster… ▽ More We propose a new method to probe the Warm Hot Intergalactic Medium (WHIM) beyond the virial radius (R_200) of a cluster of galaxies, where X-ray observations are not easily achievable. In this method, we use dispersion measures (DMs) of Fast Radio Bursts (FRBs) that appear behind the cluster and the Sunyaev-Zel'dovich (SZ) effect towards the cluster. The DMs reflect the density of the intracluster medium (ICM) including the WHIM. If we observe a sufficient number of FRBs in the direction of the cluster, we can derive the density profile from the DMs. Similarly, we can derive the pressure profile from the SZ effect. By combining the density and the pressure profiles, the temperature profile can be obtained. Based on mock observations of nearby clusters, we find that the density of the WHIM can be determined even at > 2 R_200 from the cluster center when FRB observations with the Square Kilometre Array (SKA) become available. The temperature can be derived out to r~ 1.5 R_200, and the radius is limited by the current sensitivity of SZ observations. △ Less

Submitted 7 November, 2016; v1 submitted 12 September, 2016; originally announced September 2016.

Comments: Accepted for publication in ApJ

Sagittarius A* as an Origin of the Galactic PeV Cosmic Rays?

Authors: Yutaka Fujita, Kohta Murase, Shigeo S. Kimura

Abstract: Supernova remnants (SNRs) have commonly been considered as a source of the observed PeV cosmic rays (CRs) or a Galactic PeV particle accelerator ("Pevatron"). In this work, we study Sagittarius A* (Sgr A*), which is the low-luminosity active galactic nucleus of the Milky Way Galaxy, as another possible canditate of the Pevatron, because it sometimes became very active in the past. We assume that a… ▽ More Supernova remnants (SNRs) have commonly been considered as a source of the observed PeV cosmic rays (CRs) or a Galactic PeV particle accelerator ("Pevatron"). In this work, we study Sagittarius A* (Sgr A*), which is the low-luminosity active galactic nucleus of the Milky Way Galaxy, as another possible canditate of the Pevatron, because it sometimes became very active in the past. We assume that a large number of PeV CRs were injected by Sgr A* at the outburst about 10^7 yr ago when the Fermi bubbles were created. We constrain the diffusion coefficient for the CRs in the Galactic halo on the condition that the CRs have arrived on the Earth by now, while a fairly large fraction of them have escaped from the halo. Based on a diffusion-halo model, we solve a diffusion equation for the CRs and compare the results with the CR spectrum on the Earth. The observed small anisotropy of the arrival directions of CRs may be explained if the diffusion coefficient in the Galactic disk is smaller than that in the halo. Our model predicts that a boron-to-carbon ratio should be energy-independent around the knee, where the CRs from Sgr A* become dominant. It is unlikely that the spectrum of the CRs accelerated at the outburst is represented by a power-law similar to the one for those responsible for the gamma-ray emission from the central molecular zone (CMZ) around the Galactic center. △ Less

Submitted 9 April, 2017; v1 submitted 31 March, 2016; originally announced April 2016.

Comments: Accepted for publication in JCAP

Report number: OU-TAP 384

Journal ref: JCAP 04 (2017) 037

Resolving 4-D Nature of Magnetism with Depolarization and Faraday Tomography: Japanese SKA Cosmic Magnetism Science

Authors: Takuya Akahori, Yutaka Fujita, Kiyotomo Ichiki, Shinsuke Ideguchi, Takahiro Kudoh, Yuki Kudoh, Mami Machida, Hiroyuki Nakanishi, Hiroshi Ohno, Takeaki Ozawa, Keitaro Takahashi, Motokazu Takizawa

Abstract: Magnetic fields play essential roles in various astronomical objects. Radio astronomy has revealed that magnetic fields are ubiquitous in our Universe. However, the real origin and evolution of magnetic fields is poorly proven. In order to advance our knowledge of cosmic magnetism in coming decades, the Square Kilometre Array (SKA) should have supreme sensitivity than ever before, which provides n… ▽ More Magnetic fields play essential roles in various astronomical objects. Radio astronomy has revealed that magnetic fields are ubiquitous in our Universe. However, the real origin and evolution of magnetic fields is poorly proven. In order to advance our knowledge of cosmic magnetism in coming decades, the Square Kilometre Array (SKA) should have supreme sensitivity than ever before, which provides numerous observation points in the cosmic space. Furthermore, the SKA should be designed to facilitate wideband polarimetry so as to allow us to examine sightline structures of magnetic fields by means of depolarization and Faraday Tomography. The SKA will be able to drive cosmic magnetism of the interstellar medium, the Milky Way, galaxies, AGN, galaxy clusters, and potentially the cosmic web which may preserve information of the primeval Universe. The Japan SKA Consortium (SKA-JP) Magnetism Science Working Group (SWG) proposes the project "Resolving 4-D Nature of Magnetism with Depolarization and Faraday Tomography", which contains ten scientific use cases. △ Less

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

Comments: 24 pages, 15 figures, Magnetism chapter of Japan SKA Consortium Science Book

Turbulent Cosmic-Ray Reacceleration and the Curved Radio Spectrum of the Radio Relic in the Sausage Cluster

Authors: Yutaka Fujita, Hiroki Akamatsu, Shigeo S. Kimura

Abstract: It has often been thought that the northern radio relic in the galaxy cluster CIZA J2242.8+5301 (the "Sausage" Cluster) is associated with cosmic ray (CR) electrons that are accelerated at a shock through the diffusive shock acceleration (DSA) mechanism. However, recent radio observations have shown that the radio spectrum is curved, which is inconsistent with the prediction of a simple DSA model.… ▽ More It has often been thought that the northern radio relic in the galaxy cluster CIZA J2242.8+5301 (the "Sausage" Cluster) is associated with cosmic ray (CR) electrons that are accelerated at a shock through the diffusive shock acceleration (DSA) mechanism. However, recent radio observations have shown that the radio spectrum is curved, which is inconsistent with the prediction of a simple DSA model. Moreover, the CR electron spectrum before being affected by radiative cooling seems to be too hard for DSA. In this study, we show that these facts are natural consequences if the electrons are reaccelerated in turbulence downstream of the shock. In this model, DSA is not the main mechanism for generating high-energy electrons. We find that the mean free path of the electrons should be much shorter than the Coulomb mean free path for efficient reacceleration. The scale of the turbulent eddies must be smaller than the width of the relic. We also predict hard X-ray spectra of inverse Compton scattering of photons. △ Less

Submitted 27 April, 2016; v1 submitted 23 February, 2016; originally announced February 2016.

Comments: Replaced with published version

Suzaku X-ray Observations of the Nearest Non-Cool Core Cluster, Antlia: Dynamically Young but with Remarkably Relaxed Outskirts

Authors: Ka-Wah Wong, Jimmy A. Irwin, Daniel R. Wik, Ming Sun, Craig L. Sarazin, Yutaka Fujita, Thomas H. Reiprich

Abstract: We present the results of seven Suzaku mosaic observations (>200 ks in total) of the nearest non-cool core cluster, the Antlia Cluster (or Group), beyond its degree-scale virial radius in its eastern direction. The temperature is consistent with the scaled profiles of many other clusters. Its pressure follows the universal profile. The density slope in its outskirts is significantly steeper than t… ▽ More We present the results of seven Suzaku mosaic observations (>200 ks in total) of the nearest non-cool core cluster, the Antlia Cluster (or Group), beyond its degree-scale virial radius in its eastern direction. The temperature is consistent with the scaled profiles of many other clusters. Its pressure follows the universal profile. The density slope in its outskirts is significantly steeper than that of the nearest cool core cluster (Virgo) with a similar temperature as Antlia, but shallower than those of the massive clusters. The entropy increases all the way out to R_200, which is consistent in value with the baseline model predicted by a gravity heating-only mechanism in the outskirts. Antlia is quite relaxed in this direction. However, the entropy inside ~R_500 is significantly higher than the baseline model, which is similar to many other nearby low mass clusters or groups. The enclosed gas-mass fraction does not exceed the cosmic value out to 1.3R_200. Thus, there is no evidence of significant gas clumping, electron-ion non-equipartition, or departure from the hydrostatic equilibrium approximation that are suggested to explain the entropy and gas fraction anomalies found in the outskirts of some massive clusters. We also present scaling relations for the gas fraction (f_gas,200), entropy (K_200), and temperature (T_500) using 22 groups and clusters with published data in the literature. The enclosed baryon fraction at R_200 is broadly consistent with the cosmic value. The power law slope of the K_200-T_500 relation is 0.638+/-0.205. The entropy deficit at R_200 cannot be fully accounted for by the bias or deviation in the gas fraction. △ Less

Submitted 25 September, 2016; v1 submitted 22 February, 2016; originally announced February 2016.

Comments: 24 pages, 16 figures, minor modifications to match the published version in ApJ

AGN Jet Power, Formation of X-ray Cavities, and FR I/II Dichotomy in Galaxy Clusters

Authors: Yutaka Fujita, Nozomu Kawakatu, Isaac Shlosman

Abstract: We investigate the ability of jets in active galactic nuclei (AGNs) to break out of the ambient gas with sufficiently large advance velocities. Using observationally estimated jet power, we analyze 28 bright elliptical galaxies in nearby galaxy clusters. Because the gas density profiles in the innermost regions of galaxies have not been resolved so far, we consider two extreme cases for temperatur… ▽ More We investigate the ability of jets in active galactic nuclei (AGNs) to break out of the ambient gas with sufficiently large advance velocities. Using observationally estimated jet power, we analyze 28 bright elliptical galaxies in nearby galaxy clusters. Because the gas density profiles in the innermost regions of galaxies have not been resolved so far, we consider two extreme cases for temperature and density profiles. We also follow two types of evolution for the jet cocoons: being driven by the pressure inside the cocoon (Fanaroff-Riley [FR] I type), and being driven by the jet momentum (FR II type). Our main result is that regardless of the assumed form of density profiles, jets with observed powers of <~ 10^{44} erg s^{-1} are not powerful enough to evolve as FR II sources. Instead, they evolve as FR I sources and appear to be decelerated below the buoyant velocities of the cocoons when jets were propagating through the central dense regions of the host galaxies. This explains the reason why FR I sources are more frequent than FR II sources in clusters. Furthermore, we predict the sizes of X-ray cavities from the observed jet powers and compare them with the observed ones --- they are consistent within a factor of two if the FR I type evolution is realized. Finally, we find that the jets with a power >~ 10^{44} erg s^{-1} are less affected by the ambient medium, and some of them, but not all, could serve as precursors of the FR II sources. △ Less

Submitted 26 February, 2016; v1 submitted 14 January, 2016; originally announced January 2016.

Comments: Replaced with published version

Suzaku observation of a high entropy cluster Abell 548W

Authors: Kazuhiro Nakazawa, Yuichi Kato, Liyi Gu, Madoka Kawaharada, Motokazu Takizawa, Yutaka Fujita, Kazuo Makishima

Abstract: Abell 548W, one of the galaxy clusters located in the Abell 548 region, has about an order of magnitude lower X-ray luminosity compared to ordinal clusters in view of the well known intracluster medium (ICM) temperature vs X-ray luminosity (kT-L_X) relation. The cluster hosts a pair of diffuse radio sources to the north west and north, both about 10' apart from the cluster center. They are candida… ▽ More Abell 548W, one of the galaxy clusters located in the Abell 548 region, has about an order of magnitude lower X-ray luminosity compared to ordinal clusters in view of the well known intracluster medium (ICM) temperature vs X-ray luminosity (kT-L_X) relation. The cluster hosts a pair of diffuse radio sources to the north west and north, both about 10' apart from the cluster center. They are candidate radio relics, frequently associated with merging clusters. A Suzaku deep observation with exposure of 84.4 ks was performed to search signatures for merging in this cluster. The XIS detectors successfully detected the ICM emission out to 16' from the cluster center. The temperature is ~3.6 keV around its center, and ~2 keV at the outermost regions. The hot region (~6 keV) aside the relic candidates shifted to the cluster center reported by XMM-Newton was not seen in the Suzaku data, although its temperature of 3.6 keV itself is higher than the average temperature of 2.5 keV around the radio sources. In addition, a signature of a cool (kT ~0.9 keV) component was found around the north west source. A marginal temperature jump at its outer-edge was also found, consistent with the canonical idea of shock acceleration origin of the radio relics. The cluster has among the highest central entropy of ~400 keV cm^2 and is one of the so-called low surface brightness clusters. Taking into account the fact that its shape itself is relatively circular and smooth and also its temperature structure is nearly flat, possible scenarios for merging is discussed. △ Less

Submitted 10 February, 2016; v1 submitted 3 December, 2015; originally announced December 2015.

Comments: 10 pages, 7 figures, accepted to be published in PASJ

Turbulent Cosmic Ray Reacceleration at Radio Relics and Halos in Clusters of Galaxies

Authors: Yutaka Fujita, Motokazu Takizawa, Ryo Yamazaki, Hiroki Akamatsu, Hiroshi Ohno

Abstract: Radio relics are synchrotron emission found on the periphery of galaxy clusters. From the position and the morphology, it is often believed that the relics are generated by cosmic-ray (CR) electrons accelerated at shocks through a diffusive shock acceleration (DSA) mechanism. However, some radio relics have harder spectra than the prediction of the standard DSA model. One example is observed in th… ▽ More Radio relics are synchrotron emission found on the periphery of galaxy clusters. From the position and the morphology, it is often believed that the relics are generated by cosmic-ray (CR) electrons accelerated at shocks through a diffusive shock acceleration (DSA) mechanism. However, some radio relics have harder spectra than the prediction of the standard DSA model. One example is observed in the cluster 1RXS J0603.3+4214, which is often called the ``Toothbrush Cluster''. Interestingly, the position of the relic is shifted from that of a possible shock. In this study, we show that these discrepancies in the spectrum and the position can be solved if turbulent (re)acceleration is very effective behind the shock. This means that for some relics turbulent reacceleration may be the main mechanism to produce high-energy electrons, contrary to the common belief that it is the DSA. Moreover, we show that for efficient reacceleration, the effective mean free path of the electrons has to be much smaller than their Coulomb mean free path. We also study the merging cluster 1E 0657-56 or the ``Bullet Cluster'', in which a radio relic has not been found at the position of the prominent shock ahead of the bullet. We indicate that a possible relic at the shock is obscured by the observed large radio halo that is generated by strong turbulence behind the shock. We propose a simple explanation if the morphological differences of radio emission among the Toothbrush, the Bullet, and the Sausage (CIZA J2242.8+5301) Clusters. △ Less

Submitted 27 December, 2015; v1 submitted 5 November, 2015; originally announced November 2015.

Comments: Replaced with published version

Report number: OU-TAP 377

Journal ref: ApJ, 815, 116 (2015)