Showing 1–50 of 176 results for author: Inoue, T

Information content in anisotropic cosmological fields: Impact of different multipole expansion scheme for galaxy density and ellipticity correlations

Authors: Takuya Inoue, Teppei Okumura, Shohei Saga, Atsushi Taruya

Abstract: Multipole expansions have been often used for extracting cosmological information from anisotropic quantities in observation. However, which basis of the expansion is best suited to quantify the anisotropies is not a trivial question in any summary statistics. In this paper, using the Fisher matrix formalism, we investigate the information content in multipole moments of the power spectra of galax… ▽ More Multipole expansions have been often used for extracting cosmological information from anisotropic quantities in observation. However, which basis of the expansion is best suited to quantify the anisotropies is not a trivial question in any summary statistics. In this paper, using the Fisher matrix formalism, we investigate the information content in multipole moments of the power spectra of galaxy density and intrinsic ellipticity fields in linear theory from the Alcock-Paczynski effect and redshift-space distortions (RSD). We consider two expansion schemes, the associated Legendre basis as well as the standard Legendre basis conventionally used in literature. We find that the information in the multipoles of the intrinsic alignment (IA) power spectra in the associated Legendre basis converges more slowly to that in the full 2D power spectra than in the Legendre basis. This trend is particularly significant when we consider a high number density sample. In this high number density case, we show that the errors on the Hubble parameter obtained from the multipoles of the IA cross- and auto-power spectra in the associated Legendre basis are respectively about $6 \%$ and $10 \%$ larger than the full 2D case even when we use multipoles up to $\ell = 6$. Our results demonstrate that the choice of basis is arbitrary but changes the information content encoded in multipoles depending on the sample and statistics under consideration. △ Less

Submitted 28 June, 2024; originally announced June 2024.

Comments: 19 pages, 7 figures

Report number: YITP-24-77

The Impact of Stellar Radiative Feedback on Formation of Young Massive Clusters via Fast HI Gas Collisions

Authors: Ryunosuke Maeda, Tsuyoshi Inoue, Kazuyuki Omukai, Yasuo Fukui, Kisetsu Tsuge

Abstract: Young massive clusters (YMCs) are dense aggregates of young stars and are often speculated as potential precursors to globular clusters. However, the formation mechanism of massive and compact gas clumps that precede YMCs remains unknown. In this paper, we study the formation of such massive clumps via fast HI gas collisions (~100 km/s) as suggested by recent observations and their subsequent evol… ▽ More Young massive clusters (YMCs) are dense aggregates of young stars and are often speculated as potential precursors to globular clusters. However, the formation mechanism of massive and compact gas clumps that precede YMCs remains unknown. In this paper, we study the formation of such massive clumps via fast HI gas collisions (~100 km/s) as suggested by recent observations and their subsequent evolution into YMCs by using three-dimensional magnetohydrodynamics simulations involving self-gravity and detailed thermal/chemical processes. In particular, the impact of ionization feedback from stellar radiation is included in an approximate fashion where the temperature within the HII regions is elevated to 10,000 K, while supernova feedback is not included. We examine whether the resulting massive clumps can survive this ionization feedback and evolve into YMCs. Our simulations reveal the emergence of gas clumps that do not only possess substantial mass (~10^5 M_sun) but also sufficient compactness (~5 pc). Notably, these clumps exhibit significantly higher escape velocities compared to the sound speed of the HII region, indicating effective gravitational retention of gas against feedback-induced evaporation. Consequently, these conditions foster efficient star formation within the massive gas clumps, ultimately leading to their evolution into YMCs. We also perform simulations involving lower-velocity gas collisions, approximately 15 km/s, typical shock velocities induced by galactic superbubbles. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: 15 pages, 13 figures, and 1 table; Accepted for publication in ApJ

Systematic comparison of neural networks used in discovering strong gravitational lenses

Authors: Anupreeta More, Raoul Canameras, Anton T. Jaelani, Yiping Shu, Yuichiro Ishida, Kenneth C. Wong, Kaiki Taro Inoue, Stefan Schuldt, Alessandro Sonnenfeld

Abstract: Efficient algorithms are being developed to search for strong gravitational lens systems owing to increasing large imaging surveys. Neural networks have been successfully used to discover galaxy-scale lens systems in imaging surveys such as the Kilo Degree Survey, Hyper-Suprime Cam (HSC) Survey and Dark Energy Survey over the last few years. Thus, it has become imperative to understand how some of… ▽ More Efficient algorithms are being developed to search for strong gravitational lens systems owing to increasing large imaging surveys. Neural networks have been successfully used to discover galaxy-scale lens systems in imaging surveys such as the Kilo Degree Survey, Hyper-Suprime Cam (HSC) Survey and Dark Energy Survey over the last few years. Thus, it has become imperative to understand how some of these networks compare, their strengths and the role of the training datasets as most of the networks make use of supervised learning algorithms. In this work, we present the first-of-its-kind systematic comparison and benchmarking of networks from four teams that have analysed the HSC Survey data. Each team has designed their training samples and developed neural networks independently but coordinated apriori in reserving specific datasets strictly for test purposes. The test sample consists of mock lenses, real (candidate) lenses and real non-lenses gathered from various sources to benchmark and characterise the performance of each of the network. While each team's network performed much better on their own constructed test samples compared to those from others, all networks performed comparable on the test sample with real (candidate) lenses and non-lenses. We also investigate the impact of swapping the training samples amongst the teams while retaining the same network architecture. We find that this resulted in improved performance for some networks. These results have direct implications on measures to be taken for lens searches with upcoming imaging surveys such as the Rubin-Legacy Survey of Space and Time, Roman and Euclid. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: 13 pages, 8 figures

High-mass star formation in the Large Magellanic Cloud triggered by colliding HI flows

Authors: K. Tsuge, H. Sano, K. Tachihara, K. Bekki, K. Tokuda, T. Inoue, N. Mizuno, A. Kawamura, T. Onishi, Y. Fukui

Abstract: The galactic tidal interaction is a possible mechanism to trigger the active star formation in galaxies. The recent analyses using the HI data in the Large Magellanic Cloud (LMC) proposed that the tidally driven HI flow, the L-component, is colliding with the LMC disk, the D-component, and is triggering high-mass star formation toward the active star-forming regions R136 and N44. In order to explo… ▽ More The galactic tidal interaction is a possible mechanism to trigger the active star formation in galaxies. The recent analyses using the HI data in the Large Magellanic Cloud (LMC) proposed that the tidally driven HI flow, the L-component, is colliding with the LMC disk, the D-component, and is triggering high-mass star formation toward the active star-forming regions R136 and N44. In order to explore the role of the collision over the entire LMC disk, we investigated the I-component, the collision-compressed gas between the L- and D-components, over the LMC disk, and found that 74% of the O/WR stars are located toward the I-component, suggesting their formation in the colliding gas. We compared four star-forming regions (R136, N44, N11, N77-N79-N83 complex). We found a positive correlation between the number of high-mass stars and the compressed gas pressure generated by collisions, suggesting that the pressure may be a key parameter in star formation. △ Less

Submitted 8 May, 2024; originally announced May 2024.

Comments: 33 pages, 18 figures, 3 Tables, accepted for publication in PASJ. arXiv admin note: text overlap with arXiv:2010.08816

Internal 1000 AU-scale Structures of the R CrA Cluster-forming Cloud -- I: Filamentary Structures

Authors: Kengo Tachihara, Naofumi Fukaya, Kazuki Tokuda, Yasumasa Yamasaki, Takeru Nishioka, Daisei Abe, Tsuyoshi Inoue, Naoto Harada, Ayumu Shoshi, Shingo Nozaki, Asako Sato, Mitsuki Omura, Kakeru Fujishiro, Misato Fukagawa, Masahiro N. Machida, Takahiro Kanai, Yumiko Oasa, Toshikazu Onishi, Kazuya Saigo, Yasuo Fukui

Abstract: We report on ALMA ACA observations of a high-density region of the Corona Australis cloud forming a young star cluster, and the results of resolving internal structures. In addition to embedded Class 0/I protostars in continuum, a number of complex dense filamentary structures are detected in the C18O and SO lines by the 7m array. These are sub-structures of the molecular clump that are detected b… ▽ More We report on ALMA ACA observations of a high-density region of the Corona Australis cloud forming a young star cluster, and the results of resolving internal structures. In addition to embedded Class 0/I protostars in continuum, a number of complex dense filamentary structures are detected in the C18O and SO lines by the 7m array. These are sub-structures of the molecular clump that are detected by the TP array as the extended emission. We identify 101 and 37 filamentary structures with a few thousand AU widths in C18O and SO, respectively, called as feathers. The typical column density of the feathers in C18O is about 10^{22} cm^{-2}, and the volume density and line mass are ~ 10^5 cm^{-3}, and a few times M_{sun} pc^{-1}, respectively. This line mass is significantly smaller than the critical line mass expected for cold and dense gas. These structures have complex velocity fields, indicating a turbulent internal property. The number of feathers associated with Class 0/I protostars is only ~ 10, indicating that most of them do not form stars but rather being transient structures. The formation of feathers can be interpreted as a result of colliding gas flow as the morphology well reproduced by MHD simulations, supported by the the presence of HI shells in the vicinity. The colliding gas flows may accumulate gas and form filaments and feathers, and trigger the active star formation of the R CrA cluster. △ Less

Submitted 17 April, 2024; originally announced April 2024.

Comments: 24 pages, 13 figures; Accepted for publication in ApJ

Bell Instability-Mediated Diffusive Shock Acceleration at Supernova Blast Wave Shock Propagating in the ISM

Authors: Tsuyoshi Inoue, Alexandre Marcowith, Gwenael Giacinti

Abstract: Supernova blast wave shock is a very important site of cosmic-ray acceleration. However, the detailed physical process of acceleration, in particular, non-linear interplay between cosmic-ray streaming and magnetic field amplification has not been studied under a realistic environment. In this paper, using a unique and novel numerical method, we study cosmic-ray acceleration at supernova blast wave… ▽ More Supernova blast wave shock is a very important site of cosmic-ray acceleration. However, the detailed physical process of acceleration, in particular, non-linear interplay between cosmic-ray streaming and magnetic field amplification has not been studied under a realistic environment. In this paper, using a unique and novel numerical method, we study cosmic-ray acceleration at supernova blast wave shock propagating in the interstellar medium with well-resolved magnetic field amplification by non-resonant hybrid instability (or Bell instability). We find that the magnetic field is mildly amplified under typical ISM conditions that leads to maximum cosmic-ray energy ~30 TeV for supernova remnants with age ~1000 years consistent with gamma-ray observations. The strength of the amplified magnetic field does not reach so-called saturation level, because cosmic-ray electric current towards the shock upstream has finite spatial extent, by which Bell instability cannot experience many e-folding times. △ Less

Submitted 5 April, 2024; originally announced April 2024.

Comments: ApJ accepted

Filamentary Network and Magnetic Field Structures Revealed with BISTRO in the High-Mass Star-Forming Region NGC2264 : Global Properties and Local Magnetogravitational Configurations

Authors: Jia-Wei Wang, Patrick M. Koch, Seamus D. Clarke, Gary Fuller, Nicolas Peretto, Ya-Wen Tang, Hsi-Wei Yen, Shih-Ping Lai, Nagayoshi Ohashi, Doris Arzoumanian, Doug Johnstone, Ray Furuya, Shu-ichiro Inutsuka, Chang Won Lee, Derek Ward-Thompson, Valentin J. M. Le Gouellec, Hong-Li Liu, Lapo Fanciullo, Jihye Hwang, Kate Pattle, Frédérick Poidevin, Mehrnoosh Tahani, Takashi Onaka, Mark G. Rawlings, Eun Jung Chung , et al. (132 additional authors not shown)

Abstract: We report 850 $μ$m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B-fields In STar forming Regions Observations (BISTRO) large program on the James Clerk Maxwell Telescope (JCMT). These data reveal a well-structured non-uniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30 deg from… ▽ More We report 850 $μ$m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B-fields In STar forming Regions Observations (BISTRO) large program on the James Clerk Maxwell Telescope (JCMT). These data reveal a well-structured non-uniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30 deg from north to east. Field strengths estimates and a virial analysis for the major clumps indicate that NGC 2264C is globally dominated by gravity while in 2264D magnetic, gravitational, and kinetic energies are roughly balanced. We present an analysis scheme that utilizes the locally resolved magnetic field structures, together with the locally measured gravitational vector field and the extracted filamentary network. From this, we infer statistical trends showing that this network consists of two main groups of filaments oriented approximately perpendicular to one another. Additionally, gravity shows one dominating converging direction that is roughly perpendicular to one of the filament orientations, which is suggestive of mass accretion along this direction. Beyond these statistical trends, we identify two types of filaments. The type-I filament is perpendicular to the magnetic field with local gravity transitioning from parallel to perpendicular to the magnetic field from the outside to the filament ridge. The type-II filament is parallel to the magnetic field and local gravity. We interpret these two types of filaments as originating from the competition between radial collapsing, driven by filament self-gravity, and the longitudinal collapsing, driven by the region's global gravity. △ Less

Submitted 23 January, 2024; originally announced January 2024.

Comments: Accepted for publication in the Astrophysical Journal. 43 pages, 32 figures, and 4 tables (including Appendix)

Probing warm and mixed dark matter models using lensing shift power spectrum

Authors: Kaiki Taro Inoue, Takumi Shinohara, Teruaki Suyama, Tomo Takahashi

Abstract: We argue that the lensing power spectrum of astrometric shift (lensing shift power spectrum) is a powerful tool of the clustering property of dark matter on subgalactic scales. First we give the formalism to probe the nature of dark matter by using the lensing shift power spectrum. Then, leveraging recent measurements of the lensing shift power spectrum on an angular scale of approximately $1~$arc… ▽ More We argue that the lensing power spectrum of astrometric shift (lensing shift power spectrum) is a powerful tool of the clustering property of dark matter on subgalactic scales. First we give the formalism to probe the nature of dark matter by using the lensing shift power spectrum. Then, leveraging recent measurements of the lensing shift power spectrum on an angular scale of approximately $1~$arcsec towards the gravitationally lensed quasar MG$\,$J0414+0534 at the redshift of $z_S=2.639$, we place constraints on the mass of warm dark matter (WDM) particles $m_{\rm WDM}$ and their fraction in a mixed dark matter (MDM) model $r_{\rm WDM}$, in which WDM and cold dark matter coexist. Although the constraint derived from the above single lensing system is not as strong as the existing constraints, as we show in this paper, the lensing shift power spectrum has a great potential to obtain much tighter constraints on WDM and MDM models through future observations, highlighting the importance of well-controlled systematic error considerations for achieving enhanced precision. △ Less

Submitted 29 December, 2023; originally announced December 2023.

Comments: 24 pages, 7 figures

Survey of Gravitationally lensed Objects in HSC Imaging (SuGOHI) $-$ X. Strong Lens Finding in The HSC-SSP using Convolutional Neural Networks

Authors: Anton T. Jaelani, Anupreeta More, Kenneth C. Wong, Kaiki T. Inoue, Dani C. -Y. Chao, Premana W. Premadi, Raoul Cañameras

Abstract: We apply a novel model based on convolutional neural networks (CNNs) to identify gravitationally-lensed galaxies in multi-band imaging of the Hyper Suprime Cam Subaru Strategic Program (HSC-SSP) Survey. The trained model is applied to a parent sample of 2 350 061 galaxies selected from the $\sim$ 800 deg$^2$ Wide area of the HSC-SSP Public Data Release 2. The galaxies in HSC Wide are selected base… ▽ More We apply a novel model based on convolutional neural networks (CNNs) to identify gravitationally-lensed galaxies in multi-band imaging of the Hyper Suprime Cam Subaru Strategic Program (HSC-SSP) Survey. The trained model is applied to a parent sample of 2 350 061 galaxies selected from the $\sim$ 800 deg$^2$ Wide area of the HSC-SSP Public Data Release 2. The galaxies in HSC Wide are selected based on stringent pre-selection criteria, such as multiband magnitudes, stellar mass, star formation rate, extendedness limit, photometric redshift range, etc. Initially, the CNNs provide a total of 20 241 cutouts with a score greater than 0.9, but this number is subsequently reduced to 1 522 cutouts by removing definite non-lenses for further inspection by human eyes. We discover 43 definite and 269 probable lenses, of which 97 are completely new. In addition, out of 880 potential lenses, we recovered 289 known systems in the literature. We identify 143 candidates from the known systems that had higher confidence in previous searches. Our model can also recover 285 candidate galaxy-scale lenses from the Survey of Gravitationally lensed Objects in HSC Imaging (SuGOHI), where a single foreground galaxy acts as the deflector. Even though group-scale and cluster-scale lens systems were not included in the training, a sample of 32 SuGOHI-c (i.e., group/cluster-scale systems) lens candidates was retrieved. Our discoveries will be useful for ongoing and planned spectroscopic surveys, such as the Subaru Prime Focus Spectrograph project, to measure lens and source redshifts in order to enable detailed lens modelling. △ Less

Submitted 12 December, 2023; originally announced December 2023.

Comments: Submitted to MNRAS, 16 pages, 13 figures. Comments welcome

The Gamma Ray Origin in RXJ0852.0-4622 Quantifying the Hadronic and Leptonic Components: Further Evidence for the Cosmic Ray Acceleration in Young Shell-type SNRs

Authors: Yasuo Fukui, Maki Aruga, Hidetoshi Sano, Takahiro Hayakawa, Tsuyoshi Inoue, Gavin Rowell, Sabrina Einecke, Kengo Tachihara

Abstract: Fukui et al. (2021) quantified the hadronic and leptonic gamma rays in the young TeV gamma ray shell-type supernova remnant (SNR) RXJ1713.7-3946 (RXJ1713), and demonstrated that the gamma rays are a combination of the hadronic and leptonic gamma ray components with a ratio of $\sim 6:4$ in gamma ray counts $N_\mathrm{g}$. This discovery, which adopted a new methodology of multiple-linear gamma-ray… ▽ More Fukui et al. (2021) quantified the hadronic and leptonic gamma rays in the young TeV gamma ray shell-type supernova remnant (SNR) RXJ1713.7-3946 (RXJ1713), and demonstrated that the gamma rays are a combination of the hadronic and leptonic gamma ray components with a ratio of $\sim 6:4$ in gamma ray counts $N_\mathrm{g}$. This discovery, which adopted a new methodology of multiple-linear gamma-ray decomposition, was the first quantification of the two gamma ray components. In the present work, we applied the same methodology to another TeV gamma ray shell-type SNR RX~J0852.0$-$4622 (RXJ0852) in the 3D space characterized by [the interstellar proton column density $N_{\mathrm{p}}$]-[the nonthermal X-ray count $N_{\mathrm{x}}$]-[$N_{\mathrm{g}}$], and quantified the hadronic and leptonic gamma ray components to have a ratio of $\sim 5:5$ in $N_{\mathrm{g}}$. The present work adopted fitting of two/three flat planes in the 3D space instead of a single flat plane, which allowed to suppress fitting errors. The quantification indicates that the hadronic and leptonic gamma rays are in the same order of magnitude in these two core-collapse SNRs, verifying the significant hadronic gamma ray components. We argue that the target interstellar protons, in particular their spatial distribution, are essential in any attempts to identify type of particles responsible for the gamma-ray emission. The present results confirm that the CR energy $\lesssim 100$\,TeV is compatible with a scheme that SNRs are the dominant source of these Galactic CRs. △ Less

Submitted 19 November, 2023; originally announced November 2023.

ALMA Observations of Supernova Remnant N49 in the Large Magellanic Cloud. II. Non-LTE Analysis of Shock-heated Molecular Clouds

Authors: H. Sano, Y. Yamane, J. Th. van Loon, K. Furuya, Y. Fukui, R. Z. E. Alsaberi, A. Bamba, R. Enokiya, M. D. Filipović, R. Indebetouw, T. Inoue, A. Kawamura, M. Lakićević, C. J. Law, N. Mizuno, T. Murase, T. Onishi, S. Park, P. P. Plucinsky, J. Rho, A. M. S. Richards, G. Rowell, M. Sasaki, J. Seok, P. Sharda , et al. (6 additional authors not shown)

Abstract: We present the first compelling evidence of shock-heated molecular clouds associated with the supernova remnant (SNR) N49 in the Large Magellanic Cloud (LMC). Using $^{12}$CO($J$ = 2-1, 3-2) and $^{13}$CO($J$ = 2-1) line emission data taken with the Atacama Large Millimeter/Submillimeter Array, we derived the H$_2$ number density and kinetic temperature of eight $^{13}$CO-detected clouds using the… ▽ More We present the first compelling evidence of shock-heated molecular clouds associated with the supernova remnant (SNR) N49 in the Large Magellanic Cloud (LMC). Using $^{12}$CO($J$ = 2-1, 3-2) and $^{13}$CO($J$ = 2-1) line emission data taken with the Atacama Large Millimeter/Submillimeter Array, we derived the H$_2$ number density and kinetic temperature of eight $^{13}$CO-detected clouds using the large velocity gradient approximation at a resolution of 3.5$''$ (~0.8 pc at the LMC distance). The physical properties of the clouds are divided into two categories: three of them near the shock front show the highest temperatures of ~50 K with densities of ~500-700 cm$^{-3}$, while other clouds slightly distant from the SNR have moderate temperatures of ~20 K with densities of ~800-1300 cm$^{-3}$. The former clouds were heated by supernova shocks, but the latter were dominantly affected by the cosmic-ray heating. These findings are consistent with the efficient production of X-ray recombining plasma in N49 due to thermal conduction between the cold clouds and hot plasma. We also find that the gas pressure is roughly constant except for the three shock-engulfed clouds inside or on the SNR shell, suggesting that almost no clouds have evaporated within the short SNR age of ~4800 yr. This result is compatible with the shock-interaction model with dense and clumpy clouds inside a low-density wind bubble. △ Less

Submitted 3 November, 2023; originally announced November 2023.

Comments: 12 pages, 5 figures, 1 table, accepted for publication in The Astrophysical Journal (ApJ)

Growth of Massive Molecular Cloud Filament by Accretion Flows I: Slow Shock Instability v.s. Ambipolar Diffusion

Authors: Daisei Abe, Tsuyoshi Inoue, Shu-ichiro Inutsuka

Abstract: The Herschel Gould Belt Survey showed that stars form in dense filaments in nearby molecular clouds. Recent studies suggest that massive filaments are bound by the slow shocks caused by accretion flows onto the filaments. The slow shock is known to be unstable to the corrugation deformation of the shock front. The corrugation instability could convert the accretion flow's ram pressure into turbule… ▽ More The Herschel Gould Belt Survey showed that stars form in dense filaments in nearby molecular clouds. Recent studies suggest that massive filaments are bound by the slow shocks caused by accretion flows onto the filaments. The slow shock is known to be unstable to the corrugation deformation of the shock front. The corrugation instability could convert the accretion flow's ram pressure into turbulent pressure that influences the width of the filament, which, according to theory, determines the self-gravitational fragmentation scale and core mass. In spite of its importance, the effect of slow shock instability on star-forming filaments has not been investigated. In addition, the linear dispersion relation obtained from the ideal magnetohydrodynamics (MHD) analysis shows that the most unstable wavelength of shock corrugation is infinitesimally small (or mean free path). In the scale of dense filaments, the effect of ambipolar diffusion can suppress the instability at small scales. This study investigates the influence of ambipolar diffusion on the instability of the slow shock. We perform two-dimensional MHD simulations to examine the linear growth of the slow shock instability, considering the effect of ambipolar diffusion. The results demonstrate that the most unstable scale of slow shock instability is approximately five times the length scale of ambipolar diffusion l_AD calculated using post-shock variables, where, l_AD corresponds to the scale where the magnetic Reynolds number for ambipolar diffusivity is unity. △ Less

Submitted 25 October, 2023; originally announced October 2023.

Comments: 20 pages, 12 figures, accepted by ApJ

An improved dynamical Poisson equation solver for self-gravity

Authors: Ryunosuke Maeda, Tsuyoshi Inoue, Shu-ichiro Inutsuka

Abstract: Since self-gravity is crucial in the structure formation of the universe, many hydrodynamics simulations with the effect of self-gravity have been conducted. The multigrid method is widely used as a solver for the Poisson equation of the self-gravity; however, the parallelization efficiency of the multigrid method becomes worse when we use a massively parallel computer, and it becomes inefficient… ▽ More Since self-gravity is crucial in the structure formation of the universe, many hydrodynamics simulations with the effect of self-gravity have been conducted. The multigrid method is widely used as a solver for the Poisson equation of the self-gravity; however, the parallelization efficiency of the multigrid method becomes worse when we use a massively parallel computer, and it becomes inefficient with more than $10^4$ cores, even for highly tuned codes. To perform large-scale parallel simulations ($> 10^4$ cores), developing a new gravity solver with good parallelization efficiency is beneficial. In this article, we develop a new self-gravity solver using the telegraph equation with a damping coefficient, $κ$. Parallelization is much easier than the case of the elliptic Poisson equation since the telegraph equation is a hyperbolic partial differential equation. We analyze convergence tests of our telegraph equations solver and determine that the best non-dimensional damping coefficient of the telegraph equations is $\tildeκ \simeq 2.5$. We also show that our method can maintain high parallelization efficiency even for massively parallel computations due to the hyperbolic nature of the telegraphic equation by weak-scaling tests. If the time step of the calculation is determined by heating/cooling or chemical reactions, rather than the CFL condition, our method may provide the method for calculating self-gravity faster than other previously known methods such as the fast Fourier transform and multigrid iteration solvers because gravitational phase velocity determined by the CFL condition using these timescales is much larger than the fluid velocity plus sound speed. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: 8 pages, 13 figures, accepted to MNRAS

Global deceleration and inward movements of X-ray knots and rims of RCW 103

Authors: Hiromasa Suzuki, Takaaki Tanaka, Tsuyoshi Inoue, Hiroyuki Uchida, Takuto Narita

Abstract: Kinematics of shocks, ejecta knots, and the compact remnant of a supernova remnant gives an insight into the nature of the progenitor and surrounding environment. We report on a proper motion measurement of X-ray knots and rims of the magnetar-hosting supernova remnant RCW 103. Chandra data obtained in three epochs, 1999, 2010, and 2016 are used. We find a global deceleration of 12 knots and rims… ▽ More Kinematics of shocks, ejecta knots, and the compact remnant of a supernova remnant gives an insight into the nature of the progenitor and surrounding environment. We report on a proper motion measurement of X-ray knots and rims of the magnetar-hosting supernova remnant RCW 103. Chandra data obtained in three epochs, 1999, 2010, and 2016 are used. We find a global deceleration of 12 knots and rims both in northern and southern regions within the last $\sim 24$ yrs, even though its age is thought to be larger than 2 kyr. Some of them even changed their moving directions from outward ($\sim 1,000$ km s$^{-1}$) to inward ($\sim -2,000$ km s$^{-1}$). Our findings can be explained with a collision with a high-density medium both in the northern and southern edges of the remnant, although the remnant may still be expanding in the wind-blown cavity. The proper motion of the associated magnetar 1E161348$-$5055 is possibly detected with a velocity of $\approx 500$ km s$^{-1}$. △ Less

Submitted 6 October, 2023; originally announced October 2023.

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

Measuring Line-of-sight Distances to Haloes with Astrometric Lensing B-mode

Authors: Kaiki Taro Inoue

Abstract: Relative astrometric shifts between multiply lensed images provide a valuable tool to investigate haloes in the intergalactic space. In strong lens systems in which a single lens plays the primary role in producing multiple images, the gravitational force exerted by line-of-sight (LOS) haloes can slightly change the relative positions of multiply lensed images produced by the dominant lens. In suc… ▽ More Relative astrometric shifts between multiply lensed images provide a valuable tool to investigate haloes in the intergalactic space. In strong lens systems in which a single lens plays the primary role in producing multiple images, the gravitational force exerted by line-of-sight (LOS) haloes can slightly change the relative positions of multiply lensed images produced by the dominant lens. In such cases, a LOS halo positioned sufficiently far from the dominant lens along the LOS can create a pattern in the reduced deflection angle that corresponds to the B-mode (magnetic or divergence-free mode). By measuring both the B-mode and E-mode (electric or rotation-free mode), we can determine the LOS distance ratios, as well as the 'bare' convergence and shear perturbations in the absence of the dominant lens. However, scale variations in the distance ratio lead to mass-sheet transformations in the background lens plane, introducing some uncertainty in the distance ratio estimation. This uncertainty can be significantly reduced by measuring the time delays between the lensed images. Additionally, if we obtain the redshift values of both the dominant and perturbing haloes, along with the time delays between the multiply lensed images that are affected by the haloes, the B-mode can break the degeneracy related to mass-sheet transformations in both the foreground and background lens planes. Therefore, measuring the astrometric lensing B-mode has the potential to substantially decrease the uncertainty in determining the Hubble constant. △ Less

Submitted 19 September, 2023; originally announced September 2023.

Comments: 15 pages, 16 figures

Predicting reliable H$_2$ column density maps from molecular line data using machine learning

Authors: Yoshito Shimajiri, Yasutomo Kawanishi, Shinji Fujita, Yusuke Miyamoto, Atsushi M. Ito, Doris Arzoumanian, Philippe André, Atsushi Nishimura, Kazuki Tokuda, Hiroyuki Kaneko, Shunya Takekawa, Shota Ueda, Toshikazu Onishi, Tsuyoshi Inoue, Shimpei Nishimoto, Ryuki Yoneda

Abstract: The total mass estimate of molecular clouds suffers from the uncertainty in the H$_2$-CO conversion factor, the so-called $X_{\rm CO}$ factor, which is used to convert the $^{12}$CO (1--0) integrated intensity to the H$_2$ column density. We demonstrate the machine learning's ability to predict the H$_2$ column density from the $^{12}$CO, $^{13}$CO, and C$^{18}$O (1--0) data set of four star-formi… ▽ More The total mass estimate of molecular clouds suffers from the uncertainty in the H$_2$-CO conversion factor, the so-called $X_{\rm CO}$ factor, which is used to convert the $^{12}$CO (1--0) integrated intensity to the H$_2$ column density. We demonstrate the machine learning's ability to predict the H$_2$ column density from the $^{12}$CO, $^{13}$CO, and C$^{18}$O (1--0) data set of four star-forming molecular clouds; Orion A, Orion B, Aquila, and M17. When the training is performed on a subset of each cloud, the overall distribution of the predicted column density is consistent with that of the Herschel column density. The total column density predicted and observed is consistent within 10\%, suggesting that the machine learning prediction provides a reasonable total mass estimate of each cloud. However, the distribution of the column density for values $> \sim 2 \times 10^{22}$ cm$^{-2}$, which corresponds to the dense gas, could not be predicted well. This indicates that molecular line observations tracing the dense gas are required for the training. We also found a significant difference between the predicted and observed column density when we created the model after training the data on different clouds. This highlights the presence of different $X_{\rm CO}$ factors between the clouds, and further training in various clouds is required to correct for these variations. We also demonstrated that this method could predict the column density toward the area not observed by Herschel if the molecular line and column density maps are available for the small portion, and the molecular line data are available for the larger areas. △ Less

Submitted 13 September, 2023; originally announced September 2023.

Comments: Accepted for publication in MNRAS

An ALMA Glimpse of Dense Molecular Filaments Associated with High-mass Protostellar Systems in the Large Magellanic Cloud

Authors: Kazuki Tokuda, Naoto Harada, Kei E. I. Tanaka, Tsuyoshi Inoue, Takashi Shimonishi, Yichen Zhang, Marta Sewiło, Yuri Kunitoshi, Ayu Konishi, Yasuo Fukui, Akiko Kawamura, Toshikazu Onishi, Masahiro N. Machida

Abstract: Recent millimeter/sub-millimeter facilities have revealed the physical properties of filamentary molecular clouds in relation to high-mass star formation. A uniform survey of the nearest, face-on star-forming galaxy, the Large Magellanic Cloud (LMC), complements the Galactic knowledge. We present ALMA survey data with a spatial resolution of $\sim$0.1 pc in the 0.87 mm continuum and HCO$^{+}$(4-3)… ▽ More Recent millimeter/sub-millimeter facilities have revealed the physical properties of filamentary molecular clouds in relation to high-mass star formation. A uniform survey of the nearest, face-on star-forming galaxy, the Large Magellanic Cloud (LMC), complements the Galactic knowledge. We present ALMA survey data with a spatial resolution of $\sim$0.1 pc in the 0.87 mm continuum and HCO$^{+}$(4-3) emission toward 30 protostellar objects with luminosities of 10$^4$-10$^{5.5}$ $L_{\odot}$ in the LMC. The spatial distributions of the HCO$^{+}$(4-3) line and thermal dust emission are well correlated, indicating that the line effectively traces dense, filamentary gas with an H$_2$ volume density of $\gtrsim$10$^5$ cm$^{-3}$ and a line mass of $\sim$10$^3$-10$^{4}$ $M_{\odot}$ pc$^{-1}$. Furthermore, we obtain an increase in the velocity linewidths of filamentary clouds, which follows a power-law dependence on their H$_2$ column densities with an exponent of $\sim$0.5. This trend is consistent with observations toward filamentary clouds in nearby star-forming regions withiin $ \lesssim$1 kpc from us and suggests enhanced internal turbulence within the filaments owing to surrounding gas accretion. Among the 30 sources, we find that 14 are associated with hub-filamentary structures, and these complex structures predominantly appear in protostellar luminosities exceeding $\sim$5 $\times$10$^4$ $L_{\odot}$. The hub-filament systems tend to appear in the latest stages of their natal cloud evolution, often linked to prominent H$\;${\sc ii} regions and numerous stellar clusters. Our preliminary statistics suggest that the massive filaments accompanied by hub-type complex features may be a necessary intermediate product in forming extremely luminous high-mass stellar systems capable of ultimately dispersing the parent cloud. △ Less

Submitted 10 August, 2023; originally announced August 2023.

Comments: 21 pages, 8 figures, 4 tables, accepted for publication in ApJ

Metallicity Dependence of Molecular Cloud Hierarchical Structure at Early Evolutionary Stages

Authors: Masato I. N. Kobayashi, Kazunari Iwasaki, Kengo Tomida, Tsuyoshi Inoue, Kazuyuki Omukai, Kazuki Tokuda

Abstract: The formation of molecular clouds out of HI gas is the first step toward star formation. Its metallicity dependence plays a key role to determine star formation through the cosmic history. Previous theoretical studies with detailed chemical networks calculate thermal equilibrium states and/or thermal evolution under one-zone collapsing background. The molecular cloud formation in reality, however,… ▽ More The formation of molecular clouds out of HI gas is the first step toward star formation. Its metallicity dependence plays a key role to determine star formation through the cosmic history. Previous theoretical studies with detailed chemical networks calculate thermal equilibrium states and/or thermal evolution under one-zone collapsing background. The molecular cloud formation in reality, however, involves supersonic flows, and thus resolving the cloud internal turbulence/density structure in three dimension is still essential. We here perform magnetohydrodynamics simulations of 20 km s^-1 converging flows of Warm Neutral Medium (WNM) with 1 micro Gauss mean magnetic field in the metallicity range from the Solar (1.0 Zsun) to 0.2 Zsun environment. The Cold Neutral Medium (CNM) clumps form faster with higher metallicity due to more efficient cooling. Meanwhile, their mass functions commonly follow dn/dm proportional to m^-1.7 at three cooling times regardless of the metallicity. Their total turbulence power also commonly shows the Kolmogorov spectrum with its 80 percent in the solenoidal mode, while the CNM volume alone indicates the transition towards the Larson's law. These similarities measured at the same time in the unit of the cooling time suggest that the molecular cloud formation directly from the WNM alone requires a longer physical time in a lower metallicity environment in the 1.0--0.2 Zsun range. To explain the rapid formation of molecular clouds and subsequent massive star formation possibly within less than 10 Myr as observed in the Large/Small Magellanic Clouds (LMC/SMC), the HI gas already contains CNM volume instead of pure WNM. △ Less

Submitted 23 July, 2023; v1 submitted 3 July, 2023; originally announced July 2023.

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

MHD Simulation in Galactic Center Region with Radiative Cooling and Heating

Authors: Kensuke Kakiuchi, Takeru. K. Suzuki, Shu-ichiro Inutsuka, Tsuyoshi Inoue, Jiro Shimoda

Abstract: We investigate the role of magnetic field on the gas dynamics in a galactic bulge region by three dimensional simulations with radiative cooling and heating. While high-temperature corona with $T>10^6\ {\rm K}$ is formed in the halo regions, the temperature near the mid-plane is $\lesssim 10^4\ {\rm K}$ following the thermal equilibrium curve determined by the radiative cooling and heating. Althou… ▽ More We investigate the role of magnetic field on the gas dynamics in a galactic bulge region by three dimensional simulations with radiative cooling and heating. While high-temperature corona with $T>10^6\ {\rm K}$ is formed in the halo regions, the temperature near the mid-plane is $\lesssim 10^4\ {\rm K}$ following the thermal equilibrium curve determined by the radiative cooling and heating. Although the thermal energy of the interstellar gas is lost by radiative cooling, the saturation level of the magnetic field strength does not significantly depend on the radiative cooling and heating. The magnetic field strength is amplified to $10\ {\rm μG}$ on average, and reaches several hundred ${\rm μG}$ locally. We find the formation of magnetically dominated regions at mid-latitudes in the case with the radiative cooling and heating, which is not seen in the case without radiative effect. The vertical thickness of the mid-latitude regions is $50-150\ {\rm pc}$ at the radial location of $0.4-0.8\ {\rm kpc}$ from the galactic center, which is comparable to the observed vertical distribution of neutral atomic gas. When we take the average of different components of energy density integrated over the galactic bulge region, the magnetic energy is comparable to the thermal energy. We conclude that the magnetic field plays a substantial role in controlling the dynamical and thermal properties of the galactic bulge region. △ Less

Submitted 4 April, 2024; v1 submitted 27 June, 2023; originally announced June 2023.

Comments: Accepted to ApJ; 24 pages, 20 figures 3 tables

HOLISMOKES -- XI. Evaluation of supervised neural networks for strong-lens searches in ground-based imaging surveys

Authors: R. Canameras, S. Schuldt, Y. Shu, S. H. Suyu, S. Taubenberger, I. T. Andika, S. Bag, K. T. Inoue, A. T. Jaelani, L. Leal-Taixe, T. Meinhardt, A. Melo, A. More

Abstract: While supervised neural networks have become state of the art for identifying the rare strong gravitational lenses from large imaging data sets, their selection remains significantly affected by the large number and diversity of nonlens contaminants. This work evaluates and compares systematically the performance of neural networks in order to move towards a rapid selection of galaxy-scale strong… ▽ More While supervised neural networks have become state of the art for identifying the rare strong gravitational lenses from large imaging data sets, their selection remains significantly affected by the large number and diversity of nonlens contaminants. This work evaluates and compares systematically the performance of neural networks in order to move towards a rapid selection of galaxy-scale strong lenses with minimal human input in the era of deep, wide-scale surveys. We used multiband images from PDR2 of the HSC Wide survey to build test sets mimicking an actual classification experiment, with 189 strong lenses previously found over the HSC footprint and 70,910 nonlens galaxies in COSMOS. Multiple networks were trained on different sets of realistic strong-lens simulations and nonlens galaxies, with various architectures and data pre-processing. The overall performances strongly depend on the construction of the ground-truth training data and they typically, but not systematically, improve using our baseline residual network architecture. Improvements are found when applying random shifts to the image centroids and square root stretches to the pixel values, adding z band, or using random viewpoints of the original images, but not when adding difference images to subtract emission from the central galaxy. The most significant gain is obtained with committees of networks trained on different data sets, and showing a moderate overlap between populations of false positives. Nearly-perfect invariance to image quality can be achieved by training networks either with large number of bands, or jointly with the PSF and science frames. Overall, we show the possibility to reach a TPR0 as high as 60% for the test sets under consideration, which opens promising perspectives for pure selection of strong lenses without human input using the Rubin Observatory and other forthcoming ground-based surveys. △ Less

Submitted 5 June, 2023; originally announced June 2023.

Comments: 21 pages, 10 figures, submitted to A&A, comments are welcome

The JCMT BISTRO Survey: Studying the Complex Magnetic Field of L43

Authors: Janik Karoly, Derek Ward-Thompson, Kate Pattle, David Berry, Anthony Whitworth, Jason Kirk, Pierre Bastien, Tao-Chung Ching, Simon Coude, Jihye Hwang, Woojin Kwon, Archana Soam, Jia-Wei Wang, Tetsuo Hasegawa, Shih-Ping Lai, Keping Qiu, Doris Arzoumanian, Tyler L. Bourke, Do-Young Byun, Huei-Ru Vivien Chen, Wen Ping Chen, Mike Chen, Zhiwei Chen, Jungyeon Cho, Minho Choi , et al. (133 additional authors not shown)

Abstract: We present observations of polarized dust emission at 850 $μ$m from the L43 molecular cloud which sits in the Ophiuchus cloud complex. The data were taken using SCUBA-2/POL-2 on the James Clerk Maxwell Telescope as a part of the BISTRO large program. L43 is a dense ($N_{\rm H_2}\sim 10^{22}$-10$^{23}$ cm$^{-2}$) complex molecular cloud with a submillimetre-bright starless core and two protostellar… ▽ More We present observations of polarized dust emission at 850 $μ$m from the L43 molecular cloud which sits in the Ophiuchus cloud complex. The data were taken using SCUBA-2/POL-2 on the James Clerk Maxwell Telescope as a part of the BISTRO large program. L43 is a dense ($N_{\rm H_2}\sim 10^{22}$-10$^{23}$ cm$^{-2}$) complex molecular cloud with a submillimetre-bright starless core and two protostellar sources. There appears to be an evolutionary gradient along the isolated filament that L43 is embedded within, with the most evolved source closest to the Sco OB2 association. One of the protostars drives a CO outflow that has created a cavity to the southeast. We see a magnetic field that appears to be aligned with the cavity walls of the outflow, suggesting interaction with the outflow. We also find a magnetic field strength of up to $\sim$160$\pm$30 $μ$G in the main starless core and up to $\sim$90$\pm$40 $μ$G in the more diffuse, extended region. These field strengths give magnetically super- and sub-critical values respectively and both are found to be roughly trans-Alfvénic. We also present a new method of data reduction for these denser but fainter objects like starless cores. △ Less

Submitted 22 May, 2023; v1 submitted 18 May, 2023; originally announced May 2023.

Comments: Accepted for publication in ApJ. 23 pages, 9 figures (7 main text, 2 appendix)

First BISTRO observations of the dark cloud Taurus L1495A-B10: the role of the magnetic field in the earliest stages of low-mass star formation

Authors: Derek Ward-Thompson, Janik Karoly, Kate Pattle, Anthony Whitworth, Jason Kirk, David Berry, Pierre Bastien, Tao-Chung Ching, Simon Coude, Jihye Hwang, Woojin Kwon, Archana Soam, Jia-Wei Wang, Tetsuo Hasegawa, Shih-Ping Lai, Keping Qiu, Doris Arzoumanian, Tyler L. Bourke, Do-Young Byun, Huei-Ru Vivien Chen, Wen Ping Chen, Mike Chen, Zhiwei Chen, Jungyeon Cho, Minho Choi , et al. (133 additional authors not shown)

Abstract: We present BISTRO Survey 850 μm dust emission polarisation observations of the L1495A-B10 region of the Taurus molecular cloud, taken at the JCMT. We observe a roughly triangular network of dense filaments. We detect 9 of the dense starless cores embedded within these filaments in polarisation, finding that the plane-of-sky orientation of the core-scale magnetic field lies roughly perpendicular to… ▽ More We present BISTRO Survey 850 μm dust emission polarisation observations of the L1495A-B10 region of the Taurus molecular cloud, taken at the JCMT. We observe a roughly triangular network of dense filaments. We detect 9 of the dense starless cores embedded within these filaments in polarisation, finding that the plane-of-sky orientation of the core-scale magnetic field lies roughly perpendicular to the filaments in almost all cases. We also find that the large-scale magnetic field orientation measured by Planck is not correlated with any of the core or filament structures, except in the case of the lowest-density core. We propose a scenario for early prestellar evolution that is both an extension to, and consistent with, previous models, introducing an additional evolutionary transitional stage between field-dominated and matter-dominated evolution, observed here for the first time. In this scenario, the cloud collapses first to a sheet-like structure. Uniquely, we appear to be seeing this sheet almost face-on. The sheet fragments into filaments, which in turn form cores. However, the material must reach a certain critical density before the evolution changes from being field-dominated to being matter-dominated. We measure the sheet surface density and the magnetic field strength at that transition for the first time and show consistency with an analytical prediction that had previously gone untested for over 50 years (Mestel 1965). △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: 14 pages, 5 figures. ApJ accepted

Study of the effect of turbulent interstellar mediums on young supernova remnant morphology

Authors: Gabriel Rigon, Tsuyoshi Inoue

Abstract: Context. Supernova remnants (SNR) are one of the main sources of galactic cosmic rays acceleration. This acceleration, believed to happen at the blast wave front, leads to an energy loss at the shock front. This results in the apparent proximity between the blast wave and the contact discontinuity. Aims. In this article, we study the effect that turbulent-like density perturbations of the interste… ▽ More Context. Supernova remnants (SNR) are one of the main sources of galactic cosmic rays acceleration. This acceleration, believed to happen at the blast wave front, leads to an energy loss at the shock front. This results in the apparent proximity between the blast wave and the contact discontinuity. Aims. In this article, we study the effect that turbulent-like density perturbations of the interstellar medium (ISM) have on the evolution of young SNRs. We focus on the impact such fluctuations have on the SNRs structure and more precisely on the resulting distance between blast wave and contact discontinuity. Since cosmic rays acceleration is necessary to explain this distance, this study indirectly put into question the cosmic rays acceleration at the blast wave front. Methods. We performed a set of purely hydrodynamic three-dimensional simulations without cosmic ray acceleration in a co-expanding frame. We randomly initialised the density variation of the interstellar medium following a Kolmogorov power law. The resulting ratios of radii between blast wave, contact discontinuity and reverse shock are then compared to the astronomical observations. Results. The addition of density perturbation doesnt significantly change the average ratio of radii. However, the simulations show a higher growth of interfacial instabilities in the presence of a turbulent ISM. The resulting deformation of the contact discontinuity could explain the proximity between contact discontinuity and blast wave. They also explain the plateau in the radial distribution of the line of sight velocity associated with the observations of Tycho. Conclusions. The density perturbation of the ISM should not be neglected in the simulation of young SNR as they have an impact comparable to the cosmic rays acceleration on the SNR structure. △ Less

Submitted 24 December, 2022; originally announced December 2022.

Journal ref: A&A 672, A194 (2023)

JCMT BISTRO Observations: Magnetic Field Morphology of Bubbles Associated with NGC 6334

Authors: Mehrnoosh Tahani, Pierre Bastien, Ray S. Furuya, Kate Pattle, Doug Johnstone, Doris Arzoumanian, Yasuo Doi, Tetsuo Hasegawa, Shu-ichiro Inutsuka, Simon Coudé, Laura Fissel, Michael Chun-Yuan Chen, Frédérick Poidevin, Sarah Sadavoy, Rachel Friesen, Patrick M. Koch, James Di Francesco, Gerald H. Moriarty-Schieven, Zhiwei Chen, Eun Jung Chung, Chakali Eswaraiah, Lapo Fanciullo, Tim Gledhill, Valentin J. M. Le Gouellec, Thiem Hoang , et al. (120 additional authors not shown)

Abstract: We study the HII regions associated with the NGC 6334 molecular cloud observed in the sub-millimeter and taken as part of the B-fields In STar-forming Region Observations (BISTRO) Survey. In particular, we investigate the polarization patterns and magnetic field morphologies associated with these HII regions. Through polarization pattern and pressure calculation analyses, several of these bubbles… ▽ More We study the HII regions associated with the NGC 6334 molecular cloud observed in the sub-millimeter and taken as part of the B-fields In STar-forming Region Observations (BISTRO) Survey. In particular, we investigate the polarization patterns and magnetic field morphologies associated with these HII regions. Through polarization pattern and pressure calculation analyses, several of these bubbles indicate that the gas and magnetic field lines have been pushed away from the bubble, toward an almost tangential (to the bubble) magnetic field morphology. In the densest part of NGC 6334, where the magnetic field morphology is similar to an hourglass, the polarization observations do not exhibit observable impact from HII regions. We detect two nested radial polarization patterns in a bubble to the south of NGC 6334 that correspond to the previously observed bipolar structure in this bubble. Finally, using the results of this study, we present steps (incorporating computer vision; circular Hough Transform) that can be used in future studies to identify bubbles that have physically impacted magnetic field lines. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: Accepted for publication in Astrophysical Journal (ApJ)

Distance determination of molecular clouds in the 1st quadrant of the Galactic plane using deep learning : I. Method and Results

Authors: Shinji Fujita, A. M. Ito, Yusuke Miyamoto, Yasutomo Kawanishi, Kazufumi Torii, Yoshito Shimajiri, Atsushi Nishimura, Kazuki Tokuda, Toshikazu Ohnishi, Hiroyuki Kaneko, Tsuyoshi Inoue, Shunya Takekawa, Mikito Kohno, Shota Ueda, Shimpei Nishimoto, Ryuki Yoneda, Kaoru Nishikawa, Daisuke Yoshida

Abstract: Machine learning has been successfully applied in varied field but whether it is a viable tool for determining the distance to molecular clouds in the Galaxy is an open question. In the Galaxy, the kinematic distance is commonly employed as the distance to a molecular cloud. However, there is a problem in that for the inner Galaxy, two different solutions, the ``Near'' solution, and the ``Far'' so… ▽ More Machine learning has been successfully applied in varied field but whether it is a viable tool for determining the distance to molecular clouds in the Galaxy is an open question. In the Galaxy, the kinematic distance is commonly employed as the distance to a molecular cloud. However, there is a problem in that for the inner Galaxy, two different solutions, the ``Near'' solution, and the ``Far'' solution, can be derived simultaneously. We attempted to construct a two-class (``Near'' or ``Far'') inference model using a Convolutional Neural Network (CNN), a form of deep learning that can capture spatial features generally. In this study, we used the CO dataset toward the 1st quadrant of the Galactic plane obtained with the Nobeyama 45-m radio telescope (l = 62-10 degree, |b| < 1 degree). In the model, we applied the three-dimensional distribution (position-position-velocity) of the 12CO (J=1-0) emissions as the main input. The dataset with ``Near'' or ``Far'' annotation was made from the HII region catalog of the infrared astronomy satellite WISE to train the model. As a result, we could construct a CNN model with a 76% accuracy rate on the training dataset. By using the model, we determined the distance to molecular clouds identified by the CLUMPFIND algorithm. We found that the mass of the molecular clouds with a distance of < 8.15 kpc identified in the 12CO data follows a power-law distribution with an index of about -2.3 in the mass range of M >10^3 Msun. Also, the detailed molecular gas distribution of the Galaxy as seen from the Galactic North pole was determined. △ Less

Submitted 12 December, 2022; originally announced December 2022.

Comments: 29 pages, 12 figures

The JCMT BISTRO-2 Survey: Magnetic Fields of the Massive DR21 Filament

Authors: Tao-Chung Ching, Keping Qiu, Di Li, Zhiyuan Ren, Shih-Ping Lai, David Berry, Kate Pattle, Ray Furuya, Derek Ward-Thompson, Doug Johnstone, Patrick M. Koch, Chang Won Lee, Thiem Hoang, Tetsuo Hasegawa, Woojin Kwon, Pierre Bastien, Chakali Eswaraiah, Jia-Wei Wang, Kyoung Hee Kim, Jihye Hwang, Archana Soam, A-Ran Lyo, Junhao Liu, Valentin J. M. Le Gouellec, Doris Arzoumanian , et al. (132 additional authors not shown)

Abstract: We present 850 $μ$m dust polarization observations of the massive DR21 filament from the B-fields In STar-forming Region Observations (BISTRO) survey, using the POL-2 polarimeter and the SCUBA-2 camera on the James Clerk Maxwell Telescope. We detect ordered magnetic fields perpendicular to the parsec-scale ridge of the DR21 main filament. In the sub-filaments, the magnetic fields are mainly parall… ▽ More We present 850 $μ$m dust polarization observations of the massive DR21 filament from the B-fields In STar-forming Region Observations (BISTRO) survey, using the POL-2 polarimeter and the SCUBA-2 camera on the James Clerk Maxwell Telescope. We detect ordered magnetic fields perpendicular to the parsec-scale ridge of the DR21 main filament. In the sub-filaments, the magnetic fields are mainly parallel to the filamentary structures and smoothly connect to the magnetic fields of the main filament. We compare the POL-2 and Planck dust polarization observations to study the magnetic field structures of the DR21 filament on 0.1--10 pc scales. The magnetic fields revealed in the Planck data are well aligned with those of the POL-2 data, indicating a smooth variation of magnetic fields from large to small scales. The plane-of-sky magnetic field strengths derived from angular dispersion functions of dust polarization are 0.6--1.0 mG in the DR21 filament and $\sim$ 0.1 mG in the surrounding ambient gas. The mass-to-flux ratios are found to be magnetically supercritical in the filament and slightly subcritical to nearly critical in the ambient gas. The alignment between column density structures and magnetic fields changes from random alignment in the low-density ambient gas probed by Planck to mostly perpendicular in the high-density main filament probed by JCMT. The magnetic field structures of the DR21 filament are in agreement with MHD simulations of a strongly magnetized medium, suggesting that magnetic fields play an important role in shaping the DR21 main filament and sub-filaments. △ Less

Submitted 4 December, 2022; originally announced December 2022.

Comments: 26 pages, 13 figures, ApJ accepted

The Effect of Shock Wave Duration on Star Formation and the Initial Condition of Massive Cluster Formation

Authors: Daisei Abe, Tsuyoshi Inoue, Rei Enokiya, Yasuo Fukui

Abstract: Stars are born in dense molecular filaments irrespective of their mass. Compression of the ISM by shocks cause filament formation in molecular clouds. Observations show that a massive star cluster formation occurs where the peak of gas column density in a cloud exceeds 10^23 cm^-2. In this study, we investigate the effect of the shock-compressed layer duration on filament/star formation and how th… ▽ More Stars are born in dense molecular filaments irrespective of their mass. Compression of the ISM by shocks cause filament formation in molecular clouds. Observations show that a massive star cluster formation occurs where the peak of gas column density in a cloud exceeds 10^23 cm^-2. In this study, we investigate the effect of the shock-compressed layer duration on filament/star formation and how the initial conditions of massive star formation are realized by performing three-dimensional (3D) isothermal magnetohydrodynamics (MHD) simulations with {gas inflow duration from the boundaries (i.e., shock wave duration)} as a controlling parameter. Filaments formed behind the shock expand after the duration time for short shock duration models, whereas long duration models lead to star formation by forming massive supercritical filaments. Moreover, when the shock duration is longer than two postshock free-fall times, the peak column density of the compressed layer exceeds 10^23 cm^-2, and {the gravitational collapse of the layer causes that} the number of OB stars expected to be formed in the shock-compressed layer reaches the order of ten (i.e., massive cluster formation). △ Less

Submitted 3 November, 2022; originally announced November 2022.

Comments: 16 pages, 8 figures

The JCMT BISTRO Survey: A Spiral Magnetic Field in a Hub-filament Structure, Monoceros R2

Authors: Jihye Hwang, Jongsoo Kim, Kate Pattle, Chang Won Lee, Patrick M. Koch, Doug Johnstone, Kohji Tomisaka, Anthony Whitworth, Ray S. Furuya, Ji-hyun Kang, A-Ran Lyo, Eun Jung Chung, Doris Arzoumanian, Geumsook Park, Woojin Kwon, Shinyoung Kim, Motohide Tamura, Jungmi Kwon, Archana Soam, Ilseung Han, Thiem Hoang, Kyoung Hee Kim, Takashi Onaka, Eswaraiah Chakali, Derek Ward-Thompson , et al. (135 additional authors not shown)

Abstract: We present and analyze observations of polarized dust emission at 850 $μ$m towards the central 1 pc $\times$ 1 pc hub-filament structure of Monoceros R2 (Mon R2). The data are obtained with SCUBA-2/POL-2 on the James Clerk Maxwell Telescope (JCMT) as part of the BISTRO (B-fields in Star-forming Region Observations) survey. The orientations of the magnetic field follow the spiral structure of Mon R… ▽ More We present and analyze observations of polarized dust emission at 850 $μ$m towards the central 1 pc $\times$ 1 pc hub-filament structure of Monoceros R2 (Mon R2). The data are obtained with SCUBA-2/POL-2 on the James Clerk Maxwell Telescope (JCMT) as part of the BISTRO (B-fields in Star-forming Region Observations) survey. The orientations of the magnetic field follow the spiral structure of Mon R2, which are well-described by an axisymmetric magnetic field model. We estimate the turbulent component of the magnetic field using the angle difference between our observations and the best-fit model of the underlying large-scale mean magnetic field. This estimate is used to calculate the magnetic field strength using the Davis-Chandrasekhar-Fermi method, for which we also obtain the distribution of volume density and velocity dispersion using a column density map derived from $Herschel$ data and the C$^{18}$O ($J$ = 3-2) data taken with HARP on the JCMT, respectively. We make maps of magnetic field strengths and mass-to-flux ratios, finding that magnetic field strengths vary from 0.02 to 3.64 mG with a mean value of 1.0 $\pm$ 0.06 mG, and the mean critical mass-to-flux ratio is 0.47 $\pm$ 0.02. Additionally, the mean Alfvén Mach number is 0.35 $\pm$ 0.01. This suggests that in Mon R2, magnetic fields provide resistance against large-scale gravitational collapse, and magnetic pressure exceeds turbulent pressure. We also investigate the properties of each filament in Mon R2. Most of the filaments are aligned along the magnetic field direction and are magnetically sub-critical. △ Less

Submitted 13 December, 2022; v1 submitted 12 October, 2022; originally announced October 2022.

Comments: This paper is accepted to the ApJ

Particle acceleration controlled by ambient density in the southwestern rim of RCW 86

Authors: Hiromasa Suzuki, Satoru Katsuda, Takaaki Tanaka, Nobuaki Sasaki, Tsuyoshi Inoue, Federico Fraschetti

Abstract: Particle acceleration physics at supernova remnant (SNR) shocks is one of the most intriguing problems in astrophysics. SNR RCW~86 provides a suitable environment for understanding the particle acceleration physics because one can extract the information of both accelerated particles and acceleration environment at the same regions through the bright X-ray emission. In this work, we study X-ray pr… ▽ More Particle acceleration physics at supernova remnant (SNR) shocks is one of the most intriguing problems in astrophysics. SNR RCW~86 provides a suitable environment for understanding the particle acceleration physics because one can extract the information of both accelerated particles and acceleration environment at the same regions through the bright X-ray emission. In this work, we study X-ray proper motions and spectral properties of the southwestern region of RCW~86. The proper motion velocities are found to be $\sim 300$--2000~km~s$^{-1}$ at a distance of 2.8~kpc. We find two inward-moving filaments, which are more likely reflected shocks rather than reverse shocks. Based on the X-ray spectroscopy, we evaluate thermal parameters such as the ambient density and temperature, and non-thermal parameters such as the power-law flux and index. From the flux decrease in time of several non-thermal filaments, we estimate the magnetic field amplitudes to be $\sim 30$--100~$μ$G. Gathering the physical parameters, we then investigate parameter correlations. We find that the synchrotron emission from thermal-dominated filaments is correlated with the ambient density $n_{\rm e}$ as $\text{(power-law flux)} \propto n_{\rm e}^{1.0 \pm 0.2}$ and $\text{(power-law index)} \propto n_{\rm e}^{0.38 \pm 0.10}$, not or only weakly with the shock velocity and shock obliquity. As an interpretation, we propose a shock-cloud interaction scenario, where locally enhanced magnetic turbulence levels have a great influence on local acceleration conditions. △ Less

Submitted 29 August, 2022; originally announced August 2022.

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

On the Potential of Faraday Tomography to Identify Shock Structures in Supernova Remnants

Authors: Shinsuke Ideguchi, Tsuyoshi Inoue, Takuya Akahori, Keitaro Takahashi

Abstract: Knowledge about the magnetic fields in supernova remnants (SNRs) is of paramount importance for constraining Galactic cosmic ray acceleration models. It could also indirectly provide information on the interstellar magnetic fields. In this paper, we predict the Faraday dispersion functions (FDFs) of SNRs for the first time. For this study, we use the results of three dimensional (3D) ideal magneto… ▽ More Knowledge about the magnetic fields in supernova remnants (SNRs) is of paramount importance for constraining Galactic cosmic ray acceleration models. It could also indirectly provide information on the interstellar magnetic fields. In this paper, we predict the Faraday dispersion functions (FDFs) of SNRs for the first time. For this study, we use the results of three dimensional (3D) ideal magnetohydrodynamic (MHD) simulations of SNRs expanding into a weak, regular magnetic field. We present the intrinsic FDFs of the shocked region of SNRs for different viewing angles. We find that the FDFs are generally Faraday complex, which implies that conventional rotation measure study is not sufficient to obtain the information on the magnetic fields in the shocked region and Faraday tomography is necessary. We also show that the FDF allows to derive the physical-depth distribution of polarization intensity when the line of sight is parallel to the initial magnetic field orientation. Furthermore, we demonstrate that the location of contact discontinuity can be identified from the radial profile of the width of the FDF with the accuracy of 0.1-0.2 pc. △ Less

Submitted 22 August, 2022; originally announced August 2022.

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 513, Issue 3, July 2022, Pages 3289-3301

Relativistic distortions in galaxy density-ellipticity correlations: gravitational redshift and peculiar velocity effects

Authors: Shohei Saga, Teppei Okumura, Atsushi Taruya, Takuya Inoue

Abstract: We study relativistic effects, arising from the light propagation in an inhomogeneous universe. We particularly investigate the effects imprinted in a cross-correlation function between galaxy positions and intrinsic galaxy shapes (GI correlation). Considering the Doppler and gravitational redshift effects as major relativistic effects, we present an analytical model of the GI correlation function… ▽ More We study relativistic effects, arising from the light propagation in an inhomogeneous universe. We particularly investigate the effects imprinted in a cross-correlation function between galaxy positions and intrinsic galaxy shapes (GI correlation). Considering the Doppler and gravitational redshift effects as major relativistic effects, we present an analytical model of the GI correlation function, from which we find that the relativistic effects induce non-vanishing odd multipole anisotropies. Focusing particularly on the dipole anisotropy, we show that the Doppler effect dominates at large scales, while the gravitational redshift effect originated from the halo potential dominates at the scales below $10$-$30\, {\rm Mpc}/h$, with the amplitude of the dipole GI correlation being positive over all the scales. Also, we newly derive the covariance matrix for the modelled GI dipole. Taking into account the full covariance, we estimate the signal-to-noise ratio and show that the GI dipole induced by the relativistic effects is detectable in future large-volume galaxy surveys. We discuss how the measurement of dipole GI correlation could be helpful to detect relativistic effects in combination with the conventional galaxy-galaxy cross correlation. △ Less

Submitted 15 August, 2023; v1 submitted 7 July, 2022; originally announced July 2022.

Comments: 16 pages, 7 figures

Report number: YITP-22-69

From Bubbles and Filaments to Cores and Disks: Gas Gathering and Growth of Structure Leading to the Formation of Stellar Systems

Authors: Jaime E. Pineda, Doris Arzoumanian, Philippe André, Rachel K. Friesen, Annie Zavagno, Seamus D. Clarke, Tsuyoshi Inoue, Che-Yu Chen, Yueh-Ning Lee, Juan D. Soler, Michael Kuffmeier

Abstract: The study of the development of structures on multiple scales in the cold interstellar medium has experienced rapid expansion in the past decade, on both the observational and the theoretical front. Spectral line studies at (sub-)millimeter wavelengths over a wide range of physical scales have provided unique probes of the kinematics of dense gas in star-forming regions, and have been complemented… ▽ More The study of the development of structures on multiple scales in the cold interstellar medium has experienced rapid expansion in the past decade, on both the observational and the theoretical front. Spectral line studies at (sub-)millimeter wavelengths over a wide range of physical scales have provided unique probes of the kinematics of dense gas in star-forming regions, and have been complemented by extensive, high dynamic range dust continuum surveys of the column density structure of molecular cloud complexes, while dust polarization maps have highlighted the role of magnetic fields. This has been accompanied by increasingly sophisticated numerical simulations including new physics (e.g., supernova driving, cosmic rays, non-ideal magneto-hydrodynamics, radiation pressure) and new techniques such as zoom-in simulations allowing multi-scale studies. Taken together, these new data have emphasized the anisotropic growth of dense structures on all scales, from giant ISM bubbles driven by stellar feedback on $\sim$50-100 pc scales through parsec-scale molecular filaments down to $<$0.1 pc dense cores and $<$1000 au protostellar disks. Combining observations and theory, we present a coherent picture for the formation and evolution of these structures and synthesize a comprehensive physical scenario for the initial conditions and early stages of star and disk formation. △ Less

Submitted 23 March, 2023; v1 submitted 8 May, 2022; originally announced May 2022.

Comments: 42 pages, 20 figures. To appear in Protostars and Planets VII, Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura. Small changes to match proofs and the updated format

An ALMA study of the massive molecular clump N159W-North in the Large Magellanic Cloud: A possible gas flow penetrating one of the most massive protocluster systems in the Local Group

Authors: Kazuki Tokuda, Taisei Minami, Yasuo Fukui, Tsuyoshi Inoue, Takeru Nishioka, Kisetsu Tsuge, Sarolta Zahorecz, Hidetoshi Sano, Ayu Konishi, C. -H. Rosie Chen, Marta Sewiło, Suzanne C. Madden, Omnarayani Nayak, Kazuya Saigo, Atsushi Nishimura, Kei E. I. Tanaka, Tsuyoshi Sawada, Remy Indebetouw, Kengo Tachihara, Akiko Kawamura, Toshikazu Onishi

Abstract: Massive dense clumps in the Large Magellanic Cloud can be an important laboratory to explore the formation of populous clusters. We report multiscale ALMA observations of the N159W-North clump, which is the most CO-intense region in the galaxy. High-resolution CO isotope and 1.3 mm continuum observations with an angular resolution of $\sim$0."25($\sim$0.07 pc) revealed more than five protostellar… ▽ More Massive dense clumps in the Large Magellanic Cloud can be an important laboratory to explore the formation of populous clusters. We report multiscale ALMA observations of the N159W-North clump, which is the most CO-intense region in the galaxy. High-resolution CO isotope and 1.3 mm continuum observations with an angular resolution of $\sim$0."25($\sim$0.07 pc) revealed more than five protostellar sources with CO outflows within the main ridge clump. One of the thermal continuum sources, MMS-2, shows especially massive/dense nature whose total H$_2$ mass and peak column density are $\sim$10$^{4}$ $M_{\odot}$ and $\sim$10$^{24}$ cm$^{-2}$, respectively, and harbors massive ($\sim$100 $M_{\odot}$) starless core candidates identified as its internal substructures. The main ridge containing this source can be categorized as one of the most massive protocluster systems in the Local Group. The CO high-resolution observations found several distinct filamentary clouds extending southward from the star-forming spots. The CO (1-0) data set with a larger field of view reveals a conical-shaped, $\sim$30 pc long complex extending toward the northern direction. These features indicate that a large-scale gas compression event may have produced the massive star-forming complex. Based on the striking similarity between the N159W-North complex and the previously reported other two high-mass star-forming clouds in the nearby regions, we propose a $"$teardrops inflow model$"$ that explains the synchronized, extreme star formation across $>$50 pc, including one of the most massive protocluster clumps in the Local Group. △ Less

Submitted 29 April, 2022; originally announced May 2022.

Comments: 25 pages, 14 figures, 6 tables, Accepted for publication in ApJ

Nature of supersonic turbulence and density distribution function in the multiphase interstellar medium

Authors: Masato I. N. Kobayashi, Tsuyoshi Inoue, Kengo Tomida, Kazunari Iwasaki, Hiroki Nakatsugawa

Abstract: Supersonic flows in the interstellar medium (ISM) are believed to be a key driver of the molecular cloud formation and evolution. Among molecular clouds' properties, the ratio between the solenoidal and compressive modes of turbulence plays important roles in determining the star formation efficiency. We use numerical simulations of supersonic converging flows of the warm neutral medium (WNM) reso… ▽ More Supersonic flows in the interstellar medium (ISM) are believed to be a key driver of the molecular cloud formation and evolution. Among molecular clouds' properties, the ratio between the solenoidal and compressive modes of turbulence plays important roles in determining the star formation efficiency. We use numerical simulations of supersonic converging flows of the warm neutral medium (WNM) resolving the thermal instability to calculate the early phase of molecular cloud formation, and investigate the turbulence structure and the density probability distribution function (density PDF) of the multiphase ISM. We find that both the solenoidal and compressive modes have their power spectrum similar to the Kolmogorov spectrum. The solenoidal (compressive) modes account for >~80% (<~20%) of the total turbulence power. When we consider both the cold neutral medium (CNM) and the thermally unstable neutral medium (UNM) up to T <~ 400 K, the density PDF follows the log-normal distribution whose width sigma_s is well explained by the known relation from the isothermal turbulence as sigma_s = ln(1 + b^2 * M^2) (where b is the parameter representing the turbulence mode ratio and M is the turbulent Mach number). The density PDF of the CNM component alone (T <= 50 K), however, exhibits a narrower sigma_s by a factor of ~ 2. These results suggest that observational estimations of b based on the CNM density PDF requires the internal turbulence within each CNM clump but not the inter-clump relative velocity, the latter of which is instead powered by the WNM/UNM turbulence. △ Less

Submitted 1 March, 2022; originally announced March 2022.

Comments: 19 pages, 21 figures, 3 tables; accepted for publication in the ApJ, included the referee comments

HOLISMOKES. VIII. High-redshift, strong-lens search in the Hyper Suprime-Cam Subaru Strategic Program

Authors: Yiping Shu, Raoul Cañameras, Stefan Schuldt, Sherry H. Suyu, Stefan Taubenberger, Kaiki Taro Inoue, Anton T. Jaelani

Abstract: We carry out a search for strong-lens systems containing high-redshift lens galaxies with the goal of extending strong-lensing-assisted galaxy evolutionary studies to earlier cosmic time. Two strong-lens classifiers are constructed from a deep residual network and trained with datasets of different lens-redshift and brightness distributions. We classify a sample of 5,356,628 pre-selected objects f… ▽ More We carry out a search for strong-lens systems containing high-redshift lens galaxies with the goal of extending strong-lensing-assisted galaxy evolutionary studies to earlier cosmic time. Two strong-lens classifiers are constructed from a deep residual network and trained with datasets of different lens-redshift and brightness distributions. We classify a sample of 5,356,628 pre-selected objects from the Wide-layer fields in the second public data release of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) by applying the two classifiers to their HSC $gri$-filter cutouts. Cutting off at thresholds that correspond to a false positive rate of $10^{-3}$ on our test set, the two classifiers identify 5,468 and 6,119 strong-lens candidates. Visually inspecting the cutouts of those candidates results in 735 grade-A or B strong-lens candidates in total, of which 277 candidates are discovered for the first time. This is the single largest set of galaxy-scale strong-lens candidates discovered with HSC data to date, and nearly half of it (331/735) contains lens galaxies with photometric redshifts above 0.6. Our discoveries will serve as a valuable target list for ongoing and scheduled spectroscopic surveys such as the Dark Energy Spectroscopic Instrument, the Subaru Prime Focus Spectrograph project, and the Maunakea Spectroscopic Explorer. △ Less

Submitted 25 May, 2022; v1 submitted 26 January, 2022; originally announced January 2022.

Comments: 15 pages, 8 figures + appendices, published on the A&A

Report number: 2022, A&A, 662, A4

Journal ref: A&A 662, A4 (2022)

X-ray Line Diagnostics of Ion Temperature at Cosmic-Ray Accelerating Collisionless Shocks

Authors: Jiro Shimoda, Yutaka Ohira, Aya Bamba, Yukikatsu Terada, Ryo Yamazaki, Tsuyoshi Inoue, Shuta J. Tanaka

Abstract: A novel collisionless shock jump condition is suggested by modeling the entropy production at the shock transition region. We also calculate downstream developments of the atomic ionization balance and the ion temperature relaxation in supernova remnants (SNRs). The injection process and subsequent acceleration of cosmic-rays (CRs) in the SNR shocks are closely related to the formation process of… ▽ More A novel collisionless shock jump condition is suggested by modeling the entropy production at the shock transition region. We also calculate downstream developments of the atomic ionization balance and the ion temperature relaxation in supernova remnants (SNRs). The injection process and subsequent acceleration of cosmic-rays (CRs) in the SNR shocks are closely related to the formation process of the collisionless shocks. The formation of the shock is caused by wave-particle interactions. Since the wave-particle interactions result in energy exchanges between electromagnetic fields and charged particles, the randomization of particles associated with the shock transition may occur with the rate given by the scalar product of the electric field and current. We find that order-of-magnitude estimates of the randomization with reasonable strength of the electromagnetic fields in the SNR constrain the amount of the CR nuclei and ion temperatures. The constrained amount of the CR nuclei can be sufficient to explain the Galactic CRs. The ion temperature becomes significantly lower than in the case of no CRs. To distinguish the case without CRs, we perform synthetic observations of atomic line emissions from the downstream region of the SNR RCW~86. Future observations by {\it XRISM} and {\it Athena} can distinguish whether the SNR shock accelerates the CRs or not from the ion temperatures. △ Less

Submitted 16 June, 2022; v1 submitted 19 January, 2022; originally announced January 2022.

Comments: 18 pages, 14 figures, 3 tables, accepted for publication in PASJ

B-fields in Star-Forming Region Observations (BISTRO): Magnetic Fields in the Filamentary Structures of Serpens Main

Authors: Woojin Kwon, Kate Pattle, Sarah Sadavoy, Charles L. H. Hull, Doug Johnstone, Derek Ward-Thompson, James Di Francesco, Patrick M. Koch, Ray Furuya, Yasuo Doi, Valentin J. M. Le Gouellec, Jihye Hwang, A-Ran Lyo, Archana Soam, Xindi Tang, Thiem Hoang, Florian Kirchschlager, Chakali Eswaraiah, Lapo Fanciullo, Kyoung Hee Kim, Takashi Onaka, Vera Könyves, Ji-hyun Kang, Chang Won Lee, Motohide Tamura , et al. (127 additional authors not shown)

Abstract: We present 850 $μ$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filament… ▽ More We present 850 $μ$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filaments with different physical properties such as density and star formation activity. Using the histogram of relative orientation (HRO) technique, we find that magnetic fields are parallel to filaments in less dense filamentary structures where $N_{H_2} < 0.93\times 10^{22}$ cm$^{-2}$ (magnetic fields perpendicular to density gradients), while being perpendicular to filaments (magnetic fields parallel to density gradients) in dense filamentary structures with star formation activity. Moreover, applying the HRO technique to denser core regions, we find that magnetic field orientations change to become perpendicular to density gradients again at $N_{H_2} \approx 4.6 \times 10^{22}$ cm$^{-2}$. This can be interpreted as a signature of core formation. At $N_{H_2} \approx 16 \times 10^{22}$ cm$^{-2}$ magnetic fields change back to being parallel to density gradients once again, which can be understood to be due to magnetic fields being dragged in by infalling material. In addition, we estimate the magnetic field strengths of the filaments ($B_{POS} = 60-300~μ$G)) using the Davis-Chandrasekhar-Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities. △ Less

Submitted 13 January, 2022; originally announced January 2022.

Comments: 18 pages, accepted for publication in ApJ

The JCMT BISTRO Survey: An 850/450$μ$m Polarization Study of NGC 2071IR in OrionB

Authors: A-Ran Lyo, Jongsoo Kim, Sarah Sadavoy, Doug Johnstone, David Berry, Kate Pattle, Woojin Kwon, Pierre Bastien, Takashi Onaka, James Di Francesco, Ji-Hyun Kang, Ray Furuya, Charles L. H. Hull, Motohide Tamura, Patrick M. Koch, Derek Ward-Thompson, Tetsuo Hasegawa, Thiem Hoang, Doris Arzoumanian, Chang Won Lee, Chin-Fei Lee, Do-Young Byun, Florian Kirchschlager, Yasuo Doi, Kee-Tae Kim , et al. (121 additional authors not shown)

Abstract: We present the results of simultaneous 450 $μ$m and 850 $μ$m polarization observations toward the massive star forming region NGC 2071IR, a target of the BISTRO (B-fields in Star-Forming Region Observations) Survey, using the POL-2 polarimeter and SCUBA-2 camera mounted on the James Clerk Maxwell Telescope. We find a pinched magnetic field morphology in the central dense core region, which could b… ▽ More We present the results of simultaneous 450 $μ$m and 850 $μ$m polarization observations toward the massive star forming region NGC 2071IR, a target of the BISTRO (B-fields in Star-Forming Region Observations) Survey, using the POL-2 polarimeter and SCUBA-2 camera mounted on the James Clerk Maxwell Telescope. We find a pinched magnetic field morphology in the central dense core region, which could be due to a rotating toroidal disk-like structure and a bipolar outflow originating from the central young stellar object, IRS 3. Using the modified Davis-Chandrasekhar-Fermi method, we obtain a plane-of-sky magnetic field strength of 563$\pm$421 $μ$G in the central $\sim$0.12 pc region from 850 $μ$m polarization data. The corresponding magnetic energy density of 2.04$\times$10$^{-8}$ erg cm$^{-3}$ is comparable to the turbulent and gravitational energy densities in the region. We find that the magnetic field direction is very well aligned with the whole of the IRS 3 bipolar outflow structure. We find that the median value of polarization fractions, 3.0 \%, at 450 $μ$m in the central 3 arcminute region, which is larger than the median value of 1.2 \% at 850 $μ$m. The trend could be due to the better alignment of warmer dust in the strong radiation environment. We also find that polarization fractions decrease with intensity at both wavelengths, with slopes, determined by fitting a Rician noise model, of $0.59 \pm 0.03$ at 450 $μ$m and $0.36 \pm 0.04$ at 850 $μ$m, respectively. We think that the shallow slope at 850 $μ$m is due to grain alignment at the center being assisted by strong radiation from the central young stellar objects. △ Less

Submitted 28 September, 2021; originally announced September 2021.

Comments: 17 pages, 11 figures, 4 tables

Journal ref: The Astrophysical Journal, 2021, 918, 85

ALMA Measurement of 10 kpc-scale Lensing Power Spectra towards the Lensed Quasar MG J0414+0534

Authors: Kaiki Taro Inoue, Takeo Minezaki, Satoki Matsushita, Koichiro Nakanishi

Abstract: The lensing power spectra for gravitational potential, astrometric shift, and convergence perturbations are powerful probes to investigate dark matter structures on small scales. We report the first lower and upper bounds of these lensing power spectra on angular scale $\sim 1$arcsec towards the anomalous quadruply lensed quasar MG$\,$J0414+0534 at a redshift $z=2.639$. To obtain the spectra, we c… ▽ More The lensing power spectra for gravitational potential, astrometric shift, and convergence perturbations are powerful probes to investigate dark matter structures on small scales. We report the first lower and upper bounds of these lensing power spectra on angular scale $\sim 1$arcsec towards the anomalous quadruply lensed quasar MG$\,$J0414+0534 at a redshift $z=2.639$. To obtain the spectra, we conducted observations of MG$\,$J0414+0534 using the Atacama Large Millimeter/submillimeter Array (ALMA) with high angular resolution (0.02 arcsec to 0.05 arcsec). We developed a new partially non-parametric method in which Fourier coefficients of potential perturbation are adjusted to minimize the difference between linear combinations of weighted mean de-lensed images. Using positions of radio jet components, extended dust emission on scales $>1\,$kpc, and mid-infrared flux ratios, the range of measured convergence, astrometric shift, and potential powers at an angular scale of $\sim 1.1$arcsec (corresponding to an angular wave number of $l=1.2\times 10^6$ or $\sim 9\,$kpc in the primary lens plane) within $1\,σ$ are $\varDelta_κ=0.021-0.028$, $\varDelta_α=7-9\,$mas, and $\varDelta_ψ=1.2-1.6\,$$\textrm{mas}^2$, respectively. Our result is consistent with the predicted abundance of halos in the line of sight and subhalos in cold dark matter models. Our partially non-parametric lens models suggest a presence of a clump in the vicinity of object Y, a possible dusty dwarf galaxy and some small clumps in the vicinity of other lensed quadruple images. Although much fainter than the previous report, we detected weak continuum emission possibly from object Y with a peak flux of $\sim 100\,μ\textrm{Jy}\, \textrm{beam} ^{-1}$ at the $\sim 4\,σ$ level. △ Less

Submitted 26 July, 2023; v1 submitted 2 September, 2021; originally announced September 2021.

Comments: 54 pages, 35 figures, version accepted in ApJ. Analyses of models with a possible dwarf galaxy, object Y and that of chi^2 fit in the visibility plane were added. A new section 7 about CDM consistency was added. The angular scale of the smallest angular bin in Figure 24 was corrected

Journal ref: The Astrophysical Journal (2023) 954, 197 (33pp) (Open Access)

Direct Numerical Simulations of Cosmic-ray Acceleration at Dense Circumstellar Medium: Magnetic Field Amplification by Bell Instability and Maximum Energy

Authors: Tsuyoshi Inoue, Alexandre Marcowith, Gwenael Giacinti, Allard Jan van Marle, Shogo Nishino

Abstract: Galactic cosmic rays are believed to be accelerated at supernova remnants. However, whether supernova remnants can be Pevatrons is still very unclear. In this work we argue that PeV cosmic rays can be accelerated during the early phase of a supernova blast wave expansion in dense red supergiant winds. We solve in spherical geometry a system combining a diffusive-convection equation which treats co… ▽ More Galactic cosmic rays are believed to be accelerated at supernova remnants. However, whether supernova remnants can be Pevatrons is still very unclear. In this work we argue that PeV cosmic rays can be accelerated during the early phase of a supernova blast wave expansion in dense red supergiant winds. We solve in spherical geometry a system combining a diffusive-convection equation which treats cosmic-ray dynamics coupled to magnetohydrodynamics to follow gas dynamics. The fast shock expanding in a dense ionized wind is able to trigger the fast non-resonant streaming instability over day timescales, and energizes cosmic-rays even under the effect of p-p losses. We find that such environments make the blast wave a Pevatron, although the maximum energy depends on various parameters such as the injection rate and mass-loss rate of the winds. Multi-PeV energies can be reached if the progenitor mass loss rates are of the order of $10^{-3}$ Msun yr$^{-1}$. It has been recently invoked that, prior to the explosion, hydrogen rich massive stars can produce enhanced mass loss rates. These enhanced rates would then favor the production of a Pevatron phase in early times after the shock breakout. △ Less

Submitted 30 August, 2021; originally announced August 2021.

Comments: 15 pages, 21 figures, accepted by the Astrophysical Journal

HOLISMOKES. VI. New galaxy-scale strong lens candidates from the HSC-SSP imaging survey

Authors: R. Canameras, S. Schuldt, Y. Shu, S. H. Suyu, S. Taubenberger, T. Meinhardt, L. Leal-Taixé, D. C. -Y. Chao, K. T. Inoue, A. T. Jaelani, A. More

Abstract: We have carried out a systematic search for galaxy-scale strong lenses in multiband imaging from the Hyper Suprime-Cam (HSC) survey. Our automated pipeline, based on realistic strong-lens simulations, deep neural network classification, and visual inspection, is aimed at efficiently selecting systems with wide image separations (Einstein radii ~1.0-3.0"), intermediate redshift lenses (z ~ 0.4-0.7)… ▽ More We have carried out a systematic search for galaxy-scale strong lenses in multiband imaging from the Hyper Suprime-Cam (HSC) survey. Our automated pipeline, based on realistic strong-lens simulations, deep neural network classification, and visual inspection, is aimed at efficiently selecting systems with wide image separations (Einstein radii ~1.0-3.0"), intermediate redshift lenses (z ~ 0.4-0.7), and bright arcs for galaxy evolution and cosmology. We classified gri images of all 62.5 million galaxies in HSC Wide with i-band Kron radius >0.8" to avoid strict pre-selections and to prepare for the upcoming era of deep, wide-scale imaging surveys with Euclid and Rubin Observatory. We obtained 206 newly-discovered candidates classified as definite or probable lenses with either spatially-resolved multiple images or extended, distorted arcs. In addition, we found 88 high-quality candidates that were assigned lower confidence in previous HSC searches, and we recovered 173 known systems in the literature. These results demonstrate that, aided by limited human input, deep learning pipelines with false positive rates as low as ~0.01% can be very powerful tools for identifying the rare strong lenses from large catalogs, and can also largely extend the samples found by traditional algorithms. We provide a ranked list of candidates for future spectroscopic confirmation. △ Less

Submitted 7 September, 2021; v1 submitted 16 July, 2021; originally announced July 2021.

Comments: 5 pages and 4 figures (plus appendix), version published in A&A

Journal ref: A&A 653, L6 (2021)

The eROSITA Final Equatorial-Depth Survey (eFEDS): A multiwavelength view of WISE mid-infrared galaxies/active galactic nuclei

Authors: Yoshiki Toba, Teng Liu, Tanya Urrutia, Mara Salvato, Junyao Li, Yoshihiro Ueda, Marcella Brusa, Naomichi Yutani, Keiichi Wada, Atsushi J. Nishizawa, Johannes Buchner, Tohru Nagao, Andrea Merloni, Masayuki Akiyama, Riccardo Arcodia, Bau-Ching Hsieh, Kohei Ichikawa, Masatoshi Imanishi, Kaiki T. Inoue, Toshihiro Kawaguchi, Georg Lamer, Kirpal Nandra, John D. Silverman, Yuichi Terashima

Abstract: We investigate the physical properties--such as the stellar mass, SFR, IR luminosity, X-ray luminosity, and hydrogen column density--of MIR galaxies and AGN at $z < 4$ in the 140 deg$^2$ field observed by SRG/eROSITA through the eFEDS survey. By cross-matching the WISE 22 $μ$m (W4)-detected sample and the eFEDS X-ray point-source catalog, we find that 692 extragalactic objects are detected by eROS… ▽ More We investigate the physical properties--such as the stellar mass, SFR, IR luminosity, X-ray luminosity, and hydrogen column density--of MIR galaxies and AGN at $z < 4$ in the 140 deg$^2$ field observed by SRG/eROSITA through the eFEDS survey. By cross-matching the WISE 22 $μ$m (W4)-detected sample and the eFEDS X-ray point-source catalog, we find that 692 extragalactic objects are detected by eROSITA. We have compiled a multiwavelength dataset. We have also performed (i) an X-ray spectral analysis, (ii) SED fitting using X-CIGALE, (iii) 2D image-decomposition analysis using Subaru HSC images, and (iv) optical spectral fitting with QSFit to investigate the AGN and host-galaxy properties. For 7,088 WISE W4 objects that are undetected by eROSITA, we have performed an X-ray stacking analysis to examine the typical physical properties of these X-ray faint and/or probably obscured objects. We find that (i) 82% of the eFEDS-W4 sources are classified as X-ray AGN with $\log\,L_{\rm X} >$ 42 erg s$^{-1}$; (ii) 67% and 24% of the objects have $\log\,(L_{\rm IR}/L_{\odot}) > 12$ and 13, respectively; (iii) the relationship between $L_{\rm X}$ and the 6 $μ$m luminosity is consistent with that reported in previous works; and (iv) the relationship between the Eddington ratio and $N_{\rm H}$ for the eFEDS-W4 sample and a comparison with a model prediction from a galaxy-merger simulation indicates that approximately 5% of the eFEDS-W4 sources in our sample are likely to be in an AGN-feedback phase, in which strong radiation pressure from the AGN blows out the surrounding material from the nuclear region. Thanks to the wide area coverage of eFEDS, we have been able to constrain the ranges of the physical properties of the WISE W4 sample of AGNs at $z < 4$, providing a benchmark for forthcoming studies on a complete census of MIR galaxies selected from the full-depth eROSITA all-sky survey. △ Less

Submitted 24 August, 2021; v1 submitted 28 June, 2021; originally announced June 2021.

Comments: 18 pages, 19 figures, and 3 tables, accepted to appear on A&A, Special Issue: The Early Data Release of eROSITA and Mikhail Pavlinsky ART-XC on the SRG Mission

Journal ref: A&A 661, A15 (2022)

Associated molecular and atomic clouds with X-ray shell of superbubble 30 Doradus C in the LMC

Authors: Y. Yamane, H. Sano, M. D. Filipovic, K. Tokuda, K. Fujii, Y. Babazaki, I. Mitsuishi, T. Inoue, F. Aharonian, T. Inaba, S. Inutsuka, N. Maxted, N. Mizuno, T. Onishi, G. Rowell, K. Tsuge, F. Voisin, S. Yoshiike, T. Fukuda, A. Kawamura, A. Bamba, K. Tachihara, Y. Fukui

Abstract: 30 Doradus C is a superbubble which emits the brightest nonthermal X- and TeV gamma-rays in the Local Group. In order to explore detailed connection between the high energy radiation and the interstellar medium, we have carried out new CO and HI observations using the Atacama Large Millimeter$/$Submillimeter Array (ALMA), Atacama Submillimeter Telescope Experiment, and the Australia Telescope Comp… ▽ More 30 Doradus C is a superbubble which emits the brightest nonthermal X- and TeV gamma-rays in the Local Group. In order to explore detailed connection between the high energy radiation and the interstellar medium, we have carried out new CO and HI observations using the Atacama Large Millimeter$/$Submillimeter Array (ALMA), Atacama Submillimeter Telescope Experiment, and the Australia Telescope Compact Array with resolutions of up to 3 pc. The ALMA data of $^{12}$CO($J$ = 1-0) emission revealed 23 molecular clouds with the typical diameters of $\sim$6-12 pc and masses of $\sim$600-10000 $M_{\odot}$. The comparison with the X-rays of $XMM$-$Newton$ at $\sim$3 pc resolution shows that X-rays are enhanced toward these clouds. The CO data were combined with the HI to estimate the total interstellar protons. Comparison of the interstellar proton column density and the X-rays revealed that the X-rays are enhanced with the total proton. These are most likely due to the shock-cloud interaction modeled by the magnetohydrodynamical simulations (Inoue et al. 2012, ApJ, 744, 71). Further, we note a trend that the X-ray photon index varies with distance from the center of the high-mass star cluster, suggesting that the cosmic-ray electrons are accelerated by one or multiple supernovae in the cluster. Based on these results we discuss the role of the interstellar medium in cosmic-ray particle acceleration. △ Less

Submitted 18 June, 2021; originally announced June 2021.

Comments: 20 pages, 14 figures, 3 tables, accepted for publication in The Astrophysical Journal

Pursuing the origin of the gamma rays in RX J1713.7$-$3946 quantifying the hadronic and leptonic components

Authors: Yasuo Fukui, Hidetoshi Sano, Yumiko Yamane, Takahiro Hayakawa, Tsuyoshi Inoue, Kengo Tachihara, Gavin Rowell, Sabrina Einecke

Abstract: We analyzed the TeV gamma-ray image of a supernova remnant RX J1713.7$-$3946 (RX J1713) through a comparison with the interstellar medium (ISM) and the non-thermal X-rays. The gamma-ray datasets at two energy bands of $>$2 TeV and $>$250-300 GeV were obtained with H.E.S.S. (H.E.S.S. Collaboration 2018; Aharonian et al. 2007) and utilized in the analysis. We employed a new methodology which assumes… ▽ More We analyzed the TeV gamma-ray image of a supernova remnant RX J1713.7$-$3946 (RX J1713) through a comparison with the interstellar medium (ISM) and the non-thermal X-rays. The gamma-ray datasets at two energy bands of $>$2 TeV and $>$250-300 GeV were obtained with H.E.S.S. (H.E.S.S. Collaboration 2018; Aharonian et al. 2007) and utilized in the analysis. We employed a new methodology which assumes that the gamma-ray counts are expressed by a linear combination of two terms; one is proportional to the ISM column density and the other proportional to the X-ray count. We then assume these represent the hadronic and leptonic components, respectively. By fitting the expression to the data pixels, we find that the gamma-ray counts are well represented by a flat plane in a 3D space of the gamma-ray counts, the ISM column density and the X-ray counts. The results using the latest H.E.S.S. data at 4.8 arcmin resolution show that the hadronic and leptonic components occupy $(67\pm8)$% and $(33\pm8)$% of the total gamma rays, respectively, where the two components have been quantified for the first time. The hadronic component is greater than the leptonic component, which reflects the massive ISM of $\sim$10$^4$ $M_{\odot}$ associated with the SNR, lending support for the acceleration of the cosmic-ray protons. There is a marginal hint that the gamma rays are suppressed at high gamma-ray counts which may be ascribed to the second order effects including the shock-cloud interaction and the penetration-depth effect. △ Less

Submitted 22 June, 2021; v1 submitted 6 May, 2021; originally announced May 2021.

Comments: 20 pages, 9 figures, 3 tables, accepted for publication in The Astrophysical Journal (ApJ)

Giant Molecular Cloud Formation at the Interface of Colliding Supershells in the Large Magellanic Cloud

Authors: Kosuke Fujii, Norikazu Mizuno, J. R. Dawson, Tsuyoshi Inoue, Kazufumi Torii, Toshikazu Onishi, Akiko Kawamura, Erik Muller, Tetsuhiro Minamidani, Kisetsu Tsuge, Yasuo Fukui

Abstract: We investigate the Hi envelope of the young, massive GMCs in the star-forming regions N48 and N49, which are located within the high column density Hi ridge between two kpc-scale supergiant shells, LMC 4 and LMC 5. New long-baseline Hi 21 cm line observations with the Australia Telescope Compact Array (ATCA) were combined with archival shorter baseline data and single dish data from the Parkes tel… ▽ More We investigate the Hi envelope of the young, massive GMCs in the star-forming regions N48 and N49, which are located within the high column density Hi ridge between two kpc-scale supergiant shells, LMC 4 and LMC 5. New long-baseline Hi 21 cm line observations with the Australia Telescope Compact Array (ATCA) were combined with archival shorter baseline data and single dish data from the Parkes telescope, for a final synthesized beam size of 24.75" by 20.48", which corresponds to a spatial resolution of ~ 6 pc in the LMC. It is newly revealed that the Hi gas is highly filamentary, and that the molecular clumps are distributed along filamentary Hi features. In total 39 filamentary features are identified and their typical width is ~ 21 (8-49) [pc]. We propose a scenario in which the GMCs were formed via gravitational instabilities in atomic gas which was initially accumulated by the two shells and then further compressed by their collision. This suggests that GMC formation involves the filamentary nature of the atomic medium. △ Less

Submitted 29 April, 2021; originally announced April 2021.

Comments: 22 pages, 4 tables, 16 figures

The eROSITA Final Equatorial-Depth Survey (eFEDS): An X-ray bright, extremely luminous infrared galaxy at z = 1.87

Authors: Yoshiki Toba, Marcella Brusa, Teng Liu, Johannes Buchner, Yuichi Terashima, Tanya Urrutia, Mara Salvato, Masayuki Akiyama, Riccardo Arcodia, Andy D. Goulding, Yuichi Higuchi, Kaiki T. Inoue, Toshihiro Kawaguchi, Georg Lamer, Andrea Merloni, Tohru Nagao, Yoshihiro Ueda, Kirpal Nandra

Abstract: In this study, we investigate the X-ray properties of WISE J090924.01+000211.1 (WISEJ0909+0002), an extremely luminous infrared (IR) galaxy (ELIRG) at $z_{\rm spec}$= 1.871 in the eROSITA final equatorial depth survey (eFEDS). WISEJ0909+0002 is a WISE 22 $μ$m source, located in the GAMA-09 field, detected by eROSITA during the performance and verification phase. The corresponding optical spectrum… ▽ More In this study, we investigate the X-ray properties of WISE J090924.01+000211.1 (WISEJ0909+0002), an extremely luminous infrared (IR) galaxy (ELIRG) at $z_{\rm spec}$= 1.871 in the eROSITA final equatorial depth survey (eFEDS). WISEJ0909+0002 is a WISE 22 $μ$m source, located in the GAMA-09 field, detected by eROSITA during the performance and verification phase. The corresponding optical spectrum indicates that this object is a type-1 active galactic nucleus (AGN). Observations from eROSITA combined with Chandra and XMM-Newton archival data indicate a very luminous ($L$ (2--10 keV) = ($2.1 \pm 0.2) \times 10^{45}$ erg s$^{-1}$) unobscured AGN with a power-law photon index of $Γ$ = 1.73$_{-0.15}^{+0.16}$, and an absorption hydrogen column density of $\log\,(N_{\rm H}/{\rm cm}^{-2}) < 21.0$. The IR luminosity was estimated to be $L_{\rm IR}$ = (1.79 $\pm$ 0.09) $\times 10^{14}\, L_{\odot}$ from spectral energy distribution modeling based on 22 photometric data (X-ray to far-IR) with X-CIGALE, which confirmed that WISEJ0909+0002 is an ELIRG. A remarkably high $L_{\rm IR}$ despite very low $N_{\rm H}$ would indicate that we are witnessing a short-lived phase in which hydrogen gas along the line of sight is blown outwards, whereas warm and hot dust heated by AGNs still exist. As a consequence of eROSITA all-sky survey, $6.8_{-5.6}^{+16}\times 10^2$ such X-ray bright ELIRGs are expected to be discovered in the entire extragalactic sky ($|b| > 10^\circ$). This can potentially be the key population to constrain the bright-end of IR luminosity functions. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Comments: 10 pages, 5 figures, and 3 tables, accepted for publication in A&A Letters (special Issue: First science highlights from SRG/eROSITA)

Journal ref: A&A 649, L11 (2021)

Revealing the diverse magnetic field morphologies in Taurus dense cores with sensitive sub-millimeter polarimetry

Authors: Chakali Eswaraiah, Di Li, Ray S. Furuya, Tetsuo Hasegawa, Derek Ward-Thompson, Keping Qiu, Nagayoshi Ohashi, Kate Pattle, Sarah Sadavoy, Charles L. H. Hull, David Berry, Yasuo Doi, Tao-Chung Ching, Shih-Ping Lai, Jia-Wei Wang, Patrick M. Koch, Jungmi Kwon, Woojin Kwon, Pierre Bastien, Doris Arzoumanian, Simon Coudé, Archana Soam, Lapo Fanciullo, Hsi-Wei Yen, Junhao Liu , et al. (120 additional authors not shown)

Abstract: We have obtained sensitive dust continuum polarization observations at 850 $μ$m in the B213 region of Taurus using POL-2 on SCUBA-2 at the James Clerk Maxwell Telescope (JCMT), as part of the BISTRO (B-fields in STar-forming Region Observations) survey. These observations allow us to probe magnetic field (B-field) at high spatial resolution ($\sim$2000 au or $\sim$0.01 pc at 140 pc) in two protost… ▽ More We have obtained sensitive dust continuum polarization observations at 850 $μ$m in the B213 region of Taurus using POL-2 on SCUBA-2 at the James Clerk Maxwell Telescope (JCMT), as part of the BISTRO (B-fields in STar-forming Region Observations) survey. These observations allow us to probe magnetic field (B-field) at high spatial resolution ($\sim$2000 au or $\sim$0.01 pc at 140 pc) in two protostellar cores (K04166 and K04169) and one prestellar core (Miz-8b) that lie within the B213 filament. Using the Davis-Chandrasekhar-Fermi method, we estimate the B-field strengths in K04166, K04169, and Miz-8b to be 38$\pm$14 $μ$G, 44$\pm$16 $μ$G, and 12$\pm$5 $μ$G, respectively. These cores show distinct mean B-field orientations. B-field in K04166 is well ordered and aligned parallel to the orientations of the core minor axis, outflows, core rotation axis, and large-scale uniform B-field, in accordance with magnetically regulated star formation via ambipolar diffusion taking place in K04166. B-field in K04169 is found to be ordered but oriented nearly perpendicular to the core minor axis and large-scale B-field, and not well-correlated with other axes. In contrast, Miz-8b exhibits disordered B-field which show no preferred alignment with the core minor axis or large-scale field. We found that only one core, K04166, retains a memory of the large-scale uniform B-field. The other two cores, K04169 and Miz-8b, are decoupled from the large-scale field. Such a complex B-field configuration could be caused by gas inflow onto the filament, even in the presence of a substantial magnetic flux. △ Less

Submitted 3 March, 2021; originally announced March 2021.

Comments: 20 pages (10 are main), 5 figures (3 are main), and 3 tables (2 are main); Accepted for its publication in ApJL

Dust polarized emission observations of NGC 6334; BISTRO reveals the details of the complex but organized magnetic field structure of the high-mass star-forming hub-filament network

Authors: D. Arzoumanian, R. Furuya, T. Hasegawa, M. Tahani, S. Sadavoy, C. L. H. Hull, D. Johnstone, P. M. Koch, S. -i. Inutsuka, Y. Doi, T. Hoang, T. Onaka, K. Iwasaki, Y. Shimajiri, T. Inoue, N. Peretto, P. André, P. Bastien, D. Berry, H. -R. V. Chen, J. Di Francesco, C. Eswaraiah, L. Fanciullo, L. M. Fissel, J. Hwang , et al. (123 additional authors not shown)

Abstract: [Abridged] Filaments and hubs have received special attention recently thanks to studies showing their role in star formation. While the column density and velocity structures of both filaments and hubs have been studied, their magnetic fields (B-field) are not yet characterized. We aim to understand the role of the B-field in the dynamical evolution of the NGC 6334 hub-filament network. We presen… ▽ More [Abridged] Filaments and hubs have received special attention recently thanks to studies showing their role in star formation. While the column density and velocity structures of both filaments and hubs have been studied, their magnetic fields (B-field) are not yet characterized. We aim to understand the role of the B-field in the dynamical evolution of the NGC 6334 hub-filament network. We present new observations of the dust polarized emission at 850$μ$m towards NGC 6334 obtained with the JCMT/POL-2. We study the distribution and dispersion of the polarized intensity ($PI$), the polarization fraction ($PF$), and the B-field angle ($θ_{B}$). We derive the power spectrum of the intensity and $θ_{B}$ along the ridge crest. Our analyses show a complex B-field structure when observed over the whole region ($\sim10$ pc), however, at smaller scales ($\sim1$ pc), $θ_{B}$ varies coherently along the filaments. The observed power spectrum of $θ_{B}$ can be well represented with a power law function with a slope $-1.33\pm0.23$, which is $\sim20\%$ shallower than that of $I$. This result is compatible with the properties of simulated filaments and may indicate the processes at play in the formation of filaments. $θ_{B}$ rotates from being mostly perpendicular to the filament crests to mostly parallel as they merge with the hubs. This variation of $θ_{B}$ may be tracing local velocity flows of matter in-falling onto the hubs. Our analysis suggests a variation of the energy balance along the crests of these filaments, from magnetically critical/supercritical at their far ends to magnetically subcritical near the hubs. We detect an increase of $PF$ towards the high-column density star cluster-forming hubs that may result from the increase of grain alignment efficiency due to stellar radiation from the newborn stars. △ Less

Submitted 23 December, 2020; originally announced December 2020.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 647, A78 (2021)

Observations of magnetic fields surrounding LkH$α$ 101 taken by the BISTRO survey with JCMT-POL-2

Authors: Nguyen Bich Ngoc, Pham Ngoc Diep, Harriet Parsons, Kate Pattle, Thiem Hoang, Derek Ward-Thompson, Le Ngoc Tram, Charles L. H. Hull, Mehrnoosh Tahani, Ray Furuya, Pierre Bastien, Keping Qiu, Tetsuo Hasegawa, Woojin Kwon, Yasuo Doi, Shih-Ping Lai, Simon Coude, David Berry, Tao-Chung Ching, Jihye Hwang, Archana Soam, Jia-Wei Wang, Doris Arzoumanian, Tyler L. Bourke, Do-Young Byun , et al. (124 additional authors not shown)

Abstract: We report the first high spatial resolution measurement of magnetic fields surrounding LkH$α$ 101, a part of the Auriga-California molecular cloud. The observations were taken with the POL-2 polarimeter on the James Clerk Maxwell Telescope within the framework of the B-fields In Star-forming Region Observations (BISTRO) survey. Observed polarization of thermal dust emission at 850 $μ$m is found to… ▽ More We report the first high spatial resolution measurement of magnetic fields surrounding LkH$α$ 101, a part of the Auriga-California molecular cloud. The observations were taken with the POL-2 polarimeter on the James Clerk Maxwell Telescope within the framework of the B-fields In Star-forming Region Observations (BISTRO) survey. Observed polarization of thermal dust emission at 850 $μ$m is found to be mostly associated with the red-shifted gas component of the cloud. The magnetic field displays a relatively complex morphology. Two variants of the Davis-Chandrasekhar-Fermi method, unsharp masking and structure function, are used to calculate the strength of magnetic fields in the plane of the sky, yielding a similar result of $B_{\rm POS}\sim 115$ $\mathrmμ$G. The mass-to-magnetic-flux ratio in critical value units, $λ\sim0.3$, is the smallest among the values obtained for other regions surveyed by POL-2. This implies that the LkH$α$ 101 region is sub-critical and the magnetic field is strong enough to prevent gravitational collapse. The inferred $δB/B_0\sim 0.3$ implies that the large scale component of the magnetic field dominates the turbulent one. The variation of the polarization fraction with total emission intensity can be fitted by a power-law with an index of $α=0.82\pm0.03$, which lies in the range previously reported for molecular clouds. We find that the polarization fraction decreases rapidly with proximity to the only early B star (LkH$α$ 101) in the region. The magnetic field tangling and the joint effect of grain alignment and rotational disruption by radiative torques are potential of explaining such a decreasing trend. △ Less

Submitted 8 December, 2020; originally announced December 2020.

Comments: 25 pages, 19 figures, Accepted for publication in The Astrophysical Journal

Classification of Filament Formation Mechanisms in Magnetized Molecular Clouds

Authors: Daisei Abe, Tsuyoshi Inoue, Shu-ichiro Inutsuka, Tomoaki Matsumoto

Abstract: Recent observations of molecular clouds show that dense filaments are the sites of present-day star formation. Thus, it is necessary to understand the filament formation process because these filaments provide the initial condition for star formation. Theoretical research suggests that shock waves in molecular clouds trigger filament formation. Since several different mechanisms have been proposed… ▽ More Recent observations of molecular clouds show that dense filaments are the sites of present-day star formation. Thus, it is necessary to understand the filament formation process because these filaments provide the initial condition for star formation. Theoretical research suggests that shock waves in molecular clouds trigger filament formation. Since several different mechanisms have been proposed for filament formation, the formation mechanism of the observed star-forming filaments requires clarification. In the present study, we perform a series of isothermal magnetohydrodynamics simulations of filament formation. We focus on the influences of shock velocity and turbulence on the formation mechanism and identified three different mechanisms for the filament formation. The results indicate that when the shock is fast, at shock velocity v_sh = 7 km/s, the gas flows driven by the curved shock wave create filaments irrespective of the presence of turbulence and self-gravity. However, at a slow shock velocity v_sh = 2.5 km/s, the compressive flow component involved in the initial turbulence induces filament formation. When both the shock velocities and turbulence are low, the self-gravity in the shock-compressed sheet becomes important for filament formation. Moreover, we analyzed the line-mass distribution of the filaments and showed that strong shock waves can naturally create high-line-mass filaments such as those observed in the massive star-forming regions in a short time. We conclude that the dominant filament formation mode changes with the velocity of the shock wave triggering the filament formation. △ Less

Submitted 13 August, 2021; v1 submitted 3 December, 2020; originally announced December 2020.

Comments: 25 pages, 21 figures, accepted to ApJ

Journal ref: The Astrophysical Journal, Volume 916, Issue 2, id.83, 22 pp., (2021)