Showing 1–50 of 253 results for author: Tanaka, I

The RAdio Galaxy Environment Reference Survey (RAGERS): Evidence of an anisotropic distribution of submillimeter galaxies in the 4C 23.56 protocluster at z=2.48

Authors: Dazhi Zhou, Thomas R. Greve, Bitten Gullberg, Minju M. Lee, Luca Di Mascolo, Simon R. Dicker, Charles E. Romero, Scott C. Chapman, Chian-Chou Chen, Thomas Cornish, Mark J. Devlin, Luis C. Ho, Kotaro Kohno, Claudia D. P. Lagos, Brian S. Mason, Tony Mroczkowski, Jeff F. W. Wagg, Q. Daniel Wang, Ran Wang, Malte. Brinch, Helmut Dannerbauer, Xue-Jian Jiang, Lynge R. B. Lauritsen, Aswin P. Vijayan, David Vizgan , et al. (19 additional authors not shown)

Abstract: High-redshift radio(-loud) galaxies (H$z$RGs) are massive galaxies with powerful radio-loud active galactic nuclei (AGNs) and serve as beacons for protocluster identification. However, the interplay between H$z$RGs and the large-scale environment remains unclear. To understand the connection between H$z$RGs and the surrounding obscured star formation, we investigated the overdensity and spatial di… ▽ More High-redshift radio(-loud) galaxies (H$z$RGs) are massive galaxies with powerful radio-loud active galactic nuclei (AGNs) and serve as beacons for protocluster identification. However, the interplay between H$z$RGs and the large-scale environment remains unclear. To understand the connection between H$z$RGs and the surrounding obscured star formation, we investigated the overdensity and spatial distribution of submillimeter-bright galaxies (SMGs) in the field of 4C\,23.56, a well-known H$z$RG at $z=2.48$. We used SCUBA-2 data ($σ\,{\sim}\,0.6$\,mJy) to estimate the $850\,{\rm μm}$ source number counts and examine the radial and azimuthal overdensities of the $850\,{\rm μm}$ sources in the vicinity of the H$z$RG. The angular distribution of SMGs is inhomogeneous around the H$z$RG 4C\,23.56, with fewer sources oriented along the radio jet. We also find a significant overdensity of bright SMGs (${\rm S}_{850\rm\,μm}\geq5\,$mJy). Faint and bright SMGs exhibit different spatial distributions. The former are concentrated in the core region, while the latter prefer the outskirts of the H$z$RG field. High-resolution observations show that the seven brightest SMGs in our sample are intrinsically bright, suggesting that the overdensity of bright SMGs is less likely due to the source multiplicity. △ Less

Submitted 4 August, 2024; originally announced August 2024.

Comments: 19 pages, 17 figures, 5 tables, accepted to A&A

Star formation activity of low-mass galaxies at the peak epoch of galaxy formation probed by deep narrow-band imaging

Authors: Kazuki Daikuhara, Tadayuki Kodama, Jose M. Pérez-Martínez, Rhythm Shimakawa, Tomoko L. Suzuki, Ken-ichi Tadaki, Yusei Koyama, Ichi Tanaka

Abstract: Low-mass galaxies at high redshifts are the building blocks of more massive galaxies at later times and are thus key populations for understanding galaxy formation and evolution. We have made deep narrow-band observations for two protoclusters and the general field in COSMOS at $z$ $\sim$ 2. In a clumpy young protocluster, USS1558$-$003, at $z$ = 2.53, we find many star-forming galaxies well above… ▽ More Low-mass galaxies at high redshifts are the building blocks of more massive galaxies at later times and are thus key populations for understanding galaxy formation and evolution. We have made deep narrow-band observations for two protoclusters and the general field in COSMOS at $z$ $\sim$ 2. In a clumpy young protocluster, USS1558$-$003, at $z$ = 2.53, we find many star-forming galaxies well above the star-forming main sequence of field galaxies at the low-mass end ($M_{\star}/\mathrm{M_{\odot}}<10^{8.9}$). This suggests that some environmental effects may be at work in low-mass galaxies in high-density regions to enhance their star formation activities. In the core of this protocluster, we also find that enhanced star formation activity of middle-mass galaxies ($10^{8.9} < M_{\star}/\mathrm{M_{\odot}} < 10^{10.2}$) while such trends are not observed in a more mature protocluster, PKS1138$-$262 at $z$ = 2.16. We expect these activities to be mainly due to galaxy mergers/interactions and differences in the amount of cold gas accretion. As one piece of evidence, we show that the star formation activity within individual galaxies in the protoclusters is more centrally concentrated than those in the field. This is probably due to the enhanced interactions between galaxies in the protocluster, which can reduce the angular momentum of the gas, drive the gas towards the galaxy center, and lead to a central starburst. △ Less

Submitted 31 May, 2024; originally announced May 2024.

Superconducting phase diagram of Sb and Se substitution in CeOBiS2 single crystals

Authors: Tatsuya Suzuki, Masanori Nagao, Yuki Maruyama, Satoshi Watauchi, Isao Tanaka

Abstract: Sb- and Se-substituted CeOBiS2 single crystals have been successfully grown using CsCl/KCl flux. Sb and Se substitution dependence of the superconductivity on CeOBiS2 was investigated and the superconducting phase diagram at above 0.3 K was described. The non-linear boundary between superconductivity and non-superconductivity was revealed. Sb- and Se-substituted CeOBiS2 single crystals have been successfully grown using CsCl/KCl flux. Sb and Se substitution dependence of the superconductivity on CeOBiS2 was investigated and the superconducting phase diagram at above 0.3 K was described. The non-linear boundary between superconductivity and non-superconductivity was revealed. △ Less

Submitted 23 April, 2024; originally announced April 2024.

Comments: 16 pages, 2 figures

Journal ref: Journal of the Physical Society of Japan 93, 054708 (2024)

INvestigations of massive Filaments ANd sTar formation (INFANT). I. Core Identification and Core Mass Function

Authors: Yu Cheng, Xing Lu, Patricio Sanhueza, Hauyu Baobab Liu, Qizhou Zhang, Roberto Galván-Madrid, Ke Wang, Fumitaka Nakamura, Tie Liu, Siyi Feng, Shanghuo Li, Sihan Jiao, Kei E. I. Tanaka, Xunchuan Liu, Pak Shing Li, Qiuyi Luo, Qilao Gu, Yuxin Lin, András E. Guzmán

Abstract: Filamentary structures are ubiquitously found in high-mass star-forming clouds. To investigate the relationship between filaments and star formation, we carry out the INFANT (INvestigations of massive Filaments ANd sTar formation) survey, a multi-scale, multi-wavelength survey of massive filamentary clouds with ALMA band 3/band 6 and VLA K band. In this first paper, we present the ALMA band 6 cont… ▽ More Filamentary structures are ubiquitously found in high-mass star-forming clouds. To investigate the relationship between filaments and star formation, we carry out the INFANT (INvestigations of massive Filaments ANd sTar formation) survey, a multi-scale, multi-wavelength survey of massive filamentary clouds with ALMA band 3/band 6 and VLA K band. In this first paper, we present the ALMA band 6 continuum observations toward a sample of 8 high-mass star forming filaments. We covered each target with approximately rectangular mosaic field of view with two 12-m array configurations, achieving an angular resolution of $\sim$0.6" (2700 AU at 4.5 kpc) and a continuum rms of $\sim$0.1 mJy/beam ($\sim$0.06 Msun in gas mass assuming 15 K). We identify cores using the getsf and astrodendro and find the former is more robust in terms of both identification and measuring flux densities. We identify in total 183 dense cores (15--36 cores in each cloud) and classify their star formation states via outflow and warm gas tracers. The protostellar cores are statistically more massive than the prestellar cores, possibly indicating further accretion onto cores after formation of protostars. For the high-mass end ($M_\text{core}$ $>$ 1.5 Msun) of the core mass function (CMF) we derive a power-law index of $-$1.15 $\pm$ 0.12 for the whole sample, and $-$1.70 $\pm$ 0.25 for the prestellar population. We also find a steepening trend in CMF with cloud evolution ($-$0.89 $\pm$ 0.15 for the young group v.s. $-$1.44 $\pm$ 0.25 for the evolved group) and discuss its implication for cluster formation. △ Less

Submitted 4 April, 2024; originally announced April 2024.

Comments: 25 pages, 8 figures, accepted for ApJ

Imaging quantum interference in a monolayer Kitaev quantum spin liquid candidate

Authors: Y. Kohsaka, S. Akutagawa, S. Omachi, Y. Iwamichi, T. Ono, I. Tanaka, S. Tateishi, H. Murayama, S. Suetsugu, K. Hashimoto, T. Shibauchi, M. O. Takahashi, M. G. Yamada, S. Nikolaev, T. Mizushima, S. Fujimoto, T. Terashima, T. Asaba, Y. Kasahara, Y. Matsuda

Abstract: Single atomic defects are prominent windows to look into host quantum states because collective responses from the host states emerge as localized states around the defects. Friedel oscillations and Kondo clouds in Fermi liquids are quintessential examples. However, the situation is quite different for quantum spin liquid (QSL), an exotic state of matter with fractionalized quasiparticles and topo… ▽ More Single atomic defects are prominent windows to look into host quantum states because collective responses from the host states emerge as localized states around the defects. Friedel oscillations and Kondo clouds in Fermi liquids are quintessential examples. However, the situation is quite different for quantum spin liquid (QSL), an exotic state of matter with fractionalized quasiparticles and topological order arising from a profound impact of quantum entanglement. Elucidating the underlying local electronic property has been challenging due to the charge neutrality of fractionalized quasiparticles and the insulating nature of QSLs. Here, using spectroscopic-imaging scanning tunneling microscopy, we report atomically resolved images of monolayer $α$-RuCl$_3$, the most promising Kitaev QSL candidate, on metallic substrates. We find quantum interference in the insulator manifesting as incommensurate and decaying spatial oscillations of the local density of states around defects with a characteristic bias dependence. The oscillation differs from any known spatial structures in its nature and does not exist in other Mott insulators, implying it is an exotic oscillation involved with excitations unique to $α$-RuCl$_3$. Numerical simulations can reproduce the observed oscillation by assuming that itinerant Majorana fermions of Kitaev QSL are scattered across the Majorana Fermi surface. The oscillation provides a new approach to exploring Kitaev QSLs through the local response against defects like Friedel oscillations in metals. △ Less

Submitted 25 March, 2024; originally announced March 2024.

JWST's first glimpse of a z > 2 forming cluster reveals a top-heavy stellar mass function

Authors: Hanwen Sun, Tao Wang, Ke Xu, Emanuele Daddi, Qing Gu, Tadayuki Kodama, Anita Zanella, David Elbaz, Ichi Tanaka, Raphael Gobat, Qi Guo, Jiaxin Han, Shiying Lu, Luwenjia Zhou

Abstract: Clusters and their progenitors (protoclusters) at z = 2-4, the peak epoch of star formation, are ideal laboratories to study the formation process of both the clusters themselves and their member galaxies. However, a complete census of their member galaxies has been challenging due to observational difficulties. Here we present new JWST/NIRCam observations targeting the distant cluster CLJ1001 at… ▽ More Clusters and their progenitors (protoclusters) at z = 2-4, the peak epoch of star formation, are ideal laboratories to study the formation process of both the clusters themselves and their member galaxies. However, a complete census of their member galaxies has been challenging due to observational difficulties. Here we present new JWST/NIRCam observations targeting the distant cluster CLJ1001 at z = 2.51 from the COSMOS-Web program, which, in combination with previous narrowband imaging targeting H-alpha emitters and deep millimeter surveys of CO emitters, provide a complete view of massive galaxy assembly in CLJ1001. In particular, JWST reveals a population of massive, extremely red cluster members in the long-wavelength bands that were invisible in previous Hubble Space Telescope (HST)/F160W imaging (HST-dark members). Based on this highly complete spectroscopic sample of member galaxies, we show that the spatial distribution of galaxies in CLJ1001 exhibits a strong central concentration, with the central galaxy density already resembling that of low-z clusters. Moreover, we reveal a "top-heavy" stellar mass function for the star-forming galaxies (SFGs), with an overabundance of massive SFGs piled up in the cluster core. These features strongly suggest that CLJ1001 is caught in a rapid transition, with many of its massive SFGs likely soon becoming quiescent. In the context of cluster formation, these findings suggest that the earliest clusters form from the inside out and top to bottom, with the massive galaxies in the core assembling first, followed by the less massive ones in the outskirts. △ Less

Submitted 29 May, 2024; v1 submitted 8 March, 2024; originally announced March 2024.

Comments: 14 pages, 9 figures, 1 table, published by ApJL

Journal ref: ApJL, 967, L34 (2024)

Predictive power of polynomial machine learning potentials for liquid states in 22 elemental systems

Authors: Hayato Wakai, Atsuto Seko, Hirosato Izuta, Takayuki Nishiyama, Isao Tanaka

Abstract: The polynomial machine learning potentials (MLPs) described by polynomial rotational invariants have been systematically developed for various systems and used in diverse applications in crystalline states. In this study, we systematically investigate the predictive power of the polynomial MLPs for liquid structural properties in 22 elemental systems with diverse chemical bonding properties, inclu… ▽ More The polynomial machine learning potentials (MLPs) described by polynomial rotational invariants have been systematically developed for various systems and used in diverse applications in crystalline states. In this study, we systematically investigate the predictive power of the polynomial MLPs for liquid structural properties in 22 elemental systems with diverse chemical bonding properties, including those showing anomalous melting behavior, such as Si, Ge, and Bi. We compare liquid structural properties obtained from molecular dynamics simulations using the density functional theory (DFT) calculation, the polynomial MLPs, and other interatomic potentials in the literature. The current results demonstrate that the polynomial MLPs consistently exhibit high predictive power for liquid structural properties with the same accuracy as that of typical DFT calculations. △ Less

Submitted 11 June, 2024; v1 submitted 26 January, 2024; originally announced January 2024.

Comments: REVTeX 4-2, 30 pages with 35 figures

Growth and Characterization of Superconducting Bulk Crystal [(SnSe)$_{1+δ}$]$_m$(NbSe$_2$) Misfit Layer Compounds

Authors: Ryufa Shu, Masanori Nagao, Chiaya Yamamoto, Keisuke Arimoto, Junji Yamanaka, Yuki Maruyama, Satoshi Watauchi, Isao Tanaka

Abstract: [(SnSe)$_{1+δ}$]$_m$(NbSe$_2$) ($m$ = 1-6, 8, and 12) highly orientated crystals 1-2 mm in size and well-defined c-planes were successfully grown using CsCl/KCl flux, including the first growth of crystals with $m = 12$. The stacked layers along the $c$ axis in the obtained crystals were directly observed by transmission electron microscopy as m alternating layers of SnSe and single layers of NbSe… ▽ More [(SnSe)$_{1+δ}$]$_m$(NbSe$_2$) ($m$ = 1-6, 8, and 12) highly orientated crystals 1-2 mm in size and well-defined c-planes were successfully grown using CsCl/KCl flux, including the first growth of crystals with $m = 12$. The stacked layers along the $c$ axis in the obtained crystals were directly observed by transmission electron microscopy as m alternating layers of SnSe and single layers of NbSe$_2$. The superconducting transition temperature of the obtained [(SnSe)$_{1+δ}$]$_m$(NbSe$_2$) crystals decreased with an increase in the number of SnSe layers per unit cell. As the superconducting anisotropy parameters increase, a significant increase is observed between $m = 4$ and 5. This indicates that the superconducting dimensionality becomes more two-dimensional with an increasing $m$. △ Less

Submitted 14 January, 2024; originally announced January 2024.

Report number: Journal of Alloys and Compounds, vol.978 (2024) 173486

Journal ref: Journal of Alloys and Compounds, vol.978 (2024) 173486

Suppression of superconductivity by Sb-substitution into CeOBiS2 single crystals

Authors: Tatsuya Suzuki, Yuji Hanada, Masanori Nagao, Yuki Maruyama, Satoshi Watauchi, Isao Tanaka

Abstract: Sb-substituted CeOBiS2 single crystals with 0.2-1.0 mm size have been successfully grown using CsCl/KCl flux. Sb substitution strongly suppressed the superconductivity in CeOBiS2, and the substitution of more than approximately 4 atomic percent in Bi-site disappeared the superconductivity at a measurement range of above 0.36 K. Furthermore, the drastic change of c-axis lattice parameter was observ… ▽ More Sb-substituted CeOBiS2 single crystals with 0.2-1.0 mm size have been successfully grown using CsCl/KCl flux. Sb substitution strongly suppressed the superconductivity in CeOBiS2, and the substitution of more than approximately 4 atomic percent in Bi-site disappeared the superconductivity at a measurement range of above 0.36 K. Furthermore, the drastic change of c-axis lattice parameter was observed in a neighborhood to Sb substitution amount, which maybe suggests the structural transition in CeOBiS2. The Ce valence state of the obtained single crystals was evaluated. However, the established correlations between Ce valence state and superconductivity could not be revealed. △ Less

Submitted 18 December, 2023; originally announced December 2023.

Journal ref: Solid State Communications, vol.378 (2024) 115414

Ownership Types for Verification of Programs with Pointer Arithmetic

Authors: Izumi Tanaka, Ken Sakayori, Naoki Kobayashi

Abstract: Toman et al. have proposed a type system for automatic verification of low-level programs, which combines ownership types and refinement types to enable strong updates of refinement types in the presence of pointer aliases. We extend their type system to support pointer arithmetic, and prove its soundness. Based on the proposed type system, we have implemented a prototype tool for automated verifi… ▽ More Toman et al. have proposed a type system for automatic verification of low-level programs, which combines ownership types and refinement types to enable strong updates of refinement types in the presence of pointer aliases. We extend their type system to support pointer arithmetic, and prove its soundness. Based on the proposed type system, we have implemented a prototype tool for automated verification of the lack of assertion errors of low-level programs with pointer arithmetic, and confirmed its effectiveness through experiments. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: An extended version of the paper to appear in Proceedings of PEPM 2024

Enhanced star formation and metallicity deficit in the USS 1558-003 forming protocluster at z=2.53

Authors: Jose Manuel Pérez-Martínez, Tadayuki Kodama, Yusei Koyama, Rhythm Shimakawa, Tomoko L. Suzuki, Kazuki Daikuhara, Kota Adachi, Masato Onodera, Ichi Tanaka

Abstract: We use K-band multi-object near-infrared spectroscopy with Keck/MOSFIRE to search for environmental imprints on the gas properties of 27 narrow-band selected H$α$ emitters (HAEs) across the three major clumps of the assembling USS1558--003 protocluster at $z=2.53$. We target the H$α$ and [NII]$λ$6584 emission lines to obtain star-formation rates (SFR) and gas-phase oxygen abundances for our source… ▽ More We use K-band multi-object near-infrared spectroscopy with Keck/MOSFIRE to search for environmental imprints on the gas properties of 27 narrow-band selected H$α$ emitters (HAEs) across the three major clumps of the assembling USS1558--003 protocluster at $z=2.53$. We target the H$α$ and [NII]$λ$6584 emission lines to obtain star-formation rates (SFR) and gas-phase oxygen abundances for our sources, confirming the membership of 23 objects. HAEs belonging to this protocluster display enhanced SFRs with respect to the main sequence of star formation at the same cosmic epoch. This effect is more prominent for low-mass galaxies ($\mathrm{\log M_*/M_\odot<10.0}$), which may be experiencing a vigorous phase of mass assembly shortly after they were formed. We compute the individual and stacked gas-phase metallicities for our sources finding a metallicity deficit for low-mass objects when compared against the field mass-metallicity relation and the massive Spiderweb protocluster at $z=2.16$. These results suggest that HAEs within USS1558--003 may be less evolved than those in the Spiderweb protocluster. Finally, we explore the gas metallicity - gas fraction relation for a small sample of five galaxies with CO(3-2) molecular gas information. Assuming our objects are in equilibrium, we obtain a relatively wide range of mass loading factors ($\mathrm{λ=0.5-2}$) matching field samples at the cosmic noon but in contrast with our previous results in the Spiderweb protocluster. We speculate that these discrepancies between protoclusters may be (partly) driven by differences in their current dynamical and mass assembly stages, hinting at the co-evolution of protoclusters and their galaxy populations at $2<z<3$. △ Less

Submitted 6 December, 2023; originally announced December 2023.

Comments: 18 pages, 10 figures. Accepted for publication in MNRAS

The Detection of Higher-Order Millimeter Hydrogen Recombination Lines in the Large Magellanic Cloud

Authors: Marta Sewiło, Kazuki Tokuda, Stan E. Kurtz, Steven B. Charnley, Thomas Möller, Jennifer Wiseman, C. -H. Rosie Chen, Remy Indebetouw, Álvaro Sánchez-Monge, Kei E. I. Tanaka, Peter Schilke, Toshikazu Onishi, Naoto Harada

Abstract: We report the first extragalactic detection of the higher-order millimeter hydrogen recombination lines ($Δn>2$). The $γ$-, $ε$-, and $η$-transitions have been detected toward the millimeter continuum source N105-1A in the star-forming region N105 in the Large Magellanic Cloud (LMC) with the Atacama Large Millimeter/submillimeter Array (ALMA). We use the H40$α$ line, the brightest of the detected… ▽ More We report the first extragalactic detection of the higher-order millimeter hydrogen recombination lines ($Δn>2$). The $γ$-, $ε$-, and $η$-transitions have been detected toward the millimeter continuum source N105-1A in the star-forming region N105 in the Large Magellanic Cloud (LMC) with the Atacama Large Millimeter/submillimeter Array (ALMA). We use the H40$α$ line, the brightest of the detected recombination lines (H40$α$, H36$β$, H50$β$, H41$γ$, H57$γ$, H49$ε$, H53$η$, and H54$η$), and/or the 3 mm free-free continuum emission to determine the physical parameters of N105-1A (the electron temperature, emission measure, electron density, and size) and study ionized gas kinematics. We compare the physical properties of N105-1A to a large sample of Galactic compact and ultracompact (UC) H II regions and conclude that N105-1A is similar to the most luminous ($L>10^5$ $L_{\odot}$) UC H II regions in the Galaxy. N105-1A is ionized by an O5.5 V star, it is deeply embedded in its natal molecular clump, and likely associated with a (proto)cluster. We incorporate high-resolution molecular line data including CS, SO, SO$_2$, and CH$_3$OH ($\sim$0.12 pc), and HCO$^{+}$ and CO ($\sim$0.087 pc) to explore the molecular environment of N105-1A. Based on the CO data, we find evidence for a cloud-cloud collision that likely triggered star formation in the region. We find no clear outflow signatures, but the presence of filaments and streamers indicates on-going accretion onto the clump hosting the UC H II region. Sulfur chemistry in N105-1A is consistent with the accretion shock model predictions. △ Less

Submitted 5 September, 2023; originally announced September 2023.

Comments: 51 pages, 30 figures, 2 tables (including appendices); accepted for publication in The Astrophysical Journal (ApJ)

Chemical evolution of some selected complex organic molecules in low-mass star-forming regions

Authors: Bratati Bhat, Rumela Kar, Suman Kumar Mondal, Rana Ghosh, Prasanta Gorai, Takashi Shimonishi, Kei E. I. Tanaka, Kenji Furuya, Ankan Das

Abstract: The destiny of complex organic molecules (COMs) in star-forming regions is interlinked with various evolutionary phases. Therefore, identifying these species in diversified environments of identical star-forming regions would help to comprehend their physical and chemical heritage. We identified multiple COMs utilizing the Large Program `Astrochemical Surveys At IRAM' (ASAI) data, dedicated to che… ▽ More The destiny of complex organic molecules (COMs) in star-forming regions is interlinked with various evolutionary phases. Therefore, identifying these species in diversified environments of identical star-forming regions would help to comprehend their physical and chemical heritage. We identified multiple COMs utilizing the Large Program `Astrochemical Surveys At IRAM' (ASAI) data, dedicated to chemical surveys in Sun-like star-forming regions with the IRAM 30 m telescope. It was an unbiased survey in the millimetre regime, covering the prestellar core, protostar, outflow region, and protoplanetary disk phase. Here, we have reported some transitions of seven COMs, namely, methanol (CH3OH), acetaldehyde (CH3CHO), methyl formate (CH3OCHO), ethanol (C2H5OH), propynal (HCCCHO), dimethyl ether (CH3OCH3), and methyl cyanide (CH3CN) in some sources L1544, B1-b, IRAS4A, and SVS13A. We found a trend among these species from the derived abundances using the rotational diagram method and MCMC fit. We have found that the abundances of all of the COMs, except for HCCCHO, increase from the L1544 (prestellar core) and peaks at IRAS16293-2422 (class 0 phase). It is noticed that the abundance of these molecules correlate with the luminosity of the sources. The obtained trend is also visible from the previous interferometric observations and considering the beam dilution effect. △ Less

Submitted 20 August, 2023; originally announced August 2023.

Comments: 44 pages, 25 figures, and 12 tables. Accepted for the publication in the Astrophysical Journal

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

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

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

Submitted 10 August, 2023; originally announced August 2023.

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

Algorithm for spin symmetry operation search

Authors: Kohei Shinohara, Atsushi Togo, Hikaru Watanabe, Takuya Nomoto, Isao Tanaka, Ryotaro Arita

Abstract: A spin space group provides a suitable way to fully exploit the symmetry of a spin arrangement with a negligible spin-orbit coupling. There has been a growing interest in applying spin symmetry analysis with the spin space group in the field of magnetism. However, there is no established algorithm to search for spin symmetry operations of the spin space group. This paper presents an exhaustive alg… ▽ More A spin space group provides a suitable way to fully exploit the symmetry of a spin arrangement with a negligible spin-orbit coupling. There has been a growing interest in applying spin symmetry analysis with the spin space group in the field of magnetism. However, there is no established algorithm to search for spin symmetry operations of the spin space group. This paper presents an exhaustive algorithm for determining spin symmetry operations of commensurate spin arrangements. The present algorithm searches for spin symmetry operations from the symmetry operations of a corresponding nonmagnetic crystal structure and determines their spin-rotation parts by solving a Procrustes problem. An implementation is distributed under a permissive free software license in spinspg v0.1.1: https://github.com/spglib/spinspg. △ Less

Submitted 21 October, 2023; v1 submitted 23 July, 2023; originally announced July 2023.

New insights into the role of AGNs in forming the cluster red sequence

Authors: Rhythm Shimakawa, José Manuel Pérez-Martínez, Yusei Koyama, Masayuki Tanaka, Ichi Tanaka, Tadayuki Kodama, Nina A. Hatch, Huub J. A. Röttgering, Helmut Dannerbauer, Jaron D. Kurk

Abstract: As a considerable investment of time from various telescope facilities were dedicated toward studying the Spiderweb protocluster at $z=2.2$, it so far remains one of the most extensively studied protocluster. We report here the latest results in this field, adding a new dimension to previous research on cluster formation at high redshift. Previous studies have reported a significant overdensity (… ▽ More As a considerable investment of time from various telescope facilities were dedicated toward studying the Spiderweb protocluster at $z=2.2$, it so far remains one of the most extensively studied protocluster. We report here the latest results in this field, adding a new dimension to previous research on cluster formation at high redshift. Previous studies have reported a significant overdensity ($δ\sim10$) of massive H$α$ (+ [Nii]) -emitting galaxies in 3700 comoving Mpc$^3$. Many of these were previously considered to be dusty, actively star-forming galaxies, given their rest-frame optical and infrared features. However, this study argues that a third of them are more likely to be "passively-evolving" galaxies with low-luminosity active galactic nuclei (AGNs) rather than star-forming galaxies, given the multi-wavelength spectral energy distribution (SED) fitting including an AGN component. For their SED-based star formation rates to be valid, bulk of their H$α$ + [Nii] emission should come from the central AGNs. This difference in interpretation between this work and past studies, including ours, is particularly supported by the recent deep Chandra X-ray observation. Furthermore, we have spectroscopically confirmed a quiescent nature for one of these AGNs, with its multiple stellar absorption lines but also low ionisation emission lines. This important update provides new insights into the role of AGNs in forming the cluster red sequence observed in the present-day universe. △ Less

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

Comments: 17 pages, 10 figures, 5 tables, accepted for publication in MNRAS

Observational signatures of forming young massive clusters: continuum emission from dense HII regions

Authors: Mutsuko Inoguchi, Takashi Hosokawa, Hajime Fukushima, Kei E. I. Tanaka, Hidenobu Yajima, Shin Mineshige

Abstract: Young massive clusters (YMCs) are the most massive star clusters forming in nearby galaxies and are thought to be a young analogue to the globular clusters. Understanding the formation process of YMCs leads to looking into very efficient star formation in high-redshift galaxies suggested by recent JWST observations. We investigate possible observational signatures of their formation stage, particu… ▽ More Young massive clusters (YMCs) are the most massive star clusters forming in nearby galaxies and are thought to be a young analogue to the globular clusters. Understanding the formation process of YMCs leads to looking into very efficient star formation in high-redshift galaxies suggested by recent JWST observations. We investigate possible observational signatures of their formation stage, particularly when the mass of a cluster is increasing via accretion from a natal molecular cloud. To this end, we study the broad-band continuum emission from ionized gas and dust enshrouding YMCs, whose formation is followed by recent radiation-hydrodynamics simulations. We perform post-process radiative transfer calculations using simulation snapshots and find characteristic spectral features at radio and far-infrared frequencies. We show that a striking feature is long-lasting, strong free-free emission from a $\sim$ 10pc-scale HII region with a large emission measure of $\gtrsim 10^7 \mathrm{cm}^{-6} \ \mathrm{pc}$, corresponding to the mean electron density of $\gtrsim 10^3~\mathrm{cm}^{-3}$. There is a turnover feature below $\sim$ 10 GHz, a signature of the optically-thick free-free emission, often found in Galactic ultra-compact HII regions. These features come from the peculiar YMC formation process, where the cluster's gravity effectively traps photoionized gas for a long duration and enables continuous star formation within the cluster. Such large and dense HII regions show distinct distribution on the density-size diagram, apart from the standard sequence of Galactic HII regions. This is consistent with the observational trend inferred for extragalactic HII regions associated with YMCs. △ Less

Submitted 25 October, 2023; v1 submitted 30 May, 2023; originally announced May 2023.

Comments: 12 pages, 10 figures, accepted for publication in MNRAS

Growth of self-integrated atomic quantum wires and junctions of a Mott semiconductor

Authors: Tomoya Asaba, Lang Peng, Takahiro Ono, Satoru Akutagawa, Ibuki Tanaka, Hinako Murayama, Shota Suetsugu, Aleksandar Razpopov, Yuichi Kasahara, Takahito Terashima, Yuhki Kohsaka, Takasada Shibauchi, Masatoshi Ichikawa, Roser Valentí, Shin-ichi Sasa, Yuji Matsuda

Abstract: Continued advances in quantum technologies rely on producing nanometer-scale wires. Although several state-of-the-art nanolithographic technologies and bottom-up synthesis processes have been used to engineer such wires, critical challenges remain in growing uniform atomic-scale crystalline wires and constructing their network structures. Here we discover a simple method to fabricate atomic-scale… ▽ More Continued advances in quantum technologies rely on producing nanometer-scale wires. Although several state-of-the-art nanolithographic technologies and bottom-up synthesis processes have been used to engineer such wires, critical challenges remain in growing uniform atomic-scale crystalline wires and constructing their network structures. Here we discover a simple method to fabricate atomic-scale wires with various arrangements, including stripes, X-, Y-junctions, and nanorings. Single-crystalline atomic-scale wires of a Mott insulator, whose band gap is comparable to those of wide-gap semiconductors, are spontaneously grown on graphite substrates \DEL{and epitaxial monolayer graphene on SiC }by pulsed-laser deposition. These wires are one-unit-cell-thick and have an exact width of two- and four-unit-cells (1.4 and 2.8\,nm) and lengths up to a few $μm$. We show that the non-equilibrium reaction-diffusion processes may play an essential role in atomic pattern formation. Our findings offer a new perspective on the non-equilibrium self-organization phenomena on an atomic scale, paving a unique way for the quantum architecture of nano-network. △ Less

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

Journal ref: Sci. Adv.9,eabq5561(2023)

Effects of equivalent composition on superconducting properties of high-entropy REOBiS$_2$ (RE = La, Ce, Pr, Nd, Sm, Gd) single crystals

Authors: Yuma Fujita, Masanori Nagao, Akira Miura, Daisuke Urushihara, Yoshikazu Mizuguchi, Yuki Maruyama, Satoshi Watauchi, Yoshihiko Takano, Isao Tanaka

Abstract: Superconductors are influenced by high-entropy alloys (HEAs); these have been investigated in various functional materials. REOBiS$_2$ (RE = La, Ce, Pr, Nd, Sm, and Gd in different combinations) single crystals with HEAs at the RE-site were successfully grown using the flux method. The obtained crystals were plate-shaped (1 mm$^2$) with a well-developed c-plane. Ce was present in both trivalent (C… ▽ More Superconductors are influenced by high-entropy alloys (HEAs); these have been investigated in various functional materials. REOBiS$_2$ (RE = La, Ce, Pr, Nd, Sm, and Gd in different combinations) single crystals with HEAs at the RE-site were successfully grown using the flux method. The obtained crystals were plate-shaped (1 mm$^2$) with a well-developed c-plane. Ce was present in both trivalent (Ce$^{3+}$) and tetravalent (Ce$^{4+}$) electronic configurations; the concentration of Ce$^{4+}$ at the RE-site was approximately 10 at% in all single crystals. The single crystals showed superconducting transition temperature with zero resistivity within 1.2-4.2 K. The superconducting transition temperature, superconducting anisotropy, electronic specific heat coefficient, and Debye temperature of the crystals were not correlated with the mixed entropy at the RE-site. Except for the electronic specific heat coefficient, the variation of these parameters as a function of mixed entropy showed different trends for equivalent and non-equivalent RE element compositions. Thus, the configuration of RE elements influences the superconducting properties of REOBiS$_2$ single crystals, alluding to a method of modulating transition temperatures. △ Less

Submitted 11 April, 2023; originally announced April 2023.

Journal ref: Physica C, vol.608 (2023) 1354254

Growth of REBa2Cu3Ox single-crystal whiskers utilizing the concept of high-entropy alloys

Authors: Yasunori Suzuki, Masanori Nagao, Yuki Maruyama, Satoshi Watauchi, Isao Tanaka

Abstract: Single-crystal whiskers of REBa2Cu3Ox (RE-123, RE: rare earth element) cuprate superconductors were successfully grown using Te-doped precursors utilizing the concept of high-entropy-alloys (HEAs) at the RE site. The obtained whiskers were 0.5-1.0 mm length with approximately 10 um thickness. The width of the flat surface, which corresponds to the ab-plane, was 10-50 um. The RE site of the grown w… ▽ More Single-crystal whiskers of REBa2Cu3Ox (RE-123, RE: rare earth element) cuprate superconductors were successfully grown using Te-doped precursors utilizing the concept of high-entropy-alloys (HEAs) at the RE site. The obtained whiskers were 0.5-1.0 mm length with approximately 10 um thickness. The width of the flat surface, which corresponds to the ab-plane, was 10-50 um. The RE site of the grown whiskers could be easily substituted by large ionic radius RE elements (Gd, Dy, and Ho). Substitution with small ionic radius RE elements (Tm, Yb, and Lu) resulted in the opposite trend. The superconducting transition temperature and superconducting anisotropy of the grown RE-123 whiskers were 89-93 K and 6-10, respectively. The mixed entropy at the RE site in the RE-123 whiskers did not affect the superconducting transition temperature and superconducting anisotropy. △ Less

Submitted 22 March, 2023; originally announced March 2023.

Comments: 27 pages, 8 figures, 2 tables

Journal ref: Jpn. J. Appl. Phys. 62 (2023) 033001

Massive Protostellar Disks as a Hot Laboratory of Silicate Grain Evolution

Authors: Ryota Yamamuro, Kei E. I. Tanaka, Satoshi Okuzumi

Abstract: Typical accretion disks around massive protostars are hot enough for water ice to sublimate. We here propose to utilize the massive protostellar disks for investigating the collisional evolution of silicate grains with no ice mantle, which is an essential process for the formation of rocky planetesimals in protoplanetary disks around lower-mass stars. We for the first time develop a model of massi… ▽ More Typical accretion disks around massive protostars are hot enough for water ice to sublimate. We here propose to utilize the massive protostellar disks for investigating the collisional evolution of silicate grains with no ice mantle, which is an essential process for the formation of rocky planetesimals in protoplanetary disks around lower-mass stars. We for the first time develop a model of massive protostellar disks that includes the coagulation, fragmentation, and radial drift of dust. We show that the maximum grain size in the disks is limited by collisional fragmentation rather than by radial drift. We derive analytic formulas that produce the radial distribution of the maximum grain size and dust surface density in the steady state. Applying the analytic formulas to the massive protostellar disk of GGD27-MM1, where the grain size is constrained from a millimeter polarimetric observation, we infer that the silicate grains in this disk fragment at collision velocities above ~ 10 m/s. The inferred fragmentation threshold velocity is lower than the maximum grain collision velocity in typical protoplanetary disks around low-mass stars, implying that coagulation alone may not lead to the formation of rocky planetesimals in those disks. With future measurements of grain sizes in massive protostellar disks, our model will provide more robust constraints on the sticking property of silicate grains. △ Less

Submitted 6 June, 2023; v1 submitted 16 March, 2023; originally announced March 2023.

Comments: 17 pages, 5 figures,accepted for publication to The Astrophysical Journal

Report number: AAS43716R2 MSC Class: 85

The Detection of Hot Molecular Cores in the Small Magellanic Cloud

Authors: Takashi Shimonishi, Kei E. I. Tanaka, Yichen Zhang, Kenji Furuya

Abstract: We report the first detection of hot molecular cores in the Small Magellanic Cloud, a nearby dwarf galaxy with 0.2 solar metallicity. We observed two high-mass young stellar objects in the SMC with ALMA, and detected emission lines of CO, HCO+, H13CO+, SiO, H2CO, CH3OH, SO, and SO2. Compact hot-core regions are traced by SO2, whose spatial extent is about 0.1 pc, and the gas temperature is higher… ▽ More We report the first detection of hot molecular cores in the Small Magellanic Cloud, a nearby dwarf galaxy with 0.2 solar metallicity. We observed two high-mass young stellar objects in the SMC with ALMA, and detected emission lines of CO, HCO+, H13CO+, SiO, H2CO, CH3OH, SO, and SO2. Compact hot-core regions are traced by SO2, whose spatial extent is about 0.1 pc, and the gas temperature is higher than 100 K based on the rotation diagram analysis. In contrast, CH3OH, a classical hot-core tracer, is dominated by extended (0.2-0.3 pc) components in both sources, and the gas temperature is estimated to be 39+-8 K for one source. Protostellar outflows are also detected from both sources as high-velocity components of CO. The metallicity-scaled abundances of SO2 in hot cores are comparable among the SMC, LMC, and Galactic sources, suggesting that the chemical reactions leading to SO2 formation would be regulated by elemental abundances. On the other hand, CH3OH shows a large abundance variation within SMC and LMC hot cores. The diversity in the initial condition of star formation (e.g., degree of shielding, local radiation field strength) may lead to the large abundance variation of organic molecules in hot cores. This work, in conjunction with previous hot-core studies in the LMC and outer/inner Galaxy, suggests that the formation of a hot core would be a common phenomenon during high-mass star formation across the metallicity range of 0.2-1 solar metallicity. High-excitation SO2 lines will be a useful hot-core tracer in the low-metallicity environments of the SMC and LMC. △ Less

Submitted 9 March, 2023; originally announced March 2023.

Comments: Accepted for publication in ApJL, 17 pages, 8 figures, 4 tables. arXiv admin note: text overlap with arXiv:2109.11123

Astrochemical Diagnostics of the Isolated Massive Protostar G28.20-0.05

Authors: Prasanta Gorai, Chi-Yan Law, Jonathan C. Tan, Yichen Zhang, Ruben Fedriani, Kei E. I. Tanaka, Melisse Bonfand, Giuliana Cosentino, Diego Mardones, Maria T. Beltran, Guido Garay

Abstract: We study the astrochemical diagnostics of the isolated massive protostar G28.20-0.05. We analyze data from ALMA 1.3~mm observations with resolution of 0.2 arcsec ($\sim$1,000 au). We detect emission from a wealth of species, including oxygen-bearing (e.g., $\rm{H_2CO}$, $\rm{CH_3OH}$, $\rm{CH_3OCH_3}$), sulfur-bearing (SO$_2$, H$_2$S) and nitrogen-bearing (e.g., HNCO, NH$_2$CHO, C$_2$H$_3$CN, C… ▽ More We study the astrochemical diagnostics of the isolated massive protostar G28.20-0.05. We analyze data from ALMA 1.3~mm observations with resolution of 0.2 arcsec ($\sim$1,000 au). We detect emission from a wealth of species, including oxygen-bearing (e.g., $\rm{H_2CO}$, $\rm{CH_3OH}$, $\rm{CH_3OCH_3}$), sulfur-bearing (SO$_2$, H$_2$S) and nitrogen-bearing (e.g., HNCO, NH$_2$CHO, C$_2$H$_3$CN, C$_2$H$_5$CN) molecules. We discuss their spatial distributions, physical conditions, correlation between different species and possible chemical origins. In the central region near the protostar, we identify three hot molecular cores (HMCs). HMC1 is part of a mm continuum ring-like structure, is closest in projection to the protostar, has the highest temperature of $\sim300\:$K, and shows the most line-rich spectra. HMC2 is on the other side of the ring, has a temperature of $\sim250\:$K, and is of intermediate chemical complexity. HMC3 is further away, $\sim3,000\:$au in projection, cooler ($\sim70\:$K) and is the least line-rich. The three HMCs have similar mass surface densities ($\sim10\:{\rm{g\:cm}}^{-2}$), number densities ($n_{\rm H}\sim10^9\:{\rm{cm}}^{-3}$) and masses of a few $M_\odot$. The total gas mass in the cores and in the region out to $3,000\:$au is $\sim 25\:M_\odot$, which is comparable to that of the central protostar. Based on spatial distributions of peak line intensities as a function of excitation energy, we infer that the HMCs are externally heated by the protostar. We estimate column densities and abundances of the detected species and discuss the implications for hot core astrochemistry. △ Less

Submitted 2 November, 2023; v1 submitted 5 March, 2023; originally announced March 2023.

Comments: Accepted for publication in ApJ

Follow-up survey for the binary black hole merger GW200224_222234 using Subaru/HSC and GTC/OSIRIS

Authors: Takayuki Ohgami, Josefa Becerra Gonzalez, Nozomu Tominaga, Tomoki Morokuma, Yousuke Utsumi, Yuu Niino, Masaomi Tanaka, Smaranika Banerjee, Frederick Poidevin, Jose Antonio Acosta-Pulido, Ismael Perez-Fournon, Teo Munoz-Darias, Hiroshi Akitaya, Kenshi Yanagisawa, Mahito Sasada, Michitoshi Yoshida, Mirko Simunovic, Ryou Ohsawa, Ichi Tanaka, Tsuyoshi Terai, Yuhei Takagi, The J-GEM collaboration

Abstract: The LIGO/Virgo detected a gravitational wave (GW) event, named GW200224_222234 (a.k.a. S200224ca) and classified as a binary-black-hole coalescence, on February 24, 2020. Given its relatively small localization skymap (71 deg$^2$ for a 90% credible region; revised to 50 deg$^2$ in GWTC-3), we performed target-of-opportunity observations using the Subaru/Hyper Suprime-Cam (HSC) in the $r2$- and… ▽ More The LIGO/Virgo detected a gravitational wave (GW) event, named GW200224_222234 (a.k.a. S200224ca) and classified as a binary-black-hole coalescence, on February 24, 2020. Given its relatively small localization skymap (71 deg$^2$ for a 90% credible region; revised to 50 deg$^2$ in GWTC-3), we performed target-of-opportunity observations using the Subaru/Hyper Suprime-Cam (HSC) in the $r2$- and $z$-bands. Observations were conducted on February 25 and 28 and March 23, 2020, with the first epoch beginning 12.3 h after the GW detection. The survey covered the highest probability sky area of 56.6 deg$^2$, corresponding to a 91% probability. This was the first deep follow-up ($m_{r}\gtrsim24, m_{z}\gtrsim23$) for a binary-black-hole merger covering $>$90% of the localization. By performing image subtraction and candidate screening including light curve fitting with transient templates and examples, we found 22 off-nucleus transients that were not ruled out as the counterparts of GW200224_222234 with only our Subaru/HSC data. We also performed GTC/OSIRIS spectroscopy of the probable host galaxies for five candidates; two are likely to be located within the 3D skymap, whereas the others are not. In conclusion, 19 transients remain as possible optical counterparts of GW200224_222234; however, we could not identify a unique promising counterpart. If there are no counterparts in the remaining candidates, the upper limits of optical luminosity are $νL_ν < 5.2^{+2.4}_{-1.9}\times 10^{41}$ erg s$^{-1}$ and $νL_ν < 1.8^{+0.8}_{-0.6}\times 10^{42}$ erg s$^{-1}$ in the $r2$- and $z$-bands, respectively, at $\sim$12 h after GW detection. We also discuss improvements in the strategies of optical follow-ups for future GW events. △ Less

Submitted 18 February, 2023; originally announced February 2023.

Comments: 22 pages, 15 figures. Accepted for publication in the Astrophysical Journal

Implementation strategies in phonopy and phono3py

Authors: Atsushi Togo, Laurent Chaput, Terumasa Tadano, Isao Tanaka

Abstract: Scientific simulation codes are public property sustained by the community. Modern technology allows anyone to join scientific software projects, from anywhere, remotely via the internet. The phonopy and phono3py codes are widely used open source phonon calculation codes. This review describes a collection of computational methods and techniques as implemented in these codes and shows their impl… ▽ More Scientific simulation codes are public property sustained by the community. Modern technology allows anyone to join scientific software projects, from anywhere, remotely via the internet. The phonopy and phono3py codes are widely used open source phonon calculation codes. This review describes a collection of computational methods and techniques as implemented in these codes and shows their implementation strategies as a whole, aiming to be useful for the community. Some of the techniques presented here are not limited to phonon calculations and may therefore be useful in other area of condensed matter physics. △ Less

Submitted 13 January, 2023; originally announced January 2023.

Disk Wind Feedback from High-mass Protostars. II. The Evolutionary Sequence

Authors: Jan E. Staff, Kei E. I. Tanaka, Jon P. Ramsey, Yichen Zhang, Jonathan C. Tan

Abstract: Star formation is ubiquitously associated with the ejection of accretion-powered outflows that carve bipolar cavities through the infalling envelope. This feedback is expected to be important for regulating the efficiency of star formation from a natal pre-stellar core. These low-extinction outflow cavities greatly affect the appearance of a protostar by allowing the escape of shorter wavelength p… ▽ More Star formation is ubiquitously associated with the ejection of accretion-powered outflows that carve bipolar cavities through the infalling envelope. This feedback is expected to be important for regulating the efficiency of star formation from a natal pre-stellar core. These low-extinction outflow cavities greatly affect the appearance of a protostar by allowing the escape of shorter wavelength photons. Doppler-shifted CO line emission from outflows is also often the most prominent manifestation of deeply embedded early-stage star formation. Here, we present 3D magneto-hydrodynamic simulations of a disk wind outflow from a protostar forming from an initially $60\:M_\odot$ core embedded in a high pressure environment typical of massive star-forming regions. We simulate the growth of the protostar from $m_*=1\:M_\odot$ to $26\:M_\odot$ over a period of $\sim$100,000 years. The outflow quickly excavates a cavity with half opening angle of $\sim10^\circ$ through the core. This angle remains relatively constant until the star reaches $4\:M_\odot$. It then grows steadily in time, reaching a value of $\sim 50^\circ$ by the end of the simulation. We estimate a lower limit to the star formation efficiency (SFE) of 0.43. However, accounting for continued accretion from a massive disk and residual infall envelope, we estimate that the final SFE may be as high as $\sim0.7$. We examine observable properties of the outflow, especially the evolution of the cavity opening angle, total mass and momentum flux, and velocity distributions of the outflowing gas, and compare with the massive protostars G35.20-0.74N and G339.88-1.26 observed by ALMA, yielding constraints on their intrinsic properties. △ Less

Submitted 20 February, 2023; v1 submitted 2 January, 2023; originally announced January 2023.

Comments: Accepted for publication in ApJ

Algorithms for magnetic symmetry operation search and identification of magnetic space group from magnetic crystal structure

Authors: Kohei Shinohara, Atsushi Togo, Isao Tanaka

Abstract: A crystal symmetry search is crucial for computational crystallography and materials science. Although algorithms and implementations for the crystal symmetry search have been developed, their extension to magnetic space groups (MSGs) remains limited. In this paper, algorithms for determining magnetic symmetry operations of magnetic crystal structures, identifying magnetic space-group types of giv… ▽ More A crystal symmetry search is crucial for computational crystallography and materials science. Although algorithms and implementations for the crystal symmetry search have been developed, their extension to magnetic space groups (MSGs) remains limited. In this paper, algorithms for determining magnetic symmetry operations of magnetic crystal structures, identifying magnetic space-group types of given MSGs, searching for transformations to a BNS setting, and symmetrizing the magnetic crystal structures using the MSGs are presented. The determination of magnetic symmetry operations is numerically stable and is implemented with minimal modifications from the existing crystal symmetry search. Magnetic space-group types and transformations to the BNS setting are identified by a two-step approach combining space-group type identification and the use of affine normalizers. Point coordinates and magnetic moments of the magnetic crystal structures are symmetrized by projection operators for the MSGs. An implementation is distributed with a permissive free software license in spglib v2.0.2: https://github.com/spglib/spglib. △ Less

Submitted 7 June, 2023; v1 submitted 27 November, 2022; originally announced November 2022.

Salt-bearing disk candidates around high-mass young stellar objects

Authors: Adam Ginsburg, Brett A. McGuire, Patricio Sanhueza, Fernando Olguin, Luke T Maud, Kei E. I. Tanaka, Yichen Zhang, Henrik Beuther, Nick Indriolo

Abstract: Molecular lines tracing the orbital motion of gas in a well-defined disk are valuable tools for inferring both the properties of the disk and the star it surrounds. Lines that arise only from a disk, and not also from the surrounding molecular cloud core that birthed the star or from the outflow it drives, are rare. Several such emission lines have recently been discovered in one example case, tho… ▽ More Molecular lines tracing the orbital motion of gas in a well-defined disk are valuable tools for inferring both the properties of the disk and the star it surrounds. Lines that arise only from a disk, and not also from the surrounding molecular cloud core that birthed the star or from the outflow it drives, are rare. Several such emission lines have recently been discovered in one example case, those from NaCl and KCl salt molecules. We studied a sample of 23 candidate high-mass young stellar objects (HMYSOs) in 17 high-mass star-forming regions to determine how frequently emission from these species is detected. We present five new detections of water, NaCl, KCl, PN, and SiS from the innermost regions around the objects, bringing the total number of known briny disk candidates to nine. Their kinematic structure is generally disk-like, though we are unable to determine whether they arise from a disk or outflow in the sources with new detections. We demonstrate that these species are spatially coincident in a few resolved cases and show that they are generally detected together, suggesting a common origin or excitation mechanism. We also show that several disks around HMYSOs clearly do not exhibit emission in these species. Salty disks are therefore neither particularly rare in high-mass disks, nor are they ubiquitous. △ Less

Submitted 4 November, 2022; originally announced November 2022.

Comments: accepted to ApJ

A Monomer Mott Insulator (BEDT-TTF)Cu[N(CN)$_2$]$_2$ as a Potential Nodal Line System

Authors: N. Yoneyama, M. K. Nuryadin, T. Tsumuraya, S. Iguchi, T. Takei, N. Kumada, M. Nagao, I. Tanaka, T. Sasaki

Abstract: We report the band structure calculations and the experimental results of resistivity and magnetic susceptibility in a spin-1/2 (BEDT-TTF)$^{\bullet +}$ monomer Mott insulator (BEDT-TTF)Cu[N(CN)$_2$]$_2$. The band calculations indicate a Dirac semimetal state with nodal lines at the Fermi level. The resistivity and the magnetic susceptibility as functions of temperature are well interpreted in ter… ▽ More We report the band structure calculations and the experimental results of resistivity and magnetic susceptibility in a spin-1/2 (BEDT-TTF)$^{\bullet +}$ monomer Mott insulator (BEDT-TTF)Cu[N(CN)$_2$]$_2$. The band calculations indicate a Dirac semimetal state with nodal lines at the Fermi level. The resistivity and the magnetic susceptibility as functions of temperature are well interpreted in terms of the monomer Mott insulating state instead of the expected semimetal state probably owing to strong electron correlation. In addition, we find that an Arrhenius-type steep reduction of the paramagnetic susceptibility appears below approximately 25 K, which indicates a spin-singlet ground state. △ Less

Submitted 28 September, 2022; v1 submitted 28 September, 2022; originally announced September 2022.

Comments: 4 pages, 5 figures. J. Phys. Soc. Jpn (2022) in press

Journal ref: J. Phys. Soc. Jpn. vol. 91 113704 (2022)

Signs of environmental effects on star-forming galaxies in the Spiderweb protocluster at z=2.16

Authors: Jose Manuel Pérez-Martínez, Helmut Dannerbauer, Tadayuki Kodama, Yusei Koyama, Rhythm Shimakawa, Tomoko L. Suzuki, Rosa Calvi, Zhengyi Chen, Kazuki Daikuhara, Nina A. Hatch, Andrés Laza-Ramos, David Sobral, John P. Stott, Ichi Tanaka

Abstract: We use multi-object near-infrared (NIR) spectroscopy with VLT/KMOS to investigate the role of the environment in the evolution of the ionized gas properties of narrow-band selected H$α$ emitters (HAEs) in the Spiderweb protocluster at $z=2.16$. Based on rest-frame optical emission lines, H$α$ and [NII]$λ$6584, we confirm the cluster membership of 39 of our targets (i.e. 93% success rate), and meas… ▽ More We use multi-object near-infrared (NIR) spectroscopy with VLT/KMOS to investigate the role of the environment in the evolution of the ionized gas properties of narrow-band selected H$α$ emitters (HAEs) in the Spiderweb protocluster at $z=2.16$. Based on rest-frame optical emission lines, H$α$ and [NII]$λ$6584, we confirm the cluster membership of 39 of our targets (i.e. 93% success rate), and measure their star-formation rates (SFR), gas-phase oxygen abundances and effective radius. We parametrize the environment where our targets reside by using local and global density indicators based on previous samples of spectroscopic and narrow-band cluster members. We find that star-forming galaxies embedded in the Spiderweb protocluster display SFRs compatible with those of the main sequence and morphologies comparable to those of late-type galaxies at $z=2.2$ in the field. We also report a mild gas-phase metallicity enhancement ($0.6\pm0.3$ dex) at intermediate stellar masses. Furthermore, we identify two UVJ-selected quiescent galaxies with residual H$α$-based star formation and find signs of extreme dust obscuration in a small sample of SMGs based on their FIR and H$α$ emission. Interestingly, the spatial distribution of these objects differs from the rest of HAEs, avoiding the protocluster core. Finally, we explore the gas fraction-gas metallicity diagram for 7 galaxies with molecular gas masses measured by ATCA using CO(1-0). In the context of the gas-regulator model, our objects are consistent with relatively low mass-loading factors, suggesting lower outflow activity than field samples at the cosmic noon and thus, hinting at the onset of environmental effects in this massive protocluster. △ Less

Submitted 26 September, 2022; originally announced September 2022.

Comments: 29 pages, 15 figures. Accepted for publication in MNRAS

Massive Protostars in a Protocluster -- A Multi-Scale ALMA View of G35.20-0.74N

Authors: Yichen Zhang, Kei E. I. Tanaka, Jonathan C. Tan, Yao-Lun Yang, Eva Greco, Maria T. Beltrán, Nami Sakai, James M. De Buizer, Viviana Rosero, Rubén Fedriani, Guido Garay

Abstract: We present a detailed study of the massive star-forming region G35.2-0.74N with ALMA 1.3 mm multi-configuration observations. At 0.2" (440 au) resolution, the continuum emission reveals several dense cores along a filamentary structure, consistent with previous ALMA 0.85 mm observations. At 0.03" (66 au) resolution, we detect 22 compact sources, most of which are associated with the filament. Four… ▽ More We present a detailed study of the massive star-forming region G35.2-0.74N with ALMA 1.3 mm multi-configuration observations. At 0.2" (440 au) resolution, the continuum emission reveals several dense cores along a filamentary structure, consistent with previous ALMA 0.85 mm observations. At 0.03" (66 au) resolution, we detect 22 compact sources, most of which are associated with the filament. Four of the sources are associated with compact centimeter continuum emission, and two of these are associated with H30α recombination line emission. The H30α line kinematics show ordered motion of the ionized gas, consistent with disk rotation and/or outflow expansion. We construct models of photoionized regions to simultaneously fit the multi-wavelength free-free fluxes and the H30α total fluxes. The derived properties suggest the presence of at least three massive young stars with nascent hypercompact Hii regions. Two of these ionized regions are surrounded by a large rotating structure that feeds two individual disks, revealed by dense gas tracers, such as SO2, H2CO, and CH3OH. In particular, the SO2 emission highlights two spiral structures in one of the disks and probes the faster-rotating inner disks. The 12CO emission from the general region reveals a complex outflow structure, with at least four outflows identified. The remaining 18 compact sources are expected to be associated with lower-mass protostars forming in the vicinity of the massive stars. We find potential evidence for disk disruption due to dynamical interactions in the inner region of this protocluster. The spatial distribution of the sources suggests a smooth overall radial density gradient without subclustering, but with tentative evidence of primordial mass segregation. △ Less

Submitted 22 July, 2022; originally announced July 2022.

Comments: 46 pages, 21 figures, 5 tables. Accepted to ApJ

The First Detection of a Protostellar CO Outflow in the Small Magellanic Cloud with ALMA

Authors: Kazuki Tokuda, Sarolta Zahorecz, Yuri Kunitoshi, Kosuke Higashino, Kei E. I. Tanaka, Ayu Konishi, Taisei Suzuki, Naoya Kitano, Naoto Harada, Takashi Shimonishi, Naslim Neelamkodan, Yasuo Fukui, Akiko Kawamura, Toshikazu Onishi, Masahiro N. Machida

Abstract: Protostellar outflows are one of the most outstanding features of star formation. Observational studies over the last several decades have successfully demonstrated that outflows are ubiquitously associated with low- and high-mass protostars in the solar-metallicity Galactic conditions. However, the environmental dependence of protostellar outflow properties is still poorly understood, particularl… ▽ More Protostellar outflows are one of the most outstanding features of star formation. Observational studies over the last several decades have successfully demonstrated that outflows are ubiquitously associated with low- and high-mass protostars in the solar-metallicity Galactic conditions. However, the environmental dependence of protostellar outflow properties is still poorly understood, particularly in the low-metallicity regime. Here we report the first detection of a molecular outflow in the Small Magellanic Cloud with 0.2 $Z_{\odot}$, using Atacama Large Millimeter/submillimeter Array observations at a spatial resolution of 0.1 pc toward the massive protostar Y246. The bipolar outflow is nicely illustrated by high-velocity wings of CO(3-2) emission at $\gtrsim$15 km s$^{-1}$. The evaluated properties of the outflow (momentum, mechanical force, etc.) are consistent with those of the Galactic counterparts. Our results suggest that the molecular outflows, i.e., the guidepost of the disk accretion at the small scale, might be universally associated with protostars across the metallicity range of $\sim$0.2-1 $Z_{\odot}$. △ Less

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

Comments: 7 pages, 2 figures, Accepted for publication in ApJL

An AGN with an ionized gas outflow in a massive quiescent galaxy in a protocluster at $\bf z=3.09$

Authors: Mariko Kubo, Hideki Umehata, Yuichi Matsuda, Masaru Kajisawa, Charles C. Steidel, Toru Yamada, Ichi Tanaka, Bunyo Hatsukade, Yoichi Tamura, Kouichiro Nakanishi, Kotaro Kohno, Kianhong Lee, Keiichi Matsuda, Yiping Ao, Tohru Nagao, Min S. Yun

Abstract: We report the detection of an ionized gas outflow from an $X$-ray active galactic nucleus (AGN) hosted in a massive quiescent galaxy in a protocluster at $z=3.09$ (J221737.29+001823.4). It is a type-2 QSO with broad ($W_{80}>1000$ km s$^{-1}$) and strong ($\log (L_{\rm [OIII]}$ / erg s$^{-1})\approx43.4$) [O {\footnotesize III}]$λλ$4959,5007 emission lines detected by slit spectroscopy in three-po… ▽ More We report the detection of an ionized gas outflow from an $X$-ray active galactic nucleus (AGN) hosted in a massive quiescent galaxy in a protocluster at $z=3.09$ (J221737.29+001823.4). It is a type-2 QSO with broad ($W_{80}>1000$ km s$^{-1}$) and strong ($\log (L_{\rm [OIII]}$ / erg s$^{-1})\approx43.4$) [O {\footnotesize III}]$λλ$4959,5007 emission lines detected by slit spectroscopy in three-position angles using Multi-Object Infra-Red Camera and Spectrograph (MOIRCS) on the Subaru telescope and the Multi-Object Spectrometer For Infra-Red Exploration (MOSFIRE) on the Keck-I telescope. In the all slit directions, [O {\footnotesize III}] emission is extended to $\sim15$ physical kpc and indicates a powerful outflow spreading over the host galaxy. The inferred ionized gas mass outflow rate is $\rm 22\pm3~M_{\odot}~yr^{-1}$. Although it is a radio source, according to the line diagnostics using H$β$, [O {\footnotesize II}], and [O {\footnotesize III}], photoionization by the central QSO is likely the dominant ionization mechanism rather than shocks caused by radio jets. On the other hand, the spectral energy distribution of the host galaxy is well characterized as a quiescent galaxy that has shut down star formation by several hundred Myr ago. Our results suggest a scenario that QSOs are powered after the shut-down of the star formation and help to complete the quenching of massive quiescent galaxies at high redshift. △ Less

Submitted 7 July, 2022; originally announced July 2022.

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

Protostellar-disc fragmentation across all metallicities

Authors: Ryoki Matsukoba, Kei E. I. Tanaka, Kazuyuki Omukai, Eduard I. Vorobyov, Takashi Hosokawa

Abstract: Cosmic metallicity evolution possibly creates the diversity of star formation modes at different epochs. Gravitational fragmentation of circumstellar discs provides an important formation channel of multiple star systems, including close binaries. We here study the nature of disc fragmentation, systematically performing a suite of two-dimensional radiation-hydrodynamic simulations, in a broad rang… ▽ More Cosmic metallicity evolution possibly creates the diversity of star formation modes at different epochs. Gravitational fragmentation of circumstellar discs provides an important formation channel of multiple star systems, including close binaries. We here study the nature of disc fragmentation, systematically performing a suite of two-dimensional radiation-hydrodynamic simulations, in a broad range of metallicities, from the primordial to the solar values. In particular, we follow relatively long-term disc evolution over 15 kyr after the disc formation, incorporating the effect of heating by the protostellar irradiation. Our results show that the disc fragmentation occurs at all metallicities $1$--$0$ $Z_{\odot}$, yielding self-gravitating clumps. Physical properties of the clumps, such as their number and mass distributions, change with the metallicity due to different gas thermal evolution. For instance, the number of clumps is the largest for the intermediate metallicity range of $10^{-2}$--$10^{-5}$ $Z_{\odot}$, where the dust cooling is effective exclusively in a dense part of the disc and causes the fragmentation of spiral arms. The disc fragmentation is more modest for $1$--$0.1$ $Z_{\odot}$ thanks to the disc stabilization by the stellar irradiation. Such metallicity dependence agrees with the observed trend that the close binary fraction increases with decreasing metallicity in the range of $1$--$10^{-3}$ $Z_{\odot}$. △ Less

Submitted 7 June, 2022; originally announced June 2022.

Comments: 17 pages, 16 figures, submitted to MNRAS

The SOFIA Massive (SOMA) Star Formation Survey. IV. Isolated Protostars

Authors: Ruben Fedriani, Jonathan C. Tan, Zoie Telkamp, Yichen Zhang, Yao-Lun Yang, Mengyao Liu, Chi-Yan Law, Maria T. Beltran, Viviana Rosero, Kei E. I. Tanaka, Giuliana Cosentino, Prasanta Gorai, Juan Farias, Jan E. Staff, James M. De Buizer, Barbara Whitney

Abstract: We present $\sim10-40\,μ$m SOFIA-FORCAST images of 11 isolated protostars as part of the SOFIA Massive (SOMA) Star Formation Survey, with this morphological classification based on 37 $μ$m imaging. We develop an automated method to define source aperture size using the gradient of its background-subtracted enclosed flux and apply this to build spectral energy distributions (SEDs). We fit the SEDs… ▽ More We present $\sim10-40\,μ$m SOFIA-FORCAST images of 11 isolated protostars as part of the SOFIA Massive (SOMA) Star Formation Survey, with this morphological classification based on 37 $μ$m imaging. We develop an automated method to define source aperture size using the gradient of its background-subtracted enclosed flux and apply this to build spectral energy distributions (SEDs). We fit the SEDs with radiative transfer models, developed within the framework of turbulent core accretion (TCA) theory, to estimate key protostellar properties. Here, we release the sedcreator python package that carries out these methods. The SEDs are generally well fitted by the TCA models, from which we infer initial core masses $M_c$ ranging from $20-430\:M_\odot$, clump mass surface densities $Σ_{\rm cl}\sim0.3-1.7\:{\rm{g\:cm}}^{-2}$ and current protostellar masses $m_*\sim3-50\:M_\odot$. From a uniform analysis of the 40 sources in the full SOMA survey to date, we find that massive protostars form across a wide range of clump mass surface density environments, placing constraints on theories that predict a minimum threshold $Σ_{\rm cl}$ for massive star formation. However, the upper end of the $m_*-Σ_{\rm cl}$ distribution follows trends predicted by models of internal protostellar feedback that find greater star formation efficiency in higher $Σ_{\rm cl}$ conditions. We also investigate protostellar far-IR variability by comparison with IRAS data, finding no significant variation over an $\sim$40 year baseline. △ Less

Submitted 3 January, 2023; v1 submitted 23 May, 2022; originally announced May 2022.

Comments: Accepted for publication in ApJ. Published version including the proofs

High-resolution ALMA study of CO (2-1) line and dust continuum emissions in cluster galaxies at z = 1.46

Authors: Ryota Ikeda, Ken-ichi Tadaki, Daisuke Iono, Tadayuki Kodama, Jeffrey C. C. Chan, Bunyo Hatsukade, Masao Hayashi, Takuma Izumi, Kotaro Kohno, Yusei Koyama, Rhythm Shimakawa, Tomoko L. Suzuki, Yoichi Tamura, Ichi Tanaka

Abstract: We present new Atacama Large Millimeter/submillimeter Array (ALMA) results obtained from spatially resolved CO $J$=2-1 line ($0.4''$ resolution) and 870 $μ$m continuum ($0.2''$ resolution) observations of cluster galaxies in XMMXCS J2215.9-1738 at $z=1.46$. Our sample comprises 17 galaxies within $\sim0.5$ Mpc ($0.6R_{200}$) of the cluster center, all of which have previously been detected in the… ▽ More We present new Atacama Large Millimeter/submillimeter Array (ALMA) results obtained from spatially resolved CO $J$=2-1 line ($0.4''$ resolution) and 870 $μ$m continuum ($0.2''$ resolution) observations of cluster galaxies in XMMXCS J2215.9-1738 at $z=1.46$. Our sample comprises 17 galaxies within $\sim0.5$ Mpc ($0.6R_{200}$) of the cluster center, all of which have previously been detected in the CO $J$=2-1 line at a lower resolution. The effective radii of both the CO $J$=2-1 line and 870 $μ$m dust continuum emissions are robustly measured for nine galaxies by modeling the visibilities. We find that the CO $J$=2-1 line emission in all of the nine galaxies is more extended than the dust continuum emission by a factor of $2.8\pm1.4$. We investigate the spatially resolved Kennicutt-Schmidt (KS) relation in two regions within the interstellar medium of the galaxies. The relation for our sample reveals that the central region ($0<r<R_{e,{\rm 870μm}}$) of galaxies tends to have a shorter gas depletion timescale, i.e., a higher star formation efficiency, compared to the extended region ($R_{e,{\rm 870μm}}<r<R_{e,{\rm CO}}$). Overall, our result suggests that star formation activities are concentrated inside the extended gas reservoir, possibly resulting in the formation of a bulge structure. We find consistency between the ALMA 870 $μ$m radii of star-forming members and the Hubble Space Telescope/1.6 $μ$m radii of passive members in a mass-size distribution, which suggests a transition from star-forming to passive members within $\sim0.5$ Gyr. In addition, no clear differences in the KS relation nor in the sizes are found between galaxies with and without a close companion. △ Less

Submitted 23 June, 2022; v1 submitted 11 May, 2022; originally announced May 2022.

Comments: To appear in ApJ, minor revisions during proofing process, 21 pages, 10 figures, 2 tables

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

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

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

Submitted 29 April, 2022; originally announced May 2022.

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

Vibrationally-excited Lines of HC$_{3}$N Associated with the Molecular Disk around the G24.78+0.08 A1 Hyper-compact H$_{\rm {II}}$ Region

Authors: Kotomi Taniguchi, Kei E. I. Tanaka, Yichen Zhang, Rubén Fedriani, Jonathan C. Tan, Shigehisa Takakuwa, Fumitaka Nakamura, Masao Saito, Liton Majumdar, Eric Herbst

Abstract: We have analyzed Atacama Large Millimeter/submillimeter Array Band 6 data of the hyper-compact H$_{\rm {II}}$ region G24.78+0.08 A1 (G24 HC H$_{\rm {II}}$) and report the detection of vibrationally-excited lines of HC$_{3}$N ($v_{7}=2$, $J=24-23$). The spatial distribution and kinematics of a vibrationally-excited line of HC$_{3}$N ($v_{7}=2$, $J=24-23$, $l=2e$) are found to be similar to the CH… ▽ More We have analyzed Atacama Large Millimeter/submillimeter Array Band 6 data of the hyper-compact H$_{\rm {II}}$ region G24.78+0.08 A1 (G24 HC H$_{\rm {II}}$) and report the detection of vibrationally-excited lines of HC$_{3}$N ($v_{7}=2$, $J=24-23$). The spatial distribution and kinematics of a vibrationally-excited line of HC$_{3}$N ($v_{7}=2$, $J=24-23$, $l=2e$) are found to be similar to the CH$_{3}$CN vibrationally-excited line ($v_{8}=1$), which indicates that the HC$_{3}$N emission is tracing the disk around the G24 HC H$_{\rm {II}}$ region previously identified by the CH$_{3}$CN lines. We derive the $^{13}$CH$_{3}$CN/HC$^{13}$CCN abundance ratios around G24 and compare them to the CH$_{3}$CN/HC$_{3}$N abundance ratios in disks around Herbig Ae and T Tauri stars. The $^{13}$CH$_{3}$CN/HC$^{13}$CCN ratios around G24 ($\sim 3.0-3.5$) are higher than the CH$_{3}$CN/HC$_{3}$N ratios in the other disks ($\sim 0.03-0.11$) by more than one order of magnitude. The higher CH$_{3}$CN/HC$_{3}$N ratios around G24 suggest that the thermal desorption of CH$_{3}$CN in the hot dense gas and efficient destruction of HC$_{3}$N in the region irradiated by the strong UV radiation are occurring. Our results indicate that the vibrationally-excited HC$_{3}$N lines can be used as a disk tracer of massive protostars at the HC H$_{\rm {II}}$ region stage, and the combination of these nitrile species will provide information of not only chemistry but also physical conditions of the disk structures. △ Less

Submitted 24 April, 2022; v1 submitted 21 April, 2022; originally announced April 2022.

Comments: 21 pages, 13 figures, 5 tables, Accepted by The Astrophysical Journal

Investigation of superconductivity in Ce-doped (La,Pr)OBiS2 single crystals

Authors: Masanori Nagao, Yuji Hanada, Akira Miura, Yuki Maruyama, Satoshi Watauchi, Yoshihiko Takano, Isao Tanaka

Abstract: Single crystals of Ce-doped (La,Pr)OBiS2 superconductors, multinary rare-earth elements substituted ROBiS2, were successfully grown. The grown crystals typically had a size of 1-2 mm and a plate-like shape with a well-developed c-plane. The c-axis lattice constants of the obtained (La,Ce,Pr)OBiS2 single crystals were approximately 13.6-13.7 A, and the superconducting transition temperature was 1.2… ▽ More Single crystals of Ce-doped (La,Pr)OBiS2 superconductors, multinary rare-earth elements substituted ROBiS2, were successfully grown. The grown crystals typically had a size of 1-2 mm and a plate-like shape with a well-developed c-plane. The c-axis lattice constants of the obtained (La,Ce,Pr)OBiS2 single crystals were approximately 13.6-13.7 A, and the superconducting transition temperature was 1.23-2.18 K. Valence fluctuations of Ce and Pr were detected through X-ray absorption spectroscopy analysis. In contrast to (Ce,Pr)OBiS2 and (La,Ce)OBiS2, the superconducting transition temperature of (La,Ce,Pr)OBiS2 increased with increasing concentrations of the tetravalent state at the R-site. △ Less

Submitted 19 April, 2022; originally announced April 2022.

Journal ref: Materials, vol.15 (2022) 2977

LO-mode phonon of KCl and NaCl at 300 K by inelastic X ray scattering measurements and first principles calculations

Authors: Atsushi Togo, Hiroyuki Hayashi, Terumasa Tadano, Satoshi Tsutsui, Isao Tanaka

Abstract: Longitudinal-optical (LO) mode phonon branches of KCl and NaCl were measured using inelastic X-ray scattering (IXS) at 300 K and calculated by the first-principles phonon calculation with the stochastic self-consistent harmonic approximation. Spectral shapes of the IXS measurements and calculated spectral functions agreed well. We analyzed the calculated spectral functions that provide higher reso… ▽ More Longitudinal-optical (LO) mode phonon branches of KCl and NaCl were measured using inelastic X-ray scattering (IXS) at 300 K and calculated by the first-principles phonon calculation with the stochastic self-consistent harmonic approximation. Spectral shapes of the IXS measurements and calculated spectral functions agreed well. We analyzed the calculated spectral functions that provide higher resolutions of the spectra than the IXS measurements. Due to strong anharmonicity, the spectral functions of these phonon branches have several peaks and the LO modes along $Γ$--L paths are disconnected. △ Less

Submitted 22 April, 2022; v1 submitted 15 April, 2022; originally announced April 2022.

Star Formation Properties of Sloan Digital Sky Survey BOSS Void Galaxies in the Hyper Suprime-Cam Survey

Authors: Hung-Yu Jian, Lihwai Lin, Bau-Ching Hsieh, Kai-Yang Lin, Keiichi Umetsu, Carlos Lopez-Coba, Yusei Koyama, Chin-Hao Hsu, Yung-Chau Su, Yu-Yen Chang, Tadayuki Kodama, Yutaka Komiyama, Surhud More, Atsushi J. Nishizawa, Masamune Oguri, Ichi Tanaka

Abstract: We utilize the Hyper Suprime-Cam (HSC) Wide Survey to explore the properties of galaxies located in the voids identified from the Baryon Oscillation Spectroscopic Survey (BOSS) up to z~0.7. The HSC reaches i~25, allowing us to characterize the void galaxies down to 10$^{9.2}$ solar mass. We find that the revised void galaxy densities, when including faint galaxies in voids defined by bright galaxi… ▽ More We utilize the Hyper Suprime-Cam (HSC) Wide Survey to explore the properties of galaxies located in the voids identified from the Baryon Oscillation Spectroscopic Survey (BOSS) up to z~0.7. The HSC reaches i~25, allowing us to characterize the void galaxies down to 10$^{9.2}$ solar mass. We find that the revised void galaxy densities, when including faint galaxies in voids defined by bright galaxies, are still underdense compared to the mean density from the entire field. In addition, we classify galaxies into star-forming, quiescent, and green valley populations, and find that void galaxies tend to have slightly higher fractions of star-forming galaxies under the mass and redshift control, although the significance of this result is only moderate (2$σ$). However, when we focus on the star-forming population, the distribution of the specific star formation rate (sSFR) of void galaxies shows little difference from that of the control galaxies. Similarly, the median sSFR of star-forming void galaxies is also in good agreement with that of the star-forming control galaxies. Moreover, the effective green valley fraction of void galaxies, defined as the number of green valley galaxies over the number of nonquiescent galaxies, is comparable to that of the control ones, supporting the suggestion that void and control galaxies evolve under similar physical processes and quenching frequencies. Our results thus favor a scenario of the galaxy assembly bias. △ Less

Submitted 10 January, 2022; originally announced January 2022.

Comments: 17 pages, 10 figures, 3 tables, ApJ accepted

Supersonic Expansion of the Bipolar Hii Region Sh2-106: A 3,500 Year-Old Explosion?

Authors: John Bally, Zen Chia, Adam Ginsburg, Bo Reipurth, Kei E. I. Tanaka, Hans Zinnecker, John Faulhaber

Abstract: Multi-epoch narrow-band HST images of the bipolar Hii region Sh2-106 reveal highly supersonic nebular proper motions which increase with projected distance from the massive young stellar object S106~IR, reaching over ~30 mas/year (~150 km/s at D=1.09 kpc) at a projected separation of ~1.4' (0.44 pc) from S106~IR. We propose that S106~IR experienced a $\sim10^{47}$ erg explosion ~3,500 years ago. T… ▽ More Multi-epoch narrow-band HST images of the bipolar Hii region Sh2-106 reveal highly supersonic nebular proper motions which increase with projected distance from the massive young stellar object S106~IR, reaching over ~30 mas/year (~150 km/s at D=1.09 kpc) at a projected separation of ~1.4' (0.44 pc) from S106~IR. We propose that S106~IR experienced a $\sim10^{47}$ erg explosion ~3,500 years ago. The explosion may be the result of a major accretion burst, a recent encounter with another star, or a consequence of the interaction of a companion with the bloated photosphere of S106~IR as it grew from ~10 through ~15 Solar masses at a high accretion rate. Near-IR images reveal fingers of molecular hydrogen emission pointing away from S106~IR and an asymmetric photon-dominated region surrounding the ionized nebula. Radio continuum and Brackett-gamma emission reveal a C-shaped bend in the plasma, either indicating motion of S106~IR toward the east, or deflection of plasma toward the west by the surrounding cloud. The Hii region bends around a ~1' diameter dark bay west of S106~IR that may be shielded from direct illumination by a dense molecular clump. Herbig-Haro (HH) and Molecular Hydrogen Objects (MHOs) tracing outflows powered by stars in the Sh2-106 proto-cluster such as the Class 0 source S106 FIR are discussed. △ Less

Submitted 11 November, 2021; originally announced November 2021.

Comments: 47 pages, 24 figures

Intrinsic suppression of topological thermal Hall effect in an exactly solvable quantum magnet

Authors: S. Suetsugu, T. Yokoi, K. Totsuka, T. Ono, I. Tanaka, S. Kasahara, Y. Kasahara, Z. Chengchao, H. Kageyama, Y. Matsuda

Abstract: In contrast to electron (fermion) systems, topological phases of charge neutral bosons have been poorly understood despite recent extensive research on insulating magnets. The most important unresolved issue is how the inevitable inter-bosonic interactions influence the topological properties. It has been proposed that the quantum magnet SrCu$_2$(BO$_3$)$_2$ with an exact ground state serves as an… ▽ More In contrast to electron (fermion) systems, topological phases of charge neutral bosons have been poorly understood despite recent extensive research on insulating magnets. The most important unresolved issue is how the inevitable inter-bosonic interactions influence the topological properties. It has been proposed that the quantum magnet SrCu$_2$(BO$_3$)$_2$ with an exact ground state serves as an ideal platform for this investigation, as the system is expected to be a magnetic analogue of a Chern insulator with electrons replaced by bosonic magnetic excitations (triplons). Here, in order to examine topologically protected triplon chiral edge modes in SrCu$_2$(BO$_3$)$_2$, we measured and calculated the thermal Hall conductivity $κ_{xy}$. Our calculations show that the sign of $κ_{xy}$ is negative, which is opposite to the previous calculations, and its magnitude is 2$π$ times smaller. No discernible $κ_{xy}$ was observed, which is at most 20-30% of our calculations if present. This implies that even relatively weak inter-particle interactions seriously influence the topological transport properties at finite temperatures. These demonstrate that, in contrast to fermionic cases, the picture of non-interacting topological quasi-particles cannot be naively applied to bosonic systems, calling special attention to the interpretation of the topological bosonic excitations reported for various insulating magnets. △ Less

Submitted 11 January, 2022; v1 submitted 26 July, 2021; originally announced July 2021.

Comments: 11 pages, 11 figures

Journal ref: Phys. Rev. B 105, 024415 (2022)

Finding well-optimized special quasirandom structures with decision diagram

Authors: Kohei Shinohara, Atsuto Seko, Takashi Horiyama, Isao Tanaka

Abstract: The advanced data structure of the zero-suppressed binary decision diagram (ZDD) enables us to efficiently enumerate nonequivalent substitutional structures. Not only can the ZDD store a vast number of structures in a compressed manner, but also can a set of structures satisfying given constraints be extracted from the ZDD efficiently. Here, we present a ZDD-based efficient algorithm for exhaustiv… ▽ More The advanced data structure of the zero-suppressed binary decision diagram (ZDD) enables us to efficiently enumerate nonequivalent substitutional structures. Not only can the ZDD store a vast number of structures in a compressed manner, but also can a set of structures satisfying given constraints be extracted from the ZDD efficiently. Here, we present a ZDD-based efficient algorithm for exhaustively searching for special quasirandom structures (SQSs) that mimic the perfectly random substitutional structure. We demonstrate that the current approach can extract only a tiny number of SQSs from a ZDD composed of many substitutional structures (>$10^{12}$). As a result, we find SQSs that are optimized better than those proposed in the literature. A series of SQSs should be helpful for estimating the properties of substitutional solid solutions. Furthermore, the present ZDD-based algorithm should be useful for applying the ZDD to the other structure enumeration problems. △ Less

Submitted 29 November, 2021; v1 submitted 15 July, 2021; originally announced July 2021.

A massive quiescent galaxy confirmed in a protocluster at z=3.09

Authors: Mariko Kubo, Hideki Umehata, Yuichi Matsuda, Masaru Kajisawa, Charles C. Steidel, Toru Yamada, Ichi Tanaka, Bunyo Hatsukade, Yoichi Tamura, Kouichiro Nakanishi, Kotaro Kohno, Chien-Feng Lee, Keiichi Matsuda

Abstract: We report a massive quiescent galaxy at $z_{\rm spec}=3.0922^{+0.008}_{-0.004}$ spectroscopically confirmed at a protocluster in the SSA22 field by detecting the Balmer and Ca {\footnotesize II} absorption features with multi-object spectrometer for infrared exploration (MOSFIRE) on the Keck I telescope. This is the most distant quiescent galaxy confirmed in a protocluster to date. We fit the opti… ▽ More We report a massive quiescent galaxy at $z_{\rm spec}=3.0922^{+0.008}_{-0.004}$ spectroscopically confirmed at a protocluster in the SSA22 field by detecting the Balmer and Ca {\footnotesize II} absorption features with multi-object spectrometer for infrared exploration (MOSFIRE) on the Keck I telescope. This is the most distant quiescent galaxy confirmed in a protocluster to date. We fit the optical to mid-infrared photometry and spectrum simultaneously with spectral energy distribution (SED) models of parametric and nonparametric star formation histories (SFH). Both models fit the observed SED well and confirm that this object is a massive quiescent galaxy with the stellar mass of $\log(\rm M_{\star}/M_{\odot}) = 11.26^{+0.03}_{-0.04}$ and $11.54^{+0.03}_{-0.00}$, and star formation rate of $\rm SFR/M_{\odot}~yr^{-1} <0.3$ and $=0.01^{+0.03}_{-0.01}$ for parametric and nonparametric models, respectively. The SFH from the former modeling is described as an instantaneous starburst while that of the latter modeling is longer-lived but both models agree with a sudden quenching of the star formation at $\sim0.6$ Gyr ago. This massive quiescent galaxy is confirmed in an extremely dense group of galaxies predicted as a progenitor of a brightest cluster galaxy formed via multiple mergers in cosmological numerical simulations. We newly find three plausible [O III]$λ$5007 emitters at $3.0791\leq z_{\rm spec}\leq3.0833$ happened to be detected around the target. Two of them just between the target and its nearest massive galaxy are possible evidence of their interactions. They suggest the future strong size and stellar mass evolution of this massive quiescent galaxy via mergers. △ Less

Submitted 20 June, 2021; originally announced June 2021.

Comments: 13 pages, 7 figures (in main text. 21 pages including appendices), accepted for publication in ApJ

Fluorine solubility and superconducting properties of Sm(O,F)BiS$_2$ single crystals

Authors: Koki Kinami, Yuji Hanada, Masanori Nagao, Akira Miura, Shigeto Hirai, Yuki Maruyama, Satoshi Watauchi, Isao Tanaka

Abstract: Sm(O,F)BiS2 superconducting single crystals with various nominal F contents have been grown using KI-KCl flux. The solid solution limit of F at the O-site in Sm(O,F)BiS2 single crystals was approximately less than 20 at%. F concentrations of Sm(O,F)BiS2 single crystals were hardly controlled by nominal F contents. It suggests the existence of intrinsic stable phase. Superconductivity was shown at… ▽ More Sm(O,F)BiS2 superconducting single crystals with various nominal F contents have been grown using KI-KCl flux. The solid solution limit of F at the O-site in Sm(O,F)BiS2 single crystals was approximately less than 20 at%. F concentrations of Sm(O,F)BiS2 single crystals were hardly controlled by nominal F contents. It suggests the existence of intrinsic stable phase. Superconductivity was shown at around 5 K. Superconducting transition temperature with zero resistivity were increased with increasing nominal F contents within less than 70 at% in O-site, despite similar F concentration and c-axis lattice constants. △ Less

Submitted 17 June, 2021; originally announced June 2021.

Journal ref: Journal of Alloys and Compounds vol.883 (2021) 160812

Nematicity Liquid in a Trimerized-Kagome Antiferromagnet

Authors: Ichiro Tanaka, Hirokazu Tsunetsugu

Abstract: We theoretically study low-temperature properties of the antiferromagnetic Heisenberg model with half-integer spin $S$ on the Kagome lattice with large trimerization. We have derived a low-energy effective model for general $S$ in terms of spin and nematicity operators in triangular units, and studied their low-temperature correlations for S=3/2 by classical Monte Carlo simulations. The previous s… ▽ More We theoretically study low-temperature properties of the antiferromagnetic Heisenberg model with half-integer spin $S$ on the Kagome lattice with large trimerization. We have derived a low-energy effective model for general $S$ in terms of spin and nematicity operators in triangular units, and studied their low-temperature correlations for S=3/2 by classical Monte Carlo simulations. The previous study for the S=1/2 case [M. Ferrero et al., Phys. Rev. B 68, 214431 (2003)] reported a spin liquid state at low temperatures driven by a glassy behavior of isolated dimers and trimers of nematicities. The results for S=3/2 show that nematicity dimers and trimers are connected by weak links to form a defective planar network. At very low temperatures, nematicities show glassy slow dynamics, and cluster spins continue to fluctuate in a nonperiodically frozen background of nematicities. The characteristic time of glassy dynamics scales with temperature following a power law. △ Less

Submitted 20 May, 2021; v1 submitted 22 April, 2021; originally announced April 2021.

Comments: 5 pages, 6 figures

TXS 1206+549: a new $γ$-ray detected narrow-line Seyfert 1 galaxy at redshift 1.34?

Authors: Suvendu Rakshit, Malte Schramm, C. S. Stalin, I. Tanaka, Vaidehi S. Paliya, Indrani Pal, Jari Kotilainen, Jaejin Shin

Abstract: Radio and $γ$-ray loud narrow-line Seyfert 1 galaxies (NLS1s) are unique objects to study the formation and evolution of relativistic jets, as they are believed to have high accretion rates and powered by low mass black holes contrary to that known for blazars. However, only about a dozen $γ$-ray detected NLS1s ($γ$-NLS1s) are known to date and all of them are at $z\le1$. Here, we report the ident… ▽ More Radio and $γ$-ray loud narrow-line Seyfert 1 galaxies (NLS1s) are unique objects to study the formation and evolution of relativistic jets, as they are believed to have high accretion rates and powered by low mass black holes contrary to that known for blazars. However, only about a dozen $γ$-ray detected NLS1s ($γ$-NLS1s) are known to date and all of them are at $z\le1$. Here, we report the identification of a new $γ$-ray emitting NLS1 TXS 1206+549 at $z=1.344$. A near-infrared spectrum taken with the Subaru telescope showed H$β$ emission line with FWHM of $1194\pm77$ km s$^{-1}$ and weak [O III] emission line but no optical Fe II lines, due to the limited wavelength coverage and poor signal-to-noise ratio. However, UV Fe II lines are present in the SDSS optical spectrum. The source is very radio-loud, unresolved, and has a flat radio spectrum. The broadband SED of the source has the typical two hump structure shown by blazars and other $γ$-NLS1s. The source exhibits strong variability at all wavelengths such as the optical, infrared, and $γ$-ray bands. All these observed characteristics show that TXS 1206+549 is the most distant $γ$-NLS1 known to date. △ Less

Submitted 31 March, 2021; v1 submitted 27 March, 2021; originally announced March 2021.

Comments: Accepted in MNRAS Letter

The Hot and Dynamic Birth of Massive Stars from the ngVLA Perspective

Authors: Kei E. I. Tanaka, Yichen Zhang, Kazuhito Motogi

Abstract: The Next Generation Very Large Array (ngVLA) has excellent capabilities to unveil various dynamical and chemical processes in massive star formation at the unexplored innermost regions. Based on the recent observations of ALMA/VLA as well as theoretical predictions, we propose several intriguing topics in massive star formation from the perspective of the ngVLA. In the disk scale of $\lesssim$ 100… ▽ More The Next Generation Very Large Array (ngVLA) has excellent capabilities to unveil various dynamical and chemical processes in massive star formation at the unexplored innermost regions. Based on the recent observations of ALMA/VLA as well as theoretical predictions, we propose several intriguing topics in massive star formation from the perspective of the ngVLA. In the disk scale of $\lesssim$ 100 au around massive protostars, dust grains are expected to be destructed/sublimated because the physical conditions of temperature, shocks, and radiation are much more intense than those in the envelopes, which are typically observed as hot cores. The high sensitivity and resolution of the ngVLA will enable us to detect the gaseous refractories released by dust destruction, e.g., SiO, NaCl, and AlO, which trace disk kinematics and give new insights into the metallic elements in star-forming regions, i.e., astromineralogy. The multi-epoch survey by the ngVLA will provide demographics of forming massive multiples with separations of $\lesssim$ 10 au with their proper motion. Combining with observations of refractory molecular lines and hydrogen recombination lines, we can reproduce the three-dimensional orbital motions of massive proto-binaries. Moreover, the 1-mas resolution of the ngVLA could possibly take the first-ever picture of the photospheric surface of an accreting protostar, if it is bloated to the au scale by the high accretion rates of mass and thermal energy. △ Less

Submitted 16 March, 2021; originally announced March 2021.

Comments: 7 pages, 4 figures, Published in the ngVLA-Japan Memo Series

Supernova Model Discrimination with Hyper-Kamiokande

Authors: Hyper-Kamiokande Collaboration, :, K. Abe, P. Adrich, H. Aihara, R. Akutsu, I. Alekseev, A. Ali, F. Ameli, I. Anghel, L. H. V. Anthony, M. Antonova, A. Araya, Y. Asaoka, Y. Ashida, V. Aushev, F. Ballester, I. Bandac, M. Barbi, G. J. Barker, G. Barr, M. Batkiewicz-Kwasniak, M. Bellato, V. Berardi, M. Bergevin , et al. (478 additional authors not shown)

Abstract: Core-collapse supernovae are among the most magnificent events in the observable universe. They produce many of the chemical elements necessary for life to exist and their remnants -- neutron stars and black holes -- are interesting astrophysical objects in their own right. However, despite millennia of observations and almost a century of astrophysical study, the explosion mechanism of core-colla… ▽ More Core-collapse supernovae are among the most magnificent events in the observable universe. They produce many of the chemical elements necessary for life to exist and their remnants -- neutron stars and black holes -- are interesting astrophysical objects in their own right. However, despite millennia of observations and almost a century of astrophysical study, the explosion mechanism of core-collapse supernovae is not yet well understood. Hyper-Kamiokande is a next-generation neutrino detector that will be able to observe the neutrino flux from the next galactic core-collapse supernova in unprecedented detail. We focus on the first 500 ms of the neutrino burst, corresponding to the accretion phase, and use a newly-developed, high-precision supernova event generator to simulate Hyper-Kamiokande's response to five different supernova models. We show that Hyper-Kamiokande will be able to distinguish between these models with high accuracy for a supernova at a distance of up to 100 kpc. Once the next galactic supernova happens, this ability will be a powerful tool for guiding simulations towards a precise reproduction of the explosion mechanism observed in nature. △ Less

Submitted 20 July, 2021; v1 submitted 13 January, 2021; originally announced January 2021.

Comments: 21 pages, 7 figures. Article based on thesis published as arXiv:2002.01649. v2: added references and some explanations in response to reviewer comments

Journal ref: Astrophys.J. 916 (2021) 15