Search

Fragmentation contributes to the formation and evolution of stars. Observationally, high-mass stars are known to form multiple-star systems, preferentially in cluster environments. Theoretically, Jeans instability has been suggested to determine characteristic fragmentation scales, and thermal or turbulent motion in the parental gas clump mainly contributes to the instability. To search for such a characteristic fragmentation scale, we have analyzed ALMA 1.33 mm continuum observations toward 30 high-mass star-forming clumps taken by the Digging into the Interior of Hot Cores with ALMA (DIHCA) survey. We have identified 573 cores using the dendrogram algorithm and measured the separation of cores by using the Minimum Spanning Tree (MST) technique. The core separation corrected by projection effects has a distribution peaked around 5800 au. In order to remove biases produced by different distances and sensitivities, we further smooth the images to a common physical scale and perform completeness tests. Our careful analysis finds a characteristic fragmentation scale of $\sim$7000 au, comparable to the thermal Jeans length of the clumps. We conclude that thermal Jeans fragmentation plays a dominant role in determining the clump fragmentation in high-mass star-forming regions, without the need of invoking turbulent Jeans fragmentation., Comment: 30 pages, 18 figures, Accepted in ApJ

For efficient spectroscopic redshift identification of early galaxies in the northern hemisphere, we aim to combine the Large Millimeter Telescope (LMT) with a wide-band heterodyne receiver, FINER, which will cover radio frequencies of 120--360 GHz and offer a 3--21 GHz intermediate frequency (IF) per sideband and polarization. To take full advantage of such wide IFs, we present a novel 10.24-GHz-wide digital spectrometer, DRS4 (Elecs Industry Co., Ltd.). It incorporates 20.48 Gsps samplers with an FPGA-based digital signal processing module. To mitigate the noise contamination from the image sideband, it is equipped with a digital sideband separation function to improve the sideband rejection up to 25 dB. Laboratory performance evaluations show that it exhibits an Allan time of at least ~100 s and a total power dynamic range of at least 7 dB. These results demonstrate its capability of instantaneously wide-band spectroscopy toward high-redshift galaxies with position-switching observations., Comment: 10 pages, 7 figures, Proceedings paper presented in SPIE Astronomical Telescope and Instrumentation 2024

Unveiling the emergence and prevalence of massive/bright galaxies during the epoch of reionization and beyond, within the first 600 million years of the Universe, stands as a pivotal pursuit in astronomy. Remarkable progress has been made by JWST in identifying an immense population of bright galaxies, which hints at exceptionally efficient galaxy assembly processes. However, the underlying physical mechanisms propelling their rapid growth remain unclear. With this in mind, millimeter and submillimeter-wave spectroscopic observations of redshifted far-infrared spectral lines, particularly the [O III] 88 micron and [C II] 158 micron lines, offers a crucial pathway to address this fundamental query. To this end, we develop a dual-polarization sideband-separating superconductor-insulator-superconductor (SIS) mixer receiver, FINER, for the Large Millimeter Telescope (LMT) situated in Mexico. Harnessing advancements from ALMA's wideband sensitivity upgrade (WSU) technology, FINER covers radio frequencies spanning 120-360 GHz, delivering an instantaneous intermediate frequency (IF) of 3-21 GHz per sideband per polarization, which is followed by a set of 10.24 GHz-wide digital spectrometers. At 40% of ALMA's light-collecting area, the LMT's similar atmospheric transmittance and FINER's 5 times wider bandwidth compared to ALMA culminate in an unparalleled spectral scanning capability in the northern hemisphere, paving the way for finer spectral-resolution detection of distant galaxies., Comment: 12 pages, 8 figures, and 3 tables. Proceedings paper presented in SPIE Astronomical Telescope and Instrumentation 2024

Spectroscopic observations of the far-infrared [O III] and [C II] lines present a pathway to explore the mechanisms of the emergence of massive galaxies in the epoch of reionization and beyond, which is one of the most fundamental questions in astronomy. To address this question, the Far-Infrared Nebular Emission Receiver (FINER) project is developing two wideband dual-polarization sideband-separating heterodyne receivers at 120--210 GHz and 210--360 GHz for the Large Millimeter Telescope (LMT) in Mexico. Compared with Atacama Large Millimeter/submillimeter Array (ALMA), LMT provides 40% of ALMA's light-collecting area and a similar atmospheric transmittance, but FINER plans to have an instantaneous intermediate frequency (IF) of 3--21 GHz per sideband per polarization which is five times wider than current ALMA's bandwidth. Therefore, FINER is going to offer cutting-edge spectral scanning capability in the next several years. The project is currently in an active development phase. In this proceeding, the latest development status for FINER, including the optics, wideband waveguide components as well as low-noise superconductor-insulator-superconductor (SIS) mixers is reported., Comment: 14 pages, 8 figures, 1 table. Proceeding paper presented in SPIE Astronomical Telescope and Instrumentation 2024

The exploration of outflows in protobinary systems presents a challenging yet crucial endeavour, offering valuable insights into the dynamic interplay between protostars and their evolution. In this study, we examine the morphology and dynamics of jets and outflows within the IRAS\,4A protobinary system. This analysis is based on ALMA observations of SiO(5--4), H$_2$CO(3$_{0,3}$--2$_{0,3}$), and HDCO(4$_{1,4}$--3$_{1,3}$) with a spatial resolution of $\sim$150\,au. Leveraging an astrochemical approach involving the use of diverse tracers beyond traditional ones has enabled the identification of novel features and a comprehensive understanding of the broader outflow dynamics. Our analysis reveals the presence of two jets in the redshifted emission, emanating from IRAS\,4A1 and IRAS\,4A2, respectively. Furthermore, we identify four distinct outflows in the region for the first time, with each protostar, 4A1 and 4A2, contributing to two of them. We characterise the morphology and orientation of each outflow, challenging previous suggestions of bends in their trajectories. The outflow cavities of IRAS\,4A1 exhibit extensions of 10$''$ and 13$''$ with position angles (PA) of 0$^{\circ}$ and -12$^{\circ}$, respectively, while those of IRAS\,4A2 are more extended, spanning 18$''$ and 25$''$ with PAs of 29$^{\circ}$ and 26$^{\circ}$. We propose that the misalignment of the cavities is due to a jet precession in each protostar, a notion supported by the observation that the more extended cavities of the same source exhibit lower velocities, indicating they may stem from older ejection events.

Fragmentation during the early stages of high-mass star formation is crucial for understanding the formation of high-mass clusters. We investigated fragmentation within thirty-nine high-mass star-forming clumps as part of the Atacama Large Millimeter/submillimeter Array (ALMA) Survey of 70 $\mu$m Dark High-mass Clumps in Early Stages (ASHES). Considering projection effects, we have estimated core separations for 839 cores identified from the continuum emission and found mean values between 0.08 and 0.32 pc within each clump. We find compatibility of the observed core separations and masses with the thermal Jeans length and mass, respectively. We also present sub-clump structures revealed by the 7 m-array continuum emission. Comparison of the Jeans parameters using clump and sub-clump densities with the separation and masses of gravitationally bound cores suggests that they can be explained by clump fragmentation, implying the simultaneous formation of sub-clumps and cores within rather than a step-by-step hierarchical fragmentation. The number of cores in each clump positively correlates with the clump surface density and the number expected from the thermal Jeans fragmentation. We also find that the higher the fraction of protostellar cores, the larger the dynamic range of the core mass, implying that the cores are growing in mass as the clump evolves. The ASHES sample exhibits various fragmentation patterns: aligned, scattered, clustered, and sub-clustered. Using the Q-parameter, which can help to distinguish between centrally condensed and subclustered spatial core distributions, we finally find that in the early evolutionary stages of high-mass star formation, cores tend to follow a subclustered distribution., Comment: Accepted for Publication in ApJ. 19 pages, 7 figures, 4 tables

The dominant mechanism forming multiple stellar systems in the high-mass regime (M$_\ast \gtrsim $ 8 $M_{\odot}$) remained unknown because direct imaging of multiple protostellar systems at early phases of high-mass star formation is very challenging. High-mass stars are expected to form in clustered environments containing binaries and higher-order multiplicity systems. So far only a few high-mass protobinary systems, and no definitive higher-order multiples, have been detected. Here we report the discovery of one quintuple, one quadruple, one triple and four binary protostellar systems simultaneously forming in a single high-mass protocluster, G333.23--0.06, using Atacama Large Millimeter/submillimeter Array high-resolution observations. We present a new example of a group of gravitationally bound binary and higher-order multiples during their early formation phases in a protocluster. This provides the clearest direct measurement of the initial configuration of primordial high-order multiple systems, with implications for the in situ multiplicity and its origin. We find that the binary and higher-order multiple systems, and their parent cores, show no obvious sign of disk-like kinematic structure. We conclude that the observed fragmentation into binary and higher-order multiple systems can be explained by core fragmentation, indicating its crucial role in establishing the multiplicity during high-mass star cluster formation., Comment: Published in Nature Astronomy at https://www.nature.com/articles/s41550-023-02181-9. 27 pages, 12 figures, 1 table

Massive infrared dark clouds (IRDCs) are considered to host the earliest stages of high-mass star formation. In particular, 70 $\mu$m dark IRDCs are the colder and more quiescent clouds. At a scale of about 5000 au using formaldehyde (H2CO) emission, we investigate the kinetic temperature of dense cores in 12 IRDCs obtained from the pilot ALMA Survey of 70 $\mu$m dark High-mass clumps in Early Stages (ASHES). Compared to 1.3 mm dust continuum and other molecular lines, such as C18O and deuterated species, we find that H2CO is mainly sensitive to low-velocity outflow components rather than to quiescent gas expected in the early phases of star formation. The kinetic temperatures of these components range from 26 to 300 K. The Mach number reaches about 15 with an average value of about 4, suggesting that the velocity distribution of gas traced by H2CO is significantly influenced by a supersonic non-thermal component. In addition, we detect warm line emission from HC3N and OCS in 14 protostellar cores, which requires high excitation temperatures (Eu/k ~ 100 K). These results show that some of the embedded cores in the ASHES fields are in an advanced evolutionary stage, previously unexpected for 70 $\mu$m dark IRDCs., Comment: Accepted for Publication in ApJ. 39 pages, 22 figures, 6 tables

Context. Deuterium in H-bearing species is enhanced during the early stages of star formation, however, only a small number of high spatial resolution deuteration studies exist towards protostellar objects, leaving the small-scale structures unrevealed and understudied. Aims. We aim to constrain the deuterium fractionation ratios in a Class 0/I protostellar object in formaldehyde (H2CO), which has abundant deuterated isotopologues in this environment. Methods. We observed the Class 0/I protobinary system [BHB2007] 11, whose emission components are embedded in circumstellar disks that have radii of 2-3 au, using ALMA within the context of the Large Program FAUST. The system is surrounded by a complex filamentary structure connecting to the larger circumbinary disk. In this work we present the first study of formaldehyde D-fractionation towards this source with detections of H2CO 3(0,3)-2(0,2), combined with HDCO 4(2,2)-3(2,1), HDCO 4(1,4)-3(1,3) and D2CO 4(0,4)-3(0,3). These observations enable multiple velocity components associated with the methanol hotspots also uncovered by FAUST data, as well as the external envelope, to be resolved. In addition, based on the kinematics seen in the observations of the H2CO emission, we propose the presence of a second large scale outflow. Results. HDCO and D2CO are only found in the central regions of the core while H2CO is found more ubiquitously. From radiative transfer modelling, the column densities ranges found for H2CO, HDCO and D2CO are (3-8)x10$^{14}$ cm$^{-2}$, (0.8-2.9)x10$^{13}$ cm$^{-2}$ and (2.6-4.3)x10$^{12}$ cm$^{-2}$, respectively, yielding an average D/H ratio of 0.01-0.04. Following the results of kinematic modelling, the second large scale feature is inconsistent with a streamer-like nature and we thus tentatively conclude that the feature is an asymmetric molecular outflow launched by a wide-angle disk wind., Comment: 17 pages, 15 figures

The initial conditions found in infrared dark clouds (IRDCs) provide insights on how high-mass stars and stellar clusters form. We have conducted high-angular resolution and high-sensitivity observations toward thirty-nine massive IRDC clumps, which have been mosaicked using the 12m and 7m arrays from the Atacama Large Millimeter/submillimeter Array (ALMA). The targets are 70 $\mu$m dark massive (220-4900 $M_\odot$), dense ($>$10$^4$ cm$^{-3}$), and cold ($\sim$10-20K) clumps located at distances between 2 and 6 kpc. We identify an unprecedented number of 839 cores, with masses between 0.05 and 81 $M_\odot$ using 1.3 mm dust continuum emission. About 55% of the cores are low-mass ($<$1 $M_\odot$), whereas $\lesssim$1% (7/839) are high-mass ($\gtrsim$27 $M_\odot$). We detect no high-mass prestellar cores. The most massive cores (MMC) identified within individual clumps lack sufficient mass to form high-mass stars without additional mass feeding. We find that the mass of the MMCs is correlated with the clump surface density, implying denser clumps produce more massive cores and a larger number of cores. There is no significant mass segregation except for a few tentative detections. In contrast, most clumps show segregation once the clump density is considered instead of mass. Although the dust continuum emission resolves clumps in a network of filaments, some of which consist of hub-filament systems, the majority of the MMCs are not found in the hubs. Our analysis shows that high-mass cores and MMCs have no preferred location with respect to low-mass cores at the earliest stages of high-mass star formation., Comment: Accepted for Publication in ApJ. 54 pages, 40 figures, 5 tables

We present a study of chemistry toward 294 dense cores in 12 molecular clumps using the data obtained from the ALMA Survey of 70 $\mu \rm m$ dark High-mass clumps in Early Stages (ASHES). We identified 97 protostellar cores and 197 prestellar core candidates based on the detection of outflows and molecular transitions of high upper energy levels ($E_{u}/k > 45$ K). The detection rate of the N$_{2}$D$^{+}$ emission toward the protostellar cores is 38%, which is higher than 9% for the prestellar cores, indicating that N$_{2}$D$^{+}$ does not exclusively trace prestellar cores. The detection rates of the DCO$^{+}$ emission are 35% for the prestellar cores and 49% for the protostellar cores, which are higher than those of N$_{2}$D$^{+}$, implying that DCO$^{+}$ appears more frequently than N$_{2}$D$^{+}$ in both prestellar and protostellar cores. Both N$_{2}$D$^{+}$ and DCO$^{+}$ abundances appear to decrease from the prestellar to protostellar stage. The DCN, C$_{2}$D and $^{13}$CS emission lines are rarely seen in the dense cores of early evolutionary phases. The detection rate of the H$_{2}$CO emission toward dense cores is 52%, three times higher than that of CH$_{3}$OH (17%). In addition, the H$_{2}$CO detection rate, abundance, line intensities, and line widths increase with the core evolutionary status, suggesting that the H$_{2}$CO line emission is sensitive to protostellar activity., Comment: Accepted for Publication in ApJ. 36 pages, 15 figures, 5 tables. Table 5 available at https://github.com/ShanghuoLi/70-um-dark-high-mass-clumps

Studying the physical and chemical properties of cold and dense molecular clouds is crucial for the understanding of how stars form. Under the typical conditions of infrared dark clouds, CO is removed from the gas phase and trapped on to the surface of dust grains by the so-called depletion process. This suggests that the CO depletion factor ($f_{\rm D}$) can be a useful chemical indicator for identifying cold and dense regions (i.e., prestellar cores). We have used the 1.3 mm continuum and C$^{18}$O(2-1) data observed at the resolution of $\sim$5000 au in the ALMA Survey of 70 $\mu$m Dark High-mass Clumps in Early Stages (ASHES) to construct averaged maps of $f_{\rm D}$ in twelve clumps to characterise the earliest stages of the high-mass star formation process. The average $f_{\rm D}$ determined for 277 of the 294 ASHES cores follows an unexpected increase from the prestellar to the protostellar stage. If we exclude the temperature effect due to the slight variations in the NH$_3$ kinetic temperature among different cores, we explain this result as a dependence primarily on the average gas density, which increases in cores where protostellar conditions prevail. This shows that $f_{\rm D}$ determined in high-mass star-forming regions at the core scale is insufficient to distinguish among prestellar and protostellar conditions for the individual cores, and should be complemented by information provided by additional tracers. However, we confirm that the clump-averaged $f_{\rm D}$ values correlates with the luminosity-to-mass ratio of each source, which is known to trace the evolution of the star formation process., Comment: 20 pages, 10 figures, 2 tables. Accepted for publication in The Astrophysical Journal (ApJ)

In this study, 36 cores (30 starless and 6 protostellar) identified in Orion were surveyed to search for inward motions. We used the Nobeyama 45 m radio telescope, and mapped the cores in the $J = 1\rightarrow0$ transitions of HCO$^+$, H$^{13}$CO$^+$, N$_2$H$^+$, HNC, and HN$^{13}$C. The asymmetry parameter $\delta V$, which was the ratio of the difference between the HCO$^+$ and H$^{13}$CO$^+$ peak velocities to the H$^{13}$CO$^+$ line width, was biased toward negative values, suggesting that inward motions were more dominant than outward motions. Three starless cores (10% of all starless cores surveyed) were identified as cores with blue-skewed line profiles (asymmetric profiles with more intense blue-shifted emission), and another two starless cores (7%) were identified as candidate blue-skewed line profiles. The peak velocity difference between HCO$^+$ and H$^{13}$CO$^+$ of them was up to 0.9 km s$^{-1}$, suggesting that some inward motions exceeded the speed of sound for the quiescent gas ($\sim10-17$ K). The mean of $\delta V$ of the five aforementioned starless cores was derived to be $-$0.5$\pm$0.3. One core, G211.16$-$19.33North3, observed using the ALMA ACA in DCO$^+$ $J = 3\rightarrow2$ exhibited blue-skewed features. Velocity offset in the blue-skewed line profile with a dip in the DCO$^+$ $J = 3\rightarrow2$ line was larger ($\sim 0.5$ km s$^{-1}$) than that in HCO$^+$ $J = 1\rightarrow0$ ($\sim 0.2$ km s$^{-1}$), which may represent gravitational acceleration of inward motions. It seems that this core is at the last stage in the starless phase, judging from the chemical evolution factor version 2.0 (CEF2.0)., Comment: 45 pages, 28 figures, ApJ in press

We observed the high-mass protostellar core G335.579-0.272 ALMA1 at ${\sim}200$ au (0.05") resolution with the Atacama Large Millimeter/submillimeter Array (ALMA) at 226 GHz (with a mass sensitivity of $5\sigma=0.2$ M$_\odot$ at 10 K). We discovered that at least a binary system is forming inside this region, with an additional nearby bow-like structure (${\lesssim}1000$ au) that could add an additional member to the stellar system. These three sources are located at the center of the gravitational potential well of the ALMA1 region and the larger MM1 cluster. The emission from CH$_3$OH (and many other tracers) is extended ($>1000$ au), revealing a common envelope toward the binary system. We use CH$_2$CHCN line emission to estimate an inclination angle of the rotation axis of $26^\circ$ with respect to the line of sight based on geometric assumptions and derive a kinematic mass of the primary source (protostar+disk) of 3.0 M$_\odot$ within a radius of 230 au. Using SiO emission, we find that the primary source drives the large scale outflow revealed by previous observations. Precession of the binary system likely produces a change in orientation between the outflow at small scales observed here and large scales observed in previous works. The bow structure may have originated by entrainment of matter into the envelope due to widening or precession of the outflow, or, alternatively, an accretion streamer dominated by the gravity of the central sources. An additional third source, forming due to instabilities in the streamer, cannot be ruled out as a temperature gradient is needed to produce the observed absorption spectra., Comment: 23 pages, 16 figures, 4 tables, 2 appendices. Accepted for publication in ApJ

We report a study of the low-mass Class-0 multiple system VLA 1623AB in the Ophiuchus star-forming region, using H$^{13}$CO$^+$ ($J=3-2$), CS ($J=5-4$), and CCH ($N=3-2$) lines as part of the ALMA Large Program FAUST. The analysis of the velocity fields revealed the rotation motion in the envelope and the velocity gradients in the outflows (about 2000 au down to 50 au). We further investigated the rotation of the circum-binary VLA 1623A disk as well as the VLA 1623B disk. We found that the minor axis of the circum-binary disk of VLA 1623A is misaligned by about 12 degrees with respect to the large-scale outflow and the rotation axis of the envelope. In contrast, the minor axis of the circum-binary disk is parallel to the large-scale magnetic field according to previous dust polarization observations, suggesting that the misalignment may be caused by the different directions of the envelope rotation and the magnetic field. If the velocity gradient of the outflow is caused by rotation, the outflow has a constant angular momentum and the launching radius is estimated to be $5-16$ au, although it cannot be ruled out that the velocity gradient is driven by entrainments of the two high-velocity outflows. Furthermore, we detected for the first time a velocity gradient associated with rotation toward the VLA 16293B disk. The velocity gradient is opposite to the one from the large-scale envelope, outflow, and circum-binary disk. The origin of its opposite gradient is also discussed., Comment: 27 pages, 21 figures, 2 Tables, Accepted for publication in ApJ

We have observed the 70 $\mu$m dark infrared dark cloud (IRDC) G14.492-00.139 in the N$_2$D$^+$ $J$=3--2, DCO$^+$ $J$=3--2, DCN $J$=3--2, and C$^{18}$O $J$=2--1 lines, using the Atacama Large Millimeter/submillimeter Array (ALMA) as part of the ALMA Survey of 70 $\mu$m Dark High-mass Clumps in Early Stages (ASHES). We find that the spatial distribution is different among the observed emission from the deuterated molecular lines. The N$_2$D$^+$ emission traces relatively quiescent regions, while both the DCO$^+$ and DCN emission emanates mainly from regions with signs of active star-formation. In addition, the DCO$^+$/N$_2$D$^+$ ratio is found to be lower in several dense cores than in starless cores embedded in low-mass star-forming regions. By comparing the observational results with chemical model calculations, we discuss the origin of the low DCO$^+$/N$_2$D$^+$ ratio in this IRDC clump. The low DCO$^+$/N$_2$D$^+$ ratio can be explained if the temperature of the dense cores is in the range between the sublimation temperature of N$_2$ ($\sim$20 K) and CO ($\sim$25 K). The results suggest that the dense cores in G14.492-00.139 are warmer and denser than the dense cores in low-mass star-forming regions., Comment: accepted for publication in ApJ

With a mass of $\sim$1000 $M_\odot$ and a surface density of $\sim$0.5 g cm$^{-2}$, G023.477+0.114 also known as IRDC 18310-4 is an infrared dark cloud (IRDC) that has the potential to form high-mass stars and has been recognized as a promising prestellar clump candidate. To characterize the early stages of high-mass star formation, we have observed G023.477+0.114 as part of the ALMA Survey of 70 $\mu$m Dark High-mass Clumps in Early Stages (ASHES). We have conducted $\sim$1."2 resolution observations with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.3 mm in dust continuum and molecular line emission. We identified 11 cores, whose masses range from 1.1 $M_\odot$ to 19.0 $M_\odot$. Ignoring magnetic fields, the virial parameters of the cores are below unity, implying that the cores are gravitationally bound. However, when magnetic fields are included, the prestellar cores are close to virial equilibrium, while the protostellar cores remain sub-virialized. Star formation activity has already started in this clump. Four collimated outflows are detected in CO and SiO. H$_2$CO and CH$_3$OH emission coincide with the high-velocity components seen in the CO and SiO emission. The outflows are randomly oriented for the natal filament and the magnetic field. The position-velocity diagrams suggest that episodic mass ejection has already begun even in this very early phase of protostellar formation. The masses of the identified cores are comparable to the expected maximum stellar mass that this IRDC could form (8-19 $M_\odot$). We explore two possibilities on how IRDC G023.477+0.114 could eventually form high-mass stars in the context of theoretical scenarios., Comment: Accepted for publication in ApJ (September 1, 2021). 33 pages, 20 figures, and 5 tables

We have observed the young protostellar system NGC 2264 CMM3 in the 1.3 mm and 2.0 mm bands at a resolution of about 0.1$"$ (70 au) with ALMA. The structures of two distinct components, CMM3A and CMM3B, are resolved in the continuum images of both bands. CMM3A has an elliptical structure extending along the direction almost perpendicular to the known outflow, while CMM3B reveals a round shape. We have fitted two 2D-Gaussian components to the elliptical structure of CMM3A and CMM3B, and have separated the disk and envelope components for each source. The spectral index $\alpha$ between 2.0 mm and 0.8 mm is derived to be 2.4-2.7 and 2.4-2.6 for CMM3A and CMM3B, respectively, indicating the optically thick dust emission and/or the grain growth. A velocity gradient in the disk/envelope direction is detected for CMM3A in the CH$_3$CN, CH$_3$OH, and $^{13}$CH$_3$OH lines detected in the 1.3 mm band, which can be interpreted as the rotation of the disk/envelope system. From this result, the protostellar mass of CMM3A is roughly evaluated to be $0.1- 0.5$ $M_\odot$ by assuming Keplerian rotation. The mass accretion rate is thus estimated to be $5\times10^{-5}$ - 4 $\times$ $10^{-3}$ $M_\odot$ yr$^{-1}$, which is higher than typical mass accretion rate of low-mass protostars. The OCS emission line shows a velocity gradient in both outflow direction and disk/envelope direction. A hint of outflow rotation is found, and the specific angular momentum of the outflow is estimated to be comparable to that of the disk. These results provide us with novel information on the initial stage of a binary/multiple system., Comment: Accepted for publication in the Astrophysical Journal, 21 pages, 12 figures

Temperature and pressure are essential parameters in the synthesis, evaluation, and application of functional materials. This study proposes the addition of a heating function to a high-pressure diamond anvil cell (DAC) with in-situ measurement probes. The proposed DAC allows for simultaneous control of temperature and pressure within the sample space and can be used to synthesize functional materials under extreme conditions. The various components, namely the heater, thermometer, and measurement probes, were fabricated with a boron-doped diamond epitaxial film and could be used repeatedly. The developed DAC was used to successfully conduct the high-pressure annealing of La(O,F)BiS$_{2}$ single crystal and the high-pressure synthesis of EuFBiS$_{2}$ superconductors. The proposed technique shows promise for further exploration of superconductors to broaden the research field.

We present the results of on-the-fly mapping observations of 44 fields containing 107 SCUBA-2 cores in the emission lines of molecules, N$_2$H$^+$, HC$_3$N, and CCS at 82$-$94 GHz using the Nobeyama 45-m telescope. This study aimed at investigating the physical properties of cores that show high deuterium fractions and might be close to the onset of star formation. We found that the distributions of the N$_2$H$^+$ and HC$_3$N line emissions are approximately similar to that of 850-$\mu$m dust continuum emission, whereas the CCS line emission is often undetected or is distributed in a clumpy structure surrounding the peak position of the 850-$\mu$m dust continuum emission. Occasionally (12%), we observe the CCS emission which is an early-type gas tracer toward the young stellar object, probably due to local high excitation. Evolution toward star formation does not immediately affect nonthermal velocity dispersion., Comment: 45 pages, 15 figures, ApJS, in press

A full understanding of high-mass star formation requires the study of one of the most elusive components of the energy balance in the interstellar medium: magnetic fields. We report ALMA 1.2 mm, high-resolution (700 au) dust polarization and molecular line observations of the rotating hot molecular core embedded in the high-mass star-forming region IRAS 18089-1732. The dust continuum emission and magnetic field morphology present spiral-like features resembling a whirlpool. The velocity field traced by the H13CO+ (J=3-2) transition line reveals a complex structure with spiral filaments that are likely infalling and rotating, dragging the field with them. We have modeled the magnetic field and find that the best model corresponds to a weakly magnetized core with a mass-to-magnetic-flux ratio (lambda) of 8.38. The modeled magnetic field is dominated by a poloidal component, but with an important contribution from the toroidal component that has a magnitude of 30% of the poloidal component. Using the Davis-Chandrasekhar-Fermi method, we estimate a magnetic field strength of 3.5 mG. At the spatial scales accessible to ALMA, an analysis of the energy balance of the system indicates that gravity overwhelms turbulence, rotation, and the magnetic field. We show that high-mass star formation can occur in weakly magnetized environments, with gravity taking the dominant role., Comment: Accepted for publications in ApJL (12 pages, 5 figures, Appendix)

We present a spatio-kinematical analysis of the CO~($J$=2$\rightarrow$1) line emission, observed with the Atacama Large Millimter/submillimter Array (ALMA), of the outflow associated with the most massive core, ALMA1, in the 70 $\mu$m dark clump G010.991$-$00.082. The position-velocity (P-V) diagram of the molecular outflow exhibits a peculiar $\mathsf{S}$-shaped morphology that has not been seen in any other star forming region. We propose a spatio-kinematical model for the bipolar molecular outflow that consists of a decelerating high-velocity component surrounded by a slower component whose velocity increases with distance from the central source. The physical interpretation of the model is in terms of a jet that decelerates as it entrains material from the ambient medium, which has been predicted by calculations and numerical simulations of molecular outflows in the past. One side of the outflow is shorter and shows a stronger deceleration, suggesting that the medium through which the jet moves is significantly inhomogeneous. The age of the outflow is estimated to be $\tau$$\approx$1300 years, after correction for a mean inclination of the system of $\approx$57$^{\circ}$., Comment: 15 pages, 7 figures. Accepted for publication in ApJ

To date, about two dozen low-mass embedded protostars exhibit rich spectra with lines of complex organic molecule (COM). These protostars seem to possess different enrichment in COMs. However, the statistics of COM abundance in low-mass protostars are limited by the scarcity of observations. This study introduces the Perseus ALMA Chemistry Survey (PEACHES), which aims at unbiasedly characterizing the chemistry of COMs toward the embedded (Class 0/I) protostars in the Perseus molecular cloud. Of 50 embedded protostars surveyed, 58% of them have emission from COMs. A 56%, 32%, and 40% of the protostars have CH$_3$OH, CH$_3$OCHO, and N-bearing COMs, respectively. The detectability of COMs depends neither on the averaged continuum brightness temperature, a proxy of the H$_2$ column density, nor on the bolometric luminosity and the bolometric temperature. For the protostars with detected COMs, CH$_3$OH has a tight correlation with CH$_3$CN, spanning more than two orders of magnitude in column densities normalized by the continuum brightness temperature, suggesting a chemical relation between CH$_3$OH and CH$_3$CN and a large chemical diversity in the PEACHES samples at the same time. A similar trend with more scatter is also found between all identified COMs, hinting at a common chemistry for the sources with COMs. The correlation between COMs is insensitive to the protostellar properties, such as the bolometric luminosity and the bolometric temperature. The abundance of larger COMs (CH$_3$OCHO and CH$_3$OCH$_3$) relative to that of smaller COMs (CH$_3$OH and CH$_3$CN) increases with the inferred gas column density, hinting at an efficient production of complex species in denser envelopes., Comment: 55 pages, 9 tables, 24 figures. Updated after publication in ApJ. Figure 18 appears as a figure set in ApJ

We observed the high-mass star-forming region G335.579-0.292 with the Atacama Large Millimeter/submillimeter Array (ALMA) at 226 GHz with an angular resolution of 0.3'' ($\sim 1000$ au resolution at the source distance). G335.579-0.292 hosts one of the most massive cores in the Galaxy (G335-MM1). The continuum emission shows that G335-MM1 fragments into at least five sources, while molecular line emission is detected in two of the continuum sources (ALMA1 and ALMA3). We found evidence of large and small scale infall in ALMA1 revealed by an inverse P-Cygni profile and the presence of a blue-shifted spot at the center of the first moment map of the CH$_3$CN emission. In addition, hot gas expansion in the innermost region is unveiled by a red-shifted spot in the first moment map of HDCO and (CH$_3$)$_2$CO (both with $E_u > 1100$ K). Our modeling reveals that this expansion motion originates close to the central source, likely due to reversal of the accretion flow induced by the expansion of the HII region, while infall and rotation motions originate in the outer regions. ALMA3 shows clear signs of rotation, with a rotation axis inclination with respect to the line of sight close to $90^\circ$, and a system mass (disk + star) in the range of 10-30 M$_\odot$., Comment: 15 pages, 12 figures, 3 tables. Accepted for publication in ApJ

We have observed the very low-mass Class 0 protostar IRAS 15398-3359 at scales ranging from 50 au to 1800 au, as part of the ALMA Large Program FAUST. We uncover a linear feature, visible in H2CO, SO, and C18O line emission, which extends from the source along a direction almost perpendicular to the known active outflow. Molecular line emission from H2CO, SO, SiO, and CH3OH further reveals an arc-like structure connected to the outer end of the linear feature and separated from the protostar, IRAS 15398-3359, by 1200 au. The arc-like structure is blue-shifted with respect to the systemic velocity. A velocity gradient of 1.2 km/s over 1200 au along the linear feature seen in the H2CO emission connects the protostar and the arc-like structure kinematically. SO, SiO, and CH3OH are known to trace shocks, and we interpret the arc-like structure as a relic shock region produced by an outflow previously launched by IRAS 15398-3359. The velocity gradient along the linear structure can be explained as relic outflow motion. The origins of the newly observed arc-like structure and extended linear feature are discussed in relation to turbulent motions within the protostellar core and episodic accretion events during the earliest stage of protostellar evolution.

Recent spectroscopic observations by sensitive radio telescopes require accurate molecular spectral line frequencies to identify molecular species in a forest of lines detected. To measure rest frequencies of molecular spectral lines in the laboratory, an emission-type millimeter and submillimeter-wave spectrometer utilizing state-of-the-art radio-astronomical technologies is developed. The spectrometer is equipped with a 200 cm glass cylinder cell, a two sideband (2SB) Superconductor-Insulator-Superconductor (SIS) receiver in the 230 GHz band, and wide-band auto-correlation digital spectrometers. By using the four 2.5 GHz digital spectrometers, a total instantaneous bandwidth of the 2SB SIS receiver of 8 GHz can be covered with a frequency resolution of 88.5 kHz. Spectroscopic measurements of CH$_3$CN and HDO are carried out in the 230 GHz band so as to examine frequency accuracy, stability, sensitivity, as well as intensity calibration accuracy of our system. As for the result of CH$_3$CN, we confirm that the frequency accuracy for lines detected with sufficient signal to noise ratio is better than 1 kHz, when the high resolution spectrometer having a channel resolution of 17.7 kHz is used. In addition, we demonstrate the capability of this system by spectral scan measurement of CH$_3$OH from 216 GHz to 264 GHz. We assign 242 transitions of CH$_3$OH, 51 transitions of $^{13}$CH$_3$OH, and 21 unidentified emission lines for 295 detected lines. Consequently, our spectrometer demonstrates sufficient sensitivity, spectral resolution, and frequency accuracy for in-situ experimental-based rest frequency measurements of spectral lines on various molecular species., Comment: Accepted for publication in Publications of the Astronomical Society of Japan (PASJ) , 15 figures, 4 tables, typos corrected

We report a statistical analysis exploring the origin of the overall low star formation efficiency (SFE) of the Galactic central molecular zone (CMZ) and the SFE diversity among the CMZ clouds using a wide-field HCN $J$=4-3 map, whose optically thin critical density ($\sim10^7\,\mathrm{cm}^{-3}$) is the highest among the tracers ever used in CMZ surveys. Logistic regression is performed to empirically formulate star formation probability of 195 HCN clumps, 13 of which contain star formation signatures. The explanatory parameters in the best-fit model are reduced into the virial parameter $\alpha_{\mathrm{vir}}$ without significant contribution from other parameters, whereas the performance of the model without $\alpha_{\mathrm{vir}}$ is no better than that using randomly generated data. The threshold $\alpha_{\mathrm{vir}}$ is 6, which translates into a volume density ($n_{\mathrm{H_2}}$) of $10^{4.6}\,\mathrm{cm}^{-3}$ with the $n_{\mathrm{H_2}}$-$\alpha_{\mathrm{vir}}$ correlation. The scarcity of the low-$\alpha_{\mathrm{vir}}$ clumps, whose fraction to all HCN clumps is 0.1, can be considered as one of the immediate causes of the suppressed SFE. No correlation between the clump size or mass and star formation probability is found, implying that HCN $J$=4-3 does not immediately trace the mass of star-forming gas above a threshold density. Meanwhile, star-forming and non-star-forming clouds are degenerate in the physical parameters of the CS $\mathit{J}$=1-0 clouds, highlighting the efficacy of the HCN $\mathit{J}$=4-3 line to probe star-forming regions in the CMZ. The time scale of the high-$\alpha_{\mathrm{vir}}$ to low-$\alpha_{\mathrm{vir}}$ transition is $\lesssim2$ Myr, which is consistent with the tidal compression and X1/X2 orbit transition models but possibly does not fit the cloud-cloud collision picture., Comment: 25 pages, 7 figures, 5 tables, accepted for publication in the Astrophysical Journal

We present a study of outflows at extremely early stages of high-mass star formation obtained from the ALMA Survey of 70 $\mu \rm m$ dark High-mass clumps in Early Stages (ASHES). Twelve massive 3.6$-$70 $\mu \rm m$ dark prestellar clump candidates were observed with the Atacama Large Millimeter/submillimeter Array (ALMA) in Band 6. Forty-three outflows are identified toward 41 out of 301 dense cores using the CO and SiO emission lines, yielding a detection rate of 14%. We discover 6 episodic molecular outflows associated with low- to high-mass cores, indicating that episodic outflows (and therefore episodic accretion) begin at extremely early stages of protostellar evolution for a range of core masses. The time span between consecutive ejection events is much smaller than those found in more evolved stages, which indicates that the ejection episodicity timescale is likely not constant over time. The estimated outflow dynamical timescale appears to increase with core masses, which likely indicates that more massive cores have longer accretion timescales than less massive cores. The lower accretion rates in these 70 $\mu \rm m$ dark objects compared to the more evolved protostars indicate that the accretion rates increase with time. The total outflow energy rate is smaller than the turbulent energy dissipation rate, which suggests that outflow induced turbulence cannot sustain the internal clump turbulence at the current epoch. We often detect thermal SiO emission within these 70 $\mu \rm m$ dark clumps that is unrelated to CO outflows. This SiO emission could be produced by collisions, intersection flows, undetected protostars, or other motions., Comment: 32 pages, 9 figures, 4 tables, accepted for publication in ApJ

We present the results of a single-pointing survey of 207 dense cores embedded in Planck Galactic Cold Clumps distributed in five different environments ($\lambda$ Orionis, Orion A, B, Galactic plane, and high latitudes) to identify dense cores on the verge of star formation for the study of the initial conditions of star formation. We observed these cores in eight molecular lines at 76-94 GHz using the Nobeyama 45-m telescope. We find that early-type molecules (e.g., CCS) have low detection rates and that late-type molecules (e.g., N$_2$H$^+$, c-C$_3$H$_2$) and deuterated molecules (e.g., N$_2$D$^+$, DNC) have high detection rates, suggesting that most of the cores are chemically evolved. The deuterium fraction (D/H) is found to decrease with increasing distance, indicating that it suffers from differential beam dilution between the D/H pair of lines for distant cores ($>$1 kpc). For $\lambda$ Orionis, Orion A, and B located at similar distances, D/H is not significantly different, suggesting that there is no systematic difference in the observed chemical properties among these three regions. We identify at least eight high D/H cores in the Orion region and two at high latitudes, which are most likely to be close to the onset of star formation. There is no clear evidence of the evolutionary change in turbulence during the starless phase, suggesting that the dissipation of turbulence is not a major mechanism for the beginning of star formation as judged from observations with a beam size of 0.04 pc., Comment: 55 pages, 12 Figures, 8 Tables, accepted for publication in ApJS

We mapped two molecular cloud cores in the Orion A cloud with the ALMA ACA 7-m Array and with the Nobeyama 45-m radio telescope. These cores have bright N$_2$D$^+$ emission in single-pointing observations with the Nobeyama 45-m radio telescope, have relatively high deuterium fraction, and are thought to be close to the onset of star formation. One is a star-forming core, and the other is starless. These cores are located along filaments observed in N$_2$H$^+$, and show narrow linewidths of 0.41 km s$^{-1}$ and 0.45 km s$^{-1}$ in N$_2$D$^+$, respectively, with the Nobeyama 45-m telescope. Both cores were detected with the ALMA ACA 7m Array in the continuum and molecular lines at Band 6. The starless core G211 shows clumpy structure with several sub-cores, which in turn show chemical differences. Also, the sub-cores in G211 have internal motions that are almost purely thermal. The starless sub-core G211D, in particular, shows a hint of the inverse P Cygni profile, suggesting infall motion. The star-forming core G210 shows an interesting spatial feature of two N$_2$D$^+$ peaks of similar intensity and radial velocity located symmetrically with respect to the single dust continuum peak. One interpretation is that the two N$_2$D$^+$ peaks represent an edge-on pseudo-disk. The CO outflow lobes, however, are not directed perpendicular to the line connecting both N$_2$D$^+$ peaks., Comment: 31 pages, 17 figures, 9 tables, to appear in ApJ (accepted April 23, 2020)

We propose a new observing method for single-dish millimeter and submillimeter spectroscopy using a heterodyne receiver equipped with a frequency-modulating local oscillator (FMLO). Unlike conventional switching methods, which extract astronomical signals by subtracting the reference spectra of off-sources from those of on-sources, the FMLO method does not need to obtain any off-source spectra; rather, it estimates them from the on-source spectra themselves. The principle is a high dump-rate (10 Hz) spectroscopy with radio frequency modulation (FM) achieved by fast sweeping of a local oscillator (LO) of a heterodyne receiver: Because sky emission (i.e., off-source) fluctuates as $1/f$-type and is spectrally correlated, it can be estimated and subtracted from time-series spectra (a timestream) by principal component analysis. Meanwhile astronomical signals remain in the timestream since they are modulated to a higher time-frequency domain. The FMLO method therefore achieves (1) a remarkably high observation efficiency, (2) reduced spectral baseline wiggles, and (3) software-based sideband separation. We developed an FMLO system for the Nobeyama 45-m telescope and a data reduction procedure for it. Frequency modulation was realized by a tunable and programmable first local oscillator. With observations of Galactic sources, we demonstrate that the observation efficiency of the FMLO method is dramatically improved compared to conventional switching methods. Specifically, we find that the time to achieve the same noise level is reduced by a factor of 3.0 in single-pointed observations and by a factor of 1.2 in mapping observations. The FMLO method can be applied to observations of fainter ($\sim$mK) spectral lines and larger ($\sim$deg$^{2}$) mapping. It would offer much more efficient and baseline-stable observations compared to conventional switching methods., Comment: 22 pages, 15 figures, 2 tables, accepted for publication in Publications of the Astronomical Society of Japan

(Abridged) ASHES has been designed to systematically characterize the earliest stages and to constrain theories of high-mass star formation. A total of 12 massive (>500 $M_{\odot}$), cold (<15 K), 3.6-70 $\mu$m dark prestellar clump candidates, embedded in IRDCs, were carefully selected in the pilot survey to be observed with ALMA. We mosaiced each clump (~1 arcmin^2) in dust and line emission with the 12m/7m/TP arrays at 224 GHz, resulting in ~1.2" resolution (~4800 AU). As the first paper of the series, we concentrate on the dust emission to reveal the clump fragmentation. We detect 294 cores, from which 84 (29%) are categorized as protostellar based on outflow activity or 'warm core' line emission. The remaining 210 (71%) are considered prestellar core candidates. The number of detected cores is independent of the mass sensitivity range of the observations. On average, more massive clumps tend to form more cores. We find a large population of low-mass (<1 M) cores and no high-mass (>30 $M_{\odot}$) prestellar cores. The most massive prestellar core has a mass of 11 $M_{\odot}$. From the prestellar CMF, we derive a power law index of 1.17+-0.1, slightly shallower than Salpeter (1.35). We use the MST technique to characterize the separation between cores and their spatial distribution, and derive mass segregation ratios. While there is a range of core masses and separations detected in the sample, the mean separation and mass of cores are well explained by thermal fragmentation and are inconsistent with turbulent Jeans fragmentation. The core spatial distribution is well described by hierarchical subclustering rather than centrally peaked clustering. There is no conclusive evidence of mass segregation. We test several theoretical conditions, and conclude that overall, competitive accretion and global hierarchical collapse scenarios are favored over the turbulent core accretion scenario., Comment: Accepted for publication in ApJ (September 17, 2019). 37 pages, 19 figures, and 7 tables

We have observed the Class I protostellar source Elias 29 with Atacama Large Millimeter/submillimeter Array (ALMA). We have detected CS, SO, $^{34}$SO, SO$_2$, and SiO line emissions in a compact component concentrated near the protostar and a ridge component separated from the protostar by 4\arcsec\ ($\sim 500$ au). The former component is found to be abundant in SO and SO$_2$ but deficient in CS. The abundance ratio SO/CS is as high as $3^{+13}_{-2} \times 10^2$ at the protostar, which is even higher than that in the outflow-shocked region of L1157 B1. However, organic molecules (HCOOCH$_3$, CH$_3$OCH$_3$, CCH, and c-C$_3$H$_2$) are deficient in Elias 29. We attribute the deficiency in organic molecules and richness in SO and SO$_2$ to the evolved nature of the source or the relatively high dust temperature (\protect\raisebox{-0.7ex}{$\:\stackrel{\textstyle >}{\sim}\:$} 20 K) in the parent cloud of Elias 29. The SO and SO$_2$ emissions trace rotation around the protostar. Assuming a highly inclined configuration ($i \geq 65$\degr; 0\degr\ for a face-on configuration) and Keplerian motion for simplicity, the protostellar mass is estimated to be (0.8 -- 1.0) \Msun. The $^{34}$SO and SO$_2$ emissions are asymmetric in their spectra; the blue-shifted components are weaker than the red-shifted ones. Although this may be attributed to the asymmetric molecular distribution, other possibilities are also discussed.

An unbiased spectral line survey toward a solar-type Class 0/I protostar, IRAS04368+2557, in L1527 has been carried out in the 3 mm band with the Nobeyama 45 m telescope. L1527 is known as a warm carbon-chain chemistry (WCCC) source, which harbors abundant unsaturated organic species such as C$_n$H ($n = 3,\ 4,\ 5,\ldots$) in a warm and dense region near the protostar. The observation covers the frequency range from 80 to 116 GHz. A supplementary observation has also been conducted in the 70 GHz band to observe fundamental transitions of deuterated species. In total, 69 molecular species are identified, among which 27 species are carbon-chain species and their isomers, including their minor isotopologues. This spectral line survey provides us with a good template of the chemical composition of the WCCC source., Comment: published online in PASJ

In this study, we designed and experimentally evaluated a series-connected array of superconductor-insulator-superconductor (SIS) junctions in the 100-GHz band mixer for the multi-beam receiver FOREST on the Nobeyama 45-m millimeter-wave telescope. The construction of the junction chip comprised a waveguide probe antenna, impedance matching circuit, SIS array junction, and choke filter, which were made from a superconducting niobium planar circuit on a quartz substrate. The multi-stage impedance matching circuit between the feed point and the SIS junction was designed as a capacitively loaded transmission line, and it comprised two sections with high (~90 Ohm) and low (~10 Ohm) characteristic impedance transmission lines. The structure of this tuning line was simple and easy to fabricate, and the feed impedance matched with the SIS junction in a wide frequency range. The signal coupling efficiency was more than 92% and the expected receiver noise temperature was approximately two times the quantum limit for 75-125 GHz based on quantum theory. The array junction devices with 3-6 connected junctions were fabricated and we measured their performance in terms of the receiver noise temperature and gain compression in the laboratory. We successfully developed an array junction device with a receiver noise temperature of ~15-30 K and confirmed that the improvement in the saturation power corresponded to the number of junctions. The newly developed array junction mixer was installed in the FOREST receiver and it successfully detected the 12CO (J = 1-0) molecular line toward IRC+10216 with the Nobeyama 45-m telescope., Comment: 10 pasges, 8 figures, Accepted for publication in PASJ

We report molecular line observations of the NGC 1333 IRAS 4C outflow in the Perseus Molecular Cloud with the Atacama Large Millimeter/Submillimeter Array. The CCH and CS emission reveal an outflow cavity structure with clear signatures of rotation with respect to the outflow axis. The rotation is detected from about 120 au up to about 1400 au above the envelope/disk mid-plane. As the distance to the central source increases, the rotation velocity of the outflow decreases while the outflow radius increases, which gives a flat specific angular momentum distribution along the outflow. The mean specific angular momentum of the outflow is about 100 au km/s. Based on reasonable assumptions on the outward velocity of the outflow and the protostar mass, we estimate the range of outflow launching radii to be 5-15 au. Such a launching radius rules out that this outflow is launched as an X-wind, but rather, it is more consistent to be a slow disk wind launched from relatively large radii on the disk. The radius of the centrifugal barrier is roughly estimated, and the role of the centrifugal barrier in the outflow launching is discussed., Comment: Accepted to ApJ. 29 pages, 8 figures

We performed sensitive observations of the N15ND+(1-0) and 15NND+(1-0) lines toward the prestellar core L1544 using the IRAM 30m telescope. The lines are not detected down to 3 sigma levels in 0.2 km/s channels of around 6 mK. The non-detection provides the lower limit of the 14N/15N ratio for N2D+ of ~700-800, which is much higher than the elemental abundance ratio in the local ISM of ~200-300. The result indicates that N2 is depleted in 15N in the central part of L1544, because N2D+ preferentially traces the cold dense gas, and because it is a daughter molecule of N2. In-situ chemistry is unlikely responsible for the 15N depletion in N2; neither low-temperature gas phase chemistry nor isotope selective photodissociation of N2 explains the 15N depletion; the former prefers transferring 15N to N2, while the latter requires the penetration of interstellar FUV photons into the core center. The most likely explanation is that 15N is preferentially partitioned into ices compared to 14N via the combination of isotope selective photodissociation of N2 and grain surface chemistry in the parent cloud of L1544 or in the outer regions of L1544 which are not fully shielded from the interstellar FUV radiation. The mechanism is the most efficient at the chemical transition from atomic to molecular nitrogen. In other words, our result suggests that the gas in the central part of L1544 already went trough the transition from atomic to molecular nitrogen in the earlier evolutionary stage, and that N2 is currently the primary form of gas-phase nitrogen., Comment: 5 pages, 2 figures, 2 tables, Accepted for publication in A&A Letters

Using Galactic Plane surveys, we have selected a massive (1200 M$_\odot$), cold (14 K) 3.6-70 $\mu$m dark IRDC G331.372-00.116. This IRDC has the potential to form high-mass stars and, given the absence of current star formation signatures, it seems to represent the earliest stages of high-mass star formation. We have mapped the whole IRDC with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.1 and 1.3 mm in dust continuum and line emission. The dust continuum reveals 22 cores distributed across the IRDC. In this work, we analyze the physical properties of the most massive core, ALMA1, which has no molecular outflows detected in the CO (2-1), SiO (5-4), and H$_2$CO (3-2) lines. This core is relatively massive ($M$ = 17.6 M$_\odot$), subvirialized (virial parameter $\alpha_{vir}=M_{vir}/M=0.14$), and is barely affected by turbulence (transonic Mach number of 1.2). Using the HCO$^+$ (3-2) line, we find the first detection of infall signatures in a relatively massive, prestellar core (ALMA1) with the potential to form a high-mass star. We estimate an infall speed of 1.54 km s$^{-1}$ and a high accretion rate of 1.96 $\times$ 10$^{-3}$ M$_\odot$ yr$^{-1}$. ALMA1 is rapidly collapsing, out of virial equilibrium, more consistent with competitive accretion scenarios rather than the turbulent core accretion model. On the other hand, ALMA1 has a mass $\sim$6 times larger than the clumps Jeans mass, being in an intermediate mass regime ($M_{J}=2.7

Chemical diversity of the gas in low-mass protostellar cores is widely recognized. In order to explore its origin, a survey of chemical composition toward 36 Class 0/I protostars in the Perseus molecular cloud complex, which are selected in an unbiased way under certain physical conditions, has been conducted with IRAM 30 m and NRO 45 m telescope. Multiple lines of C2H, c-C3H2 and CH3OH have been observed to characterize the chemical composition averaged over a 1000 au scale around the protostar. The derived beam-averaged column densities show significant chemical diversity among the sources, where the column density ratios of C2H/CH3OH are spread out by 2 orders of magnitude. From previous studies, the hot corino sources have abundant CH3OH but deficient C2H, their C2H/CH3OH column density ratios being relatively low. In contrast, the warm-carbon-chain chemistry (WCCC) sources are found to reveal the high C2H/CH3OH column density ratios. We find that the majority of the sources have intermediate characters between these two distinct chemistry types. A possible trend is seen between the C2H/CH3OH ratio and the distance of the source from the edge of a molecular cloud. The sources located near cloud edges or in isolated clouds tend to have a high C2H/CH3OH ratio. On the other hand, the sources having a low C2H/CH3OH ratio tend to be located in inner regions of the molecular cloud complex. This result gives an important clue to an understanding of the origin of the chemical diversity of protostellar cores in terms of environmental effects., Comment: Accepted for publication in ApJS

This {\it supplement} paper presents the maps of HCN $J$=4-3, HNC $J$=1-0, $\mathrm{H^{13}CN}$ $J$=1-0, and HC$_3$N $J$=10-9 for the Galactic central molecular zone (CMZ), which have been obtained using the Atacama Submillimeter Telescope Experiment and Nobeyama Radio Observatory 45-m telescope. Three-dimensional maps (2-D in space and 1-D in velocity) of the gas kinetic temperature ($T_\mathrm{kin}$), hydrogen volume density ($n_\mathrm{H_2}$), and fractional abundances of eight molecules (HCN, HNC, $\mathrm{HC_3N}$, HCO$^+$, $\mathrm{H_2CO}$, SiO, CS, and $\mathrm{N_2H^+}$) have been constructed from our and archival data. We have developed a method with hierarchical Bayesian inference for this analysis, which has successfully suppressed the artificial correlations among the parameters created by systematic errors due to the deficiency in the simple one-zone excitation analysis and the calibration uncertainty. The typical values of $T_\mathrm{kin}$ and $n_\mathrm{H_2}$ are $10^{1.8}\ $K and $10^{4.2}\ \mathrm{cm}^{-3}$, respectively, and the presence of an additional cold/low-density component is also indicated. The distribution of high-temperature regions is poorly correlated with known active star-forming regions, while a few of them coincide with shocked clouds. Principal component analysis has identified two distinct groups in the eight analyzed molecules: one group with large PC1 and PC2 scores and the other with a large $T_\mathrm{kin}$ dependence, which could be explained using two regimes of shock chemistry with fast ($\gtrsim 20\ \mathrm{km\,s}^{-1}$) and slow ($\lesssim 20\ \mathrm{km\,s}^{-1}$) velocity shocks, respectively. This supports the idea that the mechanical sputtering of dust grains and the mechanical heating play primary roles in the chemical and thermal processes in CMZ clouds., Comment: 51 pages, 58 figures, Accepted for publication in Astrophysical Journal Supplement Series

We have observed complex organic molecules (COMs) and deuterated species toward a hot core/corino (HC) associated with the infrared dark cloud (IRDC) clump G34.43+00.24 MM3 with the Atacama Large Millimeter/submillimeter Array (ALMA). We have detected six normal-COMs (CH3OH, CH3CHO, CH3CH2CN, CH3OCH3, HCOOCH3, and NH2CHO), one deuterated-COM (CH2DCN), and two deuterated fundamental molecules (D2CO and DNC) toward G34.43+00.24 MM3 HC. None of these lines, except for CH3OH, are detected toward the shocked regions in our data, which suggests that COMs do not originate in shocks. The abundance of the COMs relative to CH3OH in G34.43+00.24 MM3 HC is found to be similar to those in high-mass hot cores, rather than those in hot corinos in low-mass star-forming regions. This result suggests that the physical conditions of the warm-up phase of G34.43+00.24 MM3 HC are similar to those of high-mass sources. On the other hand, the D2CO abundance relative to CH$_3$OH in G34.43+00.24 MM3 HC is higher than that of other hot cores, and seems to be comparable to that of hot corinos. The relatively high D2CO/CH3OH ratio of G34.43+00.24 MM3 HC implies a long cold starless phase of G34.43+00.24 MM3 HC., Comment: accepted to ApJ