Showing 1–50 of 279 results for author: Yamada, S

Updated Picture of the Active Galactic Nuclei with Dusty/Dust-free Gas Structures and Effects of the Radiation Pressure

Authors: Shoichiro Mizukoshi, Takeo Minezaki, Hiroaki Sameshima, Mitsuru Kokubo, Hirofumi Noda, Taiki Kawamuro, Satoshi Yamada, Takashi Horiuchi

Abstract: This study investigates the properties of two gas structures of X-ray selected active galactic nuclei (AGNs), that is, dusty and dust-free gas components, by separating them with the line-of-sight dust extinction ($A_V$) and the neutral gas column density ($N_{\mathrm{H}}$). The typical column density of the dusty and dust-free gas differs depending on the Seyfert type, indicating that both struct… ▽ More This study investigates the properties of two gas structures of X-ray selected active galactic nuclei (AGNs), that is, dusty and dust-free gas components, by separating them with the line-of-sight dust extinction ($A_V$) and the neutral gas column density ($N_{\mathrm{H}}$). The typical column density of the dusty and dust-free gas differs depending on the Seyfert type, indicating that both structures have anisotropic column density distributions. The number of targets with the dusty gas column density ($N_{\mathrm{H,d}}$) of $\log N_{\mathrm{H,d}}\ [\mathrm{cm^{-2}}]>23$ is much smaller than that with the same column density of the dust-free gas. This result indicates that the optically-thick part of the dusty gas structure is very thin. There are very few targets with a larger Eddington ratio ($f_{\mathrm{Edd}}$) than the effective Eddington limit of the dusty gas and the covering factor of the dusty gas with $22\leq \log N_{\mathrm{H,d}}\ [\mathrm{cm^{-2}}]<24$ exhibits a clear drop at the effective Eddington limit. These results support the scenario wherein the covering factor of the dusty torus decreases in a high Eddington ratio owing to the radiation-driven dusty gas outflow. The covering factor of the dust-free gas with the column density ($N_{\mathrm{H,df}}$) of $22\leq \log N_{\mathrm{H,df}}\ [\mathrm{cm^{-2}}]<24$ similarly exhibits the decrease in high Eddington ratio, although it may be owing to the dust-free gas outflow driven by certain other mechanisms than the radiation pressure. Finally, we propose an updated picture of the AGN gas structure based on our results and the literature. △ Less

Submitted 12 June, 2024; originally announced June 2024.

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

X-ray Winds In Nearby-to-distant Galaxies (X-WING) - I: Legacy Surveys of Galaxies with Ultrafast Outflows and Warm Absorbers in $z \sim 0$-$4$

Authors: Satoshi Yamada, Taiki Kawamuro, Misaki Mizumoto, Claudio Ricci, Shoji Ogawa, Hirofumi Noda, Yoshihiro Ueda, Teruaki Enoto, Mitsuru Kokubo, Takeo Minezaki, Hiroaki Sameshima, Takashi Horiuchi, Shoichiro Mizukoshi

Abstract: As an inaugural investigation under the X-ray Winds In Nearby-to-distant Galaxies (X-WING) program, we assembled a dataset comprising 132 active galactic nuclei (AGNs) spanning redshifts $z \sim 0$-$4$ characterized by blueshifted absorption lines indicative of X-ray winds. Through an exhaustive review of previous research, we compiled the outflow parameters for 573 X-ray winds, encompassing key a… ▽ More As an inaugural investigation under the X-ray Winds In Nearby-to-distant Galaxies (X-WING) program, we assembled a dataset comprising 132 active galactic nuclei (AGNs) spanning redshifts $z \sim 0$-$4$ characterized by blueshifted absorption lines indicative of X-ray winds. Through an exhaustive review of previous research, we compiled the outflow parameters for 573 X-ray winds, encompassing key attributes such as outflow velocities ($V_{\rm out}$), ionization parameters ($ξ$), and hydrogen column densities. By leveraging the parameters $V_{\rm out}$ and $ξ$, we systematically categorized the winds into three distinct groups: ultrafast outflows (UFOs), low-ionization parameter (low-IP) UFOs, and warm absorbers. Strikingly, a discernible absence of linear correlations in the outflow parameters, coupled with distributions approaching instrumental detection limits, was observed. Another notable finding was the identification of a velocity gap around $V_{\rm out} \sim 10,000~{\rm km~s^{-1}}$. This gap was particularly evident in the winds detected via absorption lines within the $\lesssim$2 keV band, indicating disparate origins for low-IP UFOs and warm absorbers. In cases involving Fe XXV/Fe XXVI lines, where the gap might be attributed to potential confusion between emission/absorption lines and the Fe K-edge, the possibility of UFOs and galactic-scale warm absorbers being disconnected is considered. An examination of the outflow and dust sublimation radii revealed a distinction: UFOs appear to consist of dust-free material, whereas warm absorbers likely comprise dusty gas. From 2024, the X-Ray Imaging and Spectroscopy Mission (XRISM) is poised to alleviate observational biases, providing insights into the authenticity of the identified gap, a pivotal question in comprehending AGN feedback from UFOs. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: 34 pages, 4 tables, 7 figures, accepted in ApJS. The CSV files of Tables B2-B4 are available from the URL: http://www.kusastro.kyoto-u.ac.jp/~styamada/yamada2024bApJS/XWING-PaperI_tables.zip

[O IV] and [Ne V]-weak AGNs Hidden by Compton-thick Material in Late Mergers

Authors: Satoshi Yamada, Yoshihiro Ueda, Taiki Kawamuro, Claudio Ricci, Yoshiki Toba, Masatoshi Imanishi, Takamitsu Miyaji, Atsushi Tanimoto, Kohei Ichikawa, Martin Herrera-Endoqui, Shoji Ogawa, Ryosuke Uematsu, Keiichi Wada

Abstract: We study "buried" active galactic nuclei (AGNs) almost fully covered by circumnuclear material in ultra-/luminous infrared galaxies (U/LIRGs), which show weak ionized lines from narrow line regions. Employing an indicator of [O IV] 25.89-um or [Ne V] 14.32-um line to 12-um AGN luminosity ratio, we find 17 buried AGN candidates that are [O IV]-weak ($L_{\rm [O\,IV]}$/$L_{\rm 12,AGN} \leq -$3.0) or… ▽ More We study "buried" active galactic nuclei (AGNs) almost fully covered by circumnuclear material in ultra-/luminous infrared galaxies (U/LIRGs), which show weak ionized lines from narrow line regions. Employing an indicator of [O IV] 25.89-um or [Ne V] 14.32-um line to 12-um AGN luminosity ratio, we find 17 buried AGN candidates that are [O IV]-weak ($L_{\rm [O\,IV]}$/$L_{\rm 12,AGN} \leq -$3.0) or [Ne V]-weak ($L_{\rm [Ne\,V]}$/$L_{\rm 12,AGN} \leq -$3.4) among 30 AGNs in local U/LIRGs. For the [O IV]-weak AGNs, we estimate their covering fractions of Compton-thick (CT; $N_{\rm H} \geq 10^{24}$ cm$^{-2}$) material with an X-ray clumpy torus model to be $f^{\rm (spec)}_{\rm CT} = 0.55\pm0.19$ on average. This value is consistent with the fraction of CT AGNs ($f^{\rm (stat)}_{\rm CT} = 53\pm12$%) among the [O IV]-weak AGNs in U/LIRGs and much larger than that in Swift/BAT AGNs ($23\pm6$%). The fraction of [O IV]-weak AGNs increases from $27^{+13}_{-10}$% (early) to $66^{+10}_{-12}$% (late mergers). Similar results are obtained with the [Ne V] line. The [O IV] or [Ne V]-weak AGNs in late mergers show larger $N_{\rm H}$ and Eddington ratios ($λ_{\rm Edd}$) than those of the Swift/BAT AGNs, and the largest $N_{\rm H}$ is $\gtrsim$10$^{25}$ cm$^{-2}$ at ${\log}λ_{\rm Edd} \sim -$1, close to the effective Eddington limit for CT material. These suggest that (1) the circumnuclear material in buried AGNs is regulated by the radiation force from high-$λ_{\rm Edd}$ AGNs on the CT obscurers, and (2) their dense material with large $f^{\rm (spec)}_{\rm CT}$ ($\sim 0.5 \pm 0.1$) in U/LIRGs is a likely cause of a unique structure of buried AGNs, whose amount of material may be maintained through merger-induced supply from their host galaxies. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: 13 pages, 10 figures, and 1 table, accepted for publication in ApJ

ALMA Lensing Cluster Survey: Full SED Analysis of z~0.5-6 Lensed Galaxies Detected with Millimeter Observations

Authors: Ryosuke Uematsu, Yoshihiro Ueda, Kotaro Kohno, Yoshiki Toba, Satoshi Yamada, Ian Smail, Hideki Umehata, Seiji Fujimoto, Bunyo Hatsukade, Yiping Ao, Franz Erik Bauer, Gabriel Brammer, Miroslava Dessauges-Zavadsky, Daniel Espada, Jean-Baptiste Jolly, Anton M. Koekemoer, Vasily Kokorev, Georgios E. Magdis, Masamune Oguri, Fengwu Sun

Abstract: Sub/millimeter galaxies are a key population for the study of galaxy evolution because the majority of star formation at high redshifts occurred in galaxies deeply embedded in dust. To search for this population, we have performed an extensive survey with ALMA, called the ALMA Lensing Cluster Survey (ALCS). This survey covers 133 arcmin^2 area and securely detects 180 sources at z~0.5-6 with a flu… ▽ More Sub/millimeter galaxies are a key population for the study of galaxy evolution because the majority of star formation at high redshifts occurred in galaxies deeply embedded in dust. To search for this population, we have performed an extensive survey with ALMA, called the ALMA Lensing Cluster Survey (ALCS). This survey covers 133 arcmin^2 area and securely detects 180 sources at z~0.5-6 with a flux limit of ~0.2 mJy at 1.2 mm (Fujimoto et al. 2023). Here we report the results of multi-wavelength spectral energy distribution (SED) analysis of the whole ALCS sample, utilizing the observed-frame UV to millimeter photometry. We find that the majority of the ALCS sources lie on the star-forming main sequence, with a smaller fraction showing intense starburst activities. The ALCS sample contains high infrared-excess sources IRX=log(Ldust/LUV)>1), including two extremely dust-obscured galaxies (IRX>5). We also confirm that the ALCS sample probes a broader range in lower dust mass than conventional SMG samples in the same redshift range. We identify six heavily obscured AGN candidates that are not detected in the archival Chandra data in addition to the three X-ray AGNs reported by Uematsu et al. (2023). The inferred AGN luminosity density shows a possible excess at z=2-3 compared with that determined from X-ray surveys below 10 keV. △ Less

Submitted 13 February, 2024; v1 submitted 8 February, 2024; originally announced February 2024.

Comments: 27 pages, 20 figures

High count rate effects in event processing for XRISM/Resolve x-ray microcalorimeter

Authors: Misaki Mizumoto, Masahiro Tsujimoto, Renata S. Cumbee, Megan E. Eckart, Yoshitaka Ishisaki, Caroline A. Kilbourne, Edmund Hodges-Kluck, Maurice A. Leutenegger, Frederick S. Porter, Makoto Sawada, Yoh Takei, Yuusuke Uchida, Shin'ya Yamada, the XRISM Resolve team

Abstract: The spectroscopic performance of x-ray instruments can be affected at high count rates. The effects and mitigation in the optical chain, such as x-ray attenuation filters or de-focusing mirrors, are widely discussed, but those in the signal chain are not. Using the Resolve x-ray microcalorimeter onboard the XRISM satellite, we discuss the effects observed during high count rate measurements and ho… ▽ More The spectroscopic performance of x-ray instruments can be affected at high count rates. The effects and mitigation in the optical chain, such as x-ray attenuation filters or de-focusing mirrors, are widely discussed, but those in the signal chain are not. Using the Resolve x-ray microcalorimeter onboard the XRISM satellite, we discuss the effects observed during high count rate measurements and how these can be modeled. We focus on three instrumental effects that impact performance at high count rate: CPU limit, pile up, and electrical cross talk. High count rate data were obtained during ground testing using the flight model instrument and a calibration x-ray source. A simulated observation of GX 13+1 is presented to illustrate how to estimate these effects based on these models for observation planning. The impact of these effects on high count rate observations is discussed. △ Less

Submitted 24 December, 2023; originally announced December 2023.

Comments: 16 pages, 20 figures, Proc. SPIE

Journal ref: Proc. SPIE 12181, Space Telescopes and Instrumentation 2022: Ultraviolet to Gamma Ray, 121815Z (31 August 2022)

Collisional and Fast Neutrino Flavor Instabilities in Two-dimensional Core-collapse Supernova Simulation with Boltzmann Neutrino Transport

Authors: Ryuichiro Akaho, Jiabao Liu, Hiroki Nagakura, Masamichi Zaizen, Shoichi Yamada

Abstract: We present a comprehensive study on the occurrences of the collisional flavor instability (CFI) and the fast flavor instability (FFI) of neutrinos based on a two-dimensional (2D) core-collapse supernova (CCSN) simulation performed with a Boltzmann radiation hydrodynamics code. We find that CFI occurs in a region with the baryon-mass density of… ▽ More We present a comprehensive study on the occurrences of the collisional flavor instability (CFI) and the fast flavor instability (FFI) of neutrinos based on a two-dimensional (2D) core-collapse supernova (CCSN) simulation performed with a Boltzmann radiation hydrodynamics code. We find that CFI occurs in a region with the baryon-mass density of $10^{10}\lesssim ρ\lesssim 10^{12}\,\mathrm{g}\,\mathrm{cm}^{-3}$, which is similar to the previous results in one-dimensional (1D) CCSN models. In contrast to 1D, however, the CFI region varies with time vigorously in the 2D model, whereas it had a quiescent structure in 1D. This is attributed to the fact that the turbulent flows advected from a gain region account for the temporal variations. Another noticeable difference from the 1D models is the appearance of resonance-like CFI where number densities of $ν_e$, $\barν_e$ nearly coincide each other. The CFI growth rate there is enhanced and can reach $\sim10^{-3}\,\mathrm{cm}^{-1}$. As for FFI, on the other hand, it appears in three different regions; (1) the region overlapped with the resonance-like CFI, (2) neutrino decoupling regions where $\barν_e$'s are strongly emitted, and (3) optically thin regions where neutral-current scatterings dominate over charged-current reactions. Although overall properties for FFI are consistent with previous studies, we find that the number of electron-neutrinos lepton number crossing (ELN crossing) temporary becomes multiple, which can be assessed accurately only by multi-angle treatments in neutrino transport. We find that the growth rate of FFI is always higher than CFI if both of them occur, which suggests that the former is dominant for the linear evolution. △ Less

Submitted 9 January, 2024; v1 submitted 19 November, 2023; originally announced November 2023.

Comments: published in PRD

BASS XLII: The relation between the covering factor of dusty gas and the Eddington ratio in nearby active galactic nuclei

Authors: C. Ricci, K. Ichikawa, M. Stalevski, T. Kawamuro, S. Yamada, Y. Ueda, R. Mushotzky, G. C. Privon, M. J. Koss, B. Trakhtenbrot, A. C. Fabian, L. C. Ho, D. Asmus, F. E. Bauer, C. S. Chang, K. K. Gupta, K. Oh, M. Powell, R. W. Pfeifle, A. Rojas, F. Ricci, M. J. Temple, Y. Toba, A. Tortosa, E. Treister , et al. (3 additional authors not shown)

Abstract: Accreting supermassive black holes (SMBHs) located at the center of galaxies are typically surrounded by large quantities of gas and dust. The structure and evolution of this circumnuclear material can be studied at different wavelengths, from the submillimeter to the X-rays. Recent X-ray studies have shown that the covering factor of the obscuring material tends to decrease with increasing Edding… ▽ More Accreting supermassive black holes (SMBHs) located at the center of galaxies are typically surrounded by large quantities of gas and dust. The structure and evolution of this circumnuclear material can be studied at different wavelengths, from the submillimeter to the X-rays. Recent X-ray studies have shown that the covering factor of the obscuring material tends to decrease with increasing Eddington ratio, likely due to radiative feedback on dusty gas. Here we study a sample of 549 nearby (z<0.1) hard X-ray (14-195 keV) selected non-blazar active galactic nuclei (AGN), and use the ratio between the AGN infrared and bolometric luminosity as a proxy of the covering factor. We find that, in agreement with what has been found by X-ray studies of the same sample, the covering factor decreases with increasing Eddington ratio. We also confirm previous findings which showed that obscured AGN typically have larger covering factors than unobscured sources. Finally, we find that the median covering factors of AGN located in different regions of the column density-Eddington ratio diagram are in good agreement with what would be expected from a radiation-regulated growth of SMBHs. △ Less

Submitted 11 November, 2023; v1 submitted 2 November, 2023; originally announced November 2023.

Comments: ApJ in press

Co-evolution and Nuclear Structure in the Dwarf Galaxy POX 52 Studied by Multi-wavelength Data From Radio to X-ray

Authors: Taiki Kawamuro, Claudio Ricci, Satoshi Yamada, Hirofumi Noda, Ruancun Li, Matthew J. Temple, Alessia Tortosa

Abstract: The nearby dwarf galaxy POX 52 at $z = 0.021$ hosts an active galactic nucleus (AGN) with a black-hole (BH) mass of $M_{\rm BH} \sim 10^{5-6} M_\odot$ and an Eddington ratio of $\sim$ 0.1-1. This object provides the rare opportunity to study both AGN and host-galaxy properties in a low-mass highly accreting system. To do so, we collected its multi-wavelength data from X-ray to radio. First, we con… ▽ More The nearby dwarf galaxy POX 52 at $z = 0.021$ hosts an active galactic nucleus (AGN) with a black-hole (BH) mass of $M_{\rm BH} \sim 10^{5-6} M_\odot$ and an Eddington ratio of $\sim$ 0.1-1. This object provides the rare opportunity to study both AGN and host-galaxy properties in a low-mass highly accreting system. To do so, we collected its multi-wavelength data from X-ray to radio. First, we construct a spectral energy distribution, and by fitting it with AGN and host-galaxy components, we constrain AGN-disk and dust-torus components. Then, while considering the AGN-disk emission, we decompose optical HST images. As a result, it is found that a classical bulge component is probably present, and its mass ($M_{\rm bulge}$) is consistent with an expected value from a local relation. Lastly, we analyze new quasi-simultaneous X-ray (0.2-30 keV) data obtained by NuSTAR and XMM-Newton. The X-ray spectrum can be reproduced by multi-color blackbody, warm and hot coronae, and disk and torus reflection components. Based on this, the spin is estimated to be $a_{\rm spin} = 0.998_{-0.814}$, which could suggest that most of the current BH mass was achieved by prolonged mass accretion. Given the presence of the bulge, POX 52 would have undergone a galaxy merger, while the $M_{\rm BH}$-$M_{\rm bulge}$ relation and the inferred prolonged accretion could suggest that AGN feedback occurred. Regarding the AGN structure, the spectral slope of the hot corona, its relative strength to the bolometric emission, and the torus structure are found to be consistent with Eddington-ratio dependencies found for nearby AGNs. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: 26 pages, 10 figures, 4 tables, accepted for publication in ApJ. The data in Table 4 are available from http://www.kusastro.kyoto-u.ac.jp/~kawamuro/arxiv/SEDplusComponents_u.txt

A comparative study on three modes of s-process nucleosynthesis in extremely metal-poor AGB stars

Authors: S. Yamada, T. Suada, Y. Komiya, M. Aikawa, M. Y. Fujimoto

Abstract: Carbon-enhanced metal-poor (CEMP) stars in the Galactic halo have a wide range of neutron-capture element abundance patterns. To identify their origin, we investigated three modes of $s$-process nucleosynthesis that have been proposed to operate in extremely metal-poor (EMP) Asymptotic Giant Branch (AGB) stars: the convective 13C burning, which occurs when hydrogen is engulfed by the helium flash… ▽ More Carbon-enhanced metal-poor (CEMP) stars in the Galactic halo have a wide range of neutron-capture element abundance patterns. To identify their origin, we investigated three modes of $s$-process nucleosynthesis that have been proposed to operate in extremely metal-poor (EMP) Asymptotic Giant Branch (AGB) stars: the convective 13C burning, which occurs when hydrogen is engulfed by the helium flash convection in low-mass AGB stars, the convective 22Ne burning, which occurs in the helium flash convection of intermediate-mass AGB stars, and the radiative 13C burning, which occurs in the $^{13}$C pocket that is formed during the inter-pulse periods. We show that the production of $s$-process elements per iron seed ($s$-process efficiency) does not depend on metallicity below $[{\rm Fe}/{\rm H}]=-2$, because 16O in the helium zone dominates the neutron poison. The convective 13C mode can produce a variety of $s$-process efficiencies for Sr, Ba and Pb, including the maxima observed among CEMP stars. The 22Ne mode only produce the lowest end of $s$-process efficiencies among CEMP models. We show that the combination of these two modes can explain the full range of observed enrichment of $s$-process elements in CEMP stars. In contrast, the 13C pocket mode can hardly explain the high level of enrichment observed in some CEMP stars, even if considering star-to-star variations of the mass of the 13C pocket. These results provide a basis for discussing the binary mass transfer origin of CEMP stars and their subgroups. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 22 pages, 13 figures, 1 table. Accepted for publication in Publications of the Astronomical Society of Japan

Universality of the Neutrino Collisional Flavor Instability in Core Collapse Supernovae

Authors: Jiabao Liu, Hiroki Nagakura, Ryuichiro Akaho, Akira Ito, Masamichi Zaizen, Shoichi Yamada

Abstract: Neutrinos are known to undergo flavor conversion processes among the three flavors. The fast flavor conversion (FFC) has been the central piece of flavor conversions taking place in core-collapse supernovae (CCSNe) due to its shorter timescale to the completion of flavor conversion compared to other types of flavor conversion. Although the ordinary collisions between neutrinos and matter were once… ▽ More Neutrinos are known to undergo flavor conversion processes among the three flavors. The fast flavor conversion (FFC) has been the central piece of flavor conversions taking place in core-collapse supernovae (CCSNe) due to its shorter timescale to the completion of flavor conversion compared to other types of flavor conversion. Although the ordinary collisions between neutrinos and matter were once thought to decohere neutrinos and thus damp flavor conversions, it was recently realized that they can also induce the flavor conversion. The linear analysis showed that the so-called collisional flavor instability or CFI occurs in the absence of FFC. In this paper, we investigate if CFI takes place in of the post-bounce core of CCSNe, using the results of spherically symmetric Boltzmann simulations of CCSNe for four progenitor models with different masses. We also provide a necessary (but not sufficient) condition of matter properties for the occurrence of CFI in optically thick and semi-transparent regions; baryon mass density ($ρ$), electron fraction ($Y_e$), and the degeneracy of electron-type neutrinos ($η_{ν_e}$) need to be $10^{10} {\rm g/cm^3} \lesssim ρ\lesssim 10^{12} {\rm g/cm^3}$, $Y_e\lesssim 0.4$, and $η_{ν_e} \lesssim 0.5$, respectively. This condition allows us to easily locate the place of possible CFI occurence without detailed stability analyses, which is useful for analyzing CFI in CCSN models phenomenologically △ Less

Submitted 19 November, 2023; v1 submitted 8 October, 2023; originally announced October 2023.

Comments: accepted in PRD

AGN number fraction in galaxy groups and clusters at z < 1.4 from the Subaru Hyper Suprime-Cam survey

Authors: Aoi Hashiguchi, Yoshiki Toba, Naomi Ota, Masamune Oguri, Nobuhiro Okabe, Yoshihiro Ueda, Masatoshi Imanishi, Satoshi Yamada, Tomotsugu Goto, Shuhei Koyama, Kianhong Lee, Ikuyuki Mitsuishi, Tohru Nagao, Atsushi J. Nishizawa, Akatoki Noboriguchi, Taira Oogi, Koki Sakuta, Malte Schramm, Mio Shibata, Yuichi Terashima, Takuji Yamashita, Anri Yanagawa, Anje Yoshimoto

Abstract: One of the key questions on active galactic nuclei (AGN) in galaxy clusters is how AGN could affect the formation and evolution of member galaxies and galaxy clusters in the history of the Universe. To address this issue, we investigate the dependence of AGN number fraction ($f_{\rm AGN}$) on cluster redshift ($z_{\rm cl}$) and distance from the cluster center ($R/R_{\rm 200}$). We focus on more t… ▽ More One of the key questions on active galactic nuclei (AGN) in galaxy clusters is how AGN could affect the formation and evolution of member galaxies and galaxy clusters in the history of the Universe. To address this issue, we investigate the dependence of AGN number fraction ($f_{\rm AGN}$) on cluster redshift ($z_{\rm cl}$) and distance from the cluster center ($R/R_{\rm 200}$). We focus on more than 27,000 galaxy groups and clusters at $0.1 < z_{\rm cl} < 1.4$ with more than 1 million member galaxies selected from the Subaru Hyper Suprime-Cam. By combining various AGN selection methods based on infrared (IR), radio, and X-ray data, we identify 2,688 AGN. We find that (i) $f_{\rm AGN}$ increases with $z_{\rm cl}$ and (ii) $f_{\rm AGN}$ decreases with $R/R_{\rm 200}$. The main contributors to the rapid increase of $f_{\rm AGN}$ towards high-$z$ and cluster center are IR- and radio-selected AGN, respectively. Those results indicate that the emergence of the AGN population depends on the environment and redshift, and galaxy groups and clusters at high-$z$ play an important role in AGN evolution. We also find that cluster-cluster mergers may not drive AGN activity in at least the cluster center, while we have tentative evidence that cluster-cluster mergers would enhance AGN activity in the outskirts of (particularly massive) galaxy clusters. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Comments: 16 pages, 21 figures, and 2 tables, accepted for publication in PASJ

Probing the origin of the two-component structure of broad line region by reverberation mapping of an extremely variable quasar

Authors: Shumpei Nagoshi, Fumihide Iwamuro, Satoshi Yamada, Yoshihiro Ueda, Yuto Oikawa, Masaaki Otsuka, Keisuke Isogai, Shin Mineshige

Abstract: The physical origins of quasar components, such as the broad line region (BLR) and dust torus, remain under debate. To gain insights into them, we focused on Changing-State Quasars (CSQs) which provide a unique perspective through structural changes associated with accretion disk state transitions. We targeted SDSS J125809.31+351943.0, an extremely variable CSQ, to study its central core structure… ▽ More The physical origins of quasar components, such as the broad line region (BLR) and dust torus, remain under debate. To gain insights into them, we focused on Changing-State Quasars (CSQs) which provide a unique perspective through structural changes associated with accretion disk state transitions. We targeted SDSS J125809.31+351943.0, an extremely variable CSQ, to study its central core structure and kinematics. We conducted reverberation mapping with optical spectroscopy to explore the structure of the BLR and estimate the black hole mass. The results from H$β$ reverberation mapping indicated a black hole mass of $10^{9.64^{+0.11}_{-0.20}}\rm{M_\odot}$. Additionally, we analyzed variations in the optical to X-ray spectral indices, $α_{\rm{ox}}$, before and after the state transition, to investigate the accretion disk. These variations in $α_{\rm{ox}}$ and the Eddington ratio (from 0.4 \% to 2.4 \%) exhibitied behavior similar to state transitions observed in X-ray binary systems. Spectral analysis of H$β$ revealed a predominantly double-peaked profile during dim periods, transitioning to include a single-peaked component as the quasar brightened, suggesting that H$β$ contains a mixture of two components. Each of these components has its distinct characteristics: the first is a double-peaked profile that remains stable despite changes in the accretion rate, while the second is a variable single-peaked profile. Using time lags from reverberation mapping, we estimated the spatial relationships between these BLR components, the accretion disk, and the dust torus. Our results suggest that the BLR consists of two distinct components, each differing in location and origin. △ Less

Submitted 5 February, 2024; v1 submitted 24 June, 2023; originally announced June 2023.

Comments: Accepted for publication in MNRAS. Comments are welcome

Richardson-Lucy deconvolution with a spatially Variant point-spread function of Chandra: Supernova Remnant Cassiopeia A as an Example

Authors: Yusuke Sakai, Shinya Yamada, Toshiki Sato, Ryota Hayakawa, Ryota Higurashi, Nao Kominato

Abstract: Richardson-Lucy (RL) deconvolution is one of the classical methods widely used in X-ray astronomy and other areas. Amid recent progress in image processing, RL deconvolution still leaves much room for improvement under a realistic situations. One direction is to include the positional dependence of a point-spread function (PSF), so-called RL deconvolution with a spatially variant PSF (RL… ▽ More Richardson-Lucy (RL) deconvolution is one of the classical methods widely used in X-ray astronomy and other areas. Amid recent progress in image processing, RL deconvolution still leaves much room for improvement under a realistic situations. One direction is to include the positional dependence of a point-spread function (PSF), so-called RL deconvolution with a spatially variant PSF (RL$_{\rm{sv}}$). Another is the method of estimating a reliable number of iterations and their associated uncertainties. We developed a practical method that incorporates the RL$_{\rm{sv}}$ algorithm and the estimation of uncertainties. As a typical example of bright and high-resolution images, the Chandra X-ray image of the supernova remnant Cassiopeia~A was used in this paper. RL$_{\rm{sv}}$ deconvolution enables us to uncover the smeared features in the forward/backward shocks and jet-like structures. We constructed a method to predict the appropriate number of iterations by using statistical fluctuation of the observed images. Furthermore, the uncertainties were estimated by error propagation from the last iteration, which was phenomenologically tested with the observed data. Thus, our method is a practically efficient framework to evaluate the time evolution of the remnants and their fine structures embedded in high-resolution X-ray images. △ Less

Submitted 26 June, 2023; v1 submitted 23 June, 2023; originally announced June 2023.

Comments: 9 pages, 8 figures

Broadband X-ray Spectral Study of Nuclear Structure in Local Obscured Radio Galaxies

Authors: Yuya Nakatani, Yoshihiro Ueda, Claudio Ricci, Koki Inaba, Shoji Ogawa, Kenta Setoguchi, Ryosuke Uematsu, Satoshi Yamada, Tomohiro Yoshitake

Abstract: Radio galaxies are a key population to understand the importance of relativistic jets in AGN feedback. We present the results of a systematic, broadband X-ray spectral analysis of hard X-ray selected radio galaxies to investigate their nuclear structures. In this study, we focus on the seven most radio-loud, X-ray obscured narrow line radio galaxies in the \textit{Swift}/BAT 70 month AGN catalog.… ▽ More Radio galaxies are a key population to understand the importance of relativistic jets in AGN feedback. We present the results of a systematic, broadband X-ray spectral analysis of hard X-ray selected radio galaxies to investigate their nuclear structures. In this study, we focus on the seven most radio-loud, X-ray obscured narrow line radio galaxies in the \textit{Swift}/BAT 70 month AGN catalog. The spectra from 0.5 keV up to 66 keV obtained with \textit{Suzaku} and \textit{NuSTAR} of six objects are newly analyzed here by utilizing the X-ray clumpy torus model (XCLUMPY), whereas we refer to Ogawa et al. (2021) for the results of Centaurus A. We find that these radio galaxies have similar torus covering fractions compared with radio quiet AGNs at the same Eddington ratios ($-3 < \log λ_{\rm Edd} < -1$). This result implies that (1) the torus structure is not an important factor that determines the presence of jets and (2) AGN jets have physically little effect on the torus. △ Less

Submitted 15 June, 2023; originally announced June 2023.

Comments: 10 pages (4 tables/14 figures in the main text), accepted in MNRAS

MASTER OT J030227.28+191754.5: a dwarf nova at a massive oxygen-neon white-dwarf system ?

Authors: Mariko Kimura, Kazumi Kashiyama, Toshikazu Shigeyama, Yusuke Tampo, Shinya Yamada, Teruaki Enoto

Abstract: We present timing and spectral analysis results of the {\it NICER} and {\it NuSTAR} observations of the dwarf nova MASTER OT J030227.28$+$191754.5 during the 2021--2022 outburst. The soft X-ray component was found to be dominated by blackbody radiation with a temperature of $\sim$30 eV and also showed prominent oxygen and neon emission lines. The blackbody luminosity exceeded 10$^{34}$ ergs s… ▽ More We present timing and spectral analysis results of the {\it NICER} and {\it NuSTAR} observations of the dwarf nova MASTER OT J030227.28$+$191754.5 during the 2021--2022 outburst. The soft X-ray component was found to be dominated by blackbody radiation with a temperature of $\sim$30 eV and also showed prominent oxygen and neon emission lines. The blackbody luminosity exceeded 10$^{34}$ ergs s$^{-1}$, which is consistent with theoretical predictions, and then decreased more than an order of magnitude in 3.5 days. The inferred abundances of oxygen and neon in the optically-thin coronal region surrounding the central white dwarf (WD) are several times higher than the respective solar values. Although inconclusive, the abundance enrichment may originate from the WD, indicating that it may be mainly composed of oxygen and neon. Assuming that the blackbody radiation comes from the belt-shaped boundary layer between the WD and the accretion disk, we estimated the WD radius to be $(2.9\pm1.1)\times10^{8}$ cm, which corresponds to the WD mass range of 1.15--1.34 $M_{\odot}$. If the accretion continues for another $\sim$Gyr, the WD may experience an accretion-induced collapse into a neutron star and form a so-called black-widow pulsar system. △ Less

Submitted 25 May, 2023; originally announced May 2023.

Comments: Accepted for publication in ApJ. 15 pages, 8 figures, 3 tables

Mechanical cryocooler noise observed in the ground testing of the Resolve X-ray microcalorimeter onboard XRISM

Authors: R. Imamura, H. Awaki, M. Tsujimoto, S. Yamada, F. S. Porter, C. A. Kilbourne, R. L. Kelley, Y. Takei

Abstract: Low-temperature detectors often use mechanical coolers as part of the cooling chain in order to reach sub-Kelvin operating temperatures. The microphonics noise caused by the mechanical coolers is a general and inherent issue for these detectors. We have observed this effect in the ground test data obtained with the Resolve instrument to be flown on the XRISM satellite. Resolve is a cryogenic X-ray… ▽ More Low-temperature detectors often use mechanical coolers as part of the cooling chain in order to reach sub-Kelvin operating temperatures. The microphonics noise caused by the mechanical coolers is a general and inherent issue for these detectors. We have observed this effect in the ground test data obtained with the Resolve instrument to be flown on the XRISM satellite. Resolve is a cryogenic X-ray microcalorimeter spectrometer with a required energy resolution of 7 eV at 6 keV. Five mechanical coolers are used to cool from ambient temperature to about 4 K: four two-stage Stirling coolers (STC) driven nominally at 15 Hz and a Joule-Thomson cooler (JTC) driven nominally at 52 Hz. In 2019, we operated the flight-model instrument for two weeks, in which we also obtained accelerometer data inside the cryostat at a low-temperature stage (He tank). X-ray detector and accelerometer data were obtained continuously while changing the JTC drive frequency, which produced a unique data set for investigating how the vibration from the cryocoolers propagates to the detector. In the detector noise spectra, we observed harmonics of both STCs and JTC. More interestingly, we also observed the low (<20 Hz) frequency beat between the 4'th JTC and 14'th STC harmonics and the 7'th JTC and the 23--24'th STC harmonics. We present here a description and interpretation of these measurements. △ Less

Submitted 2 March, 2023; originally announced March 2023.

Comments: J Low Temp Phys (2023)

A Systematic Study on the Resonance in Collisional Neutrino Flavor Instability

Authors: Jiabao Liu, Masamichi Zaizen, Shoichi Yamada

Abstract: Investigations on the resonance in the collisional flavor instability (CFI) of neutrinos, which were reported recently, are reported. We show that the resonance occurs not only for the isotropy-preserving modes as pointed out in the previous work but also for the isotropy-breaking modes and that it enhances the growth rate of CFI by orders of magnitude. Employing the linear analysis and nonlinear… ▽ More Investigations on the resonance in the collisional flavor instability (CFI) of neutrinos, which were reported recently, are reported. We show that the resonance occurs not only for the isotropy-preserving modes as pointed out in the previous work but also for the isotropy-breaking modes and that it enhances the growth rate of CFI by orders of magnitude. Employing the linear analysis and nonlinear numerical simulations in the two-flavor scheme and under the relaxation approximation for the collision term, we discuss the criterion for the resonance, its effect on the nonlinear evolution as well as the influences of homogeneity-breaking (k \ne 0) perturbations as well as of anisotropy in the background on the resonance. We will also touch the cohabitation of the resonance with the fast flavor conversion (FFC). △ Less

Submitted 31 May, 2023; v1 submitted 13 February, 2023; originally announced February 2023.

Comments: Accepted in PRD

ALMA Lensing Cluster Survey: Properties of Millimeter Galaxies Hosting X-ray Detected Active Galactic Nuclei

Authors: Ryosuke Uematsu, Yoshihiro Ueda, Kotaro Kohno, Satoshi Yamada, Yoshiki Toba, Seiji Fujimoto, Bunyo Hatsukade, Hideki Umehata, Daniel Espada, Fengwu Sun, Georgios E. Magdis, Vasily Kokorev, Yiping Ao

Abstract: We report the multi-wavelength properties of millimeter galaxies hosting X-ray detected active galactic nuclei (AGNs) from the ALMA Lensing Cluster Survey (ALCS). ALCS is an extensive survey of well-studied lensing clusters with ALMA, covering an area of 133 arcmin$^2$ over 33 clusters with a 1.2 mm flux-density limit of ${\sim}$60 $\mathrm{μJy}$ ($1σ$). Utilizing the archival data of Chandra, we… ▽ More We report the multi-wavelength properties of millimeter galaxies hosting X-ray detected active galactic nuclei (AGNs) from the ALMA Lensing Cluster Survey (ALCS). ALCS is an extensive survey of well-studied lensing clusters with ALMA, covering an area of 133 arcmin$^2$ over 33 clusters with a 1.2 mm flux-density limit of ${\sim}$60 $\mathrm{μJy}$ ($1σ$). Utilizing the archival data of Chandra, we identify three AGNs at $z=$1.06, 2.09, and 2.84 among the 180 millimeter sources securely detected in the ALCS (of which 155 are inside the coverage of Chandra). The X-ray spectral analysis shows that two AGNs are not significantly absorbed ($\log N_{\mathrm{H}}/\mathrm{cm}^{-2} < 23$), while the other shows signs of moderate absorption ($\log N_{\mathrm{H}}/\mathrm{cm}^{-2}\sim 23.5$). We also perform spectral energy distribution (SED) modelling of X-ray to millimeter photometry. We find that our X-ray AGN sample shows both high mass accretion rates (intrinsic 0.5--8 keV X-ray luminosities of ${\sim}10^{\text{44--45}}\,\mathrm{erg\ s^{-1}}$) and star-formation rates (${\gtrsim}100\,M_{\odot}\,\mathrm{yr}^{-1}$). This demonstrates that a wide-area survey with ALMA and Chandra can selectively detect intense growth of both galaxies and supermassive black holes (SMBHs) in the high-redshift universe. △ Less

Submitted 23 January, 2023; originally announced January 2023.

Comments: 14 pages, 3 figures, 2 tables

Hard X-Ray to Radio Multiwavelength SED Analysis of Local U/LIRGs in GOALS Sample with Self-consistent AGN Model Including Polar-dust Component

Authors: Satoshi Yamada, Yoshihiro Ueda, Martín Herrera-Endoqui, Yoshiki Toba, Takamitsu Miyaji, Shoji Ogawa, Ryosuke Uematsu, Atsushi Tanimoto, Masatoshi Imanishi, Claudio Ricci

Abstract: We conduct a hard X-ray to radio multiwavelength spectral energy distribution (SED) decomposition for 57 local luminous and ultraluminous infrared galaxies (U/LIRGs) observed with Nuclear Spectroscopic Telescope Array and/or Swift/Burst Alert Telescope in GOALS (Armus et al. 2009) sample. We modify the latest SED-fitting code X-CIGALE by implementing the infrared (IR) CLUMPY model, allowing the mu… ▽ More We conduct a hard X-ray to radio multiwavelength spectral energy distribution (SED) decomposition for 57 local luminous and ultraluminous infrared galaxies (U/LIRGs) observed with Nuclear Spectroscopic Telescope Array and/or Swift/Burst Alert Telescope in GOALS (Armus et al. 2009) sample. We modify the latest SED-fitting code X-CIGALE by implementing the infrared (IR) CLUMPY model, allowing the multiwavelength study with the X-ray torus model (XCLUMPY) self-consistently. Adopting the torus parameters obtained by the X-ray fitting (Yamada et al. 2021), we estimate the properties of host galaxies, active galactic nucleus (AGN) tori, and polar dust. The star formation rates (SFRs) become larger with merger stage and most of them are above the main sequence. The SFRs are correlated with radio luminosity, indicating starburst emission is dominant in the radio band. Although polar-dust extinction is much smaller than torus extinction, the UV-to-IR (mainly IR) polar dust luminosities are $\sim$2 times larger than the torus ones. The polar-dust temperature decreases while the physical size, estimated by the temperature and dust sublimation radius, increases with AGN luminosity from a few tens of parsec (early mergers) to kiloparsec scales (late mergers), where the polar dust is likely the expanding (i.e., evolving) dusty outflows. The comparison between SFRs and intrinsic AGN luminosities suggests that the starbursts occur first and AGNs arise later, and overall their growth rates follow the simultaneous coevolution relation that can establish the local galaxy-SMBH mass relation. We confirm the coexistence of intense starbursts, AGNs, and large-scale outflows in late mergers, supporting a standard AGN feedback scenario. △ Less

Submitted 9 January, 2023; originally announced January 2023.

Comments: 84 pages (5 tables/29 figures in the main text and 8 tables/18 figures in the Appendix), accepted in ApJS

Narrow Fe-K$α$ Reverberation Mapping Unveils the Deactivated Broad-Line Region in a Changing-Look Active Galactic Nucleus

Authors: Hirofumi Noda, Taisei Mineta, Takeo Minezaki, Hiroaki Sameshima, Mitsuru Kokubo, Taiki Kawamuro, Satoshi Yamada, Takashi Horiuchi, Hironori Matsumoto, Makoto Watanabe, Kumiko Morihana, Yoichi Itoh, Koji S. Kawabata, Yasushi Fukazawa

Abstract: "Changing-look active galactic nuclei" (CLAGNs) are known to change their apparent types between types 1 and 2, usually accompanied by a drastic change in their luminosity on timescales of years. However, it is still unclear whether materials in broad-line regions (BLRs) in CLAGNs appear and disappear during the type-transition or remain at the same location while the line production is simply act… ▽ More "Changing-look active galactic nuclei" (CLAGNs) are known to change their apparent types between types 1 and 2, usually accompanied by a drastic change in their luminosity on timescales of years. However, it is still unclear whether materials in broad-line regions (BLRs) in CLAGNs appear and disappear during the type-transition or remain at the same location while the line production is simply activated or deactivated. Here we present our X-ray-optical monitoring results of a CLAGN, NGC 3516, by Suzaku, Swift, and ground telescopes, with our primary focus on the narrow Fe-K$α$ emission line, which is an effective probe of the BLR materials. We detected significant variations of the narrow Fe-K$α$ line on a timescale of tens of days during the type-2 (faint) phase in 2013-2014, and conducted "narrow Fe-K$α$ reverberation mapping," comparing its flux variation with those of the X-ray continuum from a corona and $B$-band continuum from an accretion disk. We derived, as a result, a time lag of $10.1^{+5.8}_{-5.6}$ days ($1σ$ errors) for the Fe-K$α$ line behind the continuum, which is consistent with the location of the BLR determined in optical spectroscopic reverberation mapping during the type-1 (bright) phase. This finding shows that the BLR materials remained at the same location without emitting optical broad-lines during the type-2 phase. Considering the drastic decrease of the radiation during the type-transition, our result is possibly inconsistent with the hotly-discussed formation models of the BLR which propose that the radiative pressure from an accretion disk should be the main driving force. △ Less

Submitted 5 December, 2022; originally announced December 2022.

Comments: 15 pages, 7 figures, 3 tables, accepted for the publication in ApJ

Leptonic and semi-leptonic neutrino interactions with muons in the proto-neutron star cooling

Authors: Ken'ichi Sugiura, Shun Furusawa, Kohsuke Sumiyoshi, Shoichi Yamada

Abstract: It is known that muons are scarce just after the birth of a proto-neutron star via a supernova explosion but get more abundant as the proto-neutron star cools via neutrino emissions on the Kelvin-Helmholtz timescale. In this paper we evaluate all the relevant rates of the neutrino interactions with muons at different times in the proto-neutron star cooling. We are particularly interested in th… ▽ More It is known that muons are scarce just after the birth of a proto-neutron star via a supernova explosion but get more abundant as the proto-neutron star cools via neutrino emissions on the Kelvin-Helmholtz timescale. In this paper we evaluate all the relevant rates of the neutrino interactions with muons at different times in the proto-neutron star cooling. We are particularly interested in the late phase ($ t \gtrsim 10 \operatorname{s}$), which will be accessible in the next Galactic supernova but has not been studied well so far. We calculate both leptonic and semi-leptonic processes, for the latter of which we pay attention also to the form factors with their dependence on the transferred momentum as well as to the modification of the dispersion relations for nucleons on the mean field level. We find that the flavor-exchange reactions $ν_e + μ^- \rightarrow ν_μ + e^-$ and $\barν_μ + μ^- \rightarrow \barν_e + e^-$ can be dominant, particularly at low energies, over the capture of $ν_e$ on neutron and the scatterings of $\barν_μ$ on nucleons as the opacity sources for these species and that the inverse muon decay $ \barν_e + ν_μ + e^- \leftrightarrows μ^- $ can overwhelm the scatterings of $\barν_e$ and $ν_μ$ on nucleons again at low energies. At high energies, on the other hand, the corrections in the semi-leptonic processes mentioned above are more important. We also show the non-trivial energy- and angular dependences of the flavor-exchange reactions and the inverse muon decay. In the study of the diffusion coefficients from these reactions, we find that $\barν_μ$ is most affected. These pieces of information are indispensable for numerical computations and the interpretation of results thereof for the proto-neutron star cooling particularly at the very late phase. △ Less

Submitted 13 November, 2022; originally announced November 2022.

Comments: 53 pages, 19 Figures

Journal ref: Progress of Theoretical and Experimental Physics, Volume 2022, Issue 11, November 2022, 113E01

A Lagrangian construction of rotating stars

Authors: Misa Ogata, Hirotada Okawa, Kotaro Fujisawa, Nobutoshi Yasutake, Yu Yamamoto, Shoichi Yamada

Abstract: We present a new formulation for numerically obtaining axisymmetric equilibrium structures of rotating stars in two spatial dimensions. With a view to apply it to the secular evolution of rotating stars, we base it on the Lagrangian description, i.e., we solve the force-balance equations to find the spatial positions of fluid elements endowed individually with a mass, specific entropy and angular… ▽ More We present a new formulation for numerically obtaining axisymmetric equilibrium structures of rotating stars in two spatial dimensions. With a view to apply it to the secular evolution of rotating stars, we base it on the Lagrangian description, i.e., we solve the force-balance equations to find the spatial positions of fluid elements endowed individually with a mass, specific entropy and angular momentum. The system of nonlinear equations obtained by finite-differencing the basic equations are solved with the W4 method, which is a new multi-dimensional root-finding scheme of our own devising. We augment it with a remapping scheme to avoid distortions of the Lagrangian coordinates. In this first one of a series of papers, we will give a detailed description of these methods initially. We then present the results of some test calculations, which include the construction of both rapidly rotating barotropic and baroclinic equilibrium states. We gauge their accuracies quantitatively with some diagnostic quantities as well as via comparisons with the counterparts obtained with an Eulerian code. For a demonstrative purpose, we apply the code to a toy-model cooling calculation of a rotating white dwarf. △ Less

Submitted 17 October, 2022; originally announced October 2022.

Comments: 15pages, 14 figures, submitted to MNRAS

Broadband X-ray Spectral Analysis of the Dual AGN System Mrk 739

Authors: Koki Inaba, Yoshihiro Ueda, Satoshi Yamada, Shoji Ogawa, Ryosuke Uematsu, Atsushi Tanimoto, Claudio Ricci

Abstract: We present the result of a broadband (0.5-70 keV) X-ray spectral analysis of the late-merger galaxy Mrk 739, which contains a dual active galactic nucleus (AGN), Mrk 739E and Mrk 739W, with a separation of $\sim$3.4 kpc. The spectra obtained with NuSTAR, Chandra, XMM-Newton and Swift/BAT are simultaneously analyzed by separating the contributions from the two AGNs and extended emission with the Ch… ▽ More We present the result of a broadband (0.5-70 keV) X-ray spectral analysis of the late-merger galaxy Mrk 739, which contains a dual active galactic nucleus (AGN), Mrk 739E and Mrk 739W, with a separation of $\sim$3.4 kpc. The spectra obtained with NuSTAR, Chandra, XMM-Newton and Swift/BAT are simultaneously analyzed by separating the contributions from the two AGNs and extended emission with the Chandra data. To evaluate the reflection components from the AGN tori, we consider two models, a phenomenological one (pexrav and zgauss) and a more physically motivated one (XCLUMPY; Tanimoto et al. 2019). On the basis of the results with XCLUMPY, we find that the AGNs in Mrk 739E and Mrk 739W have intrinsic 2-10 keV luminosities of $1.0 \times 10^{43}$ and $7.5 \times 10^{41}\ \rm{erg}\ \rm{s}^{-1}$ absorbed by hydrogen column densities of $N_{\rm{H}} < 6.5 \times 10^{19}\ \rm{cm}^{-2}$ and $N_{\rm{H}} = 6.9^{+3.2}_{-1.7} \times 10^{21}\ \rm{cm}^{-2}$, respectively. The torus covering fraction of the material with $N_{\rm{H}} > 10^{22} \rm{cm}^{-2}$ in Mrk 739E, $C_{\rm{T}}^{(22)} < 0.50$ at a 90% confidence limit, is found to be smaller than those found for late-merger ultra/luminous infrared galaxies, $C_{\rm{T}}^{(22)} = 0.71\pm0.16$ (mean and standard deviation; Yamada et al. 2021). Considering the small star formation rate of Mrk 739E, we suggest that the gas-to-mass ratio of the host galaxy is an important parameter to determine the circumnuclear environment of an AGN in late merger. △ Less

Submitted 17 October, 2022; originally announced October 2022.

Comments: 14 pages, 4 figures, 3 tables, accepted for publication in ApJ

Optical IFU Observations of GOALS Sample with KOOLS-IFU on Seimei Telescope: Initial results of 9 U/LIRGs at $z <$ 0.04

Authors: Yoshiki Toba, Satoshi Yamada, Kazuya Matsubayashi, Koki Terao, Aoi Moriya, Yoshihiro Ueda, Kouji Ohta, Aoi Hashiguchi, Kazuharu G. Himoto, Hideyuki Izumiura, Kazuma Joh, Nanako Kato, Shuhei Koyama, Hiroyuki Maehara, Rana Misato, Akatoki Noboriguchi, Shoji Ogawa, Naomi Ota, Mio Shibata, Nozomu Tamada, Anri Yanagawa, Naoki Yonekura, Tohru Nagao, Masayuki Akiyama, Masaru Kajisawa , et al. (1 additional authors not shown)

Abstract: We present ionized gas properties of 9 local ultra/luminous infrared galaxies (U/LIRGs) at $z <$ 0.04 through IFU observations with KOOLS-IFU on Seimei Telescope. The observed targets are drawn from the Great Observatories All-sky LIRG Survey (GOALS), covering a wide range of merger stages. We successfully detect emission lines such as H$β$, [OIII]$λ$5007, H$α$, [NII]$λλ$6549,6583, and [SII]$λλ$67… ▽ More We present ionized gas properties of 9 local ultra/luminous infrared galaxies (U/LIRGs) at $z <$ 0.04 through IFU observations with KOOLS-IFU on Seimei Telescope. The observed targets are drawn from the Great Observatories All-sky LIRG Survey (GOALS), covering a wide range of merger stages. We successfully detect emission lines such as H$β$, [OIII]$λ$5007, H$α$, [NII]$λλ$6549,6583, and [SII]$λλ$6717,6731 with a spectral resolution of $R$ = 1500-2000, which provides (i) spatially-resolved ($\sim$200-700 pc) moment map of ionized gas and (ii) diagnostics for active galactic nucleus (AGN) within the central $\sim$3--11 kpc in diameter for our sample. We find that [OIII] outflow that is expected to be driven by AGN tends to be stronger (i) towards the galactic center and (ii) as a sequence of merger stage. In particular, the outflow strength in the late-stage (stage D) mergers is about 1.5 times stronger than that in the early-state (stage B) mergers, which indicates that galaxy mergers could induce AGN-driven outflow and play an important role in the co-evolution of galaxies and supermassive black holes. △ Less

Submitted 24 August, 2022; originally announced August 2022.

Comments: 12 pages, 8 figures, and 2 tables, accepted for publication in PASJ

Protoneutron Star Convection Simulated with a New General Relativistic Boltzmann Neutrino Radiation-Hydrodynamics Code

Authors: Ryuichiro Akaho, Akira Harada, Hiroki Nagakura, Wakana Iwakami, Hirotada Okawa, Shun Furusawa, Hideo Matsufuru, Kohsuke Sumiyoshi, Shoichi Yamada

Abstract: We investigate the protoneutron star (PNS) convection using our newly developed general relativistic Boltzmann neutrino radiation-hydrodynamics code. This is a pilot study for more comprehensive investigations later. As such, we take a snapshot of a PNS at 2.3 seconds after bounce from a 1D PNS cooling calculation and run our simulation for $\sim160\,\mathrm{ms}$ in 2D under axisymmetry. The origi… ▽ More We investigate the protoneutron star (PNS) convection using our newly developed general relativistic Boltzmann neutrino radiation-hydrodynamics code. This is a pilot study for more comprehensive investigations later. As such, we take a snapshot of a PNS at 2.3 seconds after bounce from a 1D PNS cooling calculation and run our simulation for $\sim160\,\mathrm{ms}$ in 2D under axisymmetry. The original PNS cooling calculation neglected convection entirely and the initial condition is linearly unstable to convection. We find in our 2D simulation that convection is instigated there indeed and expands inward after being full-fledged. The convection is then settled to a quasi-steady state in $\sim100\,\mathrm{ms}$, being sustained by the negative $Y_e$ gradient, which is in turn maintained by neutrino emissions. It enhances the luminosities and mean energies of all species of neutrinos compared to 1D. Taking advantage of the Boltzmann solver, we analyze the possible occurrence of the neutrino fast flavor conversion (FFC). We found that FFC is likely to occur in the regions, where $Y_e$ is lower, and that the growth rate can be as high as $\sim 10^{-1}\,{\mathrm{cm}^{-1}}$. △ Less

Submitted 20 December, 2022; v1 submitted 3 June, 2022; originally announced June 2022.

Comments: accepted to ApJ

Report number: RIKEN-iTHEMS-Report-22

Journal ref: The Astrophysical Journal, 944:60 (18pp), 2023 February 10

Detection of highly correlated optical and X-ray variations in SS Cygni with Tomo-e Gozen and NICER

Authors: Y. Nishino, M. Kimura, S. Sako, J. Beniyama, T. Enoto, T. Minezaki, N. Nakaniwa, R. Ohsawa, S. Takita, S. Yamada, K. C. Gendreau

Abstract: We report on simultaneous optical and X-ray observations of the dwarf nova SS Cyg with Tomo-e Gozen/1.05 m Kiso Schmidt and Neutron star Interior Composition ExploreR (it NICER) / International Space Station (ISS). A total of four observations were carried out in the quiescent state and highly correlated light variations between the two wavelengths were detected. We have extracted local brightness… ▽ More We report on simultaneous optical and X-ray observations of the dwarf nova SS Cyg with Tomo-e Gozen/1.05 m Kiso Schmidt and Neutron star Interior Composition ExploreR (it NICER) / International Space Station (ISS). A total of four observations were carried out in the quiescent state and highly correlated light variations between the two wavelengths were detected. We have extracted local brightness peaks in the light curves with a binning interval of 1 sec, called `shots', and have evaluated time lags between the optical and X-ray variations by using a cross-correlation function. Some shots exhibit significant optical lags to X-ray variations and most of them are positive ranging from $+$0.26 to 3.11 sec, which have never been detected. They may be ascribable to X-ray reprocessing in the accretion disk and/or the secondary star. This analysis is possible thanks to the high timing accuracy and the high throughput of NICER and the matching capability of Tomo-e Gozen. Also, it is confirmed that the correlation between the optical and X-ray variations changed in the middle of one of our observation windows and the X-ray spectrum softer than 2 keV varied in accordance with the correlation. △ Less

Submitted 18 May, 2022; originally announced May 2022.

Comments: 6 pages, 4 figures, 2 tables. Published by Publication of the Astronomical Society of Japan. Supplementary information is available online

X-ray Studies of the Inverted Ejecta Layers in the Southeast Area of Cassiopeia A

Authors: Tomoya Tsuchioka, Toshiki Sato, Shinya Yamada, Yasunobu Uchiyama

Abstract: The central strong activities in core-collapse supernovae expect to produce the overturning of the Fe- and Si/O-rich ejecta during the supernova explosion based on multi-dimensional simulations. X-ray observations of the supernova remnant Cassiopeia A have indicated that the Fe-rich ejecta lies outside the Si-rich materials in the southeastern region, which is consistent with the hypothesis on the… ▽ More The central strong activities in core-collapse supernovae expect to produce the overturning of the Fe- and Si/O-rich ejecta during the supernova explosion based on multi-dimensional simulations. X-ray observations of the supernova remnant Cassiopeia A have indicated that the Fe-rich ejecta lies outside the Si-rich materials in the southeastern region, which is consistent with the hypothesis on the inversion of the ejecta. We investigate the kinematic and nucleosynthetic properties of the inverted ejecta layers in detail to understand its formation process using the data taken by the Chandra X-ray Observatory. Three-dimensional velocities of Fe- and Si/O-rich ejecta are obtained as $>$4,500 km s$^{-1}$ and $\sim$2,000--3,000 km s$^{-1}$, respectively, by combining proper motion and line-of-sight velocities, indicating that the velocity of the Si/O-rich ejecta is slower than that of the Fe-rich ejecta since the early stage of the explosion. To constrain their burning regime, the Cr/Fe mass ratios are evaluated as $0.51_{-0.10}^{+0.11}$% in the outermost Fe-rich region and $1.24 ^{+0.19}_{-0.20}$% in the inner Fe/Si-rich region, suggesting that the complete Si burning layer is invertedly located to the incomplete Si burning layer. All the results support the ejecta overturning at the early stages of the remnant's evolution or during the supernova explosion of Cassiopeia A. △ Less

Submitted 10 May, 2022; originally announced May 2022.

Comments: ApJ accepted. 11 pages, 7 figures

Application of deep learning to the evaluation of goodness in the waveform processing of transition-edge sensor calorimeters

Authors: Y. Ichinohe, S. Yamada, R. Hayakawa, S. Okada, T. Hashimoto, H. Tatsuno, H. Suda, T. Okumura

Abstract: Optimal filtering is the crucial technique for the data analysis of transition-edge-sensor (TES) calorimeters to achieve their state-of-the-art energy resolutions. Filtering out the `bad' data from the dataset is important because it otherwise leads to the degradation of energy resolutions, while it is not a trivial task. We propose a neural network-based technique for the automatic goodness taggi… ▽ More Optimal filtering is the crucial technique for the data analysis of transition-edge-sensor (TES) calorimeters to achieve their state-of-the-art energy resolutions. Filtering out the `bad' data from the dataset is important because it otherwise leads to the degradation of energy resolutions, while it is not a trivial task. We propose a neural network-based technique for the automatic goodness tagging of TES pulses, which is fast and automatic and does not require bad data for training. △ Less

Submitted 25 April, 2022; originally announced April 2022.

Comments: Proceedings of the LTD19 conference

A novel formulation for the evolution of relativistic rotating stars

Authors: Hirotada Okawa, Kotaro Fujisawa, Nobutoshi Yasutake, Misa Ogata, Yu Yamamoto, Shoichi Yamada

Abstract: We present a new formulation to construct numerically equilibrium configurations of rotating stars in general relativity. Having in mind the application to their quasi static evolutions, we adopt a Lagrangian formulation of our own devising, in which we solve force balance equations to seek for the positions of fluid elements assigned to the grid points, instead of the ordinary Eulerian formulatio… ▽ More We present a new formulation to construct numerically equilibrium configurations of rotating stars in general relativity. Having in mind the application to their quasi static evolutions, we adopt a Lagrangian formulation of our own devising, in which we solve force balance equations to seek for the positions of fluid elements assigned to the grid points, instead of the ordinary Eulerian formulation. Unlike previous works in the literature, we do not employ the first integral of the Euler equation, which is not obtained by an analytic integration in general. We assign a mass, specific angular momentum and entropy to each fluid element in contrast to the previous methods, in which the spatial distribution of the angular velocity or angular momentum is specified. Those distributions are determined after the positions of all fluid elements (or grid points) are derived in our formulation. We solve the large system of algebraic nonlinear equations that are obtained by discretizing the time-independent Euler and Einstein equations in the finite-elements method by using our new multi-dimensional root-finding scheme, named the W4 method. To demonstrate the capability of our new formulation, we construct some rotational configurations both barotropic and baroclinic. We also solve three evolutionary sequences that mimic the cooling, mass-loss, and mass-accretion as simple toy models. △ Less

Submitted 21 April, 2022; originally announced April 2022.

Comments: 19 pages, 13 figures

Singularity-Avoiding Multi-Dimensional Root-Finder

Authors: Hirotada Okawa, Kotaro Fujisawa, Yu Yamamoto, Nobutoshi Yasutake, Misa Ogata, Shoichi Yamada

Abstract: We proposed in this paper a new method, which we named the W4 method, to solve nonlinear equation systems. It may be regarded as an extension of the Newton-Raphson~(NR) method to be used when the method fails. Indeed our method can be applied not only to ordinary problems with non-singular Jacobian matrices but also to problems with singular Jacobians, which essentially all previous methods that e… ▽ More We proposed in this paper a new method, which we named the W4 method, to solve nonlinear equation systems. It may be regarded as an extension of the Newton-Raphson~(NR) method to be used when the method fails. Indeed our method can be applied not only to ordinary problems with non-singular Jacobian matrices but also to problems with singular Jacobians, which essentially all previous methods that employ the inversion of the Jacobian matrix have failed to solve. In this article, we demonstrate that (i) our new scheme can define a non-singular iteration map even for those problems by utilizing the singular value decomposition, (ii) a series of vectors in the new iteration map converges to the right solution under a certain condition, (iii) the standard two-dimensional problems in the literature that no single method proposed so far has been able to solve completely are all solved by our new method. △ Less

Submitted 21 April, 2022; originally announced April 2022.

Comments: 13 pages, 4 figures

NuSTAR Observations of 52 Compton-thick Active Galactic Nuclei Selected by the Swift/BAT All-sky Hard X-Ray Survey

Authors: Atsushi Tanimoto, Yoshihiro Ueda, Hirokazu Odaka, Satoshi Yamada, Claudio Ricci

Abstract: We present the systematic broadband X-ray spectral analysis of 52 Compton-thick ($24 \leq \log N_{\mathrm{H}}^{\mathrm{LOS}}/\mathrm{cm}^{-2}$) active galactic nucleus (CTAGN) candidates selected by the Swift/BAT all-sky hard X-ray survey observed with Chandra, XMM--Newton, Swift/XRT, Suzaku, and NuSTAR. The XMM--Newton data of 10 objects and the NuSTAR data of 15 objects are published for the fir… ▽ More We present the systematic broadband X-ray spectral analysis of 52 Compton-thick ($24 \leq \log N_{\mathrm{H}}^{\mathrm{LOS}}/\mathrm{cm}^{-2}$) active galactic nucleus (CTAGN) candidates selected by the Swift/BAT all-sky hard X-ray survey observed with Chandra, XMM--Newton, Swift/XRT, Suzaku, and NuSTAR. The XMM--Newton data of 10 objects and the NuSTAR data of 15 objects are published for the first time. We use an X-ray spectral model from a clumpy torus (XClumpy) to determine the torus properties. As a result, the hydrogen column density along the line of sight $N_{\mathrm{H}}^{\mathrm{LOS}}$ obtained from the XClumpy model indicate that 24 objects are Compton-thin AGNs and 28 objects are Compton-thick AGNs in the 90\% confidence interval. The main reason is the difference in the torus model applied. The hydrogen column density along the equatorial direction $N_{\mathrm{H}}^{\mathrm{Equ}}$ of CTAGNs inferred from the XClumpy model is larger than that of less obscured AGNs. The Compton-thin torus covering factor $C_{22}$ obtained from the XClumpy model is consistent with that of Ricci et al. (2017) in the low Eddington ratio ($\log R_{\mathrm{Edd}} \leq -1.0$), whereas $C_{22}$ inferred from the XClumpy model is larger than that of Ricci et al. (2017) in the high Eddington ratio ($-1.0 \leq \log R_{\mathrm{Edd}}$). The average value of the Compton-thick torus covering factor $C_{24}$ obtained from the XClumpy model is $36_{-4}^{+4}$\%. This value is larger than that of Ricci et al. (2015) ($C_{24} \simeq 27_{-4}^{+4}$\%) based on the assumption that all AGNs have intrinsically the same torus structure. These results suggest that the structure of CTAGN may be intrinsically different from that of less obscured AGN. △ Less

Submitted 24 March, 2022; originally announced March 2022.

Comments: 44 pages, 16 figures, accepted for publication in the Astrophysical Journal Supplement Series

Gravitational wave signals in the deci-Hz range from neutrinos during the proto-neutron star cooling phase

Authors: Lei Fu, Shoichi Yamada

Abstract: We investigate the gravitational waves (GWs) at low frequencies produced by neutrinos that are emitted anisotropically from the proto-neutron star (PNS) during its cooling phase that lasts about a minute. We first give a formulation based on the spherical-harmonic expansion of the neutrino luminosity to obtain the gravitational waveform as well as the characteristic strain. In the absence of multi… ▽ More We investigate the gravitational waves (GWs) at low frequencies produced by neutrinos that are emitted anisotropically from the proto-neutron star (PNS) during its cooling phase that lasts about a minute. We first give a formulation based on the spherical-harmonic expansion of the neutrino luminosity to obtain the gravitational waveform as well as the characteristic strain. In the absence of multi-dimensional simulations of PNS cooling, from which we can extract reliable data on the neutrino luminosities as a function of solid angle, we construct them by hand. In the first model, the time evolution is approximated by piece-wise exponential functions (PEFs); in the second model we employ the time evolution obtained in a 1D cooling simulation for all harmonic components for simplicity; In both cases, we consider not only axisymmetric components but also non-axisymmetric ones; for the latter, in particular, we consider as the third model axisymetric neutrino emissions that are misaligned with the rotation axis and, as a result, rotate with the PNS. We find from the first model that the decay times in PEF at late phases can be inferred from the positions of bumps and dips in the characteristic strain of the GW in the case of a slow cooling, whereas it may be obtained by identifying the positions of slope change in the case of rapid cooling, which may be induced by convections in PNS. We also find from the first model that the rotation of PNS can be reflected in GW waveforms. Finally, we compared the characteristic strains with sensitivity curves GW detectors, finding that above phenomena are expected to be detected by DECIGO if the PNS is located at 10kpc. △ Less

Submitted 8 June, 2022; v1 submitted 30 January, 2022; originally announced January 2022.

Discovery of Stable Titanium at the Northeastern Jet of Cassiopeia A: Need for a Weak Jet Mechanism?

Authors: Takuma Ikeda, Yasunobu Uchiyama, Toshiki Sato, Ryota Higurashi, Tomoya Tsuchioka, Shinya Yamada

Abstract: The origin of the jet-like structures observed in Cassiopeia A is still unclear, although it seems to be related to its explosion mechanism. X-ray observations of the characteristic structures could provide us useful information on the explosive nucleosynthesis via the observation of elements, which is a unique approach to understand its origin. We here report the discovery of shocked stable Ti, w… ▽ More The origin of the jet-like structures observed in Cassiopeia A is still unclear, although it seems to be related to its explosion mechanism. X-ray observations of the characteristic structures could provide us useful information on the explosive nucleosynthesis via the observation of elements, which is a unique approach to understand its origin. We here report the discovery of shocked stable Ti, which is produced only at the inner region of exploding stars, in the northeast jet of Cassiopeia A using the 1-Ms deep observation with the Chandra X-ray observatory. The observed Ti coexists with other intermediate-mass elements (e.g. Si, S, Ar, Ca) and Fe at the tip of the X-ray jet structure. We found that its elemental composition is well explained with the production by the incomplete Si burning regime, indicating that the formation process of the jet structure was sub-energetic at the explosion (the peak temperature during the nuclear burning must be $\lesssim$ 5$\times$10$^{9}$ K at most). Thus, we conclude that the energy source that formed the jet structure was not the primary engine for the supernova explosion. Our results are useful to limit the power of the jet-structure formation process, and a weak jet mechanism with low temperature may be needed to explain it. △ Less

Submitted 10 January, 2022; originally announced January 2022.

Comments: Accepted for publication in PASJ, 9 pages, 5 figures, 2 tables

Quantum synchro-curvature masers and their application to neutron stars

Authors: Hiroko Tomoda, Tomoya Naoe, Shoichi Yamada

Abstract: We explore the possibility of synchro-curvature maser in the magnetosphere of neutron stars (NSs). Unlike previous studies, we employ relativistic quantum mechanics, solving the Dirac equation for an electron in helical magnetic fields and calculating the radiative transition rates perturbatively. Assuming that the curvature of magnetic-field lines is much larger than the Larmor radius, we utilize… ▽ More We explore the possibility of synchro-curvature maser in the magnetosphere of neutron stars (NSs). Unlike previous studies, we employ relativistic quantum mechanics, solving the Dirac equation for an electron in helical magnetic fields and calculating the radiative transition rates perturbatively. Assuming that the curvature of magnetic-field lines is much larger than the Larmor radius, we utilize adiabatic spinor rotations to obtain the wave functions of an electron. We classify the electron states further either by the spin operator projected on the magnetic field or by the helicity operator. We then evaluate numerically the true absorption rates accounting for the induced emission for some parameter values typical to the outer gaps of different types of NSs. We show that maser is indeed possible for a range of parameters. We will also present the dependence on those parameters systematically. We finally give a crude estimate of the amplification factor in the outer gap of NSs, which seems to favor millisecond pulsars as the host of maser emissions. △ Less

Submitted 21 November, 2021; originally announced November 2021.

Comments: 25 pages, 11 figures

Principal-Axis Analysis of the Eddington Tensor for the Early Post-Bounce Phase of Rotational Core-Collapse Supernovae

Authors: Wakana Iwakami, Akira Harada, Hiroki Nagakura, Ryuichiro Akaho, Hirotada Okawa, Shun Furusawa, Hideo Matsufuru, Kohsuke Sumiyoshi, Shoichi Yamada

Abstract: Using full Boltzmann neutrino transport, we performed two-dimensional (2D) core-collapse supernova simulations in axisymmetry for two progenitor models with 11.2M and 15.0M both rotational and non-rotational. We employed the results obtained in the early post-bounce phase (t < 20 ms) to assess performance under rapid rotation of some closure relations commonly employed in the truncated moment meth… ▽ More Using full Boltzmann neutrino transport, we performed two-dimensional (2D) core-collapse supernova simulations in axisymmetry for two progenitor models with 11.2M and 15.0M both rotational and non-rotational. We employed the results obtained in the early post-bounce phase (t < 20 ms) to assess performance under rapid rotation of some closure relations commonly employed in the truncated moment method. We first made a comparison in 1D under spherical symmetry, though, of the Eddington factor p defined in the fluid rest frame (FR). We confirmed that the maximum entropy closure for the Fermionic distribution (MEFD) performs better than others near the proto-neutron star surface, where p < 1/3 occurs, but does not work well even in 1D when the phase space occupancy satisfies e < 0.5 together with p < 1/3, the condition known to be not represented by MEFD. For the 2D models with the rapid rotation, we employed the principal-axis analysis of the Eddington tensor. We paid particular attention to the direction of the longest principal axis. We observed in FR that it is aligned neither with the radial direction nor with the neutrino flux in 2D, particularly so in convective and/or rapidly rotating regions, the fact not accommodated in the moment method. We repeated the same analysis in the laboratory frame (LB) and found again that the direction of the longest principal axis is not well reproduced by MEFD because the interpolation between the optically thick and thin limits is not very accurate in this frame. △ Less

Submitted 21 May, 2022; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: 28pages, 22figures, accepted for publication in the ApJ

Report number: RIKEN-iTHEMS-Report-21

A hard X-ray view of Luminous and Ultra-luminous Infrared Galaxies in GOALS: I - AGN obscuration along the merger sequence

Authors: C. Ricci, G. C. Privon, R. W. Pfeifle, L. Armus, K. Iwasawa, N. Torres-Alba, S. Satyapal, F. E. Bauer, E. Treister, L. C. Ho, S. Aalto, P. Arevalo, L. Barcos-Munoz, V. Charmandaris, T. Diaz-Santos, A. S. Evans, T. Gao, H. Inami, M. J. Koss, G. Lansbury, S. T. Linden, A. Medling, D. B. Sanders, Y. Song, D. Stern , et al. (3 additional authors not shown)

Abstract: The merger of two or more galaxies can enhance the inflow of material from galactic scales into the close environments of Active Galactic Nuclei (AGN), obscuring and feeding the supermassive black hole (SMBH). Both recent simulations and observations of AGN in mergers have confirmed that mergers are related to strong nuclear obscuration. However, it is still unclear how AGN obscuration evolves in… ▽ More The merger of two or more galaxies can enhance the inflow of material from galactic scales into the close environments of Active Galactic Nuclei (AGN), obscuring and feeding the supermassive black hole (SMBH). Both recent simulations and observations of AGN in mergers have confirmed that mergers are related to strong nuclear obscuration. However, it is still unclear how AGN obscuration evolves in the last phases of the merger process. We study a sample of 60 Luminous and Ultra-luminous IR galaxies (U/LIRGs) from the GOALS sample observed by NuSTAR. We find that the fraction of AGN that are Compton-thick (CT; $N_{\rm H}\geq 10^{24}\rm\,cm^{-2}$) peaks at $74_{-19}^{+14}\%$ at a late merger stage, prior to coalescence, when the nuclei have projected separations of $d_{\rm sep}\sim 0.4-6$ kpc. A similar peak is also observed in the median $N_{\rm H}$ [$(1.6\pm0.5)\times10^{24}\rm\,cm^{-2}$]. The vast majority ($85^{+7}_{-9}\%$) of the AGN in the final merger stages ($d_{\rm sep}\lesssim 10$ kpc) are heavily obscured ($N_{\rm H}\geq 10^{23}\rm\,cm^{-2}$), and the median $N_{\rm H}$ of the accreting SMBHs in our sample is systematically higher than that of local hard X-ray selected AGN, regardless of the merger stage. This implies that these objects have very obscured nuclear environments, with the $N_{\rm H}\geq 10^{23}\rm\,cm^{-2}$ gas almost completely covering the AGN in late mergers. CT AGN tend to have systematically higher absorption-corrected X-ray luminosities than less obscured sources. This could either be due to an evolutionary effect, with more obscured sources accreting more rapidly because they have more gas available in their surroundings, or to a selection bias. The latter scenario would imply that we are still missing a large fraction of heavily obscured, lower luminosity ($L_{2-10}\lesssim 10^{43}\rm\,erg\,s^{-1}$) AGN in U/LIRGs. △ Less

Submitted 9 August, 2021; v1 submitted 22 July, 2021; originally announced July 2021.

Comprehensive Broadband X-ray and Multiwavelength Study of Active Galactic Nuclei in Local 57 Ultra/luminous Infrared Galaxies Observed with NuSTAR and/or Swift/BAT

Authors: Satoshi Yamada, Yoshihiro Ueda, Atsushi Tanimoto, Masatoshi Imanishi, Yoshiki Toba, Claudio Ricci, George C. Privon

Abstract: We perform a systematic X-ray spectroscopic analysis of 57 local ultra/luminous infrared galaxy systems (containing 84 individual galaxies) observed with Nuclear Spectroscopic Telescope Array and/or Swift/BAT. Combining soft X-ray data obtained with Chandra, XMM-Newton, Suzaku and/or Swift/XRT, we identify 40 hard ($>$10 keV) X-ray detected active galactic nuclei (AGNs) and constrain their torus p… ▽ More We perform a systematic X-ray spectroscopic analysis of 57 local ultra/luminous infrared galaxy systems (containing 84 individual galaxies) observed with Nuclear Spectroscopic Telescope Array and/or Swift/BAT. Combining soft X-ray data obtained with Chandra, XMM-Newton, Suzaku and/or Swift/XRT, we identify 40 hard ($>$10 keV) X-ray detected active galactic nuclei (AGNs) and constrain their torus parameters with the X-ray clumpy torus model XCLUMPY (Tanimoto et al. 2019). Among the AGNs at $z < 0.03$, for which sample biases are minimized, the fraction of Compton-thick ($N_{\rm H} \geq 10^{24}$ cm$^{-2}$) AGNs reaches 64$^{+14}_{-15}$% (6/9 sources) in late mergers, while 24$^{+12}_{-10}$% (3/14 sources) in early mergers, consistent with the tendency reported by Ricci et al. (2017). We find that the bolometric AGN luminosities derived from the infrared data increase, but the X-ray to bolometric luminosity ratios decrease, with merger stage. The X-ray weak AGNs in late mergers ubiquitously show massive outflows at sub-pc to kpc scales. Among them, the most luminous AGNs ($L_{\rm bol,AGN} \sim 10^{46}$ erg s$^{-1}$) have relatively small column densities of $\lesssim$10$^{23}$ cm$^{-2}$ and almost super-Eddington ratios ($λ_{\rm Edd} \sim$ 1.0). Their torus covering factors ($C_{\rm T}^{\rm (22)} \sim 0.6$) are larger than those of Swift/BAT selected AGNs with similarly high Eddington ratios. These results suggest a scenario that, in the final stage of mergers, multiphase strong outflows are produced due to chaotic quasi-spherical inflows and the AGN becomes extremely X-ray weak and deeply buried due to obscuration by inflowing and/or outflowing material. △ Less

Submitted 22 July, 2021; originally announced July 2021.

Comments: 85 pages, 88 figures, accepted in ApJS

On the nature of the anomalous event in 2021 in the dwarf nova SS Cygni and its multi-wavelength transition

Authors: Mariko Kimura, Shinya Yamada, Nozomi Nakaniwa, Yoshihiro Makita, Hitoshi Negoro, Megumi Shidatsu, Taichi Kato, Teruaki Enoto, Keisuke Isogai, Tatehiro Mihara, Hidehiko Akazawa, Keith C. Gendreau, Franz-Josef Hambsch, Pavol A. Dubovsky, Igor Kudzej, Kiyoshi Kasai, Tamás Tordai, Elena Pavlenko, Aleksei A. Sosnovskij, Julia V. Babina, Oksana I. Antonyuk, Hiroshi Itoh, Hiroyuki Maehara

Abstract: SS Cyg has long been recognized as the prototype of a group of dwarf novae that show only outbursts. However, this object has entered a quite anomalous event in 2021, which at first appeared to be standstill, i.e., an almost constant luminosity state, observed in Z Cam-type dwarf novae. This unexpected event gives us a great opportunity to reconsider the nature of standstill in cataclysmic variabl… ▽ More SS Cyg has long been recognized as the prototype of a group of dwarf novae that show only outbursts. However, this object has entered a quite anomalous event in 2021, which at first appeared to be standstill, i.e., an almost constant luminosity state, observed in Z Cam-type dwarf novae. This unexpected event gives us a great opportunity to reconsider the nature of standstill in cataclysmic variables. We have observed this anomalous event and its forerunner, a gradual and simultaneous increase in the optical and X-ray flux during quiescence, through many optical telescopes and the X-ray telescopes NICER and NuSTAR. We have not found any amplification of the orbital hump during quiescence before the anomalous event, which suggests that the mass transfer rate did not significantly fluctuate on average. The estimated X-ray flux was not enough to explain the increment of the optical flux during quiescence via X-ray irradiation of the disk and the secondary star. It would be natural to consider that viscosity in the quiescent disk was enhanced before the anomalous event, which increased mass accretion rates in the disk and raised not only the optical flux but also the X-ray flux. We suggest that enhanced viscosity also triggered the standstill-like phenomenon in SS Cyg, which is considered to be a series of small outbursts. The inner part of the disk would always stay in the outburst state and only its outer part would be unstable against the thermal-viscous instability during this phenomenon, which is consistent with the observed optical color variations. This scenario is in line with our X-ray spectral analyses which imply that the X-ray emitting inner accretion flow became hotter than usual and vertically expanded and that it became denser and was cooled down after the onset of the standstill-like state. △ Less

Submitted 29 June, 2021; originally announced June 2021.

Comments: Accepted for publication in PASJ. 18 pages, 10 figures, and 2 tables

Deep Learning of the Eddington Tensor in the Core-collapse Supernova Simulation

Authors: Akira Harada, Shota Nishikawa, Shoichi Yamada

Abstract: We trained deep neural networks (DNNs) as a function of the neutrino energy density, flux, and the fluid velocity to reproduce the Eddington tensor for neutrinos obtained in our first-principles core-collapse supernova (CCSN) simulations. Although the moment method, which is one of the most popular approximations for neutrino transport, requires a closure relation, none of the analytical closure r… ▽ More We trained deep neural networks (DNNs) as a function of the neutrino energy density, flux, and the fluid velocity to reproduce the Eddington tensor for neutrinos obtained in our first-principles core-collapse supernova (CCSN) simulations. Although the moment method, which is one of the most popular approximations for neutrino transport, requires a closure relation, none of the analytical closure relations commonly employed in the literature captures all aspects of the neutrino angular distribution in momentum space. In this paper, we developed a closure relation by using the DNN that takes the neutrino energy density, flux, and the fluid velocity as the input and the Eddington tensor as the output. We consider two kinds of DNNs: a conventional DNN named a component-wise neural network (CWNN) and a tensor-basis neural network (TBNN). We found that the diagonal component of the Eddington tensor is reproduced better by the DNNs than the M1-closure relation especially for low to intermediate energies. For the off-diagonal component, the DNNs agree better with the Boltzmann solver than the M1 closure at large radii. In the comparison between the two DNNs, the TBNN has slightly better performance than the CWNN. With the new closure relations at hand based on the DNNs that well reproduce the Eddington tensor with much smaller costs, we opened up a new possibility for the moment method. △ Less

Submitted 15 November, 2021; v1 submitted 27 April, 2021; originally announced April 2021.

Comments: 15 pages, 13 figures, accepted for publication in the ApJ

Report number: RIKEN-iTHEMS-Report-21

On the origin of the asymmetry of the ejecta structure and explosion of G350.1$-$0.3

Authors: Tomoya Tsuchioka, Yasunobu Uchiyama, Ryota Higurashi, Hiroyoshi Iwasaki, Shumpei Otsuka, Shinya Yamada, Toshiki Sato

Abstract: We present X-ray analysis of the ejecta of supernova remnant G350.1$-$0.3 observed with Chandra and Suzaku, and clarify the ejecta's kinematics over a decade and obtain a new observational clue to understanding the origin of the asymmetric explosion. Two images of Chandra X-ray Observatory taken in 2009 and 2018 are analyzed in several methods, and enable us to measure the velocities in the plane… ▽ More We present X-ray analysis of the ejecta of supernova remnant G350.1$-$0.3 observed with Chandra and Suzaku, and clarify the ejecta's kinematics over a decade and obtain a new observational clue to understanding the origin of the asymmetric explosion. Two images of Chandra X-ray Observatory taken in 2009 and 2018 are analyzed in several methods, and enable us to measure the velocities in the plane of the sky. A maximum velocity is 4640$\pm$290 km s$^{-1}$ (0.218$\pm$0.014 arcsec yr$^{-1}$) in the eastern region in the remnant. These findings trigger us to scrutinize the Doppler effects in the spectra of the thermal emission, and the velocities in the line-of-sight direction are estimated to be a thousand km s$^{-1}$. The results are confirmed by analyzing the spectra of Suzaku. Combining the proper motions and line-of-sight velocities, the ejecta's three-dimensional velocities are $\sim$3000-5000 km s$^{-1}$. The center of the explosion is more stringently constrained by finding the optimal time to reproduce the observed spatial expansion. Our findings that the age of the SNR is estimated at most to be 655 years, and the CCO is observed as a point source object against the SNR strengthen the 'hydrodynamical kick' hypothesis on the origin of the remnant. △ Less

Submitted 7 April, 2021; originally announced April 2021.

Comments: ApJ accepted. 12 pages, 6 figures

X-ray Constraint on Location of AGN Torus in Circinus Galaxy

Authors: Ryosuke Uematsu, Yoshihiro Ueda, Atsushi Tanimoto, Taiki Kawamuro, Kenta Setoguchi, Shoji Ogawa, Satoshi Yamada, Hirokazu Odaka

Abstract: The location of the obscuring "torus" in an active galactic nucleus (AGN) is still an unresolved issue. The line widths of X-ray fluorescence lines originated from the torus, particularly Fe K$α$, carry key information on the radii of line emitting regions. Utilizing XCLUMPY (Tanimoto et al. 2019), an X-ray clumpy torus model, we develop a realistic model of emission line profiles from an AGN toru… ▽ More The location of the obscuring "torus" in an active galactic nucleus (AGN) is still an unresolved issue. The line widths of X-ray fluorescence lines originated from the torus, particularly Fe K$α$, carry key information on the radii of line emitting regions. Utilizing XCLUMPY (Tanimoto et al. 2019), an X-ray clumpy torus model, we develop a realistic model of emission line profiles from an AGN torus where we take into account line broadening due to the Keplerian motion around the black hole. Then, we apply the updated model to the best available broadband spectra (3-100 keV) of the Circinus galaxy observed with Suzaku, XMM-Newton, Nuclear Spectroscopic Telescope Array (NuSTAR), and Chandra, including 0.62 Ms Chandra/HETG data. We confirm that the torus is Compton-thick (hydrogen column-density along the equatorial plane is $N_\mathrm{H}^\mathrm{Equ}=2.16^{+0.24}_{-0.16}\times 10^{25}\ \mathrm{cm}^{-2}$), geometrically thin (torus angular width $σ=10.3^{+0.7}_{-0.3}\ \mathrm{degrees}$), viewed edge-on (inclination $i=78.3^{+0.4}_{-0.9}\ \mathrm{degrees}$), and has super-solar abundance ($1.52^{+0.04}_{-0.06}$ times solar). Simultaneously analyzing the Chandra/HETG first, second, and third order spectra with consideration of the spatial extent of the Fe K$α$ line emitting region, we constrain the inner radius of the torus to be $1.9^{+3.1}_{-0.8}\times 10^5$ times the gravitational radius, or $1.6^{+1.5}_{-0.9}\times 10^{-2}\ \mathrm{pc}$ for a black hole mass of $(1.7\pm 0.3)\times 10^6\ M_{\odot}$. This is about 3 times smaller than that estimated from the dust sublimation radius, suggesting that the inner side of the dusty region of the torus is composed of dust-free gas. △ Less

Submitted 20 March, 2021; originally announced March 2021.

Comments: 20 pages, 10 figures, 8 tables, accepted for publication in ApJ

Systematic Study of AGN Clumpy Tori with Broadband X-ray Spectroscopy: Updated Unified Picture of AGN Structure

Authors: Shoji Ogawa, Yoshihiro Ueda, Atsushi Tanimoto, Satoshi Yamada

Abstract: We present the results of a systematic, broadband X-ray spectral analysis of nearby active galactic nuclei (AGNs) with the X-ray clumpy torus model (XCLUMPY; Tanimoto et al. 2019). By adding 16 AGNs newly analyzed in this paper, we study total 28 AGNs including unabsorbed and absorbed AGNs taken from Ichikawa et al. (2015) and García-Bernete et al. (2019). This is the largest sample whose X-ray an… ▽ More We present the results of a systematic, broadband X-ray spectral analysis of nearby active galactic nuclei (AGNs) with the X-ray clumpy torus model (XCLUMPY; Tanimoto et al. 2019). By adding 16 AGNs newly analyzed in this paper, we study total 28 AGNs including unabsorbed and absorbed AGNs taken from Ichikawa et al. (2015) and García-Bernete et al. (2019). This is the largest sample whose X-ray and infrared spectra are analyzed by the clumpy torus models XCLUMPY and CLUMPY (Nenkova et al. 2008), respectively. The relation between the Eddington ratio and the torus covering factor determined from the X-ray torus parameters of each object follows the trend found by Ricci et al. (2017) based on a statistical analysis. We confirm the results by Tanimoto et al. (2020) that (1) the torus angular widths determined by the infrared data are larger than those by the X-ray data and that (2) the ratios of the hydrogen column density to V-band extinction ($N_{\rm H}/A_{\rm V}$) along the line of sight in obscured AGNs are similar to the Galactic value on average. Unobscured AGNs show apparently smaller line-of-sight $N_{\rm H}/A_{\rm V}$ ratios than the Galactic one. Our findings can be well explained by an updated unified picture of AGN structure including a dusty torus, dusty polar outflows, and dust-free gas, where the inclination determines the X-ray and optical classifications and observed torus properties in the X-ray and infrared bands. △ Less

Submitted 27 November, 2020; originally announced November 2020.

Comments: 24 pages, 6 figures, and 5 tables, accepted for publication in ApJ

Multidimensional Boltzmann Neutrino Transport Code in Full General Relativity for Core-collapse Simulations

Authors: Ryuichiro Akaho, Akira Harada, Hiroki Nagakura, Kohsuke Sumiyoshi, Wakana Iwakami, Hirotada Okawa, Shun Furusawa, Hideo Matsufuru, Shoichi Yamada

Abstract: We develop a neutrino transfer code for core-collapse simulations, that directly solves the multidimensional Boltzmann equations in full general relativity. We employ the discrete ordinate method, which discretizes the six dimensional phase space. The code is an extension of our special relativistic code coupled to a Newtonian hydrodynamics code, which is currently employed for core-collapse super… ▽ More We develop a neutrino transfer code for core-collapse simulations, that directly solves the multidimensional Boltzmann equations in full general relativity. We employ the discrete ordinate method, which discretizes the six dimensional phase space. The code is an extension of our special relativistic code coupled to a Newtonian hydrodynamics code, which is currently employed for core-collapse supernova simulations. In order to demonstrate our code's capability to treat general relativistic effects, we conduct some tests: we first compute the free streaming of neutrinos in the Schwarzschild and Kerr spacetimes and compare the results with the geodesic curves; in the Schwarzschild case we deploy not only a 1-dimensional grid in space under spherical symmetry but also a 2-dimensional spatial mesh under axisymmetry in order to assess the capability of the code to compute the spatial advection of neutrinos; secondly, we calculate the neutrino transport in a fixed matter background, which is taken from a core-collapse supernova simulation with our general relativistic but spherically symmetric Boltzmann-hydrodynamics code, to obtain a steady neutrino distribution; the results are compared with those given by the latter code. △ Less

Submitted 19 March, 2021; v1 submitted 21 October, 2020; originally announced October 2020.

Comments: 23 pages, 19 figures, accepted version

Journal ref: The Astrophysical Journal, 909:210 (17pp), 2021 March 10

Application of X-Ray Clumpy Torus Model (XCLUMPY) to 10 Obscured Active Galactic Nuclei Observed with Suzaku and NuSTAR

Authors: Atsushi Tanimoto, Yoshihiro Ueda, Hirokazu Odaka, Shoji Ogawa, Satoshi Yamada, Toshihiro Kawaguchi, Kohei Ichikawa

Abstract: We apply XCLUMPY, an X-ray spectral model from a clumpy torus in an active galactic nucleus (AGN), to the broadband X-ray spectra of 10 obscured AGNs observed with both Suzaku and NuSTAR. The infrared spectra of these AGNs were analyzed with the CLUMPY code. Since XCLUMPY adopts the same clump distribution as that in the CLUMPY, we can directly compare the torus parameters obtained from the X-ray… ▽ More We apply XCLUMPY, an X-ray spectral model from a clumpy torus in an active galactic nucleus (AGN), to the broadband X-ray spectra of 10 obscured AGNs observed with both Suzaku and NuSTAR. The infrared spectra of these AGNs were analyzed with the CLUMPY code. Since XCLUMPY adopts the same clump distribution as that in the CLUMPY, we can directly compare the torus parameters obtained from the X-ray spectra and those from the infrared ones. The torus angular widths determined from the infrared spectra ($σ_{\mathrm{IR}}$) are systematically larger than those from the X-ray data ($σ_{\mathrm{X}}$); the difference ($σ_{\mathrm{IR}}-σ_{\mathrm{X}}$) correlates with the inclination angle determined from the X-ray spectrum. These results can be explained by the contribution from dusty polar outflows to the observed infrared flux, which becomes more significant at higher inclinations (more edge-on views). The ratio of the hydrogen column density and V-band extinction in the line of sight absorber shows large scatter ($\simeq$1 dex) around the Galactic value, suggesting that a significant fraction of AGNs have dust-rich circumnuclear environments. △ Less

Submitted 26 May, 2020; originally announced May 2020.

Comments: 17 pages, 3 figures, accepted for publication in ApJ

Nature of Compton-thick Active Galactic Nuclei in "Non-merging'' Luminous Infrared Galaxies UGC 2608 and NGC 5135 Revealed with Broadband X-ray Spectroscopy

Authors: Satoshi Yamada, Yoshihiro Ueda, Atsushi Tanimoto, Saeko Oda, Masatoshi Imanishi, Yoshiki Toba, Claudio Ricci

Abstract: We have analyzed the broadband X-ray spectra of active galactic nuclei (AGNs) in two "non-merging" luminous infrared galaxies (LIRGs) UGC 2608 and NGC 5135, utilizing the data of NuSTAR, Suzaku, XMM-Newton, and Chandra. Applying the X-ray clumpy-torus model (XCLUMPY: Tanimoto et al. 2019), we find that both sources have similar spectra characterized by Compton-thick (CT) absorption (… ▽ More We have analyzed the broadband X-ray spectra of active galactic nuclei (AGNs) in two "non-merging" luminous infrared galaxies (LIRGs) UGC 2608 and NGC 5135, utilizing the data of NuSTAR, Suzaku, XMM-Newton, and Chandra. Applying the X-ray clumpy-torus model (XCLUMPY: Tanimoto et al. 2019), we find that both sources have similar spectra characterized by Compton-thick (CT) absorption ($N_{\rm H} \sim$ 5-7 $\times$ $10^{24}$ cm$^{-2}$) and small torus angular width ($σ$ $<$ 20$^{\circ}$). The intrinsic 2-10 keV luminosities are $3.9^{+2.2}_{-1.7}$ $\times$ $10^{43}$ erg s$^{-1}$ (UGC 2608) and $2.0^{+3.3}_{-1.0}$ $\times$ $10^{43}$ erg s$^{-1}$ (NGC 5135). The [O IV]-to-nuclear-12 $μ$m luminosity ratios are larger than those of typical Seyferts, which are consistent with the torus covering factors ($C_{\rm T} \lesssim$ 0.7) estimated from the torus angular widths and column densities by X-ray spectroscopy. The torus covering factors and Eddington ratios ($λ_{\rm Edd} \sim$ 0.1) follow the relation found by Ricci et al. (2017c) for local AGNs, implying that their tori become geometrically thin due to significant radiation pressure of the AGN that blows out some part of the tori. These results indicate that the CT AGNs in these "non-merger" LIRGs are just a normal AGN population seen edge-on through a large line-of-sight column density. They are in contrast to the buried CT AGNs in late-stage mergers that have large torus covering factors even at large Eddington ratios. △ Less

Submitted 18 May, 2020; originally announced May 2020.

Comments: 15 pages, 12 figures, and 3 tables, accepted for publication in ApJ

Simulations of the Early Post-Bounce Phase of Core-Collapse Supernovae in Three-Dimensional Space with Full Boltzmann Neutrino Transport

Authors: Wakana Iwakami, Hirotada Okawa, Hiroki Nagakura, Akira Harada, Shun Furusawa, Kosuke Sumiyoshi, Hideo Matsufuru, Shoichi Yamada

Abstract: We report on the core-collapse supernova simulation we conducted for a 11.2 M progenitor model in three-dimensional space up to 20 ms after bounce, using a radiation hydrodynamics code with full Boltzmann neutrino transport. We solve the six-dimensional Boltzmann equations for three neutrino species and the three-dimensional compressible Euler equations with Furusawa and Togashi's nuclear equation… ▽ More We report on the core-collapse supernova simulation we conducted for a 11.2 M progenitor model in three-dimensional space up to 20 ms after bounce, using a radiation hydrodynamics code with full Boltzmann neutrino transport. We solve the six-dimensional Boltzmann equations for three neutrino species and the three-dimensional compressible Euler equations with Furusawa and Togashi's nuclear equation of state. We focus on the prompt convection at 10 ms after bounce and investigate how neutrinos are transported in the convective matter. We apply a new analysis based on the eigenvalues and eigenvectors of the Eddington tensor and make a comparison between the Boltzmann transport results and the M1 closure approximation in the transition regime between the optically thick and thin limits. We visualize the eigenvalues and eigenvectors using an ellipsoid, in which each principal axis is parallel to one of the eigenvectors and has a length proportional to the corresponding eigenvalue. This approach enables us to understand the difference between the Eddington tensor derived directly from the Boltzmann simulation and the one given by the M1 prescription from a new perspective. We find that the longest principal axis of the ellipsoid is almost always nearly parallel to the energy flux in the M1 closure approximation whereas in the Boltzmann simulation it becomes perpendicular in some transition regions, where the mean free path is 0.1 times the radius. In three spatial dimensions, the convective motions make it difficult to predict where this happens and to possibly improve the closure relation there. △ Less

Submitted 31 August, 2020; v1 submitted 5 April, 2020; originally announced April 2020.

Comments: 29 pages, 22 figures, submitted to ApJ

The Boltzmann-radiation-hydrodynamics Simulations of Core-collapse Supernovae with Different Equations of State: the Role of Nuclear Composition and the Behavior of Neutrinos

Authors: Akira Harada, Hiroki Nagakura, Wakana Iwakami, Hirotada Okawa, Shun Furusawa, Kohsuke Sumiyoshi, Hideo Matsufuru, Shoichi Yamada

Abstract: Using the Boltzmann-radiation-hydrodynamics code, which solves the Boltzmann equation for neutrino transport, we present the results of the simulations with the nuclear equations of state (EOSs) of Lattimer and Swesty (LS) and Furusawa and Shen (FS). We extend the simulation time of the LS model and conduct thorough investigations, though our previous paper briefly reported some of the results. On… ▽ More Using the Boltzmann-radiation-hydrodynamics code, which solves the Boltzmann equation for neutrino transport, we present the results of the simulations with the nuclear equations of state (EOSs) of Lattimer and Swesty (LS) and Furusawa and Shen (FS). We extend the simulation time of the LS model and conduct thorough investigations, though our previous paper briefly reported some of the results. Only the LS model shows the shock revival. This seems to originate from the nuclear composition: the different nuclear composition results in the different energy loss by photodissociation and hence the different strength of the prompt convection and the later neutrino-driven convection. The protoneutron star seen in the FS model is more compact than that in the LS model because the existence of multinuclear species softens the EOS. For the behavior of neutrinos, we examined the flux and the Eddington tensor of neutrinos. In the optically thick region, the diffusion of neutrinos and the dragging by the motion of matter determine the flux. In the optically thin region, the free-streaming determines it. The Eddington tensor is compared with that obtained from the M1-closure relation. The M1-closure scheme overestimates the contribution from the velocity-dependent terms in the semitransparent region. △ Less

Submitted 17 November, 2020; v1 submitted 19 March, 2020; originally announced March 2020.

Comments: 30 pages, 24 figures, 1 explosion, published in ApJ

Neutrino transport with Monte Carlo method: I. Towards fully consistent implementation of nucleon recoils in core-collapse supernova simulations

Authors: Chinami Kato, Hiroki Nagakura, Yusuke Hori, Shoichi Yamada

Abstract: The small energy exchange via nucleon recoils in neutrino-nucleon scattering is now supposed to be one of the important factors for successful explosion of core-collapse supernovae (CCSNe) as they can change neutrino spectra through accumulation of a large number of scatterings. In finite-difference methods employed for neutrino transport in CCSN simulations, we normally can not afford to deploy a… ▽ More The small energy exchange via nucleon recoils in neutrino-nucleon scattering is now supposed to be one of the important factors for successful explosion of core-collapse supernovae (CCSNe) as they can change neutrino spectra through accumulation of a large number of scatterings. In finite-difference methods employed for neutrino transport in CCSN simulations, we normally can not afford to deploy a large enough number of energy bins needed to resolve this small energy exchange and sub-grid techniques are employed one way or another. In this paper we study quantitatively with the Monte Carlo (MC) method how well such a treatment performs. We first investigate the effects of nucleon recoils on the neutrino spectra and confirm that the average energy is reduced by $\sim$15% for heavy-lepton neutrinos and by much smaller quantities for other types of neutrinos in a typical post-bounce situation. It is also observed that the nucleon scattering dominates the electron scattering in the thermalization of neutrino spectra in all flavors. We then study possible artifacts that the coarse energy grid may produce in the finite-difference methods. In order to mimic the latter calculation, we re-distribute MC particles in each energy bin after a certain interval in a couple of ways and study how the results are affected and depend on the energy-resolution. We also discuss possible implications of our results for the finite-difference methods. △ Less

Submitted 29 January, 2020; originally announced January 2020.

Comments: 22 pages, 11 figures, submitted to ApJS

Are Faint Supernovae Responsible for Carbon-Enhanced Metal-Poor Stars?

Authors: Yutaka Komiya, Takuma Suda, Shimako Yamada, Masayuki Y. Fujimoto

Abstract: Mixing and fallback models in faint supernova models are supposed to reproduce the abundance patterns of observed carbon-enhanced metal-poor (CEMP) stars in the Galactic halo. A fine tuning of the model parameters for individual stars is required to reproduce the observed ratios of carbon to iron. We focus on extremely metal-poor stars formed out of the ejecta from the mixing and fallback models u… ▽ More Mixing and fallback models in faint supernova models are supposed to reproduce the abundance patterns of observed carbon-enhanced metal-poor (CEMP) stars in the Galactic halo. A fine tuning of the model parameters for individual stars is required to reproduce the observed ratios of carbon to iron. We focus on extremely metal-poor stars formed out of the ejecta from the mixing and fallback models using a chemical evolution model. Our chemical evolution models take into account the contribution of individual stars to chemical enrichment in host halos together with their evolution in the context of the hierarchical clustering. Parametrized models of mixing and fallback models for Pop. III faint supernovae are implemented in the chemical evolution models with merger trees to reproduce the observed CEMP stars. A variety of choices for model parameters on star formation and metal-pollution by faint supernovae is unable to reproduce the observed stars with [Fe/H] < -4 and [C/H] > -2, which are the majority of CEMP stars among the lowest metallicity stars. Only possible solution is to form stars from small ejecta mass, which produces an inconsistent metallicity distribution function. We conclude that not all the CEMP stars are explicable by the mixing and fallback models. We also tested the contribution of binary mass transfers from AGB stars that are also supposed to reproduce the abundances of known CEMP stars. This model reasonably reproduces the distribution of carbon and iron abundances simultaneously only if we assume that long-period binaries are favored at [Fe/H] < -3.5. △ Less

Submitted 6 January, 2020; originally announced January 2020.

Comments: 17 pages, 19 figures, accepted for publication in the Astrophysical Journal

Spatio-temporal linear instability analysis for arbitrary dispersion relations on the Lefschetz thimble in multidimensional spacetime

Authors: Taiki Morinaga, Shoichi Yamada

Abstract: In linear stability analysis of field quantities described by partial differential equations, the well-established classical theory is all but impossible to apply to concrete problems in its entirety even for uniform backgrounds when the spatial dimension is more than 1. In this study, using the Lefschetz thimble method, we develop a new formalism to give an explicit expression to the asymptotic b… ▽ More In linear stability analysis of field quantities described by partial differential equations, the well-established classical theory is all but impossible to apply to concrete problems in its entirety even for uniform backgrounds when the spatial dimension is more than 1. In this study, using the Lefschetz thimble method, we develop a new formalism to give an explicit expression to the asymptotic behavior of linear perturbations. It is not only more mathematically rigorous than the previous theory but also useful practically in its applications to realistic problems, and, as such, has an impact on broad subjects in physics. △ Less

Submitted 23 December, 2019; originally announced December 2019.

Comments: 10 pages, 18 figures, accepted for publication in Physical Review Research