Showing 1–50 of 154 results for author: Huang, K

Supermassive black holes and very high-energy neutrinos: the case of NGC 1068

Authors: P. Padovani, E. Resconi, M. Ajello, C. Bellenghi, S. Bianchi, P. Blasi, K. -Y. Huang, S. Gabici, V. Gámez Rosas, H. Niederhausen, E. Peretti, B. Eichmann, D. Guetta, A. Lamastra, T. Shimizu

Abstract: We present a comprehensive multi-messenger study of NGC 1068, the prototype Seyfert II galaxy recently associated with high-energy IceCube neutrinos. Various aspects of the source, including its nuclear activity, jet, outflow, and starburst region, are analyzed in detail using a multi-wavelength approach and relevant luminosities are derived. We then explore its gamma-ray and neutrino emissions an… ▽ More We present a comprehensive multi-messenger study of NGC 1068, the prototype Seyfert II galaxy recently associated with high-energy IceCube neutrinos. Various aspects of the source, including its nuclear activity, jet, outflow, and starburst region, are analyzed in detail using a multi-wavelength approach and relevant luminosities are derived. We then explore its gamma-ray and neutrino emissions and investigate potential mechanisms underlying these phenomena and their relations with the different astrophysical components to try to understand which one is responsible for the IceCube neutrinos. By first using simple order-of-magnitude arguments and then applying specific theoretical models, we infer that only the region close to the accretion disc around the supermassive black hole has both the right density of X-ray photons needed to provide the targets for protons to sustain neutrino production and of optical/infrared photons required to absorb the associated but unobserved gamma rays. We conclude by highlighting ongoing efforts to constrain a possible broad connection between neutrinos and active galactic nuclei, as well as future synergies between astronomical and neutrino facilities. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: 18 pages, 4 figures, review paper accepted for publication in Nature Astronomy; authors' version before editorial review

An ALCHEMI inspection of sulphur-bearing species towards the central molecular zone of NGC 253

Authors: M. Bouvier, S. Viti, E. Behrens, J. Butterworth, K. -Y. Huang, J. G. Mangum, N. Harada, S. Martín, V. M. Rivilla, S. Muller, K. Sakamoto, Y. Yoshimura, K. Tanaka, K. Nakanishi, R. Herrero-Illana, L. Colzi, M. D. Gorski, C. Henkel, P. K. Humire, D. S. Meier, P. P. van der Werf, Y. T. Yan

Abstract: Sulphur-bearing species are detected in various environments within Galactic star-forming regions and are particularly abundant in the gas phase of outflows and shocks, and photo-dissociation regions. In this work, we aim to investigate the nature of the emission from the most common sulphur-bearing species observable at millimetre wavelengths towards the nuclear starburst of the galaxy NGC 253. W… ▽ More Sulphur-bearing species are detected in various environments within Galactic star-forming regions and are particularly abundant in the gas phase of outflows and shocks, and photo-dissociation regions. In this work, we aim to investigate the nature of the emission from the most common sulphur-bearing species observable at millimetre wavelengths towards the nuclear starburst of the galaxy NGC 253. We intend to understand which type of regions are probed by sulphur-bearing species and which process(es) dominate(s) the release of sulphur into the gas phase. We used the high-angular resolution (1.6" or 27 pc) observations from the ALCHEMI ALMA Large Program to image several sulphur-bearing species towards the central molecular zone (CMZ) of NGC 253. We performed local thermodynamic equilibrium (LTE) and non-LTE large velocity gradient (LVG) analyses to derive the physical conditions of the gas in which S-bearing species are emitted, and their abundance ratios across the CMZ. Finally, we compared our results with previous ALCHEMI studies and a few selected Galactic environments. We found that not all sulphur-bearing species trace the same type of gas: strong evidence indicates that H2S and part of the emission of OCS, H2CS, and SO, are tracing shocks whilst part of SO and CS emission rather trace the dense molecular gas. For some species, such as CCS and SO2, we could not firmly conclude on their origin of emission. The present analysis indicates that the emission from most sulphur-bearing species throughout the CMZ is likely dominated by shocks associated with ongoing star formation. In the inner part of the CMZ where the presence of super star clusters was previously indicated, we could not distinguish between shocks or thermal evaporation as the main process releasing the S-bearing species. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 44 pages, 20 figures, Accepted for publication in A&A

Dynamical Model of Rotation and Orbital Coupling for Deimos

Authors: Kai Huang, Lijun Zhang, Yongzhang Yang, Mao Ye, Yuqiang Li

Abstract: This paper introduces a novel dynamical model, building upon the existing dynamical model for Deimos in the current numerical ephemerides, which only encompasses the simple libration effects of Deimos. The study comprehensively incorporates the rotational dynamics of Deimos influenced by the torque exerted by the major celestial bodies (Mars, the Sun) in the solar system within the inertial space.… ▽ More This paper introduces a novel dynamical model, building upon the existing dynamical model for Deimos in the current numerical ephemerides, which only encompasses the simple libration effects of Deimos. The study comprehensively incorporates the rotational dynamics of Deimos influenced by the torque exerted by the major celestial bodies (Mars, the Sun) in the solar system within the inertial space. Consequently, a full dynamical model is formulated to account for the complete coupling between the rotation and orbit of Deimos. Simultaneously, employing precision orbit determination methods used for artificial satellites, we develop an adjustment model for fitting data to the complete model. The 12-order Adams--Bashforth--Moulton (ABM) integration algorithm is employed to synchronously integrate the 12 state variables of the full model to obtain the orbit of Deimos. Numerical simulation results indicate that the full dynamical model and adjustment model are stable and reliable. Compared to the simple model, the polar axis of Deimos in the inertial space exhibits a more complex oscillation in the full model. This work further advances the current dynamical model for Deimos and establishes the foundational model for the generation of a new set of precise numerical ephemerides for Deimos. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Multi-transition study of methanol towards NGC 1068 with ALMA

Authors: K. -Y. Huang, D. Abbink, S. Viti, S. García-Burillo

Abstract: The outflowing molecular gas in the circumnuclear disk (CND) of the nearby (D=14 Mpc) AGN-starburst composite galaxy NGC 1068 is considered as a manifestation of ongoing AGN feedback. The large spread of velocities from the outflowing gas is likely driving various kinds of shock chemistry across the CND. We performed a multiline molecular study using CH3OH with the aim of characterizing the gas pr… ▽ More The outflowing molecular gas in the circumnuclear disk (CND) of the nearby (D=14 Mpc) AGN-starburst composite galaxy NGC 1068 is considered as a manifestation of ongoing AGN feedback. The large spread of velocities from the outflowing gas is likely driving various kinds of shock chemistry across the CND. We performed a multiline molecular study using CH3OH with the aim of characterizing the gas properties probed by CH3OH in the CND of NGC 1068, and investigating its potential association with molecular shocks. Multi-transition CH3OH were imaged at the resolution of 0.''5-0.''8 with the Atacama Large Millimeter/submillimeter Array (ALMA). We performed non-LTE radiative transfer analysis coupled with a Bayesian inference process in order to determine the gas properties such as the gas volume density and the gas kinetic temperature. The gas densities traced by CH3OH point to $\sim 10^{6}$ cm\textsuperscript{-3} across all the CND regions. The gas kinetic temperature cannot be well constrained in any of the CND regions though the inferred temperature is likely low ($\lesssim 100$ K).The low gas temperature traced by CH3OH suggests shocks and subsequent fast cooling as the origin of the observed gas-phase CH3OH abundance. We also note that the E-/A- isomer column density ratio inferred is fairly close to unity, which is interestingly different from the Galactic measurements in the literature. It remains inconclusive whether CH3OH exclusively traces slow and non-dissociative shocks, or whether the CH3OH abundance can actually be boosted in both fast and slow shocks. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Physical Properties of the Southwest Outflow Streamer in the Starburst Galaxy NGC 253 with ALCHEMI

Authors: Min Bao, Nanase Harada, Kotaro Kohno, Yuki Yoshimura, Fumi Egusa, Yuri Nishimura, Kunihiko Tanaka, Kouichiro Nakanishi, Sergio Martín, Jeffrey G. Mangum, Kazushi Sakamoto, Sébastien Muller, Mathilde Bouvier, Laura Colzi, Kimberly L. Emig, David S. Meier, Christian Henkel, Pedro Humire, Ko-Yun Huang, Víctor M. Rivilla, Paul van der Werf, Serena Viti

Abstract: The physical properties of galactic molecular outflows are important as they could constrain outflow formation mechanisms. We study the properties of the southwest (SW) outflow streamer including gas kinematics, optical depth, dense gas fraction, and shock strength in the central molecular zone of the starburst galaxy NGC 253. We image the molecular emission at a spatial resolution of $\sim$27 pc… ▽ More The physical properties of galactic molecular outflows are important as they could constrain outflow formation mechanisms. We study the properties of the southwest (SW) outflow streamer including gas kinematics, optical depth, dense gas fraction, and shock strength in the central molecular zone of the starburst galaxy NGC 253. We image the molecular emission at a spatial resolution of $\sim$27 pc based on data from the ALCHEMI program. We trace the kinematics of molecular gas with CO(1-0) line. We constrain the optical depth of CO emission with CO/$^{13}$CO(1-0) ratio, the dense gas fraction with HCN/CO(1-0) ratio, as well as the shock strength with SiO(2-1)/$^{13}$CO(1-0) ratio. The CO/$^{13}$CO(1-0) integrated intensity ratio is $\sim$21 in the SW streamer region, which approximates the C/$^{13}$C isotopic abundance ratio. The higher integrated intensity ratio compared to the disk can be attributed to the optically thinner environment for CO(1-0) emission inside the SW streamer. The HCN/CO(1-0) and SiO(2-1)/$^{13}$CO(1-0) integrated intensity ratios both approach $\sim$0.2 in three giant molecular clouds (GMCs) at the base of the outflow streamers, which implies the higher dense gas fraction and enhanced strength of fast shocks in those GMCs than in the disk. The contours of those two integrated intensity ratios are extended towards the directions of outflow streamers, which connects the enhanced dense gas fraction and shock strength with molecular outflow. Moreover, the molecular gas with enhanced dense gas fraction and shock strength located at the base of the SW streamer shares the same velocity with the outflow. These phenomena suggest that the star formation inside the GMCs can trigger the shocks and further drive the molecular outflow. △ Less

Submitted 6 April, 2024; originally announced April 2024.

Comments: Accepted for publication in A&A, 14 pages, 11 figures

Density-dependent quark mean-field model for nuclear matter and neutron stars

Authors: Kaixuan Huang, Jinniu Hu, Ying Zhang, Hong Shen

Abstract: We develop a density-dependent quark mean-field (DDQMF) model to study the properties of nuclear matter and neutron stars, where the coupling strength between $σ$ meson and nucleon is generated by the degree of freedom of quarks, while other meson coupling constants are regarded as density-dependent ones. Two values for the nucleon effective mass, $M^*_{N0}/M_N=0.556,~0.70$ at the saturation densi… ▽ More We develop a density-dependent quark mean-field (DDQMF) model to study the properties of nuclear matter and neutron stars, where the coupling strength between $σ$ meson and nucleon is generated by the degree of freedom of quarks, while other meson coupling constants are regarded as density-dependent ones. Two values for the nucleon effective mass, $M^*_{N0}/M_N=0.556,~0.70$ at the saturation density are chosen based on the consideration of the core-collapse supernova simulation and finite nuclei when the meson-nucleon coupling constants are fixed. We find that the equation of state (EOS) of nuclear matter, the symmetry energy, the mass-radius relations, and the tidal deformabilities of neutron stars with larger nucleon effective mass are more sensitive to the skewness coefficient $J_0$. The EOSs with $M^*_{N0}/M_N=0.70$ are softer when the skewness coefficient $J_0=-800$ MeV. However, the maximum masses of the neutron star can be around $2.32M_\odot$ with $J_0=400$ MeV regardless of the value of the nucleon effective mass. By manipulating the coupling strength of the isovector meson to generate different slopes of symmetry energy, we construct the neutron star EOSs that can satisfy the different variables from the simultaneous mass-radius measurements of PSR J0030+0451, PSR J0740+6620 by the NICER collaboration, the mass-radius relations of HESS J1731-347, and the radius constraints from the gravitational-wave signal GW170817 in the framework of DDQMF model. At the same time, most of these constructed EOSs can also satisfy the constraints of the tidal deformability from GW170817 event. △ Less

Submitted 3 March, 2024; originally announced March 2024.

Molecular isotopologue measurements toward super star clusters and the relation to their ages in NGC253 with ALCHEMI

Authors: J. Butterworth, S. Viti, P. P. Van der Werf, J. G. Mangum, S. Martín, N. Harada, K. L. Emig, S. Muller, K. Sakamoto, Y. Yoshimura, K. Tanaka, R. Herrero-Illana, L. Colzi, V. M. Rivilla, K. Y. Huang, M. Bouvier, E. Behrens, C. Henkel, Y. T. Yan, D. S. Meier, D. Zhou

Abstract: Determining the evolution of the CNO isotopes in the interstellar medium (ISM) of starburst galaxies can yield important constraints on the ages of superstar clusters (SSCs), or on other aspects and contributing factors of their evolution. Due to the time-dependent nature of the abundances of isotopes within the ISM as they are supplied from processes such as nucleosynthesis or chemical fractionat… ▽ More Determining the evolution of the CNO isotopes in the interstellar medium (ISM) of starburst galaxies can yield important constraints on the ages of superstar clusters (SSCs), or on other aspects and contributing factors of their evolution. Due to the time-dependent nature of the abundances of isotopes within the ISM as they are supplied from processes such as nucleosynthesis or chemical fractionation, this provides the possible opportunity to probe the ability of isotopes ratios to trace the ages of high star forming regions, such as SSCs. The goal of this study is to investigate whether the isotopic variations in SSC regions within NGC253 are correlated with their different ages as derived from stellar population modelling. We have measured abundance ratios of CO, HCN and HCO$^+$ isotopologues in six regions containing SSCs within NGC253 using high spatial resolution (1.6",$\sim 28$pc) data from the ALCHEMI (ALma Comprehensive High-resolution Extragalactic Molecular Inventory) ALMA Large program. We have then analysed these ratios using RADEX radiative transfer modelling, with the parameter space sampled using the nested sampling Monte Carlo algorithm MLFriends. These abundance ratios were then compared to ages predicted in each region via the fitting of observed star formation tracers (such as Br$γ$) to starburst stellar population evolution models. We do not find any significant trend with age for the CO and HCN isotopologue ratios on the timescales for the ages of the SSC* regions observed. The driving factors of these ratios within SSCs could be the Initial Mass Function as well as possibly fractionation effects. To further probe these effects in SSCs over time a larger sample of SSCs must be observed spanning a larger age range. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Comments: 44 pages, 43 Figures, Accepted for Publication to A&A

Journal ref: A&A 686, A31 (2024)

Nuclear mass table in deformed relativistic Hartree-Bogoliubov theory in continuum, II: Even-$Z$ nuclei

Authors: DRHBc Mass Table Collaboration, Peng Guo, Xiaojie Cao, Kangmin Chen, Zhihui Chen, Myung-Ki Cheoun, Yong-Beom Choi, Pak Chung Lam, Wenmin Deng, Jianmin Dong, Pengxiang Du, Xiaokai Du, Kangda Duan, Xiaohua Fan, Wei Gao, Lisheng Geng, Eunja Ha, Xiao-Tao He, Jinniu Hu, Jingke Huang, Kun Huang, Yanan Huang, Zidan Huang, Kim Da Hyung, Hoi Yat Chan , et al. (58 additional authors not shown)

Abstract: The mass table in the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) with the PC-PK1 density functional has been established for even-$Z$ nuclei with $8\le Z\le120$, extended from the previous work for even-even nuclei [Zhang $\it{et.~al.}$ (DRHBc Mass Table Collaboration), At. Data Nucl. Data Tables 144, 101488 (2022)]. The calculated binding energies, two-nucleon and one-ne… ▽ More The mass table in the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) with the PC-PK1 density functional has been established for even-$Z$ nuclei with $8\le Z\le120$, extended from the previous work for even-even nuclei [Zhang $\it{et.~al.}$ (DRHBc Mass Table Collaboration), At. Data Nucl. Data Tables 144, 101488 (2022)]. The calculated binding energies, two-nucleon and one-neutron separation energies, root-mean-square (rms) radii of neutron, proton, matter, and charge distributions, quadrupole deformations, and neutron and proton Fermi surfaces are tabulated and compared with available experimental data. A total of 4829 even-$Z$ nuclei are predicted to be bound, with an rms deviation of 1.477 MeV from the 1244 mass data. Good agreement with the available experimental odd-even mass differences, $α$ decay energies, and charge radii is also achieved. The description accuracy for nuclear masses and nucleon separation energies as well as the prediction for drip lines is compared with the results obtained from other relativistic and nonrelativistic density functional. The comparison shows that the DRHBc theory with PC-PK1 provides an excellent microscopic description for the masses of even-$Z$ nuclei. The systematics of the nucleon separation energies, odd-even mass differences, pairing energies, two-nucleon gaps, $α$ decay energies, rms radii, quadrupole deformations, potential energy curves, neutron density distributions, and neutron mean-field potentials are discussed. △ Less

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

Comments: 394 pages, 17 figures, 2 tables, published in Atomic Data and Nuclear Data Tables, data file in the TXT form is available for download under "Ancillary files"

Journal ref: Peng Guo, et. al. (DRHBc Mass Table Collaboration), Atomic Data and Nuclear Data Tables 158 (2024) 101661

The ALCHEMI atlas: principal component analysis reveals starburst evolution in NGC 253

Authors: Nanase Harada, David S. Meier, Sergio Martín, Sebastien Muller, Kazushi Sakamoto, Toshiki Saito, Mark D. Gorski, Christian Henkel, Kunihiko Tanaka, Jeffrey G. Mangum, Susanne Aalto, Rebeca Aladro, Mathilde Bouvier, Laura Colzi, Kimberly L. Emig, Rubén Herrero-Illana, Ko-Yun Huang, Kotaro Kohno, Sabine König, Kouichiro Nakanishi, Yuri Nishimura, Shuro Takano, Víctor M. Rivilla, Serena Viti, Yoshimasa Watanabe , et al. (2 additional authors not shown)

Abstract: Molecular lines are powerful diagnostics of the physical and chemical properties of the interstellar medium (ISM). These ISM properties, which affect future star formation, are expected to differ in starburst galaxies from those of more quiescent galaxies. We investigate the ISM properties in the central molecular zone of the nearby starburst galaxy NGC 253 using the ultra-wide millimeter spectral… ▽ More Molecular lines are powerful diagnostics of the physical and chemical properties of the interstellar medium (ISM). These ISM properties, which affect future star formation, are expected to differ in starburst galaxies from those of more quiescent galaxies. We investigate the ISM properties in the central molecular zone of the nearby starburst galaxy NGC 253 using the ultra-wide millimeter spectral scan survey from the ALMA Large Program ALCHEMI. We present an atlas of velocity-integrated images at a 1".6 resolution of 148 unblended transitions from 44 species, including the first extragalactic detection of HCNH$^+$ and the first interferometric images of C$_3$H$^+$, NO, HCS$^+$. We conduct a principal component analysis (PCA) on these images to extract correlated chemical species and to identify key groups of diagnostic transitions. To the best of our knowledge, our dataset is currently the largest astronomical set of molecular lines to which PCA has been applied. The PCA can categorize transitions coming from different physical components in NGC 253 such as i) young starburst tracers characterized by high-excitation transitions of HC$_3$N and complex organic molecules (COMs) versus tracers of on-going star formation (radio recombination lines) and high-excitation transitions of CCH and CN tracing PDRs, ii) tracers of cloud-collision-induced shocks (low-excitation transitions of CH$_3$OH, HNCO, HOCO$^+$, and OCS) versus shocks from star-formation-induced outflows (high-excitation transitions of SiO), as well as iii) outflows showing emission from HOC$^+$, CCH, H$_3$O$^+$, CO isotopologues, HCN, HCO$^+$, CS, and CN. Our findings show these intensities vary with galactic dynamics, star formation activities, and stellar feedback. △ Less

Submitted 4 January, 2024; originally announced January 2024.

Comments: 65 pages, 39 figures. Accepted for publication in ApJS

Satellite Impact on Astronomical Observations Based on Elliptical Orbit Model

Authors: Tianzhu Hu, Yong Zhang, Xiangqun Cui, Zihuang Cao, Kang Huang, Jingyi Cai, Jun Li, Tong Zhou

Abstract: Space-based and ground-based telescopes have extensively documented the impact of satellites on astronomical observations. With the proliferation of satellite mega-constellation programs, their influence on astronomical observations has become undeniable. It is crucial to quantify the impact of satellites on telescopes. To address this need, we have enhanced the circular orbit model for satellites… ▽ More Space-based and ground-based telescopes have extensively documented the impact of satellites on astronomical observations. With the proliferation of satellite mega-constellation programs, their influence on astronomical observations has become undeniable. It is crucial to quantify the impact of satellites on telescopes. To address this need, we have enhanced the circular orbit model for satellites and introduced a methodology based on two-line element (TLE) orbit data. This involves constructing a satellite probability distribution model to evaluate the impact of satellites on telescopes. Using our method, we assessed the satellite impact on global observatories. The results indicate that the regions most severely affected by satellite interference currently are those near the equator, with latitudes around 50 and 80 degrees experiencing the most significant impact from low Earth orbit satellites. Furthermore, we validated the reliability of our method using imaging data obtained from the focal surface acquisition camera of the LAMOST telescope. △ Less

Submitted 21 December, 2023; originally announced December 2023.

Comments: 9 pages, 7 figures, submitted to A&A, any comments are welcome

Journal ref: A&A 687, A135 (2024)

Volume density structure of the NGC253 CMZ through ALCHEMI excitation analysis

Authors: Kunihiko Tanaka, Jeffrey G. Mangum, Serena Viti, Sergio Martin, Nanase Harada, Kazushi Sakamoto, Sebastien Muller, Yuki Yoshimura, Kouichiro Nakanishi, Ruben Herrero Illana, Kimberly L. Emig, S. Muhle, Hiroyuki Kaneko, Tomoka Tosaki, Erica Behrens, Victor M. Rivilla, Laura Colzi, Yuri Nishimura, P. K. Humire, Mathilde Bouvier, Ko-Yun Huang, Joshua Butterworth, David S. Meier, Paul P. van der Werf

Abstract: We present a spatially-resolved excitation analysis for the central molecular zone (CMZ) of the starburst galaxy NGC 253 using the data from the ALMA Large program ALCHEMI, whereby we explore parameters distinguishing NGC 253 from the quiescent Milky Way's Galactic Center (GC). Non-LTE analyses employing a hierarchical Bayesian framework are applied to Band 3-7 transitions from nine molecular spec… ▽ More We present a spatially-resolved excitation analysis for the central molecular zone (CMZ) of the starburst galaxy NGC 253 using the data from the ALMA Large program ALCHEMI, whereby we explore parameters distinguishing NGC 253 from the quiescent Milky Way's Galactic Center (GC). Non-LTE analyses employing a hierarchical Bayesian framework are applied to Band 3-7 transitions from nine molecular species to delineate the position-position-velocity distributions of column density ($N_\mathrm{H_2}$), volume density ($n_\mathrm{H_2}$), and temperature ($T_\mathrm{kin}$) at 27 pc resolution. Two distinct components are detected: a low-density component with $(n_\mathrm{H_2},\ T_\mathrm{kin})\sim(10^{3.3}\ \mathrm{cm}^{-3}, 85 K)$ and a high-density component with $(n_\mathrm{H_2},\ T_\mathrm{kin})\sim (10^{4.4}\ \mathrm{cm}^{-3}, 110\ \mathrm{K})$, separated at $n_\mathrm{H_2}\sim10^{3.8}\ \mathrm{cm}^{-3}$. NGC 253 has $\sim10$ times the high-density gas mass and $\sim3$ times the dense-gas mass fraction of the GC. These properties are consistent with their HCN/CO ratio but cannot alone explain the factor of $\sim30$ difference in their star formation efficiencies (SFEs), contradicting the dense-gas mass to star formation rate scaling law. The $n_\mathrm{H_2}$ histogram toward NGC 253 exhibits a shallow declining slope up to $n_\mathrm{H_2}\sim10^6\ \mathrm{cm}^{-3}$, while that of the GC steeply drops in $n_\mathrm{H_2}\gtrsim10^{4.5}\ \mathrm{cm}^{-3}$ and vanishes at $10^5\ \mathrm{cm}^{-3}$. Their dense-gas mass fraction ratio becomes consistent with their SFEs when the threshold $n_\mathrm{H_2}$ for the dense gas is taken at $\sim 10^{4.2\mbox{-}4.6}\ \mathrm{cm}^{-3}$. The rich abundance of gas above this density range in the NGC 253 CMZ, or its scarcity in the GC, is likely to be the critical difference characterizing the contrasting star formation in the centers of the two galaxies. △ Less

Submitted 20 November, 2023; originally announced November 2023.

Comments: 49 pages, 27 figures, 7 tables, accepted for publication in the Astrophysical Journal

LenSiam: Self-Supervised Learning on Strong Gravitational Lens Images

Authors: Po-Wen Chang, Kuan-Wei Huang, Joshua Fagin, James Hung-Hsu Chan, Joshua Yao-Yu Lin

Abstract: Self-supervised learning has been known for learning good representations from data without the need for annotated labels. We explore the simple siamese (SimSiam) architecture for representation learning on strong gravitational lens images. Commonly used image augmentations tend to change lens properties; for example, zoom-in would affect the Einstein radius. To create image pairs representing the… ▽ More Self-supervised learning has been known for learning good representations from data without the need for annotated labels. We explore the simple siamese (SimSiam) architecture for representation learning on strong gravitational lens images. Commonly used image augmentations tend to change lens properties; for example, zoom-in would affect the Einstein radius. To create image pairs representing the same underlying lens model, we introduce a lens augmentation method to preserve lens properties by fixing the lens model while varying the source galaxies. Our research demonstrates this lens augmentation works well with SimSiam for learning the lens image representation without labels, so we name it LenSiam. We also show that a pre-trained LenSiam model can benefit downstream tasks. We open-source our code and datasets at https://github.com/kuanweih/LenSiam . △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: 5 pages, 2 figures. Accepted by NeurIPS 2023 AI for Science Workshop

Real-time Monitoring for the Next Core-Collapse Supernova in JUNO

Authors: Angel Abusleme, Thomas Adam, Shakeel Ahmad, Rizwan Ahmed, Sebastiano Aiello, Muhammad Akram, Abid Aleem, Fengpeng An, Qi An, Giuseppe Andronico, Nikolay Anfimov, Vito Antonelli, Tatiana Antoshkina, Burin Asavapibhop, João Pedro Athayde Marcondes de André, Didier Auguste, Weidong Bai, Nikita Balashov, Wander Baldini, Andrea Barresi, Davide Basilico, Eric Baussan, Marco Bellato, Marco Beretta, Antonio Bergnoli , et al. (606 additional authors not shown)

Abstract: The core-collapse supernova (CCSN) is considered one of the most energetic astrophysical events in the universe. The early and prompt detection of neutrinos before (pre-SN) and during the supernova (SN) burst presents a unique opportunity for multi-messenger observations of CCSN events. In this study, we describe the monitoring concept and present the sensitivity of the system to pre-SN and SN neu… ▽ More The core-collapse supernova (CCSN) is considered one of the most energetic astrophysical events in the universe. The early and prompt detection of neutrinos before (pre-SN) and during the supernova (SN) burst presents a unique opportunity for multi-messenger observations of CCSN events. In this study, we describe the monitoring concept and present the sensitivity of the system to pre-SN and SN neutrinos at the Jiangmen Underground Neutrino Observatory (JUNO), a 20 kton liquid scintillator detector currently under construction in South China. The real-time monitoring system is designed to ensure both prompt alert speed and comprehensive coverage of progenitor stars. It incorporates prompt monitors on the electronic board as well as online monitors at the data acquisition stage. Assuming a false alert rate of 1 per year, this monitoring system exhibits sensitivity to pre-SN neutrinos up to a distance of approximately 1.6 (0.9) kiloparsecs and SN neutrinos up to about 370 (360) kiloparsecs for a progenitor mass of 30 solar masses, considering both normal and inverted mass ordering scenarios. The pointing ability of the CCSN is evaluated by analyzing the accumulated event anisotropy of inverse beta decay interactions from pre-SN or SN neutrinos. This, along with the early alert, can play a crucial role in facilitating follow-up multi-messenger observations of the next galactic or nearby extragalactic CCSN. △ Less

Submitted 4 December, 2023; v1 submitted 13 September, 2023; originally announced September 2023.

Comments: 24 pages, 9 figures, accepted for the publication at JCAP

Tracing the chemical footprint of shocks in AGN-host and starburst galaxies with ALMA multi-line molecular studies

Authors: Ko-Yun Huang, Serena Viti

Abstract: Multi-line molecular observations are an ideal tool for a systematic study of the physico-chemical processes in the Interstellar Medium (ISM), given the wide range of critical densities associated with different molecules and their transitions, and the dependencies of chemical reactions on the energy budget of the system. Recently high spatial resolution of typical shock tracers - SiO, HNCO, and C… ▽ More Multi-line molecular observations are an ideal tool for a systematic study of the physico-chemical processes in the Interstellar Medium (ISM), given the wide range of critical densities associated with different molecules and their transitions, and the dependencies of chemical reactions on the energy budget of the system. Recently high spatial resolution of typical shock tracers - SiO, HNCO, and CH3OH - have been studied in the potentially shocked regions in two nearby galaxies: NGC 1068 (an AGN-host galaxy) (Huang et al., Astron. Astrophys., 2022, 666, A102; Huang et al., in prep.) and NGC 253 (a starburst galaxy) (K.-Y. Huang et al., arXiv, 2023, preprint, arXiv:2303.12685, DOI: 10.48550/arXiv.2303.12685). This paper is dedicated to the comparative study of these two distinctively different galaxies, with the aim of determining the differences in their energetics and understanding large-scale shocks in dfferent types of galaxies. △ Less

Submitted 6 July, 2023; originally announced July 2023.

Comments: Accepted for publication in Faraday Discussions 2023. 18 pages, 7 figures and 2 tables

Intelligence of Astronomical Optical Telescope: Present Status and Future Perspectives

Authors: Kang Huang, Tianzhu Hu, Jingyi Cai, Xiushan Pang, Yonghui Hou, Yong Zhang, Huaiqing Wang, Xiangqun Cui

Abstract: Artificial intelligence technology has been widely used in astronomy, and new artificial intelligence technologies and application scenarios are constantly emerging. There have been a large number of papers reviewing the application of artificial intelligence technology in astronomy. However, relevant articles seldom mention telescope intelligence separately, and it is difficult to understand the… ▽ More Artificial intelligence technology has been widely used in astronomy, and new artificial intelligence technologies and application scenarios are constantly emerging. There have been a large number of papers reviewing the application of artificial intelligence technology in astronomy. However, relevant articles seldom mention telescope intelligence separately, and it is difficult to understand the current development status and research hotspots of telescope intelligence from these papers. This paper combines the development history of artificial intelligence technology and the difficulties of critical technologies of telescopes, comprehensively introduces the development and research hotspots of telescope intelligence, then conducts statistical analysis on various research directions of telescope intelligence and defines the research directions' merits. All kinds of research directions are evaluated, and the research trend of each telescope's intelligence is pointed out. Finally, according to the advantages of artificial intelligence technology and the development trend of telescopes, future research hotspots of telescope intelligence are given. △ Less

Submitted 16 January, 2024; v1 submitted 29 June, 2023; originally announced June 2023.

Comments: 41 pages, 10 figure, for questions or comments, please email tzhu@niaot.ac.cn

ACM Class: J.7

Hydrodynamic simulations of the Disk of Gas Around Supermassive black holes (HDGAS) -I; Molecular Gas Dynamics

Authors: Mojtaba Raouf, Serena Viti, S. García-Burillo, Alexander J. Richings, Joop schaye, Ashley Bemis, Folkert S. J. Nobels, Matteo Guainazzi, Ko-Yun Huang, Matthieu Schaller, Violette Impellizzeri, Jon Holdship

Abstract: We present hydrodynamic simulations of the interstellar medium (ISM) within the circumnuclear disk (CND) of a typical AGN-dominated galaxy influenced by mechanical feedback from an active galactic nucleus(AGN). The simulations are coupled with the CHIMES non-equilibrium chemistry network to treat the radiative-cooling and AGN-heating. A focus is placed on the central 100 pc scale where AGN outflow… ▽ More We present hydrodynamic simulations of the interstellar medium (ISM) within the circumnuclear disk (CND) of a typical AGN-dominated galaxy influenced by mechanical feedback from an active galactic nucleus(AGN). The simulations are coupled with the CHIMES non-equilibrium chemistry network to treat the radiative-cooling and AGN-heating. A focus is placed on the central 100 pc scale where AGN outflows are coupled to the ISM and constrained by observational Seyfert-2 galaxies. AGN-feedback models are implemented with different wind-velocity and mass-loading factors. We post-process the simulation snapshots with a radiative-transfer code to obtain the molecular emission lines. We find that the inclusion of an AGN promotes the formation of CO in clumpy and dense regions surrounding supermassive-blackholes (SMBH). The CO(1-0) intensity maps ($<$6 Myr) in the CND seem to match well with observations of NGC 1068 with a best match for a model with 5000 $\rm km/s$ wind-velocity and a high mass-loading factor. We attempt to discern between competing explanations for the apparent counter-rotating gas disk in the NGC 1068 through an analysis of kinematic maps of the CO line emission. We suggest that mechanical AGN-feedback could explain the alignment-stability of position-angle across the different CND radii around the SMBH through momentum and energy loading of the wind. It is the wind-velocity that drives the disk out of alignment on a 100 pc scale for a long period of time. The position-velocity diagrams are in broad agreement with the predicted Keplerian rotation-curve in the model without-AGN, but the AGN models exhibit a larger degree of scatter, in better agreement with NGC 1068 observations. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: 16 pages, 13 figures. Accepted for publication in MNRAS

JUNO sensitivity to the annihilation of MeV dark matter in the galactic halo

Authors: JUNO Collaboration, Angel Abusleme, Thomas Adam, Shakeel Ahmad, Rizwan Ahmed, Sebastiano Aiello, Muhammad Akram, Abid Aleem, Tsagkarakis Alexandros, Fengpeng An, Qi An, Giuseppe Andronico, Nikolay Anfimov, Vito Antonelli, Tatiana Antoshkina, Burin Asavapibhop, João Pedro Athayde Marcondes de André, Didier Auguste, Weidong Bai, Nikita Balashov, Wander Baldini, Andrea Barresi, Davide Basilico, Eric Baussan, Marco Bellato , et al. (581 additional authors not shown)

Abstract: We discuss JUNO sensitivity to the annihilation of MeV dark matter in the galactic halo via detecting inverse beta decay reactions of electron anti-neutrinos resulting from the annihilation. We study possible backgrounds to the signature, including the reactor neutrinos, diffuse supernova neutrino background, charged- and neutral-current interactions of atmospheric neutrinos, backgrounds from muon… ▽ More We discuss JUNO sensitivity to the annihilation of MeV dark matter in the galactic halo via detecting inverse beta decay reactions of electron anti-neutrinos resulting from the annihilation. We study possible backgrounds to the signature, including the reactor neutrinos, diffuse supernova neutrino background, charged- and neutral-current interactions of atmospheric neutrinos, backgrounds from muon-induced fast neutrons and cosmogenic isotopes. A fiducial volume cut, as well as the pulse shape discrimination and the muon veto are applied to suppress the above backgrounds. It is shown that JUNO sensitivity to the thermally averaged dark matter annihilation rate in 10 years of exposure would be significantly better than the present-day best limit set by Super-Kamiokande and would be comparable to that expected by Hyper-Kamiokande. △ Less

Submitted 13 September, 2023; v1 submitted 15 June, 2023; originally announced June 2023.

Comments: 25 pages, 9 figures, matches the publised version

Journal ref: JCAP 09 (2023) 001

The hadronic equation of state of HESS J1731-347 from the relativistic mean-field model with tensor coupling

Authors: Kaixuan Huang, Jinniu Hu, Ying Zhang, Hong Shen

Abstract: A recent report has identified a central compact object (CCO) within the supernova remnant HESS J1731-347, with a mass and radius of $M=0.77^{+0.20}_{-0.17}M{\odot}$ and $R=10.4^{+0.86}_{-0.78}$ km, respectively. To investigate this light compact star, a density-dependent relativistic mean-field (DDRMF) model, specifically the DDVT model, has been employed. The DDVT model incorporates tensor coupl… ▽ More A recent report has identified a central compact object (CCO) within the supernova remnant HESS J1731-347, with a mass and radius of $M=0.77^{+0.20}_{-0.17}M{\odot}$ and $R=10.4^{+0.86}_{-0.78}$ km, respectively. To investigate this light compact star, a density-dependent relativistic mean-field (DDRMF) model, specifically the DDVT model, has been employed. The DDVT model incorporates tensor couplings of vector mesons, which {can} successfully describe the properties of finite nuclei, such as charge radius, binding energy, and spin-orbit splitting. The introduction of tensor coupling reduces the influence of scalar mesons and generates a softer equation of state (EOS) in the outer core of the neutron star. Moreover, it has been found that the crust segment plays a crucial role in reproducing the mass-radius relation of HESS J1731-347, indicating a preference for a soft crust EOS. By manipulating the coupling strength of the isovector meson in the DDVT parameter set, a reasonable hadronic EOS has been obtained, satisfying the constraints from the gravitational-wave signal GW170817, the simultaneous mass-radius measurements from the NICER collaboration, and the properties of finite nuclei. Notably, the mass-radius relations derived from this hadronic EOS also accurately describe the observables of HESS J1731-347. Therefore, based on our estimation, the CCO in HESS J1731-347 may represent the lightest known neutron star. △ Less

Submitted 8 June, 2023; originally announced June 2023.

Comments: 20 pages, 4 figures, 4 tables, suggestions and comments are welcome!

Machine learning-based seeing estimation and prediction using multi-layer meteorological data at Dome A, Antarctica

Authors: Xu Hou, Yi Hu, Fujia Du, Michael C. B. Ashley, Chong Pei, Zhaohui Shang, Bin Ma, Erpeng Wang, Kang Huang

Abstract: Atmospheric seeing is one of the most important parameters for evaluating and monitoring an astronomical site. Moreover, being able to predict the seeing in advance can guide observing decisions and significantly improve the efficiency of telescopes. However, it is not always easy to obtain long-term and continuous seeing measurements from a standard instrument such as differential image motion mo… ▽ More Atmospheric seeing is one of the most important parameters for evaluating and monitoring an astronomical site. Moreover, being able to predict the seeing in advance can guide observing decisions and significantly improve the efficiency of telescopes. However, it is not always easy to obtain long-term and continuous seeing measurements from a standard instrument such as differential image motion monitor (DIMM), especially for those unattended observatories with challenging environments such as Dome A, Antarctica. In this paper, we present a novel machine learning-based framework for estimating and predicting seeing at a height of 8 m at Dome A, Antarctica, using only the data from a multi-layer automated weather station (AWS). In comparison with DIMM data, our estimate has a root mean square error (RMSE) of 0.18 arcsec, and the RMSE of predictions 20 minutes in the future is 0.12 arcsec for the seeing range from 0 to 2.2 arcsec. Compared with the persistence, where the forecast is the same as the last data point, our framework reduces the RMSE by 37 percent. Our method predicts the seeing within a second of computing time, making it suitable for real-time telescope scheduling. △ Less

Submitted 7 April, 2023; originally announced April 2023.

Comments: 13 pages, 14 figures, accepted for publication in Astronomy and Computing

Reconstructing the shock history in the CMZ of NGC 253 with ALCHEMI

Authors: K. -Y. Huang, S. Viti, J. Holdship, J. G. Mangum, S. Martín, N. Harada, S. Muller, K. Sakamoto, K. Tanaka, Y. Yoshimura, R. Herrero-Illana, D. S. Meier, E. Behrens, P. P. van der Werf, C. Henkel, S. García-Burillo, V. M. Rivilla, K. L. Emig, L. Colzi, P. K. Humire, R. Aladro, M. Bouvier

Abstract: HNCO and SiO are well known shock tracers and have been observed in nearby galaxies, including the nearby (D=3.5 Mpc) starburst galaxy NGC 253. The simultaneous detection of these two species in regions where the star formation rate is high may be used to study the shock history of the gas. We perform a multi-line molecular study using these two shock tracers (SiO and HNCO) with the aim of charact… ▽ More HNCO and SiO are well known shock tracers and have been observed in nearby galaxies, including the nearby (D=3.5 Mpc) starburst galaxy NGC 253. The simultaneous detection of these two species in regions where the star formation rate is high may be used to study the shock history of the gas. We perform a multi-line molecular study using these two shock tracers (SiO and HNCO) with the aim of characterizing the gas properties. We also explore the possibility of reconstructing the shock history in NGC 253's Central Molecular Zone (CMZ). Six SiO transitions and eleven HNCO transitions were imaged at high resolution $1''.6$ (28 pc) with the Atacama Large Millimeter/submillimeter Array (ALMA) as part of the ALCHEMI Large Programme. Both non-LTE radiative transfer analysis and chemical modelling were performed in order to characterize the gas properties, and to investigate the chemical origin of the emission. The non-LTE radiative transfer analysis coupled with Bayesian inference shows clear evidence that the gas traced by SiO has different densities and temperatures than that traced by HNCO, with an indication that shocks are needed to produce both species. Chemical modelling further confirms such a scenario and suggests that fast and slow shocks are responsible for SiO and HNCO production, respectively, in most GMCs. We are also able to infer the physical characteristics of the shocks traced by SiO and HNCO for each GMC. Radiative transfer and chemical analysis of the SiO and HNCO in the CMZ of NGC 253 reveal a complex picture whereby most of the GMCs are subjected to shocks. We speculate on the possible shock scenarios responsible for the observed emission and provide potential history and timescales for each shock scenario. Higher spatial resolution observations of these two species are required in order to quantitatively differentiate between scenarios. △ Less

Submitted 22 March, 2023; originally announced March 2023.

Journal ref: A&A 675, A151 (2023)

Simultaneous Radio and Optical Polarimetry of GRB 191221B Afterglow

Authors: Yuji Urata, Kenji Toma, Stefano Covino, Klaas Wiersema, Kuiyun Huang, Jiro Shimoda, Asuka Kuwata, Sota Nagao, Keiichi Asada, Hiroshi Nagai, Satoko Takahashi, Chao-En Chung, Glen Petitpas, Kazutaka Yamaoka, Luca Izzo, Johan Fynbo, Antonio de Ugarte Postigo, Maryam Arabsalmani, Makoto Tashiro

Abstract: Gamma-ray bursts (GRBs) are the most luminous transients in the universe and are utilized as probes of early stars, gravitational wave counterparts, and collisionless shock physics. In spite of studies on polarimetry of GRBs in individual wavelengths that characterized intriguing properties of prompt emission and afterglow, no coordinated multi-wavelength measurements have yet been performed. Here… ▽ More Gamma-ray bursts (GRBs) are the most luminous transients in the universe and are utilized as probes of early stars, gravitational wave counterparts, and collisionless shock physics. In spite of studies on polarimetry of GRBs in individual wavelengths that characterized intriguing properties of prompt emission and afterglow, no coordinated multi-wavelength measurements have yet been performed. Here, we report the first coordinated simultaneous polarimetry in the optical and radio bands for the afterglow associated with the typical long GRB 191221B. Our observations successfully caught the radio emission, which is not affected by synchrotron self-absorption, and show that the emission is depolarized in the radio band compared to the optical one. Our simultaneous polarization angle measurement and temporal polarization monitoring indicate the existence of cool electrons that increase the estimate of jet kinetic energy by a factor of $>$ 4 for this GRB afterglow. Further coordinated multi-wavelength polarimetric campaigns would improve our understanding of the total jet energies and magnetic field configurations in the emission regions of various types of GRBs, which are required to comprehend the mass scales of their progenitor systems and the physics of collisionless shocks. △ Less

Submitted 9 December, 2022; originally announced December 2022.

Comments: Author's version. 41 pages, 9 figures, 4 tables

Model Independent Approach of the JUNO $^8$B Solar Neutrino Program

Authors: JUNO Collaboration, Jie Zhao, Baobiao Yue, Haoqi Lu, Yufeng Li, Jiajie Ling, Zeyuan Yu, Angel Abusleme, Thomas Adam, Shakeel Ahmad, Rizwan Ahmed, Sebastiano Aiello, Muhammad Akram, Abid Aleem, Tsagkarakis Alexandros, Fengpeng An, Qi An, Giuseppe Andronico, Nikolay Anfimov, Vito Antonelli, Tatiana Antoshkina, Burin Asavapibhop, João Pedro Athayde Marcondes de André, Didier Auguste, Weidong Bai , et al. (579 additional authors not shown)

Abstract: The physics potential of detecting $^8$B solar neutrinos will be exploited at the Jiangmen Underground Neutrino Observatory (JUNO), in a model independent manner by using three distinct channels of the charged-current (CC), neutral-current (NC) and elastic scattering (ES) interactions. Due to the largest-ever mass of $^{13}$C nuclei in the liquid-scintillator detectors and the {expected} low backg… ▽ More The physics potential of detecting $^8$B solar neutrinos will be exploited at the Jiangmen Underground Neutrino Observatory (JUNO), in a model independent manner by using three distinct channels of the charged-current (CC), neutral-current (NC) and elastic scattering (ES) interactions. Due to the largest-ever mass of $^{13}$C nuclei in the liquid-scintillator detectors and the {expected} low background level, $^8$B solar neutrinos would be observable in the CC and NC interactions on $^{13}$C for the first time. By virtue of optimized event selections and muon veto strategies, backgrounds from the accidental coincidence, muon-induced isotopes, and external backgrounds can be greatly suppressed. Excellent signal-to-background ratios can be achieved in the CC, NC and ES channels to guarantee the $^8$B solar neutrino observation. From the sensitivity studies performed in this work, we show that JUNO, with ten years of data, can reach the {1$σ$} precision levels of 5%, 8% and 20% for the $^8$B neutrino flux, $\sin^2θ_{12}$, and $Δm^2_{21}$, respectively. It would be unique and helpful to probe the details of both solar physics and neutrino physics. In addition, when combined with SNO, the world-best precision of 3% is expected for the $^8$B neutrino flux measurement. △ Less

Submitted 6 March, 2024; v1 submitted 15 October, 2022; originally announced October 2022.

Comments: 19 pages, 7 figures, accepted version to appear in The Astrophysical Journal. Yufeng Li and Jiajie Ling are corresponding authors

Journal ref: Astrophysical Journal 965 (2024) 122

Strong Gravitational Lensing Parameter Estimation with Vision Transformer

Authors: Kuan-Wei Huang, Geoff Chih-Fan Chen, Po-Wen Chang, Sheng-Chieh Lin, Chia-Jung Hsu, Vishal Thengane, Joshua Yao-Yu Lin

Abstract: Quantifying the parameters and corresponding uncertainties of hundreds of strongly lensed quasar systems holds the key to resolving one of the most important scientific questions: the Hubble constant ($H_{0}$) tension. The commonly used Markov chain Monte Carlo (MCMC) method has been too time-consuming to achieve this goal, yet recent work has shown that convolution neural networks (CNNs) can be a… ▽ More Quantifying the parameters and corresponding uncertainties of hundreds of strongly lensed quasar systems holds the key to resolving one of the most important scientific questions: the Hubble constant ($H_{0}$) tension. The commonly used Markov chain Monte Carlo (MCMC) method has been too time-consuming to achieve this goal, yet recent work has shown that convolution neural networks (CNNs) can be an alternative with seven orders of magnitude improvement in speed. With 31,200 simulated strongly lensed quasar images, we explore the usage of Vision Transformer (ViT) for simulated strong gravitational lensing for the first time. We show that ViT could reach competitive results compared with CNNs, and is specifically good at some lensing parameters, including the most important mass-related parameters such as the center of lens $θ_{1}$ and $θ_{2}$, the ellipticities $e_1$ and $e_2$, and the radial power-law slope $γ'$. With this promising preliminary result, we believe the ViT (or attention-based) network architecture can be an important tool for strong lensing science for the next generation of surveys. The open source of our code and data is in \url{https://github.com/kuanweih/strong_lensing_vit_resnet}. △ Less

Submitted 8 October, 2022; originally announced October 2022.

Comments: Accepted by ECCV 2022 AI for Space Workshop

Neutron Tagging following Atmospheric Neutrino Events in a Water Cherenkov Detector

Authors: K. Abe, Y. Haga, Y. Hayato, K. Hiraide, K. Ieki, M. Ikeda, S. Imaizumi, K. Iyogi, J. Kameda, Y. Kanemura, Y. Kataoka, Y. Kato, Y. Kishimoto, S. Miki, S. Mine, M. Miura, T. Mochizuki, S. Moriyama, Y. Nagao, M. Nakahata, T. Nakajima, Y. Nakano, S. Nakayama, T. Okada, K. Okamoto , et al. (281 additional authors not shown)

Abstract: We present the development of neutron-tagging techniques in Super-Kamiokande IV using a neural network analysis. The detection efficiency of neutron capture on hydrogen is estimated to be 26%, with a mis-tag rate of 0.016 per neutrino event. The uncertainty of the tagging efficiency is estimated to be 9.0%. Measurement of the tagging efficiency with data from an Americium-Beryllium calibration agr… ▽ More We present the development of neutron-tagging techniques in Super-Kamiokande IV using a neural network analysis. The detection efficiency of neutron capture on hydrogen is estimated to be 26%, with a mis-tag rate of 0.016 per neutrino event. The uncertainty of the tagging efficiency is estimated to be 9.0%. Measurement of the tagging efficiency with data from an Americium-Beryllium calibration agrees with this value within 10%. The tagging procedure was performed on 3,244.4 days of SK-IV atmospheric neutrino data, identifying 18,091 neutrons in 26,473 neutrino events. The fitted neutron capture lifetime was measured as 218 \pm 9 μs. △ Less

Submitted 20 September, 2022; v1 submitted 18 September, 2022; originally announced September 2022.

Journal ref: JINST 17 P10029 (2022)

Tracing Interstellar Heating: An ALCHEMI Measurement of the HCN Isomers in NGC 253

Authors: Erica Behrens, Jeffrey G. Mangum, Jonathan Holdship, Serena Viti, Nanase Harada, Sergio Martin, Kazushi Sakamoto, Sebastien Muller, Kunihiko Tanaka, Kouichiro Nakanishi, Ruben Herrero-Illana, Yuki Yoshimura, Rebeca Aladro, Laura Colzi, Kimberly L. Emig, Christian Henkel, Ko-Yun Huang, P. K. Humire, David S. Meier, Victor M. Rivilla

Abstract: We analyze HCN and HNC emission in the nearby starburst galaxy NGC 253 to investigate its effectiveness in tracing heating processes associated with star formation. This study uses multiple HCN and HNC rotational transitions observed using ALMA via the ALCHEMI Large Program. To understand the conditions and associated heating mechanisms within NGC 253's dense gas, we employ Bayesian nested samplin… ▽ More We analyze HCN and HNC emission in the nearby starburst galaxy NGC 253 to investigate its effectiveness in tracing heating processes associated with star formation. This study uses multiple HCN and HNC rotational transitions observed using ALMA via the ALCHEMI Large Program. To understand the conditions and associated heating mechanisms within NGC 253's dense gas, we employ Bayesian nested sampling techniques applied to chemical and radiative transfer models which are constrained using our HCN and HNC measurements. We find that the volume density $n_{\text{H}_{2}}$ and cosmic ray ionization rate (CRIR) $ζ$ are enhanced by about an order of magnitude in the galaxy's central regions as compared to those further from the nucleus. In NGC 253's central GMCs, where observed HCN/HNC abundance ratios are lowest, $n \sim 10^{5.5}$ cm$^{-3}$ and $ζ\sim 10^{-12}$ s$^{-1}$ (greater than $10^4$ times the average Galactic rate). We find a positive correlation in the association of both density and CRIR with the number of star formation-related heating sources (supernova remnants, HII regions, and super hot cores) located in each GMC, as well as a correlation between CRIRs and supernova rates. Additionally, we see an anticorrelation between the HCN/HNC ratio and CRIR, indicating that this ratio will be lower in regions where $ζ$ is higher. Though previous studies suggested HCN and HNC may reveal strong mechanical heating processes in NGC 253's CMZ, we find cosmic ray heating dominates the heating budget, and mechanical heating does not play a significant role in the HCN and HNC chemistry. △ Less

Submitted 8 November, 2022; v1 submitted 13 September, 2022; originally announced September 2022.

Comments: 33 pages, 23 figures, accepted for publication by the Astrophysical Journal

Methanol masers in NGC 253 with ALCHEMI

Authors: P. K. Humire, C. Henkel, A. Hernández-Gómez, S. Martín, J. Mangum, N. Harada, S. Muller, K. Sakamoto, K. Tanaka, Y. Yoshimura, K. Nakanishi, S. Mühle, R. Herrero-Illana, D. S. Meier, E. Caux, R. Aladro, R. Mauersberger, S. Viti, L. Colzi, V. M. Rivilla, M. Gorski, K. M. Menten, K. -Y. Huang, S. Aalto, P. P. van der Werf , et al. (1 additional authors not shown)

Abstract: Context: Methanol masers of Class I (collisionally-pumped) and Class II (radiatively-pumped) have been studied in great detail in our Galaxy in a variety of astrophysical environments such as shocks and star-forming regions and are helpful to analyze the properties of the dense interstellar medium. However, the study of methanol masers in external galaxies is still in its infancy. Aims: Our main g… ▽ More Context: Methanol masers of Class I (collisionally-pumped) and Class II (radiatively-pumped) have been studied in great detail in our Galaxy in a variety of astrophysical environments such as shocks and star-forming regions and are helpful to analyze the properties of the dense interstellar medium. However, the study of methanol masers in external galaxies is still in its infancy. Aims: Our main goal is to search for methanol masers in the central molecular zone (CMZ; inner 500 pc) of the nearby starburst galaxy NGC 253. Methods: Covering a frequency range between 84 and 373 GHz ($λ$ = 3.6 to 0.8 mm) at high angular (1.6"$\sim$27 pc) and spectral ($\sim$8--9 km s$^{-1}$) resolution with the ALMA large program ALCHEMI, we have probed different regions across the CMZ of NGC 253. In order to look for methanol maser candidates, we employed the rotation diagram method and a set of radiative transfer models. Results: We detect for the first time masers above 84 GHz in NGC 253, covering an ample portion of the $J_{-1}\rightarrow(J-$ 1)$_{0}-E$ line series (at 84, 132, 229, and 278 GHz) and the $J_{0}\rightarrow(J-$ 1)$_{1}-A$ series (at 95, 146, and 198 GHz). This confirms the presence of the Class I maser line at 84 GHz, already reported but now being detected in more than one location. For the $J_{-1}\rightarrow(J-$ 1)$_{0}-E$ line series, we observe a lack of Class I maser candidates in the central star-forming disk. Conclusions: The physical conditions for maser excitation in the $J_{-1}\rightarrow(J-$ 1)$_{0}-E$ line series can be weak shocks and cloud-cloud collisions as suggested by shock tracers (SiO and HNCO) in bi-symmetric shock/active regions located in the outskirts of the CMZ. On the other hand, the presence of photodissociation regions due to a high star-formation rate would be needed to explain the lack of Class I masers in the very central regions. △ Less

Submitted 6 May, 2022; originally announced May 2022.

Comments: Accepted for publication in A&A. 29 pages, 17 figures (4 in Appendix)

The chemical footprint of AGN feedback in the outflowing circumnuclear disk of NGC 1068

Authors: K. -Y. Huang, S. Viti, J. Holdship, S. García-Burillo, K. Kohno, A. Taniguchi, S. Martín, R. Aladro, A. Fuente, M. Sánchez-García

Abstract: In the nearby (D=14 Mpc) AGN-starburst composite galaxy NGC 1068, it has been found that the molecular gas in the Circum-nuclear Disk (CND) is outflowing, which is a manifestation of ongoing AGN feedback. The outflowing gas has a large spread of velocities, which likely drive different shock chemistry signatures at different locations in the CND. We perform a multi-line molecular study using two s… ▽ More In the nearby (D=14 Mpc) AGN-starburst composite galaxy NGC 1068, it has been found that the molecular gas in the Circum-nuclear Disk (CND) is outflowing, which is a manifestation of ongoing AGN feedback. The outflowing gas has a large spread of velocities, which likely drive different shock chemistry signatures at different locations in the CND. We perform a multi-line molecular study using two shock tracers, SiO and HNCO, with the aim to determine the gas properties traced by these two species, and explore the possibility of reconstructing the shock history in the CND. Five SiO transitions and three HNCO transitions were imaged at high resolution $0''.5-0''.8$ with the Atacama Large Millimeter/submillimeter Array (ALMA). We performed both LTE and non-LTE radiative transfer analysis coupled with Bayesian inference process in order to characterize the gas properties, such as molecular gas density and gas temperature. We found clear evidence of chemical differentiation between SiO and HNCO, with the SiO/HNCO ratio ranging from greater than one on the east of CND to lower than one on the west side. The non-LTE radiative transfer analysis coupled with Bayesian inference confirms that the gas traced by SiO has different densities - and possibly temperatures - than that traced by HNCO. We find that SiO traces gas affected by fast shocks while the gas traced by HNCO is either just affected by slow shocks or not shocked at all. A distinct differentiation between SiO and HNCO has been revealed in our observations and the further analysis of the gas properties traced by both species, which confirms the results from previous chemical modelings. △ Less

Submitted 10 February, 2022; originally announced February 2022.

Comments: Accepted for publication in A&A

Journal ref: A&A 666, A102 (2022)

Energy functions of fast radio bursts derived from the first CHIME/FRB catalogue

Authors: Tetsuya Hashimoto, Tomotsugu Goto, Bo Han Chen, Simon C. -C. Ho, Tiger Y. -Y. Hsiao, Yi Hang Valerie Wong, Alvina Y. L. On, Seong Jin Kim, Ece Kilerci-Eser, Kai-Chun Huang, Daryl Joe D. Santos, Shotaro Yamasaki

Abstract: Fast radio bursts (FRBs) are mysterious millisecond pulses in radio, most of which originate from distant galaxies. Revealing the origin of FRBs is becoming central in astronomy. The redshift evolution of the FRB energy function, i.e., the number density of FRB sources as a function of energy, provides important implications for the FRB progenitors. Here we show the energy functions of FRBs select… ▽ More Fast radio bursts (FRBs) are mysterious millisecond pulses in radio, most of which originate from distant galaxies. Revealing the origin of FRBs is becoming central in astronomy. The redshift evolution of the FRB energy function, i.e., the number density of FRB sources as a function of energy, provides important implications for the FRB progenitors. Here we show the energy functions of FRBs selected from the recently released Canadian Hydrogen Intensity Mapping Experiment (CHIME) catalogue using the $V_{\rm max}$ method. The $V_{\rm max}$ method allows us to measure the redshift evolution of the energy functions as it is without any prior assumption on the evolution. We use a homogeneous sample of 164 non-repeating FRB sources, which are about one order of magnitude larger than previously investigated samples. The energy functions of non-repeating FRBs show Schechter function-like shapes at $z\lesssim1$. The energy functions and volumetric rates of non-repeating FRBs decrease towards higher redshifts similar to the cosmic stellar-mass density evolution: there is no significant difference between the non-repeating FRB rate and cosmic stellar-mass density evolution with a 1\% significance threshold, whereas the cosmic star-formation rate scenario is rejected with a more than 99\% confidence level. Our results indicate that the event rate of non-repeating FRBs is likely controlled by old populations rather than young populations which are traced by the cosmic star-formation rate density. This suggests old populations such as old neutron stars and black holes as more likely progenitors of non-repeating FRBs. △ Less

Submitted 10 January, 2022; originally announced January 2022.

Comments: Accepted for publication in MNRAS. A summary video is available at https://youtu.be/qbUctrFHLAs

Identifying RR Lyrae in the ZTF DR3 dataset

Authors: Kuan-Wei Huang, Sergey E. Koposov

Abstract: We present a RR Lyrae (RRL) catalogue based on the combination of the third data release of the Zwicky Transient Facility (ZTF DR3) and \textit{Gaia} EDR3. We use a multi-step classification pipeline relying on the Fourier decomposition fitting to the multi-band ZTF light curves and random forest classification. The resulting catalogue contains 71,755 RRLs with period and light curve parameter mea… ▽ More We present a RR Lyrae (RRL) catalogue based on the combination of the third data release of the Zwicky Transient Facility (ZTF DR3) and \textit{Gaia} EDR3. We use a multi-step classification pipeline relying on the Fourier decomposition fitting to the multi-band ZTF light curves and random forest classification. The resulting catalogue contains 71,755 RRLs with period and light curve parameter measurements and has completeness of 0.92 and purity of 0.92 with respect to the SOS \textit{Gaia} DR2 RRLs. The catalogue covers the Northern sky with declination $\geq -28^\circ$, its completeness is $\gtrsim 0.8$ for heliocentric distance $\leq 80$~kpc, and the most distant RRL at 132~kpc. Compared with several other RRL catalogues covering the Northern sky, our catalogue has more RRLs around the Galactic halo and is more complete at low Galactic latitude areas. Analysing the spatial distribution of RRL in the catalogue reveals the previously known major over-densities of the Galactic halo, such as the Virgo over-density and the Hercules-Aquila Cloud, with some evidence of an association between the two. We also analyse the Oosterhoff fraction differences throughout the halo, comparing it with the density distribution, finding increasing Oosterhoff I fraction at the elliptical radii between 16 and 32 kpc and some evidence of different Oosterhoff fractions across various halo substructures. △ Less

Submitted 11 December, 2021; originally announced December 2021.

ALMA Host galaxy observation of the off-axis Gamma-Ray Burst XRF020903

Authors: Jheng-Cyun Chen, Yuji Urata, Kuiyun Huang

Abstract: We investigated the radio properties of the host galaxy of X-ray flash, XRF020903, which is the best example for investigating of the off-axis origin of gamma-ray bursts(GRBs). Dust continuum at 233 GHz and CO are observed using the Atacama Large millimeter/submillimeter array. The molecular gas mass derived by applying the metalicity-dependent CO-to-H$_{2}$ conversion factor matches the global tr… ▽ More We investigated the radio properties of the host galaxy of X-ray flash, XRF020903, which is the best example for investigating of the off-axis origin of gamma-ray bursts(GRBs). Dust continuum at 233 GHz and CO are observed using the Atacama Large millimeter/submillimeter array. The molecular gas mass derived by applying the metalicity-dependent CO-to-H$_{2}$ conversion factor matches the global trend along the redshift and stellar mass of the GRB host galaxies. The estimated gas depletion timescale (pertaining to the potential critical characteristics of GRB host galaxies) is equivalent to those of GRBs and super-luminous supernova hosts in the same redshift range. These properties of the XRF020903 host galaxy observed in radio resemble those of GRB host galaxies, thereby supporting the identical origin of XRF020903 and GRBs. △ Less

Submitted 19 June, 2021; v1 submitted 12 May, 2021; originally announced May 2021.

Comments: 14 pages, 9 figures, 1 tables, ApJ in press

The possibility of the secondary object in GW190814 as a neutron star

Authors: Kaixuan Huang, Jinniu Hu, Ying Zhang, Hong Shen

Abstract: A compact object was observed with a mass $2.50-2.67~M_\odot$ by LIGO Scientific and Virgo collaborations (LVC) in GW190814, which provides a great challenge to the investigations into the supranuclear matter. To study this object, properties of neutron star are systematically calculated within the latest density-dependent relativistic mean-field (DDRMF) parameterizations, which are determined by… ▽ More A compact object was observed with a mass $2.50-2.67~M_\odot$ by LIGO Scientific and Virgo collaborations (LVC) in GW190814, which provides a great challenge to the investigations into the supranuclear matter. To study this object, properties of neutron star are systematically calculated within the latest density-dependent relativistic mean-field (DDRMF) parameterizations, which are determined by the ground state properties of spherical nuclei. The maximum masses of neutron star calculated by DD-MEX and DD-LZ1 sets can be around $2.55~M_\odot$ with quite stiff equations of state generated by their strong repulsive contributions from vector potentials at high densities. Their maximum speeds of sound $c_s/c$ are smaller than $\sqrt{0.8}$ at the center of neutron star and the dimensionless tidal deformabilities at $1.4~M_\odot$ are less than $800$. Furthermore, the radii of $1.4 ~M_\odot$ also satisfy the constraint from the observation of mass-radius simultaneous measurements (NICER). Therefore, we conclude that one cannot exclude the possibility of the secondary object in GW190814 as a neutron star composed of hadron matter from DDRMF models. △ Less

Submitted 10 August, 2020; originally announced August 2020.

Comments: 13 pages, 8 figures, 3 tables

Mapping circumstellar magnetic fields of late-type evolved stars with the Goldreich-Kylafis effect: CARMA observations at $λ1.3$ mm of R Crt and R Leo

Authors: K. -Y. Huang, A. J. Kemball, W. H. T. Vlemmings, S. -P. Lai, L. Yang, I. Agudo

Abstract: Mapping magnetic fields is the key to resolving what remains an unclear physical picture of circumstellar magnetic fields in late-type evolved stars. Observations of linearly polarized emission from thermal molecular line transitions due to the Goldreich-Kylafis (G-K) effect provides valuable insight into the magnetic field geometry in these sources that is complementary to other key studies. In t… ▽ More Mapping magnetic fields is the key to resolving what remains an unclear physical picture of circumstellar magnetic fields in late-type evolved stars. Observations of linearly polarized emission from thermal molecular line transitions due to the Goldreich-Kylafis (G-K) effect provides valuable insight into the magnetic field geometry in these sources that is complementary to other key studies. In this paper, we present the detection of spectral-line polarization from both the thermal $J=2-1$ CO line and the $v=1, J=5-4$ SiO maser line toward two thermal-pulsating (TP-) AGB stars, R Crt and R Leo. The observed fractional linear polarization in the CO emission is measured as $m_l\sim 3.1\%$ and $m_l\sim9.7\%$ for R Crt and R Leo respectively. A circumstellar envelope (CSE) model profile and the associated parameters are estimated and used as input to a more detailed modeling of the predicted linear polarization expected from the G-K effect. The observed thermal line polarization level is consistent with the predicted results from the G-K model for R Crt; additional effects need to be considered for R Leo. △ Less

Submitted 1 July, 2020; originally announced July 2020.

Search for globular clusters associated with the Milky Way dwarf galaxies using Gaia DR2

Authors: Kuan-Wei Huang, Sergey E. Koposov

Abstract: We report the result of searching for globular clusters (GCs) around 55 Milky Way satellite dwarf galaxies within the distance of 450 kpc from the Galactic Center except for the Large and Small Magellanic Clouds and the Sagittarius dwarf. For each dwarf, we analyze the stellar distribution of sources in Gaia DR2, selected by magnitude, proper motion, and source morphology. Using the kernel density… ▽ More We report the result of searching for globular clusters (GCs) around 55 Milky Way satellite dwarf galaxies within the distance of 450 kpc from the Galactic Center except for the Large and Small Magellanic Clouds and the Sagittarius dwarf. For each dwarf, we analyze the stellar distribution of sources in Gaia DR2, selected by magnitude, proper motion, and source morphology. Using the kernel density estimation of stellar number counts, we identify eleven possible GC candidates. Crossed-matched with existing imaging data, all eleven objects are known either GCs or galaxies and only Fornax GC 1-6 among them are associated with the targeted dwarf galaxy. Using simulated GCs, we calculate the GC detection limit $M_{\rm V}^{\rm lim}$ that spans the range from $M_{\rm V}^{\rm lim} \sim -7$ for distant dwarfs to $M_{\rm V}^{\rm lim} \sim 0$ for nearby systems. Assuming a Gaussian GC luminosity function, we compute that the completeness of the GC search is above 90 percent for most dwarf galaxies. We construct the 90 percent credible intervals/upper limits on the GC specific frequency $S_{\rm N}$ of the MW dwarf galaxies: $12 < S_{\rm N} < 47$ for Fornax, $S_{\rm N} < 20$ for the dwarfs with $-12 < M_{\rm V} < -10$, $S_{\rm N} < 30$ for the dwarfs with $-10 < M_{\rm V} < -7$, and $S_{\rm N} < 90$ for the dwarfs with $M_{\rm V} > -7$. Based on $S_{\rm N}$, we derive the probability of galaxies hosting GCs given their luminosity, finding that the probability of galaxies fainter than $M_{\rm V} = -9$ to host GCs is lower than 0.1. △ Less

Submitted 2 November, 2020; v1 submitted 28 May, 2020; originally announced May 2020.

Orphan GRB afterglow searches with the Pan-STARRS1 COSMOS survey

Authors: Yun-Jing Huang, Yuji Urata, Kuiyun Huang, Kuei-sheng Lee, Meng-feng Tsai, Yuji Shirasaki, Marcin Sawicki, Stephane Arnouts, Thibaud Moutard, Stephen Gwyn, Wei-Hao Wang, Sebastien Foucaud, Keiichi Asada, Mark E. Huber, Richard Wainscoat, Kenneth C. Chambers

Abstract: We present the result of a search for orphan Gamma-Ray Burst (GRB) afterglows in the Pan-STARRS1 (PS1) COSMOS survey. There is extensive theoretical and observational evidence suggesting that GRBs are collimated jets; the direct observation of orphan GRB afterglows would further support this model. An optimal survey strategy is designed by coupling the PS1 survey with the Subaru/Hyper-Suprime-Cam… ▽ More We present the result of a search for orphan Gamma-Ray Burst (GRB) afterglows in the Pan-STARRS1 (PS1) COSMOS survey. There is extensive theoretical and observational evidence suggesting that GRBs are collimated jets; the direct observation of orphan GRB afterglows would further support this model. An optimal survey strategy is designed by coupling the PS1 survey with the Subaru/Hyper-Suprime-Cam (HSC) survey. The PS1 COSMOS survey, one of the survey fields in the PS1 Medium Deep Survey (PS1/MDS), searches a field of 7 deg$^2$ from December 2011 to January 2014, reaching a limiting magnitude R $\sim$ 23. The dense cadence of PS1/MDS is crucial for identifying transients, and the deep magnitude reached by the HSC survey (R $\sim$ 26) is important for evaluating potential GRB hosts. A transient classification method is employed to select potential orphan GRB afterglow candidates. After a thorough analysis of the transient and host galaxy properties, we conclude that there are no candidates in this survey field. The null result implies that the consideration of jet structures is essential for further orphan GRB afterglow surveys. △ Less

Submitted 12 July, 2020; v1 submitted 4 May, 2020; originally announced May 2020.

Comments: Matches published version; 16 pages, 7 figures, 4 tables

Journal ref: ApJ 897 69 (2020)

Spectroscopically Confirmed Lyman-Alpha Emitters from Redshift 5 to 7 Behind Ten Galaxy Cluster Lenses

Authors: Spencer Fuller, Brian Lemaux, Maruša Bradač, Austin Hoag, Kasper Schmidt, Kuang Huang, Victoria Strait, Charlotte Mason, Tommaso Treu, Laura Pentericci, Michele Trenti, Alaina Henry, Matt Malkan

Abstract: We present 36 spectroscopically confirmed intrinsically UV-faint Ly$α$ emitting galaxies from follow-up observations with Keck/DEIMOS of gravitationally lensed high-redshift candidates. Candidates were selected to be between $5\lesssim z \lesssim 7$ from photometric data using \textit{HST} and \textit{Spitzer} imaging surveys. We used photometric redshift information to perform an integrated photo… ▽ More We present 36 spectroscopically confirmed intrinsically UV-faint Ly$α$ emitting galaxies from follow-up observations with Keck/DEIMOS of gravitationally lensed high-redshift candidates. Candidates were selected to be between $5\lesssim z \lesssim 7$ from photometric data using \textit{HST} and \textit{Spitzer} imaging surveys. We used photometric redshift information to perform an integrated photometric redshift probability cut $>1\%$ between $5<z<7$ to construct a sample of 198 high-redshift objects. Our high-redshift sample spans intrinsic UV luminosities from a few $L^*$ down to $0.001L^*$. We identified 19 high-confidence detections of Ly$α$ and an additional 17 likely detections. We divided our sample into lower-redshift ($z\sim5.5$) and higher-redshift ($z\sim6.5$) bins and ran Monte Carlo trials, incorporating the strength of the Ly$α$ emission and the photometric redshift of the non-detections. Considering only objects where Ly$α$ could be detected at EW(Ly$α$)$>$25Å at $3σ$ at the fiducial depth of our survey, and only those galaxies with EW(Ly$α$)$>$25Å as true LAEs, and finally, only objects with $m_{AB}<26.8$, we found the LAE fraction to be flat, or modestly increase from 0.26$\pm0.04$ to 0.30$\pm0.04$. These values relative to those for lower-redshift samples are consistent with a rising LAE fraction with redshift out to $z\sim6$, but at $z\sim6.5$ there is some tension between our results and results from surveys at intrinsically brighter luminosities. We conclude intrinsically fainter galaxies have Ly$α$ emission, and there is a steep drop in the LAE fraction from our high-redshift sample at $z\sim6.5$ and from similar galaxies at $z\sim7.5$. This likely indicates we are witnessing the tail end of the epoch of reionization, as such a drop is not expected due to changes of intrinsic galaxy properties between these redshifts. △ Less

Submitted 20 February, 2020; originally announced February 2020.

Comments: 21 pages, 10 figures, submitted to ApJ

Two Component Jets of GRB160623A as Shocked Jet cocoon afterglow

Authors: Wei Ju Chen, Yuji Urata, Kuiyun Huang, Satoko Takahashi, Glen Petitpas, Keiichi Asada

Abstract: Two components of jets associated with the afterglow of the gamma-ray burst GRB 160623A were observed with multi-frequency observations including long-term monitoring in a sub-millimetre range (230 GHz) using the SMA. The observed light curves with temporal breaks suggests on the basis of the standard forward-shock synchrotron radiation model that the X-ray radiation is narrowly collimated with an… ▽ More Two components of jets associated with the afterglow of the gamma-ray burst GRB 160623A were observed with multi-frequency observations including long-term monitoring in a sub-millimetre range (230 GHz) using the SMA. The observed light curves with temporal breaks suggests on the basis of the standard forward-shock synchrotron radiation model that the X-ray radiation is narrowly collimated with an opening angle $θ_{n,j}<\sim6^{\circ}$ whereas the radio radiation originated from wider jets ($\sim27^{\circ}$). The temporal and spectral evolutions of the radio afterglow agree with those expected from a synchrotron radiation modelling with typical physical parameters except for the fact that the observed wide jet opening angle for the radio emission is significantly larger than the theoretical maximum opening angle. By contrast, the opening angle of the X-ray afterglow is consistent with the typical value of GRB jets. Since the theory of the relativistic cocoon afterglow emission is similar to that of a regular afterglow with an opening angle of $\sim30^{\circ}$, the observed radio emission can be interpreted as the shocked jet cocoon emission. This result therefore indicates that the two components of the jets observed in the GRB 160623A afterglow is caused by the jet and the shocked jet cocoon afterglows. △ Less

Submitted 16 February, 2020; originally announced February 2020.

Comments: Accepted to ApJ Letters; 8 pages, 4 figures, 1 table

RELICS: A Very Large ($θ_{E}\sim40"$) Cluster Lens -- RXC J0032.1+1808

Authors: Ana Acebron, Adi Zitrin, Dan Coe, Guillaume Mahler, Keren Sharon, Masamune Oguri, Maruša Bradač, Larry Bradley, Brenda Frye, Christine J. Forman, Victoria Strait, Yuanyuan Su, Keiichi Umetsu, Felipe Andrade-Santos, Roberto J. Avila, Daniela Carrasco, Catherine Cerny, Nicole G. Czakon, William A. Dawson, Carter Fox, Austin T. Hoag, Kuang-Han Huang, Traci L. Johnson, Shotaro Kikuchihara, Daniel Lam , et al. (16 additional authors not shown)

Abstract: Extensive surveys with the \textit{Hubble Space Telescope} (HST) over the past decade, targeting some of the most massive clusters in the sky, have uncovered dozens of galaxy-cluster strong lenses. The massive cluster strong-lens scale is typically $θ_{E}\sim10\arcsec$ to $\sim30-35\arcsec$, with only a handful of clusters known with Einstein radii $θ_{E}\sim40\arcsec$ or above (for… ▽ More Extensive surveys with the \textit{Hubble Space Telescope} (HST) over the past decade, targeting some of the most massive clusters in the sky, have uncovered dozens of galaxy-cluster strong lenses. The massive cluster strong-lens scale is typically $θ_{E}\sim10\arcsec$ to $\sim30-35\arcsec$, with only a handful of clusters known with Einstein radii $θ_{E}\sim40\arcsec$ or above (for $z_{source}=2$, nominally). Here we report another very large cluster lens, RXC J0032.1+1808 ($z=0.3956$), the second richest cluster in the redMapper cluster catalog and the 85th most massive cluster in the Planck Sunyaev-Zel'dovich catalog. With our Light-Traces-Mass and fully parametric (dPIEeNFW) approaches, we construct strong lensing models based on 18 multiple images of 5 background galaxies newly identified in the \textit{Hubble} data mainly from the \textit{Reionization Lensing Cluster Survey} (RELICS), in addition to a known sextuply imaged system in this cluster. Furthermore, we compare these models to Lenstool and GLAFIC models that were produced independently as part of the RELICS program. All models reveal a large effective Einstein radius of $θ_{E}\simeq40\arcsec$ ($z_{source}=2$), owing to the obvious concentration of substructures near the cluster center. Although RXC J0032.1+1808 has a very large critical area and high lensing strength, only three magnified high-redshift candidates are found within the field targeted by RELICS. Nevertheless, we expect many more high-redshift candidates will be seen in wider and deeper observations with \textit{Hubble} or \emph{JWST}. Finally, the comparison between several algorithms demonstrates that the total error budget is largely dominated by systematic uncertainties. △ Less

Submitted 11 May, 2020; v1 submitted 5 December, 2019; originally announced December 2019.

Comments: 23 pages, accepted for publication in ApJ

Hubble Frontier Field Photometric Catalogues of Abell 370 and RXC J2248.7-4431: Multiwavelength photometry, photometric redshifts, and stellar properties

Authors: Marusa Bradac, Kuang-Han Huang, Marco Castellano, Emiliano Merlin, Ricardo Amorin, Austin Hoag, Victoria Strait, Paola Santini, Russell Ryan, Stefano Casertano, Brian Lemaux, Lori Lubin, Kasper Schmidt, Tim Schrabback-Krahe, Tommaso Treu, Anja von der Linden, Charlotte Mason, Xin Wang

Abstract: This paper presents multiwavelength photometric catalogues of the last two Hubble Frontier Fields (HFF), the massive galaxy clusters Abell 370 and RXC J2248.7-4431. The photometry ranges from imaging performed on the Hubble Space Telescope (HST) to ground based Very Large Telescope (VLT) and Spitzer/IRAC, in collaboration with the ASTRODEEP team, and using the ASTRODEEP pipeline. While the main pu… ▽ More This paper presents multiwavelength photometric catalogues of the last two Hubble Frontier Fields (HFF), the massive galaxy clusters Abell 370 and RXC J2248.7-4431. The photometry ranges from imaging performed on the Hubble Space Telescope (HST) to ground based Very Large Telescope (VLT) and Spitzer/IRAC, in collaboration with the ASTRODEEP team, and using the ASTRODEEP pipeline. While the main purpose of this paper is to release the catalogues, we also perform, as a proof of concept, a brief analysis of z > 6 objects selected using drop-out method, as well as spectroscopically confirmed sources and multiple images in both clusters. While dropout methods yield a sample of high-z galaxies, the addition of longer wavelength data reveals that as expected the samples have substantial contamination at the ~30-45% level by dusty galaxies at lower redshifts. Furthermore, we show that spectroscopic redshifts are still required to unambiguously determine redshifts of multiply imaged systems. Finally, the now publicly available ASTRODEEP catalogues were combined for all HFFs and used to explore stellar properties of a large sample of 20,000 galaxies across a large photometric redshift range. The powerful magnification provided by the HFF clusters allows us an exploration of the properties of galaxies with intrinsic stellar masses as low as $M_* \gtrsim 10^7M_{\odot}$ and intrinsic star formation rates $\mbox{SFRs}\sim 0.1\mbox{-}1M_\odot/\mbox yr$ at z > 6. △ Less

Submitted 31 July, 2019; v1 submitted 4 June, 2019; originally announced June 2019.

Comments: Replaced to match the accepted version

The early growth of supermassive black holes in cosmological hydrodynamic simulations with constrained Gaussian realizations

Authors: Kuan-Wei Huang, Yueying Ni, Yu Feng, Tiziana Di Matteo

Abstract: The paper examines the early growth of supermassive black holes (SMBHs) in cosmological hydrodynamic simulations with different BH seeding scenarios. Employing the constrained Gaussian realization, we reconstruct the initial conditions in the large-volume BlueTides simulation and run them to $z=6$ to cross-validate that the method reproduces the first quasars and their environments. Our constraine… ▽ More The paper examines the early growth of supermassive black holes (SMBHs) in cosmological hydrodynamic simulations with different BH seeding scenarios. Employing the constrained Gaussian realization, we reconstruct the initial conditions in the large-volume BlueTides simulation and run them to $z=6$ to cross-validate that the method reproduces the first quasars and their environments. Our constrained simulations in a volume of $(15\, h^{-1}{\rm Mpc})^3$ successfully recover the evolution of large-scale structure and the stellar and BH masses in the vicinity of a $\sim10^{12}\, M_{\odot}$ halo which we identified in BlueTides at $z\sim7$ hosting a $\sim10^9\, M_{\odot}$ SMBH. Among our constrained simulations, only the ones with a low-tidal field and high-density peak in the initial conditions induce the fastest BH growth required to explain the $z>6$ quasars. We run two sets of simulations with different BH seed masses of $5\times10^3$, $5\times10^4$, and $5\times10^5\, h^{-1}M_{\odot}$, (a) with the same ratio of halo to BH seed mass and (b) with the same halo threshold mass. At $z=6$, all the SMBHs converge in mass to $\sim10^9\, M_{\odot}$ except for the one with the smallest seed in (b) undergoing critical BH growth and reaching $10^8$ -- $10^9\, M_{\odot}$, albeit with most of the growth in (b) delayed compared to set (a). The finding of eight BH mergers in the small-seed scenario (four with masses $10^4$ -- $10^6\, M_{\odot}$ at $z>12$), six in the intermediate-seed scenario, and zero in the large-seed scenario suggests that the vast BHs in the small-seed scenario merge frequently during the early phases of the growth of SMBHs. The increased BH merger rate for the low-mass BH seed and halo threshold scenario provides an exciting prospect for discriminating BH formation mechanisms with the advent of multi-messenger astrophysics and next-generation gravitational wave facilities. △ Less

Submitted 27 May, 2020; v1 submitted 1 June, 2019; originally announced June 2019.

ALMA Polarimetry of AT2018cow

Authors: K. Y. Huang, J. Shimoda, Y. Urata, K. Toma, K. Yamaoka, K. Asada, H. Nagai, S. Takahashi, G. Petitpas, M. Tashiro

Abstract: We present the first radio polarimetric observations of a fast-rising blue optical transient, AT2018cow. Two epochs of polarimetry with additional coincident photometry were performed with the Atacama Large Millimeter/submillimeter Array (ALMA). The overall photometric results based on simultaneous observations in the 100 and 230 GHz bands are consistent with the non-thermal radiation model report… ▽ More We present the first radio polarimetric observations of a fast-rising blue optical transient, AT2018cow. Two epochs of polarimetry with additional coincident photometry were performed with the Atacama Large Millimeter/submillimeter Array (ALMA). The overall photometric results based on simultaneous observations in the 100 and 230 GHz bands are consistent with the non-thermal radiation model reported by Ho et al. (2019) and indicate that the spectral peaks ($\sim110$ GHz at the first epoch and $\sim67$ GHz at the second epoch) represent the synchrotron self-absorption frequency. The non-detection of linear polarization with $<$0.15% in the 230 GHz band at the phase when the effect of synchrotron self-absorption was quite small in the band may be explained by internal Faraday depolarization with high circumburst density and strong magnetic field. This result supports the stellar explosion scenario rather than the tidal disruption model. The maximum energy of accelerating particles at the shocks of AT2018cow-like objects is also discussed. △ Less

Submitted 23 May, 2019; originally announced May 2019.

Comments: 7 pages, 3 figures, 1 table, accepted for publication in ApJL

Stellar Properties of z ~ 8 Galaxies in the Reionization Lensing Cluster Survey

Authors: Victoria Strait, Marusa Bradac, Dan Coe, Larry Bradley, Brett Salmon, Brian C. Lemaux, Kuang-Han Huang, Adi Zitrin, Keren Sharon, Ana Acebron, Felipe Andrade-Santos, Roberto J. Avila, Brenda L. Frye, Austin Hoag, Guillaume Mahler, Mario Nonino, Sara Ogaz, Masamune Oguri, Masami Ouchi, Rachel Paterno-Mahler, Debora Pelliccia

Abstract: Measurements of stellar properties of galaxies when the universe was less than one billion years old yield some of the only observational constraints of the onset of star formation. We present here the inclusion of \textit{Spitzer}/IRAC imaging in the spectral energy distribution fitting of the seven highest-redshift galaxy candidates selected from the \emph{Hubble Space Telescope} imaging of the… ▽ More Measurements of stellar properties of galaxies when the universe was less than one billion years old yield some of the only observational constraints of the onset of star formation. We present here the inclusion of \textit{Spitzer}/IRAC imaging in the spectral energy distribution fitting of the seven highest-redshift galaxy candidates selected from the \emph{Hubble Space Telescope} imaging of the Reionization Lensing Cluster Survey (RELICS). We find that for 6/8 \textit{HST}-selected $z\gtrsim8$ sources, the $z\gtrsim8$ solutions are still strongly preferred over $z\sim$1-2 solutions after the inclusion of \textit{Spitzer} fluxes, and two prefer a $z\sim 7$ solution, which we defer to a later analysis. We find a wide range of intrinsic stellar masses ($5\times10^6 M_{\odot}$ -- $4\times10^9$ $M_{\odot}$), star formation rates (0.2-14 $M_{\odot}\rm yr^{-1}$), and ages (30-600 Myr) among our sample. Of particular interest is Abell1763-1434, which shows evidence of an evolved stellar population at $z\sim8$, implying its first generation of star formation occurred just $< 100$ Myr after the Big Bang. SPT0615-JD, a spatially resolved $z\sim10$ candidate, remains at its high redshift, supported by deep \textit{Spitzer}/IRAC data, and also shows some evidence for an evolved stellar population. Even with the lensed, bright apparent magnitudes of these $z \gtrsim 8$ candidates (H = 26.1-27.8 AB mag), only the \textit{James Webb Space Telescope} will be able further confirm the presence of evolved stellar populations early in the universe. △ Less

Submitted 22 May, 2019; originally announced May 2019.

Comments: 8 pages, 3 figures, 2 tables

First Detection of Radio Linear Polarization in a Gamma Ray Burst Afterglow

Authors: Y. Urata, K. Toma, K. Huang, K. Asada, H. Nagai, S. Takahashi, G. Petitpas, M. Tashiro, K. Yamaoka

Abstract: We report the first detection of radio polarization of a GRB afterglow with the first intensive combined use of telescopes in the millimeter and submillimeter ranges for GRB171205A. The linear polarization degree in the millimeter band at the sub-percent level ($0.27 \pm 0.04\%$) is lower than those observed in late-time optical afterglows (weighted average of $\sim 1\%$). The Faraday depolarizati… ▽ More We report the first detection of radio polarization of a GRB afterglow with the first intensive combined use of telescopes in the millimeter and submillimeter ranges for GRB171205A. The linear polarization degree in the millimeter band at the sub-percent level ($0.27 \pm 0.04\%$) is lower than those observed in late-time optical afterglows (weighted average of $\sim 1\%$). The Faraday depolarization by non-accelerated, cool electrons in the shocked region is one of possible mechanisms for the low value. In this scenario, larger total energy by a factor of $\sim 10$ than ordinary estimate without considering non-accelerated electrons is required. The polarization position angle varies by at least 20 degrees across the millimeter band, which is not inconsistent with this scenario. This result indicates that polarimetry in the millimeter and submillimeter ranges is a unique tool for investigating GRB energetics, and coincident observations with multiple frequencies or bands would provide more accurate measurements of the non-accelerated electron fraction. △ Less

Submitted 29 September, 2019; v1 submitted 17 April, 2019; originally announced April 2019.

Comments: Accepted to ApJ Letters; 8 pages, 4 figures, 1 table,

RELICS: Reionization Lensing Cluster Survey

Authors: Dan Coe, Brett Salmon, Marusa Bradac, Larry D. Bradley, Keren Sharon, Adi Zitrin, Ana Acebron, Catherine Cerny, Nathalia Cibirka, Victoria Strait, Rachel Paterno-Mahler, Guillaume Mahler, Roberto J. Avila, Sara Ogaz, Kuang-Han Huang, Debora Pelliccia, Daniel P. Stark, Ramesh Mainali, Pascal A. Oesch, Michele Trenti, Daniela Carrasco, William A. Dawson, Steven A. Rodney, Louis-Gregory Strolger, Adam G. Riess , et al. (32 additional authors not shown)

Abstract: Large surveys of galaxy clusters with the Hubble and Spitzer Space Telescopes, including CLASH and the Frontier Fields, have demonstrated the power of strong gravitational lensing to efficiently deliver large samples of high-redshift galaxies. We extend this strategy through a wider, shallower survey named RELICS, the Reionization Lensing Cluster Survey. This survey, described here, was designed p… ▽ More Large surveys of galaxy clusters with the Hubble and Spitzer Space Telescopes, including CLASH and the Frontier Fields, have demonstrated the power of strong gravitational lensing to efficiently deliver large samples of high-redshift galaxies. We extend this strategy through a wider, shallower survey named RELICS, the Reionization Lensing Cluster Survey. This survey, described here, was designed primarily to deliver the best and brightest high-redshift candidates from the first billion years after the Big Bang. RELICS observed 41 massive galaxy clusters with Hubble and Spitzer at 0.4-1.7um and 3.0-5.0um, respectively. We selected 21 clusters based on Planck PSZ2 mass estimates and the other 20 based on observed or inferred lensing strength. Our 188-orbit Hubble Treasury Program obtained the first high-resolution near-infrared images of these clusters to efficiently search for lensed high-redshift galaxies. We observed 46 WFC3/IR pointings (~200 arcmin^2) with two orbits divided among four filters (F105W, F125W, F140W, and F160W) and ACS imaging as needed to achieve single-orbit depth in each of three filters (F435W, F606W, and F814W). As previously reported by Salmon et al., we discovered 322 z ~ 6 - 10 candidates, including the brightest known at z ~ 6, and the most distant spatially-resolved lensed arc known at z ~ 10. Spitzer IRAC imaging (945 hours awarded, plus 100 archival) has crucially enabled us to distinguish z ~ 10 candidates from z ~ 2 interlopers. For each cluster, two HST observing epochs were staggered by about a month, enabling us to discover 11 supernovae, including 3 lensed supernovae, which we followed up with 20 orbits from our program. We delivered reduced HST images and catalogs of all clusters to the public via MAST and reduced Spitzer images via IRSA. We have also begun delivering lens models of all clusters, to be completed before the JWST GO call for proposals. △ Less

Submitted 5 March, 2019; originally announced March 2019.

Comments: 29 pages, 6 figures, submitted to ApJ. For reduced images, catalogs, lens models, and more, see relics.stsci.edu

3D Magnetic Reconnection with a spatially confined X-line extent -- Implications for Dipolarizing Flux Bundles and the Dawn-Dusk Asymmetry

Authors: Yi-Hsin Liu, Tak Chu Li, Michael Hesse, Weijie Sun, Jiang Liu, James Burch, James A. Slavin, Kai Huang

Abstract: Using 3D particle-in-cell (PIC) simulations, we study magnetic reconnection with the x-line being spatially confined in the current direction. We include thick current layers to prevent reconnection at two ends of a thin current sheet that has a thickness on an ion inertial (di) scale. The reconnection rate and outflow speed drop significantly when the extent of the thin current sheet in the curre… ▽ More Using 3D particle-in-cell (PIC) simulations, we study magnetic reconnection with the x-line being spatially confined in the current direction. We include thick current layers to prevent reconnection at two ends of a thin current sheet that has a thickness on an ion inertial (di) scale. The reconnection rate and outflow speed drop significantly when the extent of the thin current sheet in the current direction is < O(10 di). When the thin current sheet extent is long enough, we find it consists of two distinct regions; an inactive region (on the ion-drifting side) exists adjacent to the active region where reconnection proceeds normally as in a 2D case. The extent of this inactive region is ~ O(10 di), and it suppresses reconnection when the thin current sheet extent is comparable or shorter. The time-scale of current sheet thinning toward fast reconnection can be translated into the spatial-scale of this inactive region; because electron drifts inside the ion diffusion region transport the reconnected magnetic flux, that drives outflows and furthers the current sheet thinning, away from this region. This is a consequence of the Hall effect in 3D. While this inactive region may explain the shortest possible azimuthal extent of dipolarizing flux bundles at Earth, it may also explain the dawn-dusk asymmetry observed at the magnetotail of Mercury, that has a global dawn-dusk extent much shorter than that of Earth. △ Less

Submitted 29 January, 2019; originally announced January 2019.

Comments: 9 pages, 9 figures, submitted to JGR on 01/23/2019

Constraining the neutral fraction of hydrogen in the IGM at redshift 7.5

Authors: Austin Hoag, Maruša Bradač, Kuang-Han Huang, Charlotte Mason, Tommaso Treu, Kasper B. Schmidt, Michele Trenti, Victoria Strait, Brian C. Lemaux, Emily Finney, Maia Paddock

Abstract: We present a large spectroscopic campaign with Keck/MOSFIRE targeting Lyman-alpha emission (Ly$α$) from intrinsically faint Lyman-break Galaxies (LBGs) behind 12 efficient galaxy cluster lenses. Gravitational lensing allows us to probe the more abundant faint galaxy population to sensitive Ly$α$ equivalent width limits. During the campaign we targeted 70 LBG candidates with MOSFIRE Y-band, selecte… ▽ More We present a large spectroscopic campaign with Keck/MOSFIRE targeting Lyman-alpha emission (Ly$α$) from intrinsically faint Lyman-break Galaxies (LBGs) behind 12 efficient galaxy cluster lenses. Gravitational lensing allows us to probe the more abundant faint galaxy population to sensitive Ly$α$ equivalent width limits. During the campaign we targeted 70 LBG candidates with MOSFIRE Y-band, selected photometrically to cover Ly$α$ over the range $7<z<8.2$. We detect $S/N>5$ emission lines in 2 of these galaxies and find that they are likely Ly$α$ at $z=7.148\pm0.001$ and $z=7.161\pm0.001$. We present new lens models for 4 of the galaxy clusters, using our previously published lens models for the remaining clusters to determine the magnification factors for the source galaxies. Using a Bayesian framework that employs large scale reionization simulations of the intergalactic medium (IGM) as well as realistic properties of the interstellar medium and circumgalactic medium, we infer the volume-averaged neutral hydrogen fraction, $\overline{x}_{\mathrm{HI}}$, in the IGM during reionization to be $\overline{x}_{\mathrm{HI}}=0.88^{+0.05}_{-0.10}$ at $z=7.6\pm0.6$. Our result is consistent with a late and rapid reionization scenario inferred by Planck. △ Less

Submitted 25 January, 2019; originally announced January 2019.

RELICS: Strong Lensing Analysis of MACS J0417.5-1154 and Predictions for Observing the Magnified High-Redshift Universe with JWST

Authors: Guillaume Mahler, Keren Sharon, Carter Fox, Dan Coe, Mathilde Jauzac, Victoria Strait, Alastair Edge, Ana Acebron, Felipe Andrade-Santos, Roberto J. Avila, Maruša Bradač, Larry D. Bradley, Daniela Carrasco, Catherine Cerny, Nathália Cibirka, Nicole G. Czakon, William A. Dawson, Brenda L. Frye, Austin T. Hoag, Kuang-Han Huang, Traci L. Johnson, Christine Jones, Shotaro Kikuchihara, 15 Daniel Lam, Rachael Livermore , et al. (15 additional authors not shown)

Abstract: Strong gravitational lensing by clusters of galaxies probes the mass distribution at the core of each cluster and magnifies the universe behind it. MACS J0417.5-1154 at z=0.443 is one of the most massive clusters known based on weak lensing, X-ray, and Sunyaev-Zel'dovich analyses. Here we compute a strong lens model of MACS J0417 based on Hubble Space Telescope imaging observations collected, in p… ▽ More Strong gravitational lensing by clusters of galaxies probes the mass distribution at the core of each cluster and magnifies the universe behind it. MACS J0417.5-1154 at z=0.443 is one of the most massive clusters known based on weak lensing, X-ray, and Sunyaev-Zel'dovich analyses. Here we compute a strong lens model of MACS J0417 based on Hubble Space Telescope imaging observations collected, in part, by the Reionization Lensing Cluster Survey (RELICS), and recently reported spectroscopic redshifts from the MUSE instrument on the Very Large Telescope (VLT). We measure an Einstein radius of ThetaE=36'' at z = 9 and a mass projected within 200 kpc of M(200 kpc) = 1.78+0.01-0.03x10**14Msol. Using this model, we measure a ratio between the mass attributed to cluster-member galaxy halos and the main cluster halo of order 1:100. We assess the probability to detect magnified high-redshift galaxies in the field of this cluster, both for comparison with RELICS HST results and as a prediction for the James Webb Space Telescope (JWST) Guaranteed Time Observations upcoming for this cluster. Our lensing analysis indicates that this cluster has similar lensing strength to other clusters in the RELICS program. Our lensing analysis predicts a detection of at least a few z~6-8 galaxies behind this cluster, at odds with a recent analysis that yielded no such candidates in this field. Reliable strong lensing models are crucial for accurately predicting the intrinsic properties of lensed galaxies. As part of the RELICS program, our strong lensing model produced with the Lenstool parametric method is publicly available through the Mikulski Archive for Space Telescopes (MAST). △ Less

Submitted 10 March, 2019; v1 submitted 31 October, 2018; originally announced October 2018.

Comments: 20 pages, 8 figures. Accepted to ApJ

RELICS: High-Resolution Constraints on the Inner Mass Distribution of the z=0.83 Merging Cluster RXJ0152.7-1357 from strong lensing

Authors: Ana Acebron, May Alon, Adi Zitrin, Guillaume Mahler, Dan Coe, Keren Sharon, Nathália Cibirka, Maruša Bradač, Michele Trenti, Keiichi Umetsu, Felipe Andrade-Santos, Roberto J. Avila, Larry Bradley, Daniela Carrasco, Catherine Cerny, Nicole G. Czakon, William A. Dawson, Brenda Frye, Austin T. Hoag, Kuang-Han Huang, Traci L. Johnson, Christine Jones, Shotaro Kikuchihara, Daniel Lam, Rachael C. Livermore , et al. (15 additional authors not shown)

Abstract: Strong gravitational lensing (SL) is a powerful means to map the distribution of dark matter. In this work, we perform a SL analysis of the prominent X-ray cluster RXJ0152.7-1357 (z=0.83, also known as CL 0152.7-1357) in \textit{Hubble Space Telescope} images, taken in the framework of the Reionization Lensing Cluster Survey (RELICS). On top of a previously known $z=3.93$ galaxy multiply imaged by… ▽ More Strong gravitational lensing (SL) is a powerful means to map the distribution of dark matter. In this work, we perform a SL analysis of the prominent X-ray cluster RXJ0152.7-1357 (z=0.83, also known as CL 0152.7-1357) in \textit{Hubble Space Telescope} images, taken in the framework of the Reionization Lensing Cluster Survey (RELICS). On top of a previously known $z=3.93$ galaxy multiply imaged by RXJ0152.7-1357, for which we identify an additional multiple image, guided by a light-traces-mass approach we identify seven new sets of multiply imaged background sources lensed by this cluster, spanning the redshift range [1.79-3.93]. A total of 25 multiple images are seen over a small area of ~0.4 $arcmin^2$, allowing us to put relatively high-resolution constraints on the inner matter distribution. Although modestly massive, the high degree of substructure together with its very elongated shape make RXJ0152.7-1357 a very efficient lens for its size. This cluster also comprises the third-largest sample of z~6-7 candidates in the RELICS survey. Finally, we present a comparison of our resulting mass distribution and magnification estimates with those from a Lenstool model. These models are made publicly available through the MAST archive. △ Less

Submitted 26 February, 2019; v1 submitted 18 October, 2018; originally announced October 2018.

Comments: 15 Pages, 7 Figures, 4 Tables Accepted for publication in ApJ

Mass Modeling of Frontier Fields Cluster MACS J1149.5+2223 Using Strong and Weak Lensing

Authors: Emily Quinn Finney, Maruša Bradač, Kuang-Han Huang, Austin Hoag, Takahiro Morishita, Tim Schrabback, Tommaso Treu, Kasper Borello Schmidt, Brian C. Lemaux, Xin Wang, Charlotte Mason

Abstract: We present a gravitational lensing model of MACS J1149.5+2223 using ultra-deep Hubble Frontier Fields imaging data and spectroscopic redshifts from HST grism and VLT/MUSE spectroscopic data. We create total mass maps using 38 multiple images (13 sources) and 608 weak lensing galaxies, as well as 100 multiple images of 31 star-forming regions in the galaxy that hosts Supernova Refsdal. We find good… ▽ More We present a gravitational lensing model of MACS J1149.5+2223 using ultra-deep Hubble Frontier Fields imaging data and spectroscopic redshifts from HST grism and VLT/MUSE spectroscopic data. We create total mass maps using 38 multiple images (13 sources) and 608 weak lensing galaxies, as well as 100 multiple images of 31 star-forming regions in the galaxy that hosts Supernova Refsdal. We find good agreement with a range of recent models within the HST field of view. We present a map of the ratio of projected stellar mass to total mass ($f_{\star}$), and find that the stellar mass fraction for this cluster peaks on the primary BCG. Averaging within a radius of 0.3 Mpc, we obtain a value of $\langle f_{\star} \rangle = 0.012^{+0.004}_{-0.003}$, consistent with other recent results for this ratio in cluster environments, though with a large global error (up to $δf_{\star} = 0.005$) primarily due to the choice of an IMF. We compare values of $f_{\star}$ and measures of star formation efficiency for this cluster to other Hubble Frontier Fields clusters studied in the literature, finding that MACS1149 has a higher stellar mass fraction than these other clusters, but a star formation efficiency typical of massive clusters. △ Less

Submitted 2 June, 2018; originally announced June 2018.

Comments: 15 pages, 6 figures, published in ApJ May 23, 2018

Journal ref: The Astrophysical Journal, 859, 1, 2018

RELICS: A Strong Lens Model for SPT-CLJ0615-5746, a z=0.972 Cluster

Authors: Rachel Paterno-Mahler, Keren Sharon, Dan Coe, Guillaume Mahler, Catherine Cerny, Traci Johnson, Tim Schrabback, Felipe Andrade-Santos, Roberto J. Avila, Marusa Bradac, Larry D. Bradley, Daniela Carrasco, Nicole G. Czakon, William A. Dawson, Brenda L. Frye, Austin T. Hoag, Kuang-Han Huang, Christine Jones, Daniel Lam, Rachael Livermore, Lorenzo Lovisari, Ramesh Mainali, Pascal A. Oesch, Sara Ogaz, Matthew Past , et al. (8 additional authors not shown)

Abstract: We present a lens model for the cluster SPT-CLJ0615$-$5746, which is the highest redshift ($z=0.972$) system in the Reionization of Lensing Clusters Survey (RELICS), making it the highest redshift cluster for which a full strong lens model is published. We identify three systems of multiply-imaged lensed galaxies, two of which we spectroscopically confirm at $z=1.358$ and $z=4.013$, which we use a… ▽ More We present a lens model for the cluster SPT-CLJ0615$-$5746, which is the highest redshift ($z=0.972$) system in the Reionization of Lensing Clusters Survey (RELICS), making it the highest redshift cluster for which a full strong lens model is published. We identify three systems of multiply-imaged lensed galaxies, two of which we spectroscopically confirm at $z=1.358$ and $z=4.013$, which we use as constraints for the model. We find a foreground structure at $z\sim0.4$, which we include as a second cluster-sized halo in one of our models; however two different statistical tests find the best-fit model consists of one cluster-sized halo combined with three individually optimized galaxy-sized halos, as well as contributions from the cluster galaxies themselves. We find the total projected mass density within $r=26.7"$ (the region where the strong lensing constraints exist) to be $M=2.51^{+0.15}_{-0.09}\times 10^{14}$~M$_{\odot}$. If we extrapolate out to $r_{500}$, our projected mass density is consistent with the mass inferred from weak lensing and from the Sunyaev-Zel'dovich effect ($M\sim10^{15}$~M$_{\odot}$). This cluster is lensing a previously reported $z\sim10$ galaxy, which, if spectroscopically confirmed, will be the highest-redshift strongly lensed galaxy known. △ Less

Submitted 9 July, 2018; v1 submitted 24 May, 2018; originally announced May 2018.

Comments: 15 pages, 8 figures 4 tables. ApJ Accepted

Mass and Light of Abell 370: A Strong and Weak Lensing Analysis

Authors: Victoria Strait, Marusa Bradac, Austin Hoag, Kuang-Han Huang, Tommaso Treu, Xin Wang, Ricardo Amorin, Marco Castellano, Adriano Fontana, Brian Lemaux, Emiliano Merlin, Kasper B. Schmidt, Tim Schrabback, Adam Tomczack, Michele Trenti, Benedetta Vulcani

Abstract: We present a new gravitational lens model of the Hubble Frontier Fields cluster Abell 370 ($z = 0.375$) using imaging and spectroscopy from Hubble Space Telescope and ground-based spectroscopy. We combine constraints from a catalog of 1344 weakly lensed galaxies and 39 multiply-imaged sources comprised of 114 multiple images, including a system of multiply-imaged candidates at $z=7.93 \pm 0.02$, t… ▽ More We present a new gravitational lens model of the Hubble Frontier Fields cluster Abell 370 ($z = 0.375$) using imaging and spectroscopy from Hubble Space Telescope and ground-based spectroscopy. We combine constraints from a catalog of 1344 weakly lensed galaxies and 39 multiply-imaged sources comprised of 114 multiple images, including a system of multiply-imaged candidates at $z=7.93 \pm 0.02$, to obtain a best-fit mass distribution using the cluster lens modeling code Strong and Weak Lensing United. As the only analysis of A370 using strong and weak lensing constraints from Hubble Frontier Fields data, our method provides an independent check on assumptions in other methods on the mass distribution. Convergence, shear, and magnification maps are made publicly available through the HFF website. We find that the model we produce is similar to models produced by other groups, with some exceptions due to the differences in lensing code methodology. In an effort to study how our total projected mass distribution traces light, we measure the stellar mass density distribution using Spitzer/Infrared Array Camera imaging. Comparing our total mass density to our stellar mass density in a radius of 0.3 Mpc, we find a mean projected stellar to total mass ratio of $\langle f* \rangle = 0.011 \pm 0.003$ (stat.) using the diet Salpeter initial mass function. This value is in general agreement with independent measurements of $\langle f* \rangle$ in clusters of similar total mass and redshift. △ Less

Submitted 11 October, 2018; v1 submitted 22 May, 2018; originally announced May 2018.